1 #ifndef ALIANALYSISTASKCHECKHFMCPROD
2 #define ALIANALYSISTASKCHECKHFMCPROD
4 /* Copyright(c) 1998-2012, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
9 //*************************************************************************
10 // Class AliAnalysisTaskCheckHFMCProd
11 // AliAnalysisTask to check MC production at ESD+Kine level
14 // Author: F. Prino, prino@to.infn.it
16 //*************************************************************************
28 #include "AliAnalysisTaskSE.h"
30 class AliAnalysisTaskCheckHFMCProd : public AliAnalysisTaskSE {
34 AliAnalysisTaskCheckHFMCProd();
35 virtual ~AliAnalysisTaskCheckHFMCProd();
37 virtual void UserExec(Option_t *option);
38 virtual void UserCreateOutputObjects();
39 virtual void Terminate(Option_t *option);
41 Int_t CheckD0Decay(Int_t labD0, AliStack* stack) const;
42 Int_t CheckDplusDecay(Int_t labDplus, AliStack* stack) const;
43 Int_t CheckDsDecay(Int_t labDs, AliStack* stack) const;
44 Int_t CheckDstarDecay(Int_t labDstar, AliStack* stack) const;
45 Int_t CheckLcDecay(Int_t labLc, AliStack* stack) const;
46 void SetSearchUpToQuark(Bool_t opt){fSearchUpToQuark=opt;};
49 void SetReadMC(Bool_t opt) {fReadMC=opt;}
50 void SetpPb() {fSystem=2;}
51 void SetPbPb() {fSystem=1;}
52 void Setpp() {fSystem=0;}
56 AliAnalysisTaskCheckHFMCProd(const AliAnalysisTaskCheckHFMCProd &source);
57 AliAnalysisTaskCheckHFMCProd& operator=(const AliAnalysisTaskCheckHFMCProd &source);
59 TList* fOutput; //! list of output histos
60 TH1F* fHistoNEvents; //! histo with N of events
62 TH1F* fHistoPhysPrim; //! histo of n. of physical primaries in |eta|<0.5
63 TH1F* fHistoTracks; //! histo with number of ESD tracks
64 TH1F* fHistoSelTracks; //! histo with number of TPC+ITS refit ESD tracks
65 TH1F* fHistoTracklets; //! histo with number of SPD tracklets
66 TH1F* fHistoPtPhysPrim; //! histo of pt distribution of phys primaries
67 TH1F* fHistoEtaPhysPrim; //! histo of eta distribution of physical primaries
69 TH1F* fHistoSPD3DVtxX; //! histo with distr. of x coord. of SPD3D vertex
70 TH1F* fHistoSPD3DVtxY; //! histo with distr. of y coord. of SPD3D vertex
71 TH1F* fHistoSPD3DVtxZ; //! histo with distr. of z coord. of SPD3D vertex
73 TH1F* fHistoSPDZVtxX; //! histo with distr. of x coord. of SPDZ vertex
74 TH1F* fHistoSPDZVtxY; //! histo with distr. of y coord. of SPDZ vertex
75 TH1F* fHistoSPDZVtxZ; //! histo with distr. of z coord. of SPDZ vertex
77 TH1F* fHistoTRKVtxX; //! histo with distr. of x coord. of TRK vertex
78 TH1F* fHistoTRKVtxY; //! histo with distr. of y coord. of TRK vertex
79 TH1F* fHistoTRKVtxZ; //! histo with distr. of z coord. of TRK vertex
81 TH2F* fHistoNcharmed; //! histo of D mesons vs. dN/dy
82 TH2F* fHistoNbVsNc; //! histo of n. b quarks vs. n c. quarks
84 TH2F* fHistBYPtAllDecay[5]; //! histo of y vs. pt from prompt B0, B+, B*, Bs, Lb
85 TH2F* fHistYPtAllDecay[5]; //! histo of y vs. pt from prompt D0, D+, D*, Ds, Lc, no selection on decay channel
86 TH2F* fHistYPtPromptAllDecay[5]; //! histo of y vs. pt from prompt D0, D+, D*, Ds, Lc, no selection on decay channel
87 TH2F* fHistYPtFeeddownAllDecay[5]; //! histo of y vs. pt from prompt D0, D+, D*, Ds, Lc, no selection on decay channel
88 TH2F* fHistYPtPrompt[5]; //! histo of y vs. pt from prompt D0, D+, D*, Ds, Lc
89 TH2F* fHistYPtFeeddown[5]; //! histo of y vs. pt from feeddown D0, D+, D*, Ds, Lc
90 TH2F* fHistYPtD0byDecChannel[2]; //! histo of y vs. pt for D0->Kpi and D0->Kpipipi
91 TH2F* fHistYPtDplusbyDecChannel[2]; //! histo of y vs. pt for D+->Kpipi and D+->K0*pi
92 TH2F* fHistYPtDsbyDecChannel[2]; //! histo of y vs. pt for Ds->phipi and Ds->K0*K
93 TH1F* fHistOriginPrompt; //! histo of D production point (prompt)
94 TH1F* fHistOriginFeeddown; //! histo of D production point (feeddown)
95 TH1F* fHistMotherID; //! histo of mother ID
96 TH1F* fHistDSpecies; //! histo of D hadron species
97 TH1F* fHistBSpecies; //! histo of B hadron species
98 TH2F* fHistNcollHFtype; //! histo of B hadron species
99 Bool_t fSearchUpToQuark; // c/b separation using quarks
100 Int_t fSystem; // 0=pp, 1=PbPb, 2=pPb
103 ClassDef(AliAnalysisTaskCheckHFMCProd,4);