|
ATLAS Offline Software
|
#include <AODReader.h>
|
| AODReader (const std::string &name, ISvcLocator *pSvcLocator) |
|
virtual | ~AODReader () |
|
virtual StatusCode | initialize () |
|
virtual StatusCode | execute () |
|
virtual StatusCode | finalize () |
|
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 > &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 |
|
Definition at line 14 of file AODReader.h.
◆ StoreGateSvc_t
◆ AODReader()
AODReader::AODReader |
( |
const std::string & |
name, |
|
|
ISvcLocator * |
pSvcLocator |
|
) |
| |
◆ ~AODReader()
AODReader::~AODReader |
( |
| ) |
|
|
virtual |
◆ 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);
253 return PBASE::declareProperty(
name,hndl,
doc);
◆ 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.
229 return PBASE::declareProperty(
name,hndl,
doc);
◆ 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.
338 return PBASE::declareProperty(
name, property,
doc);
◆ 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()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode AODReader::execute |
( |
| ) |
|
|
virtual |
read the AOD electron container from persistecy storage
loop over the AOD electron container, take the first one (only one?) with the correct author
Definition at line 105 of file AODReader.cxx.
111 sc=
evtStore()->retrieve( elecTES,
"ElectronAODCollection");
112 if(
sc.isFailure() || !elecTES ) {
113 msg(MSG::WARNING) <<
"No AOD electron container found in TDS" <<
endmsg;
116 msg(MSG::INFO) <<
"ElectronContainer successfully retrieved" <<
endmsg;
120 if (( mcEvtColl =
evtStore()->retrieve<const McEventCollection>() )) {
121 msg(MSG::INFO) <<
"TruthParticleContainer successfully retrieved" <<
endmsg;
123 msg(MSG::WARNING) <<
"Could not retrieve TruthParticleContainer" <<
endmsg;
132 if (mcTrPart != mcEvtColl->
end()) {
133 trPart = (*mcTrPart)->particles_size()?*
HepMC::begin(**mcTrPart):nullptr;
135 msg(MSG::WARNING) <<
"Not a single particle event. Truth information won't be available" <<
endmsg;
138 msg(MSG::WARNING) <<
"Empty TruthParticleContainer. Truth information won't be available" <<
endmsg;
147 for (; elecItr != elecItrE; ++elecItr) {
148 if (((*elecItr)->author() & (8 + 1)) > 0) {
149 primElec = (*elecItr);
203 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 & 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.
57 return Algorithm::extraOutputDeps();
◆ finalize()
StatusCode AODReader::finalize |
( |
| ) |
|
|
virtual |
Definition at line 99 of file AODReader.cxx.
102 return StatusCode::SUCCESS;
◆ initialize()
StatusCode AODReader::initialize |
( |
| ) |
|
|
virtual |
Definition at line 28 of file AODReader.cxx.
33 msg(MSG::INFO) <<
"creating n tuple" <<
endmsg;
34 NTupleFilePtr file1(
ntupleSvc(),
"/NTUPLES/FILE");
40 nt =
ntupleSvc()->book (
"/NTUPLES/FILE/AOD", CLID_ColumnWiseTuple,
"Readed info from aod container");
43 msg(MSG::INFO) <<
"booked ntuple " <<
endmsg;
48 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for eta" <<
endmsg;
50 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for erec" <<
endmsg;
52 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for pt" <<
endmsg;
54 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for e237" <<
endmsg;
56 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for e277" <<
endmsg;
58 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for weta1" <<
endmsg;
60 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for weta2" <<
endmsg;
62 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for e2tsts1" <<
endmsg;
64 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for emins1" <<
endmsg;
66 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for fracs1" <<
endmsg;
68 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for wtots1" <<
endmsg;
70 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for f1" <<
endmsg;
72 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for f1core" <<
endmsg;
74 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for ethad1" <<
endmsg;
76 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for et" <<
endmsg;
79 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for egen" <<
endmsg;
81 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for tr_px" <<
endmsg;
83 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for tr_py" <<
endmsg;
85 if( sc.isFailure() )
msg(MSG::ERROR) <<
"Can not create ntuple entry for tr_pz" <<
endmsg;
89 msg(MSG::ERROR) <<
"didn't manage to book ntuple" <<
endmsg;
90 return StatusCode::FAILURE;
96 return StatusCode::SUCCESS;
◆ 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.
◆ 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()
◆ sysInitialize()
StatusCode AthAlgorithm::sysInitialize |
( |
| ) |
|
|
overridevirtualinherited |
◆ 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_detStore
◆ m_e237
NTuple::Item<float> AODReader::m_e237 |
|
private |
◆ m_e277
NTuple::Item<float> AODReader::m_e277 |
|
private |
◆ m_e2tsts1
NTuple::Item<float> AODReader::m_e2tsts1 |
|
private |
◆ m_emins1
NTuple::Item<float> AODReader::m_emins1 |
|
private |
◆ m_energy
NTuple::Item<float> AODReader::m_energy |
|
private |
◆ m_et
NTuple::Item<float> AODReader::m_et |
|
private |
◆ m_eta
NTuple::Item<float> AODReader::m_eta |
|
private |
◆ m_ethad1
NTuple::Item<float> AODReader::m_ethad1 |
|
private |
◆ m_evtStore
◆ m_extendedExtraObjects
DataObjIDColl AthAlgorithm::m_extendedExtraObjects |
|
privateinherited |
◆ m_f1
NTuple::Item<float> AODReader::m_f1 |
|
private |
◆ m_f1core
NTuple::Item<float> AODReader::m_f1core |
|
private |
◆ m_fracs1
NTuple::Item<float> AODReader::m_fracs1 |
|
private |
◆ m_nt
NTuple::Tuple* AODReader::m_nt |
|
private |
◆ m_pt
NTuple::Item<float> AODReader::m_pt |
|
private |
◆ m_truth_energy
NTuple::Item<float> AODReader::m_truth_energy |
|
private |
◆ m_truth_px
NTuple::Item<float> AODReader::m_truth_px |
|
private |
◆ m_truth_py
NTuple::Item<float> AODReader::m_truth_py |
|
private |
◆ m_truth_pz
NTuple::Item<float> AODReader::m_truth_pz |
|
private |
◆ m_varHandleArraysDeclared
◆ m_vhka
◆ m_weta1
NTuple::Item<float> AODReader::m_weta1 |
|
private |
◆ m_weta2
NTuple::Item<float> AODReader::m_weta2 |
|
private |
◆ m_wtots1
NTuple::Item<float> AODReader::m_wtots1 |
|
private |
The documentation for this class was generated from the following files:
Const iterator class for DataVector/DataList.
virtual double pt() const
transverse momentum
double e237() const
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x7
NTuple::Item< float > m_e2tsts1
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
NTuple::Item< float > m_emins1
NTuple::Item< float > m_f1
double weta1() const
shower width using +/-3 strips around the one with the maximal energy deposit: w3 strips = sqrt{sum(E...
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
NTuple::Item< float > m_weta2
std::vector< SG::VarHandleKeyArray * > m_vhka
NTuple::Item< float > m_e277
NTuple::Item< float > m_truth_py
const T * detail(const std::string &name="", unsigned int index=0) const
retrieve eg-detail objects:
virtual void setOwner(IDataHandleHolder *o)=0
virtual double e() const
energy
NTuple::Item< float > m_e237
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysInitialize() override
Override sysInitialize.
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
NTuple::Item< float > m_ethad1
::StatusCode StatusCode
StatusCode definition for legacy code.
double e277() const
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 7x7
double fracs1() const
shower shape in the shower core : [E(+/-3)-E(+/-1)]/E(+/-1), where E(+/-n) is the energy in +- n stri...
const CaloCluster * cluster() const
pointer to CaloCluster
NTuple::Item< float > m_truth_pz
NTuple::Item< float > m_weta1
NTuple::Item< float > m_energy
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
NTuple::Item< float > m_eta
virtual double eta() const
pseudo rapidity
This defines the McEventCollection, which is really just an ObjectVector of McEvent objects.
double weta2() const
the lateral width is calculated with a window of 3x5 cells using the energy weighted sum over all cel...
NTuple::Item< float > m_et
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
double e2tsts1() const
energy of the cell corresponding to second energy maximum in the first sampling
const GenParticle * ConstGenParticlePtr
NTuple::Item< float > m_pt
NTuple::Item< float > m_fracs1
NTuple::Item< float > m_f1core
double wtots1() const
shower width is determined in a window detaxdphi = 0,0625 x~0,2, corresponding typically to 20 strips...
double emins1() const
energy reconstructed in the strip with the minimal value between the first and second maximum
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
DataObjIDColl m_extendedExtraObjects
double f1core() const
E1(3x1)/E = fraction of the energy reconstructed in the first longitudinal compartment of the electro...
#define ATH_MSG_WARNING(x)
NTuple::Item< float > m_wtots1
NTuple::Item< float > m_truth_px
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
GenEvent::particle_iterator begin(HepMC::GenEvent &e)
AthAlgorithm()
Default constructor:
virtual double et() const
transverse energy defined to be e*sin(theta)
NTuple::Item< float > m_truth_energy
double ethad1() const
transverse energy in the first sampling of the hadronic calorimeters behind the cluster calculated fr...
double f1() const
E1/E = fraction of energy reconstructed in the first sampling, where E1 is energy in all strips belon...
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.