Skip to content

Commit

Permalink
new coctails with HIJING+gamma-jet or jet-jet or EPOSLHC+gamma-jet or…
Browse files Browse the repository at this point in the history
… jet-jet
  • Loading branch information
gconesab authored and miweberSMI committed Aug 1, 2017
1 parent 05d6b82 commit 4c99aa6
Show file tree
Hide file tree
Showing 2 changed files with 154 additions and 0 deletions.
77 changes: 77 additions & 0 deletions MC/CustomGenerators/PWGGA/EPOSLHC_Pythia_GammaTriggerAndJet.C
Original file line number Diff line number Diff line change
@@ -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;
}
77 changes: 77 additions & 0 deletions MC/CustomGenerators/PWGGA/Hijing_Pythia_GammaTriggerAndJet.C
Original file line number Diff line number Diff line change
@@ -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;
}

0 comments on commit 4c99aa6

Please sign in to comment.