StCentrality.h

#ifndef __StCentrality_h__
#define __StCentrality_h__

#include <vector>
#include "TString.h"

//____________________________________________________________________________________________________
// Class StCentrality: Centrality class
//   - Centrality: multiplicity cut, and corresponding centrality range
//   - NBD: npp, k and x parameters
class StCentrality {
        public:
                StCentrality(const TString system = "AuAu_200GeV", const TString type="default"); 
                virtual ~StCentrality(); 

                static void help() ; 

                // Getters

                // mode
                //  0    default
                //  1    -5% total cross section
                //  2    +5% total cross section
                Double_t GetCentrality(const UInt_t multiplicity, const UInt_t mode=0) const ; // Get centrality

                Double_t GetNpp()         const ; 
                Double_t GetK()           const ; 
                Double_t GetX()           const ; 
                Double_t GetEfficiency()  const ; 
                Double_t GetTriggerBias() const ; 

                Double_t GetReweighting(const UInt_t multiplicity) const ; 

        private:
                // Functions
                Double_t GetNppError(const Double_t npp) const ; 

                // Initialization
                void Init_AuAu200GeV() ; 
                void Init_AuAu62GeV()  ; 
                void Init_AuAu39GeV()  ; 
                void Init_AuAu27GeV()  ; 
                void Init_AuAu19GeV()  ; 
                void Init_AuAu11GeV()  ; 
                void Init_AuAu7GeV()   ; 
                //    void Init_SmSm200GeV() ; /// Initialization of Sm + Sm 200 GeV
                void Init_dAu200GeV()  ; 

                // Need actual implementation
                void Init_CuCu200GeV() ;
                void Init_ZrZr200GeV() ;
                void Init_UU200GeV() ;
                void Init_PbPb2760GeV() ;

                // Data members
                enum {
                        mNMode = 3
                };

                const TString mType   ; 
                Double_t mNpp         ; 
                Double_t mK           ; 
                Double_t mX           ; 
                Double_t mEfficiency  ; 
                Double_t mTriggerBias ; 
                Double_t mParReweighting[2]; 
                std::vector<UInt_t> mMultiplicityCut[mNMode] ; 
                std::vector<Double_t> mCentralityMin[mNMode] ; 
                std::vector<Double_t> mCentralityMax[mNMode] ; 

                ClassDef(StCentrality, 0)
};

inline Double_t StCentrality::GetNpp()          const { return mNpp ; }
inline Double_t StCentrality::GetK()            const { return mK ; }
inline Double_t StCentrality::GetX()            const { return mX ; }
inline Double_t StCentrality::GetEfficiency()   const { return mEfficiency ; }
inline Double_t StCentrality::GetTriggerBias()  const { return mTriggerBias ; }

#endif