From 4c99aa6410aa500dcfd500ef8f0ff34d8c333889 Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Date: Mon, 31 Jul 2017 16:04:45 +0200 Subject: [PATCH] new coctails with HIJING+gamma-jet or jet-jet or EPOSLHC+gamma-jet or jet-jet --- .../PWGGA/EPOSLHC_Pythia_GammaTriggerAndJet.C | 77 +++++++++++++++++++ .../PWGGA/Hijing_Pythia_GammaTriggerAndJet.C | 77 +++++++++++++++++++ 2 files changed, 154 insertions(+) create mode 100644 MC/CustomGenerators/PWGGA/EPOSLHC_Pythia_GammaTriggerAndJet.C create mode 100644 MC/CustomGenerators/PWGGA/Hijing_Pythia_GammaTriggerAndJet.C diff --git a/MC/CustomGenerators/PWGGA/EPOSLHC_Pythia_GammaTriggerAndJet.C b/MC/CustomGenerators/PWGGA/EPOSLHC_Pythia_GammaTriggerAndJet.C new file mode 100644 index 00000000..f287ba24 --- /dev/null +++ b/MC/CustomGenerators/PWGGA/EPOSLHC_Pythia_GammaTriggerAndJet.C @@ -0,0 +1,77 @@ +/// +/// \file EPOSLHC_Pythia_GammaTriggerAndJet.C +/// \brief Configuration EPOSLHC plus a gamma-jet or jet-jet with/out decay gamma trigger +/// +/// Generate PYTHIA8 or PYTHIA6 gamma-jet (kPyDirectGamma) or jet-jet (kPyJets), with or without +/// triggering the gamma in one of the calorimeter possible acceptances on top of a EPOSLHC event +/// Options are: +/// * process and trigger : +/// * Pythia6Jets/Pythia8Jets: jet-jet events, kPyJets (optionally with jet axis restricted to some zone or not, preferred open) +/// * Pythia6GammaJet/Pythia8GammaJet: gamma-jet events, kPyDirectGamma, with gamma restricted to a calorimeter acceptance. +/// * Pythia6JetsTrg/Pythia8JetsGammaTrg: jet-jet events, kPyJets, with a decay gamma to a calorimeter acceptance above a given pT. +/// +/// * acceptance of photon : +/// * FullDetector: Minimal acceptance cuts +/// * EMCRun1 : EMCal |eta| < 0.70, 80 < phi < 180 +/// * EMCRun2 : EMCal |eta| < 0.70, 80 < phi < 187 +/// * PHSRun1 : PHOS |eta| < 0.13, 260 < phi < 320 +/// * PHSRun2 : PHOS |eta| < 0.13, 250 < phi < 320 +/// * PHSDMC : DCal+PHOS |eta| < 0.7, 250 < phi < 327 + +/// +/// Main configuration method +/// +AliGenerator * +GeneratorCustom +(TString opt = "kFullDetector") +{ + // Init cocktail + AliGenCocktail *ctl = GeneratorCocktail(Form("EPOSLHC_%s",processConfig.Data())); + + // + // EPOSLHC + // + AliGenerator *epos = GeneratorEPOSLHC(); + ctl->AddGenerator(epos, "EPOSLHC", 1.); + + // + // PYTHIA + // + // acceptance + Int_t acceptance = kCalorimeterAcceptance_FullDetector; + if (opt.EqualTo("FullDetector")) + acceptance = kCalorimeterAcceptance_FullDetector; + if (opt.EqualTo("EMCRun1")) + acceptance = kCalorimeterAcceptance_EMCRun1; + if (opt.EqualTo("PHSRun1")) + acceptance = kCalorimeterAcceptance_PHSRun1; + if (opt.EqualTo("EMCRun2")) + acceptance = kCalorimeterAcceptance_EMCRun2; + if (opt.EqualTo("PHSRun2")) + acceptance = kCalorimeterAcceptance_PHSRun2; + if (opt.EqualTo("PHSDMC")) + acceptance = kCalorimeterAcceptance_PHSDMC; + + // process + AliGenerator *gammajet = 0; + + // PYTHIA6 + if (processConfig.EqualTo("Pythia6Jets")) + gammajet = GeneratorPythia6Jets (kPythia6Tune_Perugia2011, acceptance); + if (processConfig.EqualTo("Pythia6GammaJet")) + gammajet = GeneratorPythia6GammaJet (kPythia6Tune_Perugia2011, acceptance); + if (processConfig.EqualTo("Pythia6JetsGammaTrg")) + gammajet = GeneratorPythia6JetsGammaTrg (kPythia6Tune_Perugia2011, acceptance); + + // PYTHIA8 + if (processConfig.EqualTo("Pythia8Jets")) + gammajet = GeneratorPythia8Jets (kPythia8Tune_Monash2013 , acceptance); + if (processConfig.EqualTo("Pythia8GammaJet")) + gammajet = GeneratorPythia8GammaJet (kPythia8Tune_Monash2013 , acceptance); + if (processConfig.EqualTo("Pythia8JetsGammaTrg")) + gammajet = GeneratorPythia8JetsGammaTrg (kPythia8Tune_Monash2013 , acceptance); + + ctl->AddGenerator(gammajet, processConfig, 1.); + + return ctl; +} diff --git a/MC/CustomGenerators/PWGGA/Hijing_Pythia_GammaTriggerAndJet.C b/MC/CustomGenerators/PWGGA/Hijing_Pythia_GammaTriggerAndJet.C new file mode 100644 index 00000000..23f27907 --- /dev/null +++ b/MC/CustomGenerators/PWGGA/Hijing_Pythia_GammaTriggerAndJet.C @@ -0,0 +1,77 @@ +/// +/// \file Hijing_Pythia_GammaTriggerAndJet.C +/// \brief Configuration HIJING plus a gamma-jet or jet-jet with/out decay gamma trigger +/// +/// Generate PYTHIA8 or PYTHIA6 gamma-jet (kPyDirectGamma) or jet-jet (kPyJets), with or without +/// triggering the gamma in one of the calorimeter possible acceptances on top of a HIJING event +/// Options are: +/// * process and trigger : +/// * Pythia6Jets/Pythia8Jets: jet-jet events, kPyJets (optionally with jet axis restricted to some zone or not, preferred open) +/// * Pythia6GammaJet/Pythia8GammaJet: gamma-jet events, kPyDirectGamma, with gamma restricted to a calorimeter acceptance. +/// * Pythia6JetsTrg/Pythia8JetsGammaTrg: jet-jet events, kPyJets, with a decay gamma to a calorimeter acceptance above a given pT. +/// +/// * acceptance of photon : +/// * FullDetector: Minimal acceptance cuts +/// * EMCRun1 : EMCal |eta| < 0.70, 80 < phi < 180 +/// * EMCRun2 : EMCal |eta| < 0.70, 80 < phi < 187 +/// * PHSRun1 : PHOS |eta| < 0.13, 260 < phi < 320 +/// * PHSRun2 : PHOS |eta| < 0.13, 250 < phi < 320 +/// * PHSDMC : DCal+PHOS |eta| < 0.7, 250 < phi < 327 + +/// +/// Main configuration method +/// +AliGenerator * +GeneratorCustom +(TString opt = "kFullDetector") +{ + // Init cocktail + AliGenCocktail *ctl = GeneratorCocktail(Form("HIJING_%s",processConfig.Data())); + + // + // HIJING + // + AliGenerator *hij = GeneratorHijing(); + ctl->AddGenerator(hij, "Hijing", 1.); + + // + // PYTHIA + // + // acceptance + Int_t acceptance = kCalorimeterAcceptance_FullDetector; + if (opt.EqualTo("FullDetector")) + acceptance = kCalorimeterAcceptance_FullDetector; + if (opt.EqualTo("EMCRun1")) + acceptance = kCalorimeterAcceptance_EMCRun1; + if (opt.EqualTo("PHSRun1")) + acceptance = kCalorimeterAcceptance_PHSRun1; + if (opt.EqualTo("EMCRun2")) + acceptance = kCalorimeterAcceptance_EMCRun2; + if (opt.EqualTo("PHSRun2")) + acceptance = kCalorimeterAcceptance_PHSRun2; + if (opt.EqualTo("PHSDMC")) + acceptance = kCalorimeterAcceptance_PHSDMC; + + // process + AliGenerator *gammajet = 0; + + // PYTHIA6 + if (processConfig.EqualTo("Pythia6Jets")) + gammajet = GeneratorPythia6Jets (kPythia6Tune_Perugia2011, acceptance); + if (processConfig.EqualTo("Pythia6GammaJet")) + gammajet = GeneratorPythia6GammaJet (kPythia6Tune_Perugia2011, acceptance); + if (processConfig.EqualTo("Pythia6JetsGammaTrg")) + gammajet = GeneratorPythia6JetsGammaTrg (kPythia6Tune_Perugia2011, acceptance); + + // PYTHIA8 + if (processConfig.EqualTo("Pythia8Jets")) + gammajet = GeneratorPythia8Jets (kPythia8Tune_Monash2013 , acceptance); + if (processConfig.EqualTo("Pythia8GammaJet")) + gammajet = GeneratorPythia8GammaJet (kPythia8Tune_Monash2013 , acceptance); + if (processConfig.EqualTo("Pythia8JetsGammaTrg")) + gammajet = GeneratorPythia8JetsGammaTrg (kPythia8Tune_Monash2013 , acceptance); + + ctl->AddGenerator(gammajet, processConfig, 1.); + + return ctl; +}