Increase version number
[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AliCFHeavyFlavourTaskMultiVarMultiStep.h
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
21 // Base class for HF Unfolding - agrelli@uu.nl
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;
35 class THnSparse;
36
37 class AliCFHeavyFlavourTaskMultiVarMultiStep : public AliAnalysisTaskSE {
38   public:
39
40   enum {
41     kStepGenerated       = 0,
42     kStepAcceptance      = 1,
43     kStepVertex          = 2,
44     kStepRefit           = 3,
45     kStepReconstructed   = 4,
46     kStepRecoAcceptance  = 5,
47     kStepRecoITSClusters = 6,
48     kStepRecoPPR         = 7
49   };
50
51   AliCFHeavyFlavourTaskMultiVarMultiStep();
52   AliCFHeavyFlavourTaskMultiVarMultiStep(const Char_t* name);
53   AliCFHeavyFlavourTaskMultiVarMultiStep& operator= (const AliCFHeavyFlavourTaskMultiVarMultiStep& c);
54   AliCFHeavyFlavourTaskMultiVarMultiStep(const AliCFHeavyFlavourTaskMultiVarMultiStep& c);
55   virtual ~AliCFHeavyFlavourTaskMultiVarMultiStep();
56
57   // ANALYSIS FRAMEWORK STUFF to loop on data and fill output objects
58   void     UserCreateOutputObjects();
59   void     UserExec(Option_t *option);
60   void     Terminate(Option_t *);
61
62  // UNFOLDING
63   void     SetCorrelationMatrix(THnSparse* h) {fCorrelation=h;}
64   void     SetAcceptanceUnf(Bool_t AcceptanceUnf) {fAcceptanceUnf = AcceptanceUnf;}
65   Bool_t   GetAcceptanceUnf() const {return fAcceptanceUnf;}
66
67   
68   // CORRECTION FRAMEWORK RELATED FUNCTIONS
69   void           SetCFManager(AliCFManager* io) {fCFManager = io;}   // global correction manager
70   AliCFManager * GetCFManager()                 {return fCFManager;} // get corr manager
71
72   void     SetPDG(Int_t code) {fPDG = code; }     // defines the PDG code of searched HF
73   Double_t CosThetaStar(AliAODMCParticle* mcPart, AliAODMCParticle* mcPartDaughter0, AliAODMCParticle* mcPartDaughter1) const;  // returns cos(ThetaStar) of the D0 decay
74   Double_t CT(AliAODMCParticle* mcPart, AliAODMCParticle* mcPartDaughter0, AliAODMCParticle* mcPartDaughter1) const;            // returns cT of the D0 decay
75   void     SetFillFromGenerated(Bool_t flag) {fFillFromGenerated = flag;}
76   Bool_t   GetFillFromGenerated() const {return fFillFromGenerated;}
77   Bool_t   GetGeneratedValuesFromMCParticle(AliAODMCParticle* mcPart, TClonesArray* mcArray, Double_t* vectorMC) const;
78   void     SetMinITSClusters(Int_t minITSClusters) {fMinITSClusters = minITSClusters;}
79   Int_t    GetMinITSClusters() const {return fMinITSClusters;}
80   Int_t    CheckOrigin(AliAODMCParticle* mcPart, TClonesArray* mcArray) const;
81
82  protected:
83   Int_t           fPDG;         //  PDG code of searched V0's
84   AliCFManager   *fCFManager;   //  pointer to the CF manager
85   TH1I *fHistEventsProcessed;   //! simple histo for monitoring the number of events processed
86   THnSparse* fCorrelation;      //  response matrix for unfolding
87   Int_t fCountMC;               //  MC particle found
88   Int_t fCountAcc;              //  MC particle found that satisfy acceptance cuts
89   Int_t fCountVertex;       //  Reco particle found that satisfy vertex constrained
90   Int_t fCountRefit;        //  Reco particle found that satisfy kTPCrefit and kITSrefit
91   Int_t fCountReco;             //  Reco particle found that satisfy cuts
92   Int_t fCountRecoAcc;          //  Reco particle found that satisfy cuts in requested acceptance
93   Int_t fCountRecoITSClusters;  //  Reco particle found that satisfy cuts in n. of ITS clusters
94   Int_t fCountRecoPPR;          //  Reco particle found that satisfy cuts in PPR
95   Int_t fEvents;                //  n. of events
96   Bool_t fFillFromGenerated;    //  flag to indicate whether data container should be filled 
97                                 //  with generated values also for reconstructed particles
98   Int_t fMinITSClusters;        //  min n. of ITS clusters for RecoDecay
99   Bool_t fAcceptanceUnf;        //  flag for unfolding before or after cuts.
100   
101   ClassDef(AliCFHeavyFlavourTaskMultiVarMultiStep,3); // class for HF corrections as a function of many variables
102 };
103
104 #endif