|
ATLAS Offline Software
|
Go to the documentation of this file.
26 const std::set<CaloSampling::CaloSample> samplings{
35 unsigned sampling_pattern = 0;
36 for (
auto sample : samplings) {
37 sampling_pattern |= 0x1U <<
sample;
41 for (
auto sample : samplings) {
46 if (std::abs(
eta) < 1.45) {
62 etaCaloDecor(*cluster) =
eta;
63 phiCaloDecor(*cluster) =
phi;
77 using namespace asg::CreateDummyEl;
79 setMsgLevel(MSG::INFO);
83 std::unique_ptr<xAOD::EventInfo> eventInfo =
84 std::make_unique<xAOD::EventInfo>();
91 if (!
store.
record(std::move(eventInfo),
"EventInfo").isSuccess()) {
93 return StatusCode::FAILURE;
96 size_t numel = pt_eta.size();
98 std::unique_ptr<xAOD::CaloClusterContainer>
clusters =
99 std::make_unique<xAOD::CaloClusterContainer>();
100 std::unique_ptr<xAOD::CaloClusterAuxContainer> clAux =
101 std::make_unique<xAOD::CaloClusterAuxContainer>();
104 for (
const auto&
i : pt_eta) {
105 const double pt{
i.first };
106 const double eta{
i.second };
107 const double e{
pt * cosh(eta) };
110 fill_cluster(cluster, eta, 0.0,
e);
114 !
store.
record(std::move(clAux),
"MyClustersAux.").isSuccess()) {
116 return StatusCode::FAILURE;
119 std::vector<ElementLink<xAOD::CaloClusterContainer>>
links{};
120 std::unique_ptr<xAOD::ElectronContainer>
electrons =
121 std::make_unique<xAOD::ElectronContainer>();
122 std::unique_ptr<xAOD::ElectronAuxContainer> electronsAux =
123 std::make_unique<xAOD::ElectronAuxContainer>();
126 for (
const auto&
i : pt_eta) {
127 const double pt{
i.first };
128 const double eta{
i.second };
132 links.emplace_back(
"MyClusters",
el->index());
133 el->setCaloClusterLinks(
links);
140 !
store.
record(std::move(electronsAux),
"MyElectronAux.").isSuccess()) {
142 return StatusCode::FAILURE;
145 return StatusCode::SUCCESS;
Electron_v1 Electron
Definition of the current "egamma version".
void setEventNumber(uint64_t value)
Set the current event's event number.
Scalar phi() const
phi method
Scalar eta() const
pseudorapidity method
void insertMoment(MomentType type, double value)
void setSamplingPattern(const unsigned sp, const bool clearSamplingVars=false)
Set sampling pattern (one bit per sampling.
bool setEnergy(const CaloSample sampling, const float e)
Set energy for a given sampling. Returns false if the sample isn't part of the cluster.
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
Description of a calorimeter cluster.
::StatusCode StatusCode
StatusCode definition for legacy code.
ANA_MSG_HEADER(msgSTT) ANA_MSG_SOURCE(msgSTT
StatusCode getElectrons(const std::vector< std::pair< double, double >> &pt_eta, int runNumber, xAOD::TStore &store)
void makePrivateStore()
Create a new (empty) private store for this object.
void record(const T *p, const std::string &key)
@ PHICALOFRAME
Phi in the calo frame (for egamma)
A relatively simple transient store for objects created in analysis.
void setRunNumber(uint32_t value)
Set the current event's run number.
bool setPhi(const CaloSample sampling, const float phi)
Set in a given sampling. Returns false if the sample isn't part of the cluster.
bool setEta(const CaloSample sampling, const float eta)
Set in a given sampling. Returns false if the sample isn't part of the cluster.
@ ETACALOFRAME
Eta in the calo frame (for egamma)