Electron_v1.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// $Id: Electron_v1$
 
// EDM include(s):
#include "xAODCore/AuxStoreAccessorMacros.h"
#include "TruthUtils/ParticleConstants.h"
// Local include(s):
#include "xAODEgamma/versions/Electron_v1.h"
#include "ElectronAccessors_v1.h"
 
#include <stdexcept>
 
namespace xAOD {
 
  Electron_v1::Electron_v1() : Egamma_v1() {
  }
 
  Electron_v1::Electron_v1(const Electron_v1& el) : Egamma_v1(el){
    this->makePrivateStore(el);  
  }
 
  Electron_v1& Electron_v1::operator=(const Electron_v1& el ){
 
    if (this != &el){ // protect against invalid self-assignment      
      this->Egamma_v1::operator=( el );
    }
    // by convention, always return *this
    return *this;
  }
 
  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( Electron_v1, float, charge, setCharge )
 
  Type::ObjectType Electron_v1::type() const {
    return Type::Electron;
  }
 
  double Electron_v1::m() const {
    return ParticleConstants::electronMassInMeV;
  }

  //
  //         Implementation of the TrackParticle accessor functions
  //
 
  size_t xAOD::Electron_v1::nTrackParticles() const {
 
     static const SG::AuxElement::Accessor< Electron_v1::TPELVec_t >
        trackAcc( "trackParticleLinks" );
 
     if( trackAcc.isAvailable( *this ) ) {
        return trackAcc( *this ).size();
     }
     return 0;
  }
 
  const xAOD::TrackParticle* Electron_v1::trackParticle( size_t index ) const {
 
     if( index >= nTrackParticles() ) {
        return nullptr;
     }
     const TPELVec_t& links = trackParticleLinks();
     if( ! links[ index ].isValid() ) {
        return nullptr;
     }
     return *( links[ index ] );
  }
 
  const ElementLink< TrackParticleContainer >&
  Electron_v1::trackParticleLink( size_t index ) const {
 
     if( index >= nTrackParticles() ) {
        static const ElementLink< TrackParticleContainer > dummy;
        return dummy;
     }
     return trackParticleLinks()[ index ];
  }
 
  AUXSTORE_OBJECT_SETTER_AND_GETTER( Electron_v1, Electron_v1::TPELVec_t,
                                     trackParticleLinks, setTrackParticleLinks )
 
  //

  //Implementation of the summary value accessor functions
  
  bool Electron_v1::trackCaloMatchValue( float& value, const EgammaParameters::TrackCaloMatchType information ) const {
 
    const xAOD::Electron_v1::Accessor< float >* acc = trackCaloMatchAccessorV1( information );
    if( !acc ) { 
      return false;
    }
    if(!acc->isAvailable( *this) ){
      return  false;
    }
    // Retrieve the value:
    value = ( *acc )( *this );
    return true;
  }
 
  float Electron_v1::trackCaloMatchValue( const EgammaParameters::TrackCaloMatchType information ) const {
    const xAOD::Electron_v1::Accessor< float >* acc = trackCaloMatchAccessorV1( information );
    if(!acc ) throw std::runtime_error( "Unknown/Unavailable Track to Calo Match type requested" );
    return ( *acc )( *this );
  }
 
  bool Electron_v1::setTrackCaloMatchValue( const float value, 
                                            const EgammaParameters::TrackCaloMatchType information ) {
 
    const xAOD::Electron_v1::Accessor< float >* acc = trackCaloMatchAccessorV1( information );
    if( !acc ) return false;
    // Set the value:
    ( *acc )( *this ) = value;
    return true;
  }
 
  bool Electron_v1::setTrackCaloMatchValues(
    const std::array<double, 4> &deltaEta,
    const std::array<double, 4> &deltaPhi,
    const std::array<double, 4> &deltaPhiRescaled,
    const double deltaPhiLast
  ) {
    bool ret = true;
    // It's tempting, but clang warns about chaining these together with &.
    if (!setTrackCaloMatchValue(static_cast<float> (deltaEta[0]), xAOD::EgammaParameters::deltaEta0)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaEta[1]), xAOD::EgammaParameters::deltaEta1)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaEta[2]), xAOD::EgammaParameters::deltaEta2)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaEta[3]), xAOD::EgammaParameters::deltaEta3)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaPhi[0]), xAOD::EgammaParameters::deltaPhi0)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaPhi[1]), xAOD::EgammaParameters::deltaPhi1)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaPhi[2]), xAOD::EgammaParameters::deltaPhi2)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaPhi[3]), xAOD::EgammaParameters::deltaPhi3)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaPhiRescaled[0]), xAOD::EgammaParameters::deltaPhiRescaled0)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaPhiRescaled[1]), xAOD::EgammaParameters::deltaPhiRescaled1)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaPhiRescaled[2]), xAOD::EgammaParameters::deltaPhiRescaled2)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaPhiRescaled[3]), xAOD::EgammaParameters::deltaPhiRescaled3)) ret = false;
    if (!setTrackCaloMatchValue(static_cast<float> (deltaPhiLast), xAOD::EgammaParameters::deltaPhiFromLastMeasurement)) ret = false;
    return ret;
  }
 
  bool Electron_v1::trackParticleSummaryValue( uint8_t& value, const SummaryType information, int index ) const {
 
    const xAOD::TrackParticle* tempTrackParticle = trackParticle(index);
    if (!tempTrackParticle) return false;
    return tempTrackParticle->summaryValue(value,information);
  }
 
  bool Electron_v1::trackParticleSummaryValue( float& value, const SummaryType information, int index ) const {
    
    const xAOD::TrackParticle* tempTrackParticle = trackParticle(index);
    if (!tempTrackParticle) return false;
    return tempTrackParticle->summaryValue(value,information);
  }
  
  uint8_t Electron_v1::trackParticleSummaryIntValue( const SummaryType information, int index) const {
    const xAOD::TrackParticle* tempTrackParticle = trackParticle(index);
    if (!tempTrackParticle) {
      throw std::runtime_error( "TrackParticle not available" );
    }
    uint8_t value=0;
    if(!tempTrackParticle->summaryValue(value,information)){
      throw std::runtime_error( "Unknown/Unavailable Int Track Summary type requested" );
    }
    return value;
  }
 
  float Electron_v1::trackParticleSummaryFloatValue( const SummaryType information, int index) const {
    const xAOD::TrackParticle* tempTrackParticle = trackParticle(index);
    if (!tempTrackParticle) {
      throw std::runtime_error( "TrackParticle not available" );
    }
    float value=0;
    if(!tempTrackParticle->summaryValue(value,information)){
      throw std::runtime_error( "Unknown/Unavailable Float Track Summary type requested" );
    }
    return value;
  }

 
 
} // namespace xAOD