Updated impact parameter resolution for 2010 tracks (L. Milano)

diff --git a/PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.cxx b/PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.cxx
index 37483a6f876ec18244cf653929fb468cf59ba507..d2047ecbcf66b6ae2b6dabc25f1bc002752f43cf 100644 (file)
--- a/PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.cxx
+++ b/PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.cxx
@@ -45,6 +45,7 @@
 #include "AliESDtrackCuts.h"
 
 ClassImp(AliAnalysisTaskSEITSsaSpectra)
 #include "AliESDtrackCuts.h"
 
 ClassImp(AliAnalysisTaskSEITSsaSpectra)

+/* $Id$ */

 
 //________________________________________________________________________
 AliAnalysisTaskSEITSsaSpectra::AliAnalysisTaskSEITSsaSpectra():
 
 //________________________________________________________________________
 AliAnalysisTaskSEITSsaSpectra::AliAnalysisTaskSEITSsaSpectra():


@@ -70,6 +71,7 @@ AliAnalysisTaskSE("Task CFits"),
   fNSigmaDCAz(7.),
   fLowMult(0),
   fUpMult(9999),
   fNSigmaDCAz(7.),
   fLowMult(0),
   fUpMult(9999),

+  fYear(2010),

   fMC(kFALSE), 
   fSmearMC(kFALSE),
   fSmearP(0.),
   fMC(kFALSE), 
   fSmearMC(kFALSE),
   fSmearP(0.),


@@ -145,20 +147,38 @@ Bool_t AliAnalysisTaskSEITSsaSpectra::DCAcut(Double_t impactXY, Double_t impactZ
   Double_t xyP[3];
   Double_t zP[3];
   if(optMC){
   Double_t xyP[3];
   Double_t zP[3];
   if(optMC){

-    xyP[0]=88.63;//SIMpass4a12
-    xyP[1]=19.57;
-    xyP[2]=1.65;
-    zP[0]=140.98;
-    zP[1]=62.33;
-    zP[2]=1.15;
+    if(fYear==2009){
+      xyP[0]=88.63; //MC LHC10a12
+      xyP[1]=19.57;
+      xyP[2]=1.65;
+      zP[0]=140.98;
+      zP[1]=62.33;
+      zP[2]=1.15;
+    }else{
+      xyP[0]=36.; //MC LHC10d1
+      xyP[1]=43.9;
+      xyP[2]=1.3;
+      zP[0]=111.9;
+      zP[1]=59.8;
+      zP[2]=1.2;
+    }

   }
   else{
   }
   else{

-    xyP[0]=85.28;//DATApass6
-    xyP[1]=25.78;
-    xyP[2]=1.55;
-    zP[0]=146.80;
-    zP[1]=70.07;
-    zP[2]=1.11;
+    if(fYear==2009){
+      xyP[0]=85.28;//DATA 900 GeV pass6
+      xyP[1]=25.78;
+      xyP[2]=1.55;
+      zP[0]=146.80;
+      zP[1]=70.07;
+      zP[2]=1.11;
+    }else{
+      xyP[0]=32.7;//DATA 7 TeV pass2
+      xyP[1]=44.8;
+      xyP[2]=1.3;
+      zP[0]=117.3;
+      zP[1]=66.8;
+      zP[2]=1.2;
+    }

   }
   Double_t xySigma = xyP[0] + xyP[1]/TMath::Power(TMath::Abs(pt),xyP[2]);
   Double_t xyMax = fNSigmaDCAxy*xySigma; //in micron
   }
   Double_t xySigma = xyP[0] + xyP[1]/TMath::Power(TMath::Abs(pt),xyP[2]);
   Double_t xyMax = fNSigmaDCAxy*xySigma; //in micron


@@ -365,7 +385,7 @@ void AliAnalysisTaskSEITSsaSpectra::UserCreateOutputObjects(){
     fOutput->Add(fHistNegK[i]);
     fOutput->Add(fHistNegP[i]);
     
     fOutput->Add(fHistNegK[i]);
     fOutput->Add(fHistNegP[i]);
     

-    fOutput->Add(fHistDCAPosPi[i]);//DCA distr.
+    fOutput->Add(fHistDCAPosPi[i]); //DCA distr

     fOutput->Add(fHistDCAPosK[i]);
     fOutput->Add(fHistDCAPosP[i]);
     fOutput->Add(fHistDCANegPi[i]);
     fOutput->Add(fHistDCAPosK[i]);
     fOutput->Add(fHistDCAPosP[i]);
     fOutput->Add(fHistDCANegPi[i]);


@@ -471,8 +491,13 @@ void AliAnalysisTaskSEITSsaSpectra::UserExec(Option_t *){
   esdTrackCutsMult->SetRequireITSRefit(kTRUE);
   esdTrackCutsMult->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
                                             AliESDtrackCuts::kAny);
   esdTrackCutsMult->SetRequireITSRefit(kTRUE);
   esdTrackCutsMult->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
                                             AliESDtrackCuts::kAny);

-  esdTrackCutsMult->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/TMath::Power(pt,0.9)");  
-  esdTrackCutsMult->SetMaxDCAToVertexZ(20);
+  if(fYear==2009){
+    esdTrackCutsMult->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/TMath::Power(pt,0.9)");  //2009 standard cut
+    esdTrackCutsMult->SetMaxDCAToVertexZ(20); //2009 standard cut
+  }else{
+    esdTrackCutsMult->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01"); //2010 standard cut
+    esdTrackCutsMult->SetMaxDCAToVertexZ(2); //2010 standard cut
+  }

   esdTrackCutsMult->SetDCAToVertex2D(kFALSE);
   esdTrackCutsMult->SetRequireSigmaToVertex(kFALSE);
   esdTrackCutsMult->SetEtaRange(-0.8,+0.8);
   esdTrackCutsMult->SetDCAToVertex2D(kFALSE);
   esdTrackCutsMult->SetRequireSigmaToVertex(kFALSE);
   esdTrackCutsMult->SetEtaRange(-0.8,+0.8);

diff --git a/PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.h b/PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.h
index 79fbbecc5b33db12968591fffb80ead94cd8085b..baf8334be94b55bf7f4c2776f7d89a93e981b7e2 100644 (file)
--- a/PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.h
+++ b/PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.h
@@ -12,6 +12,8 @@
 // F. Prino, prino@to.infn.it
 ///////////////////////////////////////////////////////////////////////////
 
 // F. Prino, prino@to.infn.it
 ///////////////////////////////////////////////////////////////////////////
 

+/* $Id$ */
+

 class TString;
 class TTree;
 class TH1F;
 class TString;
 class TTree;
 class TH1F;


@@ -58,6 +60,7 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
     fUpMult=UpBin;
   }
 
     fUpMult=UpBin;
   }
 

