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 */
8 //#############################################################
10 //# Class AliDielectronCF #
11 //# Dielectron Correction Framework Manager #
14 //# Anton Andronic, GSI / A.Andronic@gsi.de #
15 //# Ionut C. Arsene, GSI / I.C.Arsene@gsi.de #
16 //# Julian Book, Uni Ffm / Julian.Book@cern.ch #
17 //# Frederick Kramer, Uni Ffm, / Frederick.Kramer@cern.ch #
18 //# Magnus Mager, CERN / Magnus.Mager@cern.ch #
19 //# WooJin J. Park, GSI / W.J.Park@gsi.de #
20 //# Jens Wiechula, Uni HD / Jens.Wiechula@cern.ch #
22 //#############################################################
27 #include <TVectorDfwd.h>
28 #include "AliDielectronVarManager.h"
30 class AliAnalysisCuts;
31 class AliAnalysisFilter;
33 class AliDielectronPair;
36 class AliDielectronCF : public TNamed {
38 enum {kNmaxAddSteps=50};
41 AliDielectronCF(const char* name, const char* title);
42 virtual ~AliDielectronCF();
44 void SetStepForMCtruth(Bool_t steps=kTRUE) { fStepForMCtruth=steps; }
45 void SetStepForNoCutsMCmotherPid(Bool_t steps=kTRUE) { fStepForNoCutsMCmotherPid=steps; }
46 void SetStepForAfterAllCuts(Bool_t steps=kTRUE) { fStepForAfterAllCuts=steps; }
47 void SetStepForPreFilter(Bool_t steps=kTRUE) { fStepForPreFilter=steps; }
48 void SetStepsForEachCut(Bool_t steps=kTRUE) { fStepsForEachCut=steps; }
49 void SetStepsForCutsIncreasing(Bool_t steps=kTRUE) { fStepsForCutsIncreasing=steps; }
50 void SetStepsForSignal(Bool_t steps=kTRUE) { fStepsForSignal=steps; }
51 void SetStepsForBackground(Bool_t steps=kTRUE) { fStepsForBackground=steps; }
53 void SetPdgMother(Int_t pdg) { fPdgMother=pdg; }
55 void AddStepMask(UInt_t mask) { fStepMasks[fNStepMasks++]=mask; }
57 void AddVariable(AliDielectronVarManager::ValueTypes type, Int_t nbins,
58 Double_t min, Double_t max, Bool_t leg=kFALSE, Bool_t log=kFALSE);
59 void AddVariable(AliDielectronVarManager::ValueTypes type, const char* binLimitStr, Bool_t leg=kFALSE);
60 void AddVariable(AliDielectronVarManager::ValueTypes type, TVectorD *binLimits, Bool_t leg=kFALSE);
62 void InitialiseContainer(const AliAnalysisFilter& filter);
64 // void Fill(UInt_t mask, const TObject *particle);
65 void Fill(UInt_t mask, const AliDielectronPair *particle);
66 void FillMC(const TObject *particle);
68 AliCFContainer* GetContainer() const { return fCfContainer; }
71 UInt_t fVariables[AliDielectronVarManager::kNMaxValues]; //configured variables
72 UInt_t fVariablesLeg[AliDielectronVarManager::kNMaxValues]; //configured variables for the legs
74 Int_t fNSteps; // number of selection steps
76 Int_t fNVars; // number of variables
77 TObjArray *fVarBinLimits; // array of bin limits
79 Int_t fNVarsLeg; // number of variables for the legs
80 TObjArray *fVarBinLimitsLeg; // array of bin limits of the legs
82 Int_t fNCuts; // Number of cuts in the filter concerned
84 Double_t *fValues; //! Value array for filling the container
86 Bool_t fStepForMCtruth; //create a step for the MC truth
87 Bool_t fStepForNoCutsMCmotherPid; //create a step for before cuts, but with MC truth of the mother
88 Bool_t fStepForAfterAllCuts; //create a step for before cuts, but with MC truth of the mother
89 Bool_t fStepForPreFilter; //create a step after PreFilter application
90 Bool_t fStepsForEachCut; //create steps for each cut?
91 Bool_t fStepsForCutsIncreasing; //create steps for increasing cut combinatons?
92 //e.g. cut1&cut2, cut1&cut2&cut3 ...
93 Bool_t fStepsForSignal; //steps for pure signal
94 Bool_t fStepsForBackground; //steps for pure background
96 UInt_t fStepMasks[kNmaxAddSteps]; //steps for additional cut combinatons
97 UInt_t fNStepMasks; //number of configured step masks
99 Int_t fPdgMother; //Pdg code of MCtruth validation
100 AliCFContainer* fCfContainer; //the CF container
102 Bool_t fHasMC; //if MC info is available
103 Int_t fNAddSteps; //number of additional MC related steps per cut step
105 TVectorD* MakeLogBinning(Int_t nbinsX, Double_t xmin, Double_t xmax) const;
106 TVectorD* MakeLinBinning(Int_t nbinsX, Double_t xmin, Double_t xmax) const;
108 AliDielectronCF(const AliDielectronCF &c);
109 AliDielectronCF &operator=(const AliDielectronCF &c);
111 ClassDef(AliDielectronCF,4) //Dielectron Correction Framework handler