Changeset 223 in svn for trunk/interface/SmearUtil.h
- Timestamp:
- Feb 2, 2009, 12:39:42 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/interface/SmearUtil.h
r198 r223 19 19 #include "TLorentzVector.h" 20 20 21 #include "Utilities/ExRootAnalysis/interface/BlockClasses.h" 22 #include "Utilities/ExRootAnalysis/interface/TSimpleArray.h" 23 24 #include "Utilities/Fastjet/plugins/CDFCones/interface/PhysicsTower.hh" 25 26 #include <iostream> 27 #include <sstream> 28 #include <fstream> 29 #include <iomanip> 21 #include "BlockClasses.h" 22 #include "TSimpleArray.h" 23 #include "PhysicsTower.hh" 30 24 31 25 using namespace std; … … 34 28 35 29 public: 36 ParticleUtil(const TLorentzVector &genMomentum,const int &pid) : TLVector(genMomentum) , Pid(pid) {} 37 38 float E() {return TLVector.E();} // particle energy [GeV] 39 float Px() {return TLVector.Px();} // horizontal coordinate of momentum [GeV] 40 float Py() {return TLVector.Py();} // vertical coordinate of momentum [GeV] 41 float Pz() {return TLVector.Pz();} // longitudinal coordinate of momentum [GeV] 42 float Pt() {return TLVector.Pt();} // transverse momentum [GeV] 43 float Eta() {return TLVector.Eta();} // pseudorapidity [GeV] 44 float Phi() {return TLVector.Phi();} // azimuthal angle [GeV] 45 float PID() {return Pid;} // particle energy in [GeV] 30 ParticleUtil(const TLorentzVector &genMomentum, int pid); 31 32 float E() {return _e;} // particle energy [GeV] 33 float Px() {return _px;} // horizontal coordinate of momentum [GeV] 34 float Py() {return _py;} // vertical coordinate of momentum [GeV] 35 float Pz() {return _pz;} // longitudinal coordinate of momentum [GeV] 36 float Pt() {return _pt;} // transverse momentum [GeV] 37 float EtaCalo() {return _etaCalo;} // pseudorapidity 38 float Eta() {return _eta;} // pseudorapidity 39 float PhiCalo() {return _phiCalo;} // azimuthal angle 40 float Phi() {return _phi;} // azimuthal angle 41 int PID() {return _pid;} // particle energy in [GeV] 46 42 47 43 private: 48 TLorentzVector TLVector; 49 int Pid; 44 float _e, _px, _py, _pz, _pt; 45 float _eta, _etaCalo, _phi, _phiCalo; 46 int _pid; 50 47 }; 51 48 … … 56 53 /// Constructor 57 54 RESOLution(); 58 55 RESOLution(const RESOLution & DET); 56 RESOLution& operator=(const RESOLution& DET); 57 ~RESOLution() { delete [] TOWER_eta_edges; delete [] TOWER_dphi;}; 58 59 59 // Detector coverage 60 60 float CEN_max_tracker; // tracker pseudorapidity coverage … … 155 155 156 156 int NEvents_Frog; 157 158 159 160 161 162 157 float PT_QUARKS_MIN; // minimal pt needed for quarks to reach the tracker, in GeV 163 158 … … 169 164 170 165 /// Create the output log file 171 void Logfile( stringLogName);166 void Logfile(const string& LogName); 172 167 173 168 /// Provides the smeared TLorentzVector for the electrons … … 180 175 void SmearHadron(TLorentzVector &hadron, const float frac); 181 176 182 // *****************************fonction pour avoir les taus************************************177 /// For electromagnetic collimation in tau jets 183 178 double EnergySmallCone(const vector<PhysicsTower> &towers, const float eta, const float phi); 184 179 185 // ***************** Fonction pour avoir le nombre de traces pour les taus****************************180 /// Number of tracks in tau jet algo 186 181 unsigned int NumTracks(const vector<TLorentzVector> &tracks, const float pt_track, const float eta, const float phi); 187 182 188 // **********************fonction pour avoir les b-jets******************************183 /// b-jets 189 184 int Bjets(const TSimpleArray<TRootGenParticle> &subarray, const float eta, const float phi); 190 185 191 // ******************** retourne l'efficacite de b-tagging ******************************186 /// b-tag efficiency and misidentification 192 187 bool Btaggedjet(const TLorentzVector &JET, const TSimpleArray<TRootGenParticle> &subarray); 193 188 194 // ******************************isolation criteria**************************************195 bool Isolation( Float_t phi,Float_t eta,const vector<TLorentzVector> &tracks,float PT_TRACK2);189 /// Lepton isolation 190 bool Isolation(const float phi, const float eta,const vector<TLorentzVector> &tracks,float PT_TRACK2); 196 191 197 192 //********************* returns a segmented value for eta and phi, for calo towers ***** … … 226 221 227 222 const double speed_of_light = 299792458; // m/s 223 const float UNDEFINED=-9999.; 224 228 225 229 226 #ifndef __PI__
Note:
See TracChangeset
for help on using the changeset viewer.