VHbbDataFormats/interface/VHbbCandidateTools.h

#ifndef VHBBCANDIDATETOOLS_H
#define VHBBCANDIDATETOOLS_H

#include "VHbbAnalysis/VHbbDataFormats/interface/VHbbCandidate.h"

#include <iostream>

struct CompareJetPtMuons {
  bool operator()( const VHbbEvent::MuonInfo& j1, const  VHbbEvent::MuonInfo& j2 ) const {
    return j1.p4.Pt() > j2.p4.Pt();
  }
};
struct CompareJetPtElectrons {
  bool operator()( const VHbbEvent::ElectronInfo& j1, const  VHbbEvent::ElectronInfo& j2 ) const {
    return j1.p4.Pt() > j2.p4.Pt();
  }
};


class VHbbCandidateTools {
 public:

 VHbbCandidateTools(bool verbose = false): verbose_(verbose){}

  float deltaPhi(float in2, float in1){
    float dphi = in2-in1;
    if ( dphi > M_PI ) {
      dphi -= 2.0*M_PI;
    } else if ( dphi <= -M_PI ) {
      dphi += 2.0*M_PI;
    }
    return dphi;
  }
  
  
  VHbbCandidate getHZmumuCandidate(const VHbbCandidate & in, bool & ok, std::vector<unsigned int>& pos){
    if (verbose_){
      std::cout <<" getHZmumuCandidate input mu "<<in.V.muons.size()<<" e "<<in.V.electrons.size()<<std::endl;
    }
    ok = false;
    VHbbCandidate temp=in;
    
    //
    // change: allow for additional leptons; by definition 
    //
    if (temp.V.muons.size()<2) return in ;
    //    if (temp.V.electrons.size()!=0) return in ;
    std::vector<VHbbEvent::MuonInfo> muons_ = temp.V.muons;

    // beware: assumes already sorted!!!!

    //    CompareJetPtMuons ptComparator;
    //    std::sort(muons_.begin(), muons_.end(), ptComparator);
    if (muons_[0].p4.Pt()<20 || muons_[1].p4.Pt()<20 ) return in;
    
    //
    // now I need to ask also for the charge
    //
    int selectMu2=1;
    
    if (muons_[0].charge* muons_[selectMu2].charge > 0){
      if (muons_.size() ==2) return in;
      //
      // i need to find a proper pair
      //
      
      for (unsigned int it=2; it< muons_.size(); ++it){
        if (  muons_[it].charge * muons_[0].charge < 0) {
          selectMu2 = it;
          break;
        }
        if (selectMu2 == 1) return in;
      }  
    }
    temp.V.p4 = muons_[0].p4+muons_[selectMu2].p4;
    std::vector<VHbbEvent::MuonInfo> muons2_;
    for (std::vector<VHbbEvent::MuonInfo>::const_iterator it = muons_.begin(); it!= muons_.end(); ++it){
      if (it->p4.Pt()>20) muons2_.push_back(*it);
    }
    temp.V.muons = muons2_;
    
    // the same for electrons
    std::vector<VHbbEvent::ElectronInfo> electrons_ = temp.V.electrons;

    // beware; assumes already sorted

    //    CompareJetPtElectrons ptComparator2;
    //    std::sort(electrons_.begin(), electrons_.end(), ptComparator2);
    std::vector<VHbbEvent::ElectronInfo> electrons2_;
    for (std::vector<VHbbEvent::ElectronInfo>::const_iterator it = electrons_.begin(); it!= electrons_.end(); ++it){
      if (it->p4.Pt()>20) electrons2_.push_back(*it);
    }
    temp.V.electrons = electrons2_;
    
    //
    // consider all 
    //
    
    
    //    if (temp.V.Pt()<150 ) return in;
    //    if (temp.H.Pt()<150) return in;
    //    if (temp.H.firstJet().csv< 0.9) return in;
    //    if (temp.H.secondJet().csv<0.5) return in;
    //    if (deltaPhi(temp.V.Phi(),temp.H.Phi())<2.7) return in;
    //    if (temp.V.FourMomentum.Mass()<75 || temp.V.FourMomentum.Mass()>105) return in; 
    //    if (temp.additionalJets.size()>0) return in;
    //    if (std::Abs(deltaTheta) ????

    temp.V.firstLepton = pos[0];
    temp.V.secondLepton = pos[selectMu2];
   ok = true;
    return temp;
  }  
  VHbbCandidate getHZeeCandidate(const VHbbCandidate & in, bool & ok, std::vector<unsigned int>& pos){
    if (verbose_){
      std::cout <<" getHZeeCandidate input mu "<<in.V.electrons.size()<<" e "<<in.V.muons.size()<<std::endl;
    }
    ok = false;
    VHbbCandidate temp=in;
    
    //
    // change: allow for additional leptons; by definition 
    //
    if (temp.V.electrons.size()<2) return in ;
    //    if (temp.V.electrons.size()!=0) return in ;
    std::vector<VHbbEvent::ElectronInfo> electrons_ = temp.V.electrons;

    // beware assumes already sorted

    //    CompareJetPtElectrons ptComparator;
    //    std::sort(electrons_.begin(), electrons_.end(), ptComparator);
    if (electrons_[0].p4.Pt()<20 || electrons_[1].p4.Pt()<20 ) return in;
    //
    // now I need to ask also for the charge
    //
    int selectE2=1;
    
    if (electrons_[0].charge* electrons_[selectE2].charge > 0){
      if (electrons_.size() ==2) return in;
      //
      // i need to find a proper pair
      //
      
      for (unsigned int it=2; it< electrons_.size(); ++it){
        if (  electrons_[it].charge * electrons_[0].charge < 0) {
          selectE2 = it;
          break;
        }
        if (selectE2 == 1) return in;
      }  
    }
   temp.V.p4 = electrons_[0].p4+electrons_[selectE2].p4;
    std::vector<VHbbEvent::ElectronInfo> electrons2_;
    for (std::vector<VHbbEvent::ElectronInfo>::const_iterator it = electrons_.begin(); it!= electrons_.end(); ++it){
      if (it->p4.Pt()>20) electrons2_.push_back(*it);
    }
    temp.V.electrons = electrons2_;
  
    // the same for muonss
    std::vector<VHbbEvent::MuonInfo> muons_ = temp.V.muons;

    // beware assumes already sorted

    //    CompareJetPtMuons ptComparator2;
    //    std::sort(muons_.begin(), muons_.end(), ptComparator2);
    std::vector<VHbbEvent::MuonInfo> muons2_;
    for (std::vector<VHbbEvent::MuonInfo>::const_iterator it = muons_.begin(); it!= muons_.end(); ++it){
      if (it->p4.Pt()>20) muons2_.push_back(*it);
    }
    temp.V.muons = muons2_;

    temp.V.firstLepton = pos[0];
    temp.V.secondLepton = pos[selectE2];

    
   ok = true;
    return temp;
  }  
  VHbbCandidate getHZnnCandidate(const VHbbCandidate & in, bool & ok){
    if (verbose_){
      std::cout <<" getHZnnCandidate input mu "<<in.V.muons.size()<<" e "<<in.V.electrons.size()<<" met "<<in.V.mets.size()<<std::endl;
    }
    
    ok = false;
    VHbbCandidate temp=in;
    if (temp.V.mets.size()!=1) return in;
//always build a NuNu candidate, exclusion come from if/else series on the caller side
// this allow to still build candidates for mu+e  and same sign dilept
//    if (temp.V.muons.size()!=0) return in ;
//    if (temp.V.electrons.size()!=0) return in ;
    
    temp.V.p4 = temp.V.mets[0].p4;
    if (verbose_) {
      std::cout <<" debug met "<< temp.V.mets[0].metSig << " " <<  temp.V.mets[0].sumEt<< std::endl;
    }   
    if (temp.V.mets[0].metSig<5) return in;
    if (temp.V.mets[0].p4.Pt()<80) return in;
    //    if (temp.H.p4.Pt()<150)return in;
    //    if (temp.H.firstJet().csv< 0.9) return in;
    //    if (temp.H.secondJet().csv<0.5) return in;
    //    if (deltaPhi(temp.V.p4.Phi(),temp.H.p4.Phi())<2.95) return in;
    //    if (temp.V.electrons.size()>0 || temp.V.muons.size()>0 ) return in;
    //    if (temp.additionalJets.size()>0) return in;
    //    if (std::Abs(deltaTheta) ????
    
    ok = true;
    return temp;
  }

