class AliDielectronBtoJPSItoEleCDFfitFCN: public TNamed
Class AliDielectronBtoJPSItoEleCDFfitFCN
Definition of main function used in
unbinned log-likelihood fit for
the channel B -> JPsi + X -> e+e- + X
Origin: C.Di Giglio
Contact: Carmelo.Digiglio@ba.infn.it , Giuseppe.Bruno@ba.infn.it
Function Members (Methods)
public:
protected:
| virtual void | TObject::DoError(int level, const char* location, const char* fmt, va_list va) const |
| void | TObject::MakeZombie() |
private:
Data Members
public:
| enum TObject::EStatusBits { | kCanDelete | |
| kMustCleanup | ||
| kObjInCanvas | ||
| kIsReferenced | ||
| kHasUUID | ||
| kCannotPick | ||
| kNoContextMenu | ||
| kInvalidObject | ||
| }; | ||
| enum TObject::[unnamed] { | kIsOnHeap | |
| kNotDeleted | ||
| kZombie | ||
| kBitMask | ||
| kSingleKey | ||
| kOverwrite | ||
| kWriteDelete | ||
| }; |
protected:
| TString | TNamed::fName | object identifier |
| TString | TNamed::fTitle | object title |
private:
| Float_t**** | fBkgParams | x background parameters in pt, mass, type bins |
| Double_t | fChangeMass | constant to change the RMS of the signal mass function |
| Double_t | fChangeResolution | constant to change the RMS of the resolution function |
| Bool_t | fCrystalBallParam | Boolean to switch to Crystall Ball parameterisation |
| Bool_t | fExponentialParam | Boolean to switch to Exponential parameterisation |
| Double_t | fFMinus | slopes of the x distributions of the background |
| Double_t | fFPlus | parameter of the log-likelihood function |
| Double_t | fFSym | functions |
| TF1*** | fFunBSaved | pointers to save functions for x of non-prompt J/psi in pt, mass, type bins |
| TF1**** | fFunBkgSaved | pointers to save functions for x of background in pt, mass, type bins |
| Bool_t | fLoadFunctions | boolean to load functions saved |
| TArrayD* | fMassWindows | limits for invariant mass bins |
| Double_t | fMassWndHigh | JPSI Mass window higher limit |
| Double_t | fMassWndLow | JPSI Mass window lower limit |
| Bool_t | fMultivariate | switch-on multivariate fit |
| Double_t | fParameters[49] | par[0] = weightRes; |
| TArrayD* | fPtWindows | limits for pt bins |
| Double_t*** | fResParams | resolution function parameters in pt and type bins |
| Double_t | fShiftTemplate | to shift the MC template |
| Double_t | fSignalBinForExtrapolation | inv. mass region in which extrapolate the background shape |
| Double_t | fWeightType[3] | vector with weights of candidates types (used to draw functions) |
| Double_t*** | fWeights | weights to interpolate bkg shape in pt, mass, type bins under signal region |
| TH1F* | fhCsiMC | X distribution used as MC template for JPSI from B |
| Double_t | fintmMassBkg | integral of invariant mass distribution for the bkg |
| Double_t | fintmMassSig | integral of invariant mass distribution for the signal |
Class Charts
| Inheritance Chart: | ||||||||
|
Function documentation
AliDielectronBtoJPSItoEleCDFfitFCN(const AliDielectronBtoJPSItoEleCDFfitFCN& source)
Copy constructor
AliDielectronBtoJPSItoEleCDFfitFCN& operator=(const AliDielectronBtoJPSItoEleCDFfitFCN& source)
Assignment operator
Double_t EvaluateLikelihood(const Double_t* pseudoproperdecaytime, const Double_t* invariantmass, const Double_t* pt, const Int_t* type, Int_t ncand) const
This function evaluates the Likelihood fnction It returns the -Log(of the likelihood function)
void ComputeMassIntegral()
this function compute the integral of the likelihood function (theoretical function) in order to normalize it to unity
void SetResolutionConstants(const Double_t* resolutionConst, Int_t type)
Resolution function is parametrized as the sum of two gaussian
Double_t EvaluateCDFfunc(Double_t x, Double_t m, Double_t pt, Int_t type) const
evaluate likelihood function
printf("CDF func ---> x = %f m = %f pt = %f type = %d \n",x,m,pt,type);
Double_t EvaluateCDFfuncNorm(Double_t x, Double_t m, Double_t pt, Int_t type) const
evaluate likelihood function
Double_t EvaluateCDFfuncSignalPart(Double_t x, Double_t m, Double_t pt, Int_t type) const
evaluate psproper signal
Double_t EvaluateCDFDecayTimeSigDistr(Double_t x, Double_t pt, Int_t type) const
Implementation of the Background part of the Likelyhood function
Double_t EvaluateCDFInvMassSigDistr(Double_t m) const
Parametrization of signal part invariant mass distribution It can be either Crystal Ball function or sum of two Landau
Double_t FunB(Double_t x, Double_t pt, Int_t type) const
Parameterisation of the fit function for the x part of the Background
Double_t FunP(Double_t x, Double_t pt, Int_t type) const
Parameterisation of the Prompt part for the x distribution
Double_t CsiMC(Double_t x) const
Distribution (template) of the x distribution for the x variable for the J/psi coming from Beauty hadrons
Double_t EvaluateCDFfuncBkgPart(Double_t x, Double_t m, Double_t pt, Int_t type) const
Return the part of the likelihood function for the background hypothesis
Double_t EvaluateCDFDecayTimeBkgDistr(Double_t x, Int_t type, Double_t m = 3.09, Double_t pt = 200.) const
it returns the value of the probability to have a given x for the background in the pt, m , type correspondent range
Double_t EvaluateCDFInvMassBkgDistr(Double_t m) const
it returns the value of the probability to have a given mass for the background
Double_t FunBkgPos(Double_t x, Double_t pt, Int_t type) const
exponential with positive slopes for the background part (x)
Double_t FunBkgNeg(Double_t x, Double_t pt, Int_t type) const
exponential with negative slopes for the background part (x)
Double_t FunBkgSym(Double_t x, Double_t pt, Int_t type) const
exponential with both positive and negative slopes for the background part (x)
Double_t FunBkgSym1(Double_t x, Double_t pt, Int_t type) const
exponential with both positive and negative slopes for the background part (x)
TF1* GetCsiMC(Double_t xmin, Double_t xmax, Double_t normalization)
return the pointer to the templateMC function
TF1* GetResolutionFunc(Double_t xmin, Double_t xmax, Double_t normalization, Double_t pt, Int_t type = 2)
return the pointer to the resolution function
TF1* GetResolutionFuncAllTypes(Double_t xmin, Double_t xmax, Double_t normalization)
return the pointer to the resolution function
TF1* GetEvaluateCDFDecayTimeBkgDistr(Double_t xmin, Double_t xmax, Double_t normalization, Int_t type = 2, Double_t mass = 3.09, Double_t pt = 200., Int_t npx = 5000)
return the pointer to the background x distribution function
TF1* GetEvaluateCDFDecayTimeBkgDistrAllTypes(Double_t xmin, Double_t xmax, Double_t normalization)
return the pointer to the background x distribution function
TF1* GetEvaluateCDFDecayTimeSigDistr(Double_t xmin, Double_t xmax, Double_t normalization, Double_t type)
return the pointer to the signal x distribution function
TF1* GetEvaluateCDFInvMassBkgDistr(Double_t mMin, Double_t mMax, Double_t normalization)
return the pointer to the invariant mass distribution for the background
TF1* GetEvaluateCDFInvMassSigDistr(Double_t mMin, Double_t mMax, Double_t normalization)
return the pointer to the invariant mass distribution for the signal
TF1 * GetEvaluateCDFInvMassTotalDistr(Double_t mMin, Double_t mMax, Double_t normalization)
return the pointer to the invariant mass distribution
Double_t EvaluateCDFInvMassTotalDistr(const Double_t* x, const Double_t* par) const
evaluate invariant mass total distribution
TF1 * GetEvaluateCDFDecayTimeTotalDistr(Double_t xMin, Double_t xMax, Double_t normalization, Double_t pt = 200., Int_t type = 2)
return the pointer to the pseudoproper distribution for the background
Double_t EvaluateCDFDecayTimeTotalDistr(const Double_t* x, const Double_t* par) const
evaluate the total pseudoproper distribution for a given candidate's type. par[1] should be the candidate's type.
Double_t EvaluateCDFDecayTimeTotalDistrAllTypes(const Double_t* x, const Double_t* par) const
evaluate the total pseudoproper distribution considering all candidate's types
TF1 * GetEvaluateCDFDecayTimeTotalDistrAllTypes(Double_t xMin, Double_t xMax, Double_t normalization)
return the pointer to the pseudoproper function for the background considering all candidate's types
TF1 * GetFunB(Double_t xmin, Double_t xmax, Double_t normalization, Double_t pt, Int_t type = 2, Int_t npx = 5000)
return the pointer to the function that describe secondary jpsi pseudoproper distribution for a given candidate's type
TF1 * GetFunBAllTypes(Double_t xmin, Double_t xmax, Double_t normalization)
return the pointer to the function that describe secondary jpsi pseudoproper distribution for all types
void SetBackgroundSpecificParameters(Int_t pt, Int_t mb, Int_t tp)
methods to set specific background parameters in bins(pt, inv. mass, type)
void InitializeFunctions(Int_t ptSize, Int_t massSize)
initialize pointers to save functions for the multivariate fit
Double_t FunBsaved(Double_t x, Double_t pt, Int_t type) const
functions to describe non-prompt J/psi x distributions
void SetFunctionsSaved(Int_t npxFunB = 5000, Int_t npxFunBkg = 5000, Double_t funBLimits = 20000., Double_t funBkgLimits = 40000., Int_t signalRegion = 2)
save functions for the multivariate fit
Double_t EvaluateCDFDecayTimeBkgDistrDifferential(Double_t x, Int_t type, Double_t m = 3.09, Double_t pt = 200.) const
it returns the value of the probability to have a given x for the background in the pt, m , type correspondent range
Double_t EvaluateCDFDecayTimeBkgDistrSaved(Double_t x, Int_t type, Double_t m = 3.09, Double_t pt = 200.) const
it returns the value of the probability to have a given x for the background in the pt, m , type correspondent range
void SetMassWndHigh(Double_t limit)
{ fMassWndHigh = TDatabasePDG::Instance()->GetParticle(443)->Mass() + limit ;}void SetMassWndLow(Double_t limit)
{ fMassWndLow = TDatabasePDG::Instance()->GetParticle(443)->Mass() - limit ;}void SetWeightType(Double_t wFF, Double_t wFS, Double_t wSS)
{fWeightType[0]= wSS; fWeightType[1]= wFS; fWeightType[2]= wFF;}void ReadMCtemplates(Int_t BinNum)
void SetBkgFunction(Int_t massRange, Int_t type, Int_t ptB, TF1* histBkg)
TF1 * GetBkgFunction(Int_t massRange, Int_t ptB, Int_t type) const
{return fFunBkgSaved[ptB][massRange][type];}void SetExtrapolationRegion(Int_t extrRegion)
Double_t EvaluateCDFDecayTimeSigDistrFunc(const Double_t* x, const Double_t* par) const
{ return par[0]*EvaluateCDFDecayTimeSigDistr(x[0],par[1],(Int_t)par[2]);}Double_t EvaluateCDFInvMassSigDistrFunc(const Double_t* x, const Double_t* par) const
{return par[0]*EvaluateCDFInvMassSigDistr(x[0])/fintmMassSig;}Double_t EvaluateCDFDecayTimeBkgDistrFunc(const Double_t* x, const Double_t* par) const
{ return EvaluateCDFDecayTimeBkgDistr(x[0],(Int_t)par[1],par[2],par[3])*par[0];}Double_t EvaluateCDFDecayTimeBkgDistrFuncAllTypes(const Double_t* x, const Double_t* par) const
{return (fWeightType[2]*EvaluateCDFDecayTimeBkgDistr(x[0],2)+fWeightType[1]*EvaluateCDFDecayTimeBkgDistr(x[0],1)+fWeightType[0]*EvaluateCDFDecayTimeBkgDistr(x[0],0))*par[0];}Double_t EvaluateCDFInvMassBkgDistrFunc(const Double_t* x, const Double_t* par) const
{return par[0]*EvaluateCDFInvMassBkgDistr(x[0])/fintmMassBkg;}Double_t FunBfuncAllTypes(const Double_t* x, const Double_t* par) const
{return (fWeightType[2]*FunB(x[0],200.,2)+fWeightType[1]*FunB(x[0],200.,1)+fWeightType[0]*FunB(x[0],200.,0))*par[0];}Double_t ResolutionFuncf(const Double_t* x, const Double_t* par) const
{ return ResolutionFunc(x[0],par[1],(Int_t)par[2])*par[0];}Double_t ResolutionFuncAllTypes(const Double_t* x, const Double_t* par) const
{ return (fWeightType[2]*ResolutionFunc(x[0],200.,2)+fWeightType[1]*ResolutionFunc(x[0],200.,1)+fWeightType[0]*ResolutionFunc(x[0],200.,0))*par[0]; }