Min. bias with the default Pythia settings
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 20 Apr 2007 14:59:04 +0000 (14:59 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 20 Apr 2007 14:59:04 +0000 (14:59 +0000)
PYTHIA6/AliGenPythia.cxx
PYTHIA6/AliPythia.cxx
PYTHIA6/AliPythia.h

index bc9447f..c845edc 100644 (file)
@@ -425,6 +425,7 @@ void AliGenPythia::Init()
     case kPyJpsi:
        fParentSelect[0] = 443;
        break;
+    case kPyMbDefault:
     case kPyMb:
     case kPyMbNonDiffr:
     case kPyMbMSEL1:
@@ -586,12 +587,16 @@ void AliGenPythia::Generate()
        Int_t nc = 0;        // Total n. of selected particles
        Int_t nParents = 0;  // Selected parents
        Int_t nTkbles = 0;   // Trackable particles
-       if (fProcess != kPyMb && fProcess != kPyJets && 
+       if (fProcess != kPyMbDefault && 
+           fProcess != kPyMb && 
+           fProcess != kPyJets && 
            fProcess != kPyDirectGamma &&
            fProcess != kPyMbNonDiffr  &&
            fProcess != kPyMbMSEL1     &&
-           fProcess != kPyW && fProcess != kPyZ &&
-           fProcess != kPyCharmppMNRwmi && fProcess != kPyBeautyppMNRwmi) {
+           fProcess != kPyW && 
+           fProcess != kPyZ &&
+           fProcess != kPyCharmppMNRwmi && 
+           fProcess != kPyBeautyppMNRwmi) {
            
            for (i = 0; i < np; i++) {
                TParticle* iparticle = (TParticle *) fParticles->At(i);
@@ -899,7 +904,11 @@ Int_t  AliGenPythia::GenerateMB()
     }
 
     //Introducing child cuts in case kPyW, kPyZ, kPyMb, and kPyMbNonDiff
-    if ( (fProcess == kPyW || fProcess == kPyZ || fProcess == kPyMb || fProcess == kPyMbNonDiffr)  
+    if ( (fProcess == kPyW ||
+         fProcess == kPyZ ||
+         fProcess == kPyMbDefault ||
+         fProcess == kPyMb ||
+         fProcess == kPyMbNonDiffr)  
         && (fCutOnChild == 1) ) {
       if ( !CheckKinematicsOnChild() ) {
        delete[] pParent;
index fa015ca..6520a13 100644 (file)
@@ -237,6 +237,18 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
        SetMSUB(95,1);             // low pt production
 
        AtlasTuning();
+       break;
+    case kPyMbDefault:
+// Minimum Bias pp-Collisions
+//
+//   
+//      select Pythia min. bias model
+       SetMSEL(0);
+       SetMSUB(92,1);             // single diffraction AB-->XB
+       SetMSUB(93,1);             // single diffraction AB-->AX
+       SetMSUB(94,1);             // double diffraction
+       SetMSUB(95,1);             // low pt production
+
        break;
     case kPyLhwgMb:
 // Les Houches Working Group 05 Minimum Bias pp-Collisions: hep-ph/0604120
index 01a670f..9a18dac 100644 (file)
@@ -17,7 +17,7 @@ typedef enum
  kPyCharmppMNR, kPyCharmppMNRwmi, kPyD0ppMNR, kPyDPlusppMNR, kPyDPlusStrangeppMNR, 
  kPyBeautyppMNR, kPyBeautyppMNRwmi, kPyW, kPyZ, kPyMbMSEL1,
  kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn,
- kPyLhwgMb}
+ kPyLhwgMb,kPyMbDefault}
 Process_t;
 
 class AliFastGlauber;