add setter for dca cut parameters, initialize missing histograms in ctor
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 14 Jan 2013 10:04:23 +0000 (10:04 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 14 Jan 2013 10:04:23 +0000 (10:04 +0000)
PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.cxx
PWGGA/CaloTrackCorrelations/AliAnaChargedParticles.h

index 42ce812..f606eb3 100755 (executable)
@@ -87,6 +87,7 @@ ClassImp(AliAnaChargedParticles)
   
   for(Int_t i = 0; i < 3; i++)
   {
+    fDCACutParam          [i] = 0 ;
     fhPtDCA               [i] = 0 ;
     
     fhPtDCASPDRefit       [i] = 0 ;
@@ -104,8 +105,10 @@ ClassImp(AliAnaChargedParticles)
     fhPtDCAVtxInBC0       [i] = 0 ;
     fhPtDCAVtxOutBC0PileUp[i] = 0 ;
     fhPtDCAVtxInBC0PileUp [i] = 0 ;
-    fhPtDCAVtxOutBC0NoTOFHit      [i] = 0 ;
+    fhPtDCAVtxOutBC0NoTOFHit[i] = 0 ;
+    fhPtDCAVtxInBC0NoTOFHit [i] = 0 ;
     fhPtDCAVtxOutBC0PileUpNoTOFHit[i] = 0 ;
+    fhPtDCAVtxInBC0PileUpNoTOFHit [i] = 0 ;
   }
   
   //Initialize parameters
@@ -118,7 +121,7 @@ Bool_t  AliAnaChargedParticles::AcceptDCA(const Float_t pt, const Float_t dca)
 {
  // Accept track if DCA is smaller than function
   
-  Float_t cut = 0.0105+0.0350/TMath::Power(pt,1.1);
+  Float_t cut = fDCACutParam[0]+fDCACutParam[1]/TMath::Power(pt,fDCACutParam[2]);
   
   if(TMath::Abs(dca) < cut)
     return kTRUE;
@@ -658,6 +661,11 @@ void AliAnaChargedParticles::InitParameters()
 
   AddToHistogramsName("AnaCharged_");
   
+  // dca_xy cut = 0.0105+0.0350/TMath::Power(pt,1.1);
+  fDCACutParam[0]= 0.0105;
+  fDCACutParam[1]= 0.0350;
+  fDCACutParam[2]= 1.1;
+  
 }
 
 //____________________________________________________________
index ec50991..3536b63 100755 (executable)
@@ -48,11 +48,14 @@ class AliAnaChargedParticles : public AliAnaCaloTrackCorrBaseClass {
   void    SwitchOnRecalculateVertexBC()      { fRecalculateVertexBC     = kTRUE  ; }
   void    SwitchOffRecalculateVertexBC()     { fRecalculateVertexBC     = kFALSE ; }
 
+  void    SetDCACutParameters(Int_t i, Float_t par) { if(i >= 0 && i < 3) fDCACutParam[i] = par ; }
+  
  private:
   
-  Bool_t fFillPileUpHistograms;    // Fill pile-up related histograms
-  Bool_t fFillVertexBC0Histograms; // Fill histograms for tracks with vertex BC=0 or not related histograms
-  Bool_t fRecalculateVertexBC;     // Recalculate vertex BC for older AODs
+  Bool_t  fFillPileUpHistograms;    // Fill pile-up related histograms
+  Bool_t  fFillVertexBC0Histograms; // Fill histograms for tracks with vertex BC=0 or not related histograms
+  Bool_t  fRecalculateVertexBC;     // Recalculate vertex BC for older AODs
+  Float_t fDCACutParam[3];          // DCA cut function parameters
   
   //Histograms
   TH1F * fhNtracks;     //! track multiplicity distribution