+  void SetYear(Int_t year){fYear=year;}

   void SetReadMC(Bool_t flag = kTRUE) {fMC = flag;}
   void SetFillNtuple(Bool_t fill=kTRUE) {fFillNtuple=fill;}
   void SetSmearMC(Double_t smearp, Double_t smeardedx){
   void SetReadMC(Bool_t flag = kTRUE) {fMC = flag;}
   void SetFillNtuple(Bool_t fill=kTRUE) {fFillNtuple=fill;}
   void SetSmearMC(Double_t smearp, Double_t smeardedx){


@@ -118,7 +121,7 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
   TH1F *fHistNegK[kNbins]; //! histo with dedx distibution in the kaons hypotesis (negative)
   TH1F *fHistNegP[kNbins]; //! histo with dedx distibution in the protons hypotesis (negative)
   
   TH1F *fHistNegK[kNbins]; //! histo with dedx distibution in the kaons hypotesis (negative)
   TH1F *fHistNegP[kNbins]; //! histo with dedx distibution in the protons hypotesis (negative)
   

-  TH1F *fHistDCAPosPi[kNbins]; //! histo with DCA distibution in the pions hypotesis (positive)
+  TH1F *fHistDCAPosPi[kNbins]; //! histo with DCA distibution in the kaons hypotesis (positive)

   TH1F *fHistDCAPosK[kNbins]; //! histo with DCA distibution in the kaons hypotesis (positive)
   TH1F *fHistDCAPosP[kNbins]; //! histo with DCA distibution in the protons hypotesis (positive)
   TH1F *fHistDCANegPi[kNbins]; //! histo with DCA distibution in the pions hypotesis (negative)
   TH1F *fHistDCAPosK[kNbins]; //! histo with DCA distibution in the kaons hypotesis (positive)
   TH1F *fHistDCAPosP[kNbins]; //! histo with DCA distibution in the protons hypotesis (positive)
   TH1F *fHistDCANegPi[kNbins]; //! histo with DCA distibution in the pions hypotesis (negative)


@@ -182,17 +185,17 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
   Double_t fNSigmaDCAz;      // DCA cut along z
   Int_t fLowMult;      // Multiplicity bin
   Int_t fUpMult;      // Multiplicity bin
   Double_t fNSigmaDCAz;      // DCA cut along z
   Int_t fLowMult;      // Multiplicity bin
   Int_t fUpMult;      // Multiplicity bin

- 
+  Int_t fYear;        // Year (2009, 2010)

   Bool_t   fMC;        //flag to switch on the MC analysis for the efficiency estimation
   Bool_t   fSmearMC;   // flag to apply extra smearing on MC 
   Double_t fSmearP;    // extra relative smearing on simulated momentum
   Double_t fSmeardEdx; // extra relative smearing on simulated dE/dx
   Bool_t   fMC;        //flag to switch on the MC analysis for the efficiency estimation
   Bool_t   fSmearMC;   // flag to apply extra smearing on MC 
   Double_t fSmearP;    // extra relative smearing on simulated momentum
   Double_t fSmeardEdx; // extra relative smearing on simulated dE/dx

-  TRandom3* fRandGener;    // generator for smaring
+  TRandom3* fRandGener;    // generator for smearing

   Bool_t       fFillNtuple;      // fill ntuple  
   TNtuple     *fNtupleNSigma;//! output ntuple
   TNtuple     *fNtupleMC;//! output MC ntuple
   
   Bool_t       fFillNtuple;      // fill ntuple  
   TNtuple     *fNtupleNSigma;//! output ntuple
   TNtuple     *fNtupleMC;//! output MC ntuple
   

-  ClassDef(AliAnalysisTaskSEITSsaSpectra, 1);
+  ClassDef(AliAnalysisTaskSEITSsaSpectra, 2);

 };
 
 #endif
 };
 
 #endif