kPsiPrimeJpsiDiElectron added

Moritz Pohl

diff --git a/EVGEN/AliDecayer.h b/EVGEN/AliDecayer.h
index 46fd0bdc3193f68d2406b003b88995f1f929b74a..3fe4f37c4c535a72abe704177cb1d83a2c92c113 100644 (file)
--- a/EVGEN/AliDecayer.h
+++ b/EVGEN/AliDecayer.h
@@ -19,7 +19,7 @@ typedef enum
     kNoDecay, kHadronicD, kHadronicDWithout4Bodies, kOmega, kLambda, kPhiKK, 
     kAll, kNoDecayHeavy, kHardMuons, kBJpsi,
     kWToMuon,kWToCharm, kWToCharmToMuon, kZDiMuon, kZDiElectron, kNeutralPion, kAllMuonic,
-    kChiToJpsiGammaToMuonMuon, kChiToJpsiGammaToElectronElectron, kNoDecayBeauty
+    kChiToJpsiGammaToMuonMuon, kChiToJpsiGammaToElectronElectron, kNoDecayBeauty, kPsiPrimeJpsiDiElectron
 } Decay_t;
 #endif
 

diff --git a/EVGEN/AliGenMC.cxx b/EVGEN/AliGenMC.cxx
index 68bbef8ed9a4d50972d120529c7a99c839157833..6660454fcf0547232208b4c27aa84c53595b759f 100644 (file)
--- a/EVGEN/AliGenMC.cxx
+++ b/EVGEN/AliGenMC.cxx
@@ -162,6 +162,10 @@ void AliGenMC::Init()
        fChildSelect[0]= kProton;
         fChildSelect[1]= 211;
        break;
+    case kPsiPrimeJpsiDiElectron:
+      fChildSelect[0]= 211;
+      fChildSelect[1]= 11;
+      break;
     case kOmega:
     case kAll:
     case kAllMuonic:

diff --git a/PYTHIA6/AliDecayerPythia.cxx b/PYTHIA6/AliDecayerPythia.cxx
index 3584f44e372ae29f7821bd94d38f782766ba23b5..19e895b0dd546843b9ca277b79cffdcb43bbe93d 100644 (file)
--- a/PYTHIA6/AliDecayerPythia.cxx
+++ b/PYTHIA6/AliDecayerPythia.cxx
@@ -286,6 +286,14 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay(100553,11,2);// Upsilon'
        ForceParticleDecay(200553,11,2);// Upsilon''
        break;
+    case kPsiPrimeJpsiDiElectron:
+       products[0] =    443;
+        products[1] =    211;
+        mult[0] = 1;
+        mult[1] = 2;
+        ForceParticleDecay( 100443, products, mult, 2, 1);
+        ForceParticleDecay(    443,11,2);             
+        break; 
     case kBJpsiDiMuon:
 
        products[0] =    443;
@@ -613,7 +621,7 @@ void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t product, Int_t m
     if (norm > 0.) fBraPart[kc] /= norm;
 }
 
-void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t* products, Int_t* mult, Int_t npart)
+void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t* products, Int_t* mult, Int_t npart, Bool_t flag)
 {
 //
 //  Force decay of particle into products with multiplicity mult
@@ -637,9 +645,9 @@ void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t* products, Int_t
        for (Int_t i = 0; i < npart; i++) {
            nprod += (CountProducts(channel, products[i]) >= mult[i]);
        }
-       if (nprod) {
+       if ((nprod && !flag) || ((nprod == npart) && flag)) {
            fPythia->SetMDME(channel,1,1);
-       } else {
+       } else {                // 
            fPythia->SetMDME(channel,1,0);
            fBraPart[kc] -= fPythia->GetBRAT(channel);
        }

diff --git a/PYTHIA6/AliDecayerPythia.h b/PYTHIA6/AliDecayerPythia.h
index 6b124fbd9876a0f8a0c269e0f391a31ba4e67dcf..8869a0597b9f0f858f7d79bf30ac749ffe085f17 100644 (file)
--- a/PYTHIA6/AliDecayerPythia.h
+++ b/PYTHIA6/AliDecayerPythia.h
@@ -42,7 +42,7 @@ public AliDecayer
                    Double_t mom, Double_t theta, Double_t phi);
     Int_t    CountProducts(Int_t channel, Int_t particle);
     void     ForceParticleDecay(Int_t particle, Int_t product, Int_t mult);
-    void     ForceParticleDecay(Int_t particle, Int_t* products, Int_t* mult, Int_t npart);
+    void     ForceParticleDecay(Int_t particle, Int_t* products, Int_t* mult, Int_t npart, Bool_t flag = 0);
     void     ForceHadronicD(Int_t optUser4Bodies=1);
     void     ForceOmega();
     void     ForceLambda();

diff --git a/PYTHIA8/AliDecayerPythia8.cxx b/PYTHIA8/AliDecayerPythia8.cxx
index 1819d0aa9eddced337ac451731a67c267245e80b..2f9a4376dee9f45913716cd9084cd258ea15f506 100644 (file)
--- a/PYTHIA8/AliDecayerPythia8.cxx
+++ b/PYTHIA8/AliDecayerPythia8.cxx
@@ -540,6 +540,7 @@ void AliDecayerPythia8::ForceDecay()
     case kNoDecayHeavy:
     case kNoDecayBeauty:
     case kNeutralPion:
+    case kPsiPrimeJpsiDiElectron: 
        break;
     }
 }