Changeset 91acc76 in git
- Timestamp:
- Aug 23, 2016, 3:28:48 PM (8 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 19389b8, 7993cad
- Parents:
- e55f5b0 (diff), 98ce52a (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:
-
- 11 added
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
CMakeLists.txt
re55f5b0 r91acc76 41 41 $<TARGET_OBJECTS:tcl> 42 42 $<TARGET_OBJECTS:Hector> 43 $<TARGET_OBJECTS:PUPPI> 43 44 ) 44 45 -
Makefile
re55f5b0 r91acc76 324 324 modules/SimpleCalorimeter.h \ 325 325 modules/Calorimeter.h \ 326 modules/OldCalorimeter.h \ 326 327 modules/Isolation.h \ 327 328 modules/EnergyScale.h \ … … 732 733 external/ExRootAnalysis/ExRootFilter.h \ 733 734 external/ExRootAnalysis/ExRootClassifier.h 735 tmp/modules/OldCalorimeter.$(ObjSuf): \ 736 modules/OldCalorimeter.$(SrcSuf) \ 737 modules/OldCalorimeter.h \ 738 classes/DelphesClasses.h \ 739 classes/DelphesFactory.h \ 740 classes/DelphesFormula.h \ 741 external/ExRootAnalysis/ExRootResult.h \ 742 external/ExRootAnalysis/ExRootFilter.h \ 743 external/ExRootAnalysis/ExRootClassifier.h 734 744 tmp/modules/ParticlePropagator.$(ObjSuf): \ 735 745 modules/ParticlePropagator.$(SrcSuf) \ … … 982 992 tmp/modules/Merger.$(ObjSuf) \ 983 993 tmp/modules/MomentumSmearing.$(ObjSuf) \ 994 tmp/modules/OldCalorimeter.$(ObjSuf) \ 984 995 tmp/modules/ParticlePropagator.$(ObjSuf) \ 985 996 tmp/modules/PdgCodeFilter.$(ObjSuf) \ … … 1008 1019 1009 1020 tmp/external/PUPPI/PuppiAlgo.$(ObjSuf): \ 1010 external/PUPPI/PuppiAlgo.$(SrcSuf) \ 1011 external/fastjet/internal/base.hh 1021 external/PUPPI/PuppiAlgo.$(SrcSuf) 1012 1022 tmp/external/PUPPI/PuppiContainer.$(ObjSuf): \ 1013 1023 external/PUPPI/PuppiContainer.$(SrcSuf) \ 1014 external/fastjet/internal/base.hh \1015 1024 external/fastjet/Selector.hh 1016 1025 tmp/external/PUPPI/puppiCleanContainer.$(ObjSuf): \ … … 1370 1379 external/PUPPI/RecoObj2.hh \ 1371 1380 external/PUPPI/AlgoObj.hh \ 1381 external/PUPPI/PuppiContainer.hh \ 1382 external/fastjet/PseudoJet.hh \ 1372 1383 classes/DelphesClasses.h \ 1373 1384 classes/DelphesFactory.h \ … … 1761 1772 @touch $@ 1762 1773 1774 modules/OldCalorimeter.h: \ 1775 classes/DelphesModule.h 1776 @touch $@ 1777 1763 1778 external/fastjet/ClusterSequence1GhostPassiveArea.hh: \ 1764 1779 external/fastjet/PseudoJet.hh \ … … 1805 1820 1806 1821 modules/RunPUPPI.h: \ 1807 classes/DelphesModule.h \ 1808 external/PUPPI/PuppiContainer.hh 1822 classes/DelphesModule.h 1809 1823 @touch $@ 1810 1824 … … 1912 1926 external/fastjet/plugins/CDFCones/fastjet/CDFMidPointPlugin.hh: \ 1913 1927 external/fastjet/JetDefinition.hh 1928 @touch $@ 1929 1930 external/PUPPI/PuppiContainer.hh: \ 1931 external/fastjet/PseudoJet.hh 1914 1932 @touch $@ 1915 1933 -
external/CMakeLists.txt
re55f5b0 r91acc76 3 3 add_subdirectory(Hector) 4 4 add_subdirectory(tcl) 5 add_subdirectory(PUPPI) -
modules/ModulesLinkDef.h
re55f5b0 r91acc76 40 40 #include "modules/SimpleCalorimeter.h" 41 41 #include "modules/Calorimeter.h" 42 #include "modules/OldCalorimeter.h" 42 43 #include "modules/Isolation.h" 43 44 #include "modules/EnergyScale.h" … … 59 60 #include "modules/PdgCodeFilter.h" 60 61 #include "modules/BeamSpotFilter.h" 62 #include "modules/RecoPuFilter.h" 61 63 #include "modules/Cloner.h" 62 64 #include "modules/Weighter.h" … … 89 91 #pragma link C++ class SimpleCalorimeter+; 90 92 #pragma link C++ class Calorimeter+; 93 #pragma link C++ class OldCalorimeter+; 91 94 #pragma link C++ class Isolation+; 92 95 #pragma link C++ class EnergyScale+; … … 108 111 #pragma link C++ class PdgCodeFilter+; 109 112 #pragma link C++ class BeamSpotFilter+; 113 #pragma link C++ class RecoPuFilter+; 110 114 #pragma link C++ class Cloner+; 111 115 #pragma link C++ class Weighter+; -
modules/PhotonConversions.cc
re55f5b0 r91acc76 59 59 fItInputArray(0), fConversionMap(0), fDecayXsec(0) 60 60 { 61 fDecayXsec = new TF1 ;61 fDecayXsec = new TF1("decayXsec","1.0 - 4.0/3.0 * x * (1.0 - x)", 0.0, 1.0); 62 62 fConversionMap = new DelphesCylindricalFormula; 63 63 } … … 84 84 85 85 fConversionMap->Compile(GetString("ConversionMap", "0.0")); 86 87 #if ROOT_VERSION_CODE < ROOT_VERSION(6,04,00)88 fDecayXsec->Compile("1.0 - 4.0/3.0 * x * (1.0 - x)");89 #else90 fDecayXsec->GetFormula()->Compile("1.0 - 4.0/3.0 * x * (1.0 - x)");91 #endif92 fDecayXsec->SetRange(0.0, 1.0);93 86 94 87 // import array with output from filter/classifier module -
readers/DelphesCMSFWLite.cpp
re55f5b0 r91acc76 65 65 ExRootTreeBranch *branchEvent, ExRootTreeBranch *branchRwgt, 66 66 DelphesFactory *factory, TObjArray *allParticleOutputArray, 67 TObjArray *stableParticleOutputArray, TObjArray *partonOutputArray )67 TObjArray *stableParticleOutputArray, TObjArray *partonOutputArray, Bool_t firstEvent) 68 68 { 69 69 70 fwlite::Handle< GenEventInfoProduct > handleGenEventInfo; 70 71 71 fwlite::Handle< LHEEventProduct > handleLHEEvent; 72 73 72 fwlite::Handle< vector< reco::GenParticle > > handleParticle; 73 74 74 vector< reco::GenParticle >::const_iterator itParticle; 75 75 … … 78 78 79 79 handleGenEventInfo.getByLabel(event, "generator"); 80 handleLHEEvent.getByLabel(event, "externalLHEProducer"); 81 handleParticle.getByLabel(event, "genParticles"); 80 81 if (!((handleLHEEvent.getBranchNameFor(event, "source")).empty())) 82 { 83 handleLHEEvent.getByLabel(event, "source"); 84 } 85 else if (!((handleLHEEvent.getBranchNameFor(event, "externalLHEProducer")).empty())) 86 { 87 handleLHEEvent.getByLabel(event, "externalLHEProducer"); 88 } 89 else if (firstEvent) 90 { 91 std::cout<<"Wrong LHEEvent Label! Please, check the input file."<<std::endl; 92 } 93 94 if (!((handleParticle.getBranchNameFor(event, "genParticles")).empty())) 95 { 96 handleParticle.getByLabel(event, "genParticles"); 97 } 98 else if (!((handleParticle.getBranchNameFor(event, "prunedGenParticles")).empty())) 99 { 100 handleParticle.getByLabel(event, "prunedGenParticles"); 101 } 102 else 103 { 104 std::cout<<"Wrong GenParticle Label! Please, check the input file."<<std::endl; 105 exit(-1); 106 } 107 108 Bool_t foundLHE = !((handleLHEEvent.getBranchNameFor(event, "source")).empty()) || !((handleLHEEvent.getBranchNameFor(event, "externalLHEProducer")).empty()); 82 109 83 110 HepMCEvent *element; … … 91 118 Double_t px, py, pz, e, mass; 92 119 Double_t x, y, z; 93 94 const vector< gen::WeightsInfo > &vectorWeightsInfo = handleLHEEvent->weights();95 vector< gen::WeightsInfo >::const_iterator itWeightsInfo;96 120 97 121 element = static_cast<HepMCEvent *>(branchEvent->NewEntry()); … … 117 141 element->ProcTime = 0.0; 118 142 119 for(itWeightsInfo = vectorWeightsInfo.begin(); itWeightsInfo != vectorWeightsInfo.end(); ++itWeightsInfo) 120 { 121 weight = static_cast<Weight *>(branchRwgt->NewEntry()); 122 weight->Weight = itWeightsInfo->wgt; 143 144 if(foundLHE) 145 { 146 const vector< gen::WeightsInfo > &vectorWeightsInfo = handleLHEEvent->weights(); 147 vector< gen::WeightsInfo >::const_iterator itWeightsInfo; 148 149 for(itWeightsInfo = vectorWeightsInfo.begin(); itWeightsInfo != vectorWeightsInfo.end(); ++itWeightsInfo) 150 { 151 weight = static_cast<Weight *>(branchRwgt->NewEntry()); 152 weight->Weight = itWeightsInfo->wgt; 153 } 123 154 } 124 155 … … 207 238 Int_t i; 208 239 Long64_t eventCounter, numberOfEvents; 240 Bool_t firstEvent = kTRUE; 209 241 210 242 if(argc < 4) … … 281 313 modularDelphes->Clear(); 282 314 treeWriter->Clear(); 315 283 316 for(event.toBegin(); !event.atEnd() && !interrupted; ++event) 284 317 { 285 318 ConvertInput(event, eventCounter, branchEvent, branchRwgt, factory, 286 allParticleOutputArray, stableParticleOutputArray, partonOutputArray );319 allParticleOutputArray, stableParticleOutputArray, partonOutputArray, firstEvent); 287 320 modularDelphes->ProcessTask(); 321 322 firstEvent = kFALSE; 288 323 289 324 treeWriter->Fill();
Note:
See TracChangeset
for help on using the changeset viewer.