ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbb/interface/ntupleReader.hpp
(Generate patch)

Comparing UserCode/VHbb/interface/ntupleReader.hpp (file contents):
Revision 1.2 by bortigno, Wed May 2 19:20:07 2012 UTC vs.
Revision 1.3 by bortigno, Tue May 8 08:13:10 2012 UTC

# Line 21 | Line 21 | public :
21     virtual int CountAddLeptons();
22     virtual bool TriggerBit();
23     virtual TLorentzVector VectorBoson(); //vector boson TLorentz vector
24 +   virtual double resolution(double eta); //smearing for the jet energy resolution
25 +   virtual TLorentzVector hJet_jer( int idx, double sign );
26 +   virtual TLorentzVector H_jer( double sing );
27     virtual double hJet_PT( int idx, int sign ); //higgs jet energy correction
28     virtual double aJet_PT( int idx, int sign ); //addtional jet energy correction
29     virtual TLorentzVector Higgs( int sign ); //higgs candidate jet energy correction
# Line 69 | Line 72 | TLorentzVector ntupleReader::VectorBoson
72    return (l1+l2);
73   }
74  
75 + //JER
76 + double ntupleReader::resolution(double eta){
77 +  double inner = 0.06;
78 +  double outer = 0.1;
79 +  double eta_tracker = 1.1;
80 +  if(abs(eta) < eta_tracker) return inner; else return outer;
81 + }
82 + TLorentzVector ntupleReader::hJet_jer( int idx, double sign ){
83 +  TLorentzVector tmp;
84 +  double hJet_pt_jer = hJet_pt[idx] + sign * resolution(hJet_eta[idx])*TMath::Abs(hJet_pt[idx]-hJet_genPt[idx]);
85 +  tmp.SetPtEtaPhiE(  hJet_pt_jer,
86 +                     hJet_eta[idx],
87 +                     hJet_phi[idx],
88 +                     hJet_e[idx]*(hJet_pt_jer/hJet_pt[idx]) );
89 + return tmp;
90 + }
91 +
92 + TLorentzVector ntupleReader::H_jer( double sign ){
93 +  TLorentzVector h;
94 +  h = hJet_jer( 0, sign ) + hJet_jer( 1, sign );
95 +  return h;
96 + }
97 +
98   double ntupleReader::hJet_PT( int idx, int sign ){ return  hJet_pt[idx]*(1 + (sign)*hJet_JECUnc[idx]); }
99   double ntupleReader::aJet_PT( int idx, int sign ){ return  aJet_pt[idx]*(1 + (sign)*aJet_JECUnc[idx]); }
100   TLorentzVector ntupleReader::Higgs( int sign ){

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines