am->LoadBranch("PrimaryVertex.");
*priv = *esdin->GetPrimaryVertexTracks();
}
-#if 1
AliESDCaloCells *ecells = dynamic_cast<AliESDCaloCells*>(objsout->FindObject("EMCALCells"));
if (ecells) {
am->LoadBranch("EMCALCells.");
*ecells = *esdin->GetEMCALCells();
}
-#else
- TClonesArray *ecells = dynamic_cast<TClonesArray*>(objsout->FindObject("EMCALCells"));
- if (ecells) {
- am->LoadBranch("EMCALCells.");
- AliESDCaloCells *ein = esdin->GetEMCALCells();
- ecells->Clear();
- for(Int_t i=0;i<ein->GetNumberOfCells();++i) {
- Short_t cellNumber = 0;
- Double_t amplitude = 0;
- Double_t time = 0;
- ein->GetCell(i, cellNumber, amplitude, time);
- new ((*ecells)[i]) AliCaloCell(cellNumber, amplitude, time);
- }
- }
-#endif
-
AliESDCaloCells *pcells = dynamic_cast<AliESDCaloCells*>(objsout->FindObject("PHOSCells"));
if (pcells) {
am->LoadBranch("PHOSCells.");
}
if (fDoEmCs) {
fEvent->AddObject(new AliESDCaloCells("EMCALCells","EMCALCells"));
- //TClonesArray *arr = new TClonesArray("AliCaloCell",0);
- //arr->SetName("EMCALCells");
- //fEvent->AddObject(arr);
}
if (fDoPCs) {
fEvent->AddObject(new AliESDCaloCells("PHOSCells","PHOSCells"));
class AliESDEvent;
class AliESDtrackCuts;
-class AliCaloCell : public TObject
-{
- public:
- AliCaloCell(Short_t num=0, Double_t amp=0., Double_t time=0) :
- TObject(), fNum(num), fAmp(amp), fTime(time*1e9), fMcFr(0) {;}
-
- Bool_t IsSortable() const { return kTRUE; }
- //Int_t Compare(const TObject* obj) const;
-
- public:
- Short_t fNum; // cell number/id
- Double32_t fAmp; //[0,0,16] cell amplitude
- Double32_t fTime; //[0,0,16] cell time
- Double32_t fMcFr; //[0,1,10] fraction of MC energy (fMcAmp=fMcFr*fAmp)
-
- ClassDef(AliCaloCell,1) // Calo cell class
-};
-
class AliEsdSkimTask : public AliAnalysisTaskSE {
public:
AliEsdSkimTask(const char *opt=0);