LambdaBLambda force decay option added.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 10 May 2012 15:15:37 +0000 (15:15 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 10 May 2012 15:15:37 +0000 (15:15 +0000)
EVGEN/AliDecayer.h
EVGEN/AliGenMC.cxx
PYTHIA6/AliDecayerPythia.cxx
PYTHIA8/AliDecayerPythia8.cxx
TEvtGen/AliDecayerEvtGen.cxx

index eb0269f..8701604 100644 (file)
@@ -20,7 +20,7 @@ typedef enum
     kAll, kNoDecayHeavy, kHardMuons, kBJpsi,
     kWToMuon,kWToCharm, kWToCharmToMuon, kZDiMuon, kZDiElectron, kNeutralPion, kAllMuonic,
     kChiToJpsiGammaToMuonMuon, kChiToJpsiGammaToElectronElectron, kNoDecayBeauty, kPsiPrimeJpsiDiElectron,
-    kElectronEM, kGammaEM, kDiElectronEM
+    kElectronEM, kGammaEM, kDiElectronEM, kLambdaBLambda
 } Decay_t;
 #endif
 
index 95a6291..9d47827 100644 (file)
@@ -180,6 +180,7 @@ void AliGenMC::Init()
     case kNoDecayHeavy:
     case kNoDecayBeauty:
     case kNeutralPion:
+    case kLambdaBLambda:  
        break;
     }
 
index 0e7e5c9..999b3c9 100644 (file)
@@ -429,6 +429,10 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay(  331,22,1); // etaprime     
        ForceParticleDecay(  333,22,1); // phi     
        break;
+    case kLambdaBLambda:
+        ForceParticleDecay( 5122, 4122, 1);
+        ForceParticleDecay( 4122, 3312, 1);
+      break;
     }
 }
 
index fd3e8a7..dd957d7 100644 (file)
@@ -564,6 +564,12 @@ void AliDecayerPythia8::ForceDecay()
        AliTPythia8::Instance()->ReadString("3122:onMode = off");
        AliTPythia8::Instance()->ReadString("3122:onIfAll = 2212 211 ");
        break;
+    case kLambdaBLambda:
+        AliTPythia8::Instance()->ReadString("5122:onMode = off");
+        AliTPythia8::Instance()->ReadString("4122:onMode = off");
+       AliTPythia8::Instance()->ReadString("5122:onIfAll = 4122");
+       AliTPythia8::Instance()->ReadString("4122:onIfAll = 3122");
+       break;
     case kAll:
        break;
     case kNoDecay:
index efd158f..d3eead7 100644 (file)
@@ -233,6 +233,7 @@ void AliDecayerEvtGen::ForceDecay()
      case kNoDecayHeavy:
      case kNeutralPion:
      case kBPsiPrimeDiElectron:
+     case kLambdaBLambda:
      AliWarning(Form("Warning: case %s not implemented for this class!",(Char_t)decay));
      break;
      }