TrigMuonEFCbTrack.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
/*****************************************************************************
Name    : TrigMuonEFCbTrack.cxx
Package : offline/Trigger/TrigEvent/TrigMuonEvent
Author  : Sergio Grancagnolo
Created : November 2008
Note    : Simple Trigger AOD for EF Muon (TrigMuonEF)
 
DESCRIPTION:
 
Nov.2008: A TrigMuonEFTrack overloading that contains the match chi2.
*****************************************************************************/
 
// STL include(s):
#include <sstream>
#include <cmath>
 
// Gaudi/Athena include(s):
#include "GaudiKernel/MsgStream.h"
#include "Particle/TrackParticleContainer.h"
#include "Particle/TrackParticle.h"
 
// Local include(s):
#include "TrigMuonEvent/TrigMuonEFCbTrack.h"
 
// "Distance" used by the comparison operator(s):
static const double DELTA = 0.001;
 
TrigMuonEFCbTrack::TrigMuonEFCbTrack() :
    TrigMuonEFTrack( 0., 0., 0., 0. ),
    m_matchChi2(-999.),
    m_nIdSctHits (0),
    m_nIdPixelHits(0),
    m_nTrtHits(0),
    m_idTrackLink()
{
}
 
TrigMuonEFCbTrack::TrigMuonEFCbTrack(double theIPt, double theCotTh, double thePhi, double theM) :
    TrigMuonEFTrack(theIPt, theCotTh, thePhi ,theM ),
    m_matchChi2(-999.),
    m_nIdSctHits (0),
    m_nIdPixelHits(0),
    m_nTrtHits(0),
    m_idTrackLink()
{
}
 
TrigMuonEFCbTrack::TrigMuonEFCbTrack( const TrigMuonEFCbTrack& rhs ) :
  I4Momentum(rhs),
  P4IPtCotThPhiMBase( rhs ),
  INavigable(rhs),
  IAthenaBarCode(),
  INavigable4Momentum(rhs),
  TrigMuonEFTrack(rhs),
  m_matchChi2( rhs.m_matchChi2 ),
  m_nIdSctHits (rhs.m_nIdSctHits),
  m_nIdPixelHits(rhs.m_nIdPixelHits),
  m_nTrtHits(rhs.m_nTrtHits),
  m_idTrackLink(rhs.m_idTrackLink)
{ }
 
TrigMuonEFCbTrack::TrigMuonEFCbTrack( const TrigMuonEFTrack& rhs ) :
  I4Momentum(rhs),
  P4IPtCotThPhiMBase( rhs ),
  INavigable(rhs),
  IAthenaBarCode(),
  INavigable4Momentum(rhs),
  TrigMuonEFTrack(rhs),
  m_matchChi2(-999.),
  m_nIdSctHits (0),
  m_nIdPixelHits(0),
  m_nTrtHits(0),
  m_idTrackLink()
  { }

void TrigMuonEFCbTrack::setMatchChi2(double theMatchChi2){
  m_matchChi2 = theMatchChi2;
}

TrigMuonEFCbTrack::~TrigMuonEFCbTrack() {}

const Rec::TrackParticle* TrigMuonEFCbTrack::getIDTrackParticle() const {
  if(!m_idTrackLink.isValid()) return 0; //possible since old data has no link
  return *(m_idTrackLink);
}

void TrigMuonEFCbTrack::setIDTrack(const Rec::TrackParticleContainer* cont, unsigned int itrk) {
  m_idTrackLink.toIndexedElement(*cont, itrk);
}

void TrigMuonEFCbTrack::setIDTrackElementLink(const ElementLink<Rec::TrackParticleContainer>& link) {
  m_idTrackLink = link;
}

// helper operators

std::string str( const TrigMuonEFCbTrack& d ) {
   std::stringstream ss;
   ss << str( static_cast<const TrigMuonEFTrack&>(d))
      << "; chi2 match: " << d.matchChi2();
 
   return ss.str();
}
 
MsgStream& operator<< ( MsgStream& m, const TrigMuonEFCbTrack& d ) {
 
   return ( m << str( d ) );
}

bool operator== ( const TrigMuonEFCbTrack& a, const TrigMuonEFCbTrack& b ) {
 
   if( static_cast<const TrigMuonEFTrack&>(a) != static_cast<const TrigMuonEFTrack&>(b) )
      return false;
   if( std::abs( a.matchChi2() - b.matchChi2() ) > DELTA )
      return false;
 
   return true;
}

void diff( const TrigMuonEFCbTrack& a, const TrigMuonEFCbTrack& b,
           std::map< std::string, double >& variableChange ) {
 
   // Call the function operating on the base class:
   diff( static_cast<const TrigMuonEFTrack&>(a), static_cast<const TrigMuonEFTrack&>(b), variableChange );
 
   if( std::abs( a.matchChi2() - b.matchChi2() ) > DELTA ) {
      variableChange[ "MatchChi2" ] = a.matchChi2() - b.matchChi2();
   }
 
   return;
}