Updated impact parameter resolution for 2010 tracks (L. Milano)
authorprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 12 Nov 2010 18:18:22 +0000 (18:18 +0000)
committerprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 12 Nov 2010 18:18:22 +0000 (18:18 +0000)
PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.cxx
PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.h

index 37483a6..d2047ec 100644 (file)
@@ -45,6 +45,7 @@
 #include "AliESDtrackCuts.h"
 
 ClassImp(AliAnalysisTaskSEITSsaSpectra)
+/* $Id$ */
 
 //________________________________________________________________________
 AliAnalysisTaskSEITSsaSpectra::AliAnalysisTaskSEITSsaSpectra():
@@ -70,6 +71,7 @@ AliAnalysisTaskSE("Task CFits"),
   fNSigmaDCAz(7.),
   fLowMult(0),
   fUpMult(9999),
+  fYear(2010),
   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){
-    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{
-    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
@@ -365,7 +385,7 @@ void AliAnalysisTaskSEITSsaSpectra::UserCreateOutputObjects(){
     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]);
@@ -471,8 +491,13 @@ void AliAnalysisTaskSEITSsaSpectra::UserExec(Option_t *){
   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);
index 79fbbec..baf8334 100644 (file)
@@ -12,6 +12,8 @@
 // F. Prino, prino@to.infn.it
 ///////////////////////////////////////////////////////////////////////////
 
+/* $Id$ */
+
 class TString;
 class TTree;
 class TH1F;
@@ -58,6 +60,7 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
     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){
@@ -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 *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)
@@ -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
+  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
-  TRandom3* fRandGener;    // generator for smaring
+  TRandom3* fRandGener;    // generator for smearing
   Bool_t       fFillNtuple;      // fill ntuple  
   TNtuple     *fNtupleNSigma;//! output ntuple
   TNtuple     *fNtupleMC;//! output MC ntuple
   
-  ClassDef(AliAnalysisTaskSEITSsaSpectra, 1);
+  ClassDef(AliAnalysisTaskSEITSsaSpectra, 2);
 };
 
 #endif