  VHbbCandidate getHWmunCandidate(const VHbbCandidate & in, bool & ok , std::vector<unsigned int>& pos){
    ok = false;
    VHbbCandidate temp=in;
    // require a muon and no electrons
    if (temp.V.muons.size()!=1) return in ;
    if (temp.V.electrons.size()!=0) return in ;
    if (temp.V.mets.size()<1) return in ;
    //
    /*pT(W) > 150 GeV (pt(W) computed using lepton px,y and PF MET x and y components)
      pT(H) > 150 GeV
      best btag, CSV > 0.90
      second-best btag, CSV > 0.50
      Dphi(W,H) > 2.95
      no additional isolated leptons (pT > 15 GeV)
      
      same lepton definition as in the physics objects section 
      
      No additional ak5PFjets (pT > 30 GeV; |eta| < 2.4)
      MET>35. for the electron BDT analysis
      |cos(theta * )| (TBO)
      color flow pull angle (TBO)
      We don't cut on the transverse mass (for boosted objects cutting on the transverse mass introduces an inefficiency due to the angle between the MET and the lepton being close to 0.) 
    */
    
    temp.V.p4 = temp.V.muons[0].p4+temp.V.mets[0].p4;
    temp.V.firstLepton=pos[0];

    ok=true;
    return temp;
  }

  VHbbCandidate getHWenCandidate(const VHbbCandidate & in, bool & ok, std::vector<unsigned int>& pos){
    ok = false;
    VHbbCandidate temp=in;
    if (temp.V.electrons.size()!=1) return in ;
    if (temp.V.muons.size()!=0) return in ;
    if (temp.V.mets.size()<1) return in ;
    //
    /*pT(W) > 150 GeV (pt(W) computed using lepton px,y and PF MET x and y components)
      pT(H) > 150 GeV
      best btag, CSV > 0.90
      second-best btag, CSV > 0.50
      Dphi(W,H) > 2.95
      no additional isolated leptons (pT > 15 GeV)
      
      same lepton definition as in the physics objects section 
      
      No additional ak5PFjets (pT > 30 GeV; |eta| < 2.4)
      MET>35. for the electron BDT analysis
      |cos(theta * )| (TBO)
      color flow pull angle (TBO)
      We don't cut on the transverse mass (for boosted objects cutting on the transverse mass introduces an inefficiency due to the angle between the MET and the lepton being close to 0.) 
    */
    
    ok=true;
    temp.V.firstLepton=pos[0];
    temp.V.p4 = temp.V.electrons[0].p4+temp.V.mets[0].p4;
    return temp;
  }

 public: 
  bool verbose_;
  
 public:
  static float getDeltaTheta( const VHbbEvent::SimpleJet & j1, const VHbbEvent::SimpleJet & j2 ) {
    double deltaTheta = 1e10;
    TLorentzVector pi(0,0,0,0);
    TLorentzVector v_j1 = j1.chargedTracksFourMomentum;
    TLorentzVector v_j2 = j2.chargedTracksFourMomentum;
    
    if( v_j2.Mag() == 0 
        || v_j1.Mag() == 0 )
      return deltaTheta = 1e10;
    
    //use j1 to calculate the pull vector
    TVector2 t = j1.tVector;
    
    if( t.Mod() == 0 )
      return deltaTheta = 1e10;
    
    Double_t dphi =  v_j2.Phi()- v_j1.Phi();
    if ( dphi > M_PI ) {
      dphi -= 2.0*M_PI;
    } else if ( dphi <= -M_PI ) {
      dphi += 2.0*M_PI;
    }
    Double_t deltaeta = v_j2.Rapidity() - v_j1.Rapidity();
    TVector2 BBdir( deltaeta, dphi );
    
    deltaTheta = t.DeltaPhi(BBdir);
    
    return deltaTheta;
    
  }
  
  float getHelicity(const VHbbEvent::SimpleJet& j, TVector3 boost) const {
    double hel = 1e10;
    TLorentzVector jet = j.p4;
    jet.Boost( -boost );
    hel = TMath::Cos( jet.Vect().Angle( boost ) );
    return hel;
  }


};


#endif


