fTrackPtPhi(0),
fTrackPtPhiCut(0),
fTrackPtEta(0),
- fTrackPtEtaCut(0)
+ fTrackPtEtaCut(0),
+ fMaxCellEPhi(0)
{
// Default constructor.
for(Int_t i = 0; i < 12; i++) fGeomMatrix[i] = 0;
fTrackPtPhi(0),
fTrackPtPhiCut(0),
fTrackPtEta(0),
- fTrackPtEtaCut(0)
+ fTrackPtEtaCut(0),
+ fMaxCellEPhi(0)
{
// Constructor
fTrackPtEtaCut->Sumw2();
fQAList->Add(fTrackPtEtaCut);
+
+ fMaxCellEPhi = new TH2F("fMaxCellEPhi","Most energetic cell in event; GeV;#phi",100,-0.25,49.75,63,0,6.3);
+ fMaxCellEPhi->Sumw2();
+ fQAList->Add(fMaxCellEPhi);
+
PostData(1, fOutputList);
PostData(2, fQAList);
}
}
}
for(int ic=0;ic<nclus;ic++){
- AliVCluster *c = (AliVCluster*)fESDClusters->At(ic);
+ AliVCluster *c = dynamic_cast<AliVCluster*>(fESDClusters->At(ic));
+ //AliESDCaloCluster *c = (AliESDCaloCluster*)fESDClusters->At(ic);
if(!c)
continue;
+ if(!c->IsEMCAL())
+ continue;
Float_t clsPos[3] = {0,0,0};
c->GetPosition(clsPos);
TVector3 clsVec(clsPos);
cells = fAODCells;
if(!cells)
return;
+ Double_t maxe = 0;
+ Double_t maxphi = -10;
Int_t ncells = cells->GetNumberOfCells();
+ Double_t eta,phi;
for (Int_t i=0; i<ncells; i++) {
Short_t absid = TMath::Abs(cells->GetCellNumber(i));
Double_t e = cells->GetCellAmplitude(absid);
if(e>0.05)
fNCells50++;
+ else
+ continue;
+ fGeom->EtaPhiFromIndex(absid,eta,phi);
+ if(maxe<e){
+ maxe = e;
+ maxphi = phi;
+ }
}
+ fMaxCellEPhi->Fill(maxe,maxphi);
}
//________________________________________________________________________
TH2F *fTrackPtPhiCut; //!selected tracks pt vs. phi in "triggered event"
TH2F *fTrackPtEta; //!selected tracks pt vs. eta
TH2F *fTrackPtEtaCut; //!selected tracks pt vs. eta in "triggered event"
+ TH2F *fMaxCellEPhi; //!max cell energy vs. cell phi
AliAnalysisTaskEMCALIsoPhoton(const AliAnalysisTaskEMCALIsoPhoton&); // not implemented