WARNING: invalid set command param_card frblock 2 9.9e-01 #sinbma (failed to identify LHA information)

Asked by Sahithi

Hello experts,

I am running gridpacks in MADGRAPH5_aMC@NLO v2.6.5 for a 2HDM A -> ZH process using the model 2HDMtII_NLO and I see a couple of warnings "WARNING: invalid set command param_card frblock 1 1.00e+00 #tanbeta (failed to identify LHA information)" and "WARNING: invalid set command param_card frblock 2 9.9e-01 #sinbma (failed to identify LHA information)" when I set the tanbeta and sin(b-a) values in customize_cards.dat even though I see "frblock" in the model param card. Please let me know what I am missing. Thank you!

Customize_cards.dat:
set param_card frblock 1 1.00e+00 #tanbeta
set param_card frblock 2 9.9e-01 #sinbma
set param_card mass 25 125.0
set param_card mass 35 HMASS
set param_card mass 36 AMASS
set param_card mass 37 AMASS
set param_card width 35 5.303000e+01
set param_card width 36 2.377000e+00
set param_card width 37 1.746000e+01

proc_card.dat:
set group_subprocesses Auto
set ignore_six_quark_processes False
#set loop_color_flows False
set gauge unitary
set complex_mass_scheme False

import model 2HDMtII_NLO
define l+ = e+ mu+
define l- = e- mu-

generate p p > h3 > h2 z [QCD]
output AZHToLLtt -nojpeg

extramodels.dat:
#http://feynrules.irmp.ucl.ac.be/attachment/wiki/2HDM/2HDMtII_NLO.tar.gz
2HDMtII_NLO.tar.gz

madspin_card.dat:
set ms_dir ./madspingrid

#initialization parameters
#set Nevents_for_max_weigth 250 # number of events for the estimate of the max. weight
#set BW_cut 15 # cut on how far the particle can be off-shell
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
set max_running_process 1

define f = u c d s u~ c~ d~ s~ ve vm vt e- mu- ta- ve~ vm~ vt~ e+ mu+ ta+

# specify the decay for the final state particles
set spinmode none
decay h2 > t t~, ( t > w+ b, w+ > f f ), ( t~ > w- b~, w- > f f )
decay z > l+ l-
# running the actual code
launch
#

run_card.dat:

#*******************
# Running parameters
#*******************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Run to generate the grid pack *
#*********************************************************************
  .true. = gridpack !True = setting up the grid pack
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#*********************************************************************
  2500 = nevents ! Number of unweighted events requested
      0 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#*********************************************************************
        1 = lpp1 ! beam 1 type
        1 = lpp2 ! beam 2 type
     6500 = ebeam1 ! beam 1 total energy in GeV
     6500 = ebeam2 ! beam 2 total energy in GeV
#*********************************************************************
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#*********************************************************************
        0 = polbeam1 ! beam polarization for beam 1
        0 = polbeam2 ! beam polarization for beam 2
#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
 'lhapdf' = pdlabel ! PDF set
 $DEFAULT_PDF_SETS = lhaid
 $DEFAULT_PDF_MEMBERS = reweight_PDF ! if pdlabel=lhapdf, this is the lhapdf number
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 F = fixed_ren_scale ! if .true. use fixed ren scale
 F = fixed_fac_scale ! if .true. use fixed fac scale
 91.1880 = scale ! fixed ren scale
 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 1 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Matching - Warning! ickkw > 1 is still beta
#*********************************************************************
 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
 1 = highestmult ! for ickkw=2, highest mult group
 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
 1 = alpsfact ! scale factor for QCD emission vx
 F = chcluster ! cluster only according to channel diag
 F = pdfwgt ! for ickkw=1, perform pdf reweighting
 4 = asrwgtflavor ! highest quark flavor for a_s reweight
 T = clusinfo ! include clustering tag in output
 3.0 = lhe_version ! Change the way clustering information pass to shower.
