]>
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 -> J/psi + 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 | #include "EvtGenModels/EvtBcPsiNPi.hh" | |
21 | #include "EvtGenBase/EvtSpinType.hh" | |
22 | ||
23 | EvtBcPsiNPi::EvtBcPsiNPi() { | |
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 | ||
28 | FA0_N=5.9; FA0_c1= 0.049; FA0_c2= 0.0015; | |
29 | FAm_N=0.0; FAm_c1=0.0; FAm_c2=0.0; | |
30 | FAp_N=-0.074; FAp_c1= 0.049; FAp_c2= 0.0015; | |
31 | FV_N=0.11; FV_c1= 0.049; FV_c2= 0.0015; | |
32 | ||
33 | } | |
34 | ||
35 | EvtBcPsiNPi::~EvtBcPsiNPi() { | |
36 | ||
37 | } | |
38 | ||
39 | std::string EvtBcPsiNPi::getName() { | |
40 | ||
41 | return "BC_PSI_NPI"; | |
42 | ||
43 | } | |
44 | ||
45 | EvtDecayBase* EvtBcPsiNPi::clone() { | |
46 | ||
47 | return new EvtBcPsiNPi; | |
48 | ||
49 | } | |
50 | ||
51 | void EvtBcPsiNPi::init() { | |
52 | ||
53 | checkNArg(0); | |
54 | ||
55 | // check spins | |
56 | checkSpinParent(EvtSpinType::SCALAR); | |
57 | checkSpinDaughter(0,EvtSpinType::VECTOR); | |
58 | // the others are scalar | |
59 | for (int i=1; i<=(getNDaug()-1);i++) { | |
60 | checkSpinDaughter(i,EvtSpinType::SCALAR); | |
61 | } | |
62 | ||
63 | } | |
64 | ||
65 | void EvtBcPsiNPi::initProbMax() { | |
66 | ||
67 | setProbMax(100.); | |
68 | if( getNDaug() == 2 ) { | |
69 | setProbMax(330.); | |
70 | } else if( getNDaug() == 3 ) { | |
71 | setProbMax(11000.); // checked with 30k events | |
72 | } else if( getNDaug() == 4 ) { | |
73 | setProbMax(36000.); | |
74 | } | |
75 | ||
76 | } |