Baryon production models added by Eleazar Cuautle Flores <ecuautle@nucleares.unam.mx>
[u/mrichter/AliRoot.git] / PYTHIA6 / AliPythia.cxx
index 0909aa68ceb9fc5c4d3bebeb442c0b05b64a48d1..f5c45413dda07416cc47fb28b70ede96e0c46612 100644 (file)
@@ -92,9 +92,71 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
 //  select charm production
     switch (process) 
     {
 //  select charm production
     switch (process) 
     {
+    case kPyOldUEQ2ordered:  //Old underlying events with Q2 ordered QCD processes
+//        Multiple interactions on.
+       SetMSTP(81,1);
+// Double Gaussian matter distribution.
+       SetMSTP(82,4);
+       SetPARP(83,0.5);
+       SetPARP(84,0.4);
+//  pT0.
+       SetPARP(82,2.0);
+//  Reference energy for pT0 and energy rescaling pace.
+       SetPARP(89,1800);
+       SetPARP(90,0.25);
+//  String drawing almost completely minimizes string length.
+       SetPARP(85,0.9);
+       SetPARP(86,0.95);
+// ISR and FSR activity.
+       SetPARP(67,4);
+       SetPARP(71,4);
+// Lambda_FSR scale.
+       SetPARJ(81,0.29);
+       break;
+    case kPyOldUEQ2ordered2:   
+// Old underlying events with Q2 ordered QCD processes
+// Multiple interactions on.
+       SetMSTP(81,1);
+// Double Gaussian matter distribution.
+       SetMSTP(82,4);
+       SetPARP(83,0.5);
+       SetPARP(84,0.4);
+// pT0.
+       SetPARP(82,2.0);
+// Reference energy for pT0 and energy rescaling pace.
+       SetPARP(89,1800);
+       SetPARP(90,0.16);  // here is the difference with  kPyOldUEQ2ordered
+// String drawing almost completely minimizes string length.
+       SetPARP(85,0.9);
+       SetPARP(86,0.95);
+// ISR and FSR activity.
+       SetPARP(67,4);
+       SetPARP(71,4);
+// Lambda_FSR scale.
+       SetPARJ(81,0.29);       
+       break;
+    case kPyOldPopcorn:  
+// Old production mechanism: Old Popcorn
+       SetMSEL(1);
+       SetMSTJ(12,3); 
+// (D=2) Like MSTJ(12)=2 but added prod ofthe 1er rank baryon
+       SetMSTP(88,2); 
+// (D=1)see can be used to form  baryons (BARYON JUNCTION)
+       SetMSTJ(1,1);  
+       SetMSTP(51,kCTEQ5L);// CTEQ 5L        ! CTEQ5L pdf
+       SetMSTP(81,1);      // Multiple Interactions ON
+       SetMSTP(82,4);      // Double Gaussian Model         
+       SetPARP(82,1.8);    // [GeV]    PT_min at Ref. energy
+       SetPARP(89,1000.);  // [GeV]   Ref. energy
+       SetPARP(90,0.16);   // 2*epsilon (exponent in power law)
+       SetPARP(83,0.5);    // Core density in proton matter dist. (def.value)
+       SetPARP(84,0.5);    // Core radius
+       SetPARP(85,0.33);   // Regulates gluon prod. mechanism
+       SetPARP(86,0.66);   // Regulates gluon prod. mechanism
+       SetPARP(67,1);      // Regulate gluon prod. mechanism
+       break;
     case kPyCharm:
        SetMSEL(4);
     case kPyCharm:
        SetMSEL(4);
-//
 //  heavy quark masses
 
        SetPMAS(4,1,1.2);
 //  heavy quark masses
 
        SetPMAS(4,1,1.2);
@@ -203,18 +265,17 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
 //  QCD Jets
 //
        SetMSEL(1);
 //  QCD Jets
 //
        SetMSEL(1);
-//
-// Pythia Tune A (CDF)
-//
-       SetPARP(67,4.);            // Regulates Initial State Radiation
-       SetMSTP(82,4);             // Double Gaussian Model
-       SetPARP(82,2.0);           // [GeV]    PT_min at Ref. energy
-       SetPARP(84,0.4);           // Core radius
-       SetPARP(85,0.90) ;         // Regulates gluon prod. mechanism
-       SetPARP(86,0.95);          // Regulates gluon prod. mechanism
-       SetPARP(89,1800.);         // [GeV]   Ref. energy
-       SetPARP(90,0.25);          // 2*epsilon (exponent in power law)
-       break;
+ // Pythia Tune A (CDF)
+ //
+       SetPARP(67,4.);            // Regulates Initial State Radiation
+       SetMSTP(82,4);             // Double Gaussian Model
+       SetPARP(82,2.0);           // [GeV]    PT_min at Ref. energy
+       SetPARP(84,0.4);           // Core radius
+       SetPARP(85,0.90) ;         // Regulates gluon prod. mechanism
+       SetPARP(86,0.95);          // Regulates gluon prod. mechanism
+       SetPARP(89,1800.);         // [GeV]   Ref. energy
+       SetPARP(90,0.25);          // 2*epsilon (exponent in power law)
+       break;
     case kPyDirectGamma:
        SetMSEL(10);
        break;
     case kPyDirectGamma:
        SetMSEL(10);
        break;
@@ -1256,3 +1317,4 @@ void AliPythia::GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]
 
 }
 
 
 }
 
+