|
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
|
#include <iostream>
#include <typeinfo>
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/Utilities/interface/EDMException.h"
#include "FWCore/Utilities/interface/FriendlyName.h"
#include "GeneratorInterface/Core/interface/FortranInstance.h"
// make sure PDFSET is pulled in when linking against the
// archive lhapdf library.
extern "C" void pdfset_(void);
__attribute__((visibility("hidden"))) void dummy() { pdfset_(); }
// implementation for the Fortran callbacks from Pythia6/Herwig6
void gen::upinit_() { FortranInstance::getInstance<FortranInstance>()->upInit(); }
void gen::upevnt_() { FortranInstance::getInstance<FortranInstance>()->upEvnt(); }
void gen::upveto_(int *veto) { *veto = FortranInstance::getInstance<FortranInstance>()->upVeto(); }
// static FortranInstance members;
gen::FortranInstance *gen::FortranInstance::currentInstance = nullptr;
const std::string gen::FortranInstance::kFortranInstance = "FortranInstance";
// FortranInstance methods
gen::FortranInstance::~FortranInstance() noexcept(false) {
if (currentInstance == this) {
edm::LogWarning("ReentrancyProblem") << edm::friendlyname::friendlyName(typeid(*this).name())
<< " destroyed while it was the "
"current active instance."
<< std::endl;
currentInstance = nullptr;
}
}
// FortranInstance current instance tracking
void gen::FortranInstance::enter() {
// we should add a boost::mutex here if we care about being
// multithread-safe
if (currentInstance && currentInstance != this)
throw edm::Exception(edm::errors::LogicError) << edm::friendlyname::friendlyName(typeid(*this).name())
<< "::enter() called from a different "
"instance while an instance was already active."
<< std::endl;
if (!currentInstance && instanceNesting != 0)
throw edm::Exception(edm::errors::LogicError) << edm::friendlyname::friendlyName(typeid(*this).name())
<< "::enter() called on an empty "
"instance, but instance counter is nonzero."
<< std::endl;
currentInstance = this;
instanceNesting++;
}
void gen::FortranInstance::leave() {
if (!currentInstance)
throw edm::Exception(edm::errors::LogicError) << edm::friendlyname::friendlyName(typeid(*this).name())
<< "::leave() called without an "
"active instance."
<< std::endl;
else if (currentInstance != this)
throw edm::Exception(edm::errors::LogicError) << edm::friendlyname::friendlyName(typeid(*this).name())
<< "::leave() called from a "
"different instance."
<< std::endl;
else if (instanceNesting <= 0)
throw edm::Exception(edm::errors::LogicError) << edm::friendlyname::friendlyName(typeid(*this).name())
<< "::leave() called with a "
"nesting level of zero."
<< std::endl;
if (--instanceNesting == 0)
currentInstance = nullptr;
}
void gen::FortranInstance::throwMissingInstance() {
throw edm::Exception(edm::errors::LogicError) << "FortranInstance::getInstance() called from "
"a Fortran context, but no current instance "
"has been registered."
<< std::endl;
}
// Herwig callback stubs
void gen::FortranInstance::upInit() {
throw cms::Exception("UnimplementedCallback") << edm::friendlyname::friendlyName(typeid(*this).name())
<< "::upInit() stub called. "
"If user process needs to be generated, please derive "
"and implement the upInit() method."
<< std::endl;
}
void gen::FortranInstance::upEvnt() {
throw cms::Exception("UnimplementedCallback") << edm::friendlyname::friendlyName(typeid(*this).name())
<< "::upEvnt() stub called. "
"If user process needs to be generated, please derive "
"and implement the upEvnt() method."
<< std::endl;
}
bool gen::FortranInstance::upVeto() { return false; }
|