Force MSTP(52) = 2 for nucl. pfd
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 13 Aug 2013 14:33:15 +0000 (14:33 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 13 Aug 2013 14:33:15 +0000 (14:33 +0000)
PYTHIA6/AliPythia.cxx

index 918e67c..9b5e681 100644 (file)
@@ -20,6 +20,7 @@
 #include "AliFastGlauber.h"
 #include "AliQuenchingWeights.h"
 #include "AliOmegaDalitz.h"
+#include "AliLog.h"
 #include "TVector3.h"
 #include "TLorentzVector.h"
 #include "PyquenCommon.h"
@@ -643,8 +644,15 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
 //  Initialize PYTHIA
 //
 //  Select the tune
-    if (itune > -1) Pytune(itune);
-    
+    if (itune > -1) {
+      Pytune(itune);
+      if (GetMSTP(192) > 1 || GetMSTP(193) > 1) {
+       AliWarning(Form("Structure function for tune %5d set to %5s\n", 
+                       itune,  AliStructFuncType::PDFsetName(strucfunc).Data()));
+       SetMSTP(52,2);
+       SetMSTP(51, AliStructFuncType::PDFsetIndex(strucfunc));
+      }
+    }
 //  
     SetMSTP(41,1);   // all resonance decays switched on
     if (process == kPyJetsPWHG || process == kPyCharmPWHG || process == kPyBeautyPWHG) {