StNbdFitMaker.h

/******************************************************************************
 * $Id: StNbdFitMaker.h,v 1.2 2012/04/25 05:15:30 hmasui Exp $
 * $Log: StNbdFitMaker.h,v $
 * Revision 1.2  2012/04/25 05:15:30  hmasui
 * Added centrality calculation. Real data histogram name can be now specified in ReadData() function
 *
 ******************************************************************************/

#ifndef __StNbdFitMaker_h__
#define __StNbdFitMaker_h__

class TCanvas ;
class TGraph ;
class TH1 ;
class TH2 ;
class TLine ;
class StNegativeBinomial ;

//____________________________________________________________________________________________________
// Class StNbdFitMaker: Determine the NBD parameters from the data
class StNbdFitMaker {
        public:
                StNbdFitMaker();
                virtual ~StNbdFitMaker(); 

                void DoCentralityDetermination() ;

                void ReadData(const Char_t* data, const Char_t* glauber,
                                const Char_t* dataHistogramName = "hRefMultTpc") ;

                //TGraph* Fit(const Int_t nevents = 1000, const Char_t* outputFileName = "");//zaochen comment
                TGraph* Fit(const Int_t nevents = 1000, const TString outputFileName = "");

                Int_t Scan(const Int_t nevents,
                                const Int_t nppbin, const Double_t nppmin, const Double_t nppmax,
                                const Int_t kbin, const Double_t kmin, const Double_t kmax,
                                const Int_t xbin, const Double_t xmin, const Double_t xmax,
                                //        const Double_t x,
                                const Double_t efficiency=1.0,
                                //        const Int_t effbin, const Double_t effmin, const Double_t effmax,
                                const Double_t triggerbias=1.0, const Bool_t isConstEfficiency=kTRUE
                                );

                void SetParameters(const Double_t npp, const Double_t k, const Double_t x,
                                const Double_t efficiency, const Double_t triggerbias, const Bool_t isConstEfficiency) ;

                void SetMinimumMultiplicityCut(const Double_t cut) ;

        private:
                // Functions
                Double_t GetNormalization(const TH1& h1, const TH1& h2,
                                const Double_t min, const Double_t max) const ; // Get normalization factor in (min, max)
                Double_t CalculateChi2(const TH1& hdata, const TH1& hfunc,
                                const Double_t minimumMultiplicityCut) ; // Get chi2 from data and func

                void CalculateCentrality(const TH1& hdata, const TH1& hmc) const ;

                // Data members
                StNegativeBinomial* mNBinomial ; 
                TH1* mhRefMult     ; 
                TH1* mhRefMultSim  ; 
                TH2* mhNcoll_Npart ; 
                Int_t mNData ; 
                Double_t mMinimumMultiplicityCut ; 
                Bool_t mDoCentralityDetermination ; 

                TCanvas* mCanvas   ; 
                TLine* mCutOff[2]  ; 
                TLine* mOneLine    ; 
                TGraph* mChi2Graph ; 

                ClassDef(StNbdFitMaker, 0)
};
#endif