Changes needed to use AliGenBox with TEvtGen and to switch off individual
authormorsch <andreas.morsch@cern.ch>
Tue, 17 Jun 2014 08:22:42 +0000 (10:22 +0200)
committermorsch <andreas.morsch@cern.ch>
Tue, 17 Jun 2014 08:22:42 +0000 (10:22 +0200)
decays.
F. Fionda

EVGEN/AliGenBox.cxx
PYTHIA6/AliDecayerPythia.cxx

index 34cbeac..197a8c0 100644 (file)
@@ -153,7 +153,7 @@ void AliGenBox::GenerateN(Int_t ntimes)
              TMath::Cos(2*random[0]*TMath::Pi())*
              TMath::Sqrt(-2*TMath::Log(random[1]));
        }
-       PushTrack(fTrackIt,-1,fIpart,p,origin,polar,time,kPPrimary,nt);
+       PushTrack(fTrackIt,-1,fIpart,p,origin,polar,time,kPPrimary,nt, 1., 1);
     }
 
     AliGenEventHeader* header = new AliGenEventHeader("BOX");
index 7b2b325..7119714 100644 (file)
@@ -140,6 +140,12 @@ void AliDecayerPythia::Init()
     ForceDecay();
 }
 
+void AliDecayerPythia::SwitchOffParticle(Int_t kf)
+{
+//switch off decay for particle "kf"
+fPythia->SetMDCY(fPythia->Pycomp(kf),1,0);
+}
+
 void AliDecayerPythia::Decay(Int_t idpart, TLorentzVector* p)
 {
 //  Decay a particle