Double_t CrystalBallLowHighDivideCrystalBallLowHigh(Double_t *x,Double_t *par) {
Double_t num = 0;
num = crystalBallLowHigh(x,par);
Double_t den = 1;
den = crystalBallLowHigh(x,&par[7]);
if (den != 0) return num/den;
else return 1.;
}
double crystalBallLowHighWithRise (double* x, double* par)
{
//PG the first 7 parameters are passed to the crystalBallLowHigh
// return TMath::Log (par[7] * (x[0] - 180.)) * crystalBallLowHigh (x, par) ;
// return leftRise (x, par + 7) * crystalBallLowHigh (x, par) ;
return leftRise (x, par + 7) + crystalBallLowHigh (x, par) ;
}
//---- division of CBLowHighPlusExp with CBLowHigh ----
Double_t CrystalBallLowHighPlusExpDividedByCrystalBallLowHigh(Double_t *x,Double_t *par) {
double den = crystalBallLowHigh (x, par + 9) ; // signal only
if (den == 0) return -1. ;
double num = doubleGausCrystalBallLowHighPlusExp (x, par) ; // signal and interference
return num / den ;
}
Double_t intereferenceOverSignal (Double_t * xx, Double_t * par)
{
if (crystalBallLowHigh (xx, par + 4) < 0.0000001) return 0 ;
return doublePeakModel (xx, par) / crystalBallLowHigh (xx, par + 4) ;
}
