\r
//____________________________________________________________________\r
void AliTRDcheckESD::MakeSummary(){\r
+ //\r
+ // Draw summary plots for the ESDcheck task\r
+ //\r
TCanvas *cOut = new TCanvas(Form("summary%s1", GetName()), Form("Summary 1 for task %s", GetName()), 1024, 768);\r
cOut->cd();\r
GetRefFigure(4);\r
//____________________________________________________________________\r
Bool_t AliTRDcheckESD::GetRefFigure(Int_t ifig)\r
{\r
+ //\r
+ // Produce reference Plots during PostProcessing\r
+ //\r
if(ifig>=fNRefFigures){\r
AliWarning(Form("Ref plot %d not available. Valid only up to %d", ifig, fNRefFigures));\r
return kFALSE;\r
h = (TH2F*)fHistos->At(kNClsTrackTRD); h->Fill(esdTrack->GetP(), esdTrack->GetTRDncls());\r
// (slicePH,sliceNo) distribution and Qtot from slices\r
for(Int_t iPlane=0; iPlane<6; iPlane++) {\r
- Float_t Qtot=0;\r
+ Float_t qtot=0;\r
for(Int_t iSlice=0; iSlice<8; iSlice++) {\r
if(esdTrack->GetTRDslice(iPlane, iSlice)>20.) {\r
h = (TH2F*)fHistos->At(kPHSlice); h->Fill(iSlice, esdTrack->GetTRDslice(iPlane, iSlice));\r
- Qtot += esdTrack->GetTRDslice(iPlane, iSlice);\r
+ qtot += esdTrack->GetTRDslice(iPlane, iSlice);\r
}\r
}\r
// Qtot>100 to avoid noise\r
- if(Qtot>100.) {\r
- h = (TH2F*)fHistos->At(kQtotP); h->Fill(esdTrack->GetTRDmomentum(iPlane), fgkQs*Qtot);\r
+ if(qtot>100.) {\r
+ h = (TH2F*)fHistos->At(kQtotP); h->Fill(esdTrack->GetTRDmomentum(iPlane), fgkQs*qtot);\r
}\r
// Qtot>100 to avoid noise\r
// fgkQs*Qtot<40. so that the average will give a value close to the peak\r
- if(localCoordGood && localMomGood && Qtot>100. && fgkQs*Qtot<40.) {\r
+ if(localCoordGood && localMomGood && qtot>100. && fgkQs*qtot<40.) {\r
h2d = (TProfile2D*)fHistos->At(kTRDEtaPhiAvQtot+iPlane);\r
- h2d->Fill(eta, localSagitaPhi, fgkQs*Qtot);\r
+ h2d->Fill(eta, localSagitaPhi, fgkQs*qtot);\r
}\r
}\r
// theta distribution\r
\r
// clusters per track\r
const Int_t kNpt(30);\r
- Float_t Pt(0.2);\r
+ Float_t pt(0.2);\r
Float_t binsPt[kNpt+1];\r
- for(Int_t i=0;i<kNpt+1; i++,Pt+=(TMath::Exp(i*i*.001)-1.)) binsPt[i]=Pt;\r
+ for(Int_t i=0;i<kNpt+1; i++,pt+=(TMath::Exp(i*i*.001)-1.)) binsPt[i]=pt;\r
if(!(h = (TH2S*)gROOT->FindObject("hNCl"))){\r
h = new TH2S("hNCl", "Clusters per TRD track;N_{cl}^{TRD};SPECIES;entries", 60, 0., 180., 10, -0.5, 9.5);\r
TAxis *ay(h->GetYaxis());\r
\r
// status bits histogram\r
const Int_t kNbits(5);\r
- Float_t Bits(.5);\r
+ Float_t bits(.5);\r
Float_t binsBits[kNbits+1];\r
- for(Int_t i=0; i<kNbits+1; i++,Bits+=1.) binsBits[i]=Bits;\r
+ for(Int_t i=0; i<kNbits+1; i++,bits+=1.) binsBits[i]=bits;\r
if(!(h = (TH2I*)gROOT->FindObject("hTRDstat"))){\r
h = new TH2I("hTRDstat", "TRD status bits;p_{t} @ TRD [GeV/c];status;entries", kNpt, binsPt, kNbits, binsBits);\r
TAxis *ay(h->GetYaxis());\r
\r
// pt resolution\r
const Int_t kNdpt(100), kNspec(4*AliPID::kSPECIES);\r
- Float_t DPt(-3.), Spec(-0.5);\r
+ Float_t dpt(-3.), spec(-0.5);\r
Float_t binsDPt[kNdpt+1], binsSpec[kNspec+1];\r
- for(Int_t i=0; i<kNdpt+1; i++,DPt+=6.e-2) binsDPt[i]=DPt;\r
- for(Int_t i=0; i<kNspec+1; i++,Spec+=1.) binsSpec[i]=Spec;\r
+ for(Int_t i=0; i<kNdpt+1; i++,dpt+=6.e-2) binsDPt[i]=dpt;\r
+ for(Int_t i=0; i<kNspec+1; i++,spec+=1.) binsSpec[i]=spec;\r
if(!(h = (TH3S*)gROOT->FindObject("hPtRes"))){\r
h = new TH3S("hPtRes", "P_{t} resolution @ DCA;p_{t}^{MC} [GeV/c];#Delta p_{t}/p_{t}^{MC} [%];SPECIES", kNpt, binsPt, kNdpt, binsDPt, kNspec, binsSpec);\r
TAxis *az(h->GetZaxis());\r
}\r
\r
//____________________________________________________________________\r
-Int_t AliTRDcheckESD::Pdg2Idx(Int_t pdg)\r
+Int_t AliTRDcheckESD::Pdg2Idx(Int_t pdg) const\r
{\r
+ //\r
+ // Helper function converting PDG code into AliPID index\r
+ //\r
switch(pdg){\r
case kElectron: \r
case kPositron: return AliPID::kElectron; \r