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

Fix bug in sv operator cards #193

Merged
merged 5 commits into from
Sep 3, 2024
Merged

Fix bug in sv operator cards #193

Changes from 3 commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
eb7b355
Fix bug in sv operator cards
andreab1997 Sep 2, 2024
e73f25b
Other error of the same kind
andreab1997 Sep 2, 2024
a8501e1
Add xir*xir to the mur2_grid in the evolution
andreab1997 Sep 2, 2024
759a05e
Revert fix to xif and use a variable for ren_grid2
andreab1997 Sep 3, 2024
dc8eb9b
use ren_grid2 also for as
giacomomagni Sep 3, 2024
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
8 changes: 4 additions & 4 deletions src/pineko/evolve.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -197,7 +197,7 @@ def write_operator_card(pineappl_grid, default_card, card_path, tcard):
muf2_grid = evol_info.fac1
operators_card = copy.deepcopy(default_card)
sv_method = sv_scheme(tcard)
xif = 1.0 if sv_method is not None else tcard["XIF"]
xif = 1.0 if sv_method is None else tcard["XIF"]
# update scale variation method
operators_card["configs"]["scvar_method"] = sv_method

Expand Down Expand Up @@ -354,7 +354,7 @@ def evolve_grid(
evol_info = grid.evolve_info(order_mask)
x_grid = evol_info.x1
mur2_grid = evol_info.ren1
xif = 1.0 if operators1.operator_card.configs.scvar_method is not None else xif
xif = 1.0 if operators1.operator_card.configs.scvar_method is None else xif
andreab1997 marked this conversation as resolved.
Show resolved Hide resolved
tcard = operators1.theory_card
opcard = operators1.operator_card
# rotate the targetgrid
Expand Down Expand Up @@ -426,15 +426,15 @@ def prepare(operator):
fktable = grid.evolve_with_slice_iter2(
prepare(operators1),
prepare(operators2),
ren1=mur2_grid,
ren1=xir*xir*mur2_grid,
andreab1997 marked this conversation as resolved.
Show resolved Hide resolved
alphas=alphas_values,
xi=(xir, xif),
order_mask=order_mask,
)
else:
fktable = grid.evolve_with_slice_iter(
prepare(operators1),
ren1=mur2_grid,
ren1=xir*xir*mur2_grid,
alphas=alphas_values,
xi=(xir, xif),
order_mask=order_mask,
Expand Down
Loading