1 #ifndef ALIANALYSISTASKSED0MASS_H
2 #define ALIANALYSISTASKSED0MASS_H
4 /* Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 //*************************************************************************
8 // Class AliAnalysisTaskSED0Mass
9 // AliAnalysisTaskSE for D0 candidates invariant mass histogram
10 // and comparison to MC truth (kinematics stored in the AOD) and cut variables
12 // Authors: A.Dainese, andrea.dainese@ln.infn.it
13 // and C.Bianchin, chiara.bianchin@pd.infn.it
14 //*************************************************************************
21 #include "AliAnalysisTaskSE.h"
22 #include "AliRDHFCutsD0toKpi.h"
26 class AliAnalysisTaskSED0Mass : public AliAnalysisTaskSE
30 AliAnalysisTaskSED0Mass();
31 AliAnalysisTaskSED0Mass(const char *name,AliRDHFCutsD0toKpi* cuts);
32 virtual ~AliAnalysisTaskSED0Mass();
35 // Implementation of interface methods
36 virtual void UserCreateOutputObjects();
38 virtual void LocalInit() {Init();}
39 virtual void UserExec(Option_t *option);
40 virtual void Terminate(Option_t *option);
42 void SetArray(Int_t type=AliAnalysisTaskSED0Mass::kD0){fArray=type;}
45 void SetReadMC(Bool_t readMC=kFALSE){fReadMC=readMC;}
46 void SetCutOnDistr(Bool_t cutondistr=kFALSE){fCutOnDistr=cutondistr;}
50 AliAnalysisTaskSED0Mass(const AliAnalysisTaskSED0Mass &source);
51 AliAnalysisTaskSED0Mass& operator=(const AliAnalysisTaskSED0Mass& source);
52 void FillMassHists(AliAODRecoDecayHF2Prong *part, TClonesArray *arrMC, AliRDHFCutsD0toKpi *cuts, TList *listout);
53 void FillVarHists(AliAODEvent *aodev,AliAODRecoDecayHF2Prong *part, TClonesArray *arrMC, AliRDHFCutsD0toKpi *cuts, TList *listout);
54 AliAODVertex* GetPrimaryVtxSkipped(AliAODEvent *aodev,AliAODRecoDecayHF2Prong *d);
56 TList *fOutputMass; //! list send on output slot 1
57 TList *fDistr; //! list send on output slot 2
58 TH1F *fNentries; //! histogram with number of events on output slot 3
59 TList *fChecks; //! list send on output slot 4
60 //TList* fCutList; //!
61 AliRDHFCutsD0toKpi *fCuts; // Cuts - sent to output slot 5
62 Int_t fArray; // can be D0 or Like Sign candidates
63 Bool_t fReadMC; // flag for MC array: kTRUE = read it, kFALSE = do not read it
64 Bool_t fCutOnDistr; // Flag to decide if apply cut also on distributions: 0 no cuts, 1 looser cuts, 2 tighter cuts
65 Int_t fNPtBins; // number of pt bins
66 Int_t* fTotPosPairs; //[fNPtBins]
67 Int_t* fTotNegPairs; //[fNPtBins]
68 Double_t fLsNormalization; // normalization
71 ClassDef(AliAnalysisTaskSED0Mass,7); // AliAnalysisTaskSE for the MC association of heavy-flavour decay candidates