Skip to content

add minus_path_weights to path weights #68

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

Open
dz24 wants to merge 5 commits into
base: main
Choose a base branch
from
Open

add minus_path_weights to path weights #68

Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
bb8b36b
first
dz24 Mar 10, 2026
51d24cf
update
dz24 Mar 10, 2026
1c2ec60
update
dz24 Mar 10, 2026
6fa0ebd
push
dz24 Mar 10, 2026
b7d9f3f
Refactor minus_paths assignment for clarity
lukasbaldauf Mar 10, 2026
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
14 changes: 13 additions & 1 deletion inftools/tistools/path_weights.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10,7 +10,8 @@ def get_path_weights(
nskip: Annotated[ int, typer.Option( "-nskip", help="Skip the first nskip entries of the infretis_data.txt file",), ] = 0,
plotP: Annotated[ bool, typer.Option( "-plotP", help="If true plot the binless crossing probability"), ] = False,
outP: Annotated[ str, typer.Option( "-outP", help="Write the binless WHAM crossing probability to outP"), ] = "",
overw: Annotated[ bool, typer.Option("-O", help="Force overwriting of files") ] = False,):
overw: Annotated[ bool, typer.Option("-O", help="Force overwriting of files") ] = False,
minus: Annotated[ bool, typer.Option("-minus", help="Add the minus path weigths") ] = False):
"""Calculate the unbiased weight for paths in the plus ensembles.

The weights can be used to calculate observables as <O> = sum(wi*Oi), for
Expand Down Expand Up @@ -40,6 +41,7 @@ def get_path_weights(
data = data[nskip:]
# we only need non-zero paths
non_zero_paths = data[:, 3] == "----"
minus_paths = ~non_zero_paths
data[data == "----"] = "0.0"
D = {}

Expand Down Expand Up @@ -85,6 +87,16 @@ def get_path_weights(
)
A[j] = Q[K] * np.sum(w[j])

# add minus to weights
if minus:
minus_pn = data[minus_paths, 0:1].astype(int)
denominator = np.sum(data[minus_paths, 3].astype(float))
minus_op = data[minus_paths, 2:3].astype(float)
minus_w = data[minus_paths, 3:4].astype(float)/denominator
D["pnr"] = np.concatenate([D["pnr"], minus_pn])
D["maxop"] = np.concatenate([D["maxop"], minus_op])
A = np.concatenate([A, minus_w])

print(f"Weights saved to {out}.")
np.savetxt(
out,
Expand Down