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 30 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 69 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 32 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