// // This class contains the secondary correction and the double hit // correction used in low-flux events. // #ifndef ALIFMDCORRSECONDARYMAP_H #define ALIFMDCORRSECONDARYMAP_H #include #include #include class TH2D; /** * This class contains the secondary correction. * * The secondary correction is given by * @f[ * c_{r,v}(\eta,\varphi) = * \frac{\sum_i N_{ch,i,v,\mbox{primary}}(\eta,\varphi)}{ * \sum_i N_{ch,i,r,v,\mbox{FMD}}(\eta,\varphi)} * @f] * where @f$N_{ch,i,v,\mbox{primary}}(\eta,\varphi)@f$ is the * is the number of primary charged particles that fall within * the @f$(\eta,\varphi)@f$ bin in event @f$i@f$ with vertex @f$v@f$, * and is the total (primary and secondary) charged particles * that hit ring @f$r@f$ within @f$(\eta,\varphi)@f$ bin in event * @f$i@f$ with vertex @f$v@f$. * * These are generated from Monte-Carlo truth information. * * @ingroup pwg2_forward_corr */ class AliFMDCorrSecondaryMap : public TObject { public: /** * Default constructor */ AliFMDCorrSecondaryMap(); /** * Copy constructor * * @param o Object to copy from */ AliFMDCorrSecondaryMap(const AliFMDCorrSecondaryMap& o); /** * Destructor * */ virtual ~AliFMDCorrSecondaryMap(); /** * @{ * @name Get corrections and parameters */ /** * Assignment operator * * @param o Object to assign from * * @return Reference to this object */ AliFMDCorrSecondaryMap& operator=(const AliFMDCorrSecondaryMap& o); /** * Get the secondary correction @f$ c_{r,v}@f$ * * @param d Detector number (1-3) * @param r Ring identifier (I or O) * @param v Primary interaction point @f$z@f$ coordinate * * @return The correction @f$ c_{r,v}@f$ */ TH2D* GetCorrection(UShort_t d, Char_t r, Double_t v) const; /** * Get the secondary correction @f$ c_{r,v}@f$ * * @param d Detector number (1-3) * @param r Ring identifier (I or O) * @param b Bin corresponding to the primary interaction point * @f$z@f$ coordinate (1 based) * * @return The correction @f$ c_{r,v}@f$ */ TH2D* GetCorrection(UShort_t d, Char_t r, UShort_t b) const; /** * Get the vertex axis used * * @return vertex axis */ const TAxis& GetVertexAxis() const { return fVertexAxis; } /** * Get the eta axis used * * @return eta axis */ const TAxis& GetEtaAxis() const { return fEtaAxis; } /* @} */ /** * @{ * @name Set corrections and parameters */ /** * Set the secondary map correction @f$ c_{r,v}(\eta,\varphi)@f$. * Note, that the object takes ownership of the passed pointer. * * @param d Detector number (1-3) * @param r Ring identifier (I or O) * @param v Primary interaction point @f$z@f$ coordinate * @param h @f$ c_{r,v}(\eta,\varphi)@f$ * * @return true if operation succeeded */ Bool_t SetCorrection(UShort_t d, Char_t r, Double_t v, TH2D* h); /** * Set the secondary map correction @f$ c_{r,v}(\eta,\varphi)@f$ * Note, that the object takes ownership of the passed pointer. * * @param d Detector number (1-3) * @param r Ring identifier (I or O) * @param b Bin corresponding to the primary interaction point * @f$z@f$ coordinate (1 based) * @param h @f$ c_{r,v}(\eta,\varphi)@f$ * * @return true if operation succeeded */ Bool_t SetCorrection(UShort_t d, Char_t r, UShort_t b, TH2D* h); /** * Set the vertex axis to use * * @param axis Vertex axis */ void SetVertexAxis(const TAxis& axis); /** * Set the vertex axis to use * * @param nBins Number of bins * @param min Minimum * @param max Maximum */ void SetVertexAxis(Int_t nBins, Double_t min, Double_t max); /** * Set the eta axis to use * * @param axis Eta axis */ void SetEtaAxis(const TAxis& axis); /** * Set the eta axis to use * * @param nBins Number of bins * @param min Minimum * @param max Maximum */ void SetEtaAxis(Int_t nBins, Double_t min, Double_t max); /* @} */ /** * @{ * @name Auxiliary member functions */ /** * Declare this as a folder * * @return Always true */ Bool_t IsFolder() const { return true; } /** * Browse this object in the browser * * @param b */ void Browse(TBrowser* b); /** * Print this object * * @param option */ void Print(Option_t* option="R") const; //*MENU* /* @} */ protected: /** * Find the vertex bin that corresponds to the passed vertex * * @param vertex The interaction points @f$z@f$-coordinate * * @return Vertex bin in @f$[1,N_{\mbox{vertex}}]@f$ or negative if * out of range */ Int_t FindVertexBin(Double_t vertex) const; /** * Get the index corresponding to the given ring * * @param d Detector * @param r Ring * * @return Index (0 based) or negative in case of errors */ Int_t GetRingIndex(UShort_t d, Char_t r) const; /** * Get the ring array corresponding to the specified ring * * @param d Detector * @param r Ring * * @return Pointer to ring array, or null in case of problems */ TObjArray* GetRingArray(UShort_t d, Char_t r) const; /** * Get the ring array corresponding to the specified ring * * @param d Detector * @param r Ring * * @return Pointer to ring array, or newly created container */ TObjArray* GetOrMakeRingArray(UShort_t d, Char_t r); TObjArray fRingArray; // Array of per-ring, per-vertex 2nd map TAxis fVertexAxis; // The vertex axis TAxis fEtaAxis; // The eta axis ClassDef(AliFMDCorrSecondaryMap,1); // }; //____________________________________________________________________ inline void AliFMDCorrSecondaryMap::SetVertexAxis(Int_t nBins, Double_t min, Double_t max) { fVertexAxis.Set(nBins, min, max); } //____________________________________________________________________ inline void AliFMDCorrSecondaryMap::SetVertexAxis(const TAxis& e) { fVertexAxis.Set(e.GetNbins(), e.GetXmin(), e.GetXmax()); } //____________________________________________________________________ inline void AliFMDCorrSecondaryMap::SetEtaAxis(Int_t nBins, Double_t min, Double_t max) { fEtaAxis.Set(nBins, min, max); } //____________________________________________________________________ inline void AliFMDCorrSecondaryMap::SetEtaAxis(const TAxis& e) { fEtaAxis.Set(e.GetNbins(), e.GetXmin(), e.GetXmax()); } #endif // Local Variables: // mode: C++ // End: