Skip to content

Commit

Permalink
Move setting of weight from GeneratorConfig.C to CustomGenerator macr…
Browse files Browse the repository at this point in the history 
…o to avoid backward incompatibility
  • Loading branch information
@gconesab
gconesab committed Apr 21, 2021
1 parent 5ee8826 commit 62d035c
Show file tree
Hide file tree
Showing 2 changed files with 85 additions and 37 deletions.
74 changes: 74 additions & 0 deletions MC/CustomGenerators/PWGGA/Pythia8Weighted_GammaTriggerAndJet.C
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,74 @@
///
/// \file Pythia8Weighted_GammaTriggerAndJet.C
/// \brief Configuration of gamma-jet and jet-jet with/out decay gamma trigger
///
/// Generate PYTHIA8 gamma-jet (kPyDirectGamma) or jet-jet (kPyJets), with or without
/// triggering the gamma in one of the calorimeter possible acceptances
/// with a weight in the parton generatio to make the pT-hard distribution flatter.
/// Copy of Pythia8Weighted_GammaTriggerAndJet.C.
/// Options are:
/// * process and trigger :
/// * Pythia8Jets: jet-jet events, kPyJets (optionally with jet axis restricted to some zone or not, preferred open)
/// * Pythia8GammaJet: gamma-jet events, kPyDirectGamma, with gamma restricted to a calorimeter acceptance.
/// * 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
/// * EMCPHSDMC: EMCal+DCal+PHOS |eta| < 0.7, 250 < phi < 327 & 80 < phi < 187
///
/// This macro can only be used for AliRoot containing the setting of the weight,
/// from v5-09-54m, v5-09-55b, v5-09-56f. On master from April 2021

