ATLAS Offline Software
Public Member Functions | Static Public Attributes | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
Photospp_i Class Reference

#include <Photospp_i.h>

Inheritance diagram for Photospp_i:
Collaboration diagram for Photospp_i:

Public Member Functions

 Photospp_i (const std::string &name, ISvcLocator *pSvcLocator)
 Standard Athena algorithm constructor. More...
 
StatusCode initialize ()
 Initialise the Photospp_i algorithm and required services. More...
 
StatusCode execute ()
 Run Photos on one event Will require a pre-existing HepMC event in Storegate. More...
 
void setupPhotos ()
 Set up the Photos class This may be called in the initialize method or directly before the generation of the first event. More...
 
friend double ::phoranc_ (int *idum)
 This external fortran function is the PHOTOS++ random number generator We make it a friend so it can access our AthRNGSvc. More...
 
virtual StatusCode sysInitialize () override
 Override sysInitialize. More...
 
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies. More...
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 

Static Public Attributes

static CLHEP::HepRandomEngine * p_rndmEngine = nullptr
 

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
 
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed. More...
 

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

DataObjIDColl m_extendedExtraObjects
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Features for derived classes to use internally

ServiceHandle< IAthRNGSvcm_rndmSvc {this, "RndmSvc", "AthRNGSvc"}
 
IntegerProperty m_dsid {this, "Dsid", 999999}
 
IntegerProperty m_randomSeed {this, "RandomSeed", 1234567, "Random seed for the built-in random engine"}
 Seed for random number engine. More...
 
StringProperty m_genEventKey {this, "MCEventKey", "GEN_EVENT"}
 The GenEvent StoreGate key - FIXME should be using Read/WriteHandles here. More...
 
BooleanProperty m_exponentiation {this, "ExponentiationMode", true}
 Whether to use exponentiation mode (default = yes) More...
 
BooleanProperty m_createHistory {this, "CreateHistory", false}
 Whether to create history entries (default = yes) More...
 
BooleanProperty m_stopCritical {this, "StopCriticalErrors", false}
 Whether to stop on critical error (default = no) More...
 
BooleanProperty m_delayInitialisation {this, "DelayInitialisation", false}
 Delay initialisation until just before first event execution (default = no) More...
 
BooleanProperty m_ZMECorrection {this, "ZMECorrection", false}
 Whether to apply ME correction to Z decays (default = no, until validated) More...
 
BooleanProperty m_WMECorrection {this, "WMECorrection", false}
 Whether to apply ME correction to W decays (default = no, until validated) More...
 
BooleanProperty m_photonSplitting {this, "PhotonSplitting", false}
 Whether to include photon splitting. More...
 
DoubleProperty m_infraRedCutOff {this, "InfraRedCutOff", -1.}
 
DoubleProperty m_maxWtInterference {this, "WtInterference", 3.}
 
DoubleProperty m_alphaQED {this, "AlphaQED", 0.00729735039}
 Value of alpha_QED. More...
 
void reseedRandomEngine (const std::string &streamName, const EventContext &ctx)
 
CLHEP::HepRandomEngine * getRandomEngine (const std::string &streamName, unsigned long int randomSeedOffset, const EventContext &ctx) const
 
CLHEP::HepRandomEngine * getRandomEngineDuringInitialize (const std::string &streamName, unsigned long int randomSeedOffset, unsigned int conditionsRun=1, unsigned int lbn=1) const
 

Detailed Description

Definition at line 16 of file Photospp_i.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ Photospp_i()

Photospp_i::Photospp_i ( const std::string &  name,
ISvcLocator *  pSvcLocator 
)

Standard Athena algorithm constructor.

Definition at line 40 of file Photospp_i.cxx.

40  :
41  AthAlgorithm(name, pSvcLocator) {
42 }

Member Function Documentation

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  hndl,
const SG::VarHandleKeyArrayType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  hndl,
const SG::VarHandleKeyType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  hndl,
const SG::VarHandleType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  t,
const SG::NotHandleType  
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleBase hndl,
const std::string &  doc,
const SG::VarHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKey hndl,
const std::string &  doc,
const SG::VarHandleKeyType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType  
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
148  }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ double ::phoranc_()

Photospp_i::double ::phoranc_ ( int *  idum)

This external fortran function is the PHOTOS++ random number generator We make it a friend so it can access our AthRNGSvc.

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( ) const
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

StatusCode Photospp_i::execute ( )

Run Photos on one event Will require a pre-existing HepMC event in Storegate.

Definition at line 135 of file Photospp_i.cxx.

135  {
136 
137  // initialise if not done already
138 
140  setupPhotos();
141  m_delayInitialisation = false;
142  }
143 
144  //Re-seed the random number stream
145  const EventContext& ctx = Gaudi::Hive::currentContext();
146  reseedRandomEngine("PHOTOSPP", ctx);
147 
148  // Get the event collection
149  // FIXME should be using Read/WriteHandles here
150  McEventCollection* eventCollection = nullptr;
151  StatusCode sc = evtStore()->retrieve(eventCollection, m_genEventKey);
152  if (sc.isFailure() || eventCollection == 0) {
153  ATH_MSG_ERROR("Unable to retrieve event collection from StoreGate with key " << m_genEventKey);
154  return StatusCode::FAILURE;
155  } else {
156  ATH_MSG_DEBUG("Retrieved event collection from StoreGate with key " << m_genEventKey);
157  }
158 
159  // Get the event most recently added to the event collection
160  HepMC::GenEvent *event = eventCollection->back();
161  if(event == 0) {
162  ATH_MSG_ERROR("Photospp_i received a null HepMC event");
163  return StatusCode::FAILURE;
164  }
165  switch(event->momentum_unit()) {
166 #ifdef HEPMC3
168  Photos::setMomentumUnit(Photos::GEV);
169  Photos::forceMass(11, ParticleConstants::electronMassInMeV/1000.); // Convert MeV to GeV
170  break;
171  case HepMC3::Units::MomentumUnit::MEV:
172  Photos::setMomentumUnit(Photos::MEV);
173  Photos::forceMass(11, ParticleConstants::electronMassInMeV);
174  break;
175 #else
176  case HepMC::Units::GEV:
177  Photos::setMomentumUnit(Photos::GEV);
178  Photos::forceMass(11, ParticleConstants::electronMassInMeV/1000.); // Convert MeV to GeV
179  break;
180  case HepMC::Units::MEV:
181  Photos::setMomentumUnit(Photos::MEV);
182  Photos::forceMass(11, ParticleConstants::electronMassInMeV);
183  break;
184 #endif
185  default:
186  ATH_MSG_ERROR("Photospp_i received a event with unknown units.");
187  Photos::setMomentumUnit(Photos::DEFAULT_MOMENTUM);
188  break;
189  };
190  PhotosHepMCEvent photosEvent(event);
191  photosEvent.process();
192 
193  return StatusCode::SUCCESS;
194 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

51 {
52  // If we didn't find any symlinks to add, just return the collection
53  // from the base class. Otherwise, return the extended collection.
54  if (!m_extendedExtraObjects.empty()) {
56  }
57  return Algorithm::extraOutputDeps();
58 }

◆ getRandomEngine()

CLHEP::HepRandomEngine * Photospp_i::getRandomEngine ( const std::string &  streamName,
unsigned long int  randomSeedOffset,
const EventContext &  ctx 
) const
private

Definition at line 109 of file Photospp_i.cxx.

111 {
112  ATHRNG::RNGWrapper* rngWrapper = m_rndmSvc->getEngine(this, streamName);
113  rngWrapper->setSeed( streamName, ctx.slot(), randomSeedOffset, ctx.eventID().run_number() );
114  return rngWrapper->getEngine(ctx);
115 }

◆ getRandomEngineDuringInitialize()

CLHEP::HepRandomEngine * Photospp_i::getRandomEngineDuringInitialize ( const std::string &  streamName,
unsigned long int  randomSeedOffset,
unsigned int  conditionsRun = 1,
unsigned int  lbn = 1 
) const
private

Definition at line 118 of file Photospp_i.cxx.

119 {
120  const size_t slot=0;
121  EventContext ctx;
122  ctx.setSlot( slot );
123  ctx.setEventID (EventIDBase (conditionsRun,
124  EventIDBase::UNDEFEVT, // event
125  EventIDBase::UNDEFNUM, // timestamp
126  EventIDBase::UNDEFNUM, // timestamp ns
127  lbn));
129  Atlas::ExtendedEventContext( evtStore()->hiveProxyDict(),
130  conditionsRun) );
131  return getRandomEngine(streamName, randomSeedOffset, ctx);
132 }

◆ initialize()

StatusCode Photospp_i::initialize ( )

Initialise the Photospp_i algorithm and required services.

Definition at line 45 of file Photospp_i.cxx.

45  {
46  ATH_MSG_DEBUG("Photospp_i initializing");
47 
48  ATH_CHECK(m_rndmSvc.retrieve());
50 
52 
53  return StatusCode::SUCCESS;
54 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ msg() [1/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24  {
25  return this->msgStream();
26  }

◆ msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( const MSG::Level  lvl) const
inlineinherited

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ msgLvl()

bool AthCommonMsg< Algorithm >::msgLvl ( const MSG::Level  lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< Algorithm > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::renounceArray ( SG::VarHandleKeyArray handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ reseedRandomEngine()

void Photospp_i::reseedRandomEngine ( const std::string &  streamName,
const EventContext &  ctx 
)
private

Definition at line 101 of file Photospp_i.cxx.

102 {
103  long seeds[7];
104  ATHRNG::calculateSeedsMC21(seeds, streamName, ctx.eventID().event_number(), m_dsid, m_randomSeed);
105  p_rndmEngine->setSeeds(seeds, 0); // NOT THREAD-SAFE
106 }

◆ setupPhotos()

void Photospp_i::setupPhotos ( )

Set up the Photos class This may be called in the initialize method or directly before the generation of the first event.

Definition at line 56 of file Photospp_i.cxx.

56  {
57 
59  Photos::setAlphaQED(m_alphaQED);
60  Photos::setInterference(true);
61  Photos::setCorrectionWtForW(true);
62  Photos::maxWtInterference(m_maxWtInterference);
63  Photos::setMeCorrectionWtForW(m_WMECorrection);
64  Photos::setMeCorrectionWtForZ(m_ZMECorrection);
65  Photos::setPairEmission(m_photonSplitting);
66  Photos::forceMassFrom4Vector(true);
67  Photos::forceMassFromEventRecord(13);
68  Photos::forceMassFromEventRecord(15);
69  Photos::forceMass(11, ParticleConstants::electronMassInMeV); // The assumption that units are MEV will be checked later
70  Photos::forceMassFromEventRecord(211);
71  Photos::setTopProcessRadiation(false);
72 #ifdef HEPMC3
73  Photos::createHistoryEntries(m_createHistory, 0);
74 #else
75  Photos::createHistoryEntries(m_createHistory, 3);
76 #endif
77 
78  if(m_exponentiation) {
79  Photos::setExponentiation(true);
80  } else {
81  Photos::setInfraredCutOff(0.01);
82  Photos::setDoubleBrem(true);
83  Photos::setQuatroBrem(false);
84  phokey.iexp = 0;
85  }
86 
87  // over-ride default IR cutoff if user has set a specific value
88  if(m_infraRedCutOff > 0.) {
89  Photos::setInfraredCutOff(m_infraRedCutOff);
90  }
91 
92  Photospp::Log::LogWarning(false);
93  Photospp::Photos::setStopAtCriticalError(m_stopCritical);
94 
95  Photos::iniInfo();
96 
97  return;
98 }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

66  {
68 
69  if (sc.isFailure()) {
70  return sc;
71  }
72  ServiceHandle<ICondSvc> cs("CondSvc",name());
73  for (auto h : outputHandles()) {
74  if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
75  // do this inside the loop so we don't create the CondSvc until needed
76  if ( cs.retrieve().isFailure() ) {
77  ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
78  return StatusCode::SUCCESS;
79  }
80  if (cs->regHandle(this,*h).isFailure()) {
81  sc = StatusCode::FAILURE;
82  ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
83  << " with CondSvc");
84  }
85  }
86  }
87  return sc;
88 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase &  )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308  {
309  // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310  // << " size: " << m_vhka.size() << endmsg;
311  for (auto &a : m_vhka) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ m_alphaQED

DoubleProperty Photospp_i::m_alphaQED {this, "AlphaQED", 0.00729735039}
private

Value of alpha_QED.

Definition at line 89 of file Photospp_i.h.

◆ m_createHistory

BooleanProperty Photospp_i::m_createHistory {this, "CreateHistory", false}
private

Whether to create history entries (default = yes)

Definition at line 66 of file Photospp_i.h.

◆ m_delayInitialisation

BooleanProperty Photospp_i::m_delayInitialisation {this, "DelayInitialisation", false}
private

Delay initialisation until just before first event execution (default = no)

Definition at line 72 of file Photospp_i.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_dsid

IntegerProperty Photospp_i::m_dsid {this, "Dsid", 999999}
private

Definition at line 54 of file Photospp_i.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_exponentiation

BooleanProperty Photospp_i::m_exponentiation {this, "ExponentiationMode", true}
private

Whether to use exponentiation mode (default = yes)

Definition at line 63 of file Photospp_i.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_genEventKey

StringProperty Photospp_i::m_genEventKey {this, "MCEventKey", "GEN_EVENT"}
private

The GenEvent StoreGate key - FIXME should be using Read/WriteHandles here.

Definition at line 60 of file Photospp_i.h.

◆ m_infraRedCutOff

DoubleProperty Photospp_i::m_infraRedCutOff {this, "InfraRedCutOff", -1.}
private

Definition at line 84 of file Photospp_i.h.

◆ m_maxWtInterference

DoubleProperty Photospp_i::m_maxWtInterference {this, "WtInterference", 3.}
private

Definition at line 86 of file Photospp_i.h.

◆ m_photonSplitting

BooleanProperty Photospp_i::m_photonSplitting {this, "PhotonSplitting", false}
private

Whether to include photon splitting.

Definition at line 81 of file Photospp_i.h.

◆ m_randomSeed

IntegerProperty Photospp_i::m_randomSeed {this, "RandomSeed", 1234567, "Random seed for the built-in random engine"}
private

Seed for random number engine.

Definition at line 57 of file Photospp_i.h.

◆ m_rndmSvc

ServiceHandle<IAthRNGSvc> Photospp_i::m_rndmSvc {this, "RndmSvc", "AthRNGSvc"}
private

Definition at line 51 of file Photospp_i.h.

◆ m_stopCritical

BooleanProperty Photospp_i::m_stopCritical {this, "StopCriticalErrors", false}
private

Whether to stop on critical error (default = no)

Definition at line 69 of file Photospp_i.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_WMECorrection

BooleanProperty Photospp_i::m_WMECorrection {this, "WMECorrection", false}
private

Whether to apply ME correction to W decays (default = no, until validated)

Definition at line 78 of file Photospp_i.h.

◆ m_ZMECorrection

BooleanProperty Photospp_i::m_ZMECorrection {this, "ZMECorrection", false}
private

Whether to apply ME correction to Z decays (default = no, until validated)

Definition at line 75 of file Photospp_i.h.

◆ p_rndmEngine

CLHEP::HepRandomEngine * Photospp_i::p_rndmEngine = nullptr
static

Definition at line 39 of file Photospp_i.h.


The documentation for this class was generated from the following files:
Photospp_i::getRandomEngine
CLHEP::HepRandomEngine * getRandomEngine(const std::string &streamName, unsigned long int randomSeedOffset, const EventContext &ctx) const
Definition: Photospp_i.cxx:109
ATHRNG::RNGWrapper::setSeed
void setSeed(const std::string &algName, const EventContext &ctx)
Set the random seed using a string (e.g.
Definition: RNGWrapper.h:169
Photospp_i::m_photonSplitting
BooleanProperty m_photonSplitting
Whether to include photon splitting.
Definition: Photospp_i.h:81
Photospp_i::getRandomEngineDuringInitialize
CLHEP::HepRandomEngine * getRandomEngineDuringInitialize(const std::string &streamName, unsigned long int randomSeedOffset, unsigned int conditionsRun=1, unsigned int lbn=1) const
Definition: Photospp_i.cxx:118
initialize
void initialize()
Definition: run_EoverP.cxx:894
Photospp_i::m_exponentiation
BooleanProperty m_exponentiation
Whether to use exponentiation mode (default = yes)
Definition: Photospp_i.h:63
McEventCollection
McEventCollection
Definition: GeneratorObjectsTPCnv.cxx:60
Photospp_i::m_WMECorrection
BooleanProperty m_WMECorrection
Whether to apply ME correction to W decays (default = no, until validated)
Definition: Photospp_i.h:78
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
Photospp_i::m_rndmSvc
ServiceHandle< IAthRNGSvc > m_rndmSvc
Definition: Photospp_i.h:51
Photospp_i::m_infraRedCutOff
DoubleProperty m_infraRedCutOff
Definition: Photospp_i.h:84
Photospp_i::m_createHistory
BooleanProperty m_createHistory
Whether to create history entries (default = yes)
Definition: Photospp_i.h:66
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
AthCommonDataStore
Definition: AthCommonDataStore.h:52
AthAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthAlgorithm.cxx:66
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
event
POOL::TEvent event(POOL::TEvent::kClassAccess)
AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
Atlas::ExtendedEventContext
Definition: ExtendedEventContext.h:23
Photospp_i::m_dsid
IntegerProperty m_dsid
Definition: Photospp_i.h:54
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
GEV
#define GEV
Definition: PrintPhotonSF.cxx:25
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
Photospp_i::reseedRandomEngine
void reseedRandomEngine(const std::string &streamName, const EventContext &ctx)
Definition: Photospp_i.cxx:101
WriteHiveWithMetaData.streamName
string streamName
Definition: WriteHiveWithMetaData.py:22
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
ATHRNG::RNGWrapper
A wrapper class for event-slot-local random engines.
Definition: RNGWrapper.h:56
Photospp_i::setupPhotos
void setupPhotos()
Set up the Photos class This may be called in the initialize method or directly before the generation...
Definition: Photospp_i.cxx:56
Photospp_i::m_randomSeed
IntegerProperty m_randomSeed
Seed for random number engine.
Definition: Photospp_i.h:57
ATHRNG::RNGWrapper::getEngine
CLHEP::HepRandomEngine * getEngine(const EventContext &ctx) const
Retrieve the random engine corresponding to the provided EventContext.
Definition: RNGWrapper.h:134
Photospp_i::m_ZMECorrection
BooleanProperty m_ZMECorrection
Whether to apply ME correction to Z decays (default = no, until validated)
Definition: Photospp_i.h:75
AthAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Definition: AthAlgorithm.h:79
a
TList * a
Definition: liststreamerinfos.cxx:10
h
Photospp_i::p_rndmEngine
static CLHEP::HepRandomEngine * p_rndmEngine
Definition: Photospp_i.h:39
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
Photospp_i::m_genEventKey
StringProperty m_genEventKey
The GenEvent StoreGate key - FIXME should be using Read/WriteHandles here.
Definition: Photospp_i.h:60
Photospp_i::m_alphaQED
DoubleProperty m_alphaQED
Value of alpha_QED.
Definition: Photospp_i.h:89
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:629
Photospp_i::m_delayInitialisation
BooleanProperty m_delayInitialisation
Delay initialisation until just before first event execution (default = no)
Definition: Photospp_i.h:72
ATHRNG::calculateSeedsMC21
void calculateSeedsMC21(long *seeds, const std::string &algName, uint64_t ev, uint64_t run, uint64_t offset=0)
Set the random seed using a string (e.g.
Definition: RNGWrapper.cxx:37
AthAlgorithm::AthAlgorithm
AthAlgorithm()
Default constructor:
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:801
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
Photospp_i::m_stopCritical
BooleanProperty m_stopCritical
Whether to stop on critical error (default = no)
Definition: Photospp_i.h:69
Photospp_i::m_maxWtInterference
DoubleProperty m_maxWtInterference
Definition: Photospp_i.h:86
Atlas::setExtendedEventContext
void setExtendedEventContext(EventContext &ctx, ExtendedEventContext &&ectx)
Move an extended context into a context object.
Definition: ExtendedEventContext.cxx:50
fitman.k
k
Definition: fitman.py:528
ServiceHandle< ICondSvc >
ParticleConstants::PDG2011::electronMassInMeV
constexpr double electronMassInMeV
the mass of the electron (in MeV)
Definition: ParticleConstants.h:26