Bug fix in strangeness vs multiplicity code. Quite stupid bug but will not impact...

[u/mrichter/AliRoot.git] / PYTHIA8 / AliPythia8.cxx

diff --git a/PYTHIA8/AliPythia8.cxx b/PYTHIA8/AliPythia8.cxx
index 2e072d51694d87d846e74e707965827b8927ab0d..9477b3e89b34874f79a024820afd3fb0b1fb38df 100644 (file)
--- a/PYTHIA8/AliPythia8.cxx
+++ b/PYTHIA8/AliPythia8.cxx
@@ -104,8 +104,8 @@ void AliPythia8::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfu
     fProcess   = process;
     fEcms      = energy;
     fStrucFunc = strucfunc;
-//...Switch off decay of pi0, K0S, Lambda, Sigma+-, Xi0-, Omega-.
-    ReadString("111:mayDecay  = off");
+    ReadString("111:mayDecay  = on");
+//...Switch off decay of K0L, Lambda, Sigma+-, Xi0-, Omega-.
     ReadString("310:mayDecay  = off");
     ReadString("3122:mayDecay = off");
     ReadString("3112:mayDecay = off");
@@ -115,11 +115,9 @@ void AliPythia8::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfu
     ReadString("3322:mayDecay = off");
     ReadString("3334:mayDecay = off");
     // Select structure function 
-
-    ReadString("PDF:useLHAPDF = on");
-    ReadString(Form("PDF:LHAPDFset = %s", AliStructFuncType::PDFsetName(fStrucFunc).Data()));
-
-    // Particles produced in string fragmentation point directly to either of the two endpoints
+          ReadString("PDF:useLHAPDF = on");
+         ReadString(Form("PDF:LHAPDFset = %s", AliStructFuncType::PDFsetName(fStrucFunc).Data()));
+      // Particles produced in string fragmentation point directly to either of the two endpoints
     // of the string (depending in the side they were generated from).
 
 //    SetMSTU(16,2); // ????
@@ -243,6 +241,7 @@ void AliPythia8::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfu
        ReadString("SoftQCD:doubleDiffractive = on");
        ReadString("SoftQCD:doubleDiffractive = on");
        if (tune > -1) ReadString(Form("Tune:pp = %3d", tune));
+       ReadString("PDF:useLHAPDF = off");
        break;
     case kPyLhwgMb:
 // Les Houches Working Group 05 Minimum Bias pp-Collisions: hep-ph/0604120
@@ -511,7 +510,28 @@ void AliPythia8::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfu
        // With parton showers on we are generating "Z inclusive process"
        ReadString("PartonLevel:ISR = on");
        ReadString("PartonLevel:FSR = on");
-       break;  
+       break;
+    case kPyZgamma:
+        //Inclusive production of Z/gamma*
+        //f fbar -> Z/gamma
+        ReadString("WeakSingleBoson:ffbar2gmZ = on");
+        // Initial/final parton shower on (Pythia default)
+        // With parton showers on we are generating "Z inclusive process"
+        ReadString("PartonLevel:ISR = on");
+        ReadString("PartonLevel:FSR = on");
+        break;
+    case kPyMBRSingleDiffraction:
+      ReadString("Diffraction:PomFlux = 5");
+      ReadString("SoftQCD:singleDiffractive = on");
+      break;
+    case kPyMBRDoubleDiffraction:
+      ReadString("Diffraction:PomFlux = 5");
+      ReadString("SoftQCD:doubleDiffractive = on");
+      break;
+    case kPyMBRCentralDiffraction: 
+      ReadString("Diffraction:PomFlux = 5");
+      ReadString("SoftQCD:centralDiffractive = on");
+      break;
     case kPyMbWithDirectPhoton:
     case kPyBeautyJets:
     case kPyMbAtlasTuneMC09: