From 9fa374770f40639e2c84cf9f29276f2c630b425d Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Balbastre Date: Mon, 29 Mar 2021 10:40:39 +0200 Subject: [PATCH 1/8] select the parton checking the acceptance so that at least one of the parton/jets is within the acceptance, not the two --- MC/config/PWGGAJE/hooks/jets_hook.C | 46 ++++++++++++++++------------- 1 file changed, 26 insertions(+), 20 deletions(-) diff --git a/MC/config/PWGGAJE/hooks/jets_hook.C b/MC/config/PWGGAJE/hooks/jets_hook.C index a8e031fa3..7ada5cebe 100644 --- a/MC/config/PWGGAJE/hooks/jets_hook.C +++ b/MC/config/PWGGAJE/hooks/jets_hook.C @@ -7,8 +7,8 @@ R__ADD_INCLUDE_PATH($O2DPG_ROOT) /// \brief Select jet events within acceptance or associated parton flavor using Pythia Hooks. /// -/// Select outoging jets on the 2->2 process, at least one in a selected acceptance and -/// with a given PDG value. +/// Select outoging parton/jets on the 2->2 process, at least one in a selected acceptance and +/// optionally select the parton with a given PDG value. /// Only valid for PYTHIA8 and using Hooks /// /// \author Gustavo Conesa Balbastre (LPSC-IN2P3-CNRS) @@ -22,6 +22,14 @@ class UserHooks_jets : public Pythia8::UserHooks ~UserHooks_jets() = default; bool canVetoPartonLevel() override { return true; }; bool doVetoPartonLevel(const Pythia8::Event& event) override { + +// for (Int_t id = 0; id < 10; id++) +// { +// printf("parton %d PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n", +// id, event[id].id() , event[id].status(), event[id].mother1(), +// event[id].e() , event[id].pT(), +// event[id].eta(), event[id].phi()*TMath::RadToDeg()); +// } // Get the outgoing 2->2 partons. // The jets are in position 5 or 6. @@ -37,17 +45,6 @@ class UserHooks_jets : public Pythia8::UserHooks if ( mOutPartonPDG > 0 && TMath::Abs(event[id1].id()) != mOutPartonPDG ) okpdg1 = false; - if ( acc1 && okpdg1 ) - { - printf("--- Accepted event, jet 1 ---\n"); - printf("PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n", - event[id1].id() , event[id1].status(), event[id1].mother1(), - event[id1].e() , event[id1].pT(), - event[id1].eta(), event[id1].phi()*TMath::RadToDeg()); - - return false; - } - // Check the second jet // bool acc2 = detector_acceptance(mAcceptance, event[id2].phi(), event[id2].eta()) ; @@ -55,21 +52,30 @@ class UserHooks_jets : public Pythia8::UserHooks if ( mOutPartonPDG > 0 && TMath::Abs(event[id2].id()) != mOutPartonPDG ) okpdg2 = false; - - if ( acc2 && okpdg2 ) + + //printf("acc1 %d, acc2 %d, okpdg1 %d, okpdf2 %d\n",acc1,acc2,okpdg1,okpdg2); + + if ( (acc1 || acc2) && (okpdg1 || okpdg2) ) { - printf("--- Accepted event, jet 2 ---\n"); - printf("PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n", + printf("--- Accepted event ---\n"); + printf("\t --- jet 1 ---\n"); + printf("\t PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n", + event[id1].id() , event[id1].status(), event[id1].mother1(), + event[id1].e() , event[id1].pT(), + event[id1].eta(), event[id1].phi()*TMath::RadToDeg()); + + printf("\t --- jet 2 ---\n"); + printf("\t PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n", event[id2].id() , event[id2].status(), event[id2].mother1(), event[id2].e() , event[id2].pT(), event[id2].eta(), event[id2].phi()*TMath::RadToDeg()); - + return false; } - // Jets not found + // Jets rejected // - printf("--- Rejected event ---\n"); + printf("\t --- Rejected event ---\n"); return true; From 7e87aff4861f722619b1fea9accd82fcf964d09a Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Balbastre Date: Mon, 29 Mar 2021 10:51:42 +0200 Subject: [PATCH 2/8] add option to set predefined particle acceptance --- MC/bin/o2dpg_sim_workflow.py | 2 ++ MC/config/PWGGAJE/hooks/jets_hook.C | 4 ++-- MC/config/PWGGAJE/hooks/prompt_gamma_hook.C | 4 ++-- MC/config/PWGGAJE/trigger/decay_gamma_jets.C | 4 ++-- MC/config/PWGGAJE/trigger/prompt_gamma.C | 4 ++-- 5 files changed, 10 insertions(+), 8 deletions(-) diff --git a/MC/bin/o2dpg_sim_workflow.py b/MC/bin/o2dpg_sim_workflow.py index fc2ed7dea..35b148b20 100755 --- a/MC/bin/o2dpg_sim_workflow.py +++ b/MC/bin/o2dpg_sim_workflow.py @@ -34,6 +34,7 @@ parser.add_argument('-ptTrigMin',help='generated pT trigger minimum', default=0) parser.add_argument('-ptTrigMax',help='generated pT trigger maximum', default=-1) +parser.add_argument('-acceptance',help='select particles within predefined acceptance in ${O2DPG_ROOT}/MC/run/common/detector_acceptance.C', default=0) parser.add_argument('--embedding',action='store_true', help='With embedding into background') parser.add_argument('-nb',help='number of background events / timeframe', default=20) @@ -175,6 +176,7 @@ def getDPL_global_options(bigshm=False,nosmallrate=False): PTTRIGMIN=float(args.ptTrigMin) PTTRIGMAX=float(args.ptTrigMax) + PARTICLE_ACCEPTANCE=int(args.acceptance) ## Pt Hat productions WEIGHTPOW=int(args.weightPow) diff --git a/MC/config/PWGGAJE/hooks/jets_hook.C b/MC/config/PWGGAJE/hooks/jets_hook.C index 7ada5cebe..6f469bfdb 100644 --- a/MC/config/PWGGAJE/hooks/jets_hook.C +++ b/MC/config/PWGGAJE/hooks/jets_hook.C @@ -103,9 +103,9 @@ pythia8_userhooks_jets(int acc = 0, int pdgPar = 0) printf("Select outgoing partons with pdg = %d\n",pdgPar); } - if ( !acc && gSystem->Getenv("CONFIG_DETECTOR_ACCEPTANCE") ) + if ( !acc && gSystem->Getenv("PARTICLE_ACCEPTANCE") ) { - acc = atoi(gSystem->Getenv("CONFIG_DETECTOR_ACCEPTANCE")); + acc = atoi(gSystem->Getenv("PARTICLE_ACCEPTANCE")); printf("Requested acceptance %d\n",acc); } diff --git a/MC/config/PWGGAJE/hooks/prompt_gamma_hook.C b/MC/config/PWGGAJE/hooks/prompt_gamma_hook.C index dadf3dc75..72f6168b2 100644 --- a/MC/config/PWGGAJE/hooks/prompt_gamma_hook.C +++ b/MC/config/PWGGAJE/hooks/prompt_gamma_hook.C @@ -121,9 +121,9 @@ Pythia8::UserHooks* printf("Select outgoing partons with pdg = %d\n",pdgPar); } - if ( !acc && gSystem->Getenv("CONFIG_DETECTOR_ACCEPTANCE") ) + if ( !acc && gSystem->Getenv("PARTICLE_ACCEPTANCE") ) { - acc = atoi(gSystem->Getenv("CONFIG_DETECTOR_ACCEPTANCE")); + acc = atoi(gSystem->Getenv("PARTICLE_ACCEPTANCE")); printf("Requested acceptance %d\n",acc); } diff --git a/MC/config/PWGGAJE/trigger/decay_gamma_jets.C b/MC/config/PWGGAJE/trigger/decay_gamma_jets.C index 7b0ea58e8..587d33f5f 100644 --- a/MC/config/PWGGAJE/trigger/decay_gamma_jets.C +++ b/MC/config/PWGGAJE/trigger/decay_gamma_jets.C @@ -27,8 +27,8 @@ o2::eventgen::Trigger decay_gamma_jets( ) // Select photons within acceptance // Int_t acceptance = 0; - if ( gSystem->Getenv("CONFIG_DETECTOR_ACCEPTANCE") ) - acceptance = atoi(gSystem->Getenv("CONFIG_DETECTOR_ACCEPTANCE")); + if ( gSystem->Getenv("PARTICLE_ACCEPTANCE") ) + acceptance = atoi(gSystem->Getenv("PARTICLE_ACCEPTANCE")); //printf("Requested acceptance %d\n",acceptance); // Particle loop diff --git a/MC/config/PWGGAJE/trigger/prompt_gamma.C b/MC/config/PWGGAJE/trigger/prompt_gamma.C index 7104db36a..81050fbdf 100644 --- a/MC/config/PWGGAJE/trigger/prompt_gamma.C +++ b/MC/config/PWGGAJE/trigger/prompt_gamma.C @@ -71,8 +71,8 @@ o2::eventgen::Trigger prompt_gamma( ) // Select photons within acceptance // Int_t acceptance = 0; - if ( gSystem->Getenv("CONFIG_DETECTOR_ACCEPTANCE") ) - acceptance = atoi(gSystem->Getenv("CONFIG_DETECTOR_ACCEPTANCE")); + if ( gSystem->Getenv("PARTICLE_ACCEPTANCE") ) + acceptance = atoi(gSystem->Getenv("PARTICLE_ACCEPTANCE")); //printf("Requested acceptance %d\n",acceptance); if ( detector_acceptance(acceptance, gamma.Phi(),gamma.Eta()) ) From fd8939296a3b01230029a2f3d4695b958b1e22b6 Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Balbastre Date: Mon, 29 Mar 2021 11:13:16 +0200 Subject: [PATCH 3/8] fix Pb PDG value and add configuration for Pythia HI depending on PDG --- MC/config/common/pythia8/utils/mkpy8cfg.py | 8 ++++++++ 1 file changed, 8 insertions(+) diff --git a/MC/config/common/pythia8/utils/mkpy8cfg.py b/MC/config/common/pythia8/utils/mkpy8cfg.py index 25500bb53..49e2a5278 100755 --- a/MC/config/common/pythia8/utils/mkpy8cfg.py +++ b/MC/config/common/pythia8/utils/mkpy8cfg.py @@ -94,6 +94,14 @@ fout.write('PromptPhoton:all = on \n') fout.write('\n') +### heavy ion settings (valid for Pb-Pb 5520 only) +if args.idA==1000822080 and args.idB==1000822080: + fout.write('### heavy-ion settings (valid for Pb-Pb 5520 only) \n') + fout.write('HeavyIon:SigFitNGen = 0 \n') + fout.write('HeavyIon:SigFitDefPar = 13.88,1.84,0.22,0.0,0.0,0.0,0.0,0.0 \n') + fout.write('HeavyIon:bWidth = 14.48 \n') +fout.write('\n') + ### decays fout.write('### decays \n') fout.write('ParticleDecays:limitTau0 = on \n') ### we will need to put some parameters for these settings From 527d6a2fa38377f76fb76e2dd202aacc8e32667a Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Balbastre Date: Mon, 29 Mar 2021 11:13:41 +0200 Subject: [PATCH 4/8] fix Pb PDG value and add configuration for Pythia HI depending on PDG --- MC/bin/o2dpg_sim_workflow.py | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/MC/bin/o2dpg_sim_workflow.py b/MC/bin/o2dpg_sim_workflow.py index 35b148b20..f551eaa04 100755 --- a/MC/bin/o2dpg_sim_workflow.py +++ b/MC/bin/o2dpg_sim_workflow.py @@ -229,16 +229,16 @@ def getDPL_global_options(bigshm=False,nosmallrate=False): PDGB=2212 # proton if COLTYPE == 'PbPb': - PDGA=2212 # Pb???? #---> to be checked (seems same as pp case) - PDGB=2212 # Pb???? + PDGA=1000822080 # Pb + PDGB=1000822080 # Pb if COLTYPE == 'pPb': - PDGA=2212 # proton - PDGB=2212 # Pb???? + PDGA=2212 # proton + PDGB=1000822080 # Pb if COLTYPE == 'Pbp': - PDGA=2212 # Pb???? - PDGB=2212 # proton + PDGA=1000822080 # Pb + PDGB=2212 # proton # produce the signal configuration SGN_CONFIG_task=createTask(name='gensgnconf_'+str(tf), tf=tf, cwd=timeframeworkdir) From 5512edeeb6b4fa5f27c64d7bae5b700cfb125a51 Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Balbastre Date: Mon, 29 Mar 2021 11:54:03 +0200 Subject: [PATCH 5/8] Set default ecms to -1, so that we can automatically set the energy for Pb-Pb to 5.02, but add an error message if not set --- MC/bin/o2dpg_sim_workflow.py | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) diff --git a/MC/bin/o2dpg_sim_workflow.py b/MC/bin/o2dpg_sim_workflow.py index f551eaa04..99c50b19a 100755 --- a/MC/bin/o2dpg_sim_workflow.py +++ b/MC/bin/o2dpg_sim_workflow.py @@ -25,7 +25,7 @@ parser.add_argument('-ini',help='generator init parameters file, for example: ${O2DPG_ROOT}/MC/config/PWGHF/ini/GeneratorHF.ini', default='') parser.add_argument('-confKey',help='generator or trigger configuration key values, for example: GeneratorPythia8.config=pythia8.cfg', default='') -parser.add_argument('-eCMS',help='CMS energy', default=5200.0) +parser.add_argument('-eCMS',help='CMS energy', default=-1) parser.add_argument('-col',help='collision sytem: pp, PbPb, pPb, Pbp, ...', default='pp') parser.add_argument('-ptHatBin',help='pT hard bin number', default=-1) parser.add_argument('-ptHatMin',help='pT hard minimum when no bin requested', default=0) @@ -221,7 +221,6 @@ def getDPL_global_options(bigshm=False,nosmallrate=False): PTHATMAX=hig_edge[PTHATBIN] # translate here collision type to PDG - # not sure this is what we want to do (GCB) COLTYPE=args.col if COLTYPE == 'pp': @@ -231,6 +230,9 @@ def getDPL_global_options(bigshm=False,nosmallrate=False): if COLTYPE == 'PbPb': PDGA=1000822080 # Pb PDGB=1000822080 # Pb + if ECMS < 0: # assign 5.02 TeV to Pb-Pb + print('>>> Set CM Energy to PbPb case 5.02 TeV') + ECMS=5020.0 if COLTYPE == 'pPb': PDGA=2212 # proton @@ -240,6 +242,10 @@ def getDPL_global_options(bigshm=False,nosmallrate=False): PDGA=1000822080 # Pb PDGB=2212 # proton + if ECMS < 0: + print('Error: Collision Energy not set!!!') + exit(1) + # produce the signal configuration SGN_CONFIG_task=createTask(name='gensgnconf_'+str(tf), tf=tf, cwd=timeframeworkdir) if GENERATOR == 'pythia8': From e31d1dc3e8bc2b191b4a65dd0a5fd5ae57f0852e Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Balbastre Date: Mon, 29 Mar 2021 14:42:49 +0200 Subject: [PATCH 6/8] add beams energies as parameter --- MC/bin/o2dpg_sim_workflow.py | 23 ++++++++++++++++++---- MC/config/common/pythia8/utils/mkpy8cfg.py | 23 +++++++++++++++++----- 2 files changed, 37 insertions(+), 9 deletions(-) diff --git a/MC/bin/o2dpg_sim_workflow.py b/MC/bin/o2dpg_sim_workflow.py index 99c50b19a..ccd53722b 100755 --- a/MC/bin/o2dpg_sim_workflow.py +++ b/MC/bin/o2dpg_sim_workflow.py @@ -26,6 +26,8 @@ parser.add_argument('-confKey',help='generator or trigger configuration key values, for example: GeneratorPythia8.config=pythia8.cfg', default='') parser.add_argument('-eCMS',help='CMS energy', default=-1) +parser.add_argument('-eBeamA',help='Beam A energy', default=6499.) #6369 PbPb, 2.510 pp 5 TeV, 4 pPb +parser.add_argument('-eBeamB',help='Beam B energy', default=-1) parser.add_argument('-col',help='collision sytem: pp, PbPb, pPb, Pbp, ...', default='pp') parser.add_argument('-ptHatBin',help='pT hard bin number', default=-1) parser.add_argument('-ptHatMin',help='pT hard minimum when no bin requested', default=0) @@ -160,7 +162,9 @@ def getDPL_global_options(bigshm=False,nosmallrate=False): # function encapsulating the signal sim part # first argument is timeframe id RNDSEED=args.seed # 0 means random seed ! - ECMS=args.eCMS + ECMS=float(args.eCMS) + EBEAMA=float(args.eBeamA) + EBEAMB=float(args.eBeamB) NSIGEVENTS=args.ns GENERATOR=args.gen INIFILE='' @@ -231,7 +235,7 @@ def getDPL_global_options(bigshm=False,nosmallrate=False): PDGA=1000822080 # Pb PDGB=1000822080 # Pb if ECMS < 0: # assign 5.02 TeV to Pb-Pb - print('>>> Set CM Energy to PbPb case 5.02 TeV') + print('o2dpg_sim_workflow: Set CM Energy to PbPb case 5.02 TeV') ECMS=5020.0 if COLTYPE == 'pPb': @@ -242,8 +246,19 @@ def getDPL_global_options(bigshm=False,nosmallrate=False): PDGA=1000822080 # Pb PDGB=2212 # proton - if ECMS < 0: - print('Error: Collision Energy not set!!!') + # If not set previously, set beam energy B equal to A + if EBEAMB < 0 and ECMS < 0: + EBEAMB=EBEAMA + print('o2dpg_sim_workflow: Set beam energy same in A and B beams') + if COLTYPE=="pPb" or COLTYPE=="Pbp": + print('o2dpg_sim_workflow: Careful! both beam energies are the same') + + if ECMS > 0: + if COLTYPE=="pPb" or COLTYPE=="Pbp": + print('o2dpg_sim_workflow: Careful! ECM set for pPb/Pbp collisions!') + + if ECMS < 0 and EBEAMA < 0 and EBEAMB < 0: + print('o2dpg_sim_workflow: Error! CM or Beam Energy not set!!!') exit(1) # produce the signal configuration diff --git a/MC/config/common/pythia8/utils/mkpy8cfg.py b/MC/config/common/pythia8/utils/mkpy8cfg.py index 49e2a5278..b02fbea97 100755 --- a/MC/config/common/pythia8/utils/mkpy8cfg.py +++ b/MC/config/common/pythia8/utils/mkpy8cfg.py @@ -12,13 +12,19 @@ help='The random seed') parser.add_argument('--idA', type=int, default='2212', - help='PDG code of projectile beam') + help='PDG code of projectile beam A') parser.add_argument('--idB', type=int, default='2212', - help='PDG code of target beam') + help='PDG code of target beam B') -parser.add_argument('--eCM', type=float, default='13000.', - help='Centre-of-mass energy') +parser.add_argument('--eBeamA', type=float, default='6499.', + help='Energy of beam A') + +parser.add_argument('--eBeamB', type=float, default='6499.', + help='Energy of beam B') + +parser.add_argument('--eCM', type=float, default='-1', + help='Centre-of-mass energy (careful!, better use beam energy)') parser.add_argument('--process', default='inel', choices=['none', 'inel', 'ccbar', 'bbbar', 'heavy', 'jets', 'dirgamma'], help='Process to switch on') @@ -76,7 +82,14 @@ fout.write('### beams \n') fout.write('Beams:idA = %d \n' % (args.idA)) fout.write('Beams:idB = %d \n' % (args.idB)) -fout.write('Beams:eCM = %f \n' % (args.eCM)) +if args.eCM > 0: + fout.write('Beams:eCM = %f \n' % (args.eCM)) +elif args.eBeamA > 0 and args.eBeamB > 0: + fout.write('Beams:eA = %f \n' % (args.eBeamA)) + fout.write('Beams:eB = %f \n' % (args.eBeamB)) +else: + print('mkpy8cfg.py: Error, CM or Beam Energy not set!!!') + exit(1) fout.write('\n') ### processes From f3126d94c35b022d13fc11e237b26d4e09ac624d Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Balbastre Date: Mon, 29 Mar 2021 15:02:18 +0200 Subject: [PATCH 7/8] use same parameter naming as PYTHIA for collision and beams energy --- MC/bin/o2dpg_sim_workflow.py | 12 ++++++------ MC/config/common/pythia8/utils/mkpy8cfg.py | 10 +++++----- 2 files changed, 11 insertions(+), 11 deletions(-) diff --git a/MC/bin/o2dpg_sim_workflow.py b/MC/bin/o2dpg_sim_workflow.py index ccd53722b..472be7d54 100755 --- a/MC/bin/o2dpg_sim_workflow.py +++ b/MC/bin/o2dpg_sim_workflow.py @@ -25,9 +25,9 @@ parser.add_argument('-ini',help='generator init parameters file, for example: ${O2DPG_ROOT}/MC/config/PWGHF/ini/GeneratorHF.ini', default='') parser.add_argument('-confKey',help='generator or trigger configuration key values, for example: GeneratorPythia8.config=pythia8.cfg', default='') -parser.add_argument('-eCMS',help='CMS energy', default=-1) -parser.add_argument('-eBeamA',help='Beam A energy', default=6499.) #6369 PbPb, 2.510 pp 5 TeV, 4 pPb -parser.add_argument('-eBeamB',help='Beam B energy', default=-1) +parser.add_argument('-eCM',help='CMS energy', default=-1) +parser.add_argument('-eA',help='Beam A energy', default=6499.) #6369 PbPb, 2.510 pp 5 TeV, 4 pPb +parser.add_argument('-eB',help='Beam B energy', default=-1) parser.add_argument('-col',help='collision sytem: pp, PbPb, pPb, Pbp, ...', default='pp') parser.add_argument('-ptHatBin',help='pT hard bin number', default=-1) parser.add_argument('-ptHatMin',help='pT hard minimum when no bin requested', default=0) @@ -162,9 +162,9 @@ def getDPL_global_options(bigshm=False,nosmallrate=False): # function encapsulating the signal sim part # first argument is timeframe id RNDSEED=args.seed # 0 means random seed ! - ECMS=float(args.eCMS) - EBEAMA=float(args.eBeamA) - EBEAMB=float(args.eBeamB) + ECMS=float(args.eCM) + EBEAMA=float(args.eA) + EBEAMB=float(args.eB) NSIGEVENTS=args.ns GENERATOR=args.gen INIFILE='' diff --git a/MC/config/common/pythia8/utils/mkpy8cfg.py b/MC/config/common/pythia8/utils/mkpy8cfg.py index b02fbea97..ca9a42875 100755 --- a/MC/config/common/pythia8/utils/mkpy8cfg.py +++ b/MC/config/common/pythia8/utils/mkpy8cfg.py @@ -17,10 +17,10 @@ parser.add_argument('--idB', type=int, default='2212', help='PDG code of target beam B') -parser.add_argument('--eBeamA', type=float, default='6499.', +parser.add_argument('--eA', type=float, default='6499.', help='Energy of beam A') -parser.add_argument('--eBeamB', type=float, default='6499.', +parser.add_argument('--eB', type=float, default='6499.', help='Energy of beam B') parser.add_argument('--eCM', type=float, default='-1', @@ -84,9 +84,9 @@ fout.write('Beams:idB = %d \n' % (args.idB)) if args.eCM > 0: fout.write('Beams:eCM = %f \n' % (args.eCM)) -elif args.eBeamA > 0 and args.eBeamB > 0: - fout.write('Beams:eA = %f \n' % (args.eBeamA)) - fout.write('Beams:eB = %f \n' % (args.eBeamB)) +elif args.eA > 0 and args.eB > 0: + fout.write('Beams:eA = %f \n' % (args.eA)) + fout.write('Beams:eB = %f \n' % (args.eB)) else: print('mkpy8cfg.py: Error, CM or Beam Energy not set!!!') exit(1) From 85894c9d2a18778c52e8bd409949ca483c36a73f Mon Sep 17 00:00:00 2001 From: Gustavo Conesa Balbastre Date: Mon, 29 Mar 2021 15:20:23 +0200 Subject: [PATCH 8/8] update example comments --- MC/bin/o2dpg_sim_workflow.py | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) diff --git a/MC/bin/o2dpg_sim_workflow.py b/MC/bin/o2dpg_sim_workflow.py index 472be7d54..d7ecd3cd6 100755 --- a/MC/bin/o2dpg_sim_workflow.py +++ b/MC/bin/o2dpg_sim_workflow.py @@ -5,10 +5,16 @@ # It aims to handle the different MC possible configurations # It just creates a workflow.json txt file, to execute the workflow one must execute right after # ${O2DPG_ROOT}/MC/bin/o2_dpg_workflow_runner.py -f workflow.json +# # Execution examples: -# ./o2dpg_sim_workflow.py -e TGeant3 -nb 0 -ns 2 -j 8 -tf 1 -mod "-m TPC" -proc "jets" -ptTrigMin 3.5 -ptHatBin 3 -trigger "external" -ini "\$O2DPG_ROOT/MC/config/PWGGAJE/ini/trigger_decay_gamma.ini" --embedding False +# - pp PYTHIA jets, 2 events, triggered on high pT decay photons on EMCal acceptance, eCMS 13 TeV +# ./o2dpg_sim_workflow.py -e TGeant3 -ns 2 -j 8 -tf 1 -mod "--skipModules ZDC" -col pp -eCM 13000 \ +# -proc "jets" -ptTrigMin 3.5 -acceptance 4 -ptHatBin 3 \ +# -trigger "external" -ini "\$O2DPG_ROOT/MC/config/PWGGAJE/ini/trigger_decay_gamma.ini" # -# ./o2dpg_sim_workflow.py -e TGeant3 -nb 0 -ns 2 -j 8 -tf 1 -mod "--skipModules ZDC" -proc "ccbar" --embedding True +# - pp PYTHIA ccbar events embedded into heavy-ion environment, 2 PYTHIA events into 1 bkg event, beams energy 2.510 +# ./o2dpg_sim_workflow.py -e TGeant3 -nb 1 -ns 2 -j 8 -tf 1 -mod "--skipModules ZDC" \ +# -col pp -eA 2.510 -proc "ccbar" --embedding # import argparse