Skip to content

First ana prod from central nano bug fixes#40

Open
valeriadamante wants to merge 75 commits intocms-flaf:mainfrom
valeriadamante:firstAnaProdFromCentralNano_bugFixes
Open

First ana prod from central nano bug fixes#40
valeriadamante wants to merge 75 commits intocms-flaf:mainfrom
valeriadamante:firstAnaProdFromCentralNano_bugFixes

Conversation

@valeriadamante
Copy link
Contributor

This PR is intended to apply bug fixes during the first anaTuple production from the Central NanoAOD. Some of them were indeed not tested (as all DY samples etc) that is the reason of the existence of this PR.

@valeriadamante
Copy link
Contributor Author

@cms-flaf-bot please test:

  • H_mumu_version=PR_40

@cms-flaf-bot
Copy link

pipeline#14077617 started

@cms-flaf-bot
Copy link

pipeline#14077617 failed

@valeriadamante
Copy link
Contributor Author

@cms-flaf-bot please test

@cms-flaf-bot
Copy link

pipeline#14084955 started

@cms-flaf-bot
Copy link

pipeline#14084955 passed

category_definition:
no_cuts: "true"
baseline_WPIso_MM: "mu1_pfRelIso04_all < 0.2 && mu2_pfRelIso04_all< 0.2 && mu1_mediumId && mu2_mediumId"
mu_pt_sel: "mu1_{mu_pt_for_selection} > 20 && mu2_{mu_pt_for_selection} > 20"
Copy link
Collaborator

Choose a reason for hiding this comment

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

Where is this used?

- DY
- TT
- W_NJets
# - VV
Copy link
Collaborator

Choose a reason for hiding this comment

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

  • clean up commented processes
  • Don't you need data for WP studies?

return legs

def PassGenAcceptance(df):
def PassGenAcceptance(df):code
Copy link
Collaborator

Choose a reason for hiding this comment

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

Suggested change
def PassGenAcceptance(df):code
def PassGenAcceptance(df):

df = df.Filter("Muon_idx[Muon_B0].size()==2", "Exactly 2 muons")
df = df.Define("mu1_idx", "Muon_idx[Muon_B0][0]")
df = df.Define("mu2_idx", "Muon_idx[Muon_B0][1]")
df = df.Filter("Muon_charge[mu1_idx]*Muon_charge[mu2_idx]<0", "OS muons")
Copy link
Collaborator

Choose a reason for hiding this comment

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

Is the OS cut really needed at this level?

# # exactly two muons -- See AN/2019_185 line 118 and AN/2019_205 lines 246
df = df.Filter("Muon_idx[Muon_B0].size()==2", "No extra muons")

df = df.Filter("Muon_idx[Muon_B0].size()==2", "Exactly 2 muons")
Copy link
Collaborator

Choose a reason for hiding this comment

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

If you don't apply any muonID and iso requirement, you shouldn't be applying a veto at this stage. Applying a veto on all muons can lead to unwanted behaviour.

df = df.Define(
"JetTagSel",
"Jet_p4[Jet_NoOverlapWithMuons && Jet_btag_Veto_medium].size() < 1 && Jet_p4[Jet_NoOverlapWithMuons && Jet_btag_Veto_loose].size() < 2 && No_Jets_in_dead_Zone",
"Jet_p4[Jet_NoOverlapWithMuons && Jet_btag_Veto_medium].size() < 1 && Jet_p4[Jet_NoOverlapWithMuons && Jet_btag_Veto_loose].size() < 2 ", # && No_Jets_in_dead_Zone",
Copy link
Collaborator

Choose a reason for hiding this comment

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

Why do you remove the dead zone veto?

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.

3 participants