///
/// Main configuration method
///
AliGenerator *
GeneratorCustom
(TString opt = "kFullDetector")
{
// 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;
if (opt.EqualTo("EMCPHSDMC"))
acceptance = kCalorimeterAcceptance_EMCPHSDMC;

AliGenerator * gener = 0;

// Select the process
if (processConfig.EqualTo("Pythia8Jets")) {
gener = GeneratorPythia8Jets (kPythia8Tune_Monash2013 , acceptance); }
if (processConfig.EqualTo("Pythia8GammaJet")) {
gener = GeneratorPythia8GammaJet (kPythia8Tune_Monash2013 , acceptance);

// Rely on AliGenPythia/Plus calorimeter default acceptances or set them here
if ( acceptance == kCalorimeterAcceptance_EMCPHSDMC )
((AliGenPythiaPlus*) gener)->SetCheckBarrelCalos(kTRUE);
}
if (processConfig.EqualTo("Pythia8JetsGammaTrg")) {
gener = GeneratorPythia8JetsGammaTrg (kPythia8Tune_Monash2013 , acceptance);
// Rely on AliGenPythia/Plus calorimeter default acceptances or set them here
if ( acceptance == kCalorimeterAcceptance_EMCPHSDMC )
((AliGenPythiaPlus*) gener)->SetDecayPhotonInBarrelCalos(kTRUE);
}

// Set the desired weight
((AliGenPythiaPlus*) gener)->SetWeightPower(6.);

return gener;
}
48 changes: 11 additions & 37 deletions MC/GeneratorConfig.C
View file
Original file line number Diff line number Diff line change
Expand Up @@ -29,16 +29,10 @@ enum EGenerator_t {
kGeneratorPythia8, kGeneratorPythia8_Monash2013,
// Pythia8 jets
kGeneratorPythia8Jets, kGeneratorPythia8Jets_Monash2013,
// Pythia8 jets + ptHard weighting
kGeneratorPythia8Jets_Weighted, kGeneratorPythia8Jets_Monash2013_Weighted,
// Pythia8 jets gamma-triggered
kGeneratorPythia8JetsGammaTrg, kGeneratorPythia8JetsGammaTrg_Monash2013,
// Pythia8 jets gamma-triggered + ptHard weighting
kGeneratorPythia8JetsGammaTrg_Weighted, kGeneratorPythia8JetsGammaTrg_Monash2013_Weighted,
// Pythia8 gamma-jet
kGeneratorPythia8GammaJet, kGeneratorPythia8GammaJet_Monash2013,
// Pythia8 gamma-jet + ptHard weighting
kGeneratorPythia8GammaJet_Weighted, kGeneratorPythia8GammaJet_Monash2013_Weighted,
// Phojet
kGeneratorPhojet, kGeneratorDpmjet,
// EPOS-LHC
Expand Down Expand Up @@ -80,16 +74,10 @@ const Char_t *GeneratorName[kNGenerators] = {
"Pythia8", "Pythia8_Monash2013",
// Pythia8 jets
"Pythia8Jets", "Pythia8Jets_Monash2013",
// Pythia8 jets + ptHard weighting
"Pythia8Jets_Weighted", "Pythia8Jets_Monash2013_Weighted",
// Pythia8 jets gamma-triggered
"Pythia8JetsGammaTrg", "Pythia8JetsGammaTrg_Monash2013",
// Pythia8 jets gamma-triggered + ptHard weighting
"Pythia8JetsGammaTrg_Weighted", "Pythia8JetsGammaTrg_Monash2013_Weighted",
// Pythia8 gamma-jet
"Pythia8GammaJet", "Pythia8GammaJet_Monash2013",
// Pythia8 gamma-jet + ptHard weighting
"Pythia8GammaJet_Weighted", "Pythia8GammaJet_Monash2013_Weighted",
// Phojet / DPMjet
"Phojet", "Dpmjet",
// EPOS-LHC
Expand Down Expand Up @@ -218,9 +206,9 @@ AliGenerator *GeneratorPythia6JetsGammaTrg(Int_t tune = 0, Int_t acceptance = kC
AliGenerator *GeneratorPythia6GammaJet(Int_t tune = 0, Int_t acceptance = kCalorimeterAcceptance_FullDetector);
AliGenerator *GeneratorPythia6Heavy(Int_t process, Int_t decay, Int_t tune = 0, Bool_t HFonly = kTRUE);
AliGenerator *GeneratorPythia8(Int_t tune = 0, Int_t pdgtrig = 0, Float_t etatrig = 1.2);
AliGenerator *GeneratorPythia8Jets(Int_t tune = 0, Double_t weightPower = -1, Int_t acceptance = kCalorimeterAcceptance_FullDetector);
AliGenerator *GeneratorPythia8JetsGammaTrg(Int_t tune = 0, Double_t weightPower = -1, Int_t acceptance = kCalorimeterAcceptance_FullDetector);
AliGenerator *GeneratorPythia8GammaJet(Int_t tune = 0, Double_t weightPower = -1, Int_t acceptance = kCalorimeterAcceptance_FullDetector);
AliGenerator *GeneratorPythia8Jets(Int_t tune = 0, Int_t acceptance = kCalorimeterAcceptance_FullDetector);
AliGenerator *GeneratorPythia8JetsGammaTrg(Int_t tune = 0, Int_t acceptance = kCalorimeterAcceptance_FullDetector);
AliGenerator *GeneratorPythia8GammaJet(Int_t tune = 0, Int_t acceptance = kCalorimeterAcceptance_FullDetector);
AliGenerator *GeneratorPhojet();
AliGenerator *GeneratorEPOSLHC(Bool_t pileup = kFALSE);
AliGenerator *GeneratorHijing();
Expand Down Expand Up @@ -249,11 +237,6 @@ void GeneratorConfig(Int_t tag)

AliGenerator *gen = NULL;

double weight = -1;

// weighting exponent to produce flat pthard distribution at LHC energy
if (tag == kGeneratorPythia8Jets_Monash2013_Weighted || tag == kGeneratorPythia8JetsGammaTrg_Monash2013_Weighted || tag == kGeneratorPythia8GammaJet_Monash2013_Weighted) weight = 6;

switch (tag) {

// Default
Expand Down Expand Up @@ -292,27 +275,21 @@ void GeneratorConfig(Int_t tag)
break;

// Pythia8Jets (Monash2013)
case kGeneratorPythia8Jets:
case kGeneratorPythia8Jets_Weighted:
case kGeneratorPythia8Jets:
case kGeneratorPythia8Jets_Monash2013:
case kGeneratorPythia8Jets_Monash2013_Weighted:
gen = GeneratorPythia8Jets(kPythia8Tune_Monash2013,weight);
gen = GeneratorPythia8Jets(kPythia8Tune_Monash2013);
break;

// Pythia8JetsGammaTrg (Monash2013)
case kGeneratorPythia8JetsGammaTrg:
case kGeneratorPythia8JetsGammaTrg_Weighted:
case kGeneratorPythia8JetsGammaTrg_Monash2013:
case kGeneratorPythia8JetsGammaTrg_Monash2013_Weighted:
gen = GeneratorPythia8JetsGammaTrg(kPythia8Tune_Monash2013,weight);
gen = GeneratorPythia8JetsGammaTrg(kPythia8Tune_Monash2013);
break;

// Pythia8GammaJet (Monash2013)
case kGeneratorPythia8GammaJet:
case kGeneratorPythia8GammaJet_Weighted:
case kGeneratorPythia8GammaJet_Monash2013:
case kGeneratorPythia8GammaJet_Monash2013_Weighted:
gen = GeneratorPythia8GammaJet(kPythia8Tune_Monash2013,weight);
gen = GeneratorPythia8GammaJet(kPythia8Tune_Monash2013);
break;

// Phojet / DPMjet
Expand Down Expand Up @@ -889,7 +866,7 @@ GeneratorPythia8(Int_t tune, Int_t pdgtrig, Float_t etatrig)
/*** PYTHIA 8 JETS ***********************************************/

AliGenerator *
GeneratorPythia8Jets(Int_t tune, Double_t weight, Int_t acceptance)
GeneratorPythia8Jets(Int_t tune, Int_t acceptance)
{
//
//
Expand All @@ -900,7 +877,7 @@ GeneratorPythia8Jets(Int_t tune, Double_t weight, Int_t acceptance)
//
// jets settings
pythia->SetProcess(kPyJets);

// Careful using detector acceptance cuts and jet acceptance,
// let it open with kCalorimeterAcceptance_FullDetector
Float_t etaMax, phiMin, phiMax;
Expand All @@ -909,7 +886,6 @@ GeneratorPythia8Jets(Int_t tune, Double_t weight, Int_t acceptance)
pythia->SetJetPhiRange(phiMin, phiMax);
pythia->SetJetEtRange(0., 1000.);
pythia->SetPtHard(pthardminConfig, pthardmaxConfig); // Pt transfer of the hard scattering
if (weight > -1) pythia->SetWeightPower(weight);
pythia->SetStrucFunc(kCTEQ5L);
// quenching
pythia->SetQuench(quenchingConfig);
Expand All @@ -931,7 +907,7 @@ GeneratorPythia8Jets(Int_t tune, Double_t weight, Int_t acceptance)
/*** PYTHIA 8 JETS GAMMA-TRIGGERED ********************************/

AliGenerator *
GeneratorPythia8JetsGammaTrg(Int_t tune, Double_t weight, Int_t acceptance)
GeneratorPythia8JetsGammaTrg(Int_t tune, Int_t acceptance)
{
//
//
Expand All @@ -943,7 +919,6 @@ GeneratorPythia8JetsGammaTrg(Int_t tune, Double_t weight, Int_t acceptance)
//
// Careful with pT hard limits if triggerParticleInCalo option is on
pythia->SetTriggerParticleMinPt(pttrigminConfig);
if (weight > -1) pythia->SetWeightPower(weight);
// acceptance
Float_t etaMax, phiMin, phiMax;
GetCalorimeterAcceptance(acceptance, etaMax, phiMin, phiMax);
Expand Down Expand Up @@ -974,7 +949,7 @@ GeneratorPythia8JetsGammaTrg(Int_t tune, Double_t weight, Int_t acceptance)
/*** PYTHIA 8 GAMMA-JET ***********************************************/

AliGenerator *
GeneratorPythia8GammaJet(Int_t tune, Double_t weight, Int_t acceptance)
GeneratorPythia8GammaJet(Int_t tune, Int_t acceptance)
{
//
//
Expand All @@ -992,7 +967,6 @@ GeneratorPythia8GammaJet(Int_t tune, Double_t weight, Int_t acceptance)
// reset jets settings
pythia->SetJetEtaRange(-20., 20.); // Final state kinematic cuts
pythia->SetJetPhiRange(0., 360.);
if (weight > -1) pythia->SetWeightPower(weight);
//
return pythia;
}
Expand Down

0 comments on commit 62d035c

Please sign in to comment.