kNoDecayBeauty added.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 16 Feb 2010 22:27:15 +0000 (22:27 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 16 Feb 2010 22:27:15 +0000 (22:27 +0000)
EVGEN/AliDecayer.h
EVGEN/AliGenMC.cxx
PYTHIA6/AliDecayerPythia.cxx
PYTHIA6/AliDecayerPythia.h
PYTHIA8/AliDecayerPythia8.cxx

index b00d075..46fd0bd 100644 (file)
@@ -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
+    kChiToJpsiGammaToMuonMuon, kChiToJpsiGammaToElectronElectron, kNoDecayBeauty
 } Decay_t;
 #endif
 
index cefaa40..c35e1fd 100644 (file)
@@ -165,6 +165,7 @@ void AliGenMC::Init()
     case kAllMuonic:
     case kNoDecay:
     case kNoDecayHeavy:
+    case kNoDecayBeauty:
     case kNeutralPion:
        break;
     }
index bcfd04b..3584f44 100644 (file)
@@ -380,6 +380,9 @@ void AliDecayerPythia::ForceDecay()
     case kNoDecayHeavy:
     case kNeutralPion:
        break;
+    case kNoDecayBeauty:                                                    
+        SwitchOffBDecay(); 
+       break;
     }
 }
 
@@ -747,6 +750,14 @@ void  AliDecayerPythia::ForceLambda()
 }
 
 
+void AliDecayerPythia::SwitchOffBDecay(){
+  Int_t heavyB[]={511,521,531,5122,5132,5232,5332};
+  for(int i=0;i<4;i++)
+    {
+      fPythia->SetMDCY(fPythia->Pycomp(heavyB[i]),1,0);
+    }
+}
+
 Float_t  AliDecayerPythia::GetPartialBranchingRatio(Int_t kf)
 {
 // Get branching ratio
index 45cd884..6b124fb 100644 (file)
@@ -32,6 +32,7 @@ public AliDecayer
     virtual void    DecayLongLivedParticles()  {fLongLived    = kTRUE;}
     virtual Float_t GetPartialBranchingRatio(Int_t ipart);
     virtual Float_t GetLifetime(Int_t kf);
+    virtual void    SwitchOffBDecay();
     virtual void    WriteDecayTable();
     virtual void    ReadDecayTable();
     
index b14a1a1..1819d0a 100644 (file)
@@ -538,6 +538,7 @@ void AliDecayerPythia8::ForceDecay()
        TPythia8::Instance()->ReadString("HadronLevel:Decay = off");
        break;
     case kNoDecayHeavy:
+    case kNoDecayBeauty:
     case kNeutralPion:
        break;
     }