modified dNdPt/AlidNdPtAnalysisPbPbAOD.{h,cxx}

- added crosscheck histos for cross rows
- bug fix

diff --git a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.cxx b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.cxx
index 267e7a9027e766eb8ed71c16bf7a347c7bbb9384..2c585a1109193820c4e45d54a5418d4031c804eb 100644 (file)
--- a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.cxx
+++ b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.cxx
@@ -49,6 +49,8 @@ hMCCharge(0),
 hMCPdgPt(0),
 hMCHijingPrim(0),
 hAccNclsTPC(0),
+hAllCrossedRowsTPC(0),
+hFilterCrossedRowsTPC(0),
 hAccCrossedRowsTPC(0),
 hDCAPtAll(0),
 hDCAPtAccepted(0),
@@ -133,6 +135,8 @@ hMCCharge(0),
 hMCPdgPt(0),
 hMCHijingPrim(0),
 hAccNclsTPC(0),
+hAllCrossedRowsTPC(0),
+hFilterCrossedRowsTPC(0),
 hAccCrossedRowsTPC(0),
 hDCAPtAll(0),
 hDCAPtAccepted(0),
@@ -356,6 +360,12 @@ void AlidNdPtAnalysisPbPbAOD::UserCreateOutputObjects()
   hAccNclsTPC = new TH1F("hAccNclsTPC","hAccNclsTPC",160,0,159);
   hAccNclsTPC->GetXaxis()->SetTitle("number of clusters per track after cut");
   
+  hAllCrossedRowsTPC = new TH1F("hAllCrossedRowsTPC","hAllCrossedRowsTPC",160,0,159);
+  hAllCrossedRowsTPC->GetXaxis()->SetTitle("number of crossed rows per track for all tracks");
+  
+  hFilterCrossedRowsTPC = new TH1F("hFilterCrossedRowsTPC","hFilterCrossedRowsTPC",160,0,159);
+  hFilterCrossedRowsTPC->GetXaxis()->SetTitle("number of crossed rows per track after filter bit");
+  
   hAccCrossedRowsTPC = new TH1F("hAccCrossedRowsTPC","hAccCrossedRowsTPC",160,0,159);
   hAccCrossedRowsTPC->GetXaxis()->SetTitle("number of crossed rows per track after cut");
   
@@ -414,6 +424,8 @@ void AlidNdPtAnalysisPbPbAOD::UserCreateOutputObjects()
   fOutputList->Add(hMCPdgPt);
   fOutputList->Add(hMCHijingPrim);
   fOutputList->Add(hAccNclsTPC);
+  fOutputList->Add(hAllCrossedRowsTPC);
+  fOutputList->Add(hFilterCrossedRowsTPC);
   fOutputList->Add(hAccCrossedRowsTPC);
   fOutputList->Add(hDCAPtAll);
   fOutputList->Add(hDCAPtAccepted);
@@ -714,11 +726,15 @@ Bool_t AlidNdPtAnalysisPbPbAOD::IsTrackAccepted(AliAODTrack *tr)
   
   if(tr->Charge()==0) { return kFALSE; }
   
-  if(!(tr->TestFilterBit(AliAODTrack::kTrkGlobal)) ) { return kFALSE; }
-  
   Double_t dNClustersTPC = tr->GetTPCNcls();
   Double_t dCrossedRowsTPC = tr->GetTPCClusterInfo(2,1);
   
+  hAllCrossedRowsTPC->Fill(dCrossedRowsTPC);
+  
+  if(!(tr->TestFilterBit(AliAODTrack::kTrkGlobal)) ) { return kFALSE; }
+  
+  hFilterCrossedRowsTPC->Fill(dCrossedRowsTPC);
+  
   if(dCrossedRowsTPC < GetCutMinNCrossedRowsTPC()) { return kFALSE; }
   
   hAccNclsTPC->Fill(dNClustersTPC);

diff --git a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.h b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.h
index 2061921ab15f116180e2b79fb7bf4b5348cae42e..b506f41f3cd98e58be74880cb86c09414b832894 100644 (file)
--- a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.h
+++ b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.h
@@ -145,6 +145,8 @@ class AlidNdPtAnalysisPbPbAOD : public AliAnalysisTaskSE {
     TH2F       *hMCPdgPt; // PDGvs PT for MC Particles
     TH1F       *hMCHijingPrim; // number of particles, which are Hijing particles and primaries
     TH1F       *hAccNclsTPC; //control histo: number of clusters in TPC for accepted tracks
+    TH1F       *hAllCrossedRowsTPC; //control histo: number of crossed rows in TPC for all tracks, before Filterbit
+    TH1F       *hFilterCrossedRowsTPC; //control histo: number of crossed rows in TPC for all tracks, after Filterbit
     TH1F       *hAccCrossedRowsTPC; //control histo: number of crossed rows in TPC for accepted tracks
     THnSparseF *hDCAPtAll; //control histo: DCAz vs DCAxy vs pT for all reconstructed tracks
     THnSparseF *hDCAPtAccepted; //control histo: DCAz vs DCAxy vs pT for all accepted reco tracks