Converting TEvtGen to native cmake
[u/mrichter/AliRoot.git] / TEvtGen / EvtGen / EvtGenBase / EvtDecayIncoherent.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: EvtGen/EvtDecayIncoherent.cc
12 //
13 // Description:
14 //
15 // Modification history:
16 //
17 //    DJL/RYD     August 11, 1998         Module created
18 //
19 //------------------------------------------------------------------------
20 #include "EvtGenBase/EvtPatches.hh"
21
22 #include "EvtGenBase/EvtDecayBase.hh"
23 #include "EvtGenBase/EvtDecayIncoherent.hh"
24 #include "EvtGenBase/EvtParticle.hh"
25 #include "EvtGenBase/EvtRadCorr.hh"
26 #include "EvtGenBase/EvtReport.hh"
27 #include "EvtGenBase/EvtPDL.hh"
28
29
30 void EvtDecayIncoherent::makeDecay(EvtParticle* p, bool recursive){
31
32   //initialize this the hard way..
33   //Lange June 26, 2000
34   for (size_t i=0; i<static_cast<unsigned int>(MAX_DAUG); i++ ) { 
35     spinDensitySet[i]=0;
36   }
37
38   _daugsDecayedByParentModel=false;
39
40   decay(p);
41   p->setDecayProb(1.0);
42
43   EvtSpinDensity rho;
44
45   rho.setDiag(p->getSpinStates());
46
47   p->setSpinDensityBackward(rho);
48
49   if (getPHOTOS() || EvtRadCorr::alwaysRadCorr()) {
50     EvtRadCorr::doRadCorr(p);
51   }
52
53   if(!recursive) return;
54
55   //Now decay the daughters.
56
57   if ( !daugsDecayedByParentModel()) {
58     
59     for(size_t i=0;i<p->getNDaug();i++){
60       //Need to set the spin density of the daughters to be
61       //diagonal.
62       rho.setDiag(p->getDaug(i)->getSpinStates());
63       //if (p->getDaug(i)->getNDaug()==0){
64       //only do this if the user has not already set the 
65       //spin density matrix herself.
66       //Lange June 26, 2000
67       if ( isDaughterSpinDensitySet(i)==0 ) { 
68         p->getDaug(i)->setSpinDensityForward(rho);
69       }
70       else{
71         //report(INFO,"EvtGen") << "spinDensitymatrix already set!!!\n";
72         EvtSpinDensity temp=p->getDaug(i)->getSpinDensityForward();
73         //      report(INFO,"EvtGen") <<temp<<endl;
74       }
75       //Now decay the daughter.  Really!
76       p->getDaug(i)->decay();
77     } 
78   }
79                             
80 }
81
82
83
84
85
86
87