ROOT logo
AliRoot » PWGPP » EVCHAR » AliCentralityGlauberFit

class AliCentralityGlauberFit: public TObject


class to determine centrality percentiles from 1D distributions


Function Members (Methods)

public:
AliCentralityGlauberFit(const char* filename)
virtual~AliCentralityGlauberFit()
voidTObject::AbstractMethod(const char* method) const
voidAddHisto(TString name)
virtual voidTObject::AppendPad(Option_t* option = "")
virtual voidTObject::Browse(TBrowser* b)
static TClass*Class()
virtual const char*TObject::ClassName() const
virtual voidTObject::Clear(Option_t* = "")
virtual TObject*TObject::Clone(const char* newname = "") const
virtual Int_tTObject::Compare(const TObject* obj) const
virtual voidTObject::Copy(TObject& object) const
virtual voidTObject::Delete(Option_t* option = "")MENU
virtual Int_tTObject::DistancetoPrimitive(Int_t px, Int_t py)
virtual voidTObject::Draw(Option_t* option = "")
virtual voidTObject::DrawClass() constMENU
virtual TObject*TObject::DrawClone(Option_t* option = "") constMENU
virtual voidTObject::Dump() constMENU
virtual voidTObject::Error(const char* method, const char* msgfmt) const
virtual voidTObject::Execute(const char* method, const char* params, Int_t* error = 0)
virtual voidTObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0)
virtual voidTObject::ExecuteEvent(Int_t event, Int_t px, Int_t py)
virtual voidTObject::Fatal(const char* method, const char* msgfmt) const
virtual TObject*TObject::FindObject(const char* name) const
virtual TObject*TObject::FindObject(const TObject* obj) const
virtual Option_t*TObject::GetDrawOption() const
static Long_tTObject::GetDtorOnly()
TH1F*GetGlauberHist()
virtual const char*TObject::GetIconName() const
virtual const char*TObject::GetName() const
virtual char*TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_tTObject::GetObjectStat()
virtual Option_t*TObject::GetOption() const
TH1F*GetTempHist()
virtual const char*TObject::GetTitle() const
virtual UInt_tTObject::GetUniqueID() const
virtual Bool_tTObject::HandleTimer(TTimer* timer)
virtual ULong_tTObject::Hash() const
virtual voidTObject::Info(const char* method, const char* msgfmt) const
virtual Bool_tTObject::InheritsFrom(const char* classname) const
virtual Bool_tTObject::InheritsFrom(const TClass* cl) const
virtual voidTObject::Inspect() constMENU
voidTObject::InvertBit(UInt_t f)
virtual TClass*IsA() const
virtual Bool_tTObject::IsEqual(const TObject* obj) const
virtual Bool_tTObject::IsFolder() const
Bool_tTObject::IsOnHeap() const
virtual Bool_tTObject::IsSortable() const
Bool_tTObject::IsZombie() const
virtual voidTObject::ls(Option_t* option = "") const
voidMakeFits()
voidMakeFitsMinuitNBD(Double_t alpha = -1, Double_t mu = -1, Double_t k = -1)
voidTObject::MayNotUse(const char* method) const
virtual Bool_tTObject::Notify()
voidTObject::Obsolete(const char* method, const char* asOfVers, const char* removedFromVers) const
static voidTObject::operator delete(void* ptr)
static voidTObject::operator delete(void* ptr, void* vp)
static voidTObject::operator delete[](void* ptr)
static voidTObject::operator delete[](void* ptr, void* vp)
void*TObject::operator new(size_t sz)
void*TObject::operator new(size_t sz, void* vp)
void*TObject::operator new[](size_t sz)
void*TObject::operator new[](size_t sz, void* vp)
virtual voidTObject::Paint(Option_t* option = "")
virtual voidTObject::Pop()
virtual voidTObject::Print(Option_t* option = "") const
virtual Int_tTObject::Read(const char* name)
virtual voidTObject::RecursiveRemove(TObject* obj)
voidTObject::ResetBit(UInt_t f)
virtual voidTObject::SaveAs(const char* filename = "", Option_t* option = "") constMENU
virtual voidTObject::SavePrimitive(basic_ostream<char,char_traits<char> >& out, Option_t* option = "")
voidSetAncestorMode(Int_t f)
voidSetAncFile(TString name)
voidTObject::SetBit(UInt_t f)
voidTObject::SetBit(UInt_t f, Bool_t set)
virtual voidTObject::SetDrawOption(Option_t* option = "")MENU
static voidTObject::SetDtorOnly(void* obj)
voidSetFastFitMode(Int_t f)
voidSetGlauberParam(Int_t Nmu, Double_t mulow, Double_t muhigh, Int_t Nk, Double_t klow, Double_t khigh, Int_t Nalpha, Double_t alphalow, Double_t alphahigh)
voidSetInputFile(TString filename)
voidSetInputNtuple(TString ntuplename)
voidSetNevents(Int_t f)
voidSetNtrials(Int_t f)
static voidTObject::SetObjectStat(Bool_t stat)
voidSetOutputFile(TString filename)
voidSetOutputNtuple(TString ntuplename)
voidSetRangeToFit(Double_t multmin, Double_t multmax)
voidSetRangeToScale(Double_t scalemin)
voidSetRebin(Int_t f)
virtual voidTObject::SetUniqueID(UInt_t uid)
virtual voidShowMembers(TMemberInspector&)
virtual voidStreamer(TBuffer&)
voidStreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b)
virtual voidTObject::SysError(const char* method, const char* msgfmt) const
Bool_tTObject::TestBit(UInt_t f) const
Int_tTObject::TestBits(UInt_t f) const
voidUseAverage(Bool_t f)
voidUseChi2(Bool_t f)
virtual voidTObject::UseCurrentStyle()
virtual voidTObject::Warning(const char* method, const char* msgfmt) const
virtual Int_tTObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0)
virtual Int_tTObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const
protected:
virtual voidTObject::DoError(int level, const char* location, const char* fmt, va_list va) const
voidTObject::MakeZombie()
private:
AliCentralityGlauberFit(const AliCentralityGlauberFit&)
Double_tCalculateChi2(TH1F* hDATA, TH1F* thistGlau)
TH1F*GetTriggerEfficiencyFunction(TH1F* hist1, TH1F* hist2)
Double_tGetTriggerEfficiencyIntegral(TH1F* hist1, TH1F* hist2)
TH1F*GlauberHisto(Double_t mu, Double_t k, Double_t alpha, TH1F* hDATA, Bool_t save = kFALSE)
TH1F*MakeAncestor(Double_t alpha)
static voidMinuitFcnNBD(Int_t& npar, Double_t* gin, Double_t& f, Double_t* par, Int_t iflag)
Double_tNBD(Int_t n, Double_t mu, Double_t k) const
static Double_tNBDFunc(const Double_t* p, const Double_t* x)
TH1F*NBDhist(Double_t mu, Double_t k)
TH1F*NormalizeHisto(TString hdistributionName)
AliCentralityGlauberFit&operator=(const AliCentralityGlauberFit&)
voidSaveHisto(TH1F* hist1, TH1F* hist2, TH1F* heffi, TFile* outrootfile)

