#include "AliJetESDReaderHeader.h"
#include "AliESD.h"
#include "AliESDtrack.h"
-#include "AliEMCALGeometry.h"
+//#include "AliEMCALGeometry.h"
+#include "AliJetDummyGeo.h"
#include "AliJetFillUnitArrayTracks.h"
#include "AliJetUnitArray.h"
fArrayInitialised(0)
{
// Constructor
- SetEMCALGeometry();
}
//____________________________________________________________________________
int nMax = fChain->GetEntries();
- printf("\nTotal number of events in chain= %d",nMax);
+ printf("\nTotal number of events in chain= %d \n",nMax);
// set number of events in header
if (fReaderHeader->GetLastEvent() == -1)
fChain->SetBranchAddress("ESD", &fESD);
Int_t nMax = fChain->GetEntries();
- printf("\nTotal number of events in chain= %5d", nMax);
+ printf("\nTotal number of events in chain= %5d \n", nMax);
// set number of events in header
if (fReaderHeader->GetLastEvent() == -1)
fReaderHeader->SetLastEvent(nMax);
TVector3 p3;
// clear momentum array
- ClearArray();
- fDebug = fReaderHeader->GetDebug();
- InitParameters();
+ ClearArray();
+ fDebug = fReaderHeader->GetDebug();
// get event from chain
fChain->GetTree()->GetEntry(event);
-
+
if (!fESD) {
return kFALSE;
}
sflag[goodTrack]=0;
if (TMath::Abs(track->GetLabel()) < 10000) sflag[goodTrack]=1;
cflag[goodTrack]=0;
- if (pt > ptMin) cflag[goodTrack]=1; // pt cut
+ if (pt > ptMin) cflag[goodTrack]=1; // pt cut
goodTrack++;
}
// set the signal flags
//
//
- AliJetFillUnitArrayTracks *fillUAFromTracks = new AliJetFillUnitArrayTracks();
- fillUAFromTracks->SetReaderHeader(fReaderHeader);
- fillUAFromTracks->SetMomentumArray(fMomentumArray);
- fillUAFromTracks->SetTPCGrid(fTpcGrid);
- fillUAFromTracks->SetEMCalGrid(fEmcalGrid);
- fillUAFromTracks->SetHadCorrection(fHCorrection);
- fillUAFromTracks->SetHadCorrector(fHadCorr);
- fNeta = fillUAFromTracks->GetNeta();
- fNphi = fillUAFromTracks->GetNphi();
- fillUAFromTracks->SetActive(kFALSE);
- // TPC only or Digits+TPC or Clusters+TPC
- if(fOpt%2==!0 && fOpt!=0){
- fillUAFromTracks->SetActive(kTRUE);
- fillUAFromTracks->SetUnitArray(fUnitArray);
- fillUAFromTracks->ExecuteTask("tpc");
+ if (fTpcGrid || fEmcalGrid) {
+ SetEMCALGeometry();
+ InitParameters();
+ AliJetFillUnitArrayTracks *fillUAFromTracks = new AliJetFillUnitArrayTracks();
+ fillUAFromTracks->SetReaderHeader(fReaderHeader);
+ fillUAFromTracks->SetMomentumArray(fMomentumArray);
+ fillUAFromTracks->SetTPCGrid(fTpcGrid);
+ fillUAFromTracks->SetEMCalGrid(fEmcalGrid);
+ fillUAFromTracks->SetHadCorrection(fHCorrection);
+ fillUAFromTracks->SetHadCorrector(fHadCorr);
+ fNeta = fillUAFromTracks->GetNeta();
+ fNphi = fillUAFromTracks->GetNphi();
+ fillUAFromTracks->SetActive(kFALSE);
+ // TPC only or Digits+TPC or Clusters+TPC
+ if(fOpt%2==!0 && fOpt!=0) {
+ fillUAFromTracks->SetActive(kTRUE);
+ fillUAFromTracks->SetUnitArray(fUnitArray);
+ fillUAFromTracks->ExecuteTask("tpc");
+ }
+
+ delete fillUAFromTracks;
}
- delete fillUAFromTracks;
+
return kTRUE;
}
void AliJetESDReader::SetEMCALGeometry()
{
// Define EMCAL geometry to be able to read ESDs
- fGeom = AliEMCALGeometry::GetInstance();
+ fGeom = AliJetDummyGeo::GetInstance();
if (fGeom == 0)
- fGeom = AliEMCALGeometry::GetInstance("SHISH_77_TRD1_2X2_FINAL_110DEG","EMCAL");
+ fGeom = AliJetDummyGeo::GetInstance("SHISH_77_TRD1_2X2_FINAL_110DEG","EMCAL");
// To be setted to run some AliEMCALGeometry functions
TGeoManager::Import("geometry.root");
fGeom->GetTransformationForSM();
-
printf("\n EMCal Geometry set ! \n");
}