MELADoubleCB Class Reference

#include <MELACombinePdfs.h>

Inheritance diagram for MELADoubleCB:

Collaboration diagram for MELADoubleCB:

Public Member Functions
	MELADoubleCB ()
	MELADoubleCB (const char *name, const char *title, RooAbsReal &_x, RooAbsReal &_mean, RooAbsReal &_width, RooAbsReal &_alpha1, RooAbsReal &_n1, RooAbsReal &_alpha2, RooAbsReal &_n2)
	MELADoubleCB (const MELADoubleCB &other, const char *name=0)
virtual TObject *	clone (const char *newname) const
virtual	~MELADoubleCB ()
Int_t	getAnalyticalIntegral (RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=0) const
Double_t	analyticalIntegral (Int_t code, const char *rangeName=0) const

Protected Member Functions
Double_t	evaluate () const

Protected Attributes
RooRealProxy	x
RooRealProxy	mean
RooRealProxy	width
RooRealProxy	alpha1
RooRealProxy	n1
RooRealProxy	alpha2
RooRealProxy	n2

Detailed Description

Definition at line 102 of file MELACombinePdfs.h.

Constructor & Destructor Documentation

MELADoubleCB() [1/3]

MELADoubleCB::MELADoubleCB

(

)

Definition at line 1998 of file MELACombinePdfs.cc.

{}

MELADoubleCB() [2/3]

MELADoubleCB::MELADoubleCB	(	const char *	name,
		const char *	title,
		RooAbsReal &	_x,
		RooAbsReal &	_mean,
		RooAbsReal &	_width,
		RooAbsReal &	_alpha1,
		RooAbsReal &	_n1,
		RooAbsReal &	_alpha2,
		RooAbsReal &	_n2
	)

Definition at line 2000 of file MELACombinePdfs.cc.

    :
  RooAbsPdf(name, title),
  x("x", "x", this, _x),
  mean("mean", "mean", this, _mean),
  width("width", "width", this, _width),
  alpha1("alpha1", "alpha1", this, _alpha1),
  n1("n1", "n1", this, _n1),
  alpha2("alpha2", "alpha2", this, _alpha2),
  n2("n2", "n2", this, _n2)
{}

MELADoubleCB() [3/3]

MELADoubleCB::MELADoubleCB	(	const MELADoubleCB &	other,
		const char *	name = 0
	)

Definition at line 2020 of file MELACombinePdfs.cc.

                                                                      :
RooAbsPdf(other, name),
x("x", this, other.x),
mean("mean", this, other.mean),
width("width", this, other.width),
alpha1("alpha1", this, other.alpha1),
n1("n1", this, other.n1),
alpha2("alpha2", this, other.alpha2),
n2("n2", this, other.n2)
{}

~MELADoubleCB()

virtual MELADoubleCB::~MELADoubleCB ( )

inlinevirtual

Definition at line 117 of file MELACombinePdfs.h.

{}

Member Function Documentation

analyticalIntegral()

Double_t MELADoubleCB::analyticalIntegral	(	Int_t	code,
		const char *	rangeName = 0
	)		const

