Changeset 25291dd in git
- Timestamp:
- Oct 17, 2020, 4:05:03 AM (4 years ago)
- Children:
- 44f4a4e
- Parents:
- 4e5ed02 (diff), be1a3be (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the(diff)
links above to see all the changes relative to each parent. - Files:
-
- 23 added
- 8 deleted
- 28 edited
Legend:
- Unmodified
- Added
- Removed
-
Makefile
r4e5ed02 r25291dd 89 89 DISPLAYLIB = libDelphesDisplay.lib 90 90 91 VERSION = $(shell cat VERSION)91 VERSION = x.y.z 92 92 DISTDIR = Delphes-$(VERSION) 93 93 DISTTAR = $(DISTDIR).tar.gz … … 174 174 tmp/examples/Example1.$(ObjSuf): \ 175 175 examples/Example1.cpp \ 176 classes/DelphesClasses.h \177 external/ExRootAnalysis/ExRootResult.h \178 external/ExRootAnalysis/ExRootTreeBranch.h \179 external/ExRootAnalysis/ExRootTreeReader.h \180 external/ExRootAnalysis/ExRootTreeWriter.h \181 external/ExRootAnalysis/ExRootUtilities.h182 DelphesValidation$(ExeSuf): \183 tmp/validation/DelphesValidation.$(ObjSuf)184 185 tmp/validation/DelphesValidation.$(ObjSuf): \186 validation/DelphesValidation.cpp \187 176 classes/DelphesClasses.h \ 188 177 external/ExRootAnalysis/ExRootResult.h \ … … 199 188 stdhep2pileup$(ExeSuf) \ 200 189 CaloGrid$(ExeSuf) \ 201 Example1$(ExeSuf) \ 202 DelphesValidation$(ExeSuf) 190 Example1$(ExeSuf) 203 191 204 192 EXECUTABLE_OBJ += \ … … 210 198 tmp/converters/stdhep2pileup.$(ObjSuf) \ 211 199 tmp/examples/CaloGrid.$(ObjSuf) \ 212 tmp/examples/Example1.$(ObjSuf) \ 213 tmp/validation/DelphesValidation.$(ObjSuf) 200 tmp/examples/Example1.$(ObjSuf) 214 201 215 202 DelphesHepMC$(ExeSuf): \ … … 2434 2421 @echo ">> Building $(DISTTAR)" 2435 2422 @mkdir -p $(DISTDIR) 2436 @cp -a AUTHORS CHANGELOG CMakeLists.txt COPYING DelphesEnv.sh LICENSE NOTICE README README_4LHCb VERSIONMakefile MinBias.pileup configure cards classes converters display doc examples external modules python readers validation $(DISTDIR)2423 @cp -a AUTHORS CHANGELOG CMakeLists.txt COPYING DelphesEnv.sh LICENSE NOTICE README README_4LHCb Makefile MinBias.pileup configure cards classes converters display doc examples external modules python readers validation $(DISTDIR) 2437 2424 @find $(DISTDIR) -depth -name .\* -exec rm -rf {} \; 2438 2425 @tar -czf $(DISTTAR) $(DISTDIR) -
cards/FCC/FCChh.tcl
r4e5ed02 r25291dd 4 4 # Main authors: Michele Selvaggi (CERN) 5 5 # 6 # Released on: Dec. 1st, 2017 6 # Released on: October 14th, 2020 7 # 8 # - fix muon resolution at high pT 9 # - updated btagging, tau tagging and photon ID 10 # 7 11 # 8 12 # Configuration: FCC-hh baseline detector … … 1049 1053 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) + 1050 1054 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) + 1051 (abs(eta) <= 2.5) * (pt > 10.0) * (0.9 5) +1055 (abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + 1052 1056 1053 1057 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) + … … 1138 1142 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.01)*(1.0 - pt/15000.) + \ 1139 1143 (abs(eta) < 2.5) * (pt > 15000.0) * (0.00) + \ 1140 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.0 075) + \1141 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.0 075)*(1.0 - pt/15000.) + \1144 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.01) + \ 1145 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.01)*(1.0 - pt/15000.) + \ 1142 1146 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) + \ 1143 1147 (abs(eta) > 4.0) * (0.00)} … … 1146 1150 1147 1151 (pt <= 10.0) * (0.00) + 1148 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0. 05) + \1149 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0. 05)*(1.0 - pt/15000.) + \1152 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.15) + \ 1153 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.15)*(1.0 - pt/15000.) + \ 1150 1154 (abs(eta) < 2.5) * (pt > 15000.0) * (0.000) + \ 1151 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0. 03) + \1152 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0. 03)*(1.0 - pt/15000.) + \1155 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.10) + \ 1156 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.10)*(1.0 - pt/15000.) + \ 1153 1157 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) + \ 1154 1158 (abs(eta) > 4.0) * (0.00)} … … 1157 1161 1158 1162 (pt <= 10.0) * (0.00) + 1159 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.8 5) +1160 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.8 5)*(1.0 - pt/15000.) +1163 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.82) + 1164 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.82)*(1.0 - pt/15000.) + 1161 1165 (abs(eta) < 2.5) * (pt > 15000.0) * (0.000) + 1162 1166 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) + … … 1217 1221 1218 1222 (pt <= 10.0) * (0.00) + 1219 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.0 1) + \1220 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.0 1) *(8./9. - pt/30000.) + \1223 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.02) + \ 1224 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.02) *(8./9. - pt/30000.) + \ 1221 1225 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \ 1222 1226 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \ … … 1228 1232 1229 1233 (pt <= 10.0) * (0.00) + 1230 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.00 5) + \1231 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.00 5) *(8./9. - pt/30000.) + \1234 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.001) + \ 1235 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.001) *(8./9. - pt/30000.) + \ 1232 1236 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \ 1233 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00 375) + \1234 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00 375)*(8./9. - pt/30000.) + \1237 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.001) + \ 1238 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.001)*(8./9. - pt/30000.) + \ 1235 1239 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \ 1236 1240 (abs(eta) > 4.0) * (0.00)} -
cards/FCC/muonMomentumResolutionVsP.tcl
r4e5ed02 r25291dd 12 12 set ResolutionFormula { 13 13 14 ( abs(eta) < 1.35 ) * 15 16 (sqrt(0.0000364164 + ( 17 9.06262e-8 *pt^2* cosh( 18 eta)^2 *(2.82074e-7/sin(2*atan(exp(-abs(eta))))^2 + ( 19 504.525 *(1/400000000 + (0.117945* 1/cosh(eta)^2)/( 20 pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/ 21 sin(2*atan(exp(-abs(eta))))^2) *sin(2*atan(exp(-abs(eta))))^2)/( 22 0.00516429/sin(2*atan(exp(-abs(eta))))^2 + ( 23 96868.8 *(1/400000000 + (0.117945 * 1/cosh(eta)^2)/( 24 pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/ 25 sin(2*atan(exp(-abs(eta))))^2)) 26 27 ) + 14 15 ( abs(eta) < 1.35 ) * 16 17 (sqrt(0.0000364164 + ( 18 3*9.06262e-8 *pt^2* cosh( 19 eta)^2 *(2.82074e-7/sin(2*atan(exp(-abs(eta))))^2 + ( 20 504.525 *(1/400000000 + (0.117945* 1/cosh(eta)^2)/( 21 pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/ 22 sin(2*atan(exp(-abs(eta))))^2) *sin(2*atan(exp(-abs(eta))))^2)/( 23 0.00516429/sin(2*atan(exp(-abs(eta))))^2 + ( 24 96868.8 *(1/400000000 + 5*(0.117945 * 1/cosh(eta)^2)/( 25 pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/ 26 sin(2*atan(exp(-abs(eta))))^2)) 27 28 ) + 28 29 29 30 ( abs(eta) > 1.35 && abs(eta) < 2.00) * -
cards/ILCgen/ILCgen_BTagging_N.tcl
r4e5ed02 r25291dd 104 104 source ILCgen/ILCgen_BTagging_50.tcl 105 105 } 106 -
cards/ILCgen/ILCgen_BeamCalEfficiency.tcl
r4e5ed02 r25291dd 2 2 # Based on IDR Figure 8.8a and Moritz Hebermehl PhD Thesis Figure 4.10 3 3 # A.F.Zarnecki, June 2020 4 5 # set EfficiencyFormula {efficiency formula as a function of eta and pt}6 4 7 5 set EfficiencyFormula { -
cards/ILCgen/ILCgen_BeamCalF_Binning.tcl
r4e5ed02 r25291dd 2 2 # corrected by A.F.Zarnecki 3 3 4 set pi [expr {acos(-1)} ]5 4 set pi [expr {acos(-1)} ] 5 6 6 # BeamCal eta range 4.0 - 5.8 -> 2.099-0.347 =1.752deg 7 7 # -
cards/ILCgen/ILCgen_BeamCalR_Binning.tcl
r4e5ed02 r25291dd 2 2 # corrected by A.F.Zarnecki 3 3 4 set pi [expr {acos(-1)} ]5 4 set pi [expr {acos(-1)} ] 5 6 6 # BeamCal eta range 4.0 - 5.8 -> 2.099-0.347 =1.752deg 7 7 # -
cards/ILCgen/ILCgen_BeamCal_EnergyFractions.tcl
r4e5ed02 r25291dd 1 1 #P.Sopicki 2 # default energy fractions {abs(PDG code)} {fraction of energy deposited in BeamCal}3 2 4 3 add EnergyFraction {0} {0.0} … … 20 19 add EnergyFraction {310} {0.3} 21 20 add EnergyFraction {3122} {0.3} 22 -
cards/ILCgen/ILCgen_BeamCal_Resolution.tcl
r4e5ed02 r25291dd 1 1 # P.Sopicki, corrected by A.F.Zarnecki 2 2 # BeamCal resolution 3 4 # set ResolutionFormula {resolution formula as a function of eta and energy}5 3 6 4 set ResolutionFormula { … … 8 6 (abs(eta) > 4.8 && abs(eta) <= 5.8) * sqrt(energy^2*0.03^2 + energy*0.45^2) 9 7 } 10 -
cards/ILCgen/ILCgen_CTagging_N.tcl
r4e5ed02 r25291dd 104 104 source ILCgen/ILCgen_CTagging_20.tcl 105 105 } 106 -
cards/ILCgen/ILCgen_ChrgHadMomentumSmearing.tcl
r4e5ed02 r25291dd 6 6 7 7 set ResolutionFormula { 8 9 (abs(eta) <= 2.0) * 10 exp( -6.82012 - 0.284133*abs(eta) + 0.64585*eta*eta 8 (abs(eta) <= 2.0) * 9 exp( -6.82012 - 0.284133*abs(eta) + 0.64585*eta*eta 11 10 + (0.829675 -0.188197*abs(eta) + 0.106226 * eta*eta - 1.0) * log(pt) 12 11 + (0.0766417-0.00435667*abs(eta) + 0.00174804*eta*eta) * log(pt)*log(pt)) 13 14 + 15 16 (abs(eta) > 2.0 && abs(eta) <= 3.0) * 17 exp( -4.54609 12 + 13 (abs(eta) > 2.0 && abs(eta) <= 3.0) * 14 exp( -4.54609 18 15 + ( 0.87526 - 1.0) * log(pt) 19 16 + (-0.0800672 + 0.0891669 * abs(eta)) * log(pt)*log(pt)) 20 21 } 22 17 } -
cards/ILCgen/ILCgen_ChrgHadTrackingEff.tcl
r4e5ed02 r25291dd 4 4 # A.F.Zarnecki, June 12, 2020 5 5 6 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 7 (abs(eta) <= 1.83) * (pt > 0.1 && pt <= 0.2) * (0.70) + 8 (abs(eta) <= 1.83) * (pt > 0.2 && pt <= 0.3) * (0.93) + 9 (abs(eta) <= 1.83) * (pt > 0.3 && pt <= 0.7) * (0.995) + 10 (abs(eta) <= 1.83) * (pt > 0.7) * (1.00) + 11 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.1 && pt <= 0.2) * (0.697) + 12 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.2 && pt <= 0.3) * (0.925) + 13 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.3 && pt <= 0.7) * (0.99) + 14 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.7) * (0.995) + 15 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.1 && pt <= 0.2) * (0.665) + 16 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.2 && pt <= 0.3) * (0.884) + 17 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.3 && pt <= 0.7) * (0.945) + 18 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.7) * (0.95) + 19 (abs(eta) > 3.00) * (pt > 0.1) * (0.00) } 6 set EfficiencyFormula { 7 (pt <= 0.1) * (0.00) + 8 (abs(eta) <= 1.83) * (pt > 0.1 && pt <= 0.2) * (0.70) + 9 (abs(eta) <= 1.83) * (pt > 0.2 && pt <= 0.3) * (0.93) + 10 (abs(eta) <= 1.83) * (pt > 0.3 && pt <= 0.7) * (0.995) + 11 (abs(eta) <= 1.83) * (pt > 0.7) * (1.00) + 12 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.1 && pt <= 0.2) * (0.697) + 13 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.2 && pt <= 0.3) * (0.925) + 14 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.3 && pt <= 0.7) * (0.99) + 15 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.7) * (0.995) + 16 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.1 && pt <= 0.2) * (0.665) + 17 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.2 && pt <= 0.3) * (0.884) + 18 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.3 && pt <= 0.7) * (0.945) + 19 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.7) * (0.95) + 20 (abs(eta) > 3.00) * (pt > 0.1) * (0.00) 21 } -
cards/ILCgen/ILCgen_ECAL_EnergyFractions.tcl
r4e5ed02 r25291dd 1 1 #P.Sopicki 2 # default energy fractions {abs(PDG code)} {fraction in ECAL or LumiCal} 2 3 3 add EnergyFraction {0} {0.0} 4 4 # energy fractions for e, gamma and pi0 -
cards/ILCgen/ILCgen_ECAL_Resolution.tcl
r4e5ed02 r25291dd 2 2 # ECAL and LumiCal resolution: same formula used at the moment 3 3 # 4 # set ResolutionFormula {resolution formula as a function of eta and energy} 4 5 5 set ResolutionFormula { 6 6 (abs(eta) <= 3 ) * sqrt(energy^2*0.01^2 + energy*0.17^2)+ 7 7 (abs(eta) > 3 && abs(eta) <= 4 ) * sqrt(energy^2*0.01^2 + energy*0.17^2) 8 }8 } -
cards/ILCgen/ILCgen_ElectronEfficiency.tcl
r4e5ed02 r25291dd 1 2 # set EfficiencyFormula {efficiency formula as a function of eta and pt}3 1 4 2 # efficiency formula for electrons (ECAL and LumiCal) -
cards/ILCgen/ILCgen_ElectronMomentumSmearing.tcl
r4e5ed02 r25291dd 6 6 7 7 set ResolutionFormula { 8 9 (abs(eta) <= 2.0) * 10 exp( -6.82012 - 0.284133*abs(eta) + 0.64585*eta*eta 8 (abs(eta) <= 2.0) * 9 exp( -6.82012 - 0.284133*abs(eta) + 0.64585*eta*eta 11 10 + (0.829675 -0.188197*abs(eta) + 0.106226 * eta*eta - 1.0) * log(pt) 12 11 + (0.0766417-0.00435667*abs(eta) + 0.00174804*eta*eta) * log(pt)*log(pt)) 13 14 + 15 16 (abs(eta) > 2.0 && abs(eta) <= 3.0) * 17 exp( -4.54609 12 + 13 (abs(eta) > 2.0 && abs(eta) <= 3.0) * 14 exp( -4.54609 18 15 + ( 0.87526 - 1.0) * log(pt) 19 16 + (-0.0800672 + 0.0891669 * abs(eta)) * log(pt)*log(pt)) 20 21 } 22 17 } -
cards/ILCgen/ILCgen_ElectronTrackingEff.tcl
r4e5ed02 r25291dd 4 4 # A.F.Zarnecki, June 12, 2020 5 5 6 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 7 (abs(eta) <= 1.83) * (pt > 0.1 && pt <= 0.2) * (0.70) + 8 (abs(eta) <= 1.83) * (pt > 0.2 && pt <= 0.3) * (0.93) + 9 (abs(eta) <= 1.83) * (pt > 0.3 && pt <= 0.7) * (0.995) + 10 (abs(eta) <= 1.83) * (pt > 0.7) * (1.00) + 11 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.1 && pt <= 0.2) * (0.697) + 12 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.2 && pt <= 0.3) * (0.925) + 13 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.3 && pt <= 0.7) * (0.99) + 14 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.7) * (0.995) + 15 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.1 && pt <= 0.2) * (0.665) + 16 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.2 && pt <= 0.3) * (0.884) + 17 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.3 && pt <= 0.7) * (0.945) + 18 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.7) * (0.95) + 19 (abs(eta) > 3.00) * (pt > 0.1) * (0.00) } 6 set EfficiencyFormula { 7 (pt <= 0.1) * (0.00) + 8 (abs(eta) <= 1.83) * (pt > 0.1 && pt <= 0.2) * (0.70) + 9 (abs(eta) <= 1.83) * (pt > 0.2 && pt <= 0.3) * (0.93) + 10 (abs(eta) <= 1.83) * (pt > 0.3 && pt <= 0.7) * (0.995) + 11 (abs(eta) <= 1.83) * (pt > 0.7) * (1.00) + 12 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.1 && pt <= 0.2) * (0.697) + 13 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.2 && pt <= 0.3) * (0.925) + 14 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.3 && pt <= 0.7) * (0.99) + 15 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.7) * (0.995) + 16 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.1 && pt <= 0.2) * (0.665) + 17 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.2 && pt <= 0.3) * (0.884) + 18 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.3 && pt <= 0.7) * (0.945) + 19 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.7) * (0.95) + 20 (abs(eta) > 3.00) * (pt > 0.1) * (0.00) 21 } -
cards/ILCgen/ILCgen_HCAL_EnergyFractions.tcl
r4e5ed02 r25291dd 1 1 #P.Sopicki 2 # default energy fractions {abs(PDG code)} {fraction in HCAL or LHCAL} 2 3 3 add EnergyFraction {0} {1.0} 4 4 # energy fractions for e, gamma and pi0 -
cards/ILCgen/ILCgen_HCAL_Resolution.tcl
r4e5ed02 r25291dd 2 2 # HCAL and LHCAL resolution: same formula used at the moment 3 3 # 4 # set ResolutionFormula {resolution formula as a function of eta and energy} 4 5 5 set ResolutionFormula { 6 6 (abs(eta) <= 2.8 ) * sqrt(energy^2*0.017^2 + energy*0.45^2)+ 7 7 (abs(eta) > 2.8 && abs(eta)<=3.8 ) * sqrt(energy^2*0.017^2 + energy*0.45^2) 8 } 9 8 } -
cards/ILCgen/ILCgen_MuonMomentumSmearing.tcl
r4e5ed02 r25291dd 6 6 7 7 set ResolutionFormula { 8 9 (abs(eta) <= 2.0) * 10 exp( -6.82012 - 0.284133*abs(eta) + 0.64585*eta*eta 8 (abs(eta) <= 2.0) * 9 exp( -6.82012 - 0.284133*abs(eta) + 0.64585*eta*eta 11 10 + (0.829675 -0.188197*abs(eta) + 0.106226 * eta*eta - 1.0) * log(pt) 12 11 + (0.0766417-0.00435667*abs(eta) + 0.00174804*eta*eta) * log(pt)*log(pt)) 13 14 + 15 16 (abs(eta) > 2.0 && abs(eta) <= 3.0) * 17 exp( -4.54609 12 + 13 (abs(eta) > 2.0 && abs(eta) <= 3.0) * 14 exp( -4.54609 18 15 + ( 0.87526 - 1.0) * log(pt) 19 16 + (-0.0800672 + 0.0891669 * abs(eta)) * log(pt)*log(pt)) 20 21 } 22 17 } -
cards/ILCgen/ILCgen_MuonTrackingEff.tcl
r4e5ed02 r25291dd 4 4 # A.F.Zarnecki, June 12, 2020 5 5 6 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 7 (abs(eta) <= 1.83) * (pt > 0.1 && pt <= 0.2) * (0.70) + 8 (abs(eta) <= 1.83) * (pt > 0.2 && pt <= 0.3) * (0.93) + 9 (abs(eta) <= 1.83) * (pt > 0.3 && pt <= 0.7) * (0.995) + 10 (abs(eta) <= 1.83) * (pt > 0.7) * (1.00) + 11 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.1 && pt <= 0.2) * (0.697) + 12 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.2 && pt <= 0.3) * (0.925) + 13 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.3 && pt <= 0.7) * (0.99) + 14 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.7) * (0.995) + 15 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.1 && pt <= 0.2) * (0.665) + 16 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.2 && pt <= 0.3) * (0.884) + 17 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.3 && pt <= 0.7) * (0.945) + 18 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.7) * (0.95) + 19 (abs(eta) > 3.00) * (pt > 0.1) * (0.00) } 6 set EfficiencyFormula { 7 (pt <= 0.1) * (0.00) + 8 (abs(eta) <= 1.83) * (pt > 0.1 && pt <= 0.2) * (0.70) + 9 (abs(eta) <= 1.83) * (pt > 0.2 && pt <= 0.3) * (0.93) + 10 (abs(eta) <= 1.83) * (pt > 0.3 && pt <= 0.7) * (0.995) + 11 (abs(eta) <= 1.83) * (pt > 0.7) * (1.00) + 12 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.1 && pt <= 0.2) * (0.697) + 13 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.2 && pt <= 0.3) * (0.925) + 14 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.3 && pt <= 0.7) * (0.99) + 15 (abs(eta) > 1.83 && abs(eta) <= 2.65) * (pt > 0.7) * (0.995) + 16 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.1 && pt <= 0.2) * (0.665) + 17 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.2 && pt <= 0.3) * (0.884) + 18 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.3 && pt <= 0.7) * (0.945) + 19 (abs(eta) > 2.65 && abs(eta) <= 3.00) * (pt > 0.7) * (0.95) + 20 (abs(eta) > 3.00) * (pt > 0.1) * (0.00) 21 } -
cards/ILCgen/ILCgen_PhotonEfficiency.tcl
r4e5ed02 r25291dd 1 2 # set EfficiencyFormula {efficiency formula as a function of eta and pt} 3 4 # efficiency formula for photons (ECAL + LumiCal) 1 # set EfficiencyFormula (efficiency formula as a function of eta and pt) 2 # efficiency formula for photons (ECAL + LumiCal) 5 3 6 4 set EfficiencyFormula { -
doc/genMakefile.tcl
r4e5ed02 r25291dd 282 282 DISPLAYLIB = libDelphesDisplay.lib 283 283 284 VERSION = $(shell cat VERSION)284 VERSION = x.y.z 285 285 DISTDIR = Delphes-$(VERSION) 286 286 DISTTAR = $(DISTDIR).tar.gz … … 290 290 } 291 291 292 executableDeps {converters/*.cpp} {examples/*.cpp} {validation/*.cpp}292 executableDeps {converters/*.cpp} {examples/*.cpp} 293 293 294 294 executableDeps {readers/DelphesHepMC.cpp} {readers/DelphesLHEF.cpp} {readers/DelphesSTDHEP.cpp} {readers/DelphesROOT.cpp} … … 434 434 @echo ">> Building $(DISTTAR)" 435 435 @mkdir -p $(DISTDIR) 436 @cp -a AUTHORS CHANGELOG CMakeLists.txt COPYING DelphesEnv.sh LICENSE NOTICE README README_4LHCb VERSIONMakefile MinBias.pileup configure cards classes converters display doc examples external modules python readers validation $(DISTDIR)436 @cp -a AUTHORS CHANGELOG CMakeLists.txt COPYING DelphesEnv.sh LICENSE NOTICE README README_4LHCb Makefile MinBias.pileup configure cards classes converters display doc examples external modules python readers validation $(DISTDIR) 437 437 @find $(DISTDIR) -depth -name .\* -exec rm -rf {} \; 438 438 @tar -czf $(DISTTAR) $(DISTDIR) -
modules/DualReadoutCalorimeter.cc
r4e5ed02 r25291dd 645 645 //cout<<"Creating tower with Pt, Eta, Phi, Energy: "<<pt<<","<<eta<<","<<phi<<","<<neutralEnergy<<endl; 646 646 tower->Momentum.SetPtEtaPhiE(pt, eta, phi, neutralEnergy); 647 tower->Eem = neutralEnergy; 648 tower->Ehad = 0.0; 649 tower->PID = 22; 650 647 648 // if no hadronic energy, use ECAL resolution 649 if (fHCalTowerEnergy <= fHCalEnergyMin) 650 { 651 tower->Eem = neutralEnergy; 652 tower->Ehad = 0.0; 653 tower->PID = 22; 654 } 655 656 // if hadronic fraction > 0, use HCAL resolution 657 else 658 { 659 tower->Eem = 0; 660 tower->Ehad = neutralEnergy; 661 tower->PID = 130; 662 } 663 664 fEFlowPhotonOutputArray->Add(tower); 651 665 652 653 fEFlowPhotonOutputArray->Add(tower);654 655 666 656 667 //clone tracks -
modules/TrackCovariance.cc
r4e5ed02 r25291dd 102 102 while((candidate = static_cast<Candidate *>(fItInputArray->Next()))) 103 103 { 104 const TLorentzVector &candidatePosition = candidate->InitialPosition; 104 105 // converting to meters 106 const TLorentzVector &candidatePosition = candidate->InitialPosition*1e-03; 105 107 const TLorentzVector &candidateMomentum = candidate->Momentum; 106 108 … … 113 115 114 116 candidate->Momentum.SetVectM(track.GetObsP(), mass); 115 candidate->InitialPosition.SetXYZT(track.GetObsX().X(),track.GetObsX().Y(),track.GetObsX().Z(),candidatePosition.T());116 117 118 // converting back to mm 119 candidate->InitialPosition.SetXYZT(track.GetObsX().X()*1e03,track.GetObsX().Y()*1e03,track.GetObsX().Z()*1e03,candidatePosition.T()*1e03); 120 117 121 // save full covariance 5x5 matrix internally (D0, phi, Curvature, dz, ctg(theta)) 118 122 candidate->TrackCovariance = track.GetCov(); … … 123 127 ct = track.GetObsPar()[4]; 124 128 125 candidate->Xd = track.GetObsX().X() ;126 candidate->Yd = track.GetObsX().Y() ;127 candidate->Zd = track.GetObsX().Z() ;129 candidate->Xd = track.GetObsX().X()*1e03; 130 candidate->Yd = track.GetObsX().Y()*1e03; 131 candidate->Zd = track.GetObsX().Z()*1e03; 128 132 129 candidate->D0 = track.GetObsPar()[0];130 candidate->Phi = track.GetObsPar()[1];131 candidate->C = track.GetObsPar()[2];132 candidate->DZ = track.GetObsPar()[3];133 candidate->D0 = track.GetObsPar()[0]*1e03; 134 candidate->Phi = track.GetObsPar()[1]; 135 candidate->C = track.GetObsPar()[2]*1e03; 136 candidate->DZ = track.GetObsPar()[3]*1e03; 133 137 candidate->CtgTheta = track.GetObsPar()[4]; 134 candidate->P = track.GetObsP().Mag();135 candidate->PT = pt;136 candidate->Charge = q;138 candidate->P = track.GetObsP().Mag(); 139 candidate->PT = pt; 140 candidate->Charge = q; 137 141 138 dd0 = TMath::Sqrt(track.GetCov()(0, 0)) ;139 ddz = TMath::Sqrt(track.GetCov()(3, 3)) ;142 dd0 = TMath::Sqrt(track.GetCov()(0, 0))*1e03; 143 ddz = TMath::Sqrt(track.GetCov()(3, 3))*1e03; 140 144 dphi = TMath::Sqrt(track.GetCov()(1, 1)); 141 145 dct = TMath::Sqrt(track.GetCov()(4, 4)); 142 146 dpt = 2 * TMath::Sqrt( track.GetCov()(2, 2))*pt*pt / (0.2998*fBz); 143 147 dp = TMath::Sqrt((1.+ct*ct)*dpt*dpt + 4*pt*pt*ct*ct*dct*dct/(1.+ct*ct)/(1.+ct*ct)); 144 dC = TMath::Sqrt(track.GetCov()(2, 2)) ;148 dC = TMath::Sqrt(track.GetCov()(2, 2))*1e03; 145 149 146 150 candidate->ErrorD0 = dd0; -
readers/DelphesPythia8.cpp
r4e5ed02 r25291dd 296 296 297 297 // jet matching 298 #if PYTHIA_VERSION_INTEGER < 8300 298 299 matching = combined->getHook(*pythia); 299 300 if(!matching) … … 302 303 } 303 304 pythia->setUserHooksPtr(matching); 305 #endif 304 306 305 307 if(pythia == NULL) -
validation/DelphesValidation.cpp
r4e5ed02 r25291dd 32 32 #include "TGraph.h" 33 33 #include "TGraphErrors.h" 34 #include "TF1.h" 34 35 #include "TH1.h" 35 36 #include "TH2.h" -
validation/validation.sh
r4e5ed02 r25291dd 38 38 outputrootdir=report_${cardbase%.*}/root 39 39 cardlabel=${cardbase%.*} 40 version= $(cat VERSION)40 version=x.y.z 41 41 outpdf=$mainoutputdir/${output%.*}.pdf 42 42 cardsdir=validation/cards
Note:
See TracChangeset
for help on using the changeset viewer.