EvtBtoXsgammaFlatEnergy Class Reference

#include <EvtBtoXsgammaFlatEnergy.hh>

Inheritance diagram for EvtBtoXsgammaFlatEnergy:

Public Member Functions
	EvtBtoXsgammaFlatEnergy ()
virtual	~EvtBtoXsgammaFlatEnergy ()
void	init (int, double *)
double	GetMass (int code)
Public Member Functions inherited from EvtBtoXsgammaAbsModel
	EvtBtoXsgammaAbsModel ()
virtual	~EvtBtoXsgammaAbsModel ()

Detailed Description

Definition at line 27 of file EvtBtoXsgammaFlatEnergy.hh.

Constructor & Destructor Documentation

EvtBtoXsgammaFlatEnergy()

EvtBtoXsgammaFlatEnergy::EvtBtoXsgammaFlatEnergy ( )

inline

Definition at line 31 of file EvtBtoXsgammaFlatEnergy.hh.

{}

~EvtBtoXsgammaFlatEnergy()

EvtBtoXsgammaFlatEnergy::~EvtBtoXsgammaFlatEnergy ( )

virtual

Definition at line 30 of file EvtBtoXsgammaFlatEnergy.cc.

                                                 {
}

Member Function Documentation

GetMass()

double EvtBtoXsgammaFlatEnergy::GetMass ( int code )

virtual

Implements EvtBtoXsgammaAbsModel.

Definition at line 64 of file EvtBtoXsgammaFlatEnergy.cc.

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

init()

void EvtBtoXsgammaFlatEnergy::init ( int nArg, double * args )

virtual

Reimplemented from EvtBtoXsgammaAbsModel.

Definition at line 33 of file EvtBtoXsgammaFlatEnergy.cc.

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

---

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

• EvtBtoXsgammaFlatEnergy.hh
• EvtBtoXsgammaFlatEnergy.cc