2 // Task to analyse the AOD for for dN/deta in the forward regions
4 #ifndef ALIMCTRUTHDNDETATASK_H
5 #define ALIMCTRUTHDNDETATASK_H
7 * @file AliMCTruthdNdetaTask.h
8 * @author Christian Holm Christensen <cholm@dalsgaard.hehi.nbi.dk>
9 * @date Wed Mar 23 14:04:54 2011
14 * @ingroup pwglf_forward_dndeta
16 #include "AliBasedNdetaTask.h"
22 * Tasks to determine @f$ dN/d\eta@f$ in the forward regions
24 * @ingroup pwglf_forward_tasks_dndeta
25 * @ingroup pwglf_forward_dndeta
27 class AliMCTruthdNdetaTask : public AliBasedNdetaTask
34 AliMCTruthdNdetaTask();
38 * @param name Name of task
40 AliMCTruthdNdetaTask(const char* name);
45 virtual ~AliMCTruthdNdetaTask() {}
47 * Called at end of event processing.
49 * This is called once in the master
51 * @param option Not used
53 virtual void Terminate(Option_t* option);
58 * @param o Object to copy from
60 AliMCTruthdNdetaTask(const AliMCTruthdNdetaTask& o);
62 * Assigmement operator
64 * @return Reference to this
66 AliMCTruthdNdetaTask& operator=(const AliMCTruthdNdetaTask&) { return *this; }
69 * Retrieve the histogram
71 * @param aod AOD event
72 * @param mc Whether to get the MC histogram or not
74 * @return Retrieved histogram or null
76 TH2D* GetHistogram(const AliAODEvent* aod, Bool_t mc);
78 * Get the marker style
80 * @return Marker style
82 virtual Int_t GetMarker() const { return GetMarkerStyle(kStar); }
84 * Get the colour to use for markers (only pp - in PbPb we use a rainbow)
86 * @return Marker colour
88 virtual Int_t GetColor() const { return kGray+2; }
90 * Make a new centrality bin
92 * @param name Histogram names
96 * @return Newly allocated object (of our type)
98 AliBasedNdetaTask::CentralityBin*
99 MakeCentralityBin(const char* name, Short_t l, Short_t h) const;
102 * Class that holds data for a single centrality bin
105 class CentralityBin : public AliBasedNdetaTask::CentralityBin
111 CentralityBin() : AliBasedNdetaTask::CentralityBin(), fSumTruth(0) {}
115 * @param name Name used for histograms (e.g., Forward)
116 * @param low Lower centrality cut in percent
117 * @param high Upper centrality cut in percent
119 CentralityBin(const char* name, Short_t low, Short_t high)
120 : AliBasedNdetaTask::CentralityBin(name, low, high),
126 * @param other Object to copy from
128 CentralityBin(const CentralityBin& other)
129 : AliBasedNdetaTask::CentralityBin(other),
130 fSumTruth(other.fSumTruth)
135 virtual ~CentralityBin() {}
137 * Assignement operator
142 CentralityBin& operator=(const CentralityBin&) { return *this; }
146 * @param forward Forward data (for trigger, vertex, & centrality)
147 * @param triggerMask Trigger mask
148 * @param isZero True if this is a zero bin event
149 * @param vzMin Minimum IP z coordinate
150 * @param vzMax Maximum IP z coordinate
151 * @param data Data histogram
152 * @param mc MC histogram
154 virtual void ProcessEvent(const AliAODForwardMult* forward,
164 * @param sums List of sums
165 * @param results Output list of results
166 * @param scheme Normalisation scheme options
167 * @param shapeCorr Shape correction or nil
168 * @param trigEff Trigger efficiency
169 * @param symmetrice Whether to symmetrice the results
170 * @param rebin Whether to rebin the results
171 * @param rootProj If true, use TH2::ProjectionX
172 * @param corrEmpty Whether to correct for empty bins
173 * @param cutEdges Whether to cut edges when rebinning
174 * @param triggerMask Trigger mask
175 * @param color Marker colour
176 * @param marker Marker style
177 * @param mclist List of MC results
178 * @param truthlist List of MC truth results
180 virtual void End(TList* sums,
183 const TH2F* shapeCorr,
196 TH2D* fSumTruth; // Sum of primary histograms
197 ClassDef(CentralityBin,1); // A centrality bin
200 ClassDef(AliMCTruthdNdetaTask,1); // Determine multiplicity in forward region