Use MSTU(16) = 2 for heavy flavor. HF Hadrons point to the
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 31 May 2006 12:27:22 +0000 (12:27 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 31 May 2006 12:27:22 +0000 (12:27 +0000)
correct endpoint of the string.

PYTHIA6/AliPythia.cxx

index d3a7c43..e98836f 100644 (file)
@@ -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()