Option kBJpsi added.

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

diff --git a/PYTHIA6/AliDecayerPythia.cxx b/PYTHIA6/AliDecayerPythia.cxx
index 86790e8323f860bd6557648ca4bda86007983e3c..4ac8d197931d437378d7a021587115305e1361a5 100644 (file)
--- a/PYTHIA6/AliDecayerPythia.cxx
+++ b/PYTHIA6/AliDecayerPythia.cxx
@@ -128,6 +128,7 @@ void AliDecayerPythia::ForceDecay()
     Int_t mult[2];
     Int_t products1[3];
     Int_t mult1[3];
     Int_t mult[2];
     Int_t products1[3];
     Int_t mult1[3];

+    printf("Forcedecay %d %d \n", decay, kAll);

     
     switch (decay) 
     {
     
     switch (decay) 
     {


@@ -236,6 +237,12 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay( 5122,443,1); // Lambda_b
        ForceParticleDecay(  443,11,2);  // J/Psi    
        break;
        ForceParticleDecay( 5122,443,1); // Lambda_b
        ForceParticleDecay(  443,11,2);  // J/Psi    
        break;

+    case kBJpsi:
+       ForceParticleDecay(  511,443,1); // B0     
+       ForceParticleDecay(  521,443,1); // B+/-     
+       ForceParticleDecay(  531,443,1); // B_s     
+       ForceParticleDecay( 5122,443,1); // Lambda_b
+       break;

     case kBPsiPrimeDiElectron:
        ForceParticleDecay(  511,30443,1); // B0     
        ForceParticleDecay(  521,30443,1); // B+/-     
     case kBPsiPrimeDiElectron:
        ForceParticleDecay(  511,30443,1); // B0     
        ForceParticleDecay(  521,30443,1); // B+/-