X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=ANALYSIS%2FAliCentralitySelectionTask.h;h=d1c82aa4d88358322e09b7040fbaa3cc8958ec15;hp=f3e5c5c913f5e691638fdcd25edc5fc14dc231c1;hb=981b4d38bb8499169c666fbaa7fc8608496c0f50;hpb=5233aeb53d246b4df4af5550c93f33ebac1ef8ca diff --git a/ANALYSIS/AliCentralitySelectionTask.h b/ANALYSIS/AliCentralitySelectionTask.h index f3e5c5c..d1c82aa 100644 --- a/ANALYSIS/AliCentralitySelectionTask.h +++ b/ANALYSIS/AliCentralitySelectionTask.h @@ -35,38 +35,40 @@ class AliCentralitySelectionTask : public AliAnalysisTaskSE { virtual void UserExec(Option_t *option); virtual void Terminate(Option_t *option); - virtual void SetDebugLevel(Int_t level) {fDebug = level;} void SetInput(const char* input) {fAnalysisInput = input;} void SetMCInput() {fIsMCInput = kTRUE;} + void SetPass(Int_t pass) {fPass = pass;} void DontUseScaling() {fUseScaling=kFALSE;} - + void DontUseCleaning() {fUseCleaning=kFALSE;} void ReadCentralityHistos(TString filename); void ReadCentralityHistos2(TString filename); private: - Int_t SetupRun(AliESDEvent* esd); - Bool_t IsOutlierV0MSPD(Float_t spd, Float_t v0, Int_t cent); - Bool_t IsOutlierV0MTPC(Int_t tracks, Float_t v0, Int_t cent); - Bool_t IsOutlierV0MZDC(Float_t zdc, Float_t v0); - Float_t MyGetScaleFactor(Int_t runnumber, Int_t flag); + Int_t SetupRun(AliESDEvent* const esd); + Bool_t IsOutlierV0MSPD(Float_t spd, Float_t v0, Int_t cent) const; + Bool_t IsOutlierV0MTPC(Int_t tracks, Float_t v0, Int_t cent) const; + Bool_t IsOutlierV0MZDC(Float_t zdc, Float_t v0) const; + Bool_t IsOutlierV0MZDCECal(Float_t zdc, Float_t v0) const; + Float_t MyGetScaleFactor(Int_t runnumber, Int_t flag) const; void MyInitScaleFactor(); - Float_t MyGetScaleFactorMC(Int_t runnumber); + Float_t MyGetScaleFactorMC(Int_t runnumber) const; void MyInitScaleFactorMC(); - Int_t fDebug; // Debug flag TString fAnalysisInput; // "ESD", "AOD" Bool_t fIsMCInput; // true when input is MC + Int_t fPass; // pass of reconstruction TFile *fFile; // file that holds the centrality vs multiplicity 1d TFile *fFile2; // file that holds the centrality vs multiplicity 2d Int_t fCurrentRun; // current run number Int_t fRunNo; // reference run number Int_t fLowRunN; // first run Int_t fHighRunN; // last run - Bool_t fUseScaling; - Float_t V0MScaleFactor[2667]; // number of runs in PbPb 2010 - Float_t SPDScaleFactor[2667]; // number of runs in PbPb 2010 - Float_t TPCScaleFactor[2667]; // number of runs in PbPb 2010 - Float_t V0MScaleFactorMC[2667]; // number of runs in PbPb 2010 + Bool_t fUseScaling; // flag to use scaling + Bool_t fUseCleaning; // flag to use cleanin + Float_t fV0MScaleFactor[2667]; // number of runs in PbPb 2010 + Float_t fSPDScaleFactor[2667]; // number of runs in PbPb 2010 + Float_t fTPCScaleFactor[2667]; // number of runs in PbPb 2010 + Float_t fV0MScaleFactorMC[2667]; // number of runs in PbPb 2010 AliESDtrackCuts* fTrackCuts; //! optional track cuts @@ -105,6 +107,10 @@ class AliCentralitySelectionTask : public AliAnalysisTaskSE { TH1F *fHOutCentV0MvsFMD; //control histogram for centrality TH1F *fHOutCentTKLvsV0M; //control histogram for centrality TH1F *fHOutCentZEMvsZDC; //control histogram for centrality + TH2F *fHOutCentV0MvsCentCL1; //control histogram for centrality + TH2F *fHOutCentV0MvsCentTRK; //control histogram for centrality + TH2F *fHOutCentTRKvsCentCL1; //control histogram for centrality + TH2F *fHOutCentV0MvsCentZDC; //control histogram for centrality TH1F *fHOutMultV0M ; //control histogram for multiplicity TH1F *fHOutMultV0R ; //control histogram for multiplicity @@ -118,22 +124,29 @@ class AliCentralitySelectionTask : public AliAnalysisTaskSE { TH2F *fHOutMultZEMvsZDN; //control histogram for multiplicity TH2F *fHOutMultV0MvsZDC; //control histogram for multiplicity TH2F *fHOutMultZEMvsZDC; //control histogram for multiplicity + TH2F *fHOutMultZEMvsZDCw; //control histogram for multiplicity TH2F *fHOutMultV0MvsCL1; //control histogram for multiplicity TH2F *fHOutMultV0MvsTRK; //control histogram for multiplicity TH2F *fHOutMultTRKvsCL1; //control histogram for multiplicity - TH1F *fHOutCentV0M_qual1 ; //control histogram for centrality quality 1 - TH1F *fHOutCentTRK_qual1 ; //control histogram for centrality quality 1 - TH1F *fHOutCentCL1_qual1 ; //control histogram for centrality quality 1 + TH1F *fHOutCentV0Mqual1 ; //control histogram for centrality quality 1 + TH1F *fHOutCentTRKqual1 ; //control histogram for centrality quality 1 + TH1F *fHOutCentCL1qual1 ; //control histogram for centrality quality 1 + TH2F *fHOutMultV0MvsCL1qual1; //control histogram for multiplicity quality 1 + TH2F *fHOutMultV0MvsTRKqual1; //control histogram for multiplicity quality 1 + TH2F *fHOutMultTRKvsCL1qual1; //control histogram for multiplicity quality 1 - TH1F *fHOutCentV0M_qual2 ; //control histogram for centrality quality 2 - TH1F *fHOutCentTRK_qual2 ; //control histogram for centrality quality 2 - TH1F *fHOutCentCL1_qual2 ; //control histogram for centrality quality 2 + TH1F *fHOutCentV0Mqual2 ; //control histogram for centrality quality 2 + TH1F *fHOutCentTRKqual2 ; //control histogram for centrality quality 2 + TH1F *fHOutCentCL1qual2 ; //control histogram for centrality quality 2 + TH2F *fHOutMultV0MvsCL1qual2; //control histogram for multiplicity quality 2 + TH2F *fHOutMultV0MvsTRKqual2; //control histogram for multiplicity quality 2 + TH2F *fHOutMultTRKvsCL1qual2; //control histogram for multiplicity quality 2 TH1F *fHOutQuality ; //control histogram for quality TH1F *fHOutVertex ; //control histogram for vertex - ClassDef(AliCentralitySelectionTask, 6); + ClassDef(AliCentralitySelectionTask, 10); }; #endif