InDet::TrackStatHelper Class Reference

#include <TrackStatHelper.h>

Collaboration diagram for InDet::TrackStatHelper:

Classes
struct	Counter
struct	Counter4D

Public Member Functions
	TrackStatHelper (const std::string &, const std::string &, bool careAboutTruth=true)
	Constructor. More...
void	SetCuts (const struct cuts &)
	Sets the cuts such as the eta regions (barrel, transition,endcap) and the hit fraction fake cuts and the track matching cut. More...
void	addEvent (const TrackCollection *, std::vector< const Trk::Track * > &, const std::vector< std::pair< HepMC::ConstGenParticlePtr, int > > &, const TrackTruthCollection *, const AtlasDetectorID *const, const PixelID *, const SCT_ID *, const Trk::IExtendedTrackSummaryTool *, bool, const unsigned int *, const unsigned int *) const
	Adds hit, track and matching information for each event. More...
void	reset ()
	Resets the track collection information, called in the constructor. More...
void	print (MsgStream &out) const
	Prints all of the statistics information, calls printRegion, printTrackSummaryRegion, etc for detailed statistics. More...
void	printRegion1 (MsgStream &out, enum eta_region) const
	Prints ntracks per event,efficiencies,fake rates, and general hit information for given eta region. More...
void	printRegion2 (MsgStream &out, enum eta_region, float denominator) const
void	printSecondary (MsgStream &out) const
	Prints all of the statistics information, calls printRegion, printTrackSummaryRegion, etc for detailed statistics. More...
void	printRegionSecondary (MsgStream &out, enum eta_region, float denominator) const
	Prints ntracks per event,efficiencies,fake rates, and general hit information for given eta region. More...
bool	printTrackSummaryRegion (MsgStream &out, enum track_types, enum eta_region) const
	Sets up detailed statistics part of table, calls printTrackSummaryRegion. More...
void	printTrackSummaryAverage (MsgStream &out, enum track_types, enum eta_region, int summary_type) const
	Prints information from TrackSummaryTool. More...
bool	PassTrackCuts (const Trk::TrackParameters *para) const
	defines 'good' reco tracks More...
int	ClassifyParticle (const HepMC::ConstGenParticlePtr &particle, const double prob) const
	classifies gen particle as primary, secondary or truncated More...
const std::string &	key () const
	Returns TrackCollection Key. More...
const std::string &	Truthkey () const
	Returns Truth TrackCollection Key. More...

Static Public Member Functions
static std::string	getSummaryTypeHeader ()

Private Types
enum	ETrackCounter { kTracks_rec, kTracks_gen, kTracks_gen_signal, kNTracksCounter }
enum	EHitsCounter { kHits_rec, kHits_pri, kHits_sec, kNHitsCounter }
enum	ETrackSummaryCounter { kTrackSummarySum, kNTrackSummaryOK, kNTrackSummaryBAD, kNTrackSummaryCounter }
enum	ETrackSummaryTypes {   kNumberOfInnermostPixelLayerHits, kNumberOfPixelHits, kNumberOfPixelHoles, kNumberOfGangedPixels,   kNumberOfSCTHits, kNumberOfSCTHoles, kNumberOfSCTDoubleHoles, kNumberOfTRTHits,   kNumberOfTRTOutliers, kNumberOfTRTHighThresholdHits, kNumberOfTRTHighThresholdOutliers, kNumberOfOutliersOnTrack,   kNSummaryTypes }
using	TracksCounter = Counter< kNTracksCounter, N_TRACKTYPES, N_ETAREGIONS, int >
using	TracksCounterAtomic = Counter< kNTracksCounter, N_TRACKTYPES, N_ETAREGIONS, std::atomic< long > >
using	HitsCounter = Counter< kNHitsCounter, N_HITTYPES, N_ETAREGIONS, int >
using	HitsCounterAtomic = Counter< kNHitsCounter, N_HITTYPES, N_ETAREGIONS, std::atomic< long > >
using	TrackSummaryCounter = Counter4D< kNTrackSummaryCounter, N_TRACKTYPES, N_ETAREGIONS, kNSummaryTypes, int >
	type statistics. More...
using	TrackSummaryCounterAtomic = Counter4D< kNTrackSummaryCounter, N_TRACKTYPES, N_ETAREGIONS, kNSummaryTypes, std::atomic< long > >
typedef std::multimap< HepMcParticleLink, float >	recoToTruthMap
	map containing reco track and matched truth track barcode More...

Private Member Functions
void	setSummaryStat (track_types track_i, eta_region region_i, const Trk::TrackSummary *summary, TrackSummaryCounter &trackSummarySum) const

Private Attributes
std::string	m_TrackCollectionKey
	StoreGate Track Collection Key. More...
std::string	m_TrackTruthCollectionKey
	StoreGate Track Truth Collection Key. More...
std::atomic< bool >	m_author_found [Trk::TrackInfo::NumberOfTrackFitters] {}
	Number of tracking authors found. More...
std::atomic< long >	m_events {}
	Number of events. More...
TracksCounterAtomic m_tracks	ATLAS_THREAD_SAFE
HitsCounterAtomic m_hits	ATLAS_THREAD_SAFE
TrackSummaryCounterAtomic m_trackSummarySum	ATLAS_THREAD_SAFE
std::atomic< bool >	m_truthMissing
	Flag for if track truth is missing. More...
bool	m_careAboutTruth
struct cuts	m_cuts
std::mutex	m_authorMutex
std::bitset< Trk::TrackInfo::NumberOfTrackRecoInfo > m_recoInfo	ATLAS_THREAD_SAFE
std::bitset< Trk::TrackInfo::NumberOfTrackProperties > m_patternProperties	ATLAS_THREAD_SAFE

Static Private Attributes
static const Trk::SummaryType	s_summaryTypes [kNSummaryTypes]
	summary types for which statistics are gathered More...
static const char *const	s_summaryTypeName [kNSummaryTypes]
	table column labels for summary More...

Detailed Description

Definition at line 126 of file TrackStatHelper.h.

Member Typedef Documentation

HitsCounter

using InDet::TrackStatHelper::HitsCounter = Counter<kNHitsCounter,N_HITTYPES, N_ETAREGIONS, int>

private

Definition at line 260 of file TrackStatHelper.h.

HitsCounterAtomic

using InDet::TrackStatHelper::HitsCounterAtomic = Counter<kNHitsCounter,N_HITTYPES, N_ETAREGIONS, std::atomic<long> >

private

Definition at line 261 of file TrackStatHelper.h.

recoToTruthMap

typedef std::multimap<HepMcParticleLink,float> InDet::TrackStatHelper::recoToTruthMap

private

map containing reco track and matched truth track barcode

Definition at line 335 of file TrackStatHelper.h.

TracksCounter

using InDet::TrackStatHelper::TracksCounter = Counter<kNTracksCounter,N_TRACKTYPES,N_ETAREGIONS, int>

private

Definition at line 250 of file TrackStatHelper.h.

TracksCounterAtomic

using InDet::TrackStatHelper::TracksCounterAtomic = Counter<kNTracksCounter,N_TRACKTYPES,N_ETAREGIONS, std::atomic<long> >

private

Definition at line 251 of file TrackStatHelper.h.

TrackSummaryCounter

