TEvtGen/EvtGenModels/EvtBtoXsgammaFixedMass.cxx

   1 //--------------------------------------------------------------------------
   2 //
   3 // Environment:
   4 //      This software is part of the EvtGen package developed jointly
   5 //      for the BaBar and CLEO collaborations.  If you use all or part
   6 //      of it, please give an appropriate acknowledgement.
   7 //
   8 // Module: EvtBtoXsgammaKagan.cc
   9 //
  10 // Description:
  11 //       Routine to perform two-body B->Xs,gamma decays with a fixed hadronic
  12 //       mass. For spectrum measurements.
  13 //       The input parameters are 1: the hadronic mass
  14
  15 // Modification history:
  16 //
  17 //      Jim Libby October 11 2002
  18 //------------------------------------------------------------------------
  19 //
  20
  21 #include "EvtGenBase/EvtPatches.hh"
  22 #include <stdlib.h>
  23 #include "EvtGenModels/EvtBtoXsgamma.hh"
  24 #include "EvtGenModels/EvtBtoXsgammaFixedMass.hh"
  25 #include "EvtGenBase/EvtReport.hh"
  26 #include <fstream>
  27 using std::endl;
  28 using std::fstream;
  29
  30 EvtBtoXsgammaFixedMass::~EvtBtoXsgammaFixedMass(){
  31 }
  32
  33 void EvtBtoXsgammaFixedMass::init(int nArg, double* args){
  34
  35   if ((nArg) > 2 || (nArg > 1 && nArg <2)){
  36
  37   report(ERROR,"EvtGen") << "EvtBtoXsgamma generator model "
  38                          << "EvtBtoXsgammaFixedMass 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
  44   if(nArg == 1){
  45     _mH = 2.0;
  46   }else{
  47     _mH=args[1];
  48   }
  49 }
  50
  51 double EvtBtoXsgammaFixedMass::GetMass( int /*Xscode*/ ){
  52   return _mH;
  53 }
  54