updates to QA task, varibale pt binning, catch the case cent == 0 in spectrum task...

[u/mrichter/AliRoot.git] / PWG4 / JetTasks / AliPWG4HighPtQAMC.cxx

diff --git a/PWG4/JetTasks/AliPWG4HighPtQAMC.cxx b/PWG4/JetTasks/AliPWG4HighPtQAMC.cxx
index cccbb028d0942020bc646cacd9ea449cfe0284c6..84978cd3f14ab9700beeac6ace356dc056e6a60b 100644 (file)
--- a/PWG4/JetTasks/AliPWG4HighPtQAMC.cxx
+++ b/PWG4/JetTasks/AliPWG4HighPtQAMC.cxx
@@ -64,6 +64,7 @@ AliPWG4HighPtQAMC::AliPWG4HighPtQAMC()
   fTrackCuts(0), 
   fTrackCutsITS(0),
   fTrackType(0),
+  fSigmaConstrainedMax(5.),
   fPtMax(100.),
   fAvgTrials(1),
   fNEventAll(0),
@@ -121,6 +122,7 @@ AliPWG4HighPtQAMC::AliPWG4HighPtQAMC(const char *name):
   fTrackCuts(),
   fTrackCutsITS(),
   fTrackType(0),
+  fSigmaConstrainedMax(5.),
   fPtMax(100.),
   fAvgTrials(1),
   fNEventAll(0),
@@ -604,6 +606,14 @@ void AliPWG4HighPtQAMC::Exec(Option_t *) {
       continue;
     }
 
+    if(fTrackType==2) {
+      //Cut on chi2 of constrained fit
+      if(track->GetConstrainedChi2TPC() > fSigmaConstrainedMax*fSigmaConstrainedMax) {
+       delete track;
+       continue;
+      }
+    }
+
     Int_t label = TMath::Abs(track->GetLabel());
     if(label>=nMCtracks) {
       if (fTrackCuts->AcceptTrack(track)) { fPtSelLargeLabel->Fill(pt); }