Changeset 341014c in git for modules/ImpactParameterSmearing.cc
- Timestamp:
- Feb 12, 2019, 9:29:17 PM (6 years ago)
- Branches:
- ImprovedOutputFile, Timing, llp, master
- Children:
- 6455202
- Parents:
- 45e58be
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/ImpactParameterSmearing.cc
r45e58be r341014c 25 25 */ 26 26 27 28 27 #include "modules/ImpactParameterSmearing.h" 29 28 … … 32 31 #include "classes/DelphesFormula.h" 33 32 33 #include "ExRootAnalysis/ExRootClassifier.h" 34 #include "ExRootAnalysis/ExRootFilter.h" 34 35 #include "ExRootAnalysis/ExRootResult.h" 35 #include "ExRootAnalysis/ExRootFilter.h"36 #include "ExRootAnalysis/ExRootClassifier.h"37 36 37 #include "TDatabasePDG.h" 38 #include "TFormula.h" 39 #include "TLorentzVector.h" 38 40 #include "TMath.h" 41 #include "TObjArray.h" 42 #include "TRandom3.h" 39 43 #include "TString.h" 40 #include "TFormula.h"41 #include "TRandom3.h"42 #include "TObjArray.h"43 #include "TDatabasePDG.h"44 #include "TLorentzVector.h"45 44 46 45 #include <algorithm> 47 #include <stdexcept>48 46 #include <iostream> 49 47 #include <sstream> 48 #include <stdexcept> 50 49 51 50 using namespace std; … … 100 99 101 100 fItInputArray->Reset(); 102 while((candidate = static_cast<Candidate *>(fItInputArray->Next())))101 while((candidate = static_cast<Candidate *>(fItInputArray->Next()))) 103 102 { 104 103 105 104 // take momentum before smearing (otherwise apply double smearing on d0) 106 particle = static_cast<Candidate *>(candidate->GetCandidates()->At(0));105 particle = static_cast<Candidate *>(candidate->GetCandidates()->At(0)); 107 106 108 107 const TLorentzVector &candidateMomentum = particle->Momentum; … … 112 111 phi = candidateMomentum.Phi(); 113 112 e = candidateMomentum.E(); 114 113 115 114 px = candidateMomentum.Px(); 116 115 py = candidateMomentum.Py(); 117 116 118 117 // calculate coordinates of closest approach to track circle in transverse plane xd, yd, zd 119 xd = 120 yd = 121 zd = 118 xd = candidate->Xd; 119 yd = candidate->Yd; 120 zd = candidate->Zd; 122 121 123 122 // calculate smeared values … … 131 130 132 131 // calculate impact parameter (after-smearing) 133 d0 = (xd *py - yd*px)/pt;132 d0 = (xd * py - yd * px) / pt; 134 133 135 134 dd0 = gRandom->Gaus(0.0, fFormula->Eval(pt, eta, phi, e)); … … 138 137 mother = candidate; 139 138 140 candidate = static_cast<Candidate *>(candidate->Clone());139 candidate = static_cast<Candidate *>(candidate->Clone()); 141 140 candidate->Xd = xd; 142 141 candidate->Yd = yd;
Note:
See TracChangeset
for help on using the changeset viewer.