Changes needed for MC production of J/psi from B with radiative decays
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 18 Dec 2012 12:07:49 +0000 (12:07 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 18 Dec 2012 12:07:49 +0000 (12:07 +0000)
 Fiorella Fionda <Fiorella.Fionda@ba.infn.it>

EVGEN/AliDecayer.h
EVGEN/AliGenMC.cxx
PYTHIA6/AliDecayerPythia.cxx
TEvtGen/AliDecayerEvtGen.cxx
TEvtGen/AliGenEvtGen.cxx

index 02c4b70..0144d82 100644 (file)
@@ -17,7 +17,7 @@ typedef enum
     kBJpsiDiMuon, kBJpsiDiElectron, 
     kBPsiPrimeDiMuon, kBPsiPrimeDiElectron, kPiToMu, kKaToMu, 
     kNoDecay, kHadronicD, kHadronicDWithout4Bodies, kOmega, kLambda, kPhiKK, 
-    kAll, kNoDecayHeavy, kHardMuons, kBJpsi,
+    kAll, kNoDecayHeavy, kHardMuons, kBJpsi,  kBJpsiUndecayed,
     kWToMuon,kWToCharm, kWToCharmToMuon, kZDiMuon, kZDiElectron, kNeutralPion, kAllMuonic,
     kChiToJpsiGammaToMuonMuon, kChiToJpsiGammaToElectronElectron, kNoDecayBeauty, kPsiPrimeJpsiDiElectron,
     kElectronEM, kGammaEM, kDiElectronEM, kBeautyUpgrade
index 568c2e3..fab6219 100644 (file)
@@ -153,6 +153,7 @@ void AliGenMC::Init()
        fChildSelect[0]=kKPlus;
        break;
     case kBJpsi:
+    case kBJpsiUndecayed:
        fChildSelect[0]= 443;
        break;
    case kChiToJpsiGammaToMuonMuon:
index 7a826fe..1465fdf 100644 (file)
@@ -350,6 +350,13 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay(  531,443,1); // B_s     
        ForceParticleDecay( 5122,443,1); // Lambda_b
        break;
+    case kBJpsiUndecayed:
+        ForceParticleDecay(  511,443,1); // B0     
+        ForceParticleDecay(  521,443,1); // B+/-     
+        ForceParticleDecay(  531,443,1); // B_s     
+        ForceParticleDecay( 5122,443,1); // Lambda_b
+        fPythia->SetMDCY(fPythia->Pycomp(443),1,0); // switch-off J/psi 
+        break;
     case kBPsiPrimeDiElectron:
        ForceParticleDecay(  511,100443,1); // B0     
        ForceParticleDecay(  521,100443,1); // B+/-     
index a47bd0a..86ef8bf 100644 (file)
@@ -235,7 +235,8 @@ void AliDecayerEvtGen::ForceDecay()
      case kNeutralPion:
      case kBPsiPrimeDiElectron:
      case kBeautyUpgrade:
-       AliWarning(Form("Warning: case %d not implemented for this class!",(int)decay));
+     case kBJpsiUndecayed: 
+      AliWarning(Form("Warning: case %d not implemented for this class!",(int)decay));
      break;
      }
      ReadDecayTable();
index e3a1bfd..c1d70ad 100644 (file)
@@ -130,7 +130,7 @@ void AliGenEvtGen::Generate()
   //check if particle is already decayed by Pythia  
   if(part->GetStatusCode() != 1 || part->GetNDaughters()>0) 
     {
-    Info("AliGenEvtGen","Attention: particle %d is already decayed by Pythia!",pdg); 
+    AliDebug(1,Form("Attention: particle %d is already decayed by Pythia!",pdg));
     continue;
     }