]>
Commit | Line | Data |
---|---|---|
b2a297fa | 1 | #ifndef ALIDIELECTRONCF_H |
2 | #define ALIDIELECTRONCF_H | |
3 | /* Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | //############################################################# | |
7 | //# # | |
8 | //# Class AliDielectronCF # | |
9 | //# Dielectron Correction Framework Manager # | |
10 | //# # | |
11 | //# Authors: # | |
12 | //# Anton Andronic, GSI / A.Andronic@gsi.de # | |
13 | //# Ionut C. Arsene, GSI / I.C.Arsene@gsi.de # | |
14 | //# Julian Book, Uni Ffm / Julian.Book@cern.ch # | |
15 | //# Frederick Kramer, Uni Ffm, / Frederick.Kramer@cern.ch # | |
16 | //# Magnus Mager, CERN / Magnus.Mager@cern.ch # | |
17 | //# WooJin J. Park, GSI / W.J.Park@gsi.de # | |
18 | //# Jens Wiechula, Uni HD / Jens.Wiechula@cern.ch # | |
19 | //# # | |
20 | //############################################################# | |
21 | ||
22 | ||
23 | ||
24 | #include <TNamed.h> | |
61d106d3 | 25 | #include <TVectorDfwd.h> |
b2a297fa | 26 | #include "AliDielectronVarManager.h" |
27 | ||
28 | class AliAnalysisCuts; | |
29 | class AliAnalysisFilter; | |
30 | class AliCFContainer; | |
6551594b | 31 | class AliDielectronPair; |
61d106d3 | 32 | class TObjArray; |
b2a297fa | 33 | |
34 | class AliDielectronCF : public TNamed { | |
35 | public: | |
36 | enum {kNmaxAddSteps=50}; | |
37 | ||
38 | AliDielectronCF(); | |
39 | AliDielectronCF(const char* name, const char* title); | |
40 | virtual ~AliDielectronCF(); | |
41 | ||
e123f993 | 42 | void SetStepForMCtruth(Bool_t steps=kTRUE) { fStepForMCtruth=steps; } |
ba15fdfb | 43 | Bool_t GetStepForMCtruth() const { return fStepForMCtruth; } |
e123f993 | 44 | void SetStepForNoCutsMCmotherPid(Bool_t steps=kTRUE) { fStepForNoCutsMCmotherPid=steps; } |
45 | void SetStepForAfterAllCuts(Bool_t steps=kTRUE) { fStepForAfterAllCuts=steps; } | |
554e40f8 | 46 | void SetStepForPreFilter(Bool_t steps=kTRUE) { fStepForPreFilter=steps; } |
e123f993 | 47 | void SetStepsForEachCut(Bool_t steps=kTRUE) { fStepsForEachCut=steps; } |
48 | void SetStepsForCutsIncreasing(Bool_t steps=kTRUE) { fStepsForCutsIncreasing=steps; } | |
572b0139 | 49 | void SetStepsForSignal(Bool_t steps=kTRUE) { fStepsForSignal=steps; } |
50 | void SetStepsForBackground(Bool_t steps=kTRUE) { fStepsForBackground=steps; } | |
5720c765 | 51 | void SetStepsForMCtruthOnly(Bool_t steps=kTRUE) { fStepsForMCtruthOnly=steps; } |
b2a297fa | 52 | |
53 | void SetPdgMother(Int_t pdg) { fPdgMother=pdg; } | |
ba15fdfb | 54 | void SetSignalsMC(TObjArray* array) {fSignalsMC = array;} |
b2a297fa | 55 | |
56 | void AddStepMask(UInt_t mask) { fStepMasks[fNStepMasks++]=mask; } | |
57 | ||
3505bfad | 58 | void AddVariable(AliDielectronVarManager::ValueTypes type, Int_t nbins, |
59 | Double_t min, Double_t max, Bool_t leg=kFALSE, Bool_t log=kFALSE); | |
60 | void AddVariable(AliDielectronVarManager::ValueTypes type, const char* binLimitStr, Bool_t leg=kFALSE); | |
61 | void AddVariable(AliDielectronVarManager::ValueTypes type, TVectorD *binLimits, Bool_t leg=kFALSE); | |
61d106d3 | 62 | |
b2a297fa | 63 | void InitialiseContainer(const AliAnalysisFilter& filter); |
5720c765 | 64 | |
65 | Int_t GetNvarsPair() const {return fNVars;} | |
66 | Int_t GetNvarsLeg() const {return fNVarsLeg;} | |
67 | ||
ac390e40 | 68 | UInt_t GetVariablePair(UInt_t var) const {return (var>=(UInt_t)AliDielectronVarManager::kNMaxValues)? (UInt_t)AliDielectronVarManager::kNMaxValues+1:fVariables[var];} |
69 | UInt_t GetVariableLeg(UInt_t var) const {return (var>=(UInt_t)AliDielectronVarManager::kNMaxValues)? (UInt_t)AliDielectronVarManager::kNMaxValues+1:fVariablesLeg[var];} | |
5720c765 | 70 | |
6551594b | 71 | // void Fill(UInt_t mask, const TObject *particle); |
72 | void Fill(UInt_t mask, const AliDielectronPair *particle); | |
b2a297fa | 73 | void FillMC(const TObject *particle); |
ba15fdfb | 74 | void FillMC(Int_t label1, Int_t label2, Int_t nSignal); |
75 | ||
b2a297fa | 76 | AliCFContainer* GetContainer() const { return fCfContainer; } |
77 | ||
6551594b | 78 | private: |
b2a297fa | 79 | UInt_t fVariables[AliDielectronVarManager::kNMaxValues]; //configured variables |
6551594b | 80 | UInt_t fVariablesLeg[AliDielectronVarManager::kNMaxValues]; //configured variables for the legs |
b2a297fa | 81 | |
82 | Int_t fNSteps; // number of selection steps | |
6551594b | 83 | |
b2a297fa | 84 | Int_t fNVars; // number of variables |
3505bfad | 85 | TObjArray *fVarBinLimits; // array of bin limits |
86 | ||
87 | Int_t fNVarsLeg; // number of variables for the legs | |
88 | TObjArray *fVarBinLimitsLeg; // array of bin limits of the legs | |
6551594b | 89 | |
ac390e40 | 90 | Int_t fNCuts; // Number of cuts in the filter concerned |
b2a297fa | 91 | |
ac390e40 | 92 | Double_t *fValues; //! Value array for filling the container |
93 | Bool_t *fIsMCTruth; //! Buffer array for MC truth information | |
6551594b | 94 | |
e123f993 | 95 | Bool_t fStepForMCtruth; //create a step for the MC truth |
96 | Bool_t fStepForNoCutsMCmotherPid; //create a step for before cuts, but with MC truth of the mother | |
97 | Bool_t fStepForAfterAllCuts; //create a step for before cuts, but with MC truth of the mother | |
554e40f8 | 98 | Bool_t fStepForPreFilter; //create a step after PreFilter application |
e123f993 | 99 | Bool_t fStepsForEachCut; //create steps for each cut? |
100 | Bool_t fStepsForCutsIncreasing; //create steps for increasing cut combinatons? | |
101 | //e.g. cut1&cut2, cut1&cut2&cut3 ... | |
572b0139 | 102 | Bool_t fStepsForSignal; //steps for pure signal |
103 | Bool_t fStepsForBackground; //steps for pure background | |
5720c765 | 104 | Bool_t fStepsForMCtruthOnly; //Switch off all pair steps, allow only MC truth Class |
572b0139 | 105 | |
b2a297fa | 106 | UInt_t fStepMasks[kNmaxAddSteps]; //steps for additional cut combinatons |
107 | UInt_t fNStepMasks; //number of configured step masks | |
108 | ||
109 | Int_t fPdgMother; //Pdg code of MCtruth validation | |
ba15fdfb | 110 | TObjArray* fSignalsMC; //! array of MC signals to be studied |
b2a297fa | 111 | AliCFContainer* fCfContainer; //the CF container |
572b0139 | 112 | |
113 | Bool_t fHasMC; //if MC info is available | |
114 | Int_t fNAddSteps; //number of additional MC related steps per cut step | |
3505bfad | 115 | |
116 | TVectorD* MakeLogBinning(Int_t nbinsX, Double_t xmin, Double_t xmax) const; | |
117 | TVectorD* MakeLinBinning(Int_t nbinsX, Double_t xmin, Double_t xmax) const; | |
b2a297fa | 118 | |
119 | AliDielectronCF(const AliDielectronCF &c); | |
120 | AliDielectronCF &operator=(const AliDielectronCF &c); | |
121 | ||
3505bfad | 122 | ClassDef(AliDielectronCF,4) //Dielectron Correction Framework handler |
b2a297fa | 123 | }; |
124 | ||
125 | #endif |