Changes in display/Delphes3DGeometry.cc [f53a4d2:1fa50c2] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
display/Delphes3DGeometry.cc
rf53a4d2 r1fa50c2 17 17 */ 18 18 19 #include "display/Delphes3DGeometry.h" 19 20 #include <set> 20 21 #include <map> … … 24 25 #include <sstream> 25 26 #include <cassert> 26 27 #include "TAxis.h"28 27 #include "TGeoManager.h" 29 28 #include "TGeoVolume.h" … … 35 34 #include "TGeoCone.h" 36 35 #include "TGeoArb8.h" 36 #include "external/ExRootAnalysis/ExRootConfReader.h" 37 #include "classes/DelphesClasses.h" 37 38 #include "TF2.h" 38 #include "TFormula.h"39 39 #include "TH1F.h" 40 40 #include "TMath.h" 41 #include "TString.h"42 43 #include "display/Delphes3DGeometry.h"44 45 #include "classes/DelphesClasses.h"46 #include "external/ExRootAnalysis/ExRootConfReader.h"47 41 48 42 using namespace std; … … 96 90 tk_length_ = confReader->GetDouble(Form("%s::HalfLength",ParticlePropagator), 3.0)*100.; // tk_length 97 91 tk_Bz_ = confReader->GetDouble("ParticlePropagator::Bz", 0.0); // tk_Bz 98 99 TString buffer; 100 const char *it; 101 102 92 103 93 { 104 94 TString tkEffFormula = confReader->GetString(Form("%s::EfficiencyFormula",TrackingEfficiency),"abs(eta)<3.0"); … … 106 96 tkEffFormula.ReplaceAll("eta","y"); 107 97 tkEffFormula.ReplaceAll("phi","0."); 108 109 buffer.Clear(); 110 for(it = tkEffFormula.Data(); *it; ++it) 111 { 112 if(*it == ' ' || *it == '\t' || *it == '\r' || *it == '\n' || *it == '\\' ) continue; 113 buffer.Append(*it); 114 } 115 116 TF2* tkEffFunction = new TF2("tkEff",buffer,0,1000,-10,10); 98 TF2* tkEffFunction = new TF2("tkEff",tkEffFormula,0,1000,-10,10); 117 99 TH1F etaHisto("eta","eta",100,5.,-5.); 118 100 Double_t pt,eta; … … 136 118 muonEffFormula.ReplaceAll("eta","y"); 137 119 muonEffFormula.ReplaceAll("phi","0."); 138 139 buffer.Clear(); 140 for(it = muonEffFormula.Data(); *it; ++it) 141 { 142 if(*it == ' ' || *it == '\t' || *it == '\r' || *it == '\n' || *it == '\\' ) continue; 143 buffer.Append(*it); 144 } 145 146 TF2* muEffFunction = new TF2("muEff",buffer,0,1000,-10,10); 120 TF2* muEffFunction = new TF2("muEff",muonEffFormula,0,1000,-10,10); 147 121 TH1F etaHisto("eta2","eta2",100,5.,-5.); 148 122 Double_t pt,eta;
Note:
See TracChangeset
for help on using the changeset viewer.