fTrackPtEta(0),
fTrackPtEtaCut(0),
fMaxCellEPhi(0),
- fDetaDphiFromTM(0)
+ fDetaDphiFromTM(0),
+ fEoverPvsE(0)
{
// Default constructor.
for(Int_t i = 0; i < 12; i++) fGeomMatrix[i] = 0;
fTrackPtEta(0),
fTrackPtEtaCut(0),
fMaxCellEPhi(0),
- fDetaDphiFromTM(0)
+ fDetaDphiFromTM(0),
+ fEoverPvsE(0)
{
// Constructor
fDetaDphiFromTM->Sumw2();
fQAList->Add(fDetaDphiFromTM);
+ fEoverPvsE = new TH2F("fEoverPvsE","E^{clus}/p^{track} vs E^{clus} (80<TPCsignal<100);E^{clus} [GeV];E^{clus}/p^{track} [c^{-1}]",fNBinsPt, fPtBinLowEdge,fPtBinHighEdge,100,0,2);
+ fEoverPvsE->Sumw2();
+ fQAList->Add(fEoverPvsE);
+
PostData(1, fOutputList);
PostData(2, fQAList);
}
fTrackPtPhiCut->Fill(t->Pt(), t->Phi());
fTrackPtEtaCut->Fill(t->Pt(), t->Eta());
}
+ if(t->GetTPCsignal()<80 || t->GetTPCsignal()>100)
+ continue;
+ if(t->GetEMCALcluster()<=0 || t->GetEMCALcluster()>nclus)
+ continue;
+ AliVCluster *c = dynamic_cast<AliVCluster*>(clusters->At(t->GetEMCALcluster()));
+ if(!c)
+ continue;
+ fEoverPvsE->Fill(c->E(),c->E()/t->P());
}
for(int ic=0;ic<nclus;ic++){
AliVCluster *c = dynamic_cast<AliVCluster*>(clusters->At(ic));
TH2F *fTrackPtEtaCut; //!selected tracks pt vs. eta in "triggered event"
TH2F *fMaxCellEPhi; //!max cell energy vs. cell phi
TH2F *fDetaDphiFromTM; //!dphi vs deta of track->GetEMCALcluster() clusters
+ TH2F *fEoverPvsE; //!E/p for tracks with 80<TPCsignal<100 vs cluster E (check material)
AliAnalysisTaskEMCALIsoPhoton(const AliAnalysisTaskEMCALIsoPhoton&); // not implemented