CMS 3D CMS Logo

Public Member Functions | Protected Member Functions | Protected Attributes

Roo4lMasses2D Class Reference

#include <HZZ4LRooPdfs.h>

List of all members.

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

Detailed Description

Definition at line 310 of file HZZ4LRooPdfs.h.


Constructor & Destructor Documentation

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.

{ }

Member Function Documentation

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;
}

Member Data Documentation

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().