MuonSimHit_V1.cxx

Go to the documentation of this file.

/*
   Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
// EDM include(s):
#include <xAODCore/AuxStoreAccessorMacros.h>
 
#include <xAODMuonSimHit/versions/MuonSimHit_v1.h>
 
namespace {
   static const std::string preFixStr{"MuSim_"};
   static const SG::AuxElement::Accessor<Identifier::value_type> acc_Identifier{preFixStr+"identifier"};
   static const xAOD::PosAccessor<3> acc_localPos{"localPositionDim3"};
   static const xAOD::PosAccessor<3> acc_localDir{"localDirectionDim3"};
 
   static const SG::AuxElement::Accessor<unsigned short> acc_mcEventIndex{preFixStr+"mcEventIndex"};
   static const SG::AuxElement::Accessor<unsigned int> acc_uniqueID{preFixStr+"uniqueID"};
}
 
#define IMPLEMENT_SETTER_GETTER( DTYPE, GETTER, SETTER)                          \
      DTYPE MuonSimHit_v1::GETTER() const {                                      \
         static const SG::AuxElement::Accessor<DTYPE> acc{preFixStr + #GETTER};  \
         return acc(*this);                                                      \
      }                                                                          \
                                                                                 \
      void MuonSimHit_v1::SETTER(DTYPE value) {                        \
         static const SG::AuxElement::Accessor<DTYPE> acc{preFixStr + #GETTER};  \
         acc(*this) = value;                                                     \
      }
          
namespace xAOD {
 
MuonSimHit_v1& MuonSimHit_v1::operator=(const MuonSimHit_v1& other) {
   if (this != &other) {
      static_cast<SG::AuxElement&>(*this) = other;
#ifndef __CLING__
      releaseParticleLink();
#endif
   }
   return (*this);
}
Identifier MuonSimHit_v1::identify() const { return Identifier{acc_Identifier(*this)}; }    
void MuonSimHit_v1::setIdentifier(const Identifier& id) { acc_Identifier(*this) = id.get_compact(); }
 
IMPLEMENT_SETTER_GETTER(float, globalTime, setGlobalTime)
IMPLEMENT_SETTER_GETTER(int, pdgId, setPdgId)
IMPLEMENT_SETTER_GETTER(float, energyDeposit, setEnergyDeposit)
IMPLEMENT_SETTER_GETTER(float, kineticEnergy, setKineticEnergy)
IMPLEMENT_SETTER_GETTER(float, mass, setMass)
IMPLEMENT_SETTER_GETTER(float, stepLength, setStepLength)
 
float MuonSimHit_v1::beta() const{
   const float e = kineticEnergy();
   return std::sqrt(1. - std::min(std::pow(mass() / (e > std::numeric_limits<float>::epsilon() ? e : 1.f), 2),1.));
}
void MuonSimHit_v1::setLocalPosition(MeasVector<3> vec) {
   VectorMap<3> lPos{acc_localPos(*this).data()};
   lPos = std::move(vec);   
}
ConstVectorMap<3> MuonSimHit_v1::localPosition() const { return ConstVectorMap<3>{acc_localPos(*this).data()};}
 
void MuonSimHit_v1::setLocalDirection(MeasVector<3> vec) {
   VectorMap<3> lPos{acc_localDir(*this).data()};
   lPos = std::move(vec);   
}
ConstVectorMap<3> MuonSimHit_v1::localDirection() const { return ConstVectorMap<3>{acc_localDir(*this).data()};}
 
const HepMcParticleLink& MuonSimHit_v1::genParticleLink() const {
   if (!m_hepMCLink) {
      auto link = std::make_unique<HepMcParticleLink>(acc_uniqueID(*this),
                                                      acc_mcEventIndex(*this),
                                                      HepMcParticleLink::IS_POSITION, 
                                                      HepMcParticleLink::IS_ID);
      return *m_hepMCLink.set(std::move(link));
   }
   return (*m_hepMCLink);
}
void MuonSimHit_v1::setGenParticleLink(const HepMcParticleLink& link) {
   releaseParticleLink();
   acc_mcEventIndex(*this) = link.getEventPositionInCollection(SG::CurrentEventStore::store());
   acc_uniqueID(*this) = link.id();
}
 
std::unique_ptr<const HepMcParticleLink> MuonSimHit_v1::releaseParticleLink() {
   return m_hepMCLink.release();
}
 
}
#undef IMPLEMENT_SETTER_GETTER