[2] | 1 | #ifndef _SMEARUTIL_H_
|
---|
| 2 | #define _SMEARUTIL_H_
|
---|
| 3 |
|
---|
[260] | 4 | /***********************************************************************
|
---|
| 5 | ** **
|
---|
| 6 | ** /----------------------------------------------\ **
|
---|
| 7 | ** | Delphes, a framework for the fast simulation | **
|
---|
| 8 | ** | of a generic collider experiment | **
|
---|
| 9 | ** \----------------------------------------------/ **
|
---|
| 10 | ** **
|
---|
| 11 | ** **
|
---|
| 12 | ** This package uses: **
|
---|
| 13 | ** ------------------ **
|
---|
| 14 | ** FastJet algorithm: Phys. Lett. B641 (2006) [hep-ph/0512210] **
|
---|
| 15 | ** Hector: JINST 2:P09005 (2007) [physics.acc-ph:0707.1198v2] **
|
---|
| 16 | ** FROG: [hep-ex/0901.2718v1] **
|
---|
| 17 | ** **
|
---|
| 18 | ** ------------------------------------------------------------------ **
|
---|
| 19 | ** **
|
---|
| 20 | ** Main authors: **
|
---|
| 21 | ** ------------- **
|
---|
| 22 | ** **
|
---|
| 23 | ** Severine Ovyn Xavier Rouby **
|
---|
| 24 | ** severine.ovyn@uclouvain.be xavier.rouby@cern **
|
---|
| 25 | ** **
|
---|
| 26 | ** Center for Particle Physics and Phenomenology (CP3) **
|
---|
| 27 | ** Universite catholique de Louvain (UCL) **
|
---|
| 28 | ** Louvain-la-Neuve, Belgium **
|
---|
| 29 | ** **
|
---|
| 30 | ** Copyright (C) 2008-2009, **
|
---|
| 31 | ** All rights reserved. **
|
---|
| 32 | ** **
|
---|
| 33 | ***********************************************************************/
|
---|
[2] | 34 |
|
---|
| 35 | /// \file SmearUtil.h
|
---|
| 36 | /// \brief RESOLution class, and some generic definitions
|
---|
| 37 |
|
---|
| 38 |
|
---|
| 39 | #include <vector>
|
---|
| 40 | #include "TLorentzVector.h"
|
---|
| 41 |
|
---|
[264] | 42 | #include "D_Constants.h"
|
---|
| 43 | #include "CaloUtil.h"
|
---|
[223] | 44 | #include "BlockClasses.h"
|
---|
| 45 | #include "TSimpleArray.h"
|
---|
| 46 | #include "PhysicsTower.hh"
|
---|
[2] | 47 |
|
---|
| 48 | using namespace std;
|
---|
| 49 |
|
---|
[264] | 50 | class D_Particle {
|
---|
[73] | 51 |
|
---|
| 52 | public:
|
---|
[264] | 53 | D_Particle(const TLorentzVector & p, const int pid, const float etacalo, const float phicalo) :
|
---|
| 54 | _fourmomentum(p), _pid(pid), _etaCalo(etacalo), _phiCalo(phicalo) {}
|
---|
| 55 | //D_Particle(const float e, const float eta, const float phi, const float pt, const int pid) :
|
---|
| 56 | // _pid(pid), _etaCalo(UNDEFINED), _phiCalo(UNDEFINED) { TLorentzVector p; p.SetPtEtaPhiE(pt,eta,phi,e); _fourmomentum = p; }
|
---|
| 57 | D_Particle(const float px, const float py, const float pz, const float e, const int pid) :
|
---|
| 58 | _fourmomentum(px,py,pz,e), _pid(pid), _etaCalo(UNDEFINED), _phiCalo(UNDEFINED) {}
|
---|
[73] | 59 |
|
---|
[264] | 60 | const float E() const {return _fourmomentum.E();} // particle energy [GeV]
|
---|
| 61 | const float Px() const {return _fourmomentum.Px();} // horizontal coordinate of momentum [GeV]
|
---|
| 62 | const float Py() const {return _fourmomentum.Py();} // vertical coordinate of momentum [GeV]
|
---|
| 63 | const float Pz() const {return _fourmomentum.Pz();} // longitudinal coordinate of momentum [GeV]
|
---|
| 64 | const float Pt() const {return _fourmomentum.Pt();} // transverse momentum [GeV]
|
---|
| 65 | const float EtaCalo() const {return _etaCalo;} // pseudorapidity
|
---|
| 66 | const float Eta() const {return _fourmomentum.Eta();} // pseudorapidity
|
---|
| 67 | const float PhiCalo() const {return _phiCalo;} // azimuthal angle
|
---|
| 68 | const float Phi() const {return _fourmomentum.Phi();} // azimuthal angle
|
---|
| 69 | const int PID() const {return _pid;} // particle energy in [GeV]
|
---|
| 70 | const TLorentzVector& getFourMomentum() const {return _fourmomentum;}
|
---|
[73] | 71 |
|
---|
| 72 | private:
|
---|
[264] | 73 | TLorentzVector _fourmomentum;
|
---|
[223] | 74 | int _pid;
|
---|
[264] | 75 | float _etaCalo, _phiCalo;
|
---|
[73] | 76 | };
|
---|
| 77 |
|
---|
[2] | 78 | class RESOLution
|
---|
| 79 | {
|
---|
| 80 | public:
|
---|
| 81 | /// Constructor
|
---|
| 82 | RESOLution();
|
---|
[223] | 83 | RESOLution(const RESOLution & DET);
|
---|
| 84 | RESOLution& operator=(const RESOLution& DET);
|
---|
| 85 | ~RESOLution() { delete [] TOWER_eta_edges; delete [] TOWER_dphi;};
|
---|
| 86 |
|
---|
[2] | 87 | // Detector coverage
|
---|
[94] | 88 | float CEN_max_tracker; // tracker pseudorapidity coverage
|
---|
| 89 | float CEN_max_calo_cen; // central calorimeter pseudorapidity coverage
|
---|
| 90 | float CEN_max_calo_fwd; // forward calorimeter pseudorapidity coverage
|
---|
| 91 | float CEN_max_mu; // muon chambers pseudorapidity coverage
|
---|
[2] | 92 |
|
---|
[94] | 93 | float VFD_min_calo_vfd; // very forward calorimeter pseudorapidity coverage
|
---|
| 94 | float VFD_max_calo_vfd; // very forward calorimeter pseudorapidity coverage
|
---|
| 95 | float VFD_min_zdc; // coverage for Zero Degree Calorimeter, for photons and neutrons
|
---|
| 96 | float VFD_s_zdc; // distance of the Zero Degree Calorimeter, from the Interaction poin, in [m]
|
---|
[2] | 97 |
|
---|
[264] | 98 | float RP_220_s; // distance of the RP to the IP, in meters
|
---|
| 99 | float RP_220_x; // distance of the RP to the beam, in meters
|
---|
| 100 | float RP_420_s; // distance of the RP to the IP, in meters
|
---|
| 101 | float RP_420_x; // distance of the RP to the beam, in meters
|
---|
[252] | 102 | string RP_beam1Card; //
|
---|
| 103 | string RP_beam2Card; //
|
---|
[257] | 104 | string RP_IP_name; //
|
---|
[252] | 105 | float RP_offsetEl_s;
|
---|
| 106 | float RP_offsetEl_x;
|
---|
[254] | 107 | float RP_cross_x;
|
---|
| 108 | float RP_cross_y;
|
---|
| 109 | float RP_cross_ang;
|
---|
[62] | 110 |
|
---|
[94] | 111 |
|
---|
[2] | 112 | //energy resolution for electron/photon
|
---|
| 113 | // \sigma/E = C + N/E + S/\sqrt{E}
|
---|
| 114 | float ELG_Scen; // S term for central ECAL
|
---|
| 115 | float ELG_Ncen; // N term for central ECAL
|
---|
| 116 | float ELG_Ccen; // C term for central ECAL
|
---|
| 117 | float ELG_Sfwd; // S term for forward ECAL
|
---|
| 118 | float ELG_Cfwd; // C term for forward ECAL
|
---|
| 119 | float ELG_Nfwd; // N term for central ECAL
|
---|
| 120 |
|
---|
| 121 | //energy resolution for hadrons in ecal/hcal/hf
|
---|
| 122 | // \sigma/E = C + N/E + S/\sqrt{E}
|
---|
| 123 | float HAD_Shcal; // S term for central HCAL // hadronic calorimeter
|
---|
| 124 | float HAD_Nhcal; // N term for central HCAL
|
---|
| 125 | float HAD_Chcal; // C term for central HCAL
|
---|
| 126 | float HAD_Shf; // S term for central HF // forward calorimeter
|
---|
| 127 | float HAD_Nhf; // N term for central HF
|
---|
| 128 | float HAD_Chf; // C term for central HF
|
---|
| 129 |
|
---|
| 130 | // muon smearing
|
---|
| 131 | float MU_SmearPt;
|
---|
| 132 |
|
---|
[94] | 133 | //Magnetic Field information
|
---|
| 134 | int TRACK_radius; //radius of the BField coverage
|
---|
| 135 | int TRACK_length; //length of the BField coverage
|
---|
| 136 | float TRACK_bfield_x;
|
---|
| 137 | float TRACK_bfield_y;
|
---|
| 138 | float TRACK_bfield_z;
|
---|
| 139 | float TRACK_ptmin; // minimal pt needed to reach the calorimeter, in GeV
|
---|
| 140 | int TRACK_eff; // in percent, should be an integer
|
---|
[2] | 141 |
|
---|
[72] | 142 |
|
---|
[94] | 143 | //Define Calorimetric towers
|
---|
| 144 | unsigned int TOWER_number;
|
---|
| 145 | float * TOWER_eta_edges;
|
---|
| 146 | float * TOWER_dphi;
|
---|
[43] | 147 |
|
---|
[94] | 148 | //thresholds for reconstructed objetcs
|
---|
| 149 | float PTCUT_elec;
|
---|
| 150 | float PTCUT_muon;
|
---|
| 151 | float PTCUT_jet;
|
---|
| 152 | float PTCUT_gamma;
|
---|
| 153 | float PTCUT_taujet;
|
---|
[305] | 154 |
|
---|
| 155 | float ISOL_PT; //minimal pt of tracks for isolation criteria
|
---|
| 156 | float ISOL_Cone; //Cone for isolation criteria
|
---|
[321] | 157 | float ISOL_Calo_ET; //minimal tower energy for isolation criteria
|
---|
| 158 | float ISOL_Calo_Cone;//Cone for calorimetric isolation
|
---|
| 159 | unsigned int ISOL_Calo_Grid; //Grid size (N x N) for calorimetric isolation
|
---|
[94] | 160 |
|
---|
[43] | 161 | //General jet variable
|
---|
[94] | 162 | double JET_coneradius;
|
---|
| 163 | int JET_jetalgo;
|
---|
| 164 | double JET_seed;
|
---|
| 165 | double JET_overlap;
|
---|
| 166 |
|
---|
[2] | 167 | // MidPoint algorithm definition
|
---|
[94] | 168 | double JET_M_coneareafraction;
|
---|
| 169 | int JET_M_maxpairsize;
|
---|
| 170 | int JET_M_maxiterations;
|
---|
[2] | 171 | // Define Cone algorithm.
|
---|
[94] | 172 | int JET_C_adjacencycut;
|
---|
| 173 | int JET_C_maxiterations;
|
---|
| 174 | int JET_C_iratch;
|
---|
[44] | 175 | //Define SISCone algorithm.
|
---|
[94] | 176 | int JET_S_npass;
|
---|
| 177 | double JET_S_protojet_ptmin;
|
---|
| 178 |
|
---|
| 179 | //For Tau-jet definition
|
---|
| 180 | // R = sqrt (phi^2 + eta^2)
|
---|
| 181 | float TAU_energy_scone; // radius R of the cone for tau definition, based on energy threshold
|
---|
| 182 | float TAU_track_scone; // radius R of the cone for tau definition, based on track number
|
---|
| 183 | float TAU_track_pt; // minimal pt [GeV] for tracks to be considered in tau definition
|
---|
| 184 | float TAU_energy_frac; // fraction of energy required in the central part of the cone, for tau jets
|
---|
| 185 |
|
---|
| 186 | //tagging definition
|
---|
| 187 | int BTAG_b;
|
---|
| 188 | int BTAG_mistag_c;
|
---|
| 189 | int BTAG_mistag_l;
|
---|
| 190 |
|
---|
[44] | 191 |
|
---|
[94] | 192 | //trigger flag
|
---|
| 193 | int FLAG_trigger; //flag for trigger
|
---|
| 194 | int FLAG_frog; //flag for frog display
|
---|
| 195 | int FLAG_bfield; //flag for bfield propagation
|
---|
| 196 | int FLAG_vfd; //flag for very forward detector
|
---|
[306] | 197 | int FLAG_RP; //flag for very forward detector
|
---|
[307] | 198 | int FLAG_lhco; //flag for very forward detector
|
---|
[94] | 199 |
|
---|
| 200 | int NEvents_Frog;
|
---|
| 201 | float PT_QUARKS_MIN; // minimal pt needed for quarks to reach the tracker, in GeV
|
---|
| 202 |
|
---|
[74] | 203 | // to sort a vector
|
---|
[264] | 204 | //void SortedVector(vector<ParticleUtil> &vect);
|
---|
| 205 | void SortedVector(vector<D_Particle> &vect);
|
---|
[71] | 206 |
|
---|
[2] | 207 | /// Reads the data card for the initialisation of the parameters
|
---|
| 208 | void ReadDataCard(const string datacard);
|
---|
[44] | 209 |
|
---|
| 210 | /// Create the output log file
|
---|
[223] | 211 | void Logfile(const string& LogName);
|
---|
[2] | 212 |
|
---|
| 213 | /// Provides the smeared TLorentzVector for the electrons
|
---|
| 214 | void SmearElectron(TLorentzVector &electron);
|
---|
| 215 |
|
---|
| 216 | /// Provides the smeared TLorentzVector for the muons
|
---|
| 217 | void SmearMu(TLorentzVector &muon);
|
---|
| 218 |
|
---|
| 219 | /// Provides the smeared TLorentzVector for the hadrons
|
---|
| 220 | void SmearHadron(TLorentzVector &hadron, const float frac);
|
---|
| 221 |
|
---|
[223] | 222 | /// For electromagnetic collimation in tau jets
|
---|
[2] | 223 | double EnergySmallCone(const vector<PhysicsTower> &towers, const float eta, const float phi);
|
---|
| 224 |
|
---|
[223] | 225 | /// Number of tracks in tau jet algo
|
---|
[287] | 226 | //unsigned int NumTracks(float& charge, const vector<TLorentzVector> &tracks, const float pt_track, const float eta, const float phi);
|
---|
| 227 | unsigned int NumTracks(float& charge, const vector<TRootTracks> &tracks, const float pt_track, const float eta, const float phi);
|
---|
[2] | 228 |
|
---|
[223] | 229 | /// b-jets
|
---|
[319] | 230 | int Bjets(const TSimpleArray<GenParticle> &subarray, const float& eta, const float& phi);
|
---|
[2] | 231 |
|
---|
[223] | 232 | /// b-tag efficiency and misidentification
|
---|
[319] | 233 | bool Btaggedjet(const TLorentzVector &JET, const TSimpleArray<GenParticle> &subarray);
|
---|
[2] | 234 |
|
---|
[321] | 235 | /// Lepton isolation based on tracking
|
---|
| 236 | bool Isolation(const D_Particle& part, const vector<TRootTracks> &tracks, const float& pt_second_track, const float& isolCone, float& ptiso);
|
---|
[31] | 237 |
|
---|
[321] | 238 | /// Lepton isolation based on calorimetry (optional. Default: off)
|
---|
| 239 | float CaloIsolation(const D_Particle& part, const D_CaloTowerList & towers);
|
---|
| 240 |
|
---|
[71] | 241 | //********************* returns a segmented value for eta and phi, for calo towers *****
|
---|
| 242 | void BinEtaPhi(const float phi, const float eta, float& iPhi, float& iEta);
|
---|
| 243 |
|
---|
[2] | 244 | };
|
---|
| 245 |
|
---|
| 246 | // ** returns the sign (+1 or -1) or an integer
|
---|
| 247 | int sign(const int myint);
|
---|
| 248 | int sign(const float myfloat);
|
---|
| 249 |
|
---|
| 250 | // **************************** Return the Delta Phi****************************
|
---|
| 251 | float DeltaPhi(const float phi1, const float phi2);
|
---|
| 252 |
|
---|
| 253 | // **************************** Returns the Delta R****************************
|
---|
| 254 | float DeltaR(const float phi1, const float eta1, const float phi2, const float eta2);
|
---|
| 255 |
|
---|
| 256 | //************* Returns an array of the quarks sitting within the tracker acceptance ***************
|
---|
[270] | 257 | int ChargeVal(const int pid);
|
---|
[2] | 258 |
|
---|
| 259 | #endif
|
---|