]>
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/EvtBcBsStarNPi.hh | |
10 | // | |
11 | // Description: Decay model for Bc -> Bs* + npi | |
12 | // | |
13 | // Modification history: | |
14 | // | |
15 | // A.Berezhnoy, A.Likhoded, A.Luchinsky April 2011 Module created | |
16 | // | |
17 | //------------------------------------------------------------------------ | |
18 | ||
19 | #include "EvtGenBase/EvtPatches.hh" | |
20 | ||
21 | #include "EvtGenModels/EvtBcBsStarNPi.hh" | |
22 | #include "EvtGenBase/EvtSpinType.hh" | |
23 | ||
24 | EvtBcBsStarNPi::EvtBcBsStarNPi() { | |
25 | ||
26 | _beta=-0.108; _mRho=0.775; _gammaRho=0.149; | |
27 | _mRhopr=1.364; _gammaRhopr=0.400; _mA1=1.23; _gammaA1=0.4; | |
28 | ||
29 | FA0_N=8.1; FA0_c1=0.30; FA0_c2=0.069; | |
30 | FAm_N=0.0; FAm_c1=0.0; FAm_c2=0.0; | |
31 | FAp_N=0.15; FAp_c1=0.30; FAp_c2=0.069; | |
32 | FV_N= 1.08; FV_c1=0.30; FV_c2=0.069; | |
33 | ||
34 | } | |
35 | ||
36 | EvtBcBsStarNPi::~EvtBcBsStarNPi() { | |
37 | ||
38 | } | |
39 | ||
40 | std::string EvtBcBsStarNPi::getName() { | |
41 | ||
42 | return "BC_BSSTAR_NPI"; | |
43 | ||
44 | } | |
45 | ||
46 | EvtDecayBase* EvtBcBsStarNPi::clone() { | |
47 | ||
48 | return new EvtBcBsStarNPi; | |
49 | ||
50 | } | |
51 | ||
52 | void EvtBcBsStarNPi::init() { | |
53 | ||
54 | checkNArg(0); | |
55 | ||
56 | // check spins | |
57 | checkSpinParent(EvtSpinType::SCALAR); | |
58 | checkSpinDaughter(0,EvtSpinType::VECTOR); | |
59 | // the others are scalar | |
60 | for (int i=1; i<=(getNDaug()-1);i++) { | |
61 | checkSpinDaughter(i,EvtSpinType::SCALAR); | |
62 | } | |
63 | ||
64 | } | |
65 | ||
66 | void EvtBcBsStarNPi::initProbMax() { | |
67 | ||
68 | if ( getNDaug() == 2 ) { | |
69 | setProbMax(100.); | |
70 | } else if( getNDaug() == 3 ) { | |
71 | setProbMax(40000.); | |
72 | } else if( getNDaug() == 4 ) { | |
73 | setProbMax(620.); // checked, 30k events | |
74 | } | |
75 | ||
76 | } |