cut on Dphi for tracklets
authorlmilano <lmilano@cern.ch>
Fri, 28 Mar 2014 14:26:34 +0000 (15:26 +0100)
committerlmilano <lmilano@cern.ch>
Fri, 28 Mar 2014 14:26:34 +0000 (15:26 +0100)
PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h

index 9883026..1ec8ed2 100644 (file)
@@ -142,6 +142,7 @@ fUseChargeHadrons(kFALSE),
 fParticleSpeciesTrigger(-1),
 fParticleSpeciesAssociated(-1),
 fCheckMotherPDG(kTRUE),
+fTrackletDphiCut(999.),
 fSelectCharge(0),
 fTriggerSelectCharge(0),
 fAssociatedSelectCharge(0),
@@ -337,6 +338,8 @@ void  AliAnalysisTaskPhiCorrelations::CreateOutputObjects()
   }
   if (fTriggersFromDetector == 1 || fTriggersFromDetector == 2)
     fListOfHistos->Add(new TH1F("V0SingleCells", "V0 single cell multiplicity;multiplicity;events", 100, -0.5, 99.5));
+  if (fTriggersFromDetector == 3)
+    fListOfHistos->Add(new TH1F("DphiTrklets", "tracklets Dphi;#Delta#phi,trklets;entries", 100, -0.1, 0.1));
   
   PostData(0,fListOfHistos);
   
@@ -1495,6 +1498,9 @@ TObjArray* AliAnalysisTaskPhiCorrelations::GetParticlesFromDetector(AliVEvent* i
        {
          Double_t eta=-TMath::Log(TMath::Tan(trklets->GetTheta(itrklets)/2));
          if(TMath::Abs(eta)>fTrackEtaCut)continue;
+         if(TMath::Abs(trklets->GetDeltaPhi(itrklets))>fTrackletDphiCut)continue;
+         TH1F* DphiTrklets = (TH1F*)fListOfHistos->FindObject("DphiTrklets");
+         DphiTrklets->Fill(trklets->GetDeltaPhi(itrklets));
          AliDPhiBasicParticle* particle = new AliDPhiBasicParticle(eta,trklets->GetPhi(itrklets), 1.1, 0); // fit pT = 1.1 and charge = 0
          particle->SetUniqueID(fAnalyseUE->GetEventCounter()* 100000 + itrklets);
          
index 9206221..c5882aa 100644 (file)
@@ -100,7 +100,10 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     void   SetFoundFractionCut(Double_t value) { fFoundFractionCut = value; }
     void   SetTrackStatus(UInt_t status)     { fTrackStatus = status; }
     void   SetCheckMotherPDG(Bool_t checkpdg) { fCheckMotherPDG = checkpdg; }
-    
+   
+    // track cuts
+    void   SetTrackletDphiCut( Double_t val )    { fTrackletDphiCut = val; }
+
     void   SetEventSelectionBit( UInt_t val )        { fSelectBit = val;  }
     void   SetUseChargeHadrons( Bool_t val ) { fUseChargeHadrons = val; }
     void   SetSelectParticleSpecies( Int_t trigger, Int_t associated ) { fParticleSpeciesTrigger = trigger; fParticleSpeciesAssociated = associated; }
@@ -215,6 +218,9 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     Int_t               fParticleSpeciesAssociated; // Select which particle to use for the associated [ -1 (all, default) 0 (pions) 1 (kaons) 2 (protons) 3 (others) particles ]
     Bool_t             fCheckMotherPDG;     // Check the PDG code of mother for secondaries 
 
+    // Tracklets cuts
+    Double_t           fTrackletDphiCut;    //maximum Dphi cut on tracklets
+    
     Int_t fSelectCharge;           // (un)like sign selection when building correlations: 0: no selection; 1: unlike sign; 2: like sign
     Int_t fTriggerSelectCharge;    // select charge of trigger particle: 1: positive; -1 negative
     Int_t fAssociatedSelectCharge; // select charge of associated particle: 1: positive; -1 negative
@@ -239,7 +245,7 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     
     Bool_t fFillpT;                // fill sum pT instead of number density
     
-    ClassDef(AliAnalysisTaskPhiCorrelations, 44); // Analysis task for delta phi correlations
+    ClassDef(AliAnalysisTaskPhiCorrelations, 45); // Analysis task for delta phi correlations
   };
 
 class AliDPhiBasicParticle : public AliVParticle