Page MenuHomeHEPForge
Files F19252053

EvtBtoXsgammaFlatEnergy.cpp
No OneTemporary
Actions

Size
2 KB
Referenced Files
None
Subscribers
None

EvtBtoXsgammaFlatEnergy.cpp
View Options

//--------------------------------------------------------------------------
//
// Environment:
// This software is part of the EvtGen package developed jointly
// for the BaBar and CLEO collaborations. If you use all or part
// of it, please give an appropriate acknowledgement.
//
// Module: EvtBtoXsgammaKagan.cc
//
// Description:
// Routine to perform two-body B->Xs,gamma decays with a fixed hadronic
// mass. For spectrum measurements.
// The input parameters are 1: the hadronic mass
// Modification history:
//
// Jim Libby October 11 2002
//------------------------------------------------------------------------
//
#include "EvtGenBase/EvtPatches.hh"
#include <stdlib.h>
#include "EvtGenModels/EvtBtoXsgamma.hh"
#include "EvtGenModels/EvtBtoXsgammaFlatEnergy.hh"
#include "EvtGenBase/EvtRandom.hh"
#include "EvtGenBase/EvtReport.hh"
#include <fstream>
using std::endl;
using std::fstream;
EvtBtoXsgammaFlatEnergy::~EvtBtoXsgammaFlatEnergy(){
}
void EvtBtoXsgammaFlatEnergy::init(int nArg, double* args){
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;
}
double EvtBtoXsgammaFlatEnergy::GetMass( int /*Xscode*/ ){
double eGamma = EvtRandom::Flat(_eRange)+_eMin;
double mH = sqrt(pow(_mB0,2)-2.0*_mB0*eGamma);
return mH;
}

File Metadata

Mime Type
text/x-c
Expires
Tue, Sep 30, 6:13 AM (11 h, 54 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
6490238
Default Alt Text
EvtBtoXsgammaFlatEnergy.cpp (2 KB)

Event Timeline

HEPForge.org · .