h_2DHisto_meas.Fill( FatJet.Perp(), FatJet.M(), weight )
h_2DHisto_measSD.Fill( FatJetSD.Perp(), FatJetSD.M(), weight)
# h_m_fulldist.Fill(FatJet.M(), weight)
# h_pt_fulldist.Fill(FatJet.Perp(), weight)
# h_pt_fulldist_softdrop.Fill(FatJetSD.Perp(), weight)
# h_m_fulldist_softdrop.Fill(FatJetSD.M(), weight)
igen = getMatched( FatJet, GenJets )
igenSD = getMatched(FatJetSD, GenJetsSD, dRMax=0.5)
if igen != None : # Here we have a "Fill"
valup = getJER(FatJet.Eta(), +1) #JER nominal=0, up=+1, down=-1
recopt = FatJet.Perp()
genpt = GenJets[igen].Perp()
deltapt = (recopt-genpt)*(valup-1.0)
smearup = max(0.0, (recopt+deltapt)/recopt)
valdn = getJER(FatJet.Eta(), -1) #JER nominal=0, dn=+1, down=-1
recopt = FatJet.Perp()
genpt = GenJets[igen].Perp()
deltapt = (recopt-genpt)*(valdn-1.0)
smeardn = max(0.0, (recopt+deltapt)/recopt)
response.Fill( FatJet.Perp(), FatJet.M(), GenJets[igen].Perp(), GenJets[igen].M(), weight )
response_jecup.Fill( FatJet.Perp() * FatJetCorrUp[ijet], FatJet.M()* FatJetCorrUp[ijet], GenJets[igen].Perp(), GenJets[igen].M(), weight )
response_jecdn.Fill( FatJet.Perp() * FatJetCorrDn[ijet], FatJet.M()* FatJetCorrDn[ijet], GenJets[igen].Perp(), GenJets[igen].M(), weight )
response_jerup.Fill( FatJet.Perp() * smearup, FatJet.M()* smearup, GenJets[igen].Perp(), GenJets[igen].M(), weight )
FatJetsSD.append(FatJetSD)
FatJets.append(FatJet)
h_2DHisto_meas.Fill(FatJet.M(), FatJet.Perp(), weight)
h_2DHisto_measSD.Fill(FatJetSD.M(), FatJetSD.Perp(), weight)
# h_m_fulldist.Fill(FatJet.M(), weight)
# h_pt_fulldist.Fill(FatJet.Perp(), weight)
# h_pt_fulldist_softdrop.Fill(FatJetSD.Perp(), weight)
# h_m_fulldist_softdrop.Fill(FatJetSD.M(), weight)
igen = getMatched(FatJet, GenJets)
igenSD = getMatched(FatJetSD, GenJetsSD, dRMax=0.5)
if igen != None: # Here we have a "Fill"
valup = getJER(FatJet.Eta(),
+1) #JER nominal=0, up=+1, down=-1
recopt = FatJet.Perp()
genpt = GenJets[igen].Perp()
deltapt = (recopt - genpt) * (valup - 1.0)
smearup = max(0.0, (recopt + deltapt) / recopt)
valdn = getJER(FatJet.Eta(),
-1) #JER nominal=0, dn=+1, down=-1
recopt = FatJet.Perp()
genpt = GenJets[igen].Perp()
deltapt = (recopt - genpt) * (valdn - 1.0)
smeardn = max(0.0, (recopt + deltapt) / recopt)
response.Fill(FatJet.M(), FatJet.Perp(), GenJets[igen].M(),
GenJets[igen].Perp(), weight)
response_jecup.Fill(FatJet.M() * FatJetCorrUp[ijet],
FatJetsSD.append(FatJetSD)
FatJets.append(FatJet)
h_2DHisto_meas.Fill( FatJet.Perp(), FatJet.M(), weight )
h_2DHisto_measSD.Fill( FatJetSD.Perp(), FatJetSD.M(), weight)
igen = getMatched( FatJet, GenJets )
igenSD = getMatched(FatJetSD, GenJetsSD, dRMax=0.5)
if igen != None : # Here we have a "Fill"
valup = getJER(FatJet.Eta(), +1) #JER nominal=0, up=+1, down=-1
recopt = FatJet.Perp()
genpt = GenJets[igen].Perp()
deltapt = (recopt-genpt)*(valup-1.0)
smearup = max(0.0, (recopt+deltapt)/recopt)
valdn = getJER(FatJet.Eta(), -1) #JER nominal=0, dn=+1, down=-1
recopt = FatJet.Perp()
genpt = GenJets[igen].Perp()
deltapt = (recopt-genpt)*(valdn-1.0)
smeardn = max(0.0, (recopt+deltapt)/recopt)
valnom = getJER(FatJet.Eta(), 0)
recopt = FatJet.Perp()
genpt = GenJets[igen].Perp()
deltapt = (recopt-genpt)*(valnom-1.0)