From 38a191767bb1ebc2ebdd8c55e83336ae771b060c Mon Sep 17 00:00:00 2001 From: morsch Date: Thu, 16 Jun 2011 13:46:19 +0000 Subject: [PATCH] New decay option kJpsiDimuon lmanceau@to.infn.it --- EVGEN/AliDecayer.h | 2 +- EVGEN/AliGenMC.cxx | 1 + PYTHIA6/AliDecayerPythia.cxx | 3 +++ PYTHIA8/AliDecayerPythia8.cxx | 5 +++++ 4 files changed, 10 insertions(+), 1 deletion(-) diff --git a/EVGEN/AliDecayer.h b/EVGEN/AliDecayer.h index 8f7d7522464..eb0269f256d 100644 --- a/EVGEN/AliDecayer.h +++ b/EVGEN/AliDecayer.h @@ -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, diff --git a/EVGEN/AliGenMC.cxx b/EVGEN/AliGenMC.cxx index 435d0a050b1..95a62912120 100644 --- a/EVGEN/AliGenMC.cxx +++ b/EVGEN/AliGenMC.cxx @@ -130,6 +130,7 @@ void AliGenMC::Init() case kBSemiMuonic: case kSemiMuonic: case kDiMuon: + case kJpsiDiMuon: case kBJpsiDiMuon: case kBPsiPrimeDiMuon: case kPiToMu: diff --git a/PYTHIA6/AliDecayerPythia.cxx b/PYTHIA6/AliDecayerPythia.cxx index 2b429d185fc..3e628a71b34 100644 --- a/PYTHIA6/AliDecayerPythia.cxx +++ b/PYTHIA6/AliDecayerPythia.cxx @@ -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+/- diff --git a/PYTHIA8/AliDecayerPythia8.cxx b/PYTHIA8/AliDecayerPythia8.cxx index abfde8746df..81950d2c4cc 100644 --- a/PYTHIA8/AliDecayerPythia8.cxx +++ b/PYTHIA8/AliDecayerPythia8.cxx @@ -259,6 +259,11 @@ void AliDecayerPythia8::ForceDecay() AliTPythia8::Instance()->ReadString("5332:onMode = off"); 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- -- 2.39.3