Data Members

private:
Double_tfAlphahighalpha high
Double_tfAlphalowalpha low
Int_tfAncestorIf 1 use Npart**alpha, if 2 use alpha*Npart + (1-alpha)*Ncoll
TStringfAncfilenameancestor file name
Float_tfBimpact parameter
Float_tfEffiefficiency
Int_tfFastFitIf 1 or 2 use cruder approximation to compute curve faster 1:NBD, 2:Gauss
TH1F*fGlauberHistGlauber histogram
TNtuple*fGlauntupleglauber ntuple
vector<TString>fHistnameshistogram names
TStringfInntuplenameinput Gauber ntuple
TStringfInrootfilenameinput root file
Double_tfKhighk high
Double_tfKlowk low
Double_tfMuhighmu high
Double_tfMulowmu low
Double_tfMultmaxmult max
Double_tfMultminmult min
TF1*fNBDNBD function
Int_tfNalphaalpha points
Float_tfNcollnumber collisions
Int_tfNeventsNumber of events to use in the glauber ntuple
Int_tfNkk points
Int_tfNmumu points
Float_tfNpartnumber participants
Int_tfNtrialsNumber of trials to use for the average
TStringfOutntuplenameoutput Glauber ntuple
TStringfOutrootfilenameoutput root file
Int_tfRebinFactorrebin factor
Double_tfScaleminmult min where to scale
Float_tfTaataa
TH1F*fTempHistTemporary pointer to data histo, to be used in minimization
Bool_tfUseAverageIf true, use average
Bool_tfUseChi2If true, use chi2
TH1F*fhAncestorhisto for the ancestor distribution
TH1F*fhEffiefficiency histogram

