MinBias with direct gamma (kPyMbwithDirectPhoton) introduced. (Ana Marin)
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 3 Nov 2008 10:49:30 +0000 (10:49 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 3 Nov 2008 10:49:30 +0000 (10:49 +0000)
PYTHIA6/AliGenPythia.cxx
PYTHIA6/AliGenPythiaPlus.cxx
PYTHIA6/AliPythia.cxx
PYTHIA6/AliPythia6.cxx
PYTHIA6/PythiaProcesses.h

index 039f64c..8d0f1cd 100644 (file)
@@ -444,6 +444,7 @@ void AliGenPythia::Init()
        break;
     case kPyMbDefault:
     case kPyMb:
+    case kPyMbWithDirectPhoton:
     case kPyMbNonDiffr:
     case kPyMbMSEL1:
     case kPyJets:
@@ -608,6 +609,7 @@ void AliGenPythia::Generate()
        Int_t nTkbles = 0;   // Trackable particles
        if (fProcess != kPyMbDefault && 
            fProcess != kPyMb && 
+           fProcess != kPyMbWithDirectPhoton && 
            fProcess != kPyJets && 
            fProcess != kPyDirectGamma &&
            fProcess != kPyMbNonDiffr  &&
@@ -997,6 +999,7 @@ Int_t  AliGenPythia::GenerateMB()
          fProcess == kPyZ ||
          fProcess == kPyMbDefault ||
          fProcess == kPyMb ||
+         fProcess == kPyMbWithDirectPhoton ||
          fProcess == kPyMbNonDiffr)  
         && (fCutOnChild == 1) ) {
       if ( !CheckKinematicsOnChild() ) {
index 6ac8faa..22bf4a5 100644 (file)
@@ -425,6 +425,7 @@ void AliGenPythiaPlus::Init()
        break;
     case kPyMbDefault:
     case kPyMb:
+    case kPyMbWithDirectPhoton:
     case kPyMbNonDiffr:
     case kPyMbMSEL1:
     case kPyJets:
@@ -580,6 +581,7 @@ void AliGenPythiaPlus::Generate()
        Int_t nTkbles = 0;   // Trackable particles
        if (fProcess != kPyMbDefault && 
            fProcess != kPyMb && 
+           fProcess != kPyMbWithDirectPhoton && 
            fProcess != kPyJets && 
            fProcess != kPyDirectGamma &&
            fProcess != kPyMbNonDiffr  &&
@@ -929,6 +931,7 @@ Int_t  AliGenPythiaPlus::GenerateMB()
          fProcess == kPyZ ||
          fProcess == kPyMbDefault ||
          fProcess == kPyMb ||
+         fProcess == kPyMbWithDirectPhoton ||
          fProcess == kPyMbNonDiffr)  
         && (fCutOnChild == 1) ) {
       if ( !CheckKinematicsOnChild() ) {
index 900757d..34c3bf2 100644 (file)
@@ -243,6 +243,28 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
 
        AtlasTuning();
        break;
+
+    case kPyMbWithDirectPhoton:
+// Minimum Bias pp-Collisions with direct photon processes added 
+//
+//   
+//      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
+
+       SetMSUB(14,1);             //
+       SetMSUB(18,1);             //
+       SetMSUB(29,1);             //
+       SetMSUB(114,1);            //
+       SetMSUB(115,1);            //
+
+
+       AtlasTuning();
+       break;
+
     case kPyMbDefault:
 // Minimum Bias pp-Collisions
 //
index 8c9be14..b0b126d 100644 (file)
@@ -249,6 +249,28 @@ void AliPythia6::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfu
 
        AtlasTuning();
        break;
+
+    case kPyMbWithDirectPhoton:
+// Minimum Bias pp-Collisions with direct photon processes added 
+//
+//   
+//      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
+
+       SetMSUB(14,1);             //
+       SetMSUB(18,1);             //
+       SetMSUB(29,1);             //
+       SetMSUB(114,1);            //
+       SetMSUB(115,1);            //
+
+
+       AtlasTuning();
+       break;
+
     case kPyMbDefault:
 // Minimum Bias pp-Collisions
 //
index 788913d..79519ff 100644 (file)
@@ -4,7 +4,7 @@
 typedef enum
 {
     kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced,
-    kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
+    kPyJpsi, kPyJpsiChi, kPyMb, kPyMbWithDirectPhoton, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
     kPyCharmPbPbMNR, kPyD0PbPbMNR, kPyDPlusPbPbMNR, kPyDPlusStrangePbPbMNR, kPyBeautyPbPbMNR,
     kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyDPlusStrangepPbMNR, kPyBeautypPbMNR,
     kPyCharmppMNR, kPyCharmppMNRwmi, kPyD0ppMNR, kPyDPlusppMNR, kPyDPlusStrangeppMNR,