New decay option kJpsiDimuon
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 16 Jun 2011 13:46:19 +0000 (13:46 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 16 Jun 2011 13:46:19 +0000 (13:46 +0000)
lmanceau@to.infn.it

EVGEN/AliDecayer.h
EVGEN/AliGenMC.cxx
PYTHIA6/AliDecayerPythia.cxx
PYTHIA8/AliDecayerPythia8.cxx

index 8f7d752..eb0269f 100644 (file)
@@ -13,7 +13,7 @@ typedef TVirtualMCDecayer AliDecayer;
 #if ROOT_VERSION_CODE >= 197633  //Corresponds to Root v3-04-01
 typedef enum 
 {
-    kBSemiElectronic, kSemiElectronic, kDiElectron, kBSemiMuonic, kSemiMuonic, kDiMuon,
+    kBSemiElectronic, kSemiElectronic, kDiElectron, kBSemiMuonic, kSemiMuonic, kDiMuon, kJpsiDiMuon,
     kBJpsiDiMuon, kBJpsiDiElectron, 
     kBPsiPrimeDiMuon, kBPsiPrimeDiElectron, kPiToMu, kKaToMu, 
     kNoDecay, kHadronicD, kHadronicDWithout4Bodies, kOmega, kLambda, kPhiKK, 
index 435d0a0..95a6291 100644 (file)
@@ -130,6 +130,7 @@ void AliGenMC::Init()
     case kBSemiMuonic:
     case kSemiMuonic:
     case kDiMuon:
+    case kJpsiDiMuon:
     case kBJpsiDiMuon:
     case kBPsiPrimeDiMuon:
     case kPiToMu:
index 2b429d1..3e628a7 100644 (file)
@@ -250,6 +250,9 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay(100553,13,2);// Upsilon'
        ForceParticleDecay(200553,13,2);// Upsilon''
        break;
+    case kJpsiDiMuon:
+       ForceParticleDecay(  443,13,2); // J/Psi
+       break;
     case kBSemiElectronic:
        ForceParticleDecay(  511,11,1); // B0     
        ForceParticleDecay(  521,11,1); // B+/-     
index abfde87..81950d2 100644 (file)
@@ -260,6 +260,11 @@ void AliDecayerPythia8::ForceDecay()
        AliTPythia8::Instance()->ReadString("5332:onIfAny = 13");
 
        break;
+    case kJpsiDiMuon:
+//      J/Psi-> mu+ mu-
+       AliTPythia8::Instance()->ReadString("443:onMode = off");
+       AliTPythia8::Instance()->ReadString("443:onIfAll = 13 13");
+       break;
     case kDiMuon:
 //      Rho -> mu+ mu-
        AliTPythia8::Instance()->ReadString("113:onMode = off");