#include <FPGATrackSimTrackMonitor.h>

Inheritance diagram for FPGATrackSimTrackMonitor:

Collaboration diagram for FPGATrackSimTrackMonitor:

Public Member Functions
	FPGATrackSimTrackMonitor (const std::string &, const std::string &, const IInterface *)
virtual StatusCode	initialize () override
void	fillRoad (const std::vector< FPGATrackSimRoad > &roads, const std::vector< FPGATrackSimTruthTrack > &truthTracks, size_t nLogicalLayers)
void	fillTrack (const std::vector< const FPGATrackSimTrack * > &tracks, const std::vector< FPGATrackSimTruthTrack > &truthTracks, float chi2Cut)
unsigned	getNElements () const
unsigned	getMaxNElements () const
unsigned	getTotNElements () const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ToolHandle< GenericMonitoringTool >	m_monTool {this,"MonTool", "", "Monitoring tool"}
Gaudi::Property< std::string >	m_dir {this, "dir", {"/TRACKMON/"}, "String name of output directory"}
ServiceHandle< ITHistSvc >	m_tHistSvc {this, "THistSvc", "THistSvc/THistSvc", "Histogramming service"}
unsigned	m_nElements {0}
unsigned	m_maxNElements {0}
unsigned	m_totNElements {0}
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 37 of file FPGATrackSimTrackMonitor.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ FPGATrackSimTrackMonitor()

FPGATrackSimTrackMonitor::FPGATrackSimTrackMonitor	(	const std::string &	algname,
		const std::string &	name,
		const IInterface *	ifc )

Definition at line 19 of file FPGATrackSimTrackMonitor.cxx.

                                                                                                                         :
  AthAlgTool(algname, name, ifc)
{
  ATH_MSG_INFO("constructor");
}

Member Function Documentation

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	hndl,
		const SG::VarHandleKeyType &	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95{ return m_detStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85{ return m_evtStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

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

◆ fillRoad()

void FPGATrackSimTrackMonitor::fillRoad	(	const std::vector< FPGATrackSimRoad > &	roads,
		const std::vector< FPGATrackSimTruthTrack > &	truthTracks,
		size_t	nLogicalLayers )

Definition at line 51 of file FPGATrackSimTrackMonitor.cxx.

{
  // overall number of roads
  auto nRoads = Monitored::Scalar<int>("nRoads", roads.size());
  // push them into the monitoring tool (this->MonTool set in python)
  Monitored::Group(m_monTool, nRoads);
 
  // update counter of total roads
  m_totNElements += roads.size();
 
  // loop over all roads in the event
  for (const auto& road : roads) {
    // layerIDs
    unsigned bitmask = road.getHitLayers();
    for (size_t l = 0; l < nLogicalLayers; l++) {
      if (bitmask & (1 << l)) {
          auto mon_layerIDs = Monitored::Scalar<unsigned>("layerIDs",l);
          Monitored::Group(m_monTool, mon_layerIDs);
      }
    }
    // number of layers represented in this road
    auto mon_nLayers = Monitored::Scalar<int>("nLayers", road.getNLayers());
    // number of hits
    auto mon_nHits = Monitored::Scalar<int>("nHits", road.getNHits());
    // push all of the above scalars into the monitoring tool
    Monitored::Group(m_monTool, mon_nLayers, mon_nHits);
    // roads with 4 hits
    if (road.getNHits() == 4) {
      // fit chi2 in this road
      auto mon_fitChi2_4 = Monitored::Scalar<float>("chi2_4", road.getFitChi2());
      // fit chi2Eta
      auto mon_fitChi2Eta_4 = Monitored::Scalar<float>("chi2Eta_4", road.getFitChi2Eta());
      // fit chi2Phi
      auto mon_fitChi2Phi_4 = Monitored::Scalar<float>("chi2Phi_4", road.getFitChi2Phi());
      // push all of the above scalars into the monitoring tool
      Monitored::Group(m_monTool, mon_fitChi2_4, mon_fitChi2Eta_4, mon_fitChi2Phi_4);
    } else if (road.getNHits() == 5) {
      auto mon_fitChi2_5 = Monitored::Scalar<float>("chi2_5", road.getFitChi2());
      auto mon_fitChi2Eta_5 = Monitored::Scalar<float>("chi2Eta_5", road.getFitChi2Eta());
      auto mon_fitChi2Phi_5 = Monitored::Scalar<float>("chi2Phi_5", road.getFitChi2Phi());
      Monitored::Group(m_monTool, mon_fitChi2_5, mon_fitChi2Eta_5, mon_fitChi2Phi_5);
    }
  }
 
  // only monitor efficiency for events with truth tracks
  if (truthTracks.size() > 0) {
    // update road counter
    if (roads.size() > 0) ++m_nElements;
    // update the max number of elements
    if (roads.size() > m_maxNElements) m_maxNElements=roads.size();
    // use the first truth track (the primary one)
    const auto& truthTrack = truthTracks.front();
    // efficiency monitor
    auto mon_eff = Monitored::Scalar<bool>("eff_road", (roads.size() > 0));
    auto mon_truth_pt_zoom = Monitored::Scalar<float>("pT_zoom", truthTrack.getPt()*0.001);
    auto mon_truth_pt = Monitored::Scalar<float>("pT", truthTrack.getPt()*0.001);
    auto mon_truth_eta = Monitored::Scalar<float>("eta", truthTrack.getEta());
    auto mon_truth_phi = Monitored::Scalar<float>("phi", truthTrack.getPhi());
    auto mon_truth_d0 = Monitored::Scalar<float>("d0", truthTrack.getD0());
    auto mon_truth_z0 = Monitored::Scalar<float>("z0", truthTrack.getZ0());
    Monitored::Group( m_monTool, mon_eff, mon_truth_pt_zoom, mon_truth_pt, mon_truth_eta, mon_truth_phi, mon_truth_d0, mon_truth_z0 );
  }
}

◆ fillTrack()

void FPGATrackSimTrackMonitor::fillTrack	(	const std::vector< const FPGATrackSimTrack * > &	tracks,
		const std::vector< FPGATrackSimTruthTrack > &	truthTracks,
		float	chi2Cut )

Definition at line 119 of file FPGATrackSimTrackMonitor.cxx.

{
  // monitor chi2 for all tracks and find the best
  float bestChi2 = chi2Cut;
  float chi2_best = chi2Cut;
  float chi2Eta_best = chi2Cut;
  float chi2Phi_best = chi2Cut;
  float chi2_best_4 = chi2Cut;
  float chi2Eta_best_4 = chi2Cut;
  float chi2Phi_best_4 = chi2Cut;
  float chi2_best_5 = chi2Cut;
  float chi2Eta_best_5 = chi2Cut;
  float chi2Phi_best_5 = chi2Cut;
 
  // lambda function to monitor chi2 (all) and keep track of the best one
  auto monitorChi2 = [&](const std::string& chi2_name, float chi2_value, float& best_value, const std::string& suffix="") {
    if (chi2_value < chi2Cut) {
      // monitor all chi2
      auto mon = Monitored::Scalar<float>(chi2_name + suffix, chi2_value);
      Monitored::Group(m_monTool, mon);
      // update best chi2
      if (chi2_value < best_value) best_value = chi2_value;
    }
  };
 
 
  // overall number of tracks
  auto mon_nTracks = Monitored::Scalar<int>("nTracks", tracks.size());
  // push them into the monitoring tool (this->MonTool set in python)
  Monitored::Group(m_monTool, mon_nTracks);
 
  // update counter of total tracks
  m_totNElements += tracks.size();
 
  // loop over all tracks in the event
  for (const auto& track : tracks) {
    // number of real hits
    int nValidHits = 0;
    for (const auto& hit : track->getFPGATrackSimHits()) {
      if (hit.isReal()) {
        ++nValidHits;
      }
    }
    auto mon_nHits = Monitored::Scalar<int>("nHits", nValidHits);
    Monitored::Group(m_monTool, mon_nHits);
 
    // monitor chi2 for tracks
    float chi2_all = track->getChi2ndof();
    if (chi2_all < chi2Cut) {
      auto mon_chi2_all = Monitored::Scalar<float>("chi2_all", chi2_all);
      Monitored::Group(m_monTool, mon_chi2_all);
      // update best chi2
      if (chi2_all < bestChi2) {
        bestChi2 = chi2_all;
      }
    }
    
    // monitor base chi2 values
    monitorChi2("chi2",    track->getChi2(),    chi2_best);
    monitorChi2("chi2Eta", track->getChi2Eta(), chi2Eta_best);
    monitorChi2("chi2Phi", track->getChi2Phi(), chi2Phi_best);
    
    // monitor nHits-dependent chi2 values
    if (nValidHits == 4) {
      monitorChi2("chi2",    track->getChi2(),    chi2_best_4,    "_4");
      monitorChi2("chi2Eta", track->getChi2Eta(), chi2Eta_best_4, "_4");
      monitorChi2("chi2Phi", track->getChi2Phi(), chi2Phi_best_4, "_4");
    } else if (nValidHits == 5) {
      monitorChi2("chi2",    track->getChi2(),    chi2_best_5,    "_5");
      monitorChi2("chi2Eta", track->getChi2Eta(), chi2Eta_best_5, "_5");
      monitorChi2("chi2Phi", track->getChi2Phi(), chi2Phi_best_5, "_5");
    }
  }
  // monitor the best chi2 found
  if (bestChi2 < chi2Cut) {
    auto mon_best_chi2 = Monitored::Scalar<float>("best_chi2", bestChi2);
    Monitored::Group(m_monTool, mon_best_chi2);
  }
 
  // lambda function to monitor the best chi2 value
  auto monitorBestChi2 = [&](const std::string& best_name, float best_value) {
    if (best_value < chi2Cut) {
      auto mon_best = Monitored::Scalar<float>(best_name, best_value);
      Monitored::Group(m_monTool, mon_best);
    }
  };
 
  // best chi2 values
  monitorBestChi2("chi2_best",    chi2_best);
  monitorBestChi2("chi2Eta_best", chi2Eta_best);
  monitorBestChi2("chi2Phi_best", chi2Phi_best);
  monitorBestChi2("chi2_best_4",    chi2_best_4);
  monitorBestChi2("chi2Eta_best_4", chi2Eta_best_4);
  monitorBestChi2("chi2Phi_best_4", chi2Phi_best_4);
  monitorBestChi2("chi2_best_5",    chi2_best_5);
  monitorBestChi2("chi2Eta_best_5", chi2Eta_best_5);
  monitorBestChi2("chi2Phi_best_5", chi2Phi_best_5);
 
 
  // only monitor efficiency for events with truth tracks
  if (truthTracks.size() > 0) {
    // update track counter
    if (tracks.size() > 0) ++m_nElements;
    // update the max number of elements
    if (tracks.size() > m_maxNElements) m_maxNElements=tracks.size();
    // use the first truth track (the primary one)
    const auto& truthTrack = truthTracks.front();
    // efficiency monitor
    auto mon_eff = Monitored::Scalar<bool>("eff_track", (tracks.size() > 0));
    auto mon_truth_pt_zoom = Monitored::Scalar<float>("pT_zoom", truthTrack.getPt()*0.001);
    auto mon_truth_pt = Monitored::Scalar<float>("pT", truthTrack.getPt()*0.001f);
    auto mon_truth_eta = Monitored::Scalar<float>("eta", truthTrack.getEta());
    auto mon_truth_phi = Monitored::Scalar<float>("phi", truthTrack.getPhi());
    auto mon_truth_d0 = Monitored::Scalar<float>("d0", truthTrack.getD0());
    auto mon_truth_z0 = Monitored::Scalar<float>("z0", truthTrack.getZ0());
    Monitored::Group( m_monTool, mon_eff, mon_truth_pt_zoom, mon_truth_pt, mon_truth_eta, mon_truth_phi, mon_truth_d0, mon_truth_z0 );
  }
}

◆ getMaxNElements()

unsigned FPGATrackSimTrackMonitor::getMaxNElements ( ) const

inline

Definition at line 67 of file FPGATrackSimTrackMonitor.h.

67{return m_maxNElements;}

◆ getNElements()

unsigned FPGATrackSimTrackMonitor::getNElements ( ) const

inline

Definition at line 64 of file FPGATrackSimTrackMonitor.h.

64{return m_nElements;}

◆ getTotNElements()

unsigned FPGATrackSimTrackMonitor::getTotNElements ( ) const

inline

Definition at line 70 of file FPGATrackSimTrackMonitor.h.

70{return m_totNElements;}

◆ initialize()

StatusCode FPGATrackSimTrackMonitor::initialize ( )

overridevirtual

Definition at line 27 of file FPGATrackSimTrackMonitor.cxx.

{
  // print all configured properties (debugging)
  auto props = this->getProperties();
  for( Gaudi::Details::PropertyBase* prop : props ) {
    if (prop->ownerTypeName()==this->type()) {      
      ATH_MSG_DEBUG("Property:\t" << prop->name() << "\t : \t" << prop->toString());
    }
  }
  ATH_MSG_INFO("Initialize Track Monitoring");
  ATH_CHECK(m_tHistSvc.retrieve());
 
  return StatusCode::SUCCESS;
 
}

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

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()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

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()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_dir

Gaudi::Property<std::string> FPGATrackSimTrackMonitor::m_dir {this, "dir", {"/TRACKMON/"}, "String name of output directory"}

private

Definition at line 80 of file FPGATrackSimTrackMonitor.h.

80{this, "dir", {"/TRACKMON/"}, "String name of output directory"};

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_maxNElements

unsigned FPGATrackSimTrackMonitor::m_maxNElements {0}

private

Definition at line 90 of file FPGATrackSimTrackMonitor.h.

90{0};

◆ m_monTool

ToolHandle<GenericMonitoringTool> FPGATrackSimTrackMonitor::m_monTool {this,"MonTool", "", "Monitoring tool"}

private

Definition at line 77 of file FPGATrackSimTrackMonitor.h.

77{this,"MonTool", "", "Monitoring tool"};

◆ m_nElements

unsigned FPGATrackSimTrackMonitor::m_nElements {0}

private

Definition at line 87 of file FPGATrackSimTrackMonitor.h.

87{0};

◆ m_tHistSvc

ServiceHandle<ITHistSvc> FPGATrackSimTrackMonitor::m_tHistSvc {this, "THistSvc", "THistSvc/THistSvc", "Histogramming service"}

private

Definition at line 84 of file FPGATrackSimTrackMonitor.h.

84{this, "THistSvc", "THistSvc/THistSvc", "Histogramming service"};

◆ m_totNElements

unsigned FPGATrackSimTrackMonitor::m_totNElements {0}

private

Definition at line 93 of file FPGATrackSimTrackMonitor.h.

93{0};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

The documentation for this class was generated from the following files:

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ FPGATrackSimTrackMonitor()

Member Function Documentation

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ extraDeps_update_handler()

◆ fillRoad()

◆ fillTrack()

◆ getMaxNElements()

◆ getNElements()

◆ getTotNElements()

◆ initialize()

◆ inputHandles()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_detStore

◆ m_dir

◆ m_evtStore

◆ m_maxNElements

◆ m_monTool

◆ m_nElements

◆ m_tHistSvc

◆ m_totNElements

◆ m_varHandleArraysDeclared

◆ m_vhka