|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
|
#include <string>
#include <vector>
#include "GeneratorInterface/LHEInterface/plugins/LHEProvenanceHelper.h"
#include "GeneratorInterface/LHEInterface/interface/LHERunInfo.h"
#include "DataFormats/Provenance/interface/ProcessHistory.h"
#include "DataFormats/Provenance/interface/ProcessHistoryRegistry.h"
#include "DataFormats/Provenance/interface/ProductRegistry.h"
#include "DataFormats/Provenance/interface/BranchIDList.h"
#include "DataFormats/Provenance/interface/BranchIDListHelper.h"
#include "SimDataFormats/GeneratorProducts/interface/LesHouches.h"
#include "FWCore/AbstractServices/interface/ResourceInformation.h"
#include "FWCore/Utilities/interface/TypeID.h"
#include "FWCore/Reflection/interface/TypeWithDict.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Version/interface/GetReleaseVersion.h"
namespace edm {
LHEProvenanceHelper::LHEProvenanceHelper(TypeID const& eventProductType,
TypeID const& runProductType,
ProductRegistry& productRegistry,
BranchIDListHelper& branchIDListHelper)
: eventProductProductDescription_(ProductDescription(InEvent,
"source",
"LHEFile"
// , "LHE"
,
"LHEEventProduct",
"LHEEventProduct",
"",
TypeWithDict(eventProductType.typeInfo()),
false)),
runProductProductDescription_(ProductDescription(InRun,
"source",
"LHEFile"
// , "LHE"
,
"LHERunInfoProduct",
"LHERunInfoProduct",
"",
TypeWithDict(runProductType.typeInfo()),
false)),
eventProductProvenance_(eventProductProductDescription_.branchID()),
commonProcessParameterSet_(fillCommonProcessParameterSet()),
processParameterSet_(),
branchListIndexes_(1, 0) {
// Add the products to the product registry
auto ep = eventProductProductDescription_;
ep.setIsProvenanceSetOnRead();
productRegistry.copyProduct(ep);
auto rp = runProductProductDescription_;
rp.setIsProvenanceSetOnRead();
productRegistry.copyProduct(rp);
BranchIDList bli(1UL, ep.branchID().id());
branchIDListHelper.updateFromInput({{bli}});
}
ParameterSet LHEProvenanceHelper::fillCommonProcessParameterSet() {
// We create a process parameter set for the "LHC" process.
// This function only fills parameters whose values are independent of the LHE input files.
// We don't currently use the untracked parameters, However, we make them available, just in case.
ParameterSet pset;
std::string const& moduleLabel = eventProductProductDescription_.moduleLabel();
std::string const& processName = eventProductProductDescription_.processName();
std::string const moduleName = "LHESource";
typedef std::vector<std::string> vstring;
vstring empty;
vstring modlbl;
modlbl.reserve(1);
modlbl.push_back(moduleLabel);
pset.addParameter("@all_sources", modlbl);
ParameterSet triggerPaths;
triggerPaths.addParameter<vstring>("@trigger_paths", empty);
pset.addParameter<ParameterSet>("@trigger_paths", triggerPaths);
ParameterSet pseudoInput;
pseudoInput.addParameter<std::string>("@module_edm_type", "Source");
pseudoInput.addParameter<std::string>("@module_label", moduleLabel);
pseudoInput.addParameter<std::string>("@module_type", moduleName);
pset.addParameter<ParameterSet>(moduleLabel, pseudoInput);
pset.addParameter<vstring>("@all_esmodules", empty);
pset.addParameter<vstring>("@all_esprefers", empty);
pset.addParameter<vstring>("@all_essources", empty);
pset.addParameter<vstring>("@all_loopers", empty);
pset.addParameter<vstring>("@all_modules", empty);
pset.addParameter<vstring>("@end_paths", empty);
pset.addParameter<vstring>("@paths", empty);
pset.addParameter<std::string>("@process_name", processName);
return pset;
}
void LHEProvenanceHelper::lheAugment(lhef::LHERunInfo const* runInfo) {
processParameterSet_ = commonProcessParameterSet_;
if (runInfo == nullptr)
return;
typedef std::vector<std::string> vstring;
auto const& heprup = *runInfo->getHEPRUP();
processParameterSet_.addParameter<int>("IDBMUP1", heprup.IDBMUP.first);
processParameterSet_.addParameter<int>("IDBMUP2", heprup.IDBMUP.second);
processParameterSet_.addParameter<int>("EBMUP1", heprup.EBMUP.first);
processParameterSet_.addParameter<int>("EBMUP2", heprup.EBMUP.second);
processParameterSet_.addParameter<int>("PDFGUP1", heprup.PDFGUP.first);
processParameterSet_.addParameter<int>("PDFGUP2", heprup.PDFGUP.second);
processParameterSet_.addParameter<int>("PDFSUP1", heprup.PDFSUP.first);
processParameterSet_.addParameter<int>("PDFSUP2", heprup.PDFSUP.second);
processParameterSet_.addParameter<int>("IDWTUP", heprup.IDWTUP);
for (auto const& header : runInfo->getHeaders()) {
if (!LHERunInfoProduct::isTagComparedInMerge(header.tag())) {
continue;
}
processParameterSet_.addParameter<vstring>(header.tag(), header.lines());
}
}
ProcessHistoryID LHEProvenanceHelper::lheInit(ProcessHistoryRegistry& processHistoryRegistry) {
// Now we register the process parameter set.
processParameterSet_.registerIt();
//std::cerr << processParameterSet_.dump() << std::endl;
// Insert an entry for this process in the process history registry
ProcessHistory ph;
edm::Service<edm::ResourceInformation> resourceInformationService;
edm::HardwareResourcesDescription hwResources;
if (resourceInformationService.isAvailable()) {
hwResources = resourceInformationService->hardwareResourcesDescription();
}
ph.emplace_back(
eventProductProductDescription_.processName(), processParameterSet_.id(), getReleaseVersion(), hwResources);
processHistoryRegistry.registerProcessHistory(ph);
// Save the process history ID for use every event.
return ph.id();
}
} // namespace edm
|