]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG2/SPECTRA/AliAnalysisTaskSEITSsaSpectra.h
Updated impact parameter resolution for 2010 tracks (L. Milano)
[u/mrichter/AliRoot.git] / PWG2 / SPECTRA / AliAnalysisTaskSEITSsaSpectra.h
index e6b7214e4066588b3669809fb35298c6258cbc55..baf8334be94b55bf7f4c2776f7d89a93e981b7e2 100644 (file)
@@ -12,6 +12,8 @@
 // F. Prino, prino@to.infn.it
 ///////////////////////////////////////////////////////////////////////////
 
+/* $Id$ */
+
 class TString;
 class TTree;
 class TH1F;
@@ -31,13 +33,13 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
   virtual void   UserExec(Option_t *);
   virtual void   Terminate(Option_t *);
 
-  void SetMaxChi2Clu(Double_t chi=1.){
+  void SetMaxChi2Clu(Double_t chi=2.5){
     fMaxChi2Clu=chi;
   }
   void SetRapidityRange(Double_t dy=0.5){
     fMaxY=dy;
   }
-  void SetMinNSigma(Double_t ns=3.){
+  void SetMinNSigma(Double_t ns=1.5){
     fMinNSigma=ns;
   }
   void SetMindEdx(Double_t mind=0.){
@@ -53,7 +55,12 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
     fNSigmaDCAxy=nsxy;
     fNSigmaDCAz=nsz; 
   }
+  void SetMultBin(Int_t LowBin=0,Int_t UpBin=9999){
+    fLowMult=LowBin;
+    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){
@@ -75,7 +82,10 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
   AliESDEvent *fESD; //ESD object
   TList *fOutput; //! tlist with output
   TH1F *fHistNEvents; //! histo with number of events
+  TH1F *fHistMult; //! histo with multiplicity of the events
   TH1F *fHistNTracks; //! histo with number of tracks
+  TH1F *fHistNTracksPos; //! histo with number of tracks
+  TH1F *fHistNTracksNeg; //! histo with number of tracks
   
   TH2F *fHistDEDX; //! histo with dedx versus momentum
   TH2F *fHistDEDXdouble; //! histo with dedx versus signed momentum
@@ -83,10 +93,24 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
   TH1F *fHistBeforeEvSel; //! histo with pt distribution before my event selection
   TH1F *fHistAfterEvSel; //! histo with pt distribution after my event selection
   
-  TH1F *fHistMCpos[3]; //! histo with spectra from the MC truth (positive)
-  TH1F *fHistMCneg[3]; //! histo with spectra from the MC truth (negative)
-  TH1F *fHistMCposBefEvSel[3]; //! histo with spectra from the MC truth (positive)
-  TH1F *fHistMCnegBefEvSel[3]; //! histo with spectra from the MC truth (negative)
+  TH1F *fHistPrimMCpos[3]; //! histo with spectra of primaries from the MC truth (positive)
+  TH1F *fHistPrimMCneg[3]; //! histo with spectra of primaries from the MC truth (negative)
+  TH1F *fHistSecStrMCpos[3]; //! histo with spectra of strange decays from the MC truth (positive)
+  TH1F *fHistSecStrMCneg[3]; //! histo with spectra of strange decays from the MC truth (negative)
+  TH1F *fHistSecMatMCpos[3]; //! histo with spectra of sec. from material from the MC truth (positive)
+  TH1F *fHistSecMatMCneg[3]; //! histo with spectra of sec. from material from the MC truth (negative)
+  TH1F *fHistPrimMCposBefEvSel[3]; //! histo with spectra of primaries from the MC truth (positive)
+  TH1F *fHistPrimMCnegBefEvSel[3]; //! histo with spectra of primaries from the MC truth (negative)
+  TH1F *fHistSecStrMCposBefEvSel[3]; //! histo with spectra of strange decays from the MC truth (positive)
+  TH1F *fHistSecStrMCnegBefEvSel[3]; //! histo with spectra of strange decays from the MC truth (negative)
+  TH1F *fHistSecMatMCposBefEvSel[3]; //! histo with spectra of sec. from material from the MC truth (positive)
+  TH1F *fHistSecMatMCnegBefEvSel[3]; //! histo with spectra of sec. from material from the MC truth (negative)
+  TH1F *fHistPrimMCposReco[3]; //! histo with spectra of primaries from the MC truth (positive)
+  TH1F *fHistPrimMCnegReco[3]; //! histo with spectra of primaries from the MC truth (negative)
+  TH1F *fHistSecStrMCposReco[3]; //! histo with spectra of strange decays from the MC truth (positive)
+  TH1F *fHistSecStrMCnegReco[3]; //! histo with spectra of strange decays from the MC truth (negative)
+  TH1F *fHistSecMatMCposReco[3]; //! histo with spectra of sec. from material from the MC truth (positive)
+  TH1F *fHistSecMatMCnegReco[3]; //! histo with spectra of sec. from material from the MC truth (negative)
   
   TH1F *fHistCharge[4]; //! histo with charge distribution to check the calibration 
   
@@ -97,13 +121,34 @@ 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)
   TH1F *fHistDCANegK[kNbins]; //! histo with DCA distibution in the kaons hypotesis (negative)
   TH1F *fHistDCANegP[kNbins]; //! histo with DCA distibution in the protons hypotesis (negative)
-               
+  
+  TH1F *fHistMCPrimDCAPosPi[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCPrimDCAPosK[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCPrimDCAPosP[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCPrimDCANegPi[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCPrimDCANegK[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCPrimDCANegP[kNbins]; //! histo with DCA distibution, MC truth
+  
+  TH1F *fHistMCSecStDCAPosPi[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecStDCAPosK[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecStDCAPosP[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecStDCANegPi[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecStDCANegK[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecStDCANegP[kNbins]; //! histo with DCA distibution, MC truth
+  
+  TH1F *fHistMCSecMatDCAPosPi[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecMatDCAPosK[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecMatDCAPosP[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecMatDCANegPi[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecMatDCANegK[kNbins]; //! histo with DCA distibution, MC truth 
+  TH1F *fHistMCSecMatDCANegP[kNbins]; //! histo with DCA distibution, MC truth
   TH1F *fHistMCPosPi[kNbins]; //! histo with dedx using the MC truth in the pions hypotesis (positive)
   TH1F *fHistMCPosK[kNbins]; //! histo with dedx  using the MC truth in the kaons hypotesis (positive)
   TH1F *fHistMCPosP[kNbins]; //! histo with dedx  using the MC truth in the protons hypotesis (positive)
@@ -111,10 +156,21 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
   TH1F *fHistMCNegK[kNbins]; //! histo with dedx  using the MC truth in the kaons hypotesis (negative)
   TH1F *fHistMCNegP[kNbins]; //! histo with dedx  using the MC truth in the protons hypotesis (negative)
   
+  TH1F *fHistPosNSigmaMean[3];       //! NSigma histos for 6 species
+  TH1F *fHistPosNSigmaMCMean[3];       //! NSigma histos for 6 species
+  TH1F *fHistPosNSigmaPrimMean[3];   //! NSigma histos for 6 species
+  TH1F *fHistPosNSigmaPrimMCMean[3]; //! NSigma histos for 6 species
+  TH1F *fHistNegNSigmaMean[3];       //! NSigma histos for 6 species
+  TH1F *fHistNegNSigmaMCMean[3];       //! NSigma histos for 6 species
+  TH1F *fHistNegNSigmaPrimMean[3];   //! NSigma histos for 6 species
+  TH1F *fHistNegNSigmaPrimMCMean[3]; //! NSigma histos for 6 species
+
   TH1F *fHistPosNSigma[3];       //! NSigma histos for 6 species
+  TH1F *fHistPosNSigmaMC[3];       //! NSigma histos for 6 species
   TH1F *fHistPosNSigmaPrim[3];   //! NSigma histos for 6 species
   TH1F *fHistPosNSigmaPrimMC[3]; //! NSigma histos for 6 species
   TH1F *fHistNegNSigma[3];       //! NSigma histos for 6 species
+  TH1F *fHistNegNSigmaMC[3];       //! NSigma histos for 6 species
   TH1F *fHistNegNSigmaPrim[3];   //! NSigma histos for 6 species
   TH1F *fHistNegNSigmaPrimMC[3]; //! NSigma histos for 6 species
 
@@ -127,17 +183,19 @@ class AliAnalysisTaskSEITSsaSpectra : public AliAnalysisTaskSE {
   Double_t fMaxChi2Clu;      // maximum cluster
   Double_t fNSigmaDCAxy;     // DCA cut in bend. plane
   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