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.

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

Vbf datacards systs tt exploration #35

Closed
wants to merge 13 commits into from

Conversation

andresnava1000
Copy link
Collaborator

Edits to bbVVSkimmer where we keep pileup weights (has the vbf postprocessing edits in there too. (built off of that))

@@ -452,7 +454,7 @@ def process(self, events: ak.Array):
msds = jmsr_shifted_vars["msoftdrop"][shift]
pnetms = jmsr_shifted_vars["particleNet_mass"][shift]

if self._save_all:
if self._save_all & False:
Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do we still want this?

@@ -634,6 +636,14 @@ def process(self, events: ak.Array):

# TODO: need to be careful about the sum of gen weights used for the LHE/QCDScale uncertainties
logger.debug("weights ", weights._weights.keys())

# TEMP: save each individual weight
skimmed_events["L1EcalPrefiring"] = ak.to_numpy(events.L1PreFiringWeight.Nom)
Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Will this not be saved from the "single_weight" for loop below?

skimmed_events["L1EcalPrefiring"] = ak.to_numpy(events.L1PreFiringWeight.Nom)

for key in weights._weights.keys():
print(f"single_weight_{key}")
Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think the print() can be removed here.

@@ -1014,7 +1024,7 @@ def getVBFVars(
# int list representing the number of passing vbf jets per event
vbfVars[f"nGoodVBFJets"] = n_good_vbf_jets_sorted_pt.to_numpy()

adding_bdt_vars = True
adding_bdt_vars = False
Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I changed this to false - can you confirm that we don't need these in the processor since they can be added in postprocessing?

@@ -387,6 +387,8 @@ def process(self, events: ak.Array):
)

skimmed_events["nGoodVBFJets"] = np.array(ak.sum(vbf_jet_mask, axis=1))
# Temp
skimmed_events["gen_weights"] = gen_weights
Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Will this not be saved from the "single_weight" for loop below?

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think this branch has a lot of bad "quick fixes" for saving the weight information which should not be kept otherwise. So would we just erase the pull request?

Yea, that part is redundant I think.

@cmantill cmantill closed this Nov 17, 2023
@rkansal47 rkansal47 deleted the vbf_datacards_systs_TT_exploration branch January 11, 2024 21:28
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants