Skip to content

MET Phi correction #74

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
Anusreevijay769 wants to merge 16 commits into
base: HZZ_Analysis_2l2nu
Choose a base branch
from
Open

MET Phi correction #74

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
16 commits
Select commit Hold shift + click to select a range
881a334
updates to condor script
Anusreevijay769 Mar 13, 2025
2dc4e8a
bug fix for emu event selection
Anusreevijay769 May 13, 2025
305c0b8
update for deltaR selection for emu
Anusreevijay769 May 16, 2025
9635ddf
rochester correction
Anusreevijay769 Jun 23, 2025
d64484a
MET Phi correction
Anusreevijay769 Jul 1, 2025
25faa85
MET Phi correction
Anusreevijay769 Jul 2, 2025
3c5c440
updates
Anusreevijay769 Jul 8, 2025
60d7d91
selection of genuine and PU jets
Anusreevijay769 Jul 10, 2025
a2ad97d
updates
Anusreevijay769 Jul 15, 2025
1376c63
fix to leading and subleading jet
Anusreevijay769 Aug 13, 2025
6ad277f
mll window
Anusreevijay769 Aug 18, 2025
88ccd7b
updates
Anusreevijay769 Sep 19, 2025
6efb60f
updates
Anusreevijay769 Oct 7, 2025
437085b
updates
Anusreevijay769 Oct 16, 2025
836d3eb
bug fix in jet variables
Anusreevijay769 Oct 17, 2025
126c4ca
updates to run for nanoAOD v15 samples
Anusreevijay769 Dec 14, 2025
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
65 changes: 48 additions & 17 deletions config/Input_2018.yml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
lumi: 59740.
gumi: 59740.
TriggerChannels:
- Triggers_HZZ4l
- Triggers_HZZ2l2q
Expand Down Expand Up @@ -34,27 +34,58 @@ Triggers_HZZ2l2q:
Triggers_HZZ2l2nu:
# DoubleMu
- event.HLT_Mu17_TrkIsoVVL_Mu8_TrkIsoVVL_DZ_Mass3p8
- event.HLT_Mu17_TrkIsoVVL_Mu8_TrkIsoVVL_DZ_Mass8
- event.HLT_Mu19_TrkIsoVVL_Mu9_TrkIsoVVL_DZ_Mass3p8
- event.HLT_Mu19_TrkIsoVVL_Mu9_TrkIsoVVL_DZ_Mass8
# Single Mu
- event.HLT_IsoMu20
#- event.HLT_IsoMu20
#- event.HLT_IsoMu24
#- event.HLT_IsoMu27
- event.HLT_IsoMu24
- event.HLT_IsoMu27
- event.HLT_IsoMu30
- event.HLT_Mu50
# Double Electron
- event.HLT_Ele23_Ele12_CaloIdL_TrackIdL_IsoVL
- event.HLT_DoubleEle25_CaloIdL_MW
# Single Electron
- event.HLT_Ele32_WPTight_Gsf
#- event.HLT_Ele32_WPTight_Gsf
# MuEGamma
- event.HLT_Mu8_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL_DZ
#- event.HLT_Mu8_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL_DZ
#- event.HLT_Mu23_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL_DZ
- event.HLT_Mu12_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL_DZ
- event.HLT_Mu12_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL
- event.HLT_Mu23_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL_DZ
- event.HLT_Mu23_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL
- event.HLT_Mu8_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL_DZ
- event.HLT_Mu8_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL
# EGamma
- event.HLT_Ele23_Ele12_CaloIdL_TrackIdL_IsoVL_DZ
- event.HLT_Ele23_Ele12_CaloIdL_TrackIdL_IsoVL
- event.HLT_DiEle27_WPTightCaloOnly_L1DoubleEG

#Double electron
- event.HLT_DoubleEle33_CaloIdL_MW
- event.HLT_DoubleEle25_CaloIdL_MW
- event.HLT_DoubleEle27_CaloIdL_MW
- event.HLT_DoublePhoton70
- event.HLT_Ele115_CaloIdVT_GsfTrkIdT
- event.HLT_Ele27_WPTight_Gsf
- event.HLT_Ele32_WPTight_Gsf
- event.HLT_Ele35_WPTight_Gsf
- event.HLT_Ele38_WPTight_Gsf
- event.HLT_Ele40_WPTight_Gsf
- event.HLT_Ele32_WPTight_Gsf_L1DoubleEG
- event.HLT_Photon200
# SinglePhoton
- event.HLT_Photon50_R9Id90_HE10_IsoM
- event.HLT_Photon75_R9Id90_HE10_IsoM
- event.HLT_Photon75_R9Id90_HE10_IsoM_EBOnly_CaloMJJ300_PFJetsMJJ400DEta3
- event.HLT_Photon75_R9Id90_HE10_IsoM_EBOnly_CaloMJJ400_PFJetsMJJ600DEta3
- event.HLT_Photon90_R9Id90_HE10_IsoM
- event.HLT_Photon120_R9Id90_HE10_IsoM
- event.HLT_Photon165_R9Id90_HE10_IsoM
- event.HLT_Photon300_NoHE
#- event.HLT_Photon50_R9Id90_HE10_IsoM
#- event.HLT_Photon75_R9Id90_HE10_IsoM
#- event.HLT_Photon75_R9Id90_HE10_IsoM_EBOnly_CaloMJJ300_PFJetsMJJ400DEta3
#- event.HLT_Photon75_R9Id90_HE10_IsoM_EBOnly_CaloMJJ400_PFJetsMJJ600DEta3
#- event.HLT_Photon90_R9Id90_HE10_IsoM
#- event.HLT_Photon120_R9Id90_HE10_IsoM
#- event.HLT_Photon165_R9Id90_HE10_IsoM
#- event.HLT_Photon300_NoHE


Triggers_HZZ2l2nu_SingleLep:
Expand Down Expand Up @@ -108,8 +139,8 @@ Electron:
pTcut: 25 #7
Etacut: 2.1 #2.5
#Sip3dcut: 4
Loosedxycut: 0.045 #0.5
Loosedzcut: 0.2 #1
Loosedxycut: 0.045 #0.5 #
Loosedzcut: 0.2 #1 #0.2
ttH:
WP: 0.8
BDTWP:
Expand All @@ -129,8 +160,8 @@ Muon:

Etacut: 2.4
#Sip3dcut: 4
Loosedxycut: 0.045 #0.5
Loosedzcut: 0.2 #1
Loosedxycut: 0.045 #0.5 #0.045
Loosedzcut: 0.2 #1 #0.2
Tightdxycut: 0.045 #0.2
Tightdzcut: 0.2 #0.5
TightTrackerLayercut: 5
Expand Down Expand Up @@ -166,7 +197,7 @@ HZZ2l2nu:
SubLeading_Lep_pT: 25.0
Lep_eta: 2.5
Pt_ll: 10.0
M_ll_Window: 15.0
M_ll_Window: 0.0
dPhi_jetMET: 0.0
MZLepcut:
down: 60.0
Expand Down
120 changes: 108 additions & 12 deletions include/H4LTools.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -22,6 +22,7 @@ class H4LTools
float HZZ2l2q_Leading_Lep_pT, HZZ2l2q_SubLeading_Lep_pT, HZZ2l2q_Lep_eta, HZZ2l2q_MZLepcutdown, HZZ2l2q_MZLepcutup;
float HZZ2l2nu_Leading_Lep_pT, HZZ2l2nu_SubLeading_Lep_pT, HZZ2l2nu_Lep_eta, HZZ2l2nu_Pt_ll, HZZ2l2nu_M_ll_Window, HZZ2l2nu_dPhi_jetMET, HZZ2l2nu_MZLepcutdown, HZZ2l2nu_MZLepcutup;
bool DEBUG;
H4LTools(bool isMC_);

void InitializeElecut(float elePtcut_, float eleEtacut_, float eleLoosedxycut_, float eleLoosedzcut_, float eleIsocut_, float eleBDTWPLELP_, float eleBDTWPMELP_, float eleBDTWPHELP_, float eleBDTWPLEHP_, float eleBDTWPMEHP_, float eleBDTWPHEHP_)
{
Expand Down Expand Up @@ -101,7 +102,8 @@ class H4LTools
btag_deepJet_Tight = btag_deepJet_Tight_;
}

