ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbb/plugins/Histos.h
(Generate patch)

Comparing UserCode/VHbb/plugins/Histos.h (file contents):
Revision 1.2 by bortigno, Wed May 2 19:18:19 2012 UTC vs.
Revision 1.4 by bortigno, Fri May 11 09:46:37 2012 UTC

# Line 250 | Line 250 | public:
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  
# Line 278 | Line 279 | public:
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 );
# Line 295 | Line 297 | public:
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 );
299      StH_HdR_addJets = new TH1F(("HdR_addJets"+suffix).c_str(),("Higgs deltaR addJets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
300 +    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 );
301 +    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 );
302 +    StH_HdRnorm_addJets = new TH1F(("HdRnorm_addJets"+suffix).c_str(),("Higgs normalized deltaR addJets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
303      StH_HMass_addJets = new TH1F(("HiggsMass_addJets"+suffix).c_str(),(" Higgs Mass with addJet ("+suffix+")").c_str(), bin_mass, min_mass, max_mass );
304 +    StH_HMass_leading_addJets = new TH1F(("HiggsMass_leading_addJets"+suffix).c_str(),(" Higgs Mass with leading addJet ("+suffix+")").c_str(), bin_mass, min_mass, max_mass );
305 +    StH_HMass_following_addJets = new TH1F(("HiggsMass_following_addJets"+suffix).c_str(),(" Higgs Mass with following addJet ("+suffix+")").c_str(), bin_mass, min_mass, max_mass );
306 +    StH_HMass_both_addJets = new TH1F(("HiggsMass_both_addJets"+suffix).c_str(),(" Higgs Mass with both addJet ("+suffix+")").c_str(), bin_mass, min_mass, max_mass );
307  
308      StH_massDrop = new TH1F(("MassDrop"+suffix).c_str(),(" Higgs Mass Drop ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
309      StH_y12 = new TH1F(("Y12"+suffix).c_str(),(" Higgs Y12 ("+suffix+")").c_str(), bin_btag, min_btag, max_btag*1.2 );
# Line 315 | Line 323 | public:
323      StH_nPVs->Fill(iEvent.nPVs, w);
324      StH_addJets->Fill(iEvent.CountAddJets(), w);
325  
326 <    StH_met_et->Fill( TMath::Min(iEvent.METnoPUCh_et, iEvent.METnoPU_et ), w);
326 >    StH_typeImet_et->Fill(iEvent.typeIcorrMET(), w);
327 >    //    StH_met_et->Fill( TMath::Min(iEvent.METnoPUCh_et, iEvent.METnoPU_et ), w);
328 >    StH_met_et->Fill( iEvent.MET_et, w);
329      StH_met_sig->Fill(iEvent.MET_sig, w);
330  
331      StH_vectorLepton1_pt->Fill(iEvent.vLepton_pt[0], w);
# Line 349 | Line 359 | public:
359        StH_simpleJet2_vtxMass->Fill(iEvent.hJet_vtxMass[1], w);
360      }
361      //here I fill both jets in one histo
362 +    StH_simpleJets_bTagSum->Fill(iEvent.hJet_csv[0]+iEvent.hJet_csv[1], w);
363      StH_simpleJets_bTag->Fill(iEvent.hJet_csv[0], w);
364      StH_simpleJets_bTag->Fill(iEvent.hJet_csv[1], w);
365      StH_simpleJets_vtxMass->Fill(iEvent.hJet_vtxMass[0], w);
# Line 394 | Line 405 | public:
405      }
406      if(iEvent.CountAddJets() > 0 ){
407        StH_HdR_addJets->Fill(higgs.DeltaR(addJets.at(dRmin_idx)), w);
408 <      if( iEvent.aJet_csv[dRmin_idx] > 0.7 )
408 >      StH_HdR_leading_addJets->Fill(hJet1.DeltaR(addJets.at(dRmin_idx)), w);
409 >      StH_HdR_following_addJets->Fill(hJet2.DeltaR(addJets.at(dRmin_idx)), w);
410 >      StH_HdRnorm_addJets->Fill(higgs.DeltaR(addJets.at(dRmin_idx)) / iEvent.H_dR, w);
411 >      if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  hJet1.DeltaR(addJets.at(dRmin_idx)) < 2 )
412 >        StH_HMass_leading_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
413 >      else StH_HMass_leading_addJets->Fill(higgs.M(), w);
414 >      if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  hJet2.DeltaR(addJets.at(dRmin_idx)) < 2 )
415 >        StH_HMass_following_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
416 >      else StH_HMass_following_addJets->Fill(higgs.M(), w);
417 >      if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  hJet2.DeltaR(addJets.at(dRmin_idx)) < 2  ||  hJet1.DeltaR(addJets.at(dRmin_idx)) < 2 )
418 >        StH_HMass_both_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
419 >      else StH_HMass_both_addJets->Fill(higgs.M(), w);
420 >      if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  higgs.DeltaR(addJets.at(dRmin_idx)) < 2 )
421          StH_HMass_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
422        else StH_HMass_addJets->Fill(higgs.M(), w);
423        StH_addJet1_pt->Fill(iEvent.aJet_pt[dRmin_idx], w);
# Line 427 | Line 450 | public:
450  
451    TH1F * StH_nPVs;
452    TH1F * StH_addJets;
453 +  TH1F * StH_typeImet_et;
454    TH1F * StH_met_et;
455    TH1F * StH_met_sig;
456  
# Line 455 | Line 479 | public:
479    TH1F * StH_simpleJet2_vtxMass;
480  
481    TH1F * StH_simpleJets_bTag;
482 +  TH1F * StH_simpleJets_bTagSum;
483    TH1F * StH_simpleJets_vtxMass;
484  
485    TH1F * StH_addJet1_pt;
# Line 472 | Line 497 | public:
497    TH1F * StH_HHel;
498    TH1F * StH_HPullAngle;
499    TH1F * StH_HdR_addJets;
500 +  TH1F * StH_HdR_leading_addJets;
501 +  TH1F * StH_HdR_following_addJets;
502 +
503 +  TH1F * StH_HdRnorm_addJets;
504    TH1F * StH_HMass_addJets;
505 +  TH1F * StH_HMass_leading_addJets;
506 +  TH1F * StH_HMass_following_addJets;
507 +  TH1F * StH_HMass_both_addJets;
508  
509    TH1F * StH_massDrop;
510    TH1F * StH_y12;
# Line 614 | Line 646 | public:
646      SystUPH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystUP").c_str(),(" ZH delta Phi ("+suffix+"SystUP"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
647      SystUPH_HdR_addJets = new TH1F(("HdR_addJets"+suffix+"SystUP").c_str(),("Higgs deltaR addJets ("+suffix+"SystUP"+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
648  
649 +    SystJERUPH_simpleJet1_pt = new TH1F(("SimpleJet1_pt"+suffix+"SystJERUP").c_str(),("Simple Jet1 pt ("+suffix+"SystJERUP"+")").c_str(), bin_pt, min_pt, max_pt*0.8 );
650 +    SystJERUPH_simpleJet2_pt = new TH1F(("SimpleJet2_pt"+suffix+"SystJERUP").c_str(),("Simple Jet2 pt ("+suffix+"SystJERUP"+")").c_str(), bin_pt, min_pt, max_pt*0.8 );
651 +    SystJERUPH_HMass = new TH1F(("HiggsMass"+suffix+"SystJERUP").c_str(),(" Higgs Mass ("+suffix+"SystJERUP"+")").c_str(), bin_mass, min_mass, max_mass );
652 +    SystJERUPH_HPt = new TH1F(("HiggsPt"+suffix+"SystJERUP").c_str(),(" Higgs Pt ("+suffix+"SystJERUP"+")").c_str(), bin_pt, min_pt, max_pt );
653 +    SystJERUPH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystJERUP").c_str(),(" ZH delta Phi ("+suffix+"SystJERUP"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
654 +
655      SystUPH_simpleJet1_bTag = new TH1F(("SimpleJet1_bTag"+suffix+"SystUP").c_str(),("Simple Jet1 bTag ("+suffix+"SystUP"+")").c_str(), bin_btag, min_btag, max_btag );
656      SystUPH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix+"SystUP").c_str(),("Simple Jet2 bTag ("+suffix+"SystUP"+")").c_str(), bin_btag, min_btag, max_btag );
657      SystUPH_simpleJets_bTag = new TH1F(("SimpleJets_bTag"+suffix+"SystUP").c_str(),("Simple Jets bTag ("+suffix+"SystUP"+")").c_str(), bin_btag, min_btag, max_btag );
# Line 627 | Line 665 | public:
665      SystDOWNH_HMass = new TH1F(("HiggsMass"+suffix+"SystDOWN").c_str(),(" Higgs Mass ("+suffix+"SystDOWN"+")").c_str(), bin_mass, min_mass, max_mass );
666      SystDOWNH_HPt = new TH1F(("HiggsPt"+suffix+"SystDOWN").c_str(),(" Higgs Pt ("+suffix+"SystDOWN"+")").c_str(), bin_pt, min_pt, max_pt );
667      SystDOWNH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystDOWN").c_str(),(" ZH delta Phi ("+suffix+"SystDOWN"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
668 <    SystDOWNH_HdR_addJets = new TH1F(("HdR_addJets"+suffix+"SystDOWN").c_str(),("Higgs deltaR addJets ("+suffix+"SystDOWN"+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
668 >    SystDOWNH_HdR_addJets = new TH1F(("HdR_addJets"+suffix+"SystDWON").c_str(),("Higgs deltaR addJets ("+suffix+"SystDOWN"+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
669 >
670 >    SystJERDOWNH_simpleJet1_pt = new TH1F(("SimpleJet1_pt"+suffix+"SystJERDOWN").c_str(),("Simple Jet1 pt ("+suffix+"SystJERDOWN"+")").c_str(), bin_pt, min_pt, max_pt*0.8 );
671 >    SystJERDOWNH_simpleJet2_pt = new TH1F(("SimpleJet2_pt"+suffix+"SystJERDOWN").c_str(),("Simple Jet2 pt ("+suffix+"SystJERDOWN"+")").c_str(), bin_pt, min_pt, max_pt*0.8 );
672 >    SystJERDOWNH_HMass = new TH1F(("HiggsMass"+suffix+"SystJERDOWN").c_str(),(" Higgs Mass ("+suffix+"SystJERDOWN"+")").c_str(), bin_mass, min_mass, max_mass );
673 >    SystJERDOWNH_HPt = new TH1F(("HiggsPt"+suffix+"SystJERDOWN").c_str(),(" Higgs Pt ("+suffix+"SystJERDOWN"+")").c_str(), bin_pt, min_pt, max_pt );
674 >    SystJERDOWNH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystJERDOWN").c_str(),(" ZH delta Phi ("+suffix+"SystJERDOWN"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
675  
676      SystDOWNH_simpleJet1_bTag = new TH1F(("SimpleJet1_bTag"+suffix+"SystDOWN").c_str(),("Simple Jet1 bTag ("+suffix+"SystDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
677      SystDOWNH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix+"SystDOWN").c_str(),("Simple Jet2 bTag ("+suffix+"SystDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
# Line 636 | Line 680 | public:
680      SystFDOWNH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix+"SystFDOWN").c_str(),("Simple Jet2 bTag ("+suffix+"SystFDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
681      SystFDOWNH_simpleJets_bTag = new TH1F(("SimpleJets_bTag"+suffix+"SystFDOWN").c_str(),("Simple Jets bTag ("+suffix+"SystFDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
682  
639 –
683    }
684    
685    virtual void fill(ntupleReader &iEvent,float w) {
686  
687      TLorentzVector higgs, addJet;
688        
646 –
689      //UPSCALE
690      SystUPH_addJets->Fill(iEvent.CountAddJets_jec(+1), w);
691      SystUPH_simpleJet1_pt->Fill(iEvent.hJet_pt_jecUP(0), w);
692      SystUPH_simpleJet2_pt->Fill(iEvent.hJet_pt_jecUP(1), w);
693 +    SystJERUPH_simpleJet1_pt->Fill(iEvent.hJet_jer(0,+1).Pt(), w);
694 +    SystJERUPH_simpleJet2_pt->Fill(iEvent.hJet_jer(1,+1).Pt(), w);
695 +
696      //btag sorted by btag
697      if( iEvent.hJet_csv_cUP(0) > iEvent.hJet_csv_cUP(1) ){
698        SystUPH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cUP(0), w);
# Line 675 | Line 720 | public:
720      SystUPH_HMass->Fill(iEvent.H_jecUP().M(), w);
721      SystUPH_HPt->Fill(iEvent.H_jecUP().Pt(), w);
722      SystUPH_ZH_dPhi->Fill( TMath::Abs(iEvent.H_jecUP().DeltaPhi( iEvent.VectorBoson() ) ), w);
723 +
724 +    SystJERUPH_HMass->Fill(iEvent.H_jer(+1).M(), w);
725 +    SystJERUPH_HPt->Fill(iEvent.H_jer(+1).Pt(), w);
726 +    SystJERUPH_ZH_dPhi->Fill( TMath::Abs(iEvent.H_jer(+1).DeltaPhi( iEvent.VectorBoson() ) ), w);
727 +
728      higgs=iEvent.H_jecUP();
729      addJet.SetPtEtaPhiE(iEvent.aJet_pt_jecUP(0),iEvent.aJet_eta[0],iEvent.aJet_phi[0],iEvent.aJet_e[0]);
730      SystUPH_HdR_addJets->Fill(higgs.DeltaR(addJet), w);
# Line 683 | Line 733 | public:
733      SystDOWNH_addJets->Fill(iEvent.CountAddJets_jec(-1), w);
734      SystDOWNH_simpleJet1_pt->Fill(iEvent.hJet_pt_jecDOWN(0), w);
735      SystDOWNH_simpleJet2_pt->Fill(iEvent.hJet_pt_jecDOWN(1), w);
736 +    SystJERDOWNH_simpleJet1_pt->Fill(iEvent.hJet_jer(0,-1).Pt(), w);
737 +    SystJERDOWNH_simpleJet2_pt->Fill(iEvent.hJet_jer(1,-1).Pt(), w);
738 +
739      //btag sorted by btag
740      if( iEvent.hJet_csv_cDOWN(0) > iEvent.hJet_csv_cDOWN(1) ){
741        SystDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cDOWN(0), w);
# Line 709 | Line 762 | public:
762      SystDOWNH_HMass->Fill(iEvent.H_jecDOWN().M(), w);
763      SystDOWNH_HPt->Fill(iEvent.H_jecDOWN().Pt(), w);
764      SystDOWNH_ZH_dPhi->Fill( TMath::Abs( iEvent.H_jecDOWN().DeltaPhi( iEvent.VectorBoson() ) ), w);
765 +
766 +    SystJERDOWNH_HMass->Fill(iEvent.H_jer(-1).M(), w);
767 +    SystJERDOWNH_HPt->Fill(iEvent.H_jer(-1).Pt(), w);
768 +    SystJERDOWNH_ZH_dPhi->Fill( TMath::Abs(iEvent.H_jer(-1).DeltaPhi( iEvent.VectorBoson() ) ), w);
769 +
770      higgs=iEvent.H_jecDOWN();
771      addJet.SetPtEtaPhiE(iEvent.aJet_pt_jecDOWN(0),iEvent.aJet_eta[0],iEvent.aJet_phi[0],iEvent.aJet_e[0]);
772      SystDOWNH_HdR_addJets->Fill(higgs.DeltaR(addJet), w);
773 <        
773 >
774    }
775  
776    TH1F * SystUPH_addJets;  
# Line 728 | Line 786 | public:
786    TH1F * SystFUPH_simpleJet1_bTag;
787    TH1F * SystFUPH_simpleJet2_bTag;
788    TH1F * SystFUPH_simpleJets_bTag;
789 +  TH1F * SystJERUPH_simpleJet1_pt;
790 +  TH1F * SystJERUPH_simpleJet2_pt;
791 +  TH1F * SystJERUPH_HMass;
792 +  TH1F * SystJERUPH_HPt;
793 +  TH1F * SystJERUPH_ZH_dPhi;
794  
795    TH1F * SystDOWNH_addJets;  
796    TH1F * SystDOWNH_simpleJet1_pt;
# Line 742 | Line 805 | public:
805    TH1F * SystFDOWNH_simpleJet1_bTag;
806    TH1F * SystFDOWNH_simpleJet2_bTag;
807    TH1F * SystFDOWNH_simpleJets_bTag;
808 +  TH1F * SystJERDOWNH_simpleJet1_pt;
809 +  TH1F * SystJERDOWNH_simpleJet2_pt;
810 +  TH1F * SystJERDOWNH_HMass;
811 +  TH1F * SystJERDOWNH_HPt;
812 +  TH1F * SystJERDOWNH_ZH_dPhi;
813  
814  
815   private:

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines