Changes in / [c3c9ac5:f6b6ee7] in git
- Files:
-
- 1 added
- 10 deleted
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
CMakeLists.txt
rc3c9ac5 rf6b6ee7 7 7 8 8 # Declare ROOT dependency 9 find_package(ROOT COMPONENTS EG Eve Ge nVector Hist Physics Matrix Graf RIO Tree Gpad RGL MathCore)9 find_package(ROOT COMPONENTS EG Eve Geom Gui GuiHtml GenVector Hist Physics Matrix Graf RIO Tree Gpad RGL MathCore) 10 10 include(${ROOT_USE_FILE}) 11 11 -
Makefile
rc3c9ac5 rf6b6ee7 17 17 CXXFLAGS += $(ROOTCFLAGS) -Wno-write-strings -D_FILE_OFFSET_BITS=64 -DDROP_CGAL -I. -Iexternal -Iexternal/tcl 18 18 DELPHES_LIBS = $(shell $(RC) --libs) -lEG $(SYSLIBS) 19 DISPLAY_LIBS = $(shell $(RC) --evelibs) $(SYSLIBS)19 DISPLAY_LIBS = $(shell $(RC) --evelibs) -lGuiHtml $(SYSLIBS) 20 20 21 21 ifneq ($(CMSSW_FWLITE_INCLUDE_PATH),) 22 22 HAS_CMSSW = true 23 23 CXXFLAGS += -std=c++0x -I$(subst :, -I,$(CMSSW_FWLITE_INCLUDE_PATH)) 24 DELPHES_LIBS += -L$(subst include,lib,$(subst :, -L,$(CMSSW_FWLITE_INCLUDE_PATH)))24 OPT_LIBS += -L$(subst include,lib,$(subst :, -L,$(CMSSW_FWLITE_INCLUDE_PATH))) 25 25 ifneq ($(CMSSW_RELEASE_BASE),) 26 26 CXXFLAGS += -I$(CMSSW_RELEASE_BASE)/src 27 27 endif 28 28 ifneq ($(LD_LIBRARY_PATH),) 29 DELPHES_LIBS += -L$(subst include,lib,$(subst :, -L,$(LD_LIBRARY_PATH)))30 endif 31 DELPHES_LIBS += -lGenVector -lFWCoreFWLite -lDataFormatsFWLite -lDataFormatsPatCandidates -lDataFormatsLuminosity -lSimDataFormatsGeneratorProducts -lCommonToolsUtils29 OPT_LIBS += -L$(subst include,lib,$(subst :, -L,$(LD_LIBRARY_PATH))) 30 endif 31 OPT_LIBS += -lGenVector -lFWCoreFWLite -lDataFormatsFWLite -lDataFormatsPatCandidates -lDataFormatsLuminosity -lSimDataFormatsGeneratorProducts -lCommonToolsUtils 32 32 endif 33 33 34 34 ifneq ($(PROMC),) 35 35 HAS_PROMC = true 36 CXXFLAGS += -I$(PROMC)/include 37 DELPHES_LIBS += -L$(PROMC)/lib-lprotoc -lprotobuf -lprotobuf-lite -lcbook -lz36 CXXFLAGS += -I$(PROMC)/include -I$(PROMC)/src 37 OPT_LIBS += -L$(PROMC)/lib -lpromc -lprotoc -lprotobuf -lprotobuf-lite -lcbook -lz 38 38 endif 39 39 … … 41 41 HAS_PYTHIA8 = true 42 42 CXXFLAGS += -I$(PYTHIA8)/include 43 DELPHES_LIBS += -L$(PYTHIA8)/lib -lpythia8 -lLHAPDF -lgfortran -lz43 OPT_LIBS += -L$(PYTHIA8)/lib -lpythia8 -lLHAPDF -lgfortran -lz 44 44 else 45 45 ifneq ($(PYTHIA8DATA),) 46 46 HAS_PYTHIA8 = true 47 47 CXXFLAGS += -I$(PYTHIA8DATA)/../include 48 DELPHES_LIBS += -L$(PYTHIA8DATA)/../lib -lpythia8 -lLHAPDF -lgfortran -lz 49 endif 50 endif 48 OPT_LIBS += -L$(PYTHIA8DATA)/../lib -lpythia8 -lLHAPDF -lgfortran -lz 49 endif 50 endif 51 52 DELPHES_LIBS += $(OPT_LIBS) 53 DISPLAY_LIBS += $(OPT_LIBS) 51 54 52 55 ### … … 238 241 external/ExRootAnalysis/ExRootTreeWriter.h \ 239 242 external/ExRootAnalysis/ExRootTreeBranch.h \ 240 external/ExRootAnalysis/ExRootProgressBar.h \ 241 external/ProMC/ProMCBook.h 243 external/ExRootAnalysis/ExRootProgressBar.h 242 244 EXECUTABLE += \ 243 245 DelphesProMC$(ExeSuf) … … 245 247 EXECUTABLE_OBJ += \ 246 248 tmp/readers/DelphesProMC.$(ObjSuf) 247 248 tmp/external/ProMC/ProMCBook.$(ObjSuf): \249 external/ProMC/ProMCBook.$(SrcSuf)250 tmp/external/ProMC/ProMC.pb.$(ObjSuf): \251 external/ProMC/ProMC.pb.$(SrcSuf)252 tmp/external/ProMC/ProMCStat.pb.$(ObjSuf): \253 external/ProMC/ProMCStat.pb.$(SrcSuf)254 tmp/external/ProMC/ProMCHeader.pb.$(ObjSuf): \255 external/ProMC/ProMCHeader.pb.$(SrcSuf)256 tmp/external/ProMC/ProMCDescription.pb.$(ObjSuf): \257 external/ProMC/ProMCDescription.pb.$(SrcSuf)258 DELPHES_OBJ += \259 tmp/external/ProMC/ProMCBook.$(ObjSuf) \260 tmp/external/ProMC/ProMC.pb.$(ObjSuf) \261 tmp/external/ProMC/ProMCStat.pb.$(ObjSuf) \262 tmp/external/ProMC/ProMCHeader.pb.$(ObjSuf) \263 tmp/external/ProMC/ProMCDescription.pb.$(ObjSuf)264 265 ifeq ($(HAS_PYTHIA8),true)266 DELPHES_OBJ += \267 268 endif269 249 270 250 endif … … 381 361 display/DisplayLinkDef.h \ 382 362 display/DelphesDisplay.h \ 383 display/DelphesCaloData.h 363 display/DelphesCaloData.h \ 364 display/DelphesBranchElement.h \ 365 display/Delphes3DGeometry.h \ 366 display/DelphesEventDisplay.h \ 367 display/DelphesHtmlSummary.h \ 368 display/DelphesPlotSummary.h 384 369 DisplayDict$(PcmSuf): \ 385 370 tmp/display/DisplayDict$(PcmSuf) \ … … 1235 1220 display/DelphesDisplay.$(SrcSuf) \ 1236 1221 display/DelphesDisplay.h 1222 tmp/display/DelphesPlotSummary.$(ObjSuf): \ 1223 display/DelphesPlotSummary.$(SrcSuf) \ 1224 display/DelphesPlotSummary.h 1225 tmp/display/DelphesBranchElement.$(ObjSuf): \ 1226 display/DelphesBranchElement.$(SrcSuf) \ 1227 display/DelphesBranchElement.h \ 1228 classes/DelphesClasses.h 1229 tmp/display/DelphesEventDisplay.$(ObjSuf): \ 1230 display/DelphesEventDisplay.$(SrcSuf) \ 1231 external/ExRootAnalysis/ExRootConfReader.h \ 1232 external/ExRootAnalysis/ExRootTreeReader.h \ 1233 display/DelphesCaloData.h \ 1234 display/DelphesBranchElement.h \ 1235 display/Delphes3DGeometry.h \ 1236 display/DelphesEventDisplay.h \ 1237 classes/DelphesClasses.h 1238 tmp/display/DelphesHtmlSummary.$(ObjSuf): \ 1239 display/DelphesHtmlSummary.$(SrcSuf) \ 1240 display/DelphesHtmlSummary.h 1241 tmp/display/Delphes3DGeometry.$(ObjSuf): \ 1242 display/Delphes3DGeometry.$(SrcSuf) \ 1243 display/Delphes3DGeometry.h \ 1244 external/ExRootAnalysis/ExRootConfReader.h \ 1245 classes/DelphesClasses.h 1237 1246 tmp/display/DelphesCaloData.$(ObjSuf): \ 1238 1247 display/DelphesCaloData.$(SrcSuf) \ … … 1240 1249 DISPLAY_OBJ += \ 1241 1250 tmp/display/DelphesDisplay.$(ObjSuf) \ 1251 tmp/display/DelphesPlotSummary.$(ObjSuf) \ 1252 tmp/display/DelphesBranchElement.$(ObjSuf) \ 1253 tmp/display/DelphesEventDisplay.$(ObjSuf) \ 1254 tmp/display/DelphesHtmlSummary.$(ObjSuf) \ 1255 tmp/display/Delphes3DGeometry.$(ObjSuf) \ 1242 1256 tmp/display/DelphesCaloData.$(ObjSuf) 1243 1257 … … 1461 1475 @touch $@ 1462 1476 1477 display/DelphesBranchElement.h: \ 1478 display/DelphesCaloData.h 1479 @touch $@ 1480 1463 1481 external/fastjet/internal/DynamicNearestNeighbours.hh: \ 1464 1482 external/fastjet/internal/numconsts.hh … … 1500 1518 @touch $@ 1501 1519 1520 display/DelphesEventDisplay.h: \ 1521 external/ExRootAnalysis/ExRootTreeReader.h \ 1522 display/DelphesDisplay.h \ 1523 display/Delphes3DGeometry.h \ 1524 display/DelphesHtmlSummary.h \ 1525 display/DelphesPlotSummary.h 1526 @touch $@ 1527 1502 1528 modules/TauTagging.h: \ 1503 1529 classes/DelphesModule.h … … 1584 1610 modules/Hector.h: \ 1585 1611 classes/DelphesModule.h 1612 @touch $@ 1613 1614 display/DelphesPlotSummary.h: \ 1615 external/ExRootAnalysis/ExRootTreeReader.h 1586 1616 @touch $@ 1587 1617 -
doc/genMakefile.tcl
rc3c9ac5 rf6b6ee7 208 208 CXXFLAGS += $(ROOTCFLAGS) -Wno-write-strings -D_FILE_OFFSET_BITS=64 -DDROP_CGAL -I. -Iexternal -Iexternal/tcl 209 209 DELPHES_LIBS = $(shell $(RC) --libs) -lEG $(SYSLIBS) 210 DISPLAY_LIBS = $(shell $(RC) --evelibs) $(SYSLIBS)210 DISPLAY_LIBS = $(shell $(RC) --evelibs) -lGuiHtml $(SYSLIBS) 211 211 212 212 ifneq ($(CMSSW_FWLITE_INCLUDE_PATH),) 213 213 HAS_CMSSW = true 214 214 CXXFLAGS += -std=c++0x -I$(subst :, -I,$(CMSSW_FWLITE_INCLUDE_PATH)) 215 DELPHES_LIBS += -L$(subst include,lib,$(subst :, -L,$(CMSSW_FWLITE_INCLUDE_PATH)))215 OPT_LIBS += -L$(subst include,lib,$(subst :, -L,$(CMSSW_FWLITE_INCLUDE_PATH))) 216 216 ifneq ($(CMSSW_RELEASE_BASE),) 217 217 CXXFLAGS += -I$(CMSSW_RELEASE_BASE)/src 218 218 endif 219 219 ifneq ($(LD_LIBRARY_PATH),) 220 DELPHES_LIBS += -L$(subst include,lib,$(subst :, -L,$(LD_LIBRARY_PATH)))221 endif 222 DELPHES_LIBS += -lGenVector -lFWCoreFWLite -lDataFormatsFWLite -lDataFormatsPatCandidates -lDataFormatsLuminosity -lSimDataFormatsGeneratorProducts -lCommonToolsUtils220 OPT_LIBS += -L$(subst include,lib,$(subst :, -L,$(LD_LIBRARY_PATH))) 221 endif 222 OPT_LIBS += -lGenVector -lFWCoreFWLite -lDataFormatsFWLite -lDataFormatsPatCandidates -lDataFormatsLuminosity -lSimDataFormatsGeneratorProducts -lCommonToolsUtils 223 223 endif 224 224 225 225 ifneq ($(PROMC),) 226 226 HAS_PROMC = true 227 CXXFLAGS += -I$(PROMC)/include 228 DELPHES_LIBS += -L$(PROMC)/lib-lprotoc -lprotobuf -lprotobuf-lite -lcbook -lz227 CXXFLAGS += -I$(PROMC)/include -I$(PROMC)/src 228 OPT_LIBS += -L$(PROMC)/lib -lpromc -lprotoc -lprotobuf -lprotobuf-lite -lcbook -lz 229 229 endif 230 230 … … 232 232 HAS_PYTHIA8 = true 233 233 CXXFLAGS += -I$(PYTHIA8)/include 234 DELPHES_LIBS += -L$(PYTHIA8)/lib -lpythia8 -lLHAPDF -lgfortran -lz234 OPT_LIBS += -L$(PYTHIA8)/lib -lpythia8 -lLHAPDF -lgfortran -lz 235 235 else 236 236 ifneq ($(PYTHIA8DATA),) 237 237 HAS_PYTHIA8 = true 238 238 CXXFLAGS += -I$(PYTHIA8DATA)/../include 239 DELPHES_LIBS += -L$(PYTHIA8DATA)/../lib -lpythia8 -lLHAPDF -lgfortran -lz 240 endif 241 endif 239 OPT_LIBS += -L$(PYTHIA8DATA)/../lib -lpythia8 -lLHAPDF -lgfortran -lz 240 endif 241 endif 242 243 DELPHES_LIBS += $(OPT_LIBS) 244 DISPLAY_LIBS += $(OPT_LIBS) 242 245 243 246 ### … … 268 271 puts {ifeq ($(HAS_PROMC),true)} 269 272 executableDeps {readers/DelphesProMC.cpp} 270 sourceDeps {DELPHES} {external/ProMC/*.cc}271 273 puts {endif} 272 274 puts {} -
modules/ParticlePropagator.cc
rc3c9ac5 rf6b6ee7 158 158 if(TMath::Abs(q) < 1.0E-9 || TMath::Abs(fBz) < 1.0E-9) 159 159 { 160 // solve pt2*t^2 + 2*(px*x + py*y)*t +(fRadius2 - x*x - y*y) = 0160 // solve pt2*t^2 + 2*(px*x + py*y)*t - (fRadius2 - x*x - y*y) = 0 161 161 tmp = px*y - py*x; 162 162 discr2 = pt2*fRadius2 - tmp*tmp; -
modules/PileUpMergerPythia8.cc
rc3c9ac5 rf6b6ee7 132 132 dphi = gRandom->Uniform(-TMath::Pi(), TMath::Pi()); 133 133 134 for(i = 0; i < fPythia->event.size(); ++i)134 for(i = 1; i < fPythia->event.size(); ++i) 135 135 { 136 136 Pythia8::Particle &particle = fPythia->event[i]; -
readers/DelphesProMC.cpp
rc3c9ac5 rf6b6ee7 47 47 #include "ExRootAnalysis/ExRootProgressBar.h" 48 48 49 #include "ProMC /ProMC.pb.h"50 #include "ProMC /ProMCBook.h"51 #include "ProMC /ProMCHeader.pb.h"49 #include "ProMC.pb.h" 50 #include "ProMCBook.h" 51 #include "ProMCHeader.pb.h" 52 52 53 53 using namespace std; … … 55 55 //--------------------------------------------------------------------------- 56 56 57 void ConvertInput(ProMCEvent &event, ExRootTreeBranch *branch, DelphesFactory *factory, 58 TObjArray *allParticleOutputArray, TObjArray *stableParticleOutputArray, TObjArray *partonOutputArray, 59 TStopwatch *readStopWatch, TStopwatch *procStopWatch) 57 void ConvertInput(ProMCEvent &event, double momentumUnit, double positionUnit, 58 ExRootTreeBranch *branch, DelphesFactory *factory, 59 TObjArray *allParticleOutputArray, TObjArray *stableParticleOutputArray, 60 TObjArray *partonOutputArray, TStopwatch *readStopWatch, TStopwatch *procStopWatch) 60 61 { 61 62 Int_t i; … … 107 108 pid = mutableParticles->pdg_id(i); 108 109 status = mutableParticles->status(i); 109 px = mutableParticles->px(i); py = mutableParticles->py(i); pz = mutableParticles->pz(i); mass = mutableParticles->mass(i); 110 x = mutableParticles->x(i); y = mutableParticles->y(i); z = mutableParticles->z(i); t = mutableParticles->t(i); 110 111 px = mutableParticles->px(i)/momentumUnit; 112 py = mutableParticles->py(i)/momentumUnit; 113 pz = mutableParticles->pz(i)/momentumUnit; 114 mass = mutableParticles->mass(i)/momentumUnit; 115 x = mutableParticles->x(i)/positionUnit; 116 y = mutableParticles->y(i)/positionUnit; 117 z = mutableParticles->z(i)/positionUnit; 118 t = mutableParticles->t(i)/positionUnit; 111 119 112 120 candidate = factory->NewCandidate(); … … 172 180 Int_t i; 173 181 Long64_t eventCounter, numberOfEvents; 182 double momentumUnit = 1.0, positionUnit = 1.0; 174 183 175 184 if(argc < 4) … … 224 233 inputFile = new ProMCBook(argv[i], "r"); 225 234 235 ProMCHeader header = inputFile->getHeader(); 236 237 momentumUnit = static_cast<double>(header.momentumunit()); 238 positionUnit = static_cast<double>(header.lengthunit()); 239 240 241 226 242 if(inputFile == NULL) 227 243 { … … 248 264 249 265 procStopWatch.Start(); 250 ConvertInput(event, branchEvent, factory, 251 allParticleOutputArray, stableParticleOutputArray, partonOutputArray, 252 &readStopWatch, &procStopWatch); 266 ConvertInput(event, momentumUnit, positionUnit, 267 branchEvent, factory, 268 allParticleOutputArray, stableParticleOutputArray, 269 partonOutputArray, &readStopWatch, &procStopWatch); 253 270 modularDelphes->ProcessTask(); 254 271 procStopWatch.Stop(); -
readers/DelphesProMC.py
rc3c9ac5 rf6b6ee7 134 134 ################################################################################ 135 135 136 if len(sys.argv) < 2:136 if len(sys.argv) < 4: 137 137 print " Usage: DelphesProMC.py config_file output_file input_file(s)" 138 138 sys.exit(1) -
readers/DelphesPythia8.cpp
rc3c9ac5 rf6b6ee7 89 89 pdg = TDatabasePDG::Instance(); 90 90 91 for(i = 0; i < pythia->event.size(); ++i)91 for(i = 1; i < pythia->event.size(); ++i) 92 92 { 93 93 Pythia8::Particle &particle = pythia->event[i];
Note:
See TracChangeset
for help on using the changeset viewer.