#include <EvtBtoXsgammaFlatEnergy.hh>
Definition at line 27 of file EvtBtoXsgammaFlatEnergy.hh.
◆ EvtBtoXsgammaFlatEnergy()
EvtBtoXsgammaFlatEnergy::EvtBtoXsgammaFlatEnergy |
( |
| ) |
|
|
inline |
◆ ~EvtBtoXsgammaFlatEnergy()
EvtBtoXsgammaFlatEnergy::~EvtBtoXsgammaFlatEnergy |
( |
| ) |
|
|
virtual |
◆ GetMass()
double EvtBtoXsgammaFlatEnergy::GetMass |
( |
int | code | ) |
|
|
virtual |
◆ init()
void EvtBtoXsgammaFlatEnergy::init |
( |
int | nArg, |
|
|
double * | args ) |
|
virtual |
Reimplemented from EvtBtoXsgammaAbsModel.
Definition at line 33 of file EvtBtoXsgammaFlatEnergy.cc.
33 {
34
35 if ((nArg) > 3 || (nArg > 1 && nArg <3)){
36
37 report(
ERROR,
"EvtGen") <<
"EvtBtoXsgamma generator model "
38 << "EvtBtoXsgammaFlatEnergy expected "
39 << "either 1(default config) or two arguments but found: "<<nArg<<endl;
40 report(
ERROR,
"EvtGen") <<
"Will terminate execution!"<<endl;
41 ::abort();
42 }
43 _mB0=5.2794;
45 double mK = 0.494;
46 if(nArg == 1){
47 _eMin = 1.7;
48
49
50 _eMax = (pow(_mB0,2)-pow(
mPi+mK,2))/(2.0*_mB0);
51 }else{
52 _eMin=args[1];
53 _eMax=args[2];
54 }
55 if (_eMax>(pow(_mB0,2)-pow(
mPi+mK,2))/(2.0*_mB0)){
56 report(
ERROR,
"EvtGen") <<
"Emax greater than Kinematic limit" << endl;
57 report(
ERROR,
"EvtGen") <<
"Reset to the kinematic limit" << endl;
58 report(
ERROR,
"EvtGen") <<
"(m_B**2-(m_pi+m_k)**2)/(2m_B)" << endl;
59 _eMax = (pow(_mB0,2)-pow(
mPi+mK,2))/(2.0*_mB0);
60 }
61 _eRange=_eMax-_eMin;
62}
ostream & report(Severity severity, const char *facility)
The documentation for this class was generated from the following files: