![]() |
![]() |
#include <HZZ4LRooPdfs.h>
Public Member Functions | |
virtual TObject * | clone (const char *newname) const |
Roo4lMasses2D () | |
Roo4lMasses2D (const char *name, const char *title, RooAbsReal &_mZstar, RooAbsReal &_mZ, RooAbsReal &_mZZ, RooAbsReal &_Gamma, RooAbsReal &_p0, RooAbsReal &_p1, RooAbsReal &_p2, RooAbsReal &_CBmean, RooAbsReal &_CBwidth, RooAbsReal &_CBalpha, RooAbsReal &_CBn) | |
Roo4lMasses2D (const Roo4lMasses2D &other, const char *name=0) | |
virtual | ~Roo4lMasses2D () |
Protected Member Functions | |
Double_t | evaluate () const |
Protected Attributes | |
RooRealProxy | CBalpha |
RooRealProxy | CBmean |
RooRealProxy | CBn |
RooRealProxy | CBwidth |
RooRealProxy | Gamma |
RooRealProxy | mZ |
RooRealProxy | mZstar |
RooRealProxy | mZZ |
RooRealProxy | p0 |
RooRealProxy | p1 |
RooRealProxy | p2 |
Definition at line 310 of file HZZ4LRooPdfs.h.
Roo4lMasses2D::Roo4lMasses2D | ( | ) |
Referenced by clone().
Roo4lMasses2D::Roo4lMasses2D | ( | const char * | name, |
const char * | title, | ||
RooAbsReal & | _mZstar, | ||
RooAbsReal & | _mZ, | ||
RooAbsReal & | _mZZ, | ||
RooAbsReal & | _Gamma, | ||
RooAbsReal & | _p0, | ||
RooAbsReal & | _p1, | ||
RooAbsReal & | _p2, | ||
RooAbsReal & | _CBmean, | ||
RooAbsReal & | _CBwidth, | ||
RooAbsReal & | _CBalpha, | ||
RooAbsReal & | _CBn | ||
) |
Definition at line 2705 of file HZZ4LRooPdfs.cc.
: RooAbsPdf(name,title), mZstar("mZstar","mZstar",this,_mZstar), mZ("mZ","mZ",this,_mZ), mZZ("mZZ","mZZ",this,_mZZ), Gamma("Gamma","Gamma",this,_Gamma), p0("p0","p0",this,_p0), p1("p1","p1",this,_p1), p2("p2","p2",this,_p2), CBmean("CBmean","CBmean",this,_CBmean), CBwidth("CBwidth","CBwidth",this,_CBwidth), CBalpha("CBalpha","CBalpha",this,_CBalpha), CBn("CBn","CBn",this,_CBn) { }
Roo4lMasses2D::Roo4lMasses2D | ( | const Roo4lMasses2D & | other, |
const char * | name = 0 |
||
) |
Definition at line 2735 of file HZZ4LRooPdfs.cc.
: RooAbsPdf(other,name), mZstar("mZstar",this,other.mZstar), mZ("mZ",this,other.mZ), mZZ("mZZ",this,other.mZZ), Gamma("Gamma",this,other.Gamma), p0("p0",this,other.p0), p1("p1",this,other.p1), p2("p2",this,other.p2), CBmean("CBmean" ,this,other.CBmean), CBwidth("CBwidth",this,other.CBwidth), CBalpha("CBalpha",this,other.CBalpha), CBn("CBn",this,other.CBn) { }
virtual Roo4lMasses2D::~Roo4lMasses2D | ( | ) | [inline, virtual] |
Definition at line 328 of file HZZ4LRooPdfs.h.
{ }
virtual TObject* Roo4lMasses2D::clone | ( | const char * | newname | ) | const [inline, virtual] |
Definition at line 327 of file HZZ4LRooPdfs.h.
References Roo4lMasses2D().
{ return new Roo4lMasses2D(*this,newname); }
double Roo4lMasses2D::evaluate | ( | ) | const [protected] |
*
* / CB Shape
Definition at line 2752 of file HZZ4LRooPdfs.cc.
References a, b, beta, CBalpha, CBmean, CBn, CBwidth, alignCSCRings::e, create_public_lumi_plots::exp, Gamma, mZ, mZstar, mZZ, p0, p1, p2, mathSSE::sqrt(), and lumiQTWidget::t.
{ double mZstarDistribution; double Numerator,Denominator,beta; beta = (1-(mZstar-mZ)*(mZstar-mZ)/(mZZ*mZZ))*(1-(mZstar+mZ)*(mZstar+mZ)/(mZZ*mZZ)); if(beta<0.0) return 1e-12; //Numerator = mZstar*mZstar*sqrt(beta)*(2+(mZstar*mZstar-mZZ*mZZ+mZ*mZ)*(mZstar*mZstar-mZZ*mZZ+mZ*mZ)/(4*mZstar*mZstar*mZ*mZ)); Numerator = mZstar*mZstar*sqrt(beta)*(2+(mZstar*mZstar-mZZ*mZZ+mZ*mZ)*(mZstar*mZstar-mZZ*mZZ+mZ*mZ)/(4*mZstar*mZstar*mZ*mZ)); Denominator = 2*(mZstar*mZstar-mZ*mZ)*(mZstar*mZstar-mZ*mZ)-2*mZ*mZ*Gamma*Gamma; //mZstarDistribution = Numerator/Denominator; mZstarDistribution = Numerator/Denominator; //*/ //mZstarDistribution = 1; double acceptance = p0+p1*mZstar+p2*mZstar*mZstar; //double acceptance = 1; double mZZDistribution; /* Double_t arg= mZZ - CBmean; Double_t ret =exp(-0.5*arg*arg/(CBwidth*CBwidth)) ; //cout << "gauss x = " << x << " mean = " << mean << " sigma = " << sigma << endl ; mZZDistribution = ret ; //*/ double t = (mZZ-CBmean)/CBwidth; //std:cout << "CBwidth: " << CBwidth << " CBmean: " << CBmean << " mZZ: " << mZZ << std::endl; if (CBalpha < 0) t = -t; double absAlpha = fabs((double)CBalpha); //std::cout << "t = " << t << std::endl; if (t >= -absAlpha) { mZZDistribution = exp(-0.5*t*t); //std::cout << "Gaussian!! " << -.5*t*t << " : " << mZZDistribution << std::endl; } else { double a = TMath::Power(CBn/absAlpha,CBn)*exp(-0.5*absAlpha*absAlpha); double b = CBn/absAlpha - absAlpha; //std::cout << "power!!! a: " << a << " b: " << b << endl; mZZDistribution = a/TMath::Power(b - t, CBn); } //*/ //std::cout << "t: " << t << " absAlpha: " << absAlpha << std::endl; //double final = mZstarDistribution*acceptance; double final = mZstarDistribution*acceptance*mZZDistribution; if (final <= 0) final = 1e-12; return final; }
RooRealProxy Roo4lMasses2D::CBalpha [protected] |
Definition at line 341 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::CBmean [protected] |
Definition at line 339 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::CBn [protected] |
Definition at line 342 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::CBwidth [protected] |
Definition at line 340 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::Gamma [protected] |
Definition at line 335 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::mZ [protected] |
Definition at line 333 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::mZstar [protected] |
Definition at line 332 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::mZZ [protected] |
Definition at line 334 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::p0 [protected] |
Definition at line 336 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::p1 [protected] |
Definition at line 337 of file HZZ4LRooPdfs.h.
Referenced by evaluate().
RooRealProxy Roo4lMasses2D::p2 [protected] |
Definition at line 338 of file HZZ4LRooPdfs.h.
Referenced by evaluate().