250 |
|
StH_nPVs = new TH1F(("nPVs"+suffix).c_str(),("Number of reconstructed PV ("+suffix+")").c_str(), bin_npvs, min_npvs, max_npvs ); |
251 |
|
StH_addJets = new TH1F(("addJets"+suffix).c_str(),("Number of additional AK5 ("+suffix+")").c_str(), bin_addJets, min_addJets, max_addJets ); |
252 |
|
|
253 |
+ |
StH_typeImet_et = new TH1F(("typeIMET_et"+suffix).c_str(),("Type corrected missing transverse energy ("+suffix+")").c_str(), bin_pt, min_pt, max_pt*0.5 ); |
254 |
|
StH_met_et = new TH1F(("MET_et"+suffix).c_str(),("Missing transverse energy ("+suffix+")").c_str(), bin_pt, min_pt, max_pt*0.5 ); |
255 |
|
StH_met_sig = new TH1F(("MET_sig"+suffix).c_str(),("MET significance("+suffix+")").c_str(), bin_pt, min_pt, max_pt*0.1 ); |
256 |
|
|
279 |
|
StH_simpleJet2_vtxMass = new TH1F(("SimpleJet2_vtxMass"+suffix).c_str(),("Simple Jet2 vtxMass ("+suffix+")").c_str(), bin_SVmass, min_SVmass, max_SVmass ); |
280 |
|
|
281 |
|
StH_simpleJets_bTag = new TH1F(("SimpleJets_bTag"+suffix).c_str(),("Simple Jets bTag ("+suffix+")").c_str(), bin_btag, min_btag, max_btag ); |
282 |
+ |
StH_simpleJets_bTagSum = new TH1F(("SimpleJets_bTagSum"+suffix).c_str(),("Simple Jets bTagSum ("+suffix+")").c_str(), bin_btag*2, min_btag, max_btag*2 ); |
283 |
|
StH_simpleJets_vtxMass = new TH1F(("SimpleJets_vtxMass"+suffix).c_str(),("Simple Jets vtxMass ("+suffix+")").c_str(), bin_SVmass, min_SVmass, max_SVmass ); |
284 |
|
|
285 |
|
StH_addJet1_pt = new TH1F(("AddJet1_pt"+suffix).c_str(),("Additional Jet1 pt ("+suffix+")").c_str(), bin_pt, min_pt, max_pt*0.8 ); |
295 |
|
StH_HdPhi = new TH1F(("HdPhi"+suffix).c_str(),("Higgs deltaPhi ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi ); |
296 |
|
StH_HdEta = new TH1F(("HdEta"+suffix).c_str(),("Higgs deltaEta ("+suffix+")").c_str(), bin_deltaEta, min_deltaEta, max_deltaEta ); |
297 |
|
StH_HHel = new TH1F(("HiggsHel"+suffix).c_str(),("Higgs helicity angle ("+suffix+")").c_str(), bin_hel, min_hel, max_hel ); |
298 |
< |
StH_HPullAngle = new TH1F(("HiggsPullAngle"+suffix).c_str(),("Higgs pull angle ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi ); |
298 |
> |
StH_HPullAngle1 = new TH1F(("HiggsPullAngle1"+suffix).c_str(),("Higgs pull angle1 ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi ); |
299 |
> |
StH_HPullAngle2 = new TH1F(("HiggsPullAngle2"+suffix).c_str(),("Higgs pull angle2 ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi ); |
300 |
> |
StH_HPullAK7Angle1 = new TH1F(("HiggsPullAK7Angle1"+suffix).c_str(),("Higgs pullAK7 angle1 ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi ); |
301 |
> |
StH_HPullAK7Angle2 = new TH1F(("HiggsPullAK7Angle2"+suffix).c_str(),("Higgs pullAK7 angle2 ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi ); |
302 |
|
StH_HdR_addJets = new TH1F(("HdR_addJets"+suffix).c_str(),("Higgs deltaR addJets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR ); |
303 |
|
StH_HdR_leading_addJets = new TH1F(("HdR_leading_addJets"+suffix).c_str(),("Higgs deltaR leading-add Jets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR ); |
304 |
|
StH_HdR_following_addJets = new TH1F(("HdR_following_addJets"+suffix).c_str(),("Higgs deltaR following-add Jets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR ); |
326 |
|
StH_nPVs->Fill(iEvent.nPVs, w); |
327 |
|
StH_addJets->Fill(iEvent.CountAddJets(), w); |
328 |
|
|
329 |
< |
StH_met_et->Fill( TMath::Min(iEvent.METnoPUCh_et, iEvent.METnoPU_et ), w); |
329 |
> |
StH_typeImet_et->Fill(iEvent.typeIcorrMET(), w); |
330 |
> |
// StH_met_et->Fill( TMath::Min(iEvent.METnoPUCh_et, iEvent.METnoPU_et ), w); |
331 |
> |
StH_met_et->Fill( iEvent.MET_et, w); |
332 |
|
StH_met_sig->Fill(iEvent.MET_sig, w); |
333 |
|
|
334 |
|
StH_vectorLepton1_pt->Fill(iEvent.vLepton_pt[0], w); |
362 |
|
StH_simpleJet2_vtxMass->Fill(iEvent.hJet_vtxMass[1], w); |
363 |
|
} |
364 |
|
//here I fill both jets in one histo |
365 |
+ |
StH_simpleJets_bTagSum->Fill(iEvent.hJet_csv[0]+iEvent.hJet_csv[1], w); |
366 |
|
StH_simpleJets_bTag->Fill(iEvent.hJet_csv[0], w); |
367 |
|
StH_simpleJets_bTag->Fill(iEvent.hJet_csv[1], w); |
368 |
|
StH_simpleJets_vtxMass->Fill(iEvent.hJet_vtxMass[0], w); |
375 |
|
StH_HdEta->Fill(iEvent.H_dEta, w); |
376 |
|
|
377 |
|
StH_HHel->Fill(iEvent.hJet_cosTheta[0], w); |
378 |
< |
StH_HPullAngle->Fill(iEvent.deltaPullAngle, w); |
378 |
> |
|
379 |
> |
StH_HPullAngle1->Fill(iEvent.deltaPullAngle, w); |
380 |
> |
StH_HPullAngle2->Fill(iEvent.deltaPullAngle2, w); |
381 |
> |
StH_HPullAK7Angle1->Fill(iEvent.deltaPullAngleAK7, w); |
382 |
> |
StH_HPullAK7Angle2->Fill(iEvent.deltaPullAngle2AK7, w); |
383 |
|
|
384 |
|
bool debug_ = false; |
385 |
|
TLorentzVector higgs, hJet1, hJet2, addJet; |
457 |
|
|
458 |
|
TH1F * StH_nPVs; |
459 |
|
TH1F * StH_addJets; |
460 |
+ |
TH1F * StH_typeImet_et; |
461 |
|
TH1F * StH_met_et; |
462 |
|
TH1F * StH_met_sig; |
463 |
|
|
486 |
|
TH1F * StH_simpleJet2_vtxMass; |
487 |
|
|
488 |
|
TH1F * StH_simpleJets_bTag; |
489 |
+ |
TH1F * StH_simpleJets_bTagSum; |
490 |
|
TH1F * StH_simpleJets_vtxMass; |
491 |
|
|
492 |
|
TH1F * StH_addJet1_pt; |
503 |
|
TH1F * StH_HdEta; |
504 |
|
TH1F * StH_HHel; |
505 |
|
TH1F * StH_HPullAngle; |
506 |
+ |
TH1F * StH_HPullAngle2; |
507 |
+ |
TH1F * StH_HPullAK7Angle; |
508 |
+ |
TH1F * StH_HPullAK7Angle2; |
509 |
|
TH1F * StH_HdR_addJets; |
510 |
|
TH1F * StH_HdR_leading_addJets; |
511 |
|
TH1F * StH_HdR_following_addJets; |