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;
47 void SetReadMC(Bool_t opt) {fReadMC=opt;}
48 void SetPbPb() {fPbPb=kTRUE;}
49 void Setpp() {fPbPb=kFALSE;}
53 AliAnalysisTaskCheckHFMCProd(const AliAnalysisTaskCheckHFMCProd &source);
54 AliAnalysisTaskCheckHFMCProd& operator=(const AliAnalysisTaskCheckHFMCProd &source);
56 TList* fOutput; //! list of output histos
57 TH1F* fHistoNEvents; //! histo with N of events
59 TH1F* fHistoTracks; //! histo with number of ESD tracks
60 TH1F* fHistoSelTracks; //! histo with number of TPC+ITS refit ESD tracks
61 TH1F* fHistoTracklets; //! histo with number of SPD tracklets
63 TH1F* fHistoSPD3DVtxX; //! histo with distr. of x coord. of SPD3D vertex
64 TH1F* fHistoSPD3DVtxY; //! histo with distr. of y coord. of SPD3D vertex
65 TH1F* fHistoSPD3DVtxZ; //! histo with distr. of z coord. of SPD3D vertex
67 TH1F* fHistoSPDZVtxX; //! histo with distr. of x coord. of SPDZ vertex
68 TH1F* fHistoSPDZVtxY; //! histo with distr. of y coord. of SPDZ vertex
69 TH1F* fHistoSPDZVtxZ; //! histo with distr. of z coord. of SPDZ vertex
71 TH1F* fHistoTRKVtxX; //! histo with distr. of x coord. of TRK vertex
72 TH1F* fHistoTRKVtxY; //! histo with distr. of y coord. of TRK vertex
73 TH1F* fHistoTRKVtxZ; //! histo with distr. of z coord. of TRK vertex
75 TH2F* fHistoNcharmed; //! histo of D mesons vs. dN/dy
76 TH2F* fHistoNbVsNc; //! histo of n. b quarks vs. n c. quarks
78 TH2F* fHistYPtAllDecay[5]; //! histo of y vs. pt from prompt D0, D+, D*, Ds, Lc, no selection on decay channel
79 TH2F* fHistYPtPromptAllDecay[5]; //! histo of y vs. pt from prompt D0, D+, D*, Ds, Lc, no selection on decay channel
80 TH2F* fHistYPtFeeddownAllDecay[5]; //! histo of y vs. pt from prompt D0, D+, D*, Ds, Lc, no selection on decay channel
81 TH2F* fHistYPtPrompt[5]; //! histo of y vs. pt from prompt D0, D+, D*, Ds, Lc
82 TH2F* fHistYPtFeeddown[5]; //! histo of y vs. pt from feeddown D0, D+, D*, Ds, Lc
83 TH2F* fHistYPtD0byDecChannel[2]; //! histo of y vs. pt for D0->Kpi and D0->Kpipipi
84 TH2F* fHistYPtDplusbyDecChannel[2]; //! histo of y vs. pt for D+->Kpipi and D+->K0*pi
85 TH2F* fHistYPtDsbyDecChannel[2]; //! histo of y vs. pt for Ds->phipi and Ds->K0*K
90 ClassDef(AliAnalysisTaskCheckHFMCProd,1);