using InDet::TrackStatHelper::TrackSummaryCounter = Counter4D<kNTrackSummaryCounter, N_TRACKTYPES, N_ETAREGIONS, kNSummaryTypes, int>

private

type statistics.

Definition at line 294 of file TrackStatHelper.h.

TrackSummaryCounterAtomic

using InDet::TrackStatHelper::TrackSummaryCounterAtomic = Counter4D<kNTrackSummaryCounter, N_TRACKTYPES, N_ETAREGIONS, kNSummaryTypes, std::atomic<long> >

private

Definition at line 299 of file TrackStatHelper.h.

Member Enumeration Documentation

EHitsCounter

enum InDet::TrackStatHelper::EHitsCounter

private

Enumerator
kHits_rec	number of reconstructed hits for a given type and eta region
kHits_pri	number of hits from primary tracks for a given type and eta region
kHits_sec	number of hits from secondary tracks for a given type and eta region
kNHitsCounter

Definition at line 254 of file TrackStatHelper.h.

                      {
       kHits_rec,     
       kHits_pri,     
       kHits_sec,     
       kNHitsCounter
    };

ETrackCounter

enum InDet::TrackStatHelper::ETrackCounter

private

Enumerator
kTracks_rec	number of reconstructed tracks for a given type and eta region
kTracks_gen	number of generated tracks for a given type and eta region, looping over genevents to include possible pileup
kTracks_gen_signal	number of generated tracks for a given type and eta region, just from first genevent
kNTracksCounter

Definition at line 245 of file TrackStatHelper.h.

                       {kTracks_rec,        
                        kTracks_gen,        
                        kTracks_gen_signal, 
                        kNTracksCounter};

ETrackSummaryCounter

enum InDet::TrackStatHelper::ETrackSummaryCounter

private

Enumerator
kTrackSummarySum	Track Summary Values for each track type, region and summary type.
kNTrackSummaryOK	Number of tracks with track summary OK for given type,eta,summary type.
kNTrackSummaryBAD	Number of tracks with track summary bad for given type,eta,summary type.
kNTrackSummaryCounter

Definition at line 263 of file TrackStatHelper.h.

                              {
       kTrackSummarySum,     
       kNTrackSummaryOK,     
       kNTrackSummaryBAD,    
       kNTrackSummaryCounter
    };

ETrackSummaryTypes

enum InDet::TrackStatHelper::ETrackSummaryTypes

private

Enumerator
kNumberOfInnermostPixelLayerHits
kNumberOfPixelHits
kNumberOfPixelHoles
kNumberOfGangedPixels
kNumberOfSCTHits
kNumberOfSCTHoles
kNumberOfSCTDoubleHoles
kNumberOfTRTHits
kNumberOfTRTOutliers
kNumberOfTRTHighThresholdHits
kNumberOfTRTHighThresholdOutliers
kNumberOfOutliersOnTrack
kNSummaryTypes

Definition at line 274 of file TrackStatHelper.h.

                            {
        kNumberOfInnermostPixelLayerHits,
        kNumberOfPixelHits,
        kNumberOfPixelHoles,
        kNumberOfGangedPixels,
        kNumberOfSCTHits,
        kNumberOfSCTHoles,
        kNumberOfSCTDoubleHoles,
        kNumberOfTRTHits,
        kNumberOfTRTOutliers,
        kNumberOfTRTHighThresholdHits,
        kNumberOfTRTHighThresholdOutliers,
        kNumberOfOutliersOnTrack,
        kNSummaryTypes
    };

Constructor & Destructor Documentation

TrackStatHelper()

InDet::TrackStatHelper::TrackStatHelper	(	const std::string &	TrackCollectionKey,
		const std::string &	TrackTruthCollectionKey,
		bool	careAboutTruth = true
	)

Constructor.

Takes strings of TrackCollectionKey and Track TruthCollectionKey. Assume that truth collections will be used, but let user specify if they intentionally will not be.

Definition at line 103 of file TrackStatHelper.cxx.

                                                                                                                                       :
  m_TrackCollectionKey       (std::move(TrackCollectionKey)),
  m_TrackTruthCollectionKey  (std::move(TrackTruthCollectionKey)),
  m_truthMissing             (false),
  m_careAboutTruth           (careAboutTruth)
{
  this->reset();
}

Member Function Documentation

addEvent()

void InDet::TrackStatHelper::addEvent	(	const TrackCollection *	recTracks,
		std::vector< const Trk::Track * > &	rec,
		const std::vector< std::pair< HepMC::ConstGenParticlePtr, int > > &	gen,
		const TrackTruthCollection *	truthMap,
		const AtlasDetectorID * const	idHelper,
		const PixelID *	pixelID,
		const SCT_ID *	sctID,
		const Trk::IExtendedTrackSummaryTool *	trkSummaryTool,
		bool	useTrackSummary,
		const unsigned int *	inTimeStart,
		const unsigned int *	inTimeEnd
	)		const

Adds hit, track and matching information for each event.

Called at each event

Definition at line 119 of file TrackStatHelper.cxx.

