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 039f64c067324ce91d4afe9289ef830071119038..8d0f1cdaefefb2c2a4f560a03b4ecce045085d25 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 6ac8faa176ec8f03b727ede4c8523304ef0d3967..22bf4a5fb1014bf079e851ded1fdbf65d1aa3477 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 900757de1550ac41d7dc69fedf71138d60fcc924..34c3bf29e43d7caa14bf4d4718e6dc287afb6be3 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 8c9be14f70cc3ffc468beec8b34636eef22373d4..b0b126d0e730ae00d2f4d4cd77da0b52c2fee603 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 788913dec4e436d2d68ce4bbad22a27ccfc5e070..79519ff81a89e43b944cf41d074a183a93da97a3 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,