Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

simplify quark mass effect #3778

Merged
merged 1 commit into from
Oct 10, 2024
Merged

simplify quark mass effect #3778

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
81 changes: 44 additions & 37 deletions bin/utils/request_fragment_check.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -487,6 +487,46 @@ def vbf_dipole_recoil_check(vbf_lo,vbf_nlo,data_f2,pw_gp,dn):
print("[OK] VBF POWHEG with local recoil --> SpaceShower:dipoleRecoil = 1.")
return warning_dipole, error_dipole

def powheg_gg_H_quark_mass_effects():
warning_gg_H_quark_mass_effects = []
error_gg_H_quark_mass_effects = []
#for more information on this check, see
#https://its.cern.ch/jira/browse/CMSCOMPPR-4874
#this configuration is ok at 125 GeV, but causes trouble starting at around 170:
# ncall1=50000, itmx1=5, ncall2=50000, itmx2=5, foldcsi=1, foldy=1, foldphi=1
#from mH=300 GeV to 3 TeV, this configuration seems to be fine:
# ncall1=550000, itmx1=7, ncall2=75000, itmx2=5, foldcsi=2, foldy=5, foldphi=2
#I'm printing warnings here for anything less than the second configuration.
#Smaller numbers are probably fine at low mass
desiredvalues = {
"ncall1": 550000,
"itmx1": 7,
"ncall2": 75000,
"itmx2": 5,
"foldcsi": 2,
"foldy": 5,
"foldphi": 2,
}
if et_flag == 0 and et_flag_external == 0:
with open(os.path.join(my_path, pi, "powheg.input")) as f:
content = f.read()
matches = dict((name, re.search(r"^"+name+" *([0-9]+)", content, flags=re.MULTILINE)) for name in desiredvalues)
if et_flag == 1 and et_flag_external == 0:
with open(os.path.join(my_path, pi, "external_tarball/powheg.input")) as f:
content = f.read()
matches = dict((name, re.search(r"^"+name+" *([0-9]+)", content, flags=re.MULTILINE)) for name in desiredvalues)
bad = False
for name, match in matches.items():
if match:
actualvalue = int(match.group(1))
if actualvalue < desiredvalues[name]:
bad = True
warning_gg_H_quark_mass_effects.append("{0} = {1}, should be at least {2} (may be ok if hmass < 150 GeV, please check!)".format(name, actualvalue, desiredvalues[name]))
else:
bad = True
error_gg_H_quark_mass_effects.append("didn't find "+name+" in powheg.input")
if not bad: print("[OK] integration grid setup looks ok for gg_H_quark_mass_effects")
return warning_gg_H_quark_mass_effects, error_gg_H_quark_mass_effects

if args.dev:
print("Running on McM DEV!\n")
Expand Down Expand Up @@ -1322,7 +1362,7 @@ def root_requests_from_ticket(ticket_prepid, include_docs=False):
with open(pwg_stat_file) as f_pwg_stat:
s_pwg_stat = f_pwg_stat.read()
print("-----------------------------------------------------------------")
print("Summary from pwg-stat.dat from Powheg firdpack (for experts only):")
print("Summary from pwg-stat.dat from Powheg gridpack (for experts only):")
print("-----------------------------------------------------------------")
print(s_pwg_stat)
print("-----------------------------------------------------------------")
Expand Down Expand Up @@ -1357,42 +1397,9 @@ def root_requests_from_ticket(ticket_prepid, include_docs=False):
if match:
process = match.group(2)
if process == "gg_H_quark-mass-effects":
#for more information on this check, see
#https://its.cern.ch/jira/browse/CMSCOMPPR-4874
#this configuration is ok at 125 GeV, but causes trouble starting at around 170:
# ncall1=50000, itmx1=5, ncall2=50000, itmx2=5, foldcsi=1, foldy=1, foldphi=1
#from mH=300 GeV to 3 TeV, this configuration seems to be fine:
# ncall1=550000, itmx1=7, ncall2=75000, itmx2=5, foldcsi=2, foldy=5, foldphi=2
#I'm printing warnings here for anything less than the second configuration.
#Smaller numbers are probably fine at low mass
desiredvalues = {
"ncall1": 550000,
"itmx1": 7,
"ncall2": 75000,
"itmx2": 5,
"foldcsi": 2,
"foldy": 5,
"foldphi": 2,
}
if et_flag == 0 and et_flag_external == 0:
with open(os.path.join(my_path, pi, "powheg.input")) as f:
content = f.read()
matches = dict((name, re.search(r"^"+name+" *([0-9]+)", content, flags=re.MULTILINE)) for name in desiredvalues)
if et_flag == 1 and et_flag_external == 0:
with open(os.path.join(my_path, pi, "external_tarball/powheg.input")) as f:
content = f.read()
matches = dict((name, re.search(r"^"+name+" *([0-9]+)", content, flags=re.MULTILINE)) for name in desiredvalues)
bad = False
for name, match in matches.items():
if match:
actualvalue = int(match.group(1))
if actualvalue < desiredvalues[name]:
bad = True
warnings.append("{0} = {1}, should be at least {2} (may be ok if hmass < 150 GeV, please check!)".format(name, actualvalue, desiredvalues[name]))
else:
bad = True
errors.append("didn't find "+name+" in powheg.input")
if not bad: print("[OK] integration grid setup looks ok for gg_H_quark-mass-effects")
warn_tmp , err_tmp = powheg_gg_H_quark_mass_effects()
warnings.extend(warn_tmp)
errors.extend(err_tmp)
else:
warnings.append("Didn't find powheg process in runcmsgrid.sh")

Expand Down