{
 
  recoToTruthMap rttMap;
  
  recoToTruthMap::const_iterator imap;
  
  m_events ++;
  
  float Eta = 0;
  InDet::eta_region Region = ETA_ALL;
  int Author = 0;
  int recoClassification = 0;
  bool truth_missing=false;
 
  TracksCounter       tracks;
  HitsCounter         hits;
  TrackSummaryCounter trackSummarySum;
  bool                author_found [Trk::TrackInfo::NumberOfTrackFitters];
  std::bitset<Trk::TrackInfo::NumberOfTrackRecoInfo>   reco_info;
  std::bitset<Trk::TrackInfo::NumberOfTrackProperties> pattern_properties;
  for (unsigned int i=0; i < Trk::TrackInfo::NumberOfTrackFitters; ++i) { author_found[i]=false; }
 
  // ------------ reconstructed tracks -----------------
  for (const Trk::Track* track : rec) {
 
    const Trk::TrackParameters* para = track->trackParameters()->front();    
    recoClassification = -999;
    int thisEventIndex = -999;
    
    Author = track->info().trackFitter();
    if (Author > 0 && Author < Trk::TrackInfo::NumberOfTrackFitters){
      author_found[Author]  = true;
      reco_info |= track->info().patternRecognition();
      pattern_properties |= track->info().properties();
    }
    else {
      //FIX author_problem++;
    }
    if (!para) {
      Region =  ETA_UNKNOWN;
    }
    else {
      Eta = std::abs(para->eta());
      if (Eta < m_cuts.maxEtaBarrel) Region = ETA_BARREL;
      else if  (Eta < m_cuts.maxEtaTransition) Region = ETA_TRANSITION;
      else if  (Eta < m_cuts.maxEtaEndcap) Region = ETA_ENDCAP;
      else if ((Eta > m_cuts.minEtaFORWARD) && (Eta < m_cuts.maxEtaFORWARD)) Region = ETA_FORWARD;
      else Region = ETA_OUTSIDE;
    }
    
    if ( !(Region==ETA_BARREL || Region==ETA_TRANSITION || Region==ETA_ENDCAP || Region==ETA_FORWARD) )
      {
    continue; // Only want to tally tracks that are within the detector volume.
      }
 
    TrackTruthCollection::const_iterator found;
    TrackTruth trtruth;
    double     trprob = 0;  
 
    if (!truthMap) {
      // no truthmap
      truth_missing=true;
    }
    else {
      ElementLink<TrackCollection> tracklink;
      tracklink.setElement(track);
      tracklink.setStorableObject(*recTracks);
      const ElementLink<TrackCollection> tracklink2=tracklink;
      
      found = truthMap->find(tracklink2);
      if(!(found == truthMap->end())){
    trtruth=found->second;
    HepMcParticleLink hmpl = trtruth.particleLink();
    thisEventIndex = hmpl.eventIndex();
    rttMap.insert(std::pair<HepMcParticleLink,float>(hmpl,trtruth.probability()));
    // ME: remember prob
    trprob = trtruth.probability();
      }
      
    }
    
    
    if(!PassTrackCuts(para)) continue;
    
    tracks.m_counter[kTracks_rec][TRACK_ALL][Region]++;
    tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_ALL]++;
    // signal only tracks for the denominator in signal efficiencies
    if (thisEventIndex==0)
      {
    tracks.m_counter[kTracks_rec][TRACK_ALL_SIGNAL][Region]++;
    tracks.m_counter[kTracks_rec][TRACK_ALL_SIGNAL][ETA_ALL]++;
      }
    
    // process track summary
    
    std::unique_ptr<Trk::TrackSummary> cleanup;
    const Trk::TrackSummary* summary = track->trackSummary();
    
    if (useTrackSummary) {
       if (!track->trackSummary()) {
          cleanup = trkSummaryTool->summary(*track);
          summary=cleanup.get();
       }
 
       setSummaryStat(TRACK_ALL, Region, summary, trackSummarySum);
    }
      
    // check if current reco track has more than a specified fraction of hits
    // from a generated track. [If not, it will be counted as a fake
    // track in the statistics report.] 
      
    if (!truthMap) {
      truth_missing=true;
    }
    else {// no truthmap
    
      if (found == truthMap->end()) {
    // m_truthMissing=true;
    // no truth might happen with new truth !
    tracks.m_counter[kTracks_rec][TRACK_NOHEPMCPARTICLELINK][Region]++;
    tracks.m_counter[kTracks_rec][TRACK_NOHEPMCPARTICLELINK][ETA_ALL]++;
        setSummaryStat(TRACK_NOHEPMCPARTICLELINK, Region, summary, trackSummarySum);
      }
      else{   // no link  
      
    // ME : change logic, secondaries from G4 processes are truncated
    const HepMcParticleLink& HMPL=trtruth.particleLink();
    if (! HMPL.isValid()) {
      tracks.m_counter[kTracks_rec][TRACK_NOHEPMCPARTICLELINK][Region]++;
      tracks.m_counter[kTracks_rec][TRACK_NOHEPMCPARTICLELINK][ETA_ALL]++;
          setSummaryStat(TRACK_NOHEPMCPARTICLELINK, Region, summary, trackSummarySum);
    }
    else {
      //classify track as coming from primary, secondary or truncated gen particle
#ifdef HEPMC3
      HepMC::ConstGenParticlePtr particle = HMPL.scptr();
#else
      const HepMC::GenParticle *particle = HMPL.cptr();
#endif
      recoClassification = ClassifyParticle(particle, trprob);
        
      if (trprob < m_cuts.fakeTrackCut)
        {
          tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB][Region]++;
          tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB][ETA_ALL]++;
          if (thisEventIndex==0)
        {
          tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB_SIGNAL][Region]++;
          tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB_SIGNAL][ETA_ALL]++;
        }
              setSummaryStat(TRACK_LOWTRUTHPROB, Region, summary, trackSummarySum);
        }
      if (trprob < m_cuts.fakeTrackCut2) {
        tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB2][Region]++;
        tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB2][ETA_ALL]++;
        if (thisEventIndex==0)
          {
        tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB2_SIGNAL][Region]++;
        tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB2_SIGNAL][ETA_ALL]++;
          }
            setSummaryStat(TRACK_LOWTRUTHPROB2, Region, summary, trackSummarySum);
      }
    } // end else on !isVadid
      }
    }
      
    // ------------------ hits on reconstructed tracks ---------------
      
    EHitsCounter part_type = ( recoClassification==TRACK_PRIMARY
                               ? kHits_pri
                               : ( recoClassification==TRACK_SECONDARY
                                   ? kHits_sec
                                   : kNHitsCounter));
    for (const Trk::TrackStateOnSurface* hit : *track->trackStateOnSurfaces()) {
      if(hit){ 
    const Trk::MeasurementBase* mesh =hit->measurementOnTrack();
    if (mesh) {
      const Trk::RIO_OnTrack* rio{} ;
      if ( mesh->type(Trk::MeasurementBaseType::RIO_OnTrack)){
            assert(dynamic_cast<const Trk::RIO_OnTrack*>(mesh)!=nullptr);
        rio = static_cast<const Trk::RIO_OnTrack*>(mesh);
      } else  {
        // try CompetingROT:
        if (mesh->type(Trk::MeasurementBaseType::CompetingRIOsOnTrack)) {
              assert(dynamic_cast<const Trk::CompetingRIOsOnTrack*>(mesh) !=nullptr);
          auto comprot = static_cast<const Trk::CompetingRIOsOnTrack*>(mesh);
          rio = &comprot->rioOnTrack(comprot->indexOfMaxAssignProb());
        }
      }
        
        
      if (rio) {
        // skip outliers
        if (!hit->type(Trk::TrackStateOnSurface::Measurement)) 
          continue;
          
        hits.m_counter[kHits_rec][HIT_ALL][Region]++;
        hits.m_counter[kHits_rec][HIT_ALL][ETA_ALL]++;
        if(part_type<kNHitsCounter){
          hits.m_counter[part_type][HIT_ALL][Region]++;
          hits.m_counter[part_type][HIT_ALL][ETA_ALL]++;
        }
          
            Identifier id = rio->identify();
            int HitDet    = HIT_UNKNOWN;
            int HitLayer  = N_HITTYPES;
            bool part_type_for_all=true;
            if (idHelper->is_trt(id)){
               HitDet=HIT_TRT_ALL;
            }
            else if (idHelper->is_sct(id)){
              HitDet=HIT_SCT_ALL;
              if (sctID) {
                switch (sctID->layer_disk(id)) {
                case 0: HitLayer = HIT_SCT1; break;
                case 1: HitLayer = HIT_SCT2; break;
                case 2: HitLayer = HIT_SCT3; break;
                case 3: HitLayer = HIT_SCT4; break;
                case 4: case 5: case 6: case 7: case 8: HitLayer = HIT_SCT5TO9; break;
                default: HitLayer = HIT_UNKNOWN; break;
                }
              }
            }
            else if (idHelper->is_pixel(id)){
              part_type_for_all=false;
              HitDet = HIT_PIXEL_ALL;
              if (pixelID) {
                switch (pixelID->layer_disk(id)) {
                case 0:  HitLayer = HIT_PIX1; break;
                case 1:  HitLayer = HIT_PIX2; break;
                case 2:  HitLayer = HIT_PIX3; break;
                default: HitLayer = HIT_UNKNOWN;
                }
              }
            }
 
            hits.m_counter[kHits_rec][HitDet][Region]++;
            hits.m_counter[kHits_rec][HitDet][ETA_ALL]++;
            if(part_type_for_all && part_type<kNHitsCounter){
               hits.m_counter[part_type][HitDet][Region]++;
               hits.m_counter[part_type][HitDet][ETA_ALL]++;
            }
 
            if (HitLayer<N_HITTYPES) {
               hits.m_counter[kHits_rec][HitLayer][Region]++;
               hits.m_counter[kHits_rec][HitLayer][ETA_ALL]++;
               if(part_type<kNHitsCounter){
                  hits.m_counter[part_type][HitLayer][Region]++;
                  hits.m_counter[part_type][HitLayer][ETA_ALL]++;
               }
            }
          }
    }
      }
    }
  }
  
  // ------------------------- generated tracks, including pileup -----------------------------
  
  Eta = 0;
  Region = ETA_ALL;
  int classification=-999;
  for (auto truth = gen.begin(); truth != gen.end();  ++truth) {
    classification=-999; 
    bool inTimePileup = truth->second == 0 || (truth->second >= (int)*inTimeStart && truth->second <= (int)*inTimeEnd);
    
    auto particle = truth->first;
    
    //determine eta region
    Eta = std::abs(particle->momentum().pseudoRapidity());
    if (Eta < m_cuts.maxEtaBarrel) Region = ETA_BARREL;
    else if  (Eta < m_cuts.maxEtaTransition) Region = ETA_TRANSITION;
    else if  (Eta < m_cuts.maxEtaEndcap) Region = ETA_ENDCAP;
    else if ((Eta > m_cuts.minEtaFORWARD) && (Eta < m_cuts.maxEtaFORWARD)) Region = ETA_FORWARD;
    else Region = ETA_OUTSIDE;
    if ( !(Region==ETA_BARREL || Region==ETA_TRANSITION || Region==ETA_ENDCAP || Region==ETA_FORWARD) )
      {
    continue; // Only want to tally tracks that are within the detector volume.
      }
    tracks.m_counter[kTracks_gen][TRACK_ALL][Region]++;
    tracks.m_counter[kTracks_gen][TRACK_ALL][ETA_ALL] ++;
    if (inTimePileup) tracks.m_counter[kTracks_gen][TRACK_ALL_SIGNAL][ETA_ALL]++; // "SIGNAL" is misnomer here
    if (inTimePileup) tracks.m_counter[kTracks_gen][TRACK_ALL_SIGNAL][Region]++;
 
    //classify gen as primary, secondary or truncated
    classification = ClassifyParticle(particle,1.);
    
    if(classification==TRACK_PRIMARY){
      tracks.m_counter[kTracks_gen][TRACK_PRIMARY][ETA_ALL] ++;
      tracks.m_counter[kTracks_gen][TRACK_PRIMARY][Region]++;
    }
    if(classification==TRACK_SECONDARY){
      tracks.m_counter[kTracks_gen][TRACK_SECONDARY][ETA_ALL] ++;
      tracks.m_counter[kTracks_gen][TRACK_SECONDARY][Region]++;
    }
    
    //see if gen track has at least 1 matching reco track with high enough probability
    bool matched = false;
    int nmatched = 0;
    
    HepMcParticleLink hmpl2(particle,particle->parent_event()->event_number(),HepMcParticleLink::IS_EVENTNUM);
    recoToTruthMap::iterator rttIter=rttMap.find(hmpl2);
    if(rttIter != rttMap.end()){
      for(imap = rttMap.lower_bound(hmpl2); imap !=rttMap.upper_bound(hmpl2); ++imap){
    if(imap->second > m_cuts.matchTrackCut){
      matched = true;
      nmatched++;
    }
      }
    }
    if (matched) {
      tracks.m_counter[kTracks_gen][TRACK_MATCHED][Region]++;
      tracks.m_counter[kTracks_gen][TRACK_MATCHED][ETA_ALL]++;
      if (inTimePileup)
    {
      tracks.m_counter[kTracks_gen][TRACK_MATCHED_SIGNAL][Region]++;
      tracks.m_counter[kTracks_gen][TRACK_MATCHED_SIGNAL][ETA_ALL]++;
    }
      if(nmatched > 1){
    tracks.m_counter[kTracks_gen][TRACK_MULTMATCH][Region] += nmatched-1;
    tracks.m_counter[kTracks_gen][TRACK_MULTMATCH][ETA_ALL]+= nmatched-1;
      }
      if(classification==TRACK_PRIMARY){
    tracks.m_counter[kTracks_gen][TRACK_MATCHED_PRIMARY][Region]++;
    tracks.m_counter[kTracks_gen][TRACK_MATCHED_PRIMARY][ETA_ALL]++;
    if(nmatched > 1){
      tracks.m_counter[kTracks_gen][TRACK_MULTMATCH_PRIMARY][Region]  += nmatched-1;
      tracks.m_counter[kTracks_gen][TRACK_MULTMATCH_PRIMARY][ETA_ALL] += nmatched-1;
    }
      }
      if(classification==TRACK_SECONDARY){
    tracks.m_counter[kTracks_gen][TRACK_MATCHED_SECONDARY][Region]++;
    tracks.m_counter[kTracks_gen][TRACK_MATCHED_SECONDARY][ETA_ALL]++;
    if(nmatched > 1){
      tracks.m_counter[kTracks_gen][TRACK_MULTMATCH_SECONDARY][Region]  += nmatched-1;
      tracks.m_counter[kTracks_gen][TRACK_MULTMATCH_SECONDARY][ETA_ALL] += nmatched-1;
    }
      }
    }   
  }
  
  // ------------------------- generated tracks, just signal event -----------------------------
  // For efficiencies.
  
  Eta = 0;
  Region = ETA_ALL;
  classification=-999;
  
  for (auto truth = gen.begin(); truth != gen.end();  ++truth) 
    {
      if (truth->second != 0) // only signal event GenParticles
    continue;
      
      classification=-999; 
      
      auto particle = truth->first;
      
      //determine eta region
      Eta = std::abs(particle->momentum().pseudoRapidity());
      if (Eta < m_cuts.maxEtaBarrel) Region = ETA_BARREL;
      else if  (Eta < m_cuts.maxEtaTransition) Region = ETA_TRANSITION;
      else if  (Eta < m_cuts.maxEtaEndcap) Region = ETA_ENDCAP;
      else if ((Eta > m_cuts.minEtaFORWARD) && (Eta < m_cuts.maxEtaFORWARD)) Region = ETA_FORWARD;
      else Region = ETA_OUTSIDE;
      if ( !(Region==ETA_BARREL || Region==ETA_TRANSITION || Region==ETA_ENDCAP || Region==ETA_FORWARD) )
    {
      continue; // Only want to tally tracks that are within the detector volume.
    }
      tracks.m_counter[kTracks_gen_signal][TRACK_ALL][ETA_ALL] ++;
      tracks.m_counter[kTracks_gen_signal][TRACK_ALL][Region]++;
      
      //classify gen partilce as primary, secondary or truncated
      classification = ClassifyParticle(particle,1.);
      
      if(classification==TRACK_PRIMARY){
    tracks.m_counter[kTracks_gen_signal][TRACK_PRIMARY][ETA_ALL] ++;
    tracks.m_counter[kTracks_gen_signal][TRACK_PRIMARY][Region]++;
      }
      if(classification==TRACK_SECONDARY){
    tracks.m_counter[kTracks_gen_signal][TRACK_SECONDARY][ETA_ALL] ++;
    tracks.m_counter[kTracks_gen_signal][TRACK_SECONDARY][Region]++;
      }
      
      //see if gen track has at least 1 matching reco track with high enough probability
      bool matched = false;
      int nmatched = 0;
      
      HepMcParticleLink hmpl2(particle,truth->second,HepMcParticleLink::IS_EVENTNUM); // FIXME truth->second is actually the position of the GenEvent in the McEventCollection!! See InDetRecStatisticsAlg::selectGenSignal(...) method (only client of TrackStatsHelper)
      recoToTruthMap::iterator rttIter=rttMap.find(hmpl2);
      if(rttIter != rttMap.end()){
    for(imap = rttMap.lower_bound(hmpl2); imap !=rttMap.upper_bound(hmpl2); ++imap){
      if(imap->second > m_cuts.matchTrackCut){
        matched = true;
        nmatched++;
      }
    }
      }
      if (matched) {
    tracks.m_counter[kTracks_gen_signal][TRACK_MATCHED][Region]++;
    tracks.m_counter[kTracks_gen_signal][TRACK_MATCHED][ETA_ALL]++;
    if(nmatched > 1){
      tracks.m_counter[kTracks_gen_signal][TRACK_MULTMATCH][Region] += nmatched-1;
      tracks.m_counter[kTracks_gen_signal][TRACK_MULTMATCH][ETA_ALL]+= nmatched-1;
    }
    if(classification==TRACK_PRIMARY){
      tracks.m_counter[kTracks_gen_signal][TRACK_MATCHED_PRIMARY][Region]++;
      tracks.m_counter[kTracks_gen_signal][TRACK_MATCHED_PRIMARY][ETA_ALL]++;
      if(nmatched > 1){
        tracks.m_counter[kTracks_gen_signal][TRACK_MULTMATCH_PRIMARY][Region]  += nmatched-1;
        tracks.m_counter[kTracks_gen_signal][TRACK_MULTMATCH_PRIMARY][ETA_ALL] += nmatched-1;
      }
    }
    if(classification==TRACK_SECONDARY){
      tracks.m_counter[kTracks_gen_signal][TRACK_MATCHED_SECONDARY][Region]++;
      tracks.m_counter[kTracks_gen_signal][TRACK_MATCHED_SECONDARY][ETA_ALL]++;
      if(nmatched > 1){
        tracks.m_counter[kTracks_gen_signal][TRACK_MULTMATCH_SECONDARY][Region]  += nmatched-1;
        tracks.m_counter[kTracks_gen_signal][TRACK_MULTMATCH_SECONDARY][ETA_ALL] += nmatched-1;
      }
    }
      } 
    }  
 
  if (truth_missing) {
     m_truthMissing = truth_missing;
  }
  m_tracks          += tracks;
  m_hits            += hits;
  m_trackSummarySum += trackSummarySum;
  for (unsigned int i=0; i < Trk::TrackInfo::NumberOfTrackFitters; ++i) {
     if (author_found[i]) {
        m_author_found[i] = author_found[i];
     }
  }
  {
     std::lock_guard<std::mutex> lock(m_authorMutex);
     m_recoInfo          |= reco_info;
     m_patternProperties |= pattern_properties;
  }
}

