Updating info for ACORDE and TRD

[u/mrichter/AliRoot.git] / PYTHIA6 / AliDecayerPythia.cxx

diff --git a/PYTHIA6/AliDecayerPythia.cxx b/PYTHIA6/AliDecayerPythia.cxx
index 891f9f82d7c6b863e58aa2ce3db0f48f69b29808..6b9a0237866471e51c61e084f63da6a77b3da399 100644 (file)
--- a/PYTHIA6/AliDecayerPythia.cxx
+++ b/PYTHIA6/AliDecayerPythia.cxx
@@ -61,7 +61,8 @@ AliDecayerPythia::AliDecayerPythia():
     fDecay(kAll),
     fHeavyFlavour(kTRUE),
     fLongLived(kFALSE),
-    fPatchOmegaDalitz(0)
+    fPatchOmegaDalitz(0),
+    fPi0(1)
 {
 // Default Constructor
     for (Int_t i=0; i< 501; i++) fBraPart[i]= 1.;
@@ -74,7 +75,8 @@ AliDecayerPythia::AliDecayerPythia(const AliDecayerPythia &decayer):
     fDecay(kAll),
     fHeavyFlavour(kTRUE),
     fLongLived(kFALSE),
-    fPatchOmegaDalitz(0)
+    fPatchOmegaDalitz(0),
+    fPi0(1)
 {
     // Copy Constructor
     decayer.Copy(*this);
@@ -118,8 +120,8 @@ void AliDecayerPythia::Init()
        fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
     }
 */
+    if (fPi0) fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
 
-    fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
     Int_t isw = 0;
     if (fLongLived) isw = 1;
     
@@ -692,8 +694,7 @@ void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t product, Int_t m
     Double_t norm = 0.;
     for (Int_t channel=ifirst; channel<=ilast;channel++){
       norm+=fPythia->GetBRAT(channel);
-      printf("%f ",fPythia->GetBRAT(channel));
-    }printf("\n");
+    }
     if (norm < 1.-1.e-12 || norm > 1.+1.e-12) {
         char pName[16];
         fPythia->Pyname(particle,pName);