void SetElectrons(float Electron_pt_, float Electron_eta_, float Electron_phi_, float Electron_mass_, float Electron_dxy_, float Electron_dz_,
//for v9
/*void SetElectrons(float Electron_pt_, float Electron_eta_, float Electron_phi_, float Electron_mass_, float Electron_dxy_, float Electron_dz_,
float Electron_mvaFall17V2Iso_WP90_, int Electron_pdgId_, float Electron_pfRelIso03_all_)
{
Electron_pt.push_back(Electron_pt_);
Expand All @@ -111,13 +113,30 @@ class H4LTools
Electron_dxy.push_back(Electron_dxy_);
Electron_dz.push_back(Electron_dz_);
//Electron_sip3d.push_back(Electron_sip3d_);
//Electron_mvaFall17V2Iso.push_back(Electron_mvaFall17V2Iso_);
//Electron_mvaFall17V2Iso.push_back(Electron_mvaFall17V2Iso_); // no need
Electron_mvaFall17V2Iso_WP90.push_back(Electron_mvaFall17V2Iso_WP90_);
Electron_pdgId.push_back(Electron_pdgId_);
Electron_pfRelIso03_all.push_back(Electron_pfRelIso03_all_);
}*/

//for v15
void SetElectrons(float Electron_pt_, float Electron_eta_, float Electron_phi_, float Electron_mass_, float Electron_dxy_, float Electron_dz_,
float Electron_mvaIso_WP90_, int Electron_pdgId_, float Electron_pfRelIso03_all_)
{
Electron_pt.push_back(Electron_pt_);
Electron_phi.push_back(Electron_phi_);
Electron_eta.push_back(Electron_eta_);
Electron_mass.push_back(Electron_mass_);
Electron_dxy.push_back(Electron_dxy_);
Electron_dz.push_back(Electron_dz_);
//Electron_sip3d.push_back(Electron_sip3d_);
Electron_mvaIso_WP90.push_back(Electron_mvaIso_WP90_);
Electron_pdgId.push_back(Electron_pdgId_);
Electron_pfRelIso03_all.push_back(Electron_pfRelIso03_all_);
}

void SetJets(float Jet_pt_, float Jet_eta_, float Jet_phi_, float Jet_mass_, int Jet_jetId_, float Jet_btagDeepFlavB_,
// for v9
/*void SetJets(float Jet_pt_, float Jet_eta_, float Jet_phi_, float Jet_mass_, int Jet_jetId_, float Jet_btagDeepFlavB_,
int Jet_puId_)
{
Jet_pt.push_back(Jet_pt_);
Expand All @@ -127,26 +146,61 @@ class H4LTools
Jet_btagDeepFlavB.push_back(Jet_btagDeepFlavB_);
Jet_jetId.push_back(Jet_jetId_);
Jet_puId.push_back(Jet_puId_); // 1 or 0?
}*/

//for v15
void SetJets(float Jet_pt_, float Jet_eta_, float Jet_phi_, float Jet_mass_, float Jet_btagDeepFlavB_,
float Jet_chEmEF_, float Jet_chHEF_, float Jet_neEmEF_, float Jet_neHEF_, float Jet_muEF_, int Jet_nConstituents_, int Jet_neMultiplicity_, int Jet_chMultiplicity_)
{
Jet_pt.push_back(Jet_pt_);
Jet_phi.push_back(Jet_phi_);
Jet_eta.push_back(Jet_eta_);
Jet_mass.push_back(Jet_mass_);
Jet_btagDeepFlavB.push_back(Jet_btagDeepFlavB_);
Jet_chEmEF.push_back(Jet_chEmEF_);
Jet_chHEF.push_back(Jet_chHEF_);
Jet_neEmEF.push_back(Jet_neEmEF_);
Jet_neHEF.push_back(Jet_neHEF_);
Jet_muEF.push_back(Jet_muEF_);
Jet_nConstituents.push_back(Jet_nConstituents_);
Jet_neMultiplicity.push_back(Jet_neMultiplicity_);
Jet_chMultiplicity.push_back(Jet_chMultiplicity_);
}

void SetFatJets(float Jet_pt_, float Jet_eta_, float Jet_phi_, float Jet_mass_, int Jet_jetId_, float Jet_btagDeepB_,
float Jet_PNZvsQCD_)
//void SetFatJets(float Jet_pt_, float Jet_eta_, float Jet_phi_, float Jet_mass_)

{
FatJet_pt.push_back(Jet_pt_);
FatJet_eta.push_back(Jet_eta_);
FatJet_phi.push_back(Jet_phi_);
FatJet_SDmass.push_back(Jet_mass_);
FatJet_jetId.push_back(Jet_jetId_);
FatJet_btagDeepB.push_back(Jet_btagDeepB_);
FatJet_PNZvsQCD.push_back(Jet_PNZvsQCD_); // 1 or 0?
FatJet_jetId.push_back(Jet_jetId_); //commented out for v15
//FatJet_btagDeepB.push_back(Jet_btagDeepB_);
FatJet_PNZvsQCD.push_back(Jet_PNZvsQCD_); // 1 or 0? //commented out for v15
}

void SetGenJets(float GenJet_pt_, float GenJet_eta_, float GenJet_phi_, float GenJet_mass_)
{
GenJet_pt.push_back(GenJet_pt_);
GenJet_eta.push_back(GenJet_eta_);
GenJet_phi.push_back(GenJet_phi_);
GenJet_mass.push_back(GenJet_mass_);
}

void SetMET(float MET_pt_, float MET_phi_, float MET_sumEt_)
{
MET_pt = MET_pt_;
MET_phi = MET_phi_;
MET_sumEt = MET_sumEt_;
// std::cout<<"Inside header file: MET_sumEt = " << MET_sumEt_ << "\t" << MET_sumEt << std::endl;
}

void SetPuppiMET(float MET_pt_, float MET_phi_, float MET_sumEt_)
{
PuppiMET_pt = MET_pt_;
PuppiMET_phi = MET_phi_;
PuppiMET_sumEt = MET_sumEt_;
}

void SetMuons(float Muon_pt_, float Muon_eta_, float Muon_phi_, float Muon_mass_, bool Muon_isGlobal_, bool Muon_isTracker_,
Expand Down Expand Up @@ -222,6 +276,9 @@ class H4LTools
std::vector<float> MuonFsrPhi();
std::vector<unsigned int> SelectedJets(std::vector<unsigned int> ele, std::vector<unsigned int> mu);
std::vector<unsigned int> SelectedFatJets(std::vector<unsigned int> ele, std::vector<unsigned int> mu);
bool PassJetIDv15(unsigned int i, bool isPUPPI);
//std::vector<unsigned int> SelectedJets(std::vector<unsigned int> ele, std::vector<unsigned int> mu, std::vector<unsigned int>& genuineJets, std::vector<unsigned int>& puJets);


std::vector<TLorentzVector> Zlist;
std::vector<TLorentzVector> Zlistnofsr;
Expand Down Expand Up @@ -264,6 +321,7 @@ class H4LTools
bool flag4e;
bool flag4mu;
bool flag2e2mu;
bool muonpt;

bool isBoosted2l2q;
bool flag2e;
Expand Down Expand Up @@ -338,6 +396,7 @@ class H4LTools
Electron_dz.clear();
Electron_sip3d.clear();
Electron_mvaFall17V2Iso_WP90.clear();
Electron_mvaIso_WP90.clear();
Electron_pdgId.clear();
Electron_pfRelIso03_all.clear();
Elechg.clear();
Expand Down Expand Up @@ -382,21 +441,39 @@ class H4LTools
Jet_btagDeepFlavB.clear();
Jet_jetId.clear();
Jet_puId.clear();
Jet_chEmEF.clear();
Jet_chHEF.clear();
Jet_neEmEF.clear();
Jet_neHEF.clear();
Jet_muEF.clear();
Jet_nConstituents.clear();
Jet_neMultiplicity.clear();
Jet_chMultiplicity.clear();
FatJet_pt.clear();
FatJet_phi.clear();
FatJet_eta.clear();
FatJet_SDmass.clear();
FatJet_btagDeepB.clear();
//FatJet_btagDeepB.clear();
FatJet_PNZvsQCD.clear();
FatJet_jetId.clear();
GenJet_pt.clear();
GenJet_eta.clear();
GenJet_phi.clear();
GenJet_mass.clear();
jetidx.clear();
FatJetidx.clear();


// MET related variables
MET_pt = -999.0;
MET_phi = -999.0; ////new
MET_sumEt = -999.0;

// Puppi MET related variables
PuppiMET_pt = -999.0;
PuppiMET_phi = -999.0;
PuppiMET_sumEt = -999.0;

// FsrPhoton related variables
nFsrPhoton = 0;
FsrPhoton_dROverEt2.clear();
Expand All @@ -408,6 +485,7 @@ class H4LTools
// Generator related variables
nGenPart = 0;
GenPart_pt.clear();


// Reconstructed variables
Zlist.clear();
Expand Down Expand Up @@ -467,6 +545,7 @@ class H4LTools
etaL2_emu = -999.0;
phiL2_emu = -999.0;
massL2_emu = -999.0;
DeltaRl1l2 = -999.0;

pTj1 = -999.0;
etaj1 = -999.0;
Expand Down Expand Up @@ -502,6 +581,7 @@ class H4LTools
flag2l = false;
HZZ2l2qNu_isELE = false;
HZZ2l2qNu_cutOppositeChargeFlag = false;
muonpt = false;

HZZ2l2nu_flag2e_met = false;
HZZ2l2nu_flag2l_met = false;
Expand Down Expand Up @@ -569,10 +649,11 @@ class H4LTools
int cutMETlt150;
int HZZ2l2nu_cutMETgT100;
int HZZ2l2nu_cut2l_met_m40_180, HZZ2l2nu_cut2e_met_m40_180, HZZ2l2nu_cut2mu_met_m40_180;
int cut2e, cut2mu, cut2l, cut2l1J, cut2l2j, cut2l1Jor2j;
int cut2e, cut2mu, cut2l, cut2l1J, cut2l2j, cut2l1Jor2j, cut_mu_pt, cut_mu_eta, cut_mu_mediumid, cut_mu_isglobal_istracker, cut_mu_iso;
int cut_2mu_cutOppositeCharge;
int HZZ2l2nu_cut2e_met, HZZ2l2nu_cut2mu_met, HZZ2l2qNu_cut2l;
int cut4e, cut4mu, cut2e2mu, cutZZ4e, cutZZ4mu, cutZZ2e2mu, cutm4l4e, cutm4l4mu, cutm4l2e2mu, cutghost2e2mu, cutQCD2e2mu, cutLepPt2e2mu, cutghost4e, cutQCD4e, cutLepPt4e, cutghost4mu, cutQCD4mu, cutLepPt4mu;
float pTL1, etaL1, phiL1, massL1, pTL2, etaL2, phiL2, massL2, pTL3, etaL3, phiL3, massL3, pTL4, etaL4, phiL4, massL4;
float pTL1, etaL1, phiL1, massL1, pTL2, etaL2, phiL2, massL2, pTL3, etaL3, phiL3, massL3, pTL4, etaL4, phiL4, massL4, DeltaRl1l2;
float pTL1_emu, etaL1_emu, phiL1_emu, massL1_emu, pTL2_emu, etaL2_emu, phiL2_emu, massL2_emu;
float pTj1, etaj1, phij1, mj1, pTj2, etaj2, phij2, mj2;
int HZZ2l2qNu_cutOppositeCharge;
Expand All @@ -592,18 +673,26 @@ class H4LTools
int HZZ_emuCR_cutdPhiJetMET;
int HZZ_emuCR_cutMETgT100;


private:
std::vector<float> Electron_pt, Electron_phi, Electron_eta, Electron_mass, Electron_dxy, Electron_dz, Electron_sip3d;
std::vector<float> Electron_mvaFall17V2Iso_WP90, Electron_pfRelIso03_all;
std::vector<float> Electron_mvaFall17V2Iso_WP90, Electron_pfRelIso03_all, Electron_mvaIso_WP90;
std::vector<int> Electron_pdgId;
bool isMC;

std::vector<float> Jet_pt, Jet_phi, Jet_eta, Jet_mass, Jet_btagDeepFlavB;
std::vector<int> Jet_jetId, Jet_puId;
std::vector<float> Jet_chEmEF, Jet_chHEF, Jet_neEmEF, Jet_neHEF, Jet_muEF;
std::vector<int> Jet_nConstituents, Jet_neMultiplicity, Jet_chMultiplicity;
float MET_pt, MET_phi;
float MET_sumEt, MT_2l2nu;
float PuppiMET_pt, PuppiMET_phi;
float PuppiMET_sumEt;

std::vector<float> FatJet_pt, FatJet_phi, FatJet_eta, FatJet_SDmass, FatJet_btagDeepB, FatJet_PNZvsQCD;
std::vector<float> FatJet_pt, FatJet_phi, FatJet_eta, FatJet_SDmass, FatJet_PNZvsQCD;
//std::vector<float> FatJet_btagDeepB;
std::vector<int> FatJet_jetId;
std::vector<float> GenJet_pt, GenJet_eta, GenJet_phi, GenJet_mass;

std::vector<float> Muon_pt, Muon_phi, Muon_eta, Muon_mass, Muon_dxy, Muon_dz, Muon_sip3d, Muon_ptErr, Muon_pfRelIso03_all;
std::vector<int> Muon_nTrackerLayers, Muon_genPartIdx, Muon_pdgId, Muon_charge;
Expand Down Expand Up @@ -676,6 +765,12 @@ H4LTools::H4LTools(int year, bool DEBUG_Main)
cut2e_m40_180 = 0;
cut2mu_m40_180 = 0;
cut2l_m40_180 = 0;
cut_mu_pt = 0;
cut_mu_eta = 0;
cut_mu_mediumid = 0;
cut_mu_isglobal_istracker = 0;
cut_mu_iso = 0;
cut_2mu_cutOppositeCharge = 0;

HZZ2l2nu_cut2e_met = 0;
HZZ2l2nu_cut2mu_met = 0;
Expand All @@ -692,5 +787,6 @@ H4LTools::H4LTools(int year, bool DEBUG_Main)
HZZ_emuCR_cutbtag = 0;
HZZ_emuCR_cutdPhiJetMET = 0;
HZZ_emuCR_cutMETgT100 = 0;

}
#endif
Loading