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.3 by bortigno, Tue May 8 08:13:27 2012 UTC

# Line 295 | Line 295 | public:
295      StH_HHel = new TH1F(("HiggsHel"+suffix).c_str(),("Higgs helicity angle ("+suffix+")").c_str(), bin_hel, min_hel, max_hel );
296      StH_HPullAngle = new TH1F(("HiggsPullAngle"+suffix).c_str(),("Higgs pull angle ("+suffix+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
297      StH_HdR_addJets = new TH1F(("HdR_addJets"+suffix).c_str(),("Higgs deltaR addJets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
298 +    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 );
299 +    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 );
300 +    StH_HdRnorm_addJets = new TH1F(("HdRnorm_addJets"+suffix).c_str(),("Higgs normalized deltaR addJets ("+suffix+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
301      StH_HMass_addJets = new TH1F(("HiggsMass_addJets"+suffix).c_str(),(" Higgs Mass with addJet ("+suffix+")").c_str(), bin_mass, min_mass, max_mass );
302 +    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 );
303 +    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 );
304 +    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 );
305  
306      StH_massDrop = new TH1F(("MassDrop"+suffix).c_str(),(" Higgs Mass Drop ("+suffix+")").c_str(), bin_btag, min_btag, max_btag );
307      StH_y12 = new TH1F(("Y12"+suffix).c_str(),(" Higgs Y12 ("+suffix+")").c_str(), bin_btag, min_btag, max_btag*1.2 );
# Line 394 | Line 400 | public:
400      }
401      if(iEvent.CountAddJets() > 0 ){
402        StH_HdR_addJets->Fill(higgs.DeltaR(addJets.at(dRmin_idx)), w);
403 <      if( iEvent.aJet_csv[dRmin_idx] > 0.7 )
403 >      StH_HdR_leading_addJets->Fill(hJet1.DeltaR(addJets.at(dRmin_idx)), w);
404 >      StH_HdR_following_addJets->Fill(hJet2.DeltaR(addJets.at(dRmin_idx)), w);
405 >      StH_HdRnorm_addJets->Fill(higgs.DeltaR(addJets.at(dRmin_idx)) / iEvent.H_dR, w);
406 >      if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  hJet1.DeltaR(addJets.at(dRmin_idx)) < 2 )
407 >        StH_HMass_leading_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
408 >      else StH_HMass_leading_addJets->Fill(higgs.M(), w);
409 >      if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  hJet2.DeltaR(addJets.at(dRmin_idx)) < 2 )
410 >        StH_HMass_following_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
411 >      else StH_HMass_following_addJets->Fill(higgs.M(), w);
412 >      if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  hJet2.DeltaR(addJets.at(dRmin_idx)) < 2  ||  hJet1.DeltaR(addJets.at(dRmin_idx)) < 2 )
413 >        StH_HMass_both_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
414 >      else StH_HMass_both_addJets->Fill(higgs.M(), w);
415 >      if( iEvent.aJet_csv[dRmin_idx] > 0.8 ||  higgs.DeltaR(addJets.at(dRmin_idx)) < 2 )
416          StH_HMass_addJets->Fill((higgs+addJets.at(dRmin_idx)).M(), w);
417        else StH_HMass_addJets->Fill(higgs.M(), w);
418        StH_addJet1_pt->Fill(iEvent.aJet_pt[dRmin_idx], w);
# Line 472 | Line 490 | public:
490    TH1F * StH_HHel;
491    TH1F * StH_HPullAngle;
492    TH1F * StH_HdR_addJets;
493 +  TH1F * StH_HdR_leading_addJets;
494 +  TH1F * StH_HdR_following_addJets;
495 +
496 +  TH1F * StH_HdRnorm_addJets;
497    TH1F * StH_HMass_addJets;
498 +  TH1F * StH_HMass_leading_addJets;
499 +  TH1F * StH_HMass_following_addJets;
500 +  TH1F * StH_HMass_both_addJets;
501  
502    TH1F * StH_massDrop;
503    TH1F * StH_y12;
# Line 614 | Line 639 | public:
639      SystUPH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystUP").c_str(),(" ZH delta Phi ("+suffix+"SystUP"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
640      SystUPH_HdR_addJets = new TH1F(("HdR_addJets"+suffix+"SystUP").c_str(),("Higgs deltaR addJets ("+suffix+"SystUP"+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
641  
642 +    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 );
643 +    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 );
644 +    SystJERUPH_HMass = new TH1F(("HiggsMass"+suffix+"SystJERUP").c_str(),(" Higgs Mass ("+suffix+"SystJERUP"+")").c_str(), bin_mass, min_mass, max_mass );
645 +    SystJERUPH_HPt = new TH1F(("HiggsPt"+suffix+"SystJERUP").c_str(),(" Higgs Pt ("+suffix+"SystJERUP"+")").c_str(), bin_pt, min_pt, max_pt );
646 +    SystJERUPH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystJERUP").c_str(),(" ZH delta Phi ("+suffix+"SystJERUP"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
647 +
648      SystUPH_simpleJet1_bTag = new TH1F(("SimpleJet1_bTag"+suffix+"SystUP").c_str(),("Simple Jet1 bTag ("+suffix+"SystUP"+")").c_str(), bin_btag, min_btag, max_btag );
649      SystUPH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix+"SystUP").c_str(),("Simple Jet2 bTag ("+suffix+"SystUP"+")").c_str(), bin_btag, min_btag, max_btag );
650      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 658 | public:
658      SystDOWNH_HMass = new TH1F(("HiggsMass"+suffix+"SystDOWN").c_str(),(" Higgs Mass ("+suffix+"SystDOWN"+")").c_str(), bin_mass, min_mass, max_mass );
659      SystDOWNH_HPt = new TH1F(("HiggsPt"+suffix+"SystDOWN").c_str(),(" Higgs Pt ("+suffix+"SystDOWN"+")").c_str(), bin_pt, min_pt, max_pt );
660      SystDOWNH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystDOWN").c_str(),(" ZH delta Phi ("+suffix+"SystDOWN"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
661 <    SystDOWNH_HdR_addJets = new TH1F(("HdR_addJets"+suffix+"SystDOWN").c_str(),("Higgs deltaR addJets ("+suffix+"SystDOWN"+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
661 >    SystDOWNH_HdR_addJets = new TH1F(("HdR_addJets"+suffix+"SystDWON").c_str(),("Higgs deltaR addJets ("+suffix+"SystDOWN"+")").c_str(), bin_deltaR, min_deltaR, max_deltaR );
662 >
663 >    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 );
664 >    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 );
665 >    SystJERDOWNH_HMass = new TH1F(("HiggsMass"+suffix+"SystJERDOWN").c_str(),(" Higgs Mass ("+suffix+"SystJERDOWN"+")").c_str(), bin_mass, min_mass, max_mass );
666 >    SystJERDOWNH_HPt = new TH1F(("HiggsPt"+suffix+"SystJERDOWN").c_str(),(" Higgs Pt ("+suffix+"SystJERDOWN"+")").c_str(), bin_pt, min_pt, max_pt );
667 >    SystJERDOWNH_ZH_dPhi = new TH1F(("ZH_dPhi"+suffix+"SystJERDOWN").c_str(),(" ZH delta Phi ("+suffix+"SystJERDOWN"+")").c_str(), bin_deltaPhi, min_deltaPhi, max_deltaPhi );
668  
669      SystDOWNH_simpleJet1_bTag = new TH1F(("SimpleJet1_bTag"+suffix+"SystDOWN").c_str(),("Simple Jet1 bTag ("+suffix+"SystDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
670      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 673 | public:
673      SystFDOWNH_simpleJet2_bTag = new TH1F(("SimpleJet2_bTag"+suffix+"SystFDOWN").c_str(),("Simple Jet2 bTag ("+suffix+"SystFDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
674      SystFDOWNH_simpleJets_bTag = new TH1F(("SimpleJets_bTag"+suffix+"SystFDOWN").c_str(),("Simple Jets bTag ("+suffix+"SystFDOWN"+")").c_str(), bin_btag, min_btag, max_btag );
675  
639 –
676    }
677    
678    virtual void fill(ntupleReader &iEvent,float w) {
679  
680      TLorentzVector higgs, addJet;
681        
646 –
682      //UPSCALE
683      SystUPH_addJets->Fill(iEvent.CountAddJets_jec(+1), w);
684      SystUPH_simpleJet1_pt->Fill(iEvent.hJet_pt_jecUP(0), w);
685      SystUPH_simpleJet2_pt->Fill(iEvent.hJet_pt_jecUP(1), w);
686 +    SystJERUPH_simpleJet1_pt->Fill(iEvent.hJet_jer(0,+1).Pt(), w);
687 +    SystJERUPH_simpleJet2_pt->Fill(iEvent.hJet_jer(1,+1).Pt(), w);
688 +
689      //btag sorted by btag
690      if( iEvent.hJet_csv_cUP(0) > iEvent.hJet_csv_cUP(1) ){
691        SystUPH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cUP(0), w);
# Line 675 | Line 713 | public:
713      SystUPH_HMass->Fill(iEvent.H_jecUP().M(), w);
714      SystUPH_HPt->Fill(iEvent.H_jecUP().Pt(), w);
715      SystUPH_ZH_dPhi->Fill( TMath::Abs(iEvent.H_jecUP().DeltaPhi( iEvent.VectorBoson() ) ), w);
716 +
717 +    SystJERUPH_HMass->Fill(iEvent.H_jer(+1).M(), w);
718 +    SystJERUPH_HPt->Fill(iEvent.H_jer(+1).Pt(), w);
719 +    SystJERUPH_ZH_dPhi->Fill( TMath::Abs(iEvent.H_jer(+1).DeltaPhi( iEvent.VectorBoson() ) ), w);
720 +
721      higgs=iEvent.H_jecUP();
722      addJet.SetPtEtaPhiE(iEvent.aJet_pt_jecUP(0),iEvent.aJet_eta[0],iEvent.aJet_phi[0],iEvent.aJet_e[0]);
723      SystUPH_HdR_addJets->Fill(higgs.DeltaR(addJet), w);
# Line 683 | Line 726 | public:
726      SystDOWNH_addJets->Fill(iEvent.CountAddJets_jec(-1), w);
727      SystDOWNH_simpleJet1_pt->Fill(iEvent.hJet_pt_jecDOWN(0), w);
728      SystDOWNH_simpleJet2_pt->Fill(iEvent.hJet_pt_jecDOWN(1), w);
729 +    SystJERDOWNH_simpleJet1_pt->Fill(iEvent.hJet_jer(0,-1).Pt(), w);
730 +    SystJERDOWNH_simpleJet2_pt->Fill(iEvent.hJet_jer(1,-1).Pt(), w);
731 +
732      //btag sorted by btag
733      if( iEvent.hJet_csv_cDOWN(0) > iEvent.hJet_csv_cDOWN(1) ){
734        SystDOWNH_simpleJet1_bTag->Fill(iEvent.hJet_csv_cDOWN(0), w);
# Line 709 | Line 755 | public:
755      SystDOWNH_HMass->Fill(iEvent.H_jecDOWN().M(), w);
756      SystDOWNH_HPt->Fill(iEvent.H_jecDOWN().Pt(), w);
757      SystDOWNH_ZH_dPhi->Fill( TMath::Abs( iEvent.H_jecDOWN().DeltaPhi( iEvent.VectorBoson() ) ), w);
758 +
759 +    SystJERDOWNH_HMass->Fill(iEvent.H_jer(-1).M(), w);
760 +    SystJERDOWNH_HPt->Fill(iEvent.H_jer(-1).Pt(), w);
761 +    SystJERDOWNH_ZH_dPhi->Fill( TMath::Abs(iEvent.H_jer(-1).DeltaPhi( iEvent.VectorBoson() ) ), w);
762 +
763      higgs=iEvent.H_jecDOWN();
764      addJet.SetPtEtaPhiE(iEvent.aJet_pt_jecDOWN(0),iEvent.aJet_eta[0],iEvent.aJet_phi[0],iEvent.aJet_e[0]);
765      SystDOWNH_HdR_addJets->Fill(higgs.DeltaR(addJet), w);
766 <        
766 >
767    }
768  
769    TH1F * SystUPH_addJets;  
# Line 728 | Line 779 | public:
779    TH1F * SystFUPH_simpleJet1_bTag;
780    TH1F * SystFUPH_simpleJet2_bTag;
781    TH1F * SystFUPH_simpleJets_bTag;
782 +  TH1F * SystJERUPH_simpleJet1_pt;
783 +  TH1F * SystJERUPH_simpleJet2_pt;
784 +  TH1F * SystJERUPH_HMass;
785 +  TH1F * SystJERUPH_HPt;
786 +  TH1F * SystJERUPH_ZH_dPhi;
787  
788    TH1F * SystDOWNH_addJets;  
789    TH1F * SystDOWNH_simpleJet1_pt;
# Line 742 | Line 798 | public:
798    TH1F * SystFDOWNH_simpleJet1_bTag;
799    TH1F * SystFDOWNH_simpleJet2_bTag;
800    TH1F * SystFDOWNH_simpleJets_bTag;
801 +  TH1F * SystJERDOWNH_simpleJet1_pt;
802 +  TH1F * SystJERDOWNH_simpleJet2_pt;
803 +  TH1F * SystJERDOWNH_HMass;
804 +  TH1F * SystJERDOWNH_HPt;
805 +  TH1F * SystJERDOWNH_ZH_dPhi;
806  
807  
808   private:

Diff Legend

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