Fix bug in building local list of valid files.
[u/mrichter/AliRoot.git] / TEvtGen / EvtGenModels / EvtBToPlnuBK.cpp
1 //--------------------------------------------------------------------------
2 //
3 // Environment:
4 //      This software is part of the EvtGen package developed jointly
5 //      for the BaBar and CLEO collaborations.  If you use all or part
6 //      of it, please give an appropriate acknowledgement.
7 //
8 // Copyright Information: See EvtGen/COPYRIGHT
9 //      Copyright (C) 1998      Caltech, UCSB
10 //
11 // Module: EvtGenModels/EvtBToPlnuBK.cc
12 //
13 // Description: B->Xu l nu with BK (Becirevic-Kaidalov) parametrization
14 //              Xu is a pseudoscalar (pi_plus,pi0,eta or eta_prime)
15 //
16 // Modification history:
17 //
18 //    Martin Simard, U. de Montreal, 08/01/2007    Module created
19 //
20 //------------------------------------------------------------------------
21 //
22 #include "EvtGenBase/EvtPatches.hh"
23 #include <stdlib.h>
24 #include <assert.h>
25 #include "EvtGenBase/EvtParticle.hh"
26 #include "EvtGenBase/EvtGenKine.hh"
27 #include "EvtGenBase/EvtPDL.hh"
28 #include "EvtGenBase/EvtReport.hh"
29 #include "EvtGenModels/EvtBToPlnuBK.hh"
30 #include "EvtGenModels/EvtBToPlnuBKFF.hh"
31 #include "EvtGenBase/EvtSemiLeptonicScalarAmp.hh"
32
33 using std::fstream;
34 using std::cout;
35 using std::endl;
36
37
38 EvtBToPlnuBK::EvtBToPlnuBK():
39   BKmodel(0)
40   ,calcamp(0)
41 {}
42
43 EvtBToPlnuBK::~EvtBToPlnuBK(){  
44   delete BKmodel;
45   BKmodel=0;
46   delete calcamp;
47   calcamp=0;
48 }
49
50 std::string EvtBToPlnuBK::getName(){
51   
52   return "BTOPLNUBK";
53
54 }
55
56 EvtDecayBase* EvtBToPlnuBK::clone(){
57   
58   return new EvtBToPlnuBK;
59   
60 }
61
62
63 void EvtBToPlnuBK::initProbMax(){
64   
65   EvtId parnum,mesnum,lnum,nunum;
66   
67   parnum = getParentId();
68   mesnum = getDaug(0);
69   lnum = getDaug(1);
70   nunum = getDaug(2);
71   
72   double mymaxprob = calcamp->CalcMaxProb(parnum,mesnum,
73                                           lnum,nunum,BKmodel);
74   
75   setProbMax(mymaxprob);
76
77 }
78
79 void EvtBToPlnuBK::init(){
80
81   checkNDaug(3);
82
83   //We expect the parent to be a scalar 
84   //and the daughters to be X lepton neutrino
85   checkSpinParent(EvtSpinType::SCALAR);
86
87   checkSpinDaughter(1,EvtSpinType::DIRAC);
88   checkSpinDaughter(2,EvtSpinType::NEUTRINO);
89
90   EvtSpinType::spintype d1type = EvtPDL::getSpinType(getDaug(0));
91   if ( d1type==EvtSpinType::SCALAR) {
92     checkNArg(2);
93     BKmodel = new EvtBToPlnuBKFF(getArg(0),getArg(1));
94     calcamp = new EvtSemiLeptonicScalarAmp; 
95   }
96   else{
97     report(ERROR,"EvtGen") << "BK model handles only scalar meson daughters. Sorry."<<endl;
98     ::abort();
99   }
100
101 }
102
103 void  EvtBToPlnuBK::decay( EvtParticle *p ){
104
105   p->initializePhaseSpace(getNDaug(),getDaugs());
106   calcamp->CalcAmp(p,_amp2,BKmodel);
107
108 }
109
110
111