]>
Commit | Line | Data |
---|---|---|
f5ec6c30 | 1 | #ifndef ALICFHEAVYFLAVOURTASKMULTIVARMULTISTEP_H |
2 | #define ALICFHEAVYFLAVOURTASKMULTIVARMULTISTEP_H | |
3 | /************************************************************************** | |
4 | * Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. * | |
5 | * * | |
6 | * Author: The ALICE Off-line Project. * | |
7 | * Contributors are mentioned in the code where appropriate. * | |
8 | * * | |
9 | * Permission to use, copy, modify and distribute this software and its * | |
10 | * documentation strictly for non-commercial purposes is hereby granted * | |
11 | * without fee, provided that the above copyright notice appears in all * | |
12 | * copies and that both the copyright notice and this permission notice * | |
13 | * appear in the supporting documentation. The authors make no claims * | |
14 | * about the suitability of this software for any purpose. It is * | |
15 | * provided "as is" without express or implied warranty. * | |
16 | **************************************************************************/ | |
17 | ||
18 | //----------------------------------------------------------------------- | |
19 | // Class for HF corrections as a function of many variables and step | |
20 | // Author : C. Zampolli, CERN | |
a5d92cb4 | 21 | // Base class for HF Unfolding - agrelli@uu.nl |
f5ec6c30 | 22 | //----------------------------------------------------------------------- |
23 | ||
24 | ||
25 | #include "AliAnalysisTaskSE.h" | |
26 | ||
27 | class TH1I; | |
28 | class TParticle ; | |
29 | class TFile ; | |
30 | class TClonesArray ; | |
31 | class AliCFManager; | |
32 | class AliAODRecoDecay; | |
33 | class AliAODRecoDecayHF2Prong; | |
34 | class AliAODMCParticle; | |
a5d92cb4 | 35 | class THnSparse; |
174db2b5 | 36 | class AliRDHFCutsD0toKpi; |
f5ec6c30 | 37 | |
38 | class AliCFHeavyFlavourTaskMultiVarMultiStep : public AliAnalysisTaskSE { | |
39 | public: | |
40 | ||
41 | enum { | |
174db2b5 | 42 | kStepGeneratedLimAcc = 0, |
43 | kStepGenerated = 1, | |
44 | kStepAcceptance = 2, | |
45 | kStepVertex = 3, | |
46 | kStepRefit = 4, | |
47 | kStepReconstructed = 5, | |
48 | kStepRecoAcceptance = 6, | |
49 | kStepRecoITSClusters = 7, | |
50 | kStepRecoPPR = 8, | |
51 | kStepRecoPID = 9 | |
f5ec6c30 | 52 | }; |
53 | ||
54 | AliCFHeavyFlavourTaskMultiVarMultiStep(); | |
174db2b5 | 55 | AliCFHeavyFlavourTaskMultiVarMultiStep(const Char_t* name, AliRDHFCutsD0toKpi* cuts); |
f5ec6c30 | 56 | AliCFHeavyFlavourTaskMultiVarMultiStep& operator= (const AliCFHeavyFlavourTaskMultiVarMultiStep& c); |
57 | AliCFHeavyFlavourTaskMultiVarMultiStep(const AliCFHeavyFlavourTaskMultiVarMultiStep& c); | |
58 | virtual ~AliCFHeavyFlavourTaskMultiVarMultiStep(); | |
59 | ||
60 | // ANALYSIS FRAMEWORK STUFF to loop on data and fill output objects | |
61 | void UserCreateOutputObjects(); | |
62 | void UserExec(Option_t *option); | |
63 | void Terminate(Option_t *); | |
a5d92cb4 | 64 | |
65 | // UNFOLDING | |
66 | void SetCorrelationMatrix(THnSparse* h) {fCorrelation=h;} | |
67 | void SetAcceptanceUnf(Bool_t AcceptanceUnf) {fAcceptanceUnf = AcceptanceUnf;} | |
68 | Bool_t GetAcceptanceUnf() const {return fAcceptanceUnf;} | |
69 | ||
f5ec6c30 | 70 | |
71 | // CORRECTION FRAMEWORK RELATED FUNCTIONS | |
72 | void SetCFManager(AliCFManager* io) {fCFManager = io;} // global correction manager | |
73 | AliCFManager * GetCFManager() {return fCFManager;} // get corr manager | |
74 | ||
75 | void SetPDG(Int_t code) {fPDG = code; } // defines the PDG code of searched HF | |
76 | Double_t CosThetaStar(AliAODMCParticle* mcPart, AliAODMCParticle* mcPartDaughter0, AliAODMCParticle* mcPartDaughter1) const; // returns cos(ThetaStar) of the D0 decay | |
77 | Double_t CT(AliAODMCParticle* mcPart, AliAODMCParticle* mcPartDaughter0, AliAODMCParticle* mcPartDaughter1) const; // returns cT of the D0 decay | |
78 | void SetFillFromGenerated(Bool_t flag) {fFillFromGenerated = flag;} | |
79 | Bool_t GetFillFromGenerated() const {return fFillFromGenerated;} | |
80 | Bool_t GetGeneratedValuesFromMCParticle(AliAODMCParticle* mcPart, TClonesArray* mcArray, Double_t* vectorMC) const; | |
81 | void SetMinITSClusters(Int_t minITSClusters) {fMinITSClusters = minITSClusters;} | |
82 | Int_t GetMinITSClusters() const {return fMinITSClusters;} | |
e8158e56 | 83 | Int_t CheckOrigin(AliAODMCParticle* mcPart, TClonesArray* mcArray) const; |
f5ec6c30 | 84 | |
f6d06004 | 85 | void SetKeepD0fromB(Bool_t keepD0fromB){fKeepD0fromB=keepD0fromB;} |
86 | ||
f5ec6c30 | 87 | protected: |
88 | Int_t fPDG; // PDG code of searched V0's | |
4df1e470 | 89 | AliCFManager *fCFManager; // pointer to the CF manager |
f5ec6c30 | 90 | TH1I *fHistEventsProcessed; //! simple histo for monitoring the number of events processed |
a5d92cb4 | 91 | THnSparse* fCorrelation; // response matrix for unfolding |
f5ec6c30 | 92 | Int_t fCountMC; // MC particle found |
93 | Int_t fCountAcc; // MC particle found that satisfy acceptance cuts | |
ae3e319c | 94 | Int_t fCountVertex; // Reco particle found that satisfy vertex constrained |
95 | Int_t fCountRefit; // Reco particle found that satisfy kTPCrefit and kITSrefit | |
f5ec6c30 | 96 | Int_t fCountReco; // Reco particle found that satisfy cuts |
97 | Int_t fCountRecoAcc; // Reco particle found that satisfy cuts in requested acceptance | |
98 | Int_t fCountRecoITSClusters; // Reco particle found that satisfy cuts in n. of ITS clusters | |
99 | Int_t fCountRecoPPR; // Reco particle found that satisfy cuts in PPR | |
174db2b5 | 100 | Int_t fCountRecoPID; // Reco particle found that satisfy cuts in PPR, and PID condition |
f5ec6c30 | 101 | Int_t fEvents; // n. of events |
102 | Bool_t fFillFromGenerated; // flag to indicate whether data container should be filled | |
103 | // with generated values also for reconstructed particles | |
104 | Int_t fMinITSClusters; // min n. of ITS clusters for RecoDecay | |
a5d92cb4 | 105 | Bool_t fAcceptanceUnf; // flag for unfolding before or after cuts. |
f6d06004 | 106 | |
107 | Bool_t fKeepD0fromB; // flag to consider also D0 coming from B | |
174db2b5 | 108 | AliRDHFCutsD0toKpi* fCuts; // cuts |
f5ec6c30 | 109 | |
174db2b5 | 110 | ClassDef(AliCFHeavyFlavourTaskMultiVarMultiStep,4); // class for HF corrections as a function of many variables |
f5ec6c30 | 111 | }; |
112 | ||
113 | #endif |