#*********************************************************************
#**********************************************************
#
#**********************************************************
# Automatic ptj and mjj cuts if xqcut > 0
# (turn off for VBF and single top processes)
#**********************************************************
   T = auto_ptj_mjj ! Automatic setting of ptj and mjj
#**********************************************************
#
#**********************************
# BW cutoff (M+/-bwcutoff*Gamma)
#**********************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#**********************************************************
# Apply pt/E/eta/dr/mij cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#**********************************************************
   T = cut_decays ! Cut decay products
#*************************************************************
# Number of helicities to sum per event (0 = all helicities)
# 0 gives more stable result, but longer run time (needed for
# long decay chains e.g.).
# Use >=2 if most helicities contribute, e.g. pure QCD.
#*************************************************************
   0 = nhel ! Number of helicities used per event
#*******************
# Standard Cuts
#*******************
#
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 10 = ptj ! minimum pt for the jets
  0 = ptb ! minimum pt for the b
 10 = pta ! minimum pt for the photons
  0 = ptl ! minimum pt for the charged leptons
 0 = misset ! minimum missing Et (sum of neutrino's momenta)
  0 = ptheavy ! minimum pt for one heavy final state
 -1 = ptjmax ! maximum pt for the jets
 -1 = ptbmax ! maximum pt for the b
 -1 = ptamax ! maximum pt for the photons
 -1 = ptlmax ! maximum pt for the charged leptons
 -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
#*********************************************************************
# Minimum and maximum E's (in the center of mass frame) *
#*********************************************************************
  0 = ej ! minimum E for the jets
  0 = eb ! minimum E for the b
  0 = ea ! minimum E for the photons
  0 = el ! minimum E for the charged leptons
 -1 = ejmax ! maximum E for the jets
 -1 = ebmax ! maximum E for the b
 -1 = eamax ! maximum E for the photons
 -1 = elmax ! maximum E for the charged leptons
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
   5 = etaj ! max rap for the jets
  -1 = etab ! max rap for the b
  -1 = etaa ! max rap for the photons
  -1 = etal ! max rap for the charged leptons
   0 = etajmin ! min rap for the jets
   0 = etabmin ! min rap for the b
   0 = etaamin ! min rap for the photons
   0 = etalmin ! main rap for the charged leptons
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
 0.1 = drjj ! min distance between jets
 0 = drbb ! min distance between b's
 0 = drll ! min distance between leptons
 0 = draa ! min distance between gammas
 0 = drbj ! min distance between b and jet
 0.1 = draj ! min distance between gamma and jet
 0 = drjl ! min distance between jet and lepton
 0 = drab ! min distance between gamma and b
 0 = drbl ! min distance between b and lepton
 0.1 = dral ! min distance between gamma and lepton
 -1 = drjjmax ! max distance between jets
 -1 = drbbmax ! max distance between b's
 -1 = drllmax ! max distance between leptons
 -1 = draamax ! max distance between gammas
 -1 = drbjmax ! max distance between b and jet
 -1 = drajmax ! max distance between gamma and jet
 -1 = drjlmax ! max distance between jet and lepton
 -1 = drabmax ! max distance between gamma and b
 -1 = drblmax ! max distance between b and lepton
 -1 = dralmax ! maxdistance between gamma and lepton
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#*********************************************************************
 0 = mmjj ! min invariant mass of a jet pair
 0 = mmbb ! min invariant mass of a b pair
 0 = mmaa ! min invariant mass of gamma gamma pair
 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
 -1 = mmjjmax ! max invariant mass of a jet pair
 -1 = mmbbmax ! max invariant mass of a b pair
 -1 = mmaamax ! max invariant mass of gamma gamma pair
 -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#*********************************************************************
# Minimum and maximum invariant mass for all letpons *
#*********************************************************************
  0 = mmnl ! min invariant mass for all letpons (l+- and vl)
 -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons *
#*********************************************************************
 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
 -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0 = xptj ! minimum pt for at least one jet
 0 = xptb ! minimum pt for at least one b
 0 = xpta ! minimum pt for at least one photon
 0 = xptl ! minimum pt for at least one charged lepton
#*********************************************************************
# Control the pt's of the jets sorted by pt *
#*********************************************************************
 0 = ptj1min ! minimum pt for the leading jet in pt
 0 = ptj2min ! minimum pt for the second jet in pt
 0 = ptj3min ! minimum pt for the third jet in pt
 0 = ptj4min ! minimum pt for the fourth jet in pt
 -1 = ptj1max ! maximum pt for the leading jet in pt
 -1 = ptj2max ! maximum pt for the second jet in pt
 -1 = ptj3max ! maximum pt for the third jet in pt
 -1 = ptj4max ! maximum pt for the fourth jet in pt
 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#*********************************************************************
# Control the pt's of leptons sorted by pt *
#*********************************************************************
 0 = ptl1min ! minimum pt for the leading lepton in pt
 0 = ptl2min ! minimum pt for the second lepton in pt
 0 = ptl3min ! minimum pt for the third lepton in pt
 0 = ptl4min ! minimum pt for the fourth lepton in pt
 -1 = ptl1max ! maximum pt for the leading lepton in pt
 -1 = ptl2max ! maximum pt for the second lepton in pt
 -1 = ptl3max ! maximum pt for the third lepton in pt
 -1 = ptl4max ! maximum pt for the fourth lepton in pt
#*********************************************************************
# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
 0 = htjmin ! minimum jet HT=Sum(jet pt)
 -1 = htjmax ! maximum jet HT=Sum(jet pt)
 0 = ihtmin !inclusive Ht for all partons (including b)
 -1 = ihtmax !inclusive Ht for all partons (including b)
 0 = ht2min ! minimum Ht for the two leading jets
 0 = ht3min ! minimum Ht for the three leading jets
 0 = ht4min ! minimum Ht for the four leading jets
 -1 = ht2max ! maximum Ht for the two leading jets
 -1 = ht3max ! maximum Ht for the three leading jets
 -1 = ht4max ! maximum Ht for the four leading jets
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#***********************************************************************
 0 = ptgmin ! Min photon transverse momentum
 0.4 = R0gamma ! Radius of isolation code
 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
.true. = isoEM ! isolate photons from EM energy (photons and leptons)
#*********************************************************************
# WBF cuts *
#*********************************************************************
 0 = xetamin ! minimum rapidity for two jets in the WBF case
 0 = deltaeta ! minimum rapidity for two jets in the WBF case
#*********************************************************************
# KT DURHAM CUT *
#*********************************************************************
 -1 = ktdurham
 0.4 = dparameter
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
# Jet measure cuts *
#*********************************************************************
0 = xqcut ! minimum kt jet measure between partons
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: If use_syst is T, matched Pythia output is *
# meaningful ONLY if plotted taking matchscale *
# reweighting into account! *
#*********************************************************************
   T = use_syst ! Enable systematics studies

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#1

I guess that the issue is that you run MG5aMC like
./bin/mg5_aMC < my_cmd
instead of the recommended method
./bin/mg5_aMC my_cmd

The handling of the "#" in the file my_cmd is not handle in the first case and lead to the error that you quote.
So options are
1) to ask CMS to do the parsing of the "#" in the customize_card.dat (a file not existing within MG5aMC framework)
2) to ask CMS to use the official recomended method to run script
3) to remove the comment in the customize_card.dat

Cheers,

Olivier

PS: we do not garantee retro-compatiblity of the method "./bin/mg5_aMC < my_cmd" with future version of the code but we do for
"./bin/mg5_aMC my_cmd"

Revision history for this message
Sahithi (rudra333) said :
#2

Removing the comments resolved the issue, thank you so much

Revision history for this message
Sahithi (rudra333) said :
#3

Thanks Olivier Mattelaer, that solved my question.