ClassifyParticle()

int InDet::TrackStatHelper::ClassifyParticle	(	const HepMC::ConstGenParticlePtr &	particle,
		const double	prob
	)		const

classifies gen particle as primary, secondary or truncated

Definition at line 850 of file TrackStatHelper.cxx.

                                                                                                             {
 
  int partClass=-999;
 
  // ME: only classify matched tracks
  if (prob <= m_cuts.matchTrackCut) return partClass;  
 
  //classify as primary, secondary or truncated
  bool  primary=false;
  bool  secondary=false;
  bool  truncated=false;
 
  if (particle->production_vertex()){
 
    // primary vertex inside innermost layer?
    HepGeom::Point3D<double>    startVertex(particle->production_vertex()->position().x(),
                particle->production_vertex()->position().y(),
                particle->production_vertex()->position().z());
    if ( std::abs(startVertex.perp()) < m_cuts.maxRStartPrimary 
     && std::abs(startVertex.z()) < m_cuts.maxZStartPrimary) { 
      if (particle->end_vertex() == nullptr) {  
    primary=true;
      }
      else {
    HepGeom::Point3D<double> endVertex(particle->end_vertex()->position().x(), 
                 particle->end_vertex()->position().y(),
                 particle->end_vertex()->position().z());
    if (  endVertex.perp()         > m_cuts.minREndPrimary 
          || std::abs(startVertex.z()) > m_cuts.minZEndPrimary){
      primary=true;
    }
    else { 
      truncated = true;
    }
      }
    }
    
    else if ( startVertex.perp() <  m_cuts.maxRStartSecondary && std::abs(startVertex.z()) <  m_cuts.maxZStartSecondary) {
      if (particle->end_vertex() == nullptr) {  
    secondary=true;
      }
      else {
    HepGeom::Point3D<double> endVertex(particle->end_vertex()->position().x(),
                 particle->end_vertex()->position().y(),
                 particle->end_vertex()->position().z());
 
    if (endVertex.perp()        > m_cuts.minREndSecondary
        || std::abs(endVertex.z())  > m_cuts.minZEndSecondary) {
      secondary=true;
    }
      }
    }
  }//end classification
 
  if(secondary)partClass=TRACK_SECONDARY;
  if(primary)partClass=TRACK_PRIMARY;
  if(truncated)partClass=TRACK_TRUNCATED;
 
  return partClass;
}

getSummaryTypeHeader()

std::string InDet::TrackStatHelper::getSummaryTypeHeader ( )

static

Definition at line 81 of file TrackStatHelper.cxx.

                                                     {
   std::stringstream out;
   for (unsigned int stype_i=0; stype_i < kNSummaryTypes; ++stype_i ) {
      out << std::setw(std::max(6,static_cast<int>(strlen(s_summaryTypeName[stype_i])+1))) << s_summaryTypeName[stype_i];
   }
   return out.str();
}

key()

const std::string& InDet::TrackStatHelper::key ( ) const

inline

Returns TrackCollection Key.

Definition at line 177 of file TrackStatHelper.h.

{ return m_TrackCollectionKey; }

PassTrackCuts()

bool InDet::TrackStatHelper::PassTrackCuts

(

const Trk::TrackParameters *

para

)

const

defines 'good' reco tracks

Definition at line 841 of file TrackStatHelper.cxx.

                                                                             {
  bool passed = false;
  if(para->pT() >  m_cuts.minPt && std::abs(para->eta()) < m_cuts.maxEtaFORWARD)passed = true;
 
 
  return passed;
 
}

print()

void InDet::TrackStatHelper::print

(

MsgStream &

out

)

const

Prints all of the statistics information, calls printRegion, printTrackSummaryRegion, etc for detailed statistics.

Definition at line 584 of file TrackStatHelper.cxx.

                                                     {
 
  out << " Printing Statistics for " << key();
 
  out << "TrackCollection \"" << key() << "\" " << "\n" ;
 
  if (m_events > 0) {
    out << "(TrackAuthors:";
    std::vector<std::string> author_string;
    author_string.reserve(Trk::TrackInfo::NumberOfTrackFitters);
    for (int i=0; i<Trk::TrackInfo::NumberOfTrackFitters; ++i) {
       author_string.push_back( Trk::TrackInfo(static_cast<Trk::TrackInfo::TrackFitter>(i),
                                               Trk::undefined,
                                               m_patternProperties,
                                               m_recoInfo).dumpInfo() );
    }
    for (int i=0; i<Trk::TrackInfo::NumberOfTrackFitters; i++){
      if (m_author_found[i]){
    out << " " << author_string[i];
      }
    }
    out << " )" << std::endl
          << "TrackTruthCollection \"" << Truthkey() << "\"" << std::endl;
 
    if (m_truthMissing)
      {
    if (m_careAboutTruth)
      out << " WARNING: TrackTruth missing for part of this TrackCollection, no efficiencies or fakerates included.!" << std::endl;
    else
      out << " INFO: Intentionally no TrackTruth for this TrackCollection, no efficiencies or fakerates included." << std::endl;
      }
    if (!m_truthMissing)
      {
    out << " \t\t\t  ................................tracks................................" << std::endl;
    out << " \t\t\tn/event\tSEff.pr\tSEff.sd\tSLowP1\tSLowP2\tLowP1\tLowP2\tnoLink\tmultM\t" << std::endl;
    out << " total     "; printRegion1(out,ETA_ALL);
    out << " in barrel "; printRegion1(out,ETA_BARREL);
    out << " in trans. "; printRegion1(out,ETA_TRANSITION);
    out << " in endcap "; printRegion1(out,ETA_ENDCAP);
    out << " in forwa."; printRegion1(out,ETA_FORWARD);
 
          }
    else // no truth
      {
    out << "\t" << "tracks, n/event:" << std::endl;
    out << "\t\t" << "total" << std::setiosflags(std::ios::fixed | std::ios::showpoint) <<
      std::setw(7) << std::setprecision(2) << "\t" << m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_ALL]/(float) m_events << std::endl << std::setprecision(-1);
    out << "\t\t" << "in barrel" << std::setiosflags(std::ios::fixed | std::ios::showpoint) <<
      std::setw(7) << std::setprecision(2) << "\t" << m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_BARREL]/(float) m_events << std::endl << std::setprecision(-1);
    out << "\t\t" << "in trans." << std::setiosflags(std::ios::fixed | std::ios::showpoint) <<
      std::setw(7) << std::setprecision(2) << "\t" << m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_TRANSITION]/(float) m_events << std::endl << std::setprecision(-1);
    out << "\t\t" << "in endcap" << std::setiosflags(std::ios::fixed | std::ios::showpoint) <<
      std::setw(7) << std::setprecision(2) << "\t" << m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_ENDCAP]/(float) m_events << std::endl << std::setprecision(-1);
    out << "\t\t" << "in forwa." << std::setiosflags(std::ios::fixed | std::ios::showpoint) <<
          std::setw(7) << std::setprecision(2) << "\t" << m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_FORWARD]/(float) m_events << std::endl << std::setprecision(-1);
      }
    out << " \t\t\t ....................................hits/track............................" << std::endl;
    out << " \t\t\t\ttotal\tPIX1\tPIX2\tPIX3\tSCT1\tSCT2\tSCT3\tSCT4\tSCT5to9\tStraws" << std::endl;
    out << " total     "; printRegion2(out,ETA_ALL, (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_ALL]);
    out << " in barrel "; printRegion2(out,ETA_BARREL, (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_BARREL]);
    out << " in trans. "; printRegion2(out,ETA_TRANSITION, (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_TRANSITION] );
    out << " in endcap "; printRegion2(out,ETA_ENDCAP, (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_ENDCAP] );
    out << " in forwa. "; printRegion2(out,ETA_FORWARD, (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL][ETA_FORWARD] );
 
  }
  else
    out << ": NO EVENTS PROCESSED! " << std::endl;
}

printRegion1()

void InDet::TrackStatHelper::printRegion1	(	MsgStream &	out,
		enum eta_region	region
	)		const

Prints ntracks per event,efficiencies,fake rates, and general hit information for given eta region.

Definition at line 653 of file TrackStatHelper.cxx.

                                                                                    {
  out  <<  std::setiosflags(std::ios::fixed | std::ios::showpoint) << std::setw(7) << std::setprecision(2)
         << "\t" << m_tracks.m_counter[kTracks_rec][TRACK_ALL][region]/(float) m_events;
 
  if (m_tracks.m_counter[kTracks_gen_signal][TRACK_PRIMARY][region]) {
    out << "\t" <<  std::setprecision(4)
        << m_tracks.m_counter[kTracks_gen_signal][TRACK_MATCHED_PRIMARY][region]/ (float) m_tracks.m_counter[kTracks_gen_signal][TRACK_PRIMARY][region]     // track efficiency wrt. gen primary in signal
        <<  std::setprecision(2);
  }
  else {
    out  << "\t" << "n/a" ;
  }
  if (m_tracks.m_counter[kTracks_gen_signal][TRACK_SECONDARY][region]) {
    out << "\t" <<  std::setprecision(4)
        << m_tracks.m_counter[kTracks_gen_signal][TRACK_MATCHED_SECONDARY][region]/ (float) m_tracks.m_counter[kTracks_gen_signal][TRACK_SECONDARY][region]     // track efficiency wrt. gen secondary in signal
        <<  std::setprecision(2);
  }
  else {
    out  << "\t" << "n/a" ;
  }
  
  if (m_tracks.m_counter[kTracks_rec][TRACK_ALL_SIGNAL][region]) {
    out << "\t" << 100*(m_tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB_SIGNAL][region]/ (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL_SIGNAL][region]) << "%"  // track fake rate for signal only
        << "\t" << 100*(m_tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB2_SIGNAL][region]/ (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL_SIGNAL][region]) << "%";  // track fake rate for signal only
  }
  else {
    out  << "\t" << "n/a\tn/a" ;
  }
  
  if (m_tracks.m_counter[kTracks_rec][TRACK_ALL][region]) {
    out << "\t" << 100*(m_tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB][region]/ (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL][region]) << "%"  // track fake rate
        << "\t" << 100*(m_tracks.m_counter[kTracks_rec][TRACK_LOWTRUTHPROB2][region]/ (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL][region]) << "%"  // track fake rate
        << "\t" << 100*(m_tracks.m_counter[kTracks_rec][TRACK_NOHEPMCPARTICLELINK][region]/ (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL][region])  << "%" // rate of tracks without HepMcParticleLink
        << "\t" << 100*(m_tracks.m_counter[kTracks_gen][TRACK_MULTMATCH][region]/ (float) m_tracks.m_counter[kTracks_rec][TRACK_ALL][region])  << "%"; // rate of tracks matched to the same truth track
 
  }
  else{
    out  << "\tn/a\tn/a\tn/a";
  }
  out << std::endl << std::setprecision(-1);
}

printRegion2()

void InDet::TrackStatHelper::printRegion2	(	MsgStream &	out,
		enum eta_region	region,
		float	denominator
	)		const

Definition at line 695 of file TrackStatHelper.cxx.

                                                                                                       {
  out  <<  std::setiosflags(std::ios::fixed | std::ios::showpoint) << std::setw(7) << std::setprecision(2);
  if (denominator > 0) {     
    out << std::setiosflags(std::ios::fixed | std::ios::showpoint) << std::setw(4) << std::setprecision(1)
          <<"    \t\t" << float(m_hits.m_counter[kHits_rec][HIT_ALL    ][region]/denominator) << std::setprecision(2)
          << "\t" << float(m_hits.m_counter[kHits_rec][HIT_PIX1   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_rec][HIT_PIX2   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_rec][HIT_PIX3   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_rec][HIT_SCT1   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_rec][HIT_SCT2   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_rec][HIT_SCT3   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_rec][HIT_SCT4   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_rec][HIT_SCT5TO9][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_rec][HIT_TRT_ALL][region]/denominator)
          << std::endl << std::setprecision(-1);
  }
  else {
    out << " Unable to calculate: Denominator=0." << std::endl << std::setprecision(-1);
  } 
}

printRegionSecondary()

void InDet::TrackStatHelper::printRegionSecondary	(	MsgStream &	out,
		enum eta_region	region,
		float	denominator
	)		const

Prints ntracks per event,efficiencies,fake rates, and general hit information for given eta region.

Definition at line 808 of file TrackStatHelper.cxx.

                                                                           {
  out  <<  std::setiosflags(std::ios::fixed | std::ios::showpoint) << std::setw(7) << std::setprecision(2)
         << "\t" << m_tracks.m_counter[kTracks_gen][TRACK_MATCHED_SECONDARY][region]/(float) m_events;                                              // tracks / event
  if (m_tracks.m_counter[kTracks_gen][TRACK_SECONDARY][region]) {
    out << "\t" << m_tracks.m_counter[kTracks_gen][TRACK_MATCHED_SECONDARY][region]/ (float) m_tracks.m_counter[kTracks_gen][TRACK_SECONDARY][region];                 // track efficiency
  }
  else {
    out << "\t" << "n/a" ;
  }
 
  if (denominator > 0) {     
    out << std::setiosflags(std::ios::fixed | std::ios::showpoint) << std::setw(4) << std::setprecision(1)
          <<"   " << float(m_hits.m_counter[kHits_sec][HIT_ALL    ][region]/denominator) << std::setprecision(2)
          << "\t" << float(m_hits.m_counter[kHits_sec][HIT_PIX1   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_sec][HIT_PIX2   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_sec][HIT_PIX3   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_sec][HIT_SCT1   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_sec][HIT_SCT2   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_sec][HIT_SCT3   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_sec][HIT_SCT4   ][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_sec][HIT_SCT5TO9][region]/denominator)
          << "\t" << float(m_hits.m_counter[kHits_sec][HIT_TRT_ALL][region]/denominator)
          << std::endl << std::setprecision(-1);
  }
  else {
    out << " Unable to calculate: Denominator=0." << std::endl << std::setprecision(-1);
  } 
}

printSecondary()

void InDet::TrackStatHelper::printSecondary

(

MsgStream &

out

)

const

Prints all of the statistics information, calls printRegion, printTrackSummaryRegion, etc for detailed statistics.

Definition at line 760 of file TrackStatHelper.cxx.

                                                              {
  out << "TrackCollection \"" << key() << "\" " << std::endl;
  if (m_events > 0) {
    out << "(TrackAuthors:";
    std::vector<std::string> author_string;
    author_string.reserve(Trk::TrackInfo::NumberOfTrackFitters);
    for (int i=0; i<Trk::TrackInfo::NumberOfTrackFitters; i++) {
       author_string.push_back( Trk::TrackInfo(static_cast<Trk::TrackInfo::TrackFitter>(i),
                                               Trk::undefined,
                                               m_patternProperties,
                                               m_recoInfo).dumpInfo() );
    }
    for (int i=0; i<Trk::TrackInfo::NumberOfTrackFitters; i++){
      if (m_author_found[i]){
    out << " " << author_string[i];
      }
    }
    out << " )" << std::endl
          << "TrackTruthCollection \"" << Truthkey() << "\"" << std::endl;
 
    if (m_truthMissing)
      {
    if (m_careAboutTruth)
      out << " WARNING: TrackTruth missing for part of this TrackCollection --> No secondaries information printed!" << std::endl;
    else
      out << " INFO: Intentionally no TrackTruth for this TrackCollection. (No secondaries information printed.)" << std::endl;
      }
    if (!m_truthMissing)
      {
    out << " \t\t\t\t      ......................truth mached tracks statistics....................." << std::endl;
    out << " \t\t\t\tn/event\teff.\ttotal\tPIX1\tPIX2\tPIX3\tSCT1\tSCT2\tSCT3\tSCT4\tSCT5to9\tStraws" << std::endl;
    
    
    out << " total secondaries    ";
    printRegionSecondary(out, ETA_ALL,        (float) (m_tracks.m_counter[kTracks_gen][TRACK_MATCHED_SECONDARY][ETA_ALL]       +m_tracks.m_counter[kTracks_gen][TRACK_MULTMATCH_SECONDARY][ETA_ALL]) );
    out << " secondaries in barrel ";
    printRegionSecondary(out, ETA_BARREL,     (float) (m_tracks.m_counter[kTracks_gen][TRACK_MATCHED_SECONDARY][ETA_BARREL]    +m_tracks.m_counter[kTracks_gen][TRACK_MULTMATCH_SECONDARY][ETA_BARREL]) );
    out << " secondaries in trans. ";
    printRegionSecondary(out, ETA_TRANSITION, (float) (m_tracks.m_counter[kTracks_gen][TRACK_MATCHED_SECONDARY][ETA_TRANSITION]+m_tracks.m_counter[kTracks_gen][TRACK_MULTMATCH_SECONDARY][ETA_TRANSITION]) );
    out << " secondaries in endcap ";
    printRegionSecondary(out, ETA_ENDCAP,     (float) (m_tracks.m_counter[kTracks_gen][TRACK_MATCHED_SECONDARY][ETA_ENDCAP]    +m_tracks.m_counter[kTracks_gen][TRACK_MULTMATCH_SECONDARY][ETA_ENDCAP]) );
    out << " secondaries in forwa. ";
        printRegionSecondary(out, ETA_FORWARD,     (float) (m_tracks.m_counter[kTracks_gen][TRACK_MATCHED_SECONDARY][ETA_FORWARD]    +m_tracks.m_counter[kTracks_gen][TRACK_MULTMATCH_SECONDARY][ETA_FORWARD]) );
      }
  }
  else out << ": NO EVENTS PROCESSED! " << std::endl;
}

printTrackSummaryAverage()

void InDet::TrackStatHelper::printTrackSummaryAverage	(	MsgStream &	out,
		enum track_types	track_type,
		enum eta_region	eta_region,
		int	summary_type
	)		const

Prints information from TrackSummaryTool.

Definition at line 743 of file TrackStatHelper.cxx.

{
  assert( summary_type < kNSummaryTypes);
  long denom =  m_trackSummarySum.m_counter[kNTrackSummaryOK][track_type][eta_region][summary_type]
              + m_trackSummarySum.m_counter[kNTrackSummaryBAD][track_type][eta_region][summary_type];
 
  if (denom > 0) {
    out << float ( m_trackSummarySum.m_counter[kTrackSummarySum][track_type][eta_region][summary_type]/ (float) denom);
  }
  else {
    out << "n/a";
  }
}

printTrackSummaryRegion()

