Skip to content
Snippets Groups Projects
Commit 38b4a44d authored by Klaus Rabbertz's avatar Klaus Rabbertz
Browse files
parents 40948d27 a1e965f0
Branches
No related tags found
No related merge requests found
......@@ -11,11 +11,9 @@ job_file_dir_cleanup: False
retries = 5
gfal_transfer_timeout = 300
gfal_transfer_nbstreams = 1
# GridKa examples, oldest to newest protocol; try root or davs first
#base = gsiftp://cmssrm-kit.gridka.de:2811//pnfs/gridka.de/cms/disk-only/store/user/USERNAME/YOURPATH
#base = srm://cmssrm-kit.gridka.de:8443/srm/managerv2?SFN=/pnfs/gridka.de/cms/disk-only/store/user/USERNAME/YOURPATH
#base = root://cmsxrootd-kit.gridka.de//store/user/USERNAME/YOURPATH
#base = davs://cmswebdav-kit.gridka.de:2880/pnfs/gridka.de/cms/disk-only/store/user/USERNAME/YOURPATH
# GridKa examples:
#base = davs://cmswebdav-kit-disk.gridka.de:2880/pnfs/gridka.de/cms/disk-only/store/user/CERNUSERNAME/YOURPATH
#base = root://cmsxrootd-kit-disk.gridka.de:1094//store/user/CERNUSERNAME/YOURPATH
# DESY
#base = srm://dcache-se-cms.desy.de:8443/srm/managerv2?SFN=/pnfs/desy.de/cms/tier2/store/user/USERNAME/YOURPATH
......
This diff is collapsed.
!###############################################!
!## ##!
!## _ ___ ____ ____ ____________ ##!
!## / |/ / |/ / / / __ \__ / / __/_ __/ ##!
!## / / / /__/ /_/ / // / _/ / / ##!
!## /_/|_/_/|_/____/\____/\___/___/ /_/ ##!
!## ##!
!## RUNCARD ##!
!###############################################!
PROCESS 2jetfc
collider = pp
sqrts = 8000
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@-CMS8-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 = 50.0
select jets_abs_y max = 5.0
select njets min = 2
! Leading jets
select abs_yj1 max = 3.0
select abs_yj2 max = 3.0
END_SELECTORS
HISTOGRAMS
ptavg_j12 > fnl3832-fc-v2_yb0_ys0_ptavgj12 [133,153,174,196,220,245,272,300,330,362,395,430,468,507,548,592,638,686,737,790,846,905,967,1032,1101,1172,1248,1327,1410,1497,1588,1784] grid=fnl3832-fc-v2_yb0_ys0_ptavgj12.fast
HISTOGRAM_SELECTORS
select yboost min=0.0 max=1.0
select ystar min=0.0 max=1.0
END_HISTOGRAM_SELECTORS
ptavg_j12 > fnl3832-fc-v2_yb0_ys1_ptavgj12 [133,153,174,196,220,245,272,300,330,362,395,430,468,507,548,592,638,686,737,790,846,905,967,1032,1101,1172,1248] grid=fnl3832-fc-v2_yb0_ys1_ptavgj12.fast
HISTOGRAM_SELECTORS
select yboost min=0.0 max=1.0
select ystar min=1.0 max=2.0
END_HISTOGRAM_SELECTORS
ptavg_j12 > fnl3832-fc-v2_yb0_ys2_ptavgj12 [133,153,174,196,220,245,272,300,330,362,395,430,468,507,548] grid=fnl3832-fc-v2_yb0_ys2_ptavgj12.fast
HISTOGRAM_SELECTORS
select yboost min=0.0 max=1.0
select ystar min=2.0 max=3.0
END_HISTOGRAM_SELECTORS
ptavg_j12 > fnl3832-fc-v2_yb1_ys0_ptavgj12 [133,153,174,196,220,245,272,300,330,362,395,430,468,507,548,592,638,686,737,790,846,905,967,1032] grid=fnl3832-fc-v2_yb1_ys0_ptavgj12.fast
HISTOGRAM_SELECTORS
select yboost min=1.0 max=2.0
select ystar min=0.0 max=1.0
END_HISTOGRAM_SELECTORS
ptavg_j12 > fnl3832-fc-v2_yb1_ys1_ptavgj12 [133,153,174,196,220,245,272,300,330,362,395,430,468,507,548,592,638,686] grid=fnl3832-fc-v2_yb1_ys1_ptavgj12.fast
HISTOGRAM_SELECTORS
select yboost min=1.0 max=2.0
select ystar min=1.0 max=2.0
END_HISTOGRAM_SELECTORS
ptavg_j12 > fnl3832-fc-v2_yb2_ys0_ptavgj12 [133,153,174,196,220,245,272,300,330,362,395,430] grid=fnl3832-fc-v2_yb2_ys0_ptavgj12.fast
HISTOGRAM_SELECTORS
select yboost min=2.0 max=3.0
select ystar min=0.0 max=1.0
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 ptj1**6
CHANNELS
@CHANNEL@ ! LO, R, V, RRa, RRb, RV, VV ! no region entry anymore! No coloured particle in FS (W,Z,H) --> RR;
END_CHANNELS
......@@ -50,10 +50,11 @@ ScDescript {
# (Same description for all perturbative contributions!)
#--------------------------------------------------------------------
CodeDescript {
"NNLOJET_rev5918, please cite: T. Gehrmann et al., PoS RADCOR2017, 074 (2018), arXiv:1801.06415."
"NNLOJET_rev6591, 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."
"APPLfast interface, please cite: D. Britzger et al., Eur. Phys. J. C82, 930 (2022), arXiv:2207.13735."
"APPLfast2 interface, to be published."
}
#--------------------------------------------------------------------
......
#----------------------
# 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 {
"d3sigma-dijet_dpTavg_dyb_dys_[pb_GeV]"
"CMS_Collaboration"
"Dijet_pTavg_yb0-2_ys0-2"
"anti-kT_R=0.7"
"CMS-PAPER-SMP-16-011, arXiv:1705.02628, Eur. Phys. J. C 77, 746 (2017)."
"RIVET_ID=,INSPIRE_RECORD=1598460"
"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."
}
# --------------------------------------------
# Modify the table contributions
# --------------------------------------------
#--------------------------------------------------------------------
# Modify code description
# (Same description for all perturbative contributions!)
#--------------------------------------------------------------------
CodeDescript {
"NNLOJET_rev6591, 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. C82, 930 (2022), arXiv:2207.13735."
"APPLfast2 interface, to be published."
}
#--------------------------------------------------------------------
# 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")
"ptavgj12_[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 fnl3832-fc-v2 # (def.=NNLOJET process name) No white space allowed here!
ScenarioDescription { # Replace generic one-line text by desired (multi-line) description
"d3sigma-dijet_dpTavg_dyb_dys_[pb_GeV]"
"CMS_Collaboration"
"Dijet_pTavg_yb0-2_ys0-2"
"anti-kT_R=0.7"
"CMS-PAPER-SMP-16-011, arXiv:1705.02628, Eur. Phys. J. C 77, 746 (2017)."
"INSPIRE_RECORD=1598460"
"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)
......@@ -36,7 +36,7 @@ action() {
export PATH="$gitpath/tools:$gitpath/tools/plotting:$PATH"
# NNLOJET + APPLfast (cvmfs)
applfastpath="/cvmfs/etp.kit.edu/fastnnlo-2888_bridge-0.0.46-WiP_202307_lcg102-c7-gcc11"
applfastpath="/cvmfs/etp.kit.edu/fastnnlo-3000_bridge-1.0.0_202309_lcg102-c7-gcc11"
#applfastpath="/eos/theory/project/applfast/local"
# LCG not needed here since source'd already above!
# NNLOJET & APPLfast
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment