Fix for definitions for CINT
[u/mrichter/AliRoot.git] / TEvtGen / EvtGenModels / EvtBcBsNPi.cpp
1 //--------------------------------------------------------------------------
2 //
3 // Environment:
4 //      This software is part of the EvtGen package. If you use all or part
5 //      of it, please give an appropriate acknowledgement.
6 //
7 // Copyright Information: See EvtGen/COPYRIGHT
8 //
9 // Module: EvtGenModels/EvtBcBsNPi.hh
10 //
11 // Description: Decay model for Bc -> Bs + npi
12 //
13 // Modification history:
14 //
15 //    A.Berezhnoy, A.Likhoded, A.Luchinsky  July 2011   Module created
16 //
17 //------------------------------------------------------------------------
18
19 #include "EvtGenBase/EvtPatches.hh"
20 #include "EvtGenBase/EvtSpinType.hh"
21 #include "EvtGenModels/EvtBcBsNPi.hh"
22
23 EvtBcBsNPi::EvtBcBsNPi() {
24
25   _beta=-0.108; _mRho=0.775; _gammaRho=0.149;
26   _mRhopr=1.364; _gammaRhopr=0.400; _mA1=1.23; _gammaA1=0.4;
27   //            Fp_N=1.3; Fp_c1=0.30; Fp_c2=0.069;
28   Fp_N=3*1.3; Fp_c1=0.30; Fp_c2=0.069;
29   Fm_N=0.0; Fm_c1=0.0; Fm_c2=0.0;
30
31 }
32
33 EvtBcBsNPi::~EvtBcBsNPi() {
34 }
35
36 std::string EvtBcBsNPi::getName() {
37
38   return "BC_BS_NPI";
39
40 }
41
42 EvtDecayBase* EvtBcBsNPi::clone() {
43
44   return new EvtBcBsNPi;
45
46 }
47
48 void EvtBcBsNPi::init() {
49
50   checkNArg(0);
51
52   // check spins
53   checkSpinParent(EvtSpinType::SCALAR);
54   checkSpinDaughter(0,EvtSpinType::SCALAR);
55   // the others are scalar
56   for (int i=1; i<=(getNDaug()-1);i++) {
57     checkSpinDaughter(i,EvtSpinType::SCALAR);
58   }
59
60 }
61
62 void EvtBcBsNPi::initProbMax() {
63
64   if ( getNDaug() == 2 ) {
65     setProbMax(250.);
66   } else if ( getNDaug() == 3 ) {
67     setProbMax(25000.);// checked at 30k events
68   } else if( getNDaug() == 4 ) {
69     setProbMax(45000.); // checked at 30k events
70   }
71
72 }