Fix bug in building local list of valid files.
[u/mrichter/AliRoot.git] / TEvtGen / EvtGenModels / EvtISGW.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: EvtISGW.cc
12 //
13 // Description: Routine to implement semileptonic decays according
14 //              to the model ISGW
15 //
16 // Modification history:
17 //
18 //    DJL/RYD     September 25, 1996         Module created
19 //
20 //------------------------------------------------------------------------
21 // 
22 #include "EvtGenBase/EvtPatches.hh"
23 #include <stdlib.h>
24 #include "EvtGenBase/EvtParticle.hh"
25 #include "EvtGenBase/EvtGenKine.hh"
26 #include "EvtGenBase/EvtPDL.hh"
27 #include "EvtGenBase/EvtReport.hh"
28 #include "EvtGenModels/EvtISGW.hh"
29 #include <string>
30 #include "EvtGenModels/EvtISGWFF.hh"
31 #include "EvtGenBase/EvtSemiLeptonicScalarAmp.hh"
32 #include "EvtGenBase/EvtSemiLeptonicVectorAmp.hh"
33 #include "EvtGenBase/EvtSemiLeptonicTensorAmp.hh"
34
35 EvtISGW::~EvtISGW() {}
36
37 std::string EvtISGW::getName(){
38
39   return "ISGW";     
40
41 }
42
43
44 EvtDecayBase* EvtISGW::clone(){
45
46   return new EvtISGW;
47
48 }
49
50 void EvtISGW::decay( EvtParticle *p ){
51
52
53   p->initializePhaseSpace(getNDaug(),getDaugs());
54
55   calcamp->CalcAmp(p,_amp2,isgwffmodel);
56   return;
57
58 }
59
60
61 void EvtISGW::init(){
62
63   checkNArg(0);
64   checkNDaug(3);
65
66
67   //We expect the parent to be a scalar 
68   //and the daughters to be X lepton neutrino
69
70   EvtSpinType::spintype mesontype=EvtPDL::getSpinType(getDaug(0));
71
72   checkSpinParent(EvtSpinType::SCALAR);
73   checkSpinDaughter(1,EvtSpinType::DIRAC);
74   checkSpinDaughter(2,EvtSpinType::NEUTRINO);
75
76
77   isgwffmodel = new EvtISGWFF;
78   
79   if ( mesontype==EvtSpinType::SCALAR ) { 
80     calcamp = new EvtSemiLeptonicScalarAmp; 
81   }
82   if ( mesontype==EvtSpinType::VECTOR ) { 
83     calcamp = new EvtSemiLeptonicVectorAmp; 
84   }
85   if ( mesontype==EvtSpinType::TENSOR ) { 
86     calcamp = new EvtSemiLeptonicTensorAmp; 
87   }
88   
89 }
90