changes from soh
authormcosenti <mcosenti@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Jul 2012 14:24:21 +0000 (14:24 +0000)
committermcosenti <mcosenti@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Jul 2012 14:24:21 +0000 (14:24 +0000)
PWGGA/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.cxx
PWGGA/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.h

index c1dc525..98c92a9 100644 (file)
@@ -39,7 +39,8 @@ AliAnalysisTaskSOH::AliAnalysisTaskSOH() :
   fHEventStat(0),      
   fHTrkEffParGenPt(0), 
   fHTrkEffDetGenPt(0), 
-  fHTrkEffDetRecPt(0), 
+  fHTrkEffDetRecPt(0),
+  fHTrkEffDetRecFakePt(0), 
   fHScaleFactor(0),
   fHScaleFactor100HC(0),
   fHEOverPVsPt(0x0),
@@ -79,6 +80,7 @@ AliAnalysisTaskSOH::AliAnalysisTaskSOH(const char *name) :
   fHTrkEffParGenPt(0), 
   fHTrkEffDetGenPt(0), 
   fHTrkEffDetRecPt(0), 
+  fHTrkEffDetRecFakePt(0),
   fHScaleFactor(0),
   fHScaleFactor100HC(0),
   fHEOverPVsPt(0x0),
@@ -148,6 +150,9 @@ void AliAnalysisTaskSOH::UserCreateOutputObjects()
   fHTrkEffDetRecPt = new TH1F("fHTrkEffDetRecPt", "Reconstructed track p_{T} distribution of primary charged pions;p_{T}^{rec} (GeV/c)",500,0.0,50.0);
   fOutputList->Add(fHTrkEffDetRecPt);
 
+  fHTrkEffDetRecFakePt = new TH1F("fHTrkEffDetRecFakePt", "Reconstructed fake track p_{T} distribution of pions;p_{T}^{rec} (GeV/c)",500,0.0,50.0);
+  fOutputList->Add(fHTrkEffDetRecFakePt);
+
   fHScaleFactor = new TH1F("fHScaleFactor", "Scale factor distribution without hadronic correction;Scale factor",100,0,10);
   fOutputList->Add(fHScaleFactor);
 
@@ -333,6 +338,11 @@ void  AliAnalysisTaskSOH::ProcessTrack()
         }
       }
     }
+
+ // fake and secondary tracks
+    if(newTrack->GetLabel()<0 && newTrack->GetPID()==2) fHTrkEffDetRecFakePt->Fill(newTrack->Pt());
+
+ // Track Indices
     fTrackIndices->AddAt(itr,nTracks);
     nTracks++;
   }
index b9e8924..725fc2a 100644 (file)
@@ -53,6 +53,7 @@ class AliAnalysisTaskSOH : public AliAnalysisTaskSE {
   TH1F               *fHTrkEffParGenPt;          //!mc truth pt spectrum
   TH1F               *fHTrkEffDetGenPt;          //!mc detector level pt spectrum
   TH1F               *fHTrkEffDetRecPt;          //!reconstructed detector level pt spectrum
+  TH1F               *fHTrkEffDetRecFakePt;      //!fake and secondary tracks pt spectrum
   TH1F               *fHScaleFactor;             //!scale factor spectrum
   TH1F               *fHScaleFactor100HC;        //!scale factor with 100% HC spectrum
   TH2F               *fHEOverPVsPt;              //!(cluster energy over reconstructed track p) vs. track pt