EvtBtoXsgammaFlatEnergy Class Reference

#include <EvtBtoXsgammaFlatEnergy.hh>

Inheritance diagram for EvtBtoXsgammaFlatEnergy:

List of all members.

Public Member Functions
	EvtBtoXsgammaFlatEnergy ()
virtual	~EvtBtoXsgammaFlatEnergy ()
void	init (int, double *)
double	GetMass (int code)
Private Attributes
double	_eMin
double	_eMax
double	_mB0
double	_eRange

---

Detailed Description

Definition at line 27 of file EvtBtoXsgammaFlatEnergy.hh.

---

Constructor & Destructor Documentation

EvtBtoXsgammaFlatEnergy::EvtBtoXsgammaFlatEnergy

(

)

[inline]

Definition at line 31 of file EvtBtoXsgammaFlatEnergy.hh.

{}

EvtBtoXsgammaFlatEnergy::~EvtBtoXsgammaFlatEnergy

(

)

[virtual]

Definition at line 30 of file EvtBtoXsgammaFlatEnergy.cc.

                                                 {
}

---

Member Function Documentation

double EvtBtoXsgammaFlatEnergy::GetMass

(

int

code

)

[virtual]

Implements EvtBtoXsgammaAbsModel.

Definition at line 64 of file EvtBtoXsgammaFlatEnergy.cc.

References _eMin, _eRange, _mB0, and EvtRandom::Flat().

                                                   {
  
  double eGamma = EvtRandom::Flat(_eRange)+_eMin;
  double mH = sqrt(pow(_mB0,2)-2.0*_mB0*eGamma); 
  return mH;
}

void EvtBtoXsgammaFlatEnergy::init	(	int	,
		double *
	)			[virtual]

Reimplemented from EvtBtoXsgammaAbsModel.

Definition at line 33 of file EvtBtoXsgammaFlatEnergy.cc.

References _eMax, _eMin, _eRange, _mB0, calibUtil::ERROR, mPi, and report().

                                                        {

  if ((nArg) > 3 || (nArg > 1 && nArg <3)){
  
  report(ERROR,"EvtGen") << "EvtBtoXsgamma generator model "
                         << "EvtBtoXsgammaFlatEnergy expected " 
                         << "either 1(default config) or two arguments but found: "<<nArg<<endl;
  report(ERROR,"EvtGen") << "Will terminate execution!"<<endl;
    ::abort();  
  }
  _mB0=5.2794;  
  double mPi = 0.140;
  double mK = 0.494;
  if(nArg == 1){
    _eMin = 1.7;
    //Invariant mass of Xsd must be greater the m_pi+m_K leads to 
    //Egamma < (m_B**2-(m_pi+m_k)**2)/(2m_B)
    _eMax = (pow(_mB0,2)-pow(mPi+mK,2))/(2.0*_mB0); 
  }else{
    _eMin=args[1];
    _eMax=args[2];
  }
  if (_eMax>(pow(_mB0,2)-pow(mPi+mK,2))/(2.0*_mB0)){
    report(ERROR,"EvtGen") << "Emax greater than Kinematic limit" << endl;
    report(ERROR,"EvtGen") << "Reset to the kinematic limit" << endl;
    report(ERROR,"EvtGen") << "(m_B**2-(m_pi+m_k)**2)/(2m_B)" << endl;
    _eMax = (pow(_mB0,2)-pow(mPi+mK,2))/(2.0*_mB0);
  }
  _eRange=_eMax-_eMin;
}

---

Member Data Documentation

double EvtBtoXsgammaFlatEnergy::_eMax [private]

Definition at line 42 of file EvtBtoXsgammaFlatEnergy.hh.

Referenced by init().

double EvtBtoXsgammaFlatEnergy::_eMin [private]

Definition at line 41 of file EvtBtoXsgammaFlatEnergy.hh.

Referenced by GetMass(), and init().

double EvtBtoXsgammaFlatEnergy::_eRange [private]

Definition at line 45 of file EvtBtoXsgammaFlatEnergy.hh.

Referenced by GetMass(), and init().

double EvtBtoXsgammaFlatEnergy::_mB0 [private]

Definition at line 44 of file EvtBtoXsgammaFlatEnergy.hh.

Referenced by GetMass(), and init().

---