|
ATLAS Offline Software
|
#include <Egamma1_LArStrip_Fex.h>
|
| Egamma1_LArStrip_Fex (const std::string &name, ISvcLocator *pSvcLocator) |
|
virtual StatusCode | initialize () override |
|
virtual StatusCode | execute (const EventContext &) const override |
|
virtual StatusCode | sysInitialize () override |
| Override sysInitialize. More...
|
|
virtual bool | isClonable () const override |
| Specify if the algorithm is clonable. More...
|
|
virtual unsigned int | cardinality () const override |
| Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...
|
|
virtual StatusCode | sysExecute (const EventContext &ctx) override |
| Execute an algorithm. More...
|
|
virtual const DataObjIDColl & | extraOutputDeps () const override |
| Return the list of extra output dependencies. More...
|
|
virtual bool | filterPassed (const EventContext &ctx) const |
|
virtual void | setFilterPassed (bool state, const EventContext &ctx) const |
|
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 > &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 |
|
|
StatusCode | findNeighborhoods (const std::vector< const xAOD::eFexEMRoI * > &, const std::vector< const CaloCell * > &, LArStripNeighborhoodContainer &) const |
|
StatusCode | findNeighborhood (const xAOD::eFexEMRoI *, const std::vector< const CaloCell * > &, LArStripNeighborhoodContainer &) const |
|
StatusCode | findClosestCellToRoI (const xAOD::eFexEMRoI *, const std::vector< const CaloCell * > &, const CaloCell *&) const |
|
StatusCode | dump (const xAOD::EventInfo &eventInfo, const LArStripNeighborhoodContainer &) const |
|
StatusCode | dumpTerse (const xAOD::EventInfo &eventInfo, const LArStripNeighborhoodContainer &) const |
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &) |
| specialization for handling Gaudi::Property<SG::VarHandleKey> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &) |
| specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &) |
| specialization for handling Gaudi::Property<SG::VarHandleBase> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &) |
| specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
|
|
Definition at line 26 of file Egamma1_LArStrip_Fex.h.
◆ StoreGateSvc_t
◆ Egamma1_LArStrip_Fex()
GlobalSim::Egamma1_LArStrip_Fex::Egamma1_LArStrip_Fex |
( |
const std::string & |
name, |
|
|
ISvcLocator * |
pSvcLocator |
|
) |
| |
◆ cardinality()
unsigned int AthReentrantAlgorithm::cardinality |
( |
| ) |
const |
|
overridevirtualinherited |
Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
Override this to return 0 for reentrant algorithms.
Override this to return 0 for reentrant algorithms.
Definition at line 55 of file AthReentrantAlgorithm.cxx.
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation 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.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation 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.
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation 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.
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation 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.
◆ declareProperty() [6/6]
◆ detStore()
◆ dump()
Definition at line 290 of file Egamma1_LArStrip_Fex.cxx.
293 std::ofstream
out(
name() +
"_" +
299 <<
" is simulation " << std::boolalpha
303 for (
const auto& nbhd : neighborhoods) {
304 out << *nbhd <<
'\n';
309 return StatusCode::SUCCESS;
◆ dumpTerse()
Definition at line 340 of file Egamma1_LArStrip_Fex.cxx.
343 std::ofstream
out(
name() +
"_" +
348 <<
" is simulation " << std::boolalpha
352 for (
const auto&
n : neighborhoods) {
dump_n(
n,
out);}
357 return StatusCode::SUCCESS;
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode GlobalSim::Egamma1_LArStrip_Fex::execute |
( |
const EventContext & |
ctx | ) |
const |
|
overridevirtual |
Definition at line 59 of file Egamma1_LArStrip_Fex.cxx.
66 if(!eventInfo.isValid()) {
68 return StatusCode::FAILURE;
71 std::vector<const CaloCell*>
cells;
75 std::vector<const xAOD::eFexEMRoI*>
rois;
83 auto neighborhoods = std::make_unique<LArStripNeighborhoodContainer>();
98 CHECK(h_write.record(std::move(neighborhoods)));
102 return StatusCode::SUCCESS;
◆ extraDeps_update_handler()
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 & AthReentrantAlgorithm::extraOutputDeps |
( |
| ) |
const |
|
overridevirtualinherited |
Return the list of extra output dependencies.
This list is extended to include symlinks implied by inheritance relations.
Definition at line 79 of file AthReentrantAlgorithm.cxx.
86 return Algorithm::extraOutputDeps();
◆ filterPassed()
virtual bool AthReentrantAlgorithm::filterPassed |
( |
const EventContext & |
ctx | ) |
const |
|
inlinevirtualinherited |
◆ findClosestCellToRoI()
◆ findNeighborhood()
Definition at line 119 of file Egamma1_LArStrip_Fex.cxx.
153 auto cells_near_roi = std::vector<const CaloCell*>();
156 auto phi_ind = [](
const auto&
c) {
158 std::size_t iphi = 32 +
int(std::floor(
c->phi()/dphi));
160 std::optional<std::size_t>() : std::make_optional(iphi);
164 auto roi_phi_index = *phi_ind(roi);
167 auto roi_phi_indices = *
wrap3(roi_phi_index);
168 for (
const auto&
i : roi_phi_indices) {
173 auto close = std::vector<std::vector<const CaloCell*>>(3);
174 for (
auto&
v :close) {
v.reserve(100);}
177 constexpr
double half_deta_roi{0.05};
178 constexpr
double half_deta_neigh{8.5*0.0031};
179 constexpr
double half_deta_fid{half_deta_roi + half_deta_neigh};
182 double etalim_low =
std::max(roi->
eta()-half_deta_fid, -1.4);
183 double etalim_high =
std::min(roi->
eta()+half_deta_fid, 1.4);
190 for(
const auto& iroi : roi_phi_indices) {
191 auto c_eta =
cell->eta();
192 if (iroi ==
icell and c_eta >= etalim_low and c_eta < etalim_high) {
201 for (
const auto&
v : close) {
ATH_MSG_DEBUG(
"size close " <<
v.size());}
205 auto roi_cells = std::vector<const CaloCell*> ();
206 roi_cells.reserve(close[1].
size());
209 std::back_inserter(roi_cells),
210 [&roi](
const auto&
c) {
211 return std::abs(
c->eta() - roi->
eta()) < half_deta_roi;
217 [](
const auto&
l,
const auto&
r) {
218 return l->e() <
r->e();
223 return StatusCode::FAILURE;
232 auto neigh_cells = std::vector<std::vector<const CaloCell*>>(3);
235 const auto max_cell_eta = max_cell->eta();
238 for (std::size_t iv{0
ul}; iv != close.size(); ++iv) {
241 std::back_inserter(neigh_cells[iv]),
242 [&max_cell_eta, &half_deta_neigh](
const auto&
c){
243 return abs(
c->eta()-max_cell_eta) < half_deta_neigh;
250 if (max_neigh_cell_it ==
std::end(neigh_cells[1])){
252 return StatusCode::FAILURE;
258 auto toStripData = [](
const auto& fromCells){
259 auto stripdata = std::vector<StripData>();
260 stripdata.reserve(fromCells.size());
263 back_inserter(stripdata),
265 return StripData(
c->eta(),
271 auto low = toStripData(neigh_cells[0]);
272 auto center = toStripData(neigh_cells[1]);
273 auto high = toStripData(neigh_cells[2]);
276 Coords cell_c{max_cell->eta(), max_cell->phi()};
279 neighborhoods.push_back(std::make_unique<LArStripNeighborhood>(low,
284 max_neigh_cell_pos));
286 return StatusCode::SUCCESS;
◆ findNeighborhoods()
◆ initialize()
StatusCode GlobalSim::Egamma1_LArStrip_Fex::initialize |
( |
| ) |
|
|
overridevirtual |
◆ inputHandles()
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.
◆ isClonable()
bool AthReentrantAlgorithm::isClonable |
( |
| ) |
const |
|
overridevirtualinherited |
Specify if the algorithm is clonable.
Reentrant algorithms are clonable.
Reimplemented in Simulation::BeamEffectsAlg, InDet::SiTrackerSpacePointFinder, InDet::SCT_Clusterization, InDet::SiSPSeededTrackFinder, SCTRawDataProvider, InDet::GNNSeedingTrackMaker, SCT_PrepDataToxAOD, RoIBResultToxAOD, InDet::SiSPGNNTrackMaker, SCT_CablingCondAlgFromCoraCool, SCT_ReadCalibDataTestAlg, SCT_CablingCondAlgFromText, SCT_ReadCalibChipDataTestAlg, SCT_TestCablingAlg, SCT_ConfigurationConditionsTestAlg, ITkPixelCablingAlg, ITkStripCablingAlg, SCTEventFlagWriter, SCT_ConditionsSummaryTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_TdaqEnabledTestAlg, SCT_SiliconConditionsTestAlg, SCTSiLorentzAngleTestAlg, SCT_ByteStreamErrorsTestAlg, SCT_ConditionsParameterTestAlg, SCT_FlaggedConditionTestAlg, SCT_StripVetoTestAlg, SCT_RawDataToxAOD, and SCTSiPropertiesTestAlg.
Definition at line 44 of file AthReentrantAlgorithm.cxx.
◆ msg() [1/2]
◆ msg() [2/2]
◆ msgLvl()
◆ outputHandles()
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()
◆ renounceArray()
◆ setFilterPassed()
virtual void AthReentrantAlgorithm::setFilterPassed |
( |
bool |
state, |
|
|
const EventContext & |
ctx |
|
) |
| const |
|
inlinevirtualinherited |
◆ sysExecute()
StatusCode AthReentrantAlgorithm::sysExecute |
( |
const EventContext & |
ctx | ) |
|
|
overridevirtualinherited |
Execute an algorithm.
We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.
Definition at line 67 of file AthReentrantAlgorithm.cxx.
69 return Gaudi::Algorithm::sysExecute (ctx);
◆ sysInitialize()
StatusCode AthReentrantAlgorithm::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< Gaudi::Algorithm > >.
Reimplemented in InputMakerBase, and HypoBase.
Definition at line 96 of file AthReentrantAlgorithm.cxx.
107 if ( cs.retrieve().isFailure() ) {
109 return StatusCode::SUCCESS;
111 if (cs->regHandle(
this,*
h).isFailure()) {
112 sc = StatusCode::FAILURE;
113 ATH_MSG_ERROR(
"unable to register WriteCondHandle " <<
h->fullKey()
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ updateVHKA()
◆ m_cellProducer
Initial value:{this,
"caloCellProducer",
"EMB1CellFromCaloCells",
"AlgTool to provide a vector of CaloCells"
}
Definition at line 45 of file Egamma1_LArStrip_Fex.h.
◆ m_detStore
◆ m_dump
Gaudi::Property<bool> GlobalSim::Egamma1_LArStrip_Fex::m_dump |
|
private |
◆ m_dumpTerse
Gaudi::Property<bool> GlobalSim::Egamma1_LArStrip_Fex::m_dumpTerse |
|
private |
Initial value:{
this,
"dumpTerse",
false,
"flag to enable terse dumps"}
Definition at line 65 of file Egamma1_LArStrip_Fex.h.
◆ m_eventInfoKey
◆ m_evtStore
◆ m_extendedExtraObjects
DataObjIDColl AthReentrantAlgorithm::m_extendedExtraObjects |
|
privateinherited |
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
Empty if no symlinks were found.
Definition at line 153 of file AthReentrantAlgorithm.h.
◆ m_neighKey
Initial value:{
this,
"stripNeighborhoodKey",
"stripNeighborhoodContainer",
"location to write strip neighborhoods of EFex RoIs"}
Definition at line 72 of file Egamma1_LArStrip_Fex.h.
◆ m_roiAlgTool
ToolHandle<eFexRoIAlgTool> GlobalSim::Egamma1_LArStrip_Fex::m_roiAlgTool |
|
private |
Initial value:{this,
"roiAlgTool",
"EMB1CellFromCaloCells",
"AlgTool to provide a vector<const xAOD::eFexEMRoI*>"}
Definition at line 54 of file Egamma1_LArStrip_Fex.h.
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
std::optional< std::vector< std::size_t > > wrap3(std::size_t center)
std::string find(const std::string &s)
return a remapped string
uint64_t eventNumber() const
The current event's event number.
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
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)
ToolHandle< ICaloCellsProducer > m_cellProducer
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
std::pair< double, double > Coords
StatusCode dumpTerse(const xAOD::EventInfo &eventInfo, const LArStripNeighborhoodContainer &) const
StatusCode findNeighborhood(const xAOD::eFexEMRoI *, const std::vector< const CaloCell * > &, LArStripNeighborhoodContainer &) const
@ IS_SIMULATION
true: simulation, false: data
uint32_t runNumber() const
The current event's run number.
virtual void setOwner(IDataHandleHolder *o)=0
SG::WriteHandleKey< LArStripNeighborhoodContainer > m_neighKey
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Gaudi::Property< bool > m_dump
AthReentrantAlgorithm()
Default constructor:
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
::StatusCode StatusCode
StatusCode definition for legacy code.
SG::ReadHandleKey< xAOD::EventInfo > m_eventInfoKey
Amg::Vector3D transform(Amg::Vector3D &v, Amg::Transform3D &tr)
Transform a point from a Trasformation3D.
Out copy_if(In first, const In &last, Out res, const Pred &p)
#define CHECK(...)
Evaluate an expression and check for errors.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
std::string to_string(const DetectorType &type)
virtual StatusCode sysInitialize() override
Override sysInitialize.
float eta() const
setter for the above
Gaudi::Property< bool > m_dumpTerse
Data object for each calorimeter readout cell.
#define ATH_MSG_WARNING(x)
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
StatusCode findNeighborhoods(const std::vector< const xAOD::eFexEMRoI * > &, const std::vector< const CaloCell * > &, LArStripNeighborhoodContainer &) const
StatusCode dump(const xAOD::EventInfo &eventInfo, const LArStripNeighborhoodContainer &) const
float mcEventWeight(size_t i=0) const
The weight of one specific MC event used in the simulation.
float phi() const
Seed supercell index within central tower (0 -> 3)
ToolHandle< eFexRoIAlgTool > m_roiAlgTool
float distance(const Amg::Vector3D &p1, const Amg::Vector3D &p2)
calculates the distance between two point in 3D space
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
void dump_n(const LArStripNeighborhood *n, std::ostream &os)
bool eventType(EventType type) const
Check for one particular bitmask value.