Skip to content
Snippets Groups Projects
Commit e5729bd4 authored by Klaus Rabbertz's avatar Klaus Rabbertz
Browse files

Add setup of 2jetfc.fnl2412eff-fc-v2 production

parent 0c195248
Branches
No related tags found
No related merge requests found
[core]
no_lock = True
# Set local scheduler
local_scheduler = True
#default-scheduler-host = condorcentral.etp.kit.edu
#default-scheduler-port = 8082
[scheduler]
retry_count = 0
retry_delay = 86400
[worker]
keep_alive = False
ping_interval = 20
wait_interval = 20
max_reschedules = 1
[DEFAULT]
# Switch my debugging on/off (not yet implemented in most tasks)
#my_debug = True
# Name of your analysis
name = fnl2412eff-fc-v2
# If known, set proper Rivet_Id for data comparison plots
rivetid = CMS_2013_I1208923
# NNLOJET process and job name (the jobname is $channel-$jobnameext), and technical cutoff with leading 'y'
process = 2jetfc
jobnameext = CMS7-ak07
cutoff = y1.00E-08
multi_channel = .false. .false. .false. 5 5 5 5
# NNLOJET channels
# Complete LIST to be used in production; serves to create dictionary with settings per channel; settings must have same lengths!
allchannels = ["LO", "R", "V", "RRa", "RRb", "RV", "VV"]
# Default space-separated list of channels to be used in specific tasks; can be used as command line option.
# With coloured particles in final state (jets): 7 channels
channels = LO R V RRa RRb RV VV
# No coloured particle in final state (W,Z,H): 6 channels
# channels = LO R V RR RV VV
# Grid merging setup as defined in combine.ini; check compatibility
final_tables = {
"LO": ["LO"],
# Use the following lines when NLO is ready to be merged in as well (don't forget to add/remove a comma above)
"NLO": ["LO", "R", "V"],
"NLO_only": ["R", "V"],
# Use the following lines when NNLO is ready to be merged in as well (don't forget to add/remove a comma above)
"NNLO_only": ["RRa", "RRb", "RV", "VV"],
"NNLO": ["LO", "R", "V", "RRa", "RRb", "RV", "VV"]
}
# List of all observables (APPLfast grid names for NNLOJET histograms)
observables = [
"fnl2412eff-fc-v2_ym0_m12",
"fnl2412eff-fc-v2_ym1_m12",
"fnl2412eff-fc-v2_ym2_m12",
"fnl2412eff-fc-v2_ym3_m12",
"fnl2412eff-fc-v2_ym4_m12"
]
# Default htcondor job submission configuration (modifiable for each task)
htcondor_accounting_group = cms.jet
htcondor_requirements = (TARGET.ProvidesCPU==true)
htcondor_remote_job = True
#htcondor_user_proxy = /tmp/x509up_u12010
htcondor_user_proxy = /home/rabbertz/.globus/x509_proxy
# Time in seconds
htcondor_walltime = 84000
htcondor_request_cpus = 1
htcondor_use_cpus = 1 1 1 1 1 1 1
# For all cores in total
htcondor_request_memory = 4096
htcondor_universe = docker
htcondor_docker_image = mschnepf/slc7-condocker
# Disk usage on WN (def = 1000000)
htcondor_request_disk = 2000000
# Create log files in htcondor jobs
transfer_logs = True
# Set tolerance for workflow success with failed branches
tolerance = 0
# Submit only missing htcondor workflow branches (should always be true)
only_missing = True
# Do not resubmit automatically (default is 5); important for debugging!
#retries = 3
retries = 2
# Bootstrap file to be sourced at beginning of htcondor jobs
bootstrap_file = bootstrap_NNLOJET_modules2_rev6591_bridge-1.0.0_202309.sh
# Local directories, defaults are 'grids', 'plots', 'tags', 'warmups'
# - unpacking and merging of grids in local directories $merge_dir/$name/[$channel|Combined]
#merge_dir = grids
# - production of plots in local directories $plots_dir/$name
#plots_dir = plots
# - storing tags to mark tasks as finished in separate local directories $tags_dir/[Combine|FnloCppread]
#tags_dir = tags
# - unpacking of NNLOJET and fastNLO warmups into local directories $warm_dir/[warmup|rewarmup|fastwarm|mergefastwarm]/$name
#warm_dir = warmups
# fastNLO cppread options for grid/table comparison to NNLOJET original
# Confer the NNLOJET run file
# PDF choice for production
pdf = NNPDF31_nnlo_as_0118
# Central scale choice for grid closure, cf. NNLOJET run file & nnlobridge setup [kScale1, kScale2, scale12, scale21, kProd = scale1*scale2]
# Either one scale for all or one per observable
scale = kScale1
# Central scale choice(s) for flex-table evaluation, cf. fastNLO enumerator EScaleFunctionalForm: "0"=kScale1, "1"=kScale2, "2"=kQuadraticSum, ..., "9"=kProd, ...
# For now one number for both, mur and muf; fixed-scale tables must have 0 here
scaleindex = 0
# Central plus scale variations, cf. NNLOJET run file & nnlobridge setup
# None: --> 1; fixed-scale: --> 7; flex-scale low pT: --> -6; flex-scale high pT: --> -12
# Either one number for all or one per observable
scalecombs = 7 7 7 7 7
ascode = LHAPDF
norm = no
# Included in analysis.tar.gz
# Path to base runcard file
base_runcard = 2jetfc.fnl2412eff-fc-v2.run
# Local directory with all steering files
steer_dir = steeringfiles
#
# START of the ACTION
#
[Warmup]
# produced @ grid storage under Warmup
# Override some defaults for this task
# htcondor config
# If BWFORCLUSTER is too full
#htcondor_requirements = (TARGET.ProvidesCPU==true)
#htcondor_requirements = (TARGET.CloudSite=="blade")
#htcondor_request_cpus = 1
#htcondor_use_cpus = 1 1 1 1 1 1 1
#htcondor_requirements = (TARGET.CloudSite=="topas")
#htcondor_request_cpus = 42
#htcondor_use_cpus = 28 28 28 28 28 28 28
htcondor_requirements = (TARGET.CloudSite=="BWFORCLUSTER")
#htcondor_request_cpus = 20
#htcondor_use_cpus = 20 20 20 20 20 20 20
# for all cores in total
#htcondor_request_memory = 8000
# Set multicore usage per channel
htcondor_request_cpus = 4
htcondor_use_cpus = 4 4 4 4 4 4 4
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# NNLOJET event count and integration steps for every channel
# MUST be of the same length as channels!
#starting_seed = 0
# 1st test setup, <~ 5 minutes
# (main integration may have run separately ...;
# see also luigi_warmup_nemo and luigi_warmup_topas, if existing, configs for settings)
# first seed for first branch, counting upwards
warmup_events = 10000 10000 10000 1000 1000 1000 1000
warmup_iterations = 3 3 3 3 3 3 3
# time in seconds
htcondor_walltime = 3600
# Make warmup unusable without explicit editing to avoid reproducing 'lost' warmups on-the-fly!
# Such jobs would either run for a long time or be of very low quality.
# With seeds starting below zero, NNLOJET reports an error, but unfortunately with return code=0.
# This is captured in the Warmup task.
starting_seed = -9
# Real warmup is from git@github.com:neorobin88/APPLfast2.git: ProcessRuncards/CMS7TeVJets1212.6660/mjjR07/modules2/FC/warmupgridwithColourSample
[CopyWarmup]
# untar'ed to local storage under $warm_dir/warmup/$name
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# Up to NLO only
#channels = LO R V
[PlotVegasGrids]
# produced @ local storage under $plots_dir
# gnuplotted into local storage $warm_dir/warmup/$name
# Needs gnuplot version > 4.6 --> does not work on Centos7 portal machines
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# Up to NLO only
#channels = LO R V
[Rewarmup]
# produced @ grid storage under Rewarmup
# Override some defaults for this task
# htcondor config
#htcondor_requirements = (TARGET.CloudSite=="blade")
htcondor_request_cpus = 1
htcondor_use_cpus = 1 1 1 1 1 1 1
#htcondor_requirements = (TARGET.CloudSite=="topas")
#htcondor_request_cpus = 42
#htcondor_use_cpus = 28 28 28 28 28 28 28
htcondor_requirements = (TARGET.CloudSite=="BWFORCLUSTER")
#htcondor_request_cpus = 20
#htcondor_use_cpus = 20 20 20 20 20 20 20
# for all cores in total
htcondor_request_memory = 8000
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# NNLOJET event count and integration steps for every channel
# MUST be of the same length as channels!
# Test setup
# first seed for first branch, counting upwards
starting_seed = 10
warmup_events = 1000 1000 1000 1000 1000 1000 1000
warmup_iterations = 3 3 3 3 3 3 3
# time in seconds
htcondor_walltime = 1800
# No rewarmup was run!
[CopyRewarmup]
# untar'ed to local storage under $warm_dir/rewarmup/$name
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
[FastWarm]
# produced @ grid storage under FastWarm
# override some defaults for this task
# accept as complete at 95%
acceptance = 0.95
# htcondor config
htcondor_request_memory = 2000
# If BWFORCLUSTER is too full
#htcondor_requirements = (TARGET.ProvidesCPU==true)
# If ETP blades are too slow
htcondor_requirements = (TARGET.CloudSite=="BWFORCLUSTER")
# Everything fast
#htcondor_requirements = ((TARGET.CloudSite=="BWFORCLUSTER")||(TARGET.CloudSite=="schnepf")||(TARGET.CloudSite=="topas"))
# Up to 6h
#htcondor_requirements = ((TARGET.CloudSite=="schnepf")||(TARGET.CloudSite=="topas")||(TARGET.CloudSite=="forHLR2"))
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# NNLOJET event count and number of jobs for each channel
starting_seeds = 1000 2000 3000 4000 5000 6000 7000
# 1st test setup, <~ 5 minutes
#fastwarm_events = 1000 200 1000 100 250 200 1000
#fastwarm_events = 1000000 200000 1000000 25000 250000 200000 1000000
#fastwarm_jobs = 1 1 1 1 1 1 1
#fastwarm_jobs = 2 2 2 2 2 2 2
# time in seconds
#htcondor_walltime = 180
# 2nd test setup, <~ 20-30 minutes (18, 20, 17, 17, 19, 20, 17 min)
#fastwarm_events = 20000000 4000000 20000000 500000 5000000 4000000 20000000
#fastwarm_jobs = 2 2 2 2 2 2 2
# time in seconds
#htcondor_walltime = 1800
# Final setup, ~22 hours (15, 21, 15, 18, 19, 22, 15)
fastwarm_events = 1000000000 250000000 1000000000 30000000 300000000 250000000 1000000000
fastwarm_jobs = 10 10 10 10 10 10 10
# time in seconds
htcondor_walltime = 86400
[CopyFastWarm]
# untar'ed to local storage under $warm_dir/fastwarm/$name
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
[MergeFastWarm]
# produced @ grid storage under MergeFastWarm
# Reuse merged fast warmup from fnl2412e for fix- and flex-scale tables
[CopyMergeFastWarm]
# untar'ed to local storage under $warm_dir/mergefastwarm/$name
[FastProd]
# produced @ grid storage under FastProd
# Uncomment to define subdirectory per allchannel to delimit file numbers; default '.'
# fastprod_subdirs = LO R V ...
# override some defaults for this task
# accept as complete at 95%
acceptance = 0.95
# If BWFORCLUSTER is too full
#htcondor_requirements = (TARGET.ProvidesCPU==true)
#htcondor_requirements = (TARGET.Cloudsite=="blade")
#htcondor_requirements = (TARGET.CloudSite=="topas")
# If ETP blades are too slow
#htcondor_requirements = (TARGET.CloudSite=="BWFORCLUSTER")
# Everything fast
htcondor_requirements = ((TARGET.CloudSite=="BWFORCLUSTER")||(TARGET.CloudSite=="topas"))
#htcondor_requirements = ((TARGET.CloudSite=="BWFORCLUSTER")||(TARGET.CloudSite=="schnepf")||(TARGET.CloudSite=="topas"))
#htcondor_request_memory = 7500
htcondor_request_memory = 4000
# Set multicore usage per channel
#htcondor_request_cpus = 1
#htcondor_use_cpus = 1 1 1 1 1 1 1
htcondor_request_cpus = 4
htcondor_use_cpus = 4 4 4 4 4 4 4
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
#fastprod_subdirs = LO R . . . . .
# NNLOJET seeds, event count, and number of jobs for each channel
starting_seeds = 10000 20000 30000 40000 50000 60000 70000
# 1st test setup, 1 core, ~20 min (37, 9.3h, 11h, 4.0h, 4.0h, 2.0h, 45 min), allcpu
# 1st test setup, 4 cores, ~20 min (18, 4.0h, 4.5h, 1.0h, 1.5h, 30, 15 min), topas
#fastprod_events = 20000000 1760000 6500000 16000 170000 24000 650000
#fastprod_jobs = 2 2 2 2 2 2 2
# time in seconds
#htcondor_walltime = 1800
# 2nd test setup, ~3-5 h ( h)
#fastprod_events =
#fastprod_jobs = 3 3 3 3 3 3 3
# time in seconds
#htcondor_walltime = 40000
# Final setup, ~22 h, 1,1,1,4,4,4,4 cores (22, xx, 25, 28, 23, xx, 18) h
#fastprod_events = 700000000 4000000 13000000 350000 2500000 1000000 55000000
# Final setup, ~22 h, 1,1,1,4,4,4,4 cores (--, 11, 18, --, --, 9, --) h
#fastprod_events = 700000000 2000000 11000000 280000 2300000 500000 63000000
# Final setup, ~22 h, 1,1,1,4,4,4,4 cores (24, 28, 24, 22, 18, 22, 20) h
fastprod_events = 700000000 4000000 12000000 280000 2300000 1000000 63000000
#fastprod_jobs = 10 10 10 0 0 0 0
#fastprod_jobs = 0 0 0 10 10 10 10
#fastprod_jobs = 25 100 50 10 10 10 10
#fastprod_jobs = 10 10 10 250 100 250 100
fastprod_jobs = 25 100 50 500 100 500 100
htcondor_walltime = 86000
[CopyTables]
# copied to local storage under $merge_dir/$name/CHANNEL, see above
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# Up to LO/NLO only
#channels = LO
#channels = LO R V
# accept as complete at 95%
acceptance = 0.95
# Uncomment the following line to have this task ignore the level
# of completeness of the previous required task (FastProd).
# In that case this task works on whatever is ready to be treated,
# but interrupts the chain of requirements in the workflow.
force = True
[FnloCppread]
# produced @ local storage under $merge_dir/$name/CHANNEL
# ATTENTION: Combine now requires FnloCppread instead of CopyTables such that all grids have been checked for readability beforehand!
# ATTENTION: Incomplete table availability with unlimited branches and mutiple workers leads to unending task
# --> Needs to be stopped/killed manually
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# Up to LO/NLO only
#channels = LO
#channels = LO R V
# Accept as complete at 95%
acceptance = 0.95
# Uncomment the following line to have this task ignore the level
# of completeness of the previous required task (CopyTables).
# In that case this task works on whatever is ready to be treated,
# but interrupts the chain of requirements in the workflow.
force = True
[Combine]
# produced @ local storage under $merge_dir/$name/Combined
# Always consumes only one worker, but uses multiple threads inside the script by NNLOJET as set below by number of cores
# Always needs all channels up to some order LO, NLO, or NNLO
channels = LO R V RRa RRb RV VV
#channels = LO R V
#channels = LO
# Path to combine.ini config; check compatibility to desired order
combine_ini = combine.ini
#combine_ini = combine-LO.ini
#combine_ini = combine-NLO.ini
# Update of only part of the contributions
#combine_ini = combine-Real.ini
#combine_ini = combine-RRaRV.ini
# Number of cores for NNLOJET combine script
cores = 20
[Recombine]
# produced @ local storage under $merge_dir/$name/Combined
# Produces the files for adding stat. uncertainties to the grids; profits from multiple workers
# Needs all input channels up to some order plus combined channels as defined in combine_ini
channels = LO R V RRa RRb RV VV NLO NLO_only NNLO NNLO_only
#channels = LO R V
#channels = LO
# Path to combine.ini config; check compatibility to desired order
combine_ini = combine.ini
#combine_ini = combine-LO.ini
#combine_ini = combine-NLO.ini
# Update of only part of the contributions
#combine_ini = combine-Real.ini
#combine_ini = combine-RRaRV.ini
# Number of cores for NNLOJET combine script
cores = 20
[MergeFastProd]
# produced @ local storage under $merge_dir/$name/Combined/Final
# Always needs all channels up to some order LO, NLO, or NNLO
#channels = LO R V RRa RRb RV VV
# Up to LO/NLO only
#channels = LO
#channels = LO R V
# NNLOJET weight file to use for merging (LO, NLO, or NNLO)
weightorder = NNLO
# Accept as complete at 95%
acceptance = 0.95
# Execute workflow as local workflow instead of htcondor workflow (useful for merging small amount of grids, to be removed later)
workflow = local
[MergeFinal]
# produced @ local storage under $merge_dir/$name/Combined/Final
# Channels are set via final_tables dictionary
final_tables = {
"LO": ["LO"],
# Use the following lines when NLO is ready to be merged in as well (don't forget to add/remove a comma above)
"NLO": ["LO", "R", "V"],
"NLO_only": ["R", "V"],
# Use the following lines when NNLO is ready to be merged in as well (don't forget to add/remove a comma above)
"NNLO_only": ["RRa", "RRb", "RV", "VV"],
"NNLO": ["LO", "R", "V", "RRa", "RRb", "RV", "VV"]
}
[AddStatunc]
# produced @ local storage under $merge_dir/$name/Combined/Final
# Requires Recombine and MergeFinal to have run
# ATTENTION: Requires fastNLO v2.5!
# Should be added for each order
orders = [ "LO", "NLO", "NNLO" ]
#orders = [ "LO", "NLO" ]
#orders = [ "LO" ]
# Rivet_Id must have been defined in default section.
# Make sure SteerModify.str contains a scenario description with RIVET_ID=
# Then provide here the histogram identifiers for data comparison plots if possible
histoids = ["d02-x01-y01", "d02-x01-y02", "d02-x01-y03", "d02-x01-y04", "d02-x01-y05", "d02-x01-y01", "d02-x01-y02", "d02-x01-y03", "d02-x01-y04", "d02-x01-y05" ]
[FnloCppreadFinal]
# produced @ local storage under $merge_dir/$name/Combined/Final
# Needs all input channels up to some order plus combined channels as defined in combine_ini
channels = LO R V RRa RRb RV VV NLO NLO_only NNLO NNLO_only
#channels = LO R V NLO NLO_only
#channels = LO
#######################
# New plotting setup! #
#######################
[PlotRuntime]
# produced @ local storage under $plots_dir
# Only requires CopyTables to have run
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# Up to LO/NLO only
#channels = LO
#channels = LO R V
# Uncomment the following line to have this task ignore the level
# of completeness of the previous required task (CopyTables).
# In that case this task works on whatever is ready to be treated,
# but interrupts the chain of requirements in the workflow.
force = True
[SingleGridClosure]
# produced @ local storage under $plots_dir
# Requires FnloCppread to have run
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# Up to LO/NLO only
#channels = LO
#channels = LO R V
# Uncomment the following line to have this task ignore the level
# of completeness of the previous required task (FnloCppread).
# In that case this task works on whatever is ready to be treated,
# but interrupts the chain of requirements in the workflow.
force = True
[MultiGridClosure]
# produced @ local storage under $plots_dir
# Requires FnloCppread to have run
# Uncomment the following line to have only channels changing with jet size R to complement other analysis
#channels = R RRa RRb RV
# Up to LO/NLO only
#channels = LO
#channels = LO R V
# Needs correct fastNLO 'scalecombs' parameter to be set, see above
# Uncomment the following line to have this task ignore the level
# of completeness of the previous required task (FnloCppread).
# In that case this task works on whatever is ready to be treated,
# but interrupts the chain of requirements in the workflow.
force = True
[MergedGridClosure]
# produced @ local storage under $plots_dir
# Requires FnloCppreadFinal to have run
# Needs all input channels up to some order plus combined channels as defined in combine_ini
channels = LO R V RRa RRb RV VV NLO NLO_only NNLO NNLO_only
# Up to LO/NLO only
#channels = LO
#channels = LO R V
# Only subset available
#channels = NNLO
[StatUncertainty]
# produced @ local storage under $plots_dir
# Highest available order (if NNLO too imprecise, maybe just up to NLO as well)
orders = [ "NNLO" ]
# Only needs one grid to read scenario info; numbers are from NNLOJET dat files
# Uncomment to show statistical uncertainties only for summed up channels per order (LO, NLO_only, NNLO_only)
#channellist = sum
[Absolute]
# produced @ local storage under $plots_dir
# Needs all input channels up to some order plus combined channels as defined in combine_ini
channels = LO R V RRa RRb RV VV NLO NLO_only NNLO NNLO_only
#channels = LO R V NLO NLO_only
#channels = LO
# Only subset available
#channels = NNLO
[KfacComparison]
# produced @ local storage under $plots_dir
# Highest available order
orders = ["NNLO"]
[ScaleUncertainty]
# produced @ local storage under $plots_dir
# Highest available order (if NNLO too imprecise, maybe just up to NLO as well)
orders = [ "NNLO" ]
# Needs correct fastNLO 'scaleindex' parameter to be set
# Only one scale choice in this production
scaleindex = 0
# Show statistical uncertainties from corresponding NNLOJET dat files directly, "", or re-evaluated from fastNLO grid, "0"
# Note: '0' needs '_v25' tables (output of task AddStatunc)
#datfiles = [ "", "0" ]
datfiles = [ "" ]
[PdfUncertainty]
# produced @ local storage under $plots_dir
# Might be interesting for each order (Attention: ME and PDF orders are not matched!)
orders = [ "NNLO" ]
pdfsets = CT14nnlo,MMHT2014nnlo68cl,NNPDF31_nnlo_as_0118,ABMP16als118_5_nnlo
# Needs correct fastNLO 'scaleindex' parameter to be set
# One scale choice fully sufficient here
scaleindex = 0
[DataComparison]
# produced @ local storage under $plots_dir
# Requires AddStatunc to have run
# Produces Rivet-like comparison plots to data, if available.
# Rivet_Id must have been defined in default section.
# ATTENTION: Requires fastNLO + YODA + Rivet setup with Python2 because of malfunction in Rivet with Python3
# Show comparison to all orders
orders = [ "LO", "NLO", "NNLO" ]
#orders = [ "LO", "NLO" ]
#orders = [ "LO" ]
# Needs correct fastNLO 'scalename' parameter to be set
scalenames = [ "kScale1" ]
# Select theory uncertainty to show, default is 6P scale uncertainty, alternatives are:
# NN, 2P, 6P, HS, HA, HP, HC, MC, L6, AS, ST (see `fnlo-tk-yodaout -h` for details)
#uncertainty = ST
# Rescale uncertainties to be shown by factor, default=1
#dxsrescale = 10
!###############################################!
!## ##!
!## _ ___ ____ ____ ____________ ##!
!## / |/ / |/ / / / __ \__ / / __/_ __/ ##!
!## / / / /__/ /_/ / // / _/ / / ##!
!## /_/|_/_/|_/____/\____/\___/___/ /_/ ##!
!## ##!
!## RUNCARD ##!
!###############################################!
PROCESS 2jetfc
collider = pp
sqrts = 7000
jet = antikt[0.7]
jet_exclusive = .false. ! inclusive or exclusive Njet
jet_recomb = V4 ! recombination scheme: V4 is 4-vector addition
decay_type = 0 ! 0 no decay?
END_PROCESS
RUN @CHANNEL@-CMS7-ak07 ! job name
PDF = NNPDF31_nnlo_as_0118[0] ! PDFset[member]
tcut = 1d-8 ! y0, appears in filenames of warmup
iseed = @SEED@ ! seed index
cache_kinematics = .true. ! caching for kinematics; normally true
reset_vegas_grid = .true. ! normally true
lips_reduce = .true. ! optimisation for identical particles in FS; normally true
angular_average = .true. ! azimuthal rotation; only relevant for NNLO; can always be true; might be false for RR
multi_channel = @MULTI_CHANNEL@ ! color sampling; true faster, less precise; true/false must be consistent between warmup & production; can be set to 5 instead of true
@RUNMODE@ = @EVENTS@[@ITERATIONS@] ! Normally #iterations is 1 for production
END_RUN
PARAMETERS
MASS[H] = 125
WIDTH[H] = 0.004029643852
MASS[Z] = 91.1876
WIDTH[Z] = 2.4952
MASS[W] = 80.398
WIDTH[W] = 2.1054
MASS[t] = 173.2
WIDTH[t] = 1.41
MASS[b] = 4.18
MASS[c] = 1.275
MASS[tau] = 1.777
END_PARAMETERS
@UNIT_PHASE@
SELECTORS
! Base event selection
select jets_pt min = 30.0
select jets_abs_y max = 3.0
select njets min = 2
! Leading jets
select ptj1 min = 60.0
! select ptj2 min = 30.0 ! redundant
select abs_yj1 max = 2.5
select abs_yj2 max = 2.5
select m12 min = 197.0 max = 5058.0
! select max_abs_y_j12 max = 2.5 ! redundant
END_SELECTORS
HISTOGRAMS
m12 > fnl2412eff-fc-v2_ym0_m12 [197.0,296.0,419.0,565.0,740.0,944.0,1181.0,1455.0,1770.0,2132.0,2546.0,3019.0,3416.0,4010.0] grid=fnl2412eff-fc-v2_ym0_m12.fast ! output_type = 1 (like modules 1, default 0 means no channel breakdown), normally 0 now for modules 2
HISTOGRAM_SELECTORS
select max_abs_y_j12 min=0.0 max=0.5
END_HISTOGRAM_SELECTORS
m12 > fnl2412eff-fc-v2_ym1_m12 [270.0,386.0,526.0,693.0,890.0,1118.0,1383.0,1687.0,2037.0,2438.0,2895.0,3416.0,4010.0] grid=fnl2412eff-fc-v2_ym1_m12.fast
HISTOGRAM_SELECTORS
select max_abs_y_j12 min=0.5 max=1.0
END_HISTOGRAM_SELECTORS
m12 > fnl2412eff-fc-v2_ym2_m12 [419.0,565.0,740.0,944.0,1181.0,1455.0,1770.0,2132.0,2546.0,3147.0,3854.0,4509.0] grid=fnl2412eff-fc-v2_ym2_m12.fast
HISTOGRAM_SELECTORS
select max_abs_y_j12 min=1.0 max=1.5
END_HISTOGRAM_SELECTORS
m12 > fnl2412eff-fc-v2_ym3_m12 [565.0,740.0,944.0,1181.0,1455.0,1770.0,2132.0,2546.0,3019.0,3558.0,5058.0] grid=fnl2412eff-fc-v2_ym3_m12.fast
HISTOGRAM_SELECTORS
select max_abs_y_j12 min=1.5 max=2.0
END_HISTOGRAM_SELECTORS
m12 > fnl2412eff-fc-v2_ym4_m12 [1000.0,1246.0,1530.0,1856.0,2231.0,2659.0,3147.0,3704.0,5058.0] grid=fnl2412eff-fc-v2_ym4_m12.fast
HISTOGRAM_SELECTORS
select max_abs_y_j12 min=2.0 max=2.5
END_HISTOGRAM_SELECTORS
END_HISTOGRAMS
SCALES
! Central scale followed by 6 scale-factor variations
! To correspond to fastNLO, order must be:
muf = 1.0 * m12 mur = 1.0 * m12
muf = 0.5 * m12 mur = 0.5 * m12
muf = 2.0 * m12 mur = 2.0 * m12
muf = 1.0 * m12 mur = 0.5 * m12
muf = 0.5 * m12 mur = 1.0 * m12
muf = 2.0 * m12 mur = 1.0 * m12
muf = 1.0 * m12 mur = 2.0 * m12
END_SCALES
REWEIGHT m12**4
CHANNELS
@CHANNEL@ ! LO, R, V, RRa, RRb, RV, VV ! no region entry anymore! No coloured particle in FS (W,Z,H) --> RR;
END_CHANNELS
#----------------------
# Input and output
#----------------------
#InTable TableToModify.tab.gz
#OutTable ModifiedTable.tab.gz
# or define via command line
#--------------------
# Info print out
#--------------------
Verbosity INFO # Set print out verbosity
PrintSteeringCard true # Print steering provided on command line or within this file
PrintInputA1 true # Print table header of input table
PrintInputA2 true # Print table scenario of input table
#PrintOutputA1 true # Print table header of output table
#PrintOutputA2 true # Print table scenario of output table
#-------------------------------------------
# Modify parameters of the table header
# Leave value empty for no modification
#-------------------------------------------
Itabversion 25000 # Change table version (Expert option!)
#ScenName fnlabcd # Change scenario name
# --------------------------------------------
# Modify parameters of the table scenario
# Leave value empty for no modification
# --------------------------------------------
AttachScDescription false # Append to (true) or replace (false) existing scenario description
ScDescript {
"d2sigma-dijet_dm12_dym_[pb_GeV]"
"CMS_Collaboration"
"Dijet_m12_ym0-4"
"anti-kT_R=0.7"
"CMS-PAPER-QCD-11-004, arXiv:1212.6660, Phys. Rev. D87, 112002 (2013)."
"RIVET_ID=,INSPIRE_RECORD=1208923"
"THEORY_LABEL=NNLOJET"
"provided by:"
"fastNLO_2.5.0"
"If you use this table, please cite:"
" D. Britzger, K. Rabbertz, F. Stober, M. Wobisch, Proc. DIS 2012, 217 (2012), arXiv:1208.3641."
}
# --------------------------------------------
# Modify the table contributions
# --------------------------------------------
#--------------------------------------------------------------------
# Modify code description
# (Same description for all perturbative contributions!)
#--------------------------------------------------------------------
CodeDescript {
"NNLOJET_rev5918, please cite: T. Gehrmann et al., PoS RADCOR2017, 074 (2018), arXiv:1801.06415."
"Dijet production, please cite: J. Currie et al., Phys. Rev. Lett. 119, 152001 (2017), arXiv:1705.10271."
"Dijet scale choice, please cite: A. Gehrmann-de Ridder et al., Phys. Rev. Lett. 123, 102001 (2019), arXiv:1905.09047."
"APPLfast interface, please cite: D. Britzger et al., Eur. Phys. J. C79, 845 (2019), arXiv:1906.05303."
}
#--------------------------------------------------------------------
# Add InfoBlocks with statistical/numerical uncertainty
#--------------------------------------------------------------------
InfoBlockFlag1 0 # Default: 0 Statistical/numerical uncertainty
InfoBlockFlag2 1 # Default: 1 = quadratic addition; else: 0 = linear addition
InfoBlockDescr { # Same description for each contribution
"Default statistical uncertainties from NNLOJET."
}
# For NNLOJET dat files also add:
#
#InfoBlockFiles { # For each order separately: NNLOJET dat files, fastNLO statunc log files, txt files using 1 or 2 columns
# *.LO.*.dat
# *.NLO_only.*.dat
# *.NNLO_only.*.dat
#}
#InfoBlockOrders { # Attribution of dat files to respective orders
# "LO"
# "NLO"
# "NNLO"
#}
# -*-sh-*-
# ==================================================================== #
#
# A steering file for creating a fastNLO table/grid (with NNLOJET)
# This is an example for a global steering for ALL tables/grids of a scenario
#
# The grid production is performed in three steps:
# I. A warmup run for NNLOJET to perform essential integrations
# See NNLOJET for details. fastNLO is not used or called at this stage.
# II. A warmup run for fastNLO that determines the accessed phase space
# in momentum fraction x and process scale mu_r, mu_f.
# Only a few fastNLO settings are used here, see below.
# The output is written into a warmup file named after the
# the NNLOJET process name, e.g. 'ZJ', the NNLOJET job name, e.g.
# 'LO-CMS13-ak07', and the NNLOJET grid name, e.g.
# ptz > ZJtriple_yb0_ystar0_ptz [30, 100, 200, 300, 400, 1000] grid = ZJtriple_yb0_ystar0_ptz.fast
# NNLOJET histo fastNLO grid
# ==> ZJ.LO-CMS13-ak07.ZJtriple_yb0_ystar0_ptz.wrm
# This warmup file stores some settings for fastNLO in addition to
# the accessed phase space (with number of events per bin!) for
# each order of NNLOJET run (LO, R, V, RRa, RRb, RV, VV).
# For a later combination ALL production runs in step III below
# MUST use the SAME warmup file. The latter should be the largest
# accessed phase space as combined from any order in NNLOJET!
# ==> the expected filename for this unique warmup file is
# ZJ.ZJtriple_yb0_ystar0_ptz.wrm
# without the job name info.
# III. Many production runs filling fastNLO grids for various orders
# of a process from NNLOJET that later must be combined.
# See NNLOJET for details.
# Most fastNLO settings are taken into account here.
#
# This steering file has to be used in combination with a NNLOJET
# run card, where most of the essential seetings like the
# observable and its binning are defined.
# Modifications or additional settings respect the following hierarchy:
# 1. From the interface package nnlo-bridge default settings for
# fastNLO with NNLOJET are made and information from NNLOJET is
# passed to fastNLO.
# 2. A steering file named after the NNLOJET grid name for fastNLO,
# e.g. from the example above
# ==> ZJ.ZJtriple_yb0_ystar0_ptz.str
# can be used to set steer parameters individually for each fastNLO grid.
# Settings from step 1 can be overwritten. Be careful not to make
# changes inconsistent with what is defined in the NNLOJET run card
# or what is mandatory for the interface to work!
# 3a. In production runs the previously derived warmup file is read.
# 3b. A global steering file named after only the first part of
# the NNLOJET grid name for fastNLO, 'ZJtriple', which is also used
# as fastNLO scenario name,
# ==> ZJ.ZJtriple.str
# can be used to set steer parameters globally for all fastNLO grids.
# Settings from step 2 are not overwritten, but lead to warnings.
# 4. In principal, parameters can also be set on the command line, but
# this has not been tested yet.
#
##### TODO: The special scale-factor settings are not necessary anymore with the modules2 interface
##### TODO: The method how to set two different scales in flex-scale grids needs to be tested
#
##### BEGIN OF DESCRIPTION NEEDING UPDATE
#
# Special remarks on scale settings:
# ----------------------------------
# Lacking sufficient scale choice information from NNLOJET the following
# settings are assumed for flexible-scale tables that allow two central
# scale choices to be stored simultaneously:
# DIS pp
# - the 1st scale set for muf defines scale1 scale2
# - the 1st scale set for mur defines scale2 scale1
# To check the closure between NNLOJET and fastNLO the proper
# ASYMMETRIC central scale choice must be used when evaluating a
# fastNLO grid with fnlo-tk-cppread, either scale12 or scale21.
#
# It is recommended to set ScaleDescriptionScale1 and
# ScaleDescriptionScale2 accordingly.
#
# If multiple entries per event are requested using NNLOJET COMPOSITE
# histos AND a scale definition varying with each entry is desired,
# then one central scale, e.g. ptmax=ptj1, must be defined for both
# muf and mur, and the InclusiveJets flag must be set in the fastNLO
# steering file.
#
# DIS pp
# - the 1st scale set for muf&mur defines scale1 scale2
# - the stored observable defines scale2 scale1
# To check the closure between NNLOJET and fastNLO the proper
# SYMMETRIC central scale choice must be used when evaluating a
# fastNLO grid with fnlo-tk-cppread, either scale1 or scale2.
#
# DIS(?) pp
# - the 1st scale set for muf&mur defines scale1 scale2=mu/obs(!)
# - the stored observable defines scale2 scale1
# To check the closure between NNLOJET and fastNLO the proper
# SYMMETRIC central scale choice must be used when evaluating a
# fastNLO grid with fnlo-tk-cppread, i.e. for pp scale1*scale2.
#
##### END OF DESCRIPTION NEEDING UPDATE
#
# ==================================================================== #
# ==================================================================== #
#
# Scenario specific settings
#
# Either set individually for each fastNLO grid or globally for all
#
# ==================================================================== #
# -------------------------------------------------------------------- #
# Settings to be modified already for Warmup
# (These are stored in the warmup file!)
# -------------------------------------------------------------------- #
#CheckScaleLimitsAgainstBins true # (def.=true) Set limits for scale nodes to bin borders, if possible. Good if scale equals the binned observable.
# Scales and scale factors must be set in the NNLOJET run card
ScaleDescriptionScale1 "m12_[GeV]" # (def.='scale1') Reset the 1st scale name and unit, e.g. "<pT_1,2>_[GeV]" (Note: The 1st scale must always be in units of [GeV]!)
ScaleDescriptionScale2 "m12_[GeV]" # (def.='scale2') Reset the 2nd scale name and unit (ONLY for flexible-scale tables)
#DifferentialDimension 1 # (must be 1) So far the interface to NNLOJET supports only 1-dim histograms. DO NOT CHANGE!
DimensionLabels { # Labels (symbol and unit) for the measurement dimension (from outer to inner "loop")
"m12_[GeV]" # The default following the example above would be: "ptz"
}
#DimensionIsDifferential { # (must be 2) Specify for each dimension whether
# 2 # 0 : the cross section is NOT differential, i.e. there are two bin borders,
#} # but NO division (normalization) by bin width
# 1 : the cross section is point-wise differential, i.e. only one point is given
# 2 : the cross section is bin-wise differential, i.e. there are two bin borders
# Since NNLOJET provides bin-wise differential distributions use option 2. DO NOT CHANGE!
# -------------------------------------------------------------------- #
# Description of scenario
#
# (Write each line of description in quotation marks)
# -------------------------------------------------------------------- #
ScenarioName fnl2412eff-fc-v2 # (def.=NNLOJET process name) No white space allowed here!
ScenarioDescription { # Replace generic one-line text by desired (multi-line) description
"d2sigma-dijet_dm12_dym_[pb_GeV]"
"CMS_Collaboration"
"Dijet_m12_ym0-ym4"
"anti-kT_R=0.7"
"CMS-PAPER-QCD-11-004, arXiv:1212.6660, Phys. Rev. D 87 (2013) 112002."
"RIVET_ID=CMS_2013_I1208923/D06-x01-y01,SPIRES_ID=9901922,INSPIRE_RECORD=1208923"
"THEORY_LABEL=NNLOJET modules2 FC"
"provided by:"
"fastNLO_2.5.1"
"If you use this table, please cite:"
" D. Britzger, K. Rabbertz, F. Stober, M. Wobisch, Proc. DIS 2012, 217 (2012), arXiv:1208.3641."
}
PublicationUnits 12 # (def.=12) Cross section unit to be used for output (negative power of 10, e.g. 12->pb, 15->fb)
# For comparison, e.g. via Rivet, should correspond to published data cross sections
# -------------------------------------------------------------------- #
# Cuts
# -------------------------------------------------------------------- #
# To bet set via NNLOJET run card
# -------------------------------------------------------------------- #
# Observable normalization
# -------------------------------------------------------------------- #
#CalculateBinSize true # (def.=true) Calculate bin width from lower and upper bin boundaries
#BinSizeFactor 1. # (def.=1.) Possibility to provide additional normalization factor, e.g. of 2. for bins in |y|
#BinSize { x1 x2 x3 ... } # (def. not set) If 'CalculateBinSize' is 'false' provide table with bin widths 'by hand' for normalization
# -------------------------------------------------------------------- #
# Observable binning
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# Use either 'SingleDifferentialBinning' or
# 'DoubleDifferentialBinning' or
# 'TripleDifferentialBinning'
# in accord with 'DifferentialDimension' above
# -------------------------------------------------------------------- #
# To bet set via the NNLOJET run card
# ==================================================================== #
#
# Process specific and therefore at least partially theory-code
# dependent settings
#
# ==================================================================== #
# To bet set directly from NNLOJET interface and via the NNLOJET run card
# ==================================================================== #
#
# A few special settings only used with NNLOJET so far
#
# ==================================================================== #
#InclusiveJets true # (def.=false) Flag to vary the mur, muf scale per entry, not just per event
# # Must be set in accordance with NNLOJET run card!
# # Only possible with flexible-scale grids
# ==================================================================== #
#
# General fastNLO settings
#
# ==================================================================== #
GlobalVerbosity INFO # (def.=WARNING) Global output verbosity of fastNLO toolkit,
# (and of user scenario code if implemented).
# Possible values are: DEBUG, MANUAL, INFO, WARNING, ERROR, SILENT
FlexibleScaleTable true # (def.=false) Create table fully flexible in mu_f (larger size), true, or
# table with fixed number of mu_f scale factors
CacheType 2 # Default caching settings
CacheMax 30
CacheCompare 10
# ==================================================================== #
#
# Specific fastNLO settings (normally do not need to be changed)
# Only change when you know what you are doing ...!
#
# ==================================================================== #
#OutputFilename fastNLO.tab # Overwrites default filename of fastNLO output table, e.g. ZJ.LO-CMS13.vBa.ZJtriple_yb0_ystar0_ptz.s92394.tab.gz
#OutputPrecision 8 # (def.=8) Number of decimal digits to store in output table
#OutputCompression true # (def.=true) Write out fastNLO table/grid in gzipped format (requires zlib)
#ScaleVariationFactors { # (def.=1.0 0.5 2.0) Must be set in accordance with NNLOJET run card! Factorization scale variations (only needed for fixed-scale tables)
# 1.0 0.5 2.0 # List of scale factors must include factor '1.0'
#} # Scale factors will be ordered according to fastNLO convention: (1, min, ... , max)
#ReadBinningFromSteering false # (must be true) Specify where the binning is defined.
#ApplyPDFReweighting true # (def.=true) Apply reweighting of PDFs for an optimized interpolation
# -------------------------------------------------------------------- #
# Choose fastNLO interpolation kernels and distance measures
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# Currently implemented interpolation kernels
# Catmull
# Lagrange
# OneNode
# Linear
#
# Currently implemented distance measures
# linear
# loglog025 eq. to (log(log(4*x)))
# log10
# sqrtlog10 eq. to sqrt(log_10(x))
#
# For a constant scale like M_Z use e.g.
# Mu1_Kernel OneNode
# Mu1_DistanceMeasure "log10"
# Mu1_NNodes 1
# -------------------------------------------------------------------- #
#X_Kernel Lagrange # (def.= pp: Lagrange; DIS: Catmull)
#X_DistanceMeasure sqrtlog10 # (def.= pp: sqrtlog10; DIS: log10)
#X_NNodes 20 # (def.= pp: 20; DIS: 18)
#X_NNodeCounting NodesPerBin # (def.=NodesPerBin)
#Mu1_Kernel Lagrange # (def.=Lagrange)
#Mu1_DistanceMeasure "loglog025" # (def.=loglog025)
#Mu1_NNodes 6 # (def.=6)
# Scale2 not used for fixed-scale tables
Mu2_Kernel OneNode # (def.=Lagrange)
Mu2_DistanceMeasure "loglog025" # (def.=loglog025)
Mu2_NNodes 1 # (def.=6)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment