#include <InDetRecStatisticsAlg.h>

Inheritance diagram for InDet::InDetRecStatisticsAlg:

Classes
struct	Counter

Public Member Functions
	InDetRecStatisticsAlg (const std::string &name, ISvcLocator *pSvcLocator)
	Default Constructor.
StatusCode	initialize ()
	Initialization of services, track collections, creates TrackStatHelper for each Track Collection.
StatusCode	execute (const EventContext &ctx) const
	Calculation of statistics.
StatusCode	finalize ()
	Statistics table printed, collection cleared.
StatusCode	getServices ()
	Get various services such as StoreGate, dictionaries, detector managers etc.
StatusCode	resetStatistics ()
	Clear statistics counters, called before each track collection is processed.
void	printStatistics ()
	Print tracking statistics calculated with TrackStatHelper.
void	printTrackSummary (MsgStream &out, enum eta_region)
	Print track statistics for all and low proability tracks.
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
virtual void	setFilterPassed (bool state, const EventContext &ctx) 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	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
enum	ECounter { kN_rec_tracks_without_perigee , kN_unknown_hits , kN_rec_tracks_processed , kN_gen_tracks_processed , kN_spacepoints_processed , kNCounter }
using	CounterLocal = Counter<int>
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	selectGenSignal (const McEventCollection *, std::vector< std::pair< HepMC::ConstGenParticlePtr, int > > &, unsigned int, unsigned int, CounterLocal &counter) const
	Select charged,stable particles which pass pt and eta cuts for analysis.
void	selectRecSignal (const TrackCollection , std::vector< const Trk::Track > &, std::vector< const Trk::Track * > &, CounterLocal &counter) const
	Select for analysis reconstructed tracks passing Pt and eta cuts.
const Trk::TrackParameters *	getUnbiasedTrackParameters (const Trk::TrackParameters , const Trk::MeasurementBase )
	Get Unbiased Track Parameters.
Identifier	getIdentifier (const Trk::MeasurementBase *measurement)
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Static Private Member Functions
static float	calculatePull (const float, const float, const float)
	Calculate pull from residual, track and hit error.

Private Attributes
std::vector< class TrackStatHelper * >	m_SignalCounters
	Vector of TrackStatHelper objects, one for each track collection.
const AtlasDetectorID *	m_idHelper
	Used to find out the sub-det from PRD->identify().
const PixelID *	m_pixelID
	get pixel layer from hit ID
const SCT_ID *	m_sctID
	get sct layer from hit ID
const TRT_ID *	m_trtID
	get trt layer from hit ID
const IdDictManager *	m_idDictMgr
ToolHandle< Trk::ITruthToTrack >	m_truthToTrack
	tool to create track parameters from a gen particle
ToolHandle< Trk::IExtendedTrackSummaryTool >	m_trkSummaryTool
	tool to get track summary information from track
ToolHandle< Trk::IUpdator >	m_updatorHandle
	Tool handle of updator for unbiased states.
Trk::IUpdator *	m_updator
	updator for unbiased states
ToolHandle< Trk::IResidualPullCalculator >	m_residualPullCalculator
	The residual and pull calculator tool handle.
SG::ReadHandleKey< McEventCollection >	m_McTrackCollection_key
ToolHandle< Trk::ITrackSelectorTool >	m_trackSelectorTool
bool	m_UseTrackSummary
	Flag to print detailed statistics for each track collection.
bool	m_printSecondary
	Flag to print hit information for secondary tracks.
float	m_minPt
	Minimum Pt cut for tracks used by the algorithm.
float	m_maxEta
	Maximum Eta cut for tracks used by the algorithm.
float	m_maxEtaBarrel
	define max eta of barrel region
float	m_maxEtaTransition
	define max eta of transition region
float	m_maxEtaEndcap
	define max eta of eta region
float	m_fakeTrackCut
	Maximum probability for which a track will be considered a fake.
float	m_fakeTrackCut2
	Second definition of maximum probability for which a track will be considered a fake.
float	m_matchTrackCut
	Minimum number of hits from a truth track to be considered a matched reco track.
float	m_maxRStartPrimary
	Maximum R of start vertex to be considered primary.
float	m_maxRStartSecondary
	Maximum R of start vertex to be considered secondary.
float	m_maxZStartPrimary
	Maximum Z of start vertex to be considered primary.
float	m_maxZStartSecondary
	Maximum Z of start vertex to be considered secondary.
float	m_minREndPrimary
	If track has end vertex, this is min R of end vertex to be considered primary.
float	m_minREndSecondary
	If track has end vertex, this is min R of end vertex to be considered secondary.
float	m_minZEndPrimary
	If track has end vertex, this is min Z of end vertex to be considered primary.
float	m_minZEndSecondary
	If track has end vertex, this is min Z of end vertex to be considered secondary.
bool	m_useTrackSelection
	Use track selector tool.
bool	m_doTruth
	Use truth information.
float	m_minEtaFORWARD
float	m_maxEtaFORWARD
SG::ReadHandleKeyArray< TrackCollection >	m_RecTrackCollection_keys
SG::ReadHandleKeyArray< TrackTruthCollection >	m_TrackTruthCollection_keys
std::atomic< bool >	m_pullWarning {}
	warn only once, if pull cannot be calculated
std::atomic< bool >	m_UpdatorWarning {}
	warn only once, if unbiased track states can not be calculated
std::atomic< int >	m_isUnbiased
	if can get unbiased residuals
std::atomic< long >	m_events_processed
	number of events processed
Counter< std::atomic< long > > m_counter	ATLAS_THREAD_SAFE
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
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 46 of file InDetRecStatisticsAlg.h.

Member Typedef Documentation

◆ CounterLocal

using InDet::InDetRecStatisticsAlg::CounterLocal = Counter<int>

private

Definition at line 79 of file InDetRecStatisticsAlg.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ ECounter

enum InDet::InDetRecStatisticsAlg::ECounter

private

Enumerator
kN_rec_tracks_without_perigee	number of tracks w/o perigee
kN_unknown_hits	number of hits without track
kN_rec_tracks_processed	number of reconstructed tracks processed
kN_gen_tracks_processed	number of generated tracks processed
kN_spacepoints_processed	number of space points processed
kNCounter

Definition at line 147 of file InDetRecStatisticsAlg.h.

                    {
         kN_rec_tracks_without_perigee, 
         kN_unknown_hits,               
         kN_rec_tracks_processed,       
         kN_gen_tracks_processed,       
         kN_spacepoints_processed,      
         kNCounter
      };

Constructor & Destructor Documentation

◆ InDetRecStatisticsAlg()

InDet::InDetRecStatisticsAlg::InDetRecStatisticsAlg	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Default Constructor.

Definition at line 74 of file InDetRecStatisticsAlg.cxx.

                                                                                                 :
  AthReentrantAlgorithm(name, pSvcLocator),
  m_trtID                      (nullptr),
  m_idDictMgr                  (nullptr),
  m_truthToTrack               ("Trk::TruthToTrack"),
  m_trkSummaryTool             ("Trk::TrackSummaryTool/InDetTrackSummaryTool"),
  m_updatorHandle              ("Trk::KalmanUpdator/TrkKalmanUpdator"),
  m_updator                    (nullptr),
  m_residualPullCalculator     ("Trk::ResidualPullCalculator/ResidualPullCalculator"),
  m_McTrackCollection_key      ("TruthEvent"),
  m_trackSelectorTool          ("InDet::InDetDetailedTrackSelectorTool"),
  m_UseTrackSummary            (true),
  m_printSecondary             (false),
  m_minPt                      (1000),
  m_maxEta                     (4.2),
  m_maxEtaBarrel               (0.8),
  m_maxEtaTransition           (1.6),
  m_maxEtaEndcap               (2.5),
  m_fakeTrackCut               (0.9),
  m_fakeTrackCut2              (0.7),
  m_matchTrackCut              (0.5),
  m_maxRStartPrimary           (  25.0*CLHEP::mm),
  m_maxRStartSecondary         ( 360.0*CLHEP::mm),
  m_maxZStartPrimary           ( 200.0*CLHEP::mm),
  m_maxZStartSecondary         (2000.0*CLHEP::mm),
  m_minREndPrimary             ( 400.0*CLHEP::mm),
  m_minREndSecondary           (1000.0*CLHEP::mm),
  m_minZEndPrimary             (2300.0*CLHEP::mm),
  //m_maxZIndet                (),
  m_minZEndSecondary           (3200.0*CLHEP::mm),
  m_useTrackSelection          (false),
  m_doTruth                    (true),
  m_minEtaFORWARD              (2.5),
  m_maxEtaFORWARD              (4.2),
  m_isUnbiased                 (0),
  m_events_processed           (0)
{
  // m_RecTrackCollection_keys.push_back(std::string("Tracks"));
  // m_TrackTruthCollection_keys.push_back(std::string("TrackTruthCollection"));
 
  // Algorithm properties
  declareProperty("SummaryTool",                m_trkSummaryTool);
  declareProperty("TruthToTrackTool",           m_truthToTrack);
  declareProperty("UpdatorTool",                m_updatorHandle,
          "Measurement updator to calculate unbiased track states");
  declareProperty("ResidualPullCalculatorTool", m_residualPullCalculator,
          "Tool to calculate residuals and pulls");
  declareProperty("TrackCollectionKeys",        m_RecTrackCollection_keys);
  declareProperty("McTrackCollectionKey",       m_McTrackCollection_key);
  declareProperty("TrackTruthCollectionKeys",   m_TrackTruthCollection_keys);
  declareProperty("UseTrackSelection"       ,   m_useTrackSelection);
  declareProperty("DoTruth"                 ,   m_doTruth);
  declareProperty("TrackSelectorTool"       ,   m_trackSelectorTool);
  declareProperty("UseTrackSummary",            m_UseTrackSummary);
  declareProperty("PrintSecondary",             m_printSecondary);
  declareProperty("minPt",               m_minPt);
  declareProperty("maxEta",              m_maxEta);
  declareProperty("maxEtaBarrel",               m_maxEtaBarrel  );
  declareProperty("maxEtaTransition",           m_maxEtaTransition);
  declareProperty("maxEtaEndcap",               m_maxEtaEndcap);
  declareProperty("maxEtaFORWARD", m_maxEtaFORWARD);
  declareProperty("minEtaFORWARD", m_minEtaFORWARD);
  declareProperty("fakeTrackCut",               m_fakeTrackCut);
  declareProperty("fakeTrackCut2",              m_fakeTrackCut2);
  declareProperty("matchTrackCut",              m_matchTrackCut);
  declareProperty("maxRStartPrimary",            m_maxRStartPrimary);
  declareProperty("maxRStartSecondary",          m_maxRStartSecondary);
  declareProperty("maxZStartPrimary",            m_maxZStartPrimary);
  declareProperty("maxZStartSecondary",         m_maxZStartSecondary);
  declareProperty("minREndPrimary",      m_minREndPrimary);
  declareProperty("minREndSecondary",        m_minREndSecondary);
  declareProperty("minZEndPrimary",      m_minZEndPrimary);
  declareProperty("minZEndSecondary",            m_minZEndSecondary);
  m_idHelper  = nullptr;
  m_pixelID   = nullptr;
  m_sctID     = nullptr;
  m_UpdatorWarning = false;
  m_pullWarning = false;
}

Member Function Documentation

◆ calculatePull()

float InDet::InDetRecStatisticsAlg::calculatePull	(	const float	residual,
		const float	trkErr,
		const float	hitErr )

staticprivate

Calculate pull from residual, track and hit error.

Definition at line 771 of file InDetRecStatisticsAlg.cxx.

                                                 {
  double ErrorSum;
  ErrorSum = sqrt(pow(trkErr, 2) + pow(hitErr, 2));
  if (ErrorSum != 0) { return residual/ErrorSum; }
  else { return 0; }
}

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const

overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

{
  return 0;
}

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ execute()

StatusCode InDet::InDetRecStatisticsAlg::execute ( const EventContext & ctx ) const

Calculation of statistics.

Definition at line 265 of file InDetRecStatisticsAlg.cxx.

                                                                             {
 
    ATH_MSG_DEBUG("entering execute()");
 
    // Get reconstructed tracks , generated tracks, and truth from storegate
 
    SG::ReadHandle<McEventCollection> SimTracks;
 
    if (m_doTruth) {
      SimTracks=SG::ReadHandle<McEventCollection>(m_McTrackCollection_key,ctx);
      if (!SimTracks.isValid()) {
        // @TODO warning ?
        ATH_MSG_WARNING("Error retrieving collections !");
        return StatusCode::SUCCESS;
      }
    }
 
    // Doesn't take account of pileup:
    //m_gen_tracks_processed += (*(SimTracks->begin()))->particles_size();
    m_events_processed ++;
    CounterLocal counter;
 
    // select charged and stable generated tracks
    // apply pt, eta etc cuts to generated tracks
    // devide generated tracks into primary, truncated, secondary
 
    std::vector <std::pair<HepMC::ConstGenParticlePtr,int> > GenSignal;
    //     GenSignalPrimary, GenSignalTruncated, GenSignalSecondary;
    unsigned int inTimeStart = 0;
    unsigned int inTimeEnd   = 0;
    if (m_doTruth) selectGenSignal ((SimTracks.isValid() ? &(*SimTracks) : nullptr), GenSignal, inTimeStart, inTimeEnd, counter);
 
    // step through the various reconstructed TrackCollections and
    // corresponding TrackTruthCollections and produce statistics for each
 
    if (m_SignalCounters.empty()) {
      ATH_MSG_ERROR("No reco track collection specified! Aborting.");
      return StatusCode::FAILURE;
    }
 
    std::vector< SG::ReadHandle<TrackCollection> > rec_track_collections = m_RecTrackCollection_keys.makeHandles(ctx);
    std::vector< SG::ReadHandle<TrackTruthCollection> > truth_track_collections;
    if (m_doTruth && !m_TrackTruthCollection_keys.empty()) {
      truth_track_collections = m_TrackTruthCollection_keys.makeHandles(ctx);
      if (truth_track_collections.size() != rec_track_collections.size()) {
        ATH_MSG_ERROR("Different number of reco and truth track collections (" << rec_track_collections.size() << "!=" << truth_track_collections.size() << ")" );
      }
    }
    if (m_SignalCounters.size() != rec_track_collections.size()) {
        ATH_MSG_ERROR("Number expected reco track collections does not match the actual number of such collections ("
                      << m_SignalCounters.size() << "!=" << rec_track_collections.size() << ")" );
    }
 
    std::vector< SG::ReadHandle<TrackCollection> >::iterator rec_track_collections_iter = rec_track_collections.begin();
    std::vector< SG::ReadHandle<TrackTruthCollection> >::iterator truth_track_collections_iter = truth_track_collections.begin();
    for (std::vector <class TrackStatHelper *>::const_iterator statHelper
     =  m_SignalCounters.begin();
     statHelper !=  m_SignalCounters.end();
         ++statHelper, ++rec_track_collections_iter) {
      assert( rec_track_collections_iter != rec_track_collections.end());
 
      ATH_MSG_DEBUG("Acessing TrackCollection " <<  m_RecTrackCollection_keys.at(rec_track_collections_iter - rec_track_collections.begin()).key());
      const TrackCollection       * RecCollection = &(**rec_track_collections_iter);
      const TrackTruthCollection  * TruthMap  = nullptr;
 
      if (RecCollection)  ATH_MSG_DEBUG("Retrieved " << RecCollection->size() << " reconstructed tracks from storegate");
 
      if (m_doTruth) {
        ATH_MSG_DEBUG("Acessing TrackTruthCollection " <<  m_TrackTruthCollection_keys.at(truth_track_collections_iter - truth_track_collections.begin()).key());
        assert( truth_track_collections_iter != truth_track_collections.end());
        TruthMap = &(**truth_track_collections_iter);
        if (TruthMap)   ATH_MSG_DEBUG("Retrieved " << TruthMap->size() << " TrackTruth elements from storegate");
        ++truth_track_collections_iter;
      }
 
      //start process of getting correct track summary
 
      std::vector <const Trk::Track *>          RecTracks, RecSignal;
      selectRecSignal                     (RecCollection, RecTracks,RecSignal,counter);
 
      ATH_MSG_DEBUG(
            "  RecTracks.size()="          << RecTracks.size()
            << ", GenSignal.size()="          << GenSignal.size());
 
      ATH_MSG_DEBUG("Accumulating Statistics...");
      (*statHelper)->addEvent    (RecCollection,
                  RecTracks,
                  GenSignal,
                  TruthMap,
                  m_idHelper,
                  m_pixelID,
                  m_sctID,
                  m_trkSummaryTool.operator->(),
                  m_UseTrackSummary,
                  &inTimeStart,
                  &inTimeEnd);
 
      counter.m_counter[kN_rec_tracks_processed] += RecCollection->size();
 
      for (  TrackCollection::const_iterator it = RecCollection->begin() ;
         it < RecCollection->end(); ++ it){
    std::vector<const Trk::RIO_OnTrack*> rioOnTracks;
    Trk::RoT_Extractor::extract( rioOnTracks,
                     (*it)->measurementsOnTrack()->stdcont() );
    counter.m_counter[kN_spacepoints_processed] += rioOnTracks.size();
      }
 
    }
    m_counter += counter;
 
    ATH_MSG_DEBUG("leaving execute()");
    return StatusCode::SUCCESS;
}

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

◆ filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

◆ finalize()

StatusCode InDet::InDetRecStatisticsAlg::finalize ( )

Statistics table printed, collection cleared.

Definition at line 382 of file InDetRecStatisticsAlg.cxx.

                                                    {
 
  // Part 1: Get the messaging service, print where you are
  ATH_MSG_DEBUG("finalize()");
 
  printStatistics();
 
  for (std::vector <class TrackStatHelper *>::const_iterator collection =
       m_SignalCounters.begin(); collection !=  m_SignalCounters.end();
       ++collection) {
    ATH_MSG_DEBUG(s_linestr2);
    delete (*collection);
  }
  m_SignalCounters.clear();
  return StatusCode::SUCCESS;
}

◆ getIdentifier()

Identifier InDet::InDetRecStatisticsAlg::getIdentifier ( const Trk::MeasurementBase * measurement )

private

Definition at line 821 of file InDetRecStatisticsAlg.cxx.

                                                                                           {
  Identifier id;
  const Trk::CompetingRIOsOnTrack *comprot = nullptr;
  // identify by ROT:
  const Trk::RIO_OnTrack *rot =
    dynamic_cast<const Trk::RIO_OnTrack*>(measurement);
  if (rot) {
    id = rot->identify();
  } else {
    // identify by CompetingROT:
    comprot = dynamic_cast<const Trk::CompetingRIOsOnTrack*>(measurement);
    if (comprot) {
      rot = &comprot->rioOnTrack(comprot->indexOfMaxAssignProb());
      id = rot->identify();
    } else {
      ATH_MSG_DEBUG("measurement is neither ROT nor competingROT:"
      <<" can not determine detector type");
      id.clear();
    }
  }
  delete comprot;
  return id;
}

◆ getServices()

StatusCode InDet::InDetRecStatisticsAlg::getServices ( )

Get various services such as StoreGate, dictionaries, detector managers etc.

Definition at line 400 of file InDetRecStatisticsAlg.cxx.

{
    //Set up ATLAS ID helper to be able to identify the RIO's det-subsystem.
 
    // Get the dictionary manager from the detector store
    const IdDictManager*  idDictMgr = nullptr;
    StatusCode sc = detStore()->retrieve(idDictMgr, "IdDict");
    if (sc.isFailure()) {
      ATH_MSG_FATAL("Could not get IdDictManager !");
      return StatusCode::FAILURE;
    }
 
    // Initialize the helper with the dictionary information.
    sc = detStore()->retrieve(m_idHelper, "AtlasID");
    if (sc.isFailure()) {
      ATH_MSG_FATAL("Could not get AtlasDetectorID helper.");
      return StatusCode::FAILURE;
    }
 
   //get Pixel, SCT, TRT managers and helpers
 
   if (detStore()->retrieve(m_pixelID, "PixelID").isFailure()) {
     msg(MSG::FATAL) << "Could not get Pixel ID helper" << endmsg;
     return StatusCode::FAILURE;
   }
   if (detStore()->retrieve(m_sctID, "SCT_ID").isFailure()) {
     msg(MSG::FATAL) << "Could not get SCT ID helper" << endmsg;
     return StatusCode::FAILURE;
   }
 
   //retrieve the TRT helper only if not-SLHC layout used
   sc = detStore()->retrieve(m_idDictMgr, "IdDict");
   if (sc.isFailure()) {
     ATH_MSG_FATAL("Could not get IdDictManager !");
     return StatusCode::FAILURE;
   }
   const IdDictDictionary* dict = m_idDictMgr->manager()->find_dictionary("InnerDetector");
   if(!dict) {
     ATH_MSG_FATAL(" Cannot access InnerDetector dictionary ");
     return StatusCode::FAILURE;
   }
 
   bool isSLHC = false;
   if (dict->file_name().find("SLHC")!=std::string::npos) isSLHC=true;
 
   if(!isSLHC){
     if (detStore()->retrieve(m_trtID, "TRT_ID").isFailure()) {
       msg(MSG::FATAL) << "Could not get TRT ID helper" << endmsg;
       return StatusCode::FAILURE;
     }
   }
   //
 
   if (m_UseTrackSummary) {
     if (m_trkSummaryTool.retrieve().isFailure() ) {
       ATH_MSG_FATAL("Failed to retrieve tool "
       << m_trkSummaryTool);
       return StatusCode::FAILURE;
     } else {
       ATH_MSG_INFO("Retrieved tool " << m_trkSummaryTool);
     }
   } else {
     m_trkSummaryTool.disable();
   }
 
   // AG: init truthToTrack
   if (m_doTruth) {
     if (m_truthToTrack.retrieve().isFailure() ) {
       ATH_MSG_FATAL("Failed to retrieve tool " << m_truthToTrack);
       return StatusCode::FAILURE;
     } else {
       ATH_MSG_INFO("Retrieved tool " << m_truthToTrack);
     }
   } else {
     m_truthToTrack.disable();
   }
 
   //adding track selector tool
   if(m_useTrackSelection){
     if ( m_trackSelectorTool.retrieve().isFailure() ) {
       ATH_MSG_FATAL("Failed to retrieve tool " << m_trackSelectorTool);
       return StatusCode::FAILURE;
     } else {
       ATH_MSG_INFO("Retrieved tool " << m_trackSelectorTool);
     }
   } else {
     m_trackSelectorTool.disable();
   }
   return StatusCode :: SUCCESS;
}

◆ getUnbiasedTrackParameters()

const Trk::TrackParameters * InDet::InDetRecStatisticsAlg::getUnbiasedTrackParameters	(	const Trk::TrackParameters *	trkParameters,
		const Trk::MeasurementBase *	measurement )

private

Get Unbiased Track Parameters.

Definition at line 780 of file InDetRecStatisticsAlg.cxx.

                                                                                                                                                                 {
 
 
  const Trk::TrackParameters *unbiasedTrkParameters = nullptr;
 
 // -----------------------------------------
  // use unbiased track states or normal ones?
  // unbiased track parameters are tried to retrieve if the updator tool
  //    is available and if unbiased track states could be produced before
  //    for the current track (ie. if one trial to get unbiased track states
  //    fail
 
  if (m_updator && (m_isUnbiased==1) ) {
    if ( trkParameters->covariance() ) {
      // Get unbiased state
      ATH_MSG_VERBOSE(" getting unbiased params");
      unbiasedTrkParameters =
    m_updator->removeFromState( *trkParameters,
                    measurement->localParameters(),
                    measurement->localCovariance()).release();
 
      if (!unbiasedTrkParameters) {
    ATH_MSG_WARNING("Could not get unbiased track parameters, "
        <<"use normal parameters");
    m_isUnbiased = 0;
      }
    } else if(!m_UpdatorWarning) {
      // warn only once!
      ATH_MSG_WARNING("TrackParameters contain no covariance: "
      <<"Unbiased track states can not be calculated "
      <<"(ie. pulls and residuals will be too small)");
      m_UpdatorWarning = true;
      m_isUnbiased = 0;
    } else {
      m_isUnbiased = 0;
    }
  } // end if no measured track parameter
  return unbiasedTrkParameters;
}

◆ initialize()

StatusCode InDet::InDetRecStatisticsAlg::initialize ( )

Initialization of services, track collections, creates TrackStatHelper for each Track Collection.

Definition at line 156 of file InDetRecStatisticsAlg.cxx.

                                                 {
 
  // Part 1: Get the messaging service, print where you are
  ATH_MSG_DEBUG("initialize()");
 
  StatusCode sc1 = getServices();           // retrieve store gate service etc
  if (sc1.isFailure()) {
    ATH_MSG_FATAL("Error retrieving services !");
    return StatusCode::FAILURE;
  }
 
  if (m_RecTrackCollection_keys.empty()) {
    ATH_MSG_ERROR("No reco track collection specified! Aborting.");
    return StatusCode::FAILURE;
  }
 
  if (m_doTruth && m_RecTrackCollection_keys.size() != m_TrackTruthCollection_keys.size()) {
    ATH_MSG_ERROR("You have specified "
          << m_RecTrackCollection_keys.size()
          << " TrackCollection keys, and " <<  m_TrackTruthCollection_keys.size()
          << " TrackTruthCollection keys."
          << " You have to specify one TrackTruthCollection for each"
          << " TrackCollection! Exiting."
          );
    return StatusCode::FAILURE;
  }
 
  // ----------------------------------
  // use updator to get unbiased states
  if ( ! m_updatorHandle.empty() ) {
    if (m_updatorHandle.retrieve().isFailure()) {
      ATH_MSG_FATAL("Could not retrieve measurement updator tool: "
      << m_updatorHandle);
      return StatusCode::FAILURE;
    }
    m_updator = &(*m_updatorHandle);
  } else {
    ATH_MSG_DEBUG(
      "No Updator for unbiased track states given, use normal states!");
    m_updator = nullptr;
  }
 
 
  //get residual and pull calculator
  if (m_residualPullCalculator.empty()) {
    ATH_MSG_INFO(
      "No residual/pull calculator for general hit residuals configured."
    );
    ATH_MSG_INFO(
      "It is recommended to give R/P calculators to the det-specific tool"
    << " handle lists then.");
  } else if (m_residualPullCalculator.retrieve().isFailure()) {
    ATH_MSG_FATAL("Could not retrieve "<< m_residualPullCalculator
    <<" (to calculate residuals and pulls) ");
 
  } else {
    ATH_MSG_INFO("Generic hit residuals&pulls will be calculated in one or both "
         << "available local coordinates");
  }
 
  // create one TrackStatHelper object of each trackCollection --- this is used to accumulate track and hit statistics
 
  struct cuts ct;
  ct.maxEtaBarrel=  m_maxEtaBarrel;
  ct.maxEtaTransition= m_maxEtaTransition;
  ct.maxEtaEndcap=        m_maxEtaEndcap;
  ct.fakeTrackCut=      m_fakeTrackCut;
  ct.fakeTrackCut2=     m_fakeTrackCut2;
  ct.matchTrackCut    =  m_matchTrackCut;
  ct.maxRStartPrimary  = m_maxRStartPrimary;
  ct.maxRStartSecondary    = m_maxRStartSecondary;
  ct.maxZStartPrimary  = m_maxZStartPrimary;
  ct.maxZStartSecondary    = m_maxZStartSecondary;
  ct.minREndPrimary    = m_minREndPrimary;
  ct.minREndSecondary  = m_minREndSecondary;
  ct.minZEndPrimary    = m_minZEndPrimary;
  ct.minZEndSecondary      = m_minZEndSecondary;
  ct.minPt               = m_minPt;
  ct.minEtaFORWARD = m_minEtaFORWARD;
  ct.maxEtaFORWARD = m_maxEtaFORWARD;
 
  unsigned int nCollections = 0;
  for (SG::ReadHandleKeyArray<TrackCollection>::const_iterator
       it = m_RecTrackCollection_keys.begin();
       it < m_RecTrackCollection_keys.end(); ++ it) {
    InDet::TrackStatHelper * collection =
      new TrackStatHelper(it->key(),(m_doTruth ? m_TrackTruthCollection_keys[nCollections].key() : ""), m_doTruth);
    nCollections ++;
    collection->SetCuts(ct);
    m_SignalCounters.push_back(collection);
  }
 
  StatusCode sc3 = resetStatistics();     // reset all statistic counters
  if (sc3.isFailure()) {
    ATH_MSG_FATAL("Error in resetStatistics !");
    return StatusCode::FAILURE;
  }
 
  ATH_CHECK( m_RecTrackCollection_keys.initialize() );
  ATH_CHECK( m_McTrackCollection_key.initialize(m_doTruth && !m_McTrackCollection_key.key().empty()) );
  ATH_CHECK( m_TrackTruthCollection_keys.initialize() );
 
  return StatusCode :: SUCCESS;
 
}

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

◆ isClonable()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ITkPixelCablingAlg, ITkStripCablingAlg, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, and Simulation::BeamEffectsAlg.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

{
  // Reentrant algorithms are clonable.
  return true;
}

◆ msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::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< Gaudi::Algorithm > >::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.

◆ printStatistics()

void InDet::InDetRecStatisticsAlg::printStatistics ( )

Print tracking statistics calculated with TrackStatHelper.

Definition at line 583 of file InDetRecStatisticsAlg.cxx.

                                                     {
  if (!msgLvl(MSG::INFO)) return;
 
  ATH_MSG_INFO(" ********** Beginning InDetRecStatistics Statistics Table ***********");
  ATH_MSG_INFO("For documentation see https://twiki.cern.ch/twiki/bin/view/Atlas/InDetRecStatistics");
  ATH_MSG_INFO("(or for guaranteed latest version: http://atlas-sw.cern.ch/cgi-bin/viewcvs-atlas.cgi/offline/InnerDetector/InDetValidation/InDetRecStatistics/doc/mainpage.h?&view=markup )");
  ATH_MSG_INFO(" ********************************************************************");
 
  std::stringstream outstr;
  int def_precision(outstr.precision());
  outstr << "\n"
         << MSG::INFO
     << std::setiosflags(std::ios::fixed | std::ios::showpoint)
     << std::setw(7) << std::setprecision(2)
     << s_linestr << "\n"
         << "Summary" << "\n"
         << "\tProcessed              : " << m_events_processed
     << " events, "  << m_counter.m_counter[kN_rec_tracks_processed]
     << " reconstructed tracks with " << m_counter.m_counter[kN_spacepoints_processed]
     << " hits, and "  << m_counter.m_counter[kN_gen_tracks_processed]
     << " truth particles" << "\n"
         << "\tProblem objects        : " <<  m_counter.m_counter[kN_rec_tracks_without_perigee]
     << " tracks without perigee, "
     << m_counter.m_counter[kN_unknown_hits] << " unknown hits" << "\n"
         << "\t" << "Reco  TrackCollections : ";
  bool first = true;
  for (std::vector <class TrackStatHelper *>::const_iterator collection =
     m_SignalCounters.begin();
       collection !=  m_SignalCounters.end(); ++collection)
    {
      if (first) {
        first = false;
      }
      else {
    outstr << ", ";
      }
      outstr << "\"" << (*collection)->key() << "\"";
    }
  ATH_MSG_INFO(outstr.str());
  outstr.str("");
 
  if (m_doTruth)
    {
      outstr.str("");
      outstr << "\n"
             << "\t" << "TrackTruthCollections  : ";
      first = true;
      for (std::vector <class TrackStatHelper *>::const_iterator collection =  m_SignalCounters.begin();
       collection !=  m_SignalCounters.end(); ++collection)
    {
      if (first) {
        first = false;
      }
      else {
        outstr << ", ";
      }
      outstr << "\"" << (*collection)->Truthkey() << "\"";
    }
      ATH_MSG_INFO(outstr.str());
      outstr.str("");
    }
   outstr.str("");
   outstr << "\n"
      << s_linestr2 << "\n"
      << "Cuts and Settings for Statistics Table" << "\n"
      << "\t" << "TrackSummary Statistics" << "\t"
      << (m_UseTrackSummary     ? "YES" : "NO") << "\n"
      << "\t" << "Signal                \t" << "pT > "
      << m_minPt/1000 << " GeV/c, |eta| < " << m_maxEta << "\t\t"
      << "\t" << "Primary track start   \t" << "R < "
      << m_maxRStartPrimary << "mm and |z| < "
      << m_maxZStartPrimary << "mm" << "\n"
      << "\t" << "Barrel                \t" << 0.0
      << "< |eta| < " << m_maxEtaBarrel     << "\t\t\t"
      << "\t" << "Primary track end     \t" << "R > "
      << m_minREndPrimary   << "mm or |z| > " << m_minZEndPrimary
      << "mm" << "\n"
      << "\t" << "Transition Region     \t" << m_maxEtaBarrel
      << "< |eta| < " << m_maxEtaTransition << "\t\t\t"
      << "\t" << "Secondary (non-Primary) start \t"
      << " R < "    << m_maxRStartSecondary << "mm and"
      << " |z| < "  << m_maxZStartSecondary << " mm" << "\n"
      << "\t" << "Endcap                \t" << m_maxEtaTransition
      << "< |eta| < " << m_maxEtaEndcap     << "\t\t\t"
      << "\t" << "Secondary (non-primary) end   \t"
      << " R > "    << m_minREndSecondary   << "mm or"
      << " |z| > "  << m_minREndSecondary   << "mm" << "\n"
          << "\t" << "Forward                \t"
          << "|eta| > " << m_minEtaFORWARD     << "\n"
      << "\t" << "Low prob tracks #1    \t" << "< "
      << m_fakeTrackCut  << " of hits from single Truth Track "
      << "\n"
      << "\t" << "Low prob tracks #2    \t" << "< "
      << m_fakeTrackCut2 << " of hits from single Truth Track "
      << "\n"
      << "\t" << "No link tracks        \t  Track has no link associated to an HepMC Particle" << "\n"
      << "\t" << "Good reco tracks      \t" << "> "
      << m_matchTrackCut << " of hits from single Truth Track + a link !";
  ATH_MSG_INFO(outstr.str());
  outstr.str("");
 
  MsgStream &out = msg(MSG::INFO);
  {
  RestoreStream<MsgStream> restore(out);
  out << "\n" << s_linestr2 << "\n";
  m_SignalCounters.back()->print(out);
 
  if (m_UseTrackSummary) {
    std::string track_stummary_type_header = TrackStatHelper::getSummaryTypeHeader();
    out << "\n"
        << s_linestr2 << "\n"
        << "Detailed Statistics for Hits on Reconstructed tracks, using TrackSummary: (Preselection of tracks as described above.)" << "\n"
        << s_linestr2 << "\n"
        << "----------------------------------------------------------------------------------------------------------------------------------------------------" << "\n"
        << "  Reco Tracks                           .........................................hits/track.......................................................  " << "\n"
        << "----------------------------------------------------------------------------------------------------------------------------------------------------" << "\n"
        << "  in BARREL                tracks/event " << track_stummary_type_header << "\n"
        << "----------------------------------------------------------------------------------------------------------------------------------------------------" << "\n";
    printTrackSummary (out, ETA_BARREL);
 
    out<< "\n"
       << "----------------------------------------------------------------------------------------------------------------------------------------------------" << "\n"
       << "  in TRANSITION region     tracks/event  " << track_stummary_type_header << "\n"
       << "----------------------------------------------------------------------------------------------------------------------------------------------------"
       << "\n";
    printTrackSummary (out, ETA_TRANSITION);
 
    out << "\n"
        << "----------------------------------------------------------------------------------------------------------------------------------------------------" << "\n"
        << "  in ENDCAP                tracks/event " << track_stummary_type_header << "\n"
        << "----------------------------------------------------------------------------------------------------------------------------------------------------" << "\n";
    printTrackSummary (out, ETA_ENDCAP);
 
    out << "\n"
        << "----------------------------------------------------------------------------------------------------------------------------------------------------" << "\n"
        << "  in FORWARD region        tracks/event  " << track_stummary_type_header << "\n"
        << "----------------------------------------------------------------------------------------------------------------------------------------------------" << "\n";
    printTrackSummary (out, ETA_FORWARD);
  }
 
  if(m_printSecondary){
    outstr.str("");
    outstr << "\n" << std::setprecision(def_precision)
           <<s_linestr<<"\n"
           <<"Statistics for Secondaries (non-Primaries)"<<"\n"
           << "\t" << "Secondary track start \t"
       << " R < "   << m_maxRStartSecondary << "mm and"
       << " |z| < " << m_maxZStartSecondary << " mm" << "\n"
           << "\t" << "Secondary track end   \t"
       << " R > "    << m_minREndSecondary << "mm or"
           << " |z| > "  << m_minZEndSecondary << "mm";
    ATH_MSG_INFO(outstr.str());
    outstr.str("");
    out << "\n" << s_linestr2 << "\n";
    m_SignalCounters.back()->printSecondary(out);
 
  }
  }
  out << endmsg;
 
  ATH_MSG_INFO(" ********** Ending InDetRecStatistics Statistics Table ***********");
  ATH_MSG_INFO( "\n"
               << s_linestr );
}

◆ printTrackSummary()

void InDet::InDetRecStatisticsAlg::printTrackSummary	(	MsgStream &	out,
		enum eta_region	eta_reg )

Print track statistics for all and low proability tracks.

Definition at line 749 of file InDetRecStatisticsAlg.cxx.

{
  bool printed = m_SignalCounters.back()->printTrackSummaryRegion(out, TRACK_ALL, eta_reg);
 
  if (printed) {
     out <<  "\n"
         << "----------------------------------------------------------------------------------------------------------------------------------------------" << "\n";
  }
 
  printed = m_SignalCounters.back()->printTrackSummaryRegion(out, TRACK_LOWTRUTHPROB, eta_reg);
  if (printed) {
     out << "\n"
         << "----------------------------------------------------------------------------------------------------------------------------------------------" << "\n";
  }
 
  m_SignalCounters.back()->printTrackSummaryRegion(out, TRACK_LOWTRUTHPROB2, eta_reg);
 
}

◆ 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< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

◆ resetStatistics()

StatusCode InDet::InDetRecStatisticsAlg::resetStatistics ( )

Clear statistics counters, called before each track collection is processed.

Definition at line 491 of file InDetRecStatisticsAlg.cxx.

                                                           {
    m_counter.reset();
    m_events_processed           = 0;
 
    for (std::vector<InDet::TrackStatHelper *>::const_iterator counter =
       m_SignalCounters.begin();
     counter != m_SignalCounters.end(); ++ counter) {
      (*counter)->reset();
    }
    return StatusCode :: SUCCESS;
}

◆ selectGenSignal()

void InDet::InDetRecStatisticsAlg::selectGenSignal	(	const McEventCollection *	SimTracks,
		std::vector< std::pair< HepMC::ConstGenParticlePtr, int > > &	GenSignal,
		unsigned int	,
		unsigned int	,
		CounterLocal &	counter ) const

private

Select charged,stable particles which pass pt and eta cuts for analysis.

Classify truth tracks as primary, secondary or truncated

Definition at line 529 of file InDetRecStatisticsAlg.cxx.

{
  if (! SimTracks) return;
 
  unsigned int nb_mc_event = SimTracks->size();
  std::unique_ptr<PileUpType>  put = std::make_unique<PileUpType>(SimTracks);
 
  McEventCollection::const_iterator inTimeMBend;
  McEventCollection::const_iterator inTimeMBbegin;
 
  if (put)
    {
      inTimeMBbegin = put->in_time_minimum_bias_event_begin();
      inTimeMBend = put->in_time_minimum_bias_event_end();
    }
 
  for(unsigned int ievt=0; ievt<nb_mc_event; ++ievt)
    {
      const HepMC::GenEvent* genEvent = SimTracks->at(ievt);
      counter.m_counter[kN_gen_tracks_processed] += genEvent->particles_size();
      for (const auto& particle: *genEvent){
      // require stable particle from generation or simulation
      if (!MC::isStable(particle)) continue;
      int   pdgCode = particle->pdg_id();
          if (MC::isNucleus(pdgCode)) continue; // ignore nuclei from hadronic interactions
      float charge = MC::charge(pdgCode);
      if (std::abs(charge)<0.5) continue;
      if (std::abs(particle->momentum().perp()) >  m_minPt  &&
          std::abs(particle->momentum().pseudoRapidity()) < m_maxEta ) {
        std::pair<HepMC::ConstGenParticlePtr,int> thisPair(particle,ievt);
        GenSignal.push_back(thisPair);
      }
    } // End of a particle iteration
    } // End of one GenEvent iteration
  }

◆ selectRecSignal()

void InDet::InDetRecStatisticsAlg::selectRecSignal	(	const TrackCollection *	RecCollection,
		std::vector< const Trk::Track * > &	RecTracks,
		std::vector< const Trk::Track * > &	RecSignal,
		CounterLocal &	counter ) const

private

Select for analysis reconstructed tracks passing Pt and eta cuts.

Definition at line 503 of file InDetRecStatisticsAlg.cxx.

                                                                                                          {
 
  for (  TrackCollection::const_iterator it = RecCollection->begin() ;
     it != RecCollection->end(); ++ it){
    RecTracks.push_back(*it);
    const DataVector<const Trk::TrackParameters>* trackpara =
      (*it)->trackParameters();
 
    if(!trackpara->empty()){
      const Trk::TrackParameters* para = trackpara->front();
      if (para){
    if (para->pT() >  m_minPt && std::abs(para->eta()) < m_maxEta) {
      RecSignal.push_back(*it);
    }
      }
    }
    else {
      counter.m_counter[kN_rec_tracks_without_perigee] ++;
    }
  }
 }

◆ setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed	(	bool	state,
		const EventContext &	ctx ) const

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

◆ sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

◆ sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

                                                               {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }      
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;  
}

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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

◆ ATLAS_THREAD_SAFE

Counter<std::atomic<long> > m_counter InDet::InDetRecStatisticsAlg::ATLAS_THREAD_SAFE

mutableprivate

Definition at line 170 of file InDetRecStatisticsAlg.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_doTruth

bool InDet::InDetRecStatisticsAlg::m_doTruth

private

Use truth information.

Definition at line 135 of file InDetRecStatisticsAlg.h.

◆ m_events_processed

std::atomic<long> InDet::InDetRecStatisticsAlg::m_events_processed

mutableprivate

number of events processed

Definition at line 145 of file InDetRecStatisticsAlg.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_fakeTrackCut

float InDet::InDetRecStatisticsAlg::m_fakeTrackCut

private

Maximum probability for which a track will be considered a fake.

Definition at line 122 of file InDetRecStatisticsAlg.h.

◆ m_fakeTrackCut2

float InDet::InDetRecStatisticsAlg::m_fakeTrackCut2

private

Second definition of maximum probability for which a track will be considered a fake.

Definition at line 123 of file InDetRecStatisticsAlg.h.

◆ m_idDictMgr

const IdDictManager* InDet::InDetRecStatisticsAlg::m_idDictMgr

private

Definition at line 103 of file InDetRecStatisticsAlg.h.

◆ m_idHelper

const AtlasDetectorID* InDet::InDetRecStatisticsAlg::m_idHelper

private

Used to find out the sub-det from PRD->identify().

Definition at line 98 of file InDetRecStatisticsAlg.h.

◆ m_isUnbiased

std::atomic<int> InDet::InDetRecStatisticsAlg::m_isUnbiased

mutableprivate

if can get unbiased residuals

Definition at line 143 of file InDetRecStatisticsAlg.h.

◆ m_matchTrackCut

float InDet::InDetRecStatisticsAlg::m_matchTrackCut

private

Minimum number of hits from a truth track to be considered a matched reco track.

Definition at line 124 of file InDetRecStatisticsAlg.h.

◆ m_maxEta

float InDet::InDetRecStatisticsAlg::m_maxEta

private

Maximum Eta cut for tracks used by the algorithm.

Definition at line 118 of file InDetRecStatisticsAlg.h.

◆ m_maxEtaBarrel

float InDet::InDetRecStatisticsAlg::m_maxEtaBarrel

private

define max eta of barrel region

Definition at line 119 of file InDetRecStatisticsAlg.h.

◆ m_maxEtaEndcap

float InDet::InDetRecStatisticsAlg::m_maxEtaEndcap

private

define max eta of eta region

Definition at line 121 of file InDetRecStatisticsAlg.h.

◆ m_maxEtaFORWARD

float InDet::InDetRecStatisticsAlg::m_maxEtaFORWARD

private

Definition at line 137 of file InDetRecStatisticsAlg.h.

◆ m_maxEtaTransition

float InDet::InDetRecStatisticsAlg::m_maxEtaTransition

private

define max eta of transition region

Definition at line 120 of file InDetRecStatisticsAlg.h.

◆ m_maxRStartPrimary

float InDet::InDetRecStatisticsAlg::m_maxRStartPrimary

private

Maximum R of start vertex to be considered primary.

Definition at line 126 of file InDetRecStatisticsAlg.h.

◆ m_maxRStartSecondary

float InDet::InDetRecStatisticsAlg::m_maxRStartSecondary

private

Maximum R of start vertex to be considered secondary.

Definition at line 127 of file InDetRecStatisticsAlg.h.

◆ m_maxZStartPrimary

float InDet::InDetRecStatisticsAlg::m_maxZStartPrimary

private

Maximum Z of start vertex to be considered primary.

Definition at line 128 of file InDetRecStatisticsAlg.h.

◆ m_maxZStartSecondary

float InDet::InDetRecStatisticsAlg::m_maxZStartSecondary

private

Maximum Z of start vertex to be considered secondary.

Definition at line 129 of file InDetRecStatisticsAlg.h.

◆ m_McTrackCollection_key

SG::ReadHandleKey<McEventCollection> InDet::InDetRecStatisticsAlg::m_McTrackCollection_key

private

Definition at line 112 of file InDetRecStatisticsAlg.h.

◆ m_minEtaFORWARD

float InDet::InDetRecStatisticsAlg::m_minEtaFORWARD

private

Definition at line 136 of file InDetRecStatisticsAlg.h.

◆ m_minPt

float InDet::InDetRecStatisticsAlg::m_minPt

private

Minimum Pt cut for tracks used by the algorithm.

Definition at line 117 of file InDetRecStatisticsAlg.h.

◆ m_minREndPrimary

float InDet::InDetRecStatisticsAlg::m_minREndPrimary

private

If track has end vertex, this is min R of end vertex to be considered primary.

Definition at line 130 of file InDetRecStatisticsAlg.h.

◆ m_minREndSecondary

float InDet::InDetRecStatisticsAlg::m_minREndSecondary

private

If track has end vertex, this is min R of end vertex to be considered secondary.

Definition at line 131 of file InDetRecStatisticsAlg.h.

◆ m_minZEndPrimary

float InDet::InDetRecStatisticsAlg::m_minZEndPrimary

private

If track has end vertex, this is min Z of end vertex to be considered primary.

Definition at line 132 of file InDetRecStatisticsAlg.h.

◆ m_minZEndSecondary

float InDet::InDetRecStatisticsAlg::m_minZEndSecondary

private

If track has end vertex, this is min Z of end vertex to be considered secondary.

Definition at line 133 of file InDetRecStatisticsAlg.h.

◆ m_pixelID

const PixelID* InDet::InDetRecStatisticsAlg::m_pixelID

private

get pixel layer from hit ID

Definition at line 99 of file InDetRecStatisticsAlg.h.

◆ m_printSecondary

bool InDet::InDetRecStatisticsAlg::m_printSecondary

private

Flag to print hit information for secondary tracks.

Definition at line 116 of file InDetRecStatisticsAlg.h.

◆ m_pullWarning

std::atomic<bool> InDet::InDetRecStatisticsAlg::m_pullWarning {}

mutableprivate

warn only once, if pull cannot be calculated

Definition at line 141 of file InDetRecStatisticsAlg.h.

141{};

◆ m_RecTrackCollection_keys

SG::ReadHandleKeyArray<TrackCollection> InDet::InDetRecStatisticsAlg::m_RecTrackCollection_keys

private

Definition at line 138 of file InDetRecStatisticsAlg.h.

◆ m_residualPullCalculator

ToolHandle<Trk::IResidualPullCalculator> InDet::InDetRecStatisticsAlg::m_residualPullCalculator

private

The residual and pull calculator tool handle.

Definition at line 109 of file InDetRecStatisticsAlg.h.

◆ m_sctID

const SCT_ID* InDet::InDetRecStatisticsAlg::m_sctID

private

get sct layer from hit ID

Definition at line 100 of file InDetRecStatisticsAlg.h.

◆ m_SignalCounters

std::vector<class TrackStatHelper *> InDet::InDetRecStatisticsAlg::m_SignalCounters

private

Vector of TrackStatHelper objects, one for each track collection.

Definition at line 97 of file InDetRecStatisticsAlg.h.

◆ m_trackSelectorTool

ToolHandle<Trk::ITrackSelectorTool> InDet::InDetRecStatisticsAlg::m_trackSelectorTool

private

Definition at line 113 of file InDetRecStatisticsAlg.h.

◆ m_TrackTruthCollection_keys

SG::ReadHandleKeyArray<TrackTruthCollection> InDet::InDetRecStatisticsAlg::m_TrackTruthCollection_keys

private

Definition at line 139 of file InDetRecStatisticsAlg.h.

◆ m_trkSummaryTool

ToolHandle<Trk::IExtendedTrackSummaryTool> InDet::InDetRecStatisticsAlg::m_trkSummaryTool

private

tool to get track summary information from track

Definition at line 105 of file InDetRecStatisticsAlg.h.

◆ m_trtID

const TRT_ID* InDet::InDetRecStatisticsAlg::m_trtID

private

get trt layer from hit ID

Definition at line 101 of file InDetRecStatisticsAlg.h.

◆ m_truthToTrack

ToolHandle<Trk::ITruthToTrack> InDet::InDetRecStatisticsAlg::m_truthToTrack

private

tool to create track parameters from a gen particle

Definition at line 104 of file InDetRecStatisticsAlg.h.

◆ m_updator

Trk::IUpdator* InDet::InDetRecStatisticsAlg::m_updator

private

updator for unbiased states

Definition at line 108 of file InDetRecStatisticsAlg.h.

◆ m_updatorHandle

ToolHandle<Trk::IUpdator> InDet::InDetRecStatisticsAlg::m_updatorHandle

private

Tool handle of updator for unbiased states.

Definition at line 107 of file InDetRecStatisticsAlg.h.

◆ m_UpdatorWarning

std::atomic<bool> InDet::InDetRecStatisticsAlg::m_UpdatorWarning {}

mutableprivate

warn only once, if unbiased track states can not be calculated

Definition at line 142 of file InDetRecStatisticsAlg.h.

142{};

◆ m_useTrackSelection

bool InDet::InDetRecStatisticsAlg::m_useTrackSelection

private

Use track selector tool.

Definition at line 134 of file InDetRecStatisticsAlg.h.

◆ m_UseTrackSummary

bool InDet::InDetRecStatisticsAlg::m_UseTrackSummary

private

Flag to print detailed statistics for each track collection.

Definition at line 115 of file InDetRecStatisticsAlg.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Classes

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ CounterLocal

◆ StoreGateSvc_t

Member Enumeration Documentation

◆ ECounter

Constructor & Destructor Documentation

◆ InDetRecStatisticsAlg()

Member Function Documentation

◆ calculatePull()

◆ cardinality()

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ filterPassed()

◆ finalize()

◆ getIdentifier()

◆ getServices()

◆ getUnbiasedTrackParameters()

◆ initialize()

◆ inputHandles()

◆ isClonable()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ printStatistics()

◆ printTrackSummary()

◆ renounce()

◆ renounceArray()

◆ resetStatistics()

◆ selectGenSignal()

◆ selectRecSignal()

◆ setFilterPassed()

◆ sysExecute()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_detStore

◆ m_doTruth

◆ m_events_processed

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_fakeTrackCut

◆ m_fakeTrackCut2

◆ m_idDictMgr

◆ m_idHelper

◆ m_isUnbiased

◆ m_matchTrackCut

◆ m_maxEta

◆ m_maxEtaBarrel

◆ m_maxEtaEndcap

◆ m_maxEtaFORWARD

◆ m_maxEtaTransition

◆ m_maxRStartPrimary

◆ m_maxRStartSecondary

◆ m_maxZStartPrimary

◆ m_maxZStartSecondary

◆ m_McTrackCollection_key

◆ m_minEtaFORWARD

◆ m_minPt

◆ m_minREndPrimary

◆ m_minREndSecondary

◆ m_minZEndPrimary

◆ m_minZEndSecondary

◆ m_pixelID

◆ m_printSecondary

◆ m_pullWarning