Add TOF-track related histograms

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.cxx b/PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.cxx
index d5c4a4ef987fa445a1784415c4d884ffc9da3fc5..c9d2e761975e4dc2ee3b20d2043199e0294f8767 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.cxx
+++ b/PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.cxx
@@ -39,22 +39,26 @@
 ClassImp(AliAnaChargedParticles)
   
 //__________________________________________________
-  AliAnaChargedParticles::AliAnaChargedParticles() : 
+  AliAnaChargedParticles::AliAnaChargedParticles() :
     AliAnaCaloTrackCorrBaseClass(),
     fPdg(0), 
-    fhNtracks(0),   fhPt(0),
-    fhPhiNeg(0),    fhEtaNeg(0), 
-    fhPhiPos(0),    fhEtaPos(0), 
-    fhEtaPhiPos(0), fhEtaPhiNeg(0),
+    fhNtracks(0),     fhPt(0),
+    fhPhiNeg(0),      fhEtaNeg(0), 
+    fhPhiPos(0),      fhEtaPos(0), 
+    fhEtaPhiPos(0),   fhEtaPhiNeg(0),
     //MC
-    fhPtPion(0),    fhPhiPion(0),     fhEtaPion(0),
-    fhPtProton(0),  fhPhiProton(0),   fhEtaProton(0),
-    fhPtElectron(0),fhPhiElectron(0), fhEtaElectron(0),
-    fhPtKaon(0),    fhPhiKaon(0),     fhEtaKaon(0),
-    fhPtUnknown(0), fhPhiUnknown(0),  fhEtaUnknown(0)
+    fhPtPion(0),      fhPhiPion(0),      fhEtaPion(0),
+    fhPtProton(0),    fhPhiProton(0),    fhEtaProton(0),
+    fhPtElectron(0),  fhPhiElectron(0),  fhEtaElectron(0),
+    fhPtKaon(0),      fhPhiKaon(0),      fhEtaKaon(0),
+    fhPtUnknown(0),   fhPhiUnknown(0),   fhEtaUnknown(0),
+    fhTOFSignal(0),   fhTOFSignalPtCut(0),
+    fhPtTOFSignal(0), fhPtTOFStatus0(0), fhEtaPhiTOFStatus0(0)
 {
   //Default Ctor
 
+  for(Int_t i = 0; i < 7; i++) fhPtTOFSignalPileUp[i] = 0;
+  
   //Initialize parameters
   InitParameters();
 
@@ -106,7 +110,7 @@ TList *  AliAnaChargedParticles::GetCreateOutputObjects()
   fhEtaPos->SetXTitle("p_{T} (GeV/c)");
   outputContainer->Add(fhEtaPos);
   
-  fhEtaPhiPos  = new TH2F ("hEtaPhiPositive","pt/eta/phi of positive charge",netabins,etamin,etamax, nphibins,phimin,phimax); 
+  fhEtaPhiPos  = new TH2F ("hEtaPhiPositive","pt/eta/phi of positive charge",netabins,etamin,etamax, nphibins,phimin,phimax);
   fhEtaPhiPos->SetXTitle("#eta ");
   fhEtaPhiPos->SetYTitle("#phi (rad)");  
   outputContainer->Add(fhEtaPhiPos);
@@ -116,6 +120,46 @@ TList *  AliAnaChargedParticles::GetCreateOutputObjects()
   fhEtaPhiNeg->SetYTitle("#phi (rad)");  
   outputContainer->Add(fhEtaPhiNeg);
   
+  Int_t ntofbins = 1000;
+  Int_t mintof = -500;
+  Int_t maxtof =  500;
+
+  fhTOFSignal  = new TH1F ("hTOFSignal","TOF signal", ntofbins,mintof,maxtof);
+  fhTOFSignal->SetXTitle("TOF signal #times (ns?)");
+  outputContainer->Add(fhTOFSignal);
+
+  fhTOFSignalPtCut  = new TH1F ("hTOFSignalPtCut","TOF signal", ntofbins,mintof,maxtof);
+  fhTOFSignalPtCut->SetXTitle("TOF signal (ns?)");
+  outputContainer->Add(fhTOFSignalPtCut);
+
+  fhPtTOFSignal  = new TH2F ("hPtTOFSignal","TOF signal", nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
+  fhPtTOFSignal->SetYTitle("TOF signal (ns?)");
+  fhPtTOFSignal->SetXTitle("p_{T} (GeV/c)");
+  outputContainer->Add(fhPtTOFSignal);
+  
+  TString pileUpName[] = {"SPD","EMCAL","SPDOrEMCAL","SPDAndEMCAL","SPDAndNotEMCAL","EMCALAndNotSPD","NotSPDAndNotEMCAL"} ;
+  
+  for(Int_t i = 0 ; i < 7 ; i++)
+  {
+    fhPtTOFSignalPileUp[i]  = new TH2F(Form("hPtTOFSignalPileUp%s",pileUpName[i].Data()),
+                                       Form("Track TOF vs p_{T} distribution, %s Pile-Up event",pileUpName[i].Data()),
+                                       nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
+    fhPtTOFSignalPileUp[i]->SetXTitle("p_{T} (GeV/c)");
+    fhPtTOFSignalPileUp[i]->SetXTitle("TOF signal (ns?)");
+    outputContainer->Add(fhPtTOFSignalPileUp[i]);
+  }
+  
+  fhPtTOFStatus0  = new TH1F ("hPtTOFStatus0","p_T distribution of tracks not hitting TOF", nptbins,ptmin,ptmax);
+  fhPtTOFStatus0->SetXTitle("p_{T} (GeV/c)");
+  outputContainer->Add(fhPtTOFStatus0);
+  
+  
+  fhEtaPhiTOFStatus0  = new TH2F ("hEtaPhiTOFStatus0","pt/eta/phi for tracks without hit on TOF",netabins,etamin,etamax, nphibins,phimin,phimax);
+  fhEtaPhiTOFStatus0->SetXTitle("#eta ");
+  fhEtaPhiTOFStatus0->SetYTitle("#phi (rad)");
+  outputContainer->Add(fhEtaPhiTOFStatus0);
+
+  
   if(IsDataMC()){
     
     fhPtPion  = new TH1F ("hPtMCPion","p_T distribution from #pi", nptbins,ptmin,ptmax); 
@@ -245,10 +289,30 @@ void  AliAnaChargedParticles::MakeAnalysisFillAOD()
     
     AliVTrack * track =  (AliVTrack*) (GetCTSTracks()->At(i));
     
-    //Fill AODParticle after some selection       
+     //Fill AODParticle after some selection       
     Double_t mom[3] = {track->Px(),track->Py(),track->Pz()};
     p3.SetXYZ(mom[0],mom[1],mom[2]);
     
+    //TOF
+    ULong_t status = track->GetStatus();
+    Bool_t okTOF = ( (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ) && ( (status & AliVTrack::kTIME) == AliVTrack::kTIME );
+    Double32_t tof = track->GetTOFsignal()*1e-3;
+    //if( tof < 0) printf("TOF Signal %e, status %d, pt %f\n", tof,status,status2,p3.Pt());
+    
+    if(okTOF)
+    {
+      fhTOFSignal  ->Fill(tof);
+      fhPtTOFSignal->Fill(p3.Pt(), tof);
+            
+      if(GetReader()->IsPileUpFromSPD())               fhPtTOFSignalPileUp[0]->Fill(p3.Pt(), tof); 
+      if(GetReader()->IsPileUpFromEMCal())             fhPtTOFSignalPileUp[1]->Fill(p3.Pt(), tof);
+      if(GetReader()->IsPileUpFromSPDOrEMCal())        fhPtTOFSignalPileUp[2]->Fill(p3.Pt(), tof);
+      if(GetReader()->IsPileUpFromSPDAndEMCal())       fhPtTOFSignalPileUp[3]->Fill(p3.Pt(), tof);
+      if(GetReader()->IsPileUpFromSPDAndNotEMCal())    fhPtTOFSignalPileUp[4]->Fill(p3.Pt(), tof);
+      if(GetReader()->IsPileUpFromEMCalAndNotSPD())    fhPtTOFSignalPileUp[5]->Fill(p3.Pt(), tof);
+      if(GetReader()->IsPileUpFromNotSPDAndNotEMCal()) fhPtTOFSignalPileUp[6]->Fill(p3.Pt(), tof);
+    }
+    
     Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,"CTS") ;
     
     if(GetDebug() > 1) 
@@ -258,16 +322,25 @@ void  AliAnaChargedParticles::MakeAnalysisFillAOD()
     if(IsFiducialCutOn() && ! in ) continue ;
     
     // Momentum selection
-    if(p3.Pt() < GetMinPt() || p3.Pt() > GetMaxPt()) continue;
+    if(track->Pt() < GetMinPt() || track->Pt() > GetMaxPt()) continue;
+    
+    if(okTOF) fhTOFSignalPtCut->Fill(tof); 
+    else
+    {
+      fhPtTOFStatus0    ->Fill(track->Pt());
+      fhEtaPhiTOFStatus0->Fill(track->Eta(),track->Phi());
+    }
     
     //Keep only particles identified with fPdg
     //Selection not done for the moment
     //Should be done here.
     
     // Mixed event
-    if (GetMixedEvent()){
+    if (GetMixedEvent())
+    {
       evtIndex = GetMixedEvent()->EventIndex(track->GetID()) ;
-    } 
+    }
+    
     GetVertex(vert,evtIndex); 
     if(TMath::Abs(vert[2])> GetZvertexCut()) return; 
         

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.h b/PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.h
index 9e50decd1e5b93aa34769460a903c50a6fd1d403..ec82480a75e131bb74ac796ef5457c61426ab94a 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.h
+++ b/PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.h
@@ -50,8 +50,8 @@ class AliAnaChargedParticles : public AliAnaCaloTrackCorrBaseClass {
   TH2F * fhPhiPos;      //! phi distribution vs pT, positive
   TH2F * fhEtaPos;      //! eta distribution vs pT, positive
   TH2F * fhEtaPhiPos;   //! eta vs phi distribution of positive charge  
-  TH2F * fhEtaPhiNeg;   //! eta vs phi distribution of negative charge  
-
+  TH2F * fhEtaPhiNeg;   //! eta vs phi distribution of negative charge
+  
   //MC
   TH1F * fhPtPion;      //! pT distribution
   TH2F * fhPhiPion;     //! phi distribution vs pT
@@ -73,10 +73,18 @@ class AliAnaChargedParticles : public AliAnaCaloTrackCorrBaseClass {
   TH2F * fhPhiUnknown;  //! phi distribution vs pT
   TH2F * fhEtaUnknown;  //! eta distribution vs pT
   
+  // TOF
+  TH1F * fhTOFSignal;        //! TOF signal, good status
+  TH1F * fhTOFSignalPtCut;   //! TOF signal, good status, pt and acceptance cut
+  TH2F * fhPtTOFSignal;      //! TOF signal vs track pT, good status
+  TH2F * fhPtTOFSignalPileUp[7]; //! TOF signal vs track pT, good status, pile-up
+  TH1F * fhPtTOFStatus0;     //! pT of tracks not passing TOF status selection
+  TH2F * fhEtaPhiTOFStatus0; //! eta/phi of tracks not passing TOF status selection
+
   AliAnaChargedParticles(              const AliAnaChargedParticles & ch) ; // cpy ctor
   AliAnaChargedParticles & operator = (const AliAnaChargedParticles & ch) ; // cpy assignment
   
-  ClassDef(AliAnaChargedParticles,3)
+  ClassDef(AliAnaChargedParticles,4)
 
 } ;