myEcalShower.cxx

Go to the documentation of this file.

// Authors: M.Duranti - INFN di Perugia
#include "myEcalShower.h"
#include "debug.h"
#include "TClass.h"

using namespace std;

//--------------------------------------------------------------------

ClassImp(myEcalShower);
#ifdef _WITHGBATCH_
ClassImp(myEcalShowerFiller);
#endif //#ifdef _WITHGBATCH_

//--------------------------------------------------------------------

myEcalShower::myEcalShower(){
#ifdef PDEBUG
  printf("In myEcalShower::myEcalShower\n");
#endif
  PRINTDEBUG;
  init();
  PRINTDEBUG;
}

myEcalShower::~myEcalShower(){
#ifdef PDEBUG
  printf("In myEcalShower::~myEcalShower\n");
#endif
  PRINTDEBUG;
  Clear();
}

void myEcalShower::Clear(Option_t* option){
#ifdef PDEBUG
  printf("In myEcalShower::Clear\n");
#endif
  PRINTDEBUG;
  Status = -999999;
  Nhits = -999999;
  /*
  memset(Dir, 0, 3*sizeof(Dir[0]));
  memset(CofG, 0, 3*sizeof(CofG[0]));
  memset(Entry, 0, 3*sizeof(Entry[0]));
  memset(Exit, 0, 3*sizeof(Exit[0]));
  */
  fill_n(Dir, 3, 0.);
  fill_n(CofG, 3, 0.);
  fill_n(Exit, 3, 0.);
  fill_n(Entry, 3, 0.);

  EnergyD = -0.0000001;
  EnergyC = -0.0000001;
  EnergyE = -0.0000001;
  EnergyE2014 = -0.0000001;
  EnergyA = -0.0000001;
  EnergyF = -0.0000001;
  DispLat = -999999;
  DispLatX = -999999;
  DispLatY = -999999;
  S35RaX = -999999;
  S35RaY = -999999;
  MaxS = -9999999;
  CorrectedEnergyElectrons = -9999999;
  CorrectedEnergyPhotons = -9999999;
  PRINTDEBUG;
  return;
}

void myEcalShower::init(){
#ifdef PDEBUG
  printf("In myEcalShower::init\n");
#endif
  PRINTDEBUG;
  Clear();
  PRINTDEBUG;
  return;
}

float myEcalShower::GetCorrectedEnergy(int partid, int method){

  if (method==2 && partid==2) return CorrectedEnergyElectrons;
  else if (method==2 && partid==1) return CorrectedEnergyPhotons;
  else {
    printf("Only result of method 2 and for partid 1 and 2 has been saved, returning 0.0...\n");
  }

  return 0.0;
}

//----------------------------------------------------------------------------------------

#ifdef _WITHGBATCH_

myEcalShowerFiller::myEcalShowerFiller(){
#ifdef PDEBUG
  printf("In myEcalShowerFiller::myEcalShowerFiller\n");
#endif
  PRINTDEBUG;
  init();
  PRINTDEBUG;
}

myEcalShowerFiller::~myEcalShowerFiller(){
#ifdef PDEBUG
  printf("In myEcalShowerFiller::~myEcalShowerFiller\n");
#endif
  PRINTDEBUG;
  Clear();
}

void myEcalShowerFiller::Clear(Option_t* option){
#ifdef PDEBUG
  printf("In myEcalShowerFiller::Clear\n");
#endif
  PRINTDEBUG;
  shower=0;
  PRINTDEBUG;
  return;
}

void myEcalShowerFiller::init(){
#ifdef PDEBUG
  printf("In myEcalShowerFiller::init\n");
#endif
  PRINTDEBUG;
  Clear();
  PRINTDEBUG;
  return;
}

void myEcalShowerFiller::Fill(EcalShowerR* _shower){
#ifdef PDEBUG
  printf("In myEventFiller::Fill\n");
#endif
  PRINTDEBUG;

  shower = _shower;

  if (shower) {
    /*
       memcpy(Dir, shower->Dir, 3*sizeof(Dir[0]));
       memcpy(CofG, shower->CofG, 3*sizeof(CofG[0]));
       memcpy(Entry, shower->Entry, 3*sizeof(Entry[0]));
       memcpy(Exit, shower->Exit, 3*sizeof(Exit[0]));
       */
    copy_n(shower->Dir, 3, Dir);
    copy_n(shower->CofG, 3, CofG);
    copy_n(shower->DirCR, 3, DirCR);
    copy_n(shower->CofGCR, 3, CofGCR);
    copy_n(shower->Entry, 3, Entry);
    copy_n(shower->Exit, 3, Exit);

    Status=shower->Status;
    Nhits=shower->Nhits;

    EnergyD=shower->EnergyD;
    EnergyC=shower->EnergyC;
    EnergyE=shower->EnergyE;//Pisa merda, W Livorno!
#if !defined _B524_ && !defined _B530_ && !defined _B538_ && !defined _B550_ && !defined _B572_ && !defined _B580_ && !defined _B584_ && !defined _B594_ && !defined _B598_ && !defined _B610_ && !defined _B620_ && !defined _B630_
    EnergyE2014=shower->EnergyE2014;//Pisa merda, W Livorno!
#endif //from B700
    EnergyA=shower->EnergyA;
    EnergyF=shower->EnergyF;

    DispLat = shower->ShowerLatDisp[0]/18;
    DispLatX = shower->ShowerLatDisp[1]/8;
    DispLatY = shower->ShowerLatDisp[2]/10;
    //    printf("DispLat: XY = %f, X = %f, Y = %f\n", DispLat, DispLatX, DispLatY);

    S35RaX=shower->S35Ra[0];
    S35RaY=shower->S35Ra[1];

    MaxS=shower->ParProfile[1];

    CorrectedEnergyElectrons = shower->GetCorrectedEnergy(2,2);
    CorrectedEnergyPhotons = shower->GetCorrectedEnergy(1,2);

  }

  PRINTDEBUG;

  return;
}

#endif //#ifdef _WITHGBATCH_