Class Charts

Inheritance Chart:
TObject
AliCentralityGlauberFit

Function documentation

AliCentralityGlauberFit(const char* filename)
 Standard constructor.
void SetRangeToFit(Double_t multmin, Double_t multmax)
 Set fit range.
void SetRangeToScale(Double_t scalemin)
 Set range where to scale simulated histo to data
void SetGlauberParam(Int_t Nmu, Double_t mulow, Double_t muhigh, Int_t Nk, Double_t klow, Double_t khigh, Int_t Nalpha, Double_t alphalow, Double_t alphahigh)
 Set Glauber parameters.
void MakeFits()
 Make fits.
void MakeFitsMinuitNBD(Double_t alpha = -1, Double_t mu = -1, Double_t k = -1)
 Make fits using Minuit.
TH1F * GlauberHisto(Double_t mu, Double_t k, Double_t alpha, TH1F* hDATA, Bool_t save = kFALSE)
 Get Glauber histogram.
Double_t CalculateChi2(TH1F* hDATA, TH1F* thistGlau)
 Note that for different values of neff the mc distribution is identical, just re-scaled below
 normalize the two histogram, scale MC up but leave data alone - that preserves error bars !!!!
TH1F * GetTriggerEfficiencyFunction(TH1F* hist1, TH1F* hist2)
 Get efficiency.
Double_t GetTriggerEfficiencyIntegral(TH1F* hist1, TH1F* hist2)
 Get eff integral.
void SaveHisto(TH1F* hist1, TH1F* hist2, TH1F* heffi, TFile* outrootfile)
 Save histograms.
Double_t NBD(Int_t n, Double_t mu, Double_t k) const
 Compute NBD.
TH1F * NBDhist(Double_t mu, Double_t k)
 Interface for TH1F.
void MinuitFcnNBD(Int_t& npar, Double_t* gin, Double_t& f, Double_t* par, Int_t iflag)
 FCN for minuit.
Double_t NBDFunc(const Double_t* p, const Double_t* x)
 TF1  interface.
TH1F * MakeAncestor(Double_t alpha)
 Make ancestor histogram.
AliCentralityGlauberFit(const char* filename)
virtual ~AliCentralityGlauberFit()
{}
void AddHisto(TString name)
{ fHistnames.push_back(name); }
TH1F * GetGlauberHist()
{ return fGlauberHist; }
TH1F * GetTempHist()
{ return fTempHist; }
void SetAncestorMode(Int_t f)
{ fAncestor=f; }
void SetFastFitMode(Int_t f)
{ fFastFit=f; }
void SetInputFile(TString filename)
{ fInrootfilename = filename; }
void SetInputNtuple(TString ntuplename)
{ fInntuplename = ntuplename; }
void SetAncFile(TString name)
{ fAncfilename = name; }
void SetNevents(Int_t f)
{ fNevents=f; }
void SetNtrials(Int_t f)
{ fNtrials=f; }
void SetOutputFile(TString filename)
{ fOutrootfilename = filename; }
void SetOutputNtuple(TString ntuplename)
{ fOutntuplename = ntuplename;}
void SetRebin(Int_t f)
void UseAverage(Bool_t f)
void UseChi2(Bool_t f)
{ fUseChi2=f; }
TH1F * NormalizeHisto(TString hdistributionName)
AliCentralityGlauberFit & operator=(const AliCentralityGlauberFit& )