The code we used for the analysis of the first PbPb mult paper

[u/mrichter/AliRoot.git] / PWG2 / EVCHAR / AliAnalysisTaskSPDdNdEta.h

diff --git a/PWG2/EVCHAR/AliAnalysisTaskSPDdNdEta.h b/PWG2/EVCHAR/AliAnalysisTaskSPDdNdEta.h
index ce6009bcbd546a0486d55f8f0c523f5e75d6f303..db306fd7014fba7f5ffc9ad9eb87fbea6dad9c8b 100644 (file)
--- a/PWG2/EVCHAR/AliAnalysisTaskSPDdNdEta.h
+++ b/PWG2/EVCHAR/AliAnalysisTaskSPDdNdEta.h
@@ -38,8 +38,13 @@ class AliAnalysisTaskSPDdNdEta : public AliAnalysisTaskSE {
   void SetMCCentralityBin(MCCentralityBin mccentrbin) {fMCCentralityBin=mccentrbin;}
   void SetCentralityLowLim(Float_t centrlowlim) {fCentrLowLim=centrlowlim;}
   void SetCentralityUpLim(Float_t centruplim) {fCentrUpLim=centruplim;}
-  void SetCentralityEst(TString centrest) {fCentrEst=centrest;}
+//  void SetCentralityEst(TString centrest) {fCentrEst=centrest;}
+  void SetCentralityEst(Bool_t centrest) {fCentrEst=centrest;}
   void SetMinClusterMultLay2(Int_t minClMultLay2=0) {fMinClMultLay2=minClMultLay2;}
+  void SetMaxClusterMultLay2(Int_t maxClMultLay2=0) {fMaxClMultLay2=maxClMultLay2;}
+  void SetMinV0Mult(Int_t minV0Mult=0) {fMinMultV0=minV0Mult;}
+  void SetVertexRange(Float_t vl=7.) {fVtxLim=vl;}
+  void SetPartSpecies(Bool_t partsp = kFALSE)  {fPartSpecies=partsp;}
 
   void SetPhiWindow(Float_t w=0.08) {fPhiWindow=w;}
   void SetThetaWindow(Float_t w=0.025) {fThetaWindow=w;}
@@ -58,7 +63,6 @@ class AliAnalysisTaskSPDdNdEta : public AliAnalysisTaskSE {
   TList *fOutput;                  // ! output list send on output slot 1 
 
   Bool_t fUseMC;                   // flag to enable the calculation of correction histograms
-  Bool_t fPbPb;                    // flag to analyze PbPb data 
   AliTriggerAnalysis::Trigger fTrigger;  
   Bool_t fTR;                      // to read track references and calculate factors of track to particle correction 
   Bool_t fRecoTracklets;           // flag to recostruct tracklets
@@ -66,9 +70,15 @@ class AliAnalysisTaskSPDdNdEta : public AliAnalysisTaskSE {
   MCCentralityBin fMCCentralityBin; // to select MC centrality bin in which corrections are calculated
   Float_t fCentrLowLim;             // to select centrality bin on data
   Float_t fCentrUpLim;              // to select centrality bin on data
-  TString fCentrEst;                // to select centrality estimator
-  Int_t fMinClMultLay2;             // to select multiplicity class
+//  TString fCentrEst;                 // to select centrality estimator
+  Bool_t fCentrEst;                 // to select centrality estimator
 
+  Int_t fMinClMultLay2;             // to select multiplicity class
+  Int_t fMaxClMultLay2;             // to select multiplicity class
+  Int_t fMinMultV0;                 // to select centrality class 
+  Float_t fVtxLim;                  // to select vertex range
+  Bool_t fPartSpecies;              // to fill correction matrices for each part species (for syst studies)
+ 
   Float_t       fPhiWindow;                    // Search window in phi
   Float_t       fThetaWindow;                  // Search window in theta
   Float_t       fPhiShift;                     // Phi shift reference value (at 0.5 T)
@@ -80,6 +90,8 @@ class AliAnalysisTaskSPDdNdEta : public AliAnalysisTaskSE {
 
   AliTrackletAlg *fMultReco;       // tracklet reconstruction class
 
+
+  TH1F        *fV0Ampl;                     // ! V0 amplitudes to cut on centrality
   TH2F        *fHistSPDRAWMultvsZ;          // ! data to be corrected 
   TH2F        *fHistSPDRAWMultvsZTriggCentrEvts; // ! data to be corrected
   TH2F        *fHistSPDRAWMultvsZCentrEvts; // ! data to be corrected
@@ -91,8 +103,8 @@ class AliAnalysisTaskSPDdNdEta : public AliAnalysisTaskSE {
   TH1F        *fHistSPDmultcl1;             // ! cluster inner layer and tracklet check histos
   TH1F        *fHistSPDmultcl2;             // ! cluster inner layer and tracklet check histos
   TH2F        *fHistSPDmultcl1vscl2;        // ! cluster inner layer and tracklet check histos
-  TH2F        *fHistSPDmultvscl1;        // ! cluster inner layer and tracklet check histos
-  TH2F        *fHistSPDmultvscl2;        // ! cluster inner layer and tracklet check histos
+  TH2F        *fHistSPDmultvscl1;           // ! cluster inner layer and tracklet check histos
+  TH2F        *fHistSPDmultvscl2;           // ! cluster inner layer and tracklet check histos
 
   TH1F        *fHistSPDeta;                 // ! cluster inner layer and tracklet check histos
   TH1F        *fHistSPDphi;                 // ! cluster inner layer and tracklet check histos
@@ -105,16 +117,27 @@ class AliAnalysisTaskSPDdNdEta : public AliAnalysisTaskSE {
   TH1F        *fHistSPDvtxAnalysis;         // ! SPD vertex distributions
   TH2F        *fHistSPDdePhideTheta;        // ! histogram for combinatorial background studies
 
-  TH1F        *fHistSPDphicl1;                 // ! cluster inner layer and tracklet check histos
-  TH1F        *fHistSPDphicl2;                 // ! cluster inner layer and tracklet check histos
+  TH1F        *fHistSPDphicl1;              // ! cluster inner layer and tracklet check histos
+  TH1F        *fHistSPDphicl2;              // ! cluster inner layer and tracklet check histos
+
+  TH2F* fHistBkgCorrDen;                    // ! track level correction histograms
+  TH2F* fHistReconstructedProtons;            // ! track level correction histograms
+  TH2F* fHistReconstructedKaons;            // ! track level correction histograms
+  TH2F* fHistReconstructedPions;            // ! track level correction histograms
+  TH2F* fHistReconstructedOthers;           // ! track level correction histograms
+  TH2F* fHistReconstructedSec;        // ! track level correction histograms
 
-  TH2F* fHistBkgCorrDen;             // ! track level correction histograms
   TH2F* fHistBkgCorrDenPrimGen;      // ! track level correction histograms
+  TH2F* fHistBkgCombLabels;          // ! track level correction histograms
   TH2F* fHistBkgCorrNum;             // ! track level correction histograms
   TH2F* fHistAlgEffNum;              // ! track level correction histograms
   TH2F* fHistNonDetectableCorrDen;   // ! track level correction histograms
 
   TH2F* fHistNonDetectableCorrNum;   // ! track level correction histograms
+  TH2F* fHistProtons;                // ! track level correction histograms
+  TH2F* fHistKaons;                  // ! track level correction histograms
+  TH2F* fHistPions;                  // ! track level correction histograms
+  TH2F* fHistOthers;                 // ! track level correction histograms
   TH2F* fHistAllPrimaries;           // ! track level correction histograms
   TH2F* fHistTrackCentrEvts;         // ! track level correction histograms
   TH2F* fHistTrackTrigCentrEvts;     // ! track level correction histograms
@@ -124,10 +147,10 @@ class AliAnalysisTaskSPDdNdEta : public AliAnalysisTaskSE {
   TH2F* fHistTrigCentrEvts;          // ! event level correction histograms
   TH2F* fHistSelEvts;                // ! event level correction histograms
 
-  TH1F* fHistMCmultEtacut;           // ! MC distributions
+  TH1F* fHistMCmultEtacut;                // ! MC distributions
   TH2F* fHistMCmultEtacutvsSPDmultEtacut; // ! MC distributions
-  TH2F* fHistMCmultEtacutvsSPDmultcl1;  // ! MC distributions
-  TH2F* fHistMCmultEtacutvsSPDmultcl2;  // ! MC distributions
+  TH2F* fHistMCmultEtacutvsSPDmultcl1;    // ! MC distributions
+  TH2F* fHistMCmultEtacutvsSPDmultcl2;    // ! MC distributions
 
   TH1F* fHistMCvtxx;                 // ! MC vertex
   TH1F* fHistMCvtxy;                 // ! MC vertex
@@ -136,9 +159,9 @@ class AliAnalysisTaskSPDdNdEta : public AliAnalysisTaskSE {
   TH2F* fHistRecvsGenImpactPar;      // ! impact parameter correlation (ZDC est vs gen)
   TH1F* fHistMCNpart;                // ! distribution of number of participants from MC 
 
-  TH1F* fHistdPhiPP;                 // ! tracklet check histo
-  TH1F* fHistdPhiSS;                 // ! tracklet check histo
-  TH1F* fHistdPhiComb;               // ! tracklet check histo
+  TH2F* fHistdPhidThetaPP;           // ! tracklet check histo
+  TH2F* fHistdPhidThetaSS;           // ! tracklet check histo
+  TH2F* fHistdPhidThetaComb;         // ! tracklet check histo
 
   TH2F* fHistDeVtx;                  // ! check histo
 
@@ -146,7 +169,7 @@ class AliAnalysisTaskSPDdNdEta : public AliAnalysisTaskSE {
   AliAnalysisTaskSPDdNdEta(const AliAnalysisTaskSPDdNdEta&); // not implemented
   AliAnalysisTaskSPDdNdEta& operator=(const AliAnalysisTaskSPDdNdEta&); // not implemented 
   
-  ClassDef(AliAnalysisTaskSPDdNdEta, 3);  
+  ClassDef(AliAnalysisTaskSPDdNdEta, 5);  
 };
 
-#endif
+#endif
\ No newline at end of file