From 0c0f2010032e9162100599e7359b2c40561d3406 Mon Sep 17 00:00:00 2001
From: morsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 31 May 2006 12:27:22 +0000
Subject: [PATCH] Use MSTU(16) = 2 for heavy flavor. HF Hadrons point to the
 correct endpoint of the string.

---
 PYTHIA6/AliPythia.cxx | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/PYTHIA6/AliPythia.cxx b/PYTHIA6/AliPythia.cxx
index d3a7c439870..e98836f31a7 100644
--- 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()
-- 
2.39.3