2 // This class calculates the exclusive charged particle density
3 // in each for the 5 FMD rings.
5 #ifndef ALIFMDMCCORRECTOR_H
6 #define ALIFMDMCCORRECTOR_H
7 #include "AliFMDCorrector.h"
13 * @defgroup pwg2_forward_mc Monte-carlo code
15 * @ingroup pwg2_forward
18 * This class calculates the exclusive charged particle density
19 * in each for the 5 FMD rings.
22 * - 5 RingHistos objects - each with a number of vertex dependent
23 * 2D histograms of the inclusive charge particle density
26 * - 5 RingHistos objects - each with a number of vertex dependent
27 * 2D histograms of the exclusive charge particle density
29 * @par Corrections used:
30 * - AliFMDCorrSecondaryMap;
31 * - AliFMDCorrVertexBias
32 * - AliFMDCorrMergingEfficiency
34 * @ingroup pwg2_forward_algo
35 * @ingroup pwg2_forward_mc
37 class AliFMDMCCorrector : public AliFMDCorrector
55 * @param name Name of object
57 AliFMDMCCorrector(const char* name)
58 : AliFMDCorrector(name),
69 * @param o Object to copy from
71 AliFMDMCCorrector(const AliFMDMCCorrector& o)
83 virtual ~AliFMDMCCorrector();
85 * Assignement operator
87 * @param o Object to assign from
89 * @return Reference to this object
91 AliFMDMCCorrector& operator=(const AliFMDMCCorrector&);
93 * Initialize this object
95 * @param etaAxis Eta axis to use
97 void Init(const TAxis& etaAxis);
101 * @param hists Cache of histograms
102 * @param vtxBin Vertex bin
104 * @return true on successs
106 virtual Bool_t CorrectMC(AliForwardUtil::Histos& hists, UShort_t vtxBin);
108 * Compare the result of analysing the ESD for
109 * the inclusive charged particle density to analysing
117 virtual Bool_t CompareResults(AliForwardUtil::Histos& esd,
118 AliForwardUtil::Histos& mc);
120 * Output diagnostic histograms to directory
122 * @param dir List to write in
124 void DefineOutput(TList* dir);
127 * MAke comparison profiles
131 * @param axis Eta axis
133 * @return Newly allocated profile object
135 TProfile2D* Make(UShort_t d, Char_t r, const TAxis& axis) const;
137 * Fill comparison profiles
141 * @param esd ESD histogram
142 * @param mc MC histogram
144 void Fill(UShort_t d, Char_t r, TH2* esd, TH2* mc);
146 TProfile2D* fFMD1i; // Comparison
147 TProfile2D* fFMD2i; // Comparison
148 TProfile2D* fFMD2o; // Comparison
149 TProfile2D* fFMD3i; // Comparison
150 TProfile2D* fFMD3o; // Comparison
151 TList* fComps; // List of comparisons
153 ClassDef(AliFMDMCCorrector,1); // Calculate Nch density