Use MSTU(16) = 2 for heavy flavor. HF Hadrons point to the

correct endpoint of the string.

diff --git a/PYTHIA6/AliPythia.cxx b/PYTHIA6/AliPythia.cxx
index d3a7c439870509346fbd6514b67e1f19accfc2aa..e98836f31a74c5f3e2ffb04ad88d8d12e27574c1 100644 (file)
--- a/PYTHIA6/AliPythia.cxx
+++ b/PYTHIA6/AliPythia.cxx
@@ -84,6 +84,7 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
     //  select structure function 
     SetMSTP(52,2);
     SetMSTP(51,strucfunc);
+    
 //
 // Pythia initialisation for selected processes//
 //
@@ -184,6 +185,7 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
        break;
     case kPyCharmUnforced:
        SetMSEL(0);
+       SetMSTU(16,2);
 // gq->qg   
        SetMSUB(28,1);
 // gg->qq
@@ -192,6 +194,7 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
        SetMSUB(68,1);
        break;
     case kPyBeautyUnforced:
+       SetMSTU(16,2);
        SetMSEL(0);
 // gq->qg   
        SetMSUB(28,1);
@@ -1181,6 +1184,7 @@ void AliPythia::ConfigHeavyFlavor()
     // QCD scales
     SetMSTP(32,2);
     SetPARP(34,1.0);
+    SetMSTU(16,2);
 }
 
 void AliPythia::AtlasTuning()