]>
Commit | Line | Data |
---|---|---|
0ca57c2f | 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 | } |