From 65f2626cfb3f4345fec0eab080f93d7de0bf7f6a Mon Sep 17 00:00:00 2001 From: morsch Date: Tue, 20 Sep 2005 06:02:47 +0000 Subject: [PATCH] Baryon production models added by Eleazar Cuautle Flores --- PYTHIA6/AliGenPythia.cxx | 4 ++ PYTHIA6/AliPythia.cxx | 88 ++++++++++++++++++++++++++++++++++------ PYTHIA6/AliPythia.h | 5 +-- 3 files changed, 81 insertions(+), 16 deletions(-) diff --git a/PYTHIA6/AliGenPythia.cxx b/PYTHIA6/AliGenPythia.cxx index b21c95303d7..c6405e2b906 100644 --- a/PYTHIA6/AliGenPythia.cxx +++ b/PYTHIA6/AliGenPythia.cxx @@ -269,6 +269,10 @@ void AliGenPythia::Init() // Parent and Children Selection switch (fProcess) { + case kPyOldUEQ2ordered: + case kPyOldUEQ2ordered2: + case kPyOldPopcorn: + break; case kPyCharm: case kPyCharmUnforced: case kPyCharmPbPbMNR: diff --git a/PYTHIA6/AliPythia.cxx b/PYTHIA6/AliPythia.cxx index 0909aa68ceb..f5c45413dda 100644 --- a/PYTHIA6/AliPythia.cxx +++ b/PYTHIA6/AliPythia.cxx @@ -92,9 +92,71 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun // 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); -// // 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); -// -// 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; @@ -1256,3 +1317,4 @@ void AliPythia::GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4] } + diff --git a/PYTHIA6/AliPythia.h b/PYTHIA6/AliPythia.h index c71ffd6fee7..0f7c2c5d375 100644 --- a/PYTHIA6/AliPythia.h +++ b/PYTHIA6/AliPythia.h @@ -13,7 +13,8 @@ typedef enum kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma, kPyCharmPbPbMNR, kPyD0PbPbMNR, kPyDPlusPbPbMNR, kPyBeautyPbPbMNR, kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyBeautypPbMNR, - kPyCharmppMNR, kPyD0ppMNR, kPyDPlusppMNR, kPyBeautyppMNR, kPyW} + kPyCharmppMNR, kPyD0ppMNR, kPyDPlusppMNR, kPyBeautyppMNR, kPyW, + kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn} Process_t; /* typedef enum @@ -59,8 +60,6 @@ class AliPythia : public TPythia6, public AliRndm virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t k, Int_t iECMethod); virtual void Pyquen(Double_t a, Int_t ibf, Double_t b); virtual void GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]); - - // return instance of the singleton static AliPythia* Instance(); virtual void Quench(); -- 2.31.1