Revision:	1.11
Committed:	Wed Sep 14 07:46:59 2011 UTC (13 years, 7 months ago) by tboccali
Content type:	text/plain
Branch:	MAIN
CVS Tags:	EdmV9Sept2011, Sept19th2011_2, Sept19th2011, Sept19th, VHNtupleV9_AR1, VHSept15_AR1, Sept14th2011_2, Sept14th2011
Changes since 1.10:	+1 -1 lines
Log Message:	fix et cut for znunu and exclusivity of candidates
#	Content
1	#ifndef VHBBCANDIDATETOOLS_H
2	#define VHBBCANDIDATETOOLS_H
3
4	#include "VHbbAnalysis/VHbbDataFormats/interface/VHbbCandidate.h"
5
6	#include <iostream>
7
8	struct CompareJetPtMuons {
9	bool operator()( const VHbbEvent::MuonInfo& j1, const VHbbEvent::MuonInfo& j2 ) const {
10	return j1.p4.Pt() > j2.p4.Pt();
11	}
12	};
13	struct CompareJetPtElectrons {
14	bool operator()( const VHbbEvent::ElectronInfo& j1, const VHbbEvent::ElectronInfo& j2 ) const {
15	return j1.p4.Pt() > j2.p4.Pt();
16	}
17	};
18
19
20
21	class VHbbCandidateTools {
22	public:
23
24	VHbbCandidateTools(bool verbose = false): verbose_(verbose){}
25
26	float deltaPhi(float in2, float in1){
27	float dphi = in2-in1;
28	if ( dphi > M_PI ) {
29	dphi -= 2.0*M_PI;
30	} else if ( dphi <= -M_PI ) {
31	dphi += 2.0*M_PI;
32	}
33	return dphi;
34	}
35
36
37
38
39	VHbbCandidate getHZmumuCandidate(const VHbbCandidate & in, bool & ok, std::vector<unsigned int>& pos){
40	if (verbose_){
41	std::cout <<" getHZmumuCandidate input mu "<<in.V.muons.size()<<" e "<<in.V.electrons.size()<<std::endl;
42	}
43	ok = false;
44	VHbbCandidate temp=in;
45
46	//
47	// change: allow for additional leptons; by definition
48	//
49	if (temp.V.muons.size()<2) return in ;
50	// if (temp.V.electrons.size()!=0) return in ;
51	std::vector<VHbbEvent::MuonInfo> muons_ = temp.V.muons;
52
53	// beware: assumes already sorted!!!!
54
55	// CompareJetPtMuons ptComparator;
56	// std::sort(muons_.begin(), muons_.end(), ptComparator);
57	if (muons_[0].p4.Pt()<20 \|\| muons_[1].p4.Pt()<20 ) return in;
58
59	//
60	// now I need to ask also for the charge
61	//
62	int selectMu2=1;
63
64	if (muons_[0].charge* muons_[selectMu2].charge > 0){
65	if (muons_.size() ==2) return in;
66	//
67	// i need to find a proper pair
68	//
69
70	for (unsigned int it=2; it< muons_.size(); ++it){
71	if ( muons_[it].charge * muons_[0].charge < 0) {
72	selectMu2 = it;
73	break;
74	}
75	if (selectMu2 == 1) return in;
76	}
77	}
78	temp.V.p4 = muons_[0].p4+muons_[selectMu2].p4;
79	std::vector<VHbbEvent::MuonInfo> muons2_;
80	for (std::vector<VHbbEvent::MuonInfo>::const_iterator it = muons_.begin(); it!= muons_.end(); ++it){
81	if (it->p4.Pt()>20) muons2_.push_back(*it);
82	}
83	temp.V.muons = muons2_;
84
85	// the same for electrons
86	std::vector<VHbbEvent::ElectronInfo> electrons_ = temp.V.electrons;
87
88	// beware; assumes already sorted
89
90	// CompareJetPtElectrons ptComparator2;
91	// std::sort(electrons_.begin(), electrons_.end(), ptComparator2);
92	std::vector<VHbbEvent::ElectronInfo> electrons2_;
93	for (std::vector<VHbbEvent::ElectronInfo>::const_iterator it = electrons_.begin(); it!= electrons_.end(); ++it){
94	if (it->p4.Pt()>20) electrons2_.push_back(*it);
95	}
96	temp.V.electrons = electrons2_;
97
98	//
99	// consider all
100	//
101
102
103	// if (temp.V.Pt()<150 ) return in;
104	// if (temp.H.Pt()<150) return in;
105	// if (temp.H.firstJet().csv< 0.9) return in;
106	// if (temp.H.secondJet().csv<0.5) return in;
107	// if (deltaPhi(temp.V.Phi(),temp.H.Phi())<2.7) return in;
108	// if (temp.V.FourMomentum.Mass()<75 \|\| temp.V.FourMomentum.Mass()>105) return in;
109	// if (temp.additionalJets.size()>0) return in;
110	// if (std::Abs(deltaTheta) ????
111
112	temp.V.firstLepton = pos[0];
113	temp.V.secondLepton = pos[selectMu2];
114	ok = true;
115	return temp;
116	}
117	VHbbCandidate getHZeeCandidate(const VHbbCandidate & in, bool & ok, std::vector<unsigned int>& pos){
118	if (verbose_){
119	std::cout <<" getHZeeCandidate input mu "<<in.V.electrons.size()<<" e "<<in.V.muons.size()<<std::endl;
120	}
121	ok = false;
122	VHbbCandidate temp=in;
123
124	//
125	// change: allow for additional leptons; by definition
126	//
127	if (temp.V.electrons.size()<2) return in ;
128	// if (temp.V.electrons.size()!=0) return in ;
129	std::vector<VHbbEvent::ElectronInfo> electrons_ = temp.V.electrons;
130
131	// beware assumes already sorted
132
133	// CompareJetPtElectrons ptComparator;
134	// std::sort(electrons_.begin(), electrons_.end(), ptComparator);
135	if (electrons_[0].p4.Pt()<20 \|\| electrons_[1].p4.Pt()<20 ) return in;
136	//
137	// now I need to ask also for the charge
138	//
139	int selectE2=1;
140
141	if (electrons_[0].charge* electrons_[selectE2].charge > 0){
142	if (electrons_.size() ==2) return in;
143	//
144	// i need to find a proper pair
145	//
146
147	for (unsigned int it=2; it< electrons_.size(); ++it){
148	if ( electrons_[it].charge * electrons_[0].charge < 0) {
149	selectE2 = it;
150	break;
151	}
152	if (selectE2 == 1) return in;
153	}
154	}
155	temp.V.p4 = electrons_[0].p4+electrons_[selectE2].p4;
156	std::vector<VHbbEvent::ElectronInfo> electrons2_;
157	for (std::vector<VHbbEvent::ElectronInfo>::const_iterator it = electrons_.begin(); it!= electrons_.end(); ++it){
158	if (it->p4.Pt()>20) electrons2_.push_back(*it);
159	}
160	temp.V.electrons = electrons2_;
161
162	// the same for muonss
163	std::vector<VHbbEvent::MuonInfo> muons_ = temp.V.muons;
164
165	// beware assumes already sorted
166
167	// CompareJetPtMuons ptComparator2;
168	// std::sort(muons_.begin(), muons_.end(), ptComparator2);
169	std::vector<VHbbEvent::MuonInfo> muons2_;
170	for (std::vector<VHbbEvent::MuonInfo>::const_iterator it = muons_.begin(); it!= muons_.end(); ++it){
171	if (it->p4.Pt()>20) muons2_.push_back(*it);
172	}
173	temp.V.muons = muons2_;
174
175	temp.V.firstLepton = pos[0];
176	temp.V.secondLepton = pos[selectE2];
177
178
179	ok = true;
180	return temp;
181	}
182	VHbbCandidate getHZnnCandidate(const VHbbCandidate & in, bool & ok){
183	if (verbose_){
184	std::cout <<" getHZnnCandidate input mu "<<in.V.muons.size()<<" e "<<in.V.electrons.size()<<" met "<<in.V.mets.size()<<std::endl;
185	}
186
187	ok = false;
188	VHbbCandidate temp=in;
189	if (temp.V.mets.size()!=1) return in;
190	//always build a NuNu candidate, exclusion come from if/else series on the caller side
191	// this allow to still build candidates for mu+e and same sign dilept
192	// if (temp.V.muons.size()!=0) return in ;
193	// if (temp.V.electrons.size()!=0) return in ;
194
195	temp.V.p4 = temp.V.mets[0].p4;
196	if (verbose_) {
197	std::cout <<" debug met "<< temp.V.mets[0].metSig << " " << temp.V.mets[0].sumEt<< std::endl;
198	}
199	if (temp.V.mets[0].metSig<5) return in;
200	if (temp.V.mets[0].p4.Pt()<80) return in;
201	// if (temp.H.p4.Pt()<150)return in;
202	// if (temp.H.firstJet().csv< 0.9) return in;
203	// if (temp.H.secondJet().csv<0.5) return in;
204	// if (deltaPhi(temp.V.p4.Phi(),temp.H.p4.Phi())<2.95) return in;
205	// if (temp.V.electrons.size()>0 \|\| temp.V.muons.size()>0 ) return in;
206	// if (temp.additionalJets.size()>0) return in;
207	// if (std::Abs(deltaTheta) ????
208
209	ok = true;
210	return temp;
211	}
212
213	VHbbCandidate getHWmunCandidate(const VHbbCandidate & in, bool & ok , std::vector<unsigned int>& pos){
214	ok = false;
215	VHbbCandidate temp=in;
216	// require a muon and no electrons
217	if (temp.V.muons.size()!=1) return in ;
218	if (temp.V.electrons.size()!=0) return in ;
219	if (temp.V.mets.size()<1) return in ;
220	//
221	/*pT(W) > 150 GeV (pt(W) computed using lepton px,y and PF MET x and y components)
222	pT(H) > 150 GeV
223	best btag, CSV > 0.90
224	second-best btag, CSV > 0.50
225	Dphi(W,H) > 2.95
226	no additional isolated leptons (pT > 15 GeV)
227
228	same lepton definition as in the physics objects section
229
230	No additional ak5PFjets (pT > 30 GeV; \|eta\| < 2.4)
231	MET>35. for the electron BDT analysis
232	\|cos(theta * )\| (TBO)
233	color flow pull angle (TBO)
234	We don't cut on the transverse mass (for boosted objects cutting on the transverse mass introduces an inefficiency due to the angle between the MET and the lepton being close to 0.)
235	*/
236
237	temp.V.p4 = temp.V.muons[0].p4+temp.V.mets[0].p4;
238	temp.V.firstLepton=pos[0];
239
240	ok=true;
241	return temp;
242	}
243
244	VHbbCandidate getHWenCandidate(const VHbbCandidate & in, bool & ok, std::vector<unsigned int>& pos){
245	ok = false;
246	VHbbCandidate temp=in;
247	if (temp.V.electrons.size()!=1) return in ;
248	if (temp.V.muons.size()!=0) return in ;
249	if (temp.V.mets.size()<1) return in ;
250	//
251	/*pT(W) > 150 GeV (pt(W) computed using lepton px,y and PF MET x and y components)
252	pT(H) > 150 GeV
253	best btag, CSV > 0.90
254	second-best btag, CSV > 0.50
255	Dphi(W,H) > 2.95
256	no additional isolated leptons (pT > 15 GeV)
257
258	same lepton definition as in the physics objects section
259
260	No additional ak5PFjets (pT > 30 GeV; \|eta\| < 2.4)
261	MET>35. for the electron BDT analysis
262	\|cos(theta * )\| (TBO)
263	color flow pull angle (TBO)
264	We don't cut on the transverse mass (for boosted objects cutting on the transverse mass introduces an inefficiency due to the angle between the MET and the lepton being close to 0.)
265	*/
266
267	ok=true;
268	temp.V.firstLepton=pos[0];
269	temp.V.p4 = temp.V.electrons[0].p4+temp.V.mets[0].p4;
270	return temp;
271	}
272
273	public:
274	bool verbose_;
275
276	public:
277	static float getDeltaTheta( const VHbbEvent::SimpleJet & j1, const VHbbEvent::SimpleJet & j2 ) {
278	double deltaTheta = 1e10;
279	TLorentzVector pi(0,0,0,0);
280	TLorentzVector v_j1 = j1.chargedTracksFourMomentum;
281	TLorentzVector v_j2 = j2.chargedTracksFourMomentum;
282
283	if( v_j2.Mag() == 0
284	\|\| v_j1.Mag() == 0 )
285	return deltaTheta = 1e10;
286
287	//use j1 to calculate the pull vector
288	TVector2 t = j1.tVector;
289
290	if( t.Mod() == 0 )
291	return deltaTheta = 1e10;
292
293	Double_t dphi = v_j2.Phi()- v_j1.Phi();
294	if ( dphi > M_PI ) {
295	dphi -= 2.0*M_PI;
296	} else if ( dphi <= -M_PI ) {
297	dphi += 2.0*M_PI;
298	}
299	Double_t deltaeta = v_j2.Rapidity() - v_j1.Rapidity();
300	TVector2 BBdir( deltaeta, dphi );
301
302	deltaTheta = t.DeltaPhi(BBdir);
303
304	return deltaTheta;
305
306	}
307
308	float getHelicity(const VHbbEvent::SimpleJet& j, TVector3 boost) const {
309	double hel = 1e10;
310	TLorentzVector jet = j.p4;
311	jet.Boost( -boost );
312	hel = TMath::Cos( jet.Vect().Angle( boost ) );
313	return hel;
314	}
315
316
317	};
318
319
320
321	#endif
322
323
324
325
326
327