bool InDet::TrackStatHelper::printTrackSummaryRegion	(	MsgStream &	out,
		enum track_types	track_type,
		enum eta_region	eta_region
	)		const

Sets up detailed statistics part of table, calls printTrackSummaryRegion.

Definition at line 716 of file TrackStatHelper.cxx.

{
  long ex_denom =  m_trackSummarySum.m_counter[kNTrackSummaryOK][track_type][eta_region][kNumberOfPixelHits]
                 + m_trackSummarySum.m_counter[kNTrackSummaryBAD][track_type][eta_region][kNumberOfPixelHits];
  if (ex_denom > 0)
    {
      out << std::setiosflags(std::ios::fixed | std::ios::showpoint)
          << std::setw(8) << track_types_string[track_type]
          << std::setw(25) << key();
      
      if (m_events)
    out << std::setw(6)  << std::setprecision(2) << m_tracks.m_counter[kTracks_rec][track_type][eta_region]/(float) m_events;
      else
    out << std::setw(6) << "n/a";
      for (unsigned int stype_i=0; stype_i< kNSummaryTypes; ++stype_i) {
          out << std::setw(std::max(6,static_cast<int>(strlen(s_summaryTypeName[stype_i])+1)))
              << std::setprecision(2); printTrackSummaryAverage(out, track_type , eta_region, stype_i );
      }
      out << std::endl << std::setprecision(-1);
      return true; // yes, printed output
    }
  else
    return false; // no, didn't print output
}

reset()

void InDet::TrackStatHelper::reset

(

)

Resets the track collection information, called in the constructor.

Definition at line 571 of file TrackStatHelper.cxx.

                                {
  m_events = 0;
  m_hits.reset();
  m_tracks.reset();
  m_trackSummarySum.reset();
 
  for (int i=0; i<Trk::TrackInfo::NumberOfTrackFitters; i++) {
    m_author_found[i] = false;
  }
  m_truthMissing=false;
}

SetCuts()

void InDet::TrackStatHelper::SetCuts

(

const struct cuts &

ct

)

Sets the cuts such as the eta regions (barrel, transition,endcap) and the hit fraction fake cuts and the track matching cut.

Definition at line 113 of file TrackStatHelper.cxx.

{
  m_cuts = ct;
 
}

setSummaryStat()

void InDet::TrackStatHelper::setSummaryStat ( track_types  track_i, eta_region  region_i, const Trk::TrackSummary *  summary, TrackSummaryCounter &  trackSummarySum  ) const

inlineprivate

Definition at line 307 of file TrackStatHelper.h.

    {
       if (summary) {
          for (int stype=0; stype < kNSummaryTypes; stype++) {
             int value = summary->get(s_summaryTypes[stype]);
             //value is -1 if undefined
             if (value>0) {
                trackSummarySum.m_counter[kTrackSummarySum][track_i][region_i][stype] += value;
                trackSummarySum.m_counter[kNTrackSummaryOK][track_i][region_i][stype] ++;
             }
             else {
                trackSummarySum.m_counter[kNTrackSummaryBAD][track_i][region_i][stype] ++;
             }
          }
       }
    }

Truthkey()

const std::string& InDet::TrackStatHelper::Truthkey ( ) const

inline

Returns Truth TrackCollection Key.

Definition at line 179 of file TrackStatHelper.h.

{ return m_TrackTruthCollectionKey; }

Member Data Documentation

ATLAS_THREAD_SAFE [1/5]

TracksCounterAtomic m_tracks InDet::TrackStatHelper::ATLAS_THREAD_SAFE

mutableprivate

Definition at line 252 of file TrackStatHelper.h.

ATLAS_THREAD_SAFE [2/5]

HitsCounterAtomic m_hits InDet::TrackStatHelper::ATLAS_THREAD_SAFE

mutableprivate

Definition at line 262 of file TrackStatHelper.h.

ATLAS_THREAD_SAFE [3/5]

TrackSummaryCounterAtomic m_trackSummarySum InDet::TrackStatHelper::ATLAS_THREAD_SAFE

mutableprivate

Definition at line 304 of file TrackStatHelper.h.

ATLAS_THREAD_SAFE [4/5]

std::bitset<Trk::TrackInfo::NumberOfTrackRecoInfo> m_recoInfo InDet::TrackStatHelper::ATLAS_THREAD_SAFE

mutableprivate

Definition at line 332 of file TrackStatHelper.h.

ATLAS_THREAD_SAFE [5/5]

std::bitset<Trk::TrackInfo::NumberOfTrackProperties> m_patternProperties InDet::TrackStatHelper::ATLAS_THREAD_SAFE

mutableprivate

Definition at line 333 of file TrackStatHelper.h.

m_author_found

std::atomic<bool> InDet::TrackStatHelper::m_author_found[Trk::TrackInfo::NumberOfTrackFitters] {}

mutableprivate

Number of tracking authors found.

Definition at line 186 of file TrackStatHelper.h.

m_authorMutex

std::mutex InDet::TrackStatHelper::m_authorMutex

mutableprivate

Definition at line 331 of file TrackStatHelper.h.

m_careAboutTruth

bool InDet::TrackStatHelper::m_careAboutTruth

private

Definition at line 328 of file TrackStatHelper.h.

m_cuts

struct cuts InDet::TrackStatHelper::m_cuts

private

Definition at line 328 of file TrackStatHelper.h.

m_events

std::atomic<long> InDet::TrackStatHelper::m_events {}

mutableprivate

Number of events.

Definition at line 187 of file TrackStatHelper.h.

m_TrackCollectionKey

std::string InDet::TrackStatHelper::m_TrackCollectionKey

private

StoreGate Track Collection Key.

Definition at line 183 of file TrackStatHelper.h.

m_TrackTruthCollectionKey

std::string InDet::TrackStatHelper::m_TrackTruthCollectionKey

private

StoreGate Track Truth Collection Key.

Definition at line 184 of file TrackStatHelper.h.

m_truthMissing

std::atomic<bool> InDet::TrackStatHelper::m_truthMissing

mutableprivate

Flag for if track truth is missing.

Definition at line 327 of file TrackStatHelper.h.

s_summaryTypeName

const char *const InDet::TrackStatHelper::s_summaryTypeName

staticprivate

Initial value:

= {
   "blay",
   "pix",
   "hole",
   "gang",
   "SCT",
   "hole",
   "DHole",
   "TRT",
   "outl",
   "TRHi",
   "outl",
   "alloutl"}

table column labels for summary

Definition at line 291 of file TrackStatHelper.h.

s_summaryTypes

const Trk::SummaryType InDet::TrackStatHelper::s_summaryTypes

staticprivate

Initial value:

= {
   Trk::numberOfInnermostPixelLayerHits,
   Trk::numberOfPixelHits,
   Trk::numberOfPixelHoles,
   Trk::numberOfGangedPixels,
   Trk::numberOfSCTHits,
   Trk::numberOfSCTHoles,
   Trk::numberOfSCTDoubleHoles,
   Trk::numberOfTRTHits,
   Trk::numberOfTRTOutliers,
   Trk::numberOfTRTHighThresholdHits,
   Trk::numberOfTRTHighThresholdOutliers,
   Trk::numberOfOutliersOnTrack,
}

summary types for which statistics are gathered

Definition at line 289 of file TrackStatHelper.h.

---

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

• TrackStatHelper.h
• TrackStatHelper.cxx