Definition at line 2057 of file MELACombinePdfs.cc.

                                                                                {
  assert(code==1);
 
  double central=0;
  double left=0;
  double right=0;
 
  static const Double_t root2 = sqrt(2);
  static const Double_t rootPiBy2 = sqrt(atan2(0.0, -1.0)/2.0);
  Double_t xscale = root2*width;
 
  //compute gaussian contribution
  double central_low =max(x.min(rangeName), mean - alpha1*width);
  double central_high=min(x.max(rangeName), mean + alpha2*width);
  if (central_low < central_high) // is the gaussian part in range?
    central = rootPiBy2*width*(TMath::Erf((central_high-mean)/xscale)-TMath::Erf((central_low-mean)/xscale));
 
  //compute left tail;
  double A1 = pow(n1/fabs(alpha1), n1)*exp(-alpha1*alpha1/2);
  double B1 = n1/fabs(alpha1)-fabs(alpha1);
 
  double left_low=x.min(rangeName);
  double left_high=min(x.max(rangeName), mean - alpha1*width);
  if (left_low < left_high){ //is the left tail in range?
    if (fabs(n1-1.0)>1.e-5)
      left = A1/(-n1+1.0)*width*(pow(B1-(left_low-mean)/width, -n1+1.)-pow(B1-(left_high-mean)/width, -n1+1.));
    else
      left = A1*width*(log(B1-(left_low-mean)/width) - log(B1-(left_high-mean)/width));
  }
 
  //compute right tail;
  double A2 = pow(n2/fabs(alpha2), n2)*exp(-alpha2*alpha2/2);
  double B2 = n2/fabs(alpha2)-fabs(alpha2);
 
  double right_low=max(x.min(rangeName), mean + alpha2*width);
  double right_high=x.max(rangeName);
  if (right_low < right_high){ //is the right tail in range?
    if (fabs(n2-1.0)>1.e-5)
      right = A2/(-n2+1.0)*width*(pow(B2+(right_high-mean)/width, -n2+1.)-pow(B2+(right_low-mean)/width, -n2+1.));
    else
      right = A2*width*(log(B2+(right_high-mean)/width) - log(B2+(right_low-mean)/width));
  }
 
  return left+central+right;
}

clone()

virtual TObject* MELADoubleCB::clone ( const char *  newname ) const

inlinevirtual

Definition at line 116 of file MELACombinePdfs.h.

{ return new MELADoubleCB(*this, newname); }

evaluate()

double MELADoubleCB::evaluate ( ) const

protected

Definition at line 2031 of file MELACombinePdfs.cc.

                                   {
  double t = (x-mean)/width;
  if (t>-alpha1 && t<alpha2){
    return exp(-0.5*t*t);
  }
  else if (t<-alpha1){
    double A1 = pow(n1/fabs(alpha1), n1)*exp(-alpha1*alpha1/2);
    double B1 = n1/fabs(alpha1)-fabs(alpha1);
    return A1*pow(B1-t, -n1);
  }
  else if (t>alpha2){
    double A2 = pow(n2/fabs(alpha2), n2)*exp(-alpha2*alpha2/2);
    double B2 = n2/fabs(alpha2)-fabs(alpha2);
    return A2*pow(B2+t, -n2);
  }
  else{
    cout << "ERROR evaluating range..." << endl;
    return 99;
  }
}

getAnalyticalIntegral()

Int_t MELADoubleCB::getAnalyticalIntegral	(	RooArgSet &	allVars,
		RooArgSet &	analVars,
		const char *	rangeName = 0
	)		const

Definition at line 2052 of file MELACombinePdfs.cc.

                                                                                                         {
  if (matchArgs(allVars, analVars, x)) return 1;
  return 0;
}

Member Data Documentation

alpha1

RooRealProxy MELADoubleCB::alpha1

protected

Definition at line 126 of file MELACombinePdfs.h.

alpha2

RooRealProxy MELADoubleCB::alpha2

protected

Definition at line 128 of file MELACombinePdfs.h.

mean

RooRealProxy MELADoubleCB::mean

protected

Definition at line 124 of file MELACombinePdfs.h.

n1

RooRealProxy MELADoubleCB::n1

protected

Definition at line 127 of file MELACombinePdfs.h.

n2

RooRealProxy MELADoubleCB::n2

protected

Definition at line 129 of file MELACombinePdfs.h.

width

RooRealProxy MELADoubleCB::width

protected

Definition at line 125 of file MELACombinePdfs.h.

x

RooRealProxy MELADoubleCB::x

protected

Definition at line 123 of file MELACombinePdfs.h.

---

The documentation for this class was generated from the following files:

• MELA/interface/MELACombinePdfs.h
• MELA/src/MELACombinePdfs.cc