Fix bug in building local list of valid files.
[u/mrichter/AliRoot.git] / TEvtGen / EvtGenModels / EvtBcPsiNPi.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 -> 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 }