Settings for Pythia6 with Tune D6T added.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 11 Nov 2009 09:04:58 +0000 (09:04 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 11 Nov 2009 09:04:58 +0000 (09:04 +0000)
prod/fpprod/Config.C

index fcf7620..22656c2 100644 (file)
 
 enum PDC06Proc_t 
 {
-  kPythia6, kPhojet, kRunMax
+  kPythia6, kPythia6D6T, kPhojet, kRunMax
 };
 
 const char * pprRunName[] = {
-  "kPythia6", "kPhojet"
+    "kPythia6", "kPythia6D6T", "kPhojet" 
 };
 
 enum Mag_t
@@ -70,6 +70,7 @@ const char * pprField[] = {
 //--- Functions ---
 class AliGenPythia;
 AliGenerator *MbPythia();
+AliGenerator *MbPythiaTuneD6T();
 AliGenerator *MbPhojet();
 void ProcessEnvironmentVars();
 
@@ -100,7 +101,11 @@ void Config()
 #if defined(__CINT__)
   gSystem->Load("liblhapdf");      // Parton density functions
   gSystem->Load("libEGPythia6");   // TGenerator interface
-  gSystem->Load("libpythia6");     // Pythia
+  if (proc != kPythia6D6T) {
+      gSystem->Load("libpythia6");     // Pythia 6.2
+  } else {
+      gSystem->Load("libqpythia");     // Pythia 6.4
+  }
   gSystem->Load("libAliPythia6");  // ALICE specific implementations
   gSystem->Load("libgeant321");
 #endif
@@ -186,6 +191,8 @@ void Config()
   
   if (proc == kPythia6) {
       gener = MbPythia();
+  } else if (proc == kPythia6D6T) {
+      gener = MbPythiaTuneD6T();
   } else if (proc == kPhojet) {
       gener = MbPhojet();
   }
@@ -465,6 +472,26 @@ AliGenerator* MbPythia()
       return pythia;
 }
 
+AliGenerator* MbPythiaTuneD6T()
+{
+      comment = comment.Append(" pp at 14 TeV: Pythia low-pt");
+//
+//    Pythia
+      AliGenPythia* pythia = new AliGenPythia(-1); 
+      pythia->SetMomentumRange(0, 999999.);
+      pythia->SetThetaRange(0., 180.);
+      pythia->SetYRange(-12.,12.);
+      pythia->SetPtRange(0,1000.);
+      pythia->SetProcess(kPyMb);
+      pythia->SetEnergyCMS(energy);
+//    Tune
+//    109     D6T : Rick Field's CDF Tune D6T (NB: needs CTEQ6L pdfs externally)
+      pythia->SetTune(109);
+      pythia->SetStrucFunc(kCTEQ6l);
+//
+      return pythia;
+}
+
 AliGenerator* MbPhojet()
 {
       comment = comment.Append(" pp at 14 TeV: Pythia low-pt");