This validation tool writes basic information about Trk::TrackParticleBase into an ntuple;. More...

#include <BasicValTrkParticleNtupleTool.h>

Inheritance diagram for Trk::BasicValTrkParticleNtupleTool:

Collaboration diagram for Trk::BasicValTrkParticleNtupleTool:

Public Member Functions
	BasicValTrkParticleNtupleTool (const std::string &, const std::string &, const IInterface *)
	~BasicValTrkParticleNtupleTool ()
StatusCode	initialize ()
StatusCode	finalize ()
	finalize
virtual StatusCode	addNtupleItems (TTree *tree)
	add branches to the tree Should be called once (per track collection and tree) dunring the initialisation phase by the calling algorithm (usually Trk::TrackValidationNtupleWriter)
virtual StatusCode	writeTrackData (const Trk::Track &, const int iterationIndex, const unsigned int fitStatCode)
	fill AND write ntuple data of a given track
virtual StatusCode	writeTrackParticleData (const Trk::TrackParticleBase &)
	fill AND write ntuple data of a given track particle
virtual StatusCode	fillTrackData (const Trk::Track &, const int iterationIndex, const unsigned int fitStatCode)
	fill ntuple data of a given track without writing the record.
virtual StatusCode	fillTrackParticleData (const Trk::TrackParticleBase &)
	fill ntuple data of a given track particle without writing the record.
virtual StatusCode	fillTrackParameter (const Trk::TrackParameters *, const int iterationIndex)
	fill ntuple data of a given track without writing the record.
virtual StatusCode	fillMeasurementData (const Trk::MeasurementBase , const Trk::TrackParameters )
	fill ntuple data of given measurement and track parameters without writing the record
virtual StatusCode	fillOutlierData (const Trk::MeasurementBase , const Trk::TrackParameters , const Trk::FitQualityOnSurface *)
	fill ntuple data of an outlier measurement (without writing the record yet).
virtual StatusCode	fillHoleData (const Trk::TrackStateOnSurface &)
	fill ntuple data of holes on track without writing the record
virtual StatusCode	fillTrackTruthData (const TrackParameters *&, const TrackTruth &, const int truthIndex=-1)
	fill ntuple data of holes on track without writing the record
virtual StatusCode	writeRecord (TTree *tree)
	write the filled data into the ntuple
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()
	Interface ID, declared here, and defined below.

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

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	fillTrkParticlePerigee (const Trk::Perigee *perigee)
	fill the perigee in ntuple
StatusCode	fillTrkParticleSummary (const Trk::TrackSummary *summary)
	fill the perigee in ntuple
StatusCode	fillFitQualityData (const Trk::FitQuality *fitQuality)
void	resetVariables ()
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
bool	m_bookNewNtuple
	jobOption: book new ntuple?
std::string	m_ntupleFileName
	jobOption: Ntuple file name
std::string	m_ntupleDirName
	jobOption: Ntuple directory name
std::string	m_ntupleTreeName
	jobOption: Ntuple tree name
TTree *	m_nt
	Pointer to the NTuple tree.
int	m_TrackIDcounter
	help variable to link with event property tree
unsigned int	m_lastEventNumber
	help variable to link with event property tree
int	m_runNumber
	run number the track belongs to
int	m_eventNumber
	event number the track belongs to
unsigned char	m_TrackID
	number of the track within the current event
float	m_Rec_d0
	reconstructed track params: local d0 at perigee UNIT:mm
float	m_Rec_z0
	reconstructed track params: local z0 at perigee UNIT:mm
float	m_Rec_phi0
	reconstructed track params: angle phi at perigee UNIT:rad
float	m_Rec_eta
	reconstructed track params: pseudorapidity UNIT:1
float	m_Rec_qOverP
	reconstructed track params: curvature q/p UNIT:1/MeV
float	m_errord0
	error on local d0 at perigee UNIT:mm
float	m_errorz0
	error on local z0 at perigee UNIT:mm
float	m_errorphi0
	error on angle phi at perigee UNIT:rad
float	m_errortheta0
	error on angle theta at perigee UNIT:rad
float	m_errorqoverp
	error on curvature q/p UNIT:1/MeV
float	m_chi2
	chi sqared of the track fit normalized to DOF
unsigned char	m_numberOfPixelHits
	number of measurements in the Pixels
unsigned char	m_numberOfContribPixelLayers
	number of contributing layers of the pixel detector
unsigned char	m_numberOfPixelHoles
	number of pixel layers on track with absence of hits
unsigned char	m_numberOfPixelDeadSensors
	number of dead pixel sensors crossed
unsigned char	m_numberOfSCTHits
	number of measurements in the SCT
unsigned char	m_numberOfSCTHoles
	number of SCT holes
unsigned char	m_numberOfSCTDeadSensors
	number of dead SCT sensors crossed
unsigned char	m_numberOfTRTHits
	number of TRT hits
unsigned char	m_numberOfTRTHoles
	number of TRT holes
unsigned char	m_numberOfTRTDeadStraws
	number of dead TRT straws crossed
unsigned char	m_numberOfTRTHighThresholdHits
	number of TRT hits which pass the high threshold
unsigned long	m_idHitPattern
	bit word carrying information about hit layers
SG::ReadHandleKey< xAOD::EventInfo >	m_evt {this, "EvtInfo", "EventInfo", "EventInfo name"}
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

Static Private Attributes
static const float	s_errorEntry = -99999.0
	error entry costant

Detailed Description

This validation tool writes basic information about Trk::TrackParticleBase into an ntuple;.

It can book the ntuple, create the tree directory and add basic entries. This tool is also able to write the current record to the ntuple file and start a new entry. You can also use an existing ntuple tree (which has to have the same track-wise structure) and just use the Trk::BasicValidationNtupleTool to add some ntuple variables.

Every entry (Rec:TrackPartcile) contains the event number, a track id (numbering the tracks within one event).

Perigee parameters and track summary are filled for every track if available.

Author: Pavel.nosp@m..Ruz.nosp@m.icka@.nosp@m.cern.nosp@m..ch

Definition at line 39 of file BasicValTrkParticleNtupleTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ BasicValTrkParticleNtupleTool()

Trk::BasicValTrkParticleNtupleTool::BasicValTrkParticleNtupleTool	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p )

Definition at line 27 of file BasicValTrkParticleNtupleTool.cxx.

        :
        AthAlgTool(t,n,p),         // retrieve as private tools
        m_bookNewNtuple{},
        m_nt(nullptr),
        m_TrackIDcounter{},
        m_lastEventNumber(0),
        m_runNumber{},
        m_eventNumber{},
        m_TrackID{},
        m_Rec_d0{},     
         m_Rec_z0{},     
         m_Rec_phi0{},   
         m_Rec_eta{},    
         m_Rec_qOverP{}, 
 
         m_errord0{},    
         m_errorz0{},    
         m_errorphi0{},  
         m_errortheta0{},
         m_errorqoverp{},
         m_chi2{},
        m_numberOfPixelHits{},           
         m_numberOfContribPixelLayers{},  
         m_numberOfPixelHoles{},          
         m_numberOfPixelDeadSensors{},    
         m_numberOfSCTHits{},             
         m_numberOfSCTHoles{},            
         m_numberOfSCTDeadSensors{},      
         m_numberOfTRTHits{},             
         m_numberOfTRTHoles{},            
         m_numberOfTRTDeadStraws{},       
         m_numberOfTRTHighThresholdHits{},
         m_idHitPattern{}    
{
   
    declareInterface<IValidationNtupleTool>(this);
    // Declare the properties
    declareProperty("BookNewNtuple",                m_bookNewNtuple = false,                "Create the ntuple tree or use existing one?");
    declareProperty("NtupleFileName",               m_ntupleFileName = "/NTUPLES/FILE1",    "Ntuple file handle");
    declareProperty("NtupleDirectoryName",          m_ntupleDirName = "FitterValidation",   "dircetory name for ntuple tree");
    declareProperty("NtupleTreeName",               m_ntupleTreeName = "Validation",        "Name of the ntuple tree");
}

◆ ~BasicValTrkParticleNtupleTool()

Trk::BasicValTrkParticleNtupleTool::~BasicValTrkParticleNtupleTool ( )

default

Member Function Documentation

◆ addNtupleItems()

StatusCode Trk::BasicValTrkParticleNtupleTool::addNtupleItems ( TTree * tree )

virtual

add branches to the tree Should be called once (per track collection and tree) dunring the initialisation phase by the calling algorithm (usually Trk::TrackValidationNtupleWriter)

Implements Trk::IValidationNtupleTool.

Definition at line 123 of file BasicValTrkParticleNtupleTool.cxx.

                                                                         {
    if (!tree) return StatusCode::FAILURE;
    //-----------------
    // add items  *** Note: Documentation is in the header file, doxygen and wikis! ***
    // 
    // event info:
    tree->Branch("RunNumber",           &m_runNumber,   "run_number/I");
    tree->Branch("EventNumber",         &m_eventNumber, "event_number/I");
    tree->Branch("TrackID",             &m_TrackID,     "ID_of_track_in_event/b");
 
    // reconstructed track parameters
    tree->Branch("RecD0",               &m_Rec_d0,      "Reconstructed_d0/F");
    tree->Branch("RecZ0",               &m_Rec_z0,      "Reconstructed_z0/F");
    tree->Branch("RecPhi0",             &m_Rec_phi0,    "Reconstructed_phi0/F");
    tree->Branch("RecEta",              &m_Rec_eta,     "Reconstructed_eta/F");
    tree->Branch("RecQoverP",           &m_Rec_qOverP,  "Reconstructed_Q_over_p/F");
    
    // errors of reconstructed track parameters
    tree->Branch("RecErrD0",            &m_errord0,      "err_d0/F");
    tree->Branch("RecErrZ0",            &m_errorz0,      "err_z0/F");
    tree->Branch("RecErrPhi0",          &m_errorphi0,    "err_phi0/F");
    tree->Branch("RecErrTheta",         &m_errortheta0,  "err_theta/F");
    tree->Branch("RecErrQoverP",        &m_errorqoverp,  "err_Q_over_p/F");
   
    // fit quality parameters 
    tree->Branch("chi2",                &m_chi2,                "chi2/F");
 
    //track particle summary parameters
    tree->Branch("nPixelHits",            &m_numberOfPixelHits,            "numberOfPixelHits/b");
    tree->Branch("nContribPixelLayers",   &m_numberOfContribPixelLayers,   "numberOfContribPixelLayers/b");
    tree->Branch("nPixelHoles",           &m_numberOfPixelHoles,           "numberOfPixelHoles/b");
    tree->Branch("nPixelDeadSensors",     &m_numberOfPixelDeadSensors,     "numberOfPixelDeadSensors/b");
    tree->Branch("nSCTHits",              &m_numberOfSCTHits,              "numberOfSCTHits/b");
    tree->Branch("nSCTHoles",             &m_numberOfSCTHoles,             "numberOfSCTHoles/b");
    tree->Branch("nSCTDeadSensors"    ,   &m_numberOfSCTDeadSensors,       "numberOfSCTDeadSensors/b");
    tree->Branch("nTRTHits"  ,            &m_numberOfTRTHits,              "numberOfTRTHits/b");
    tree->Branch("nTRTHoles",             &m_numberOfTRTHoles,             "numberOfTRTHoles/b");
    tree->Branch("nTRTDeadStraws",        &m_numberOfTRTDeadStraws,        "numberOfTRTDeadStraws/b");
    tree->Branch("nTRTHighThresholdHits", &m_numberOfTRTHighThresholdHits, "numberOfTRTHighThresholdHits/b");
    
    tree->Branch("idHitPattern",          &m_idHitPattern,                 "idHitPattern/l");
 
     ATH_MSG_DEBUG ("added own branches to ntuple");
 
    return StatusCode::SUCCESS;
}

◆ declareGaudiProperty()

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareProperty()

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

95{ return m_detStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ evtStore()

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

85{ return m_evtStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ fillFitQualityData()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillFitQualityData ( const Trk::FitQuality * fitQuality )

private

Definition at line 311 of file BasicValTrkParticleNtupleTool.cxx.

                                                                                               {
  if (!fitQuality) {
        ATH_MSG_WARNING("Something is wrong - track has no fit quality data !!");
        m_chi2  = 0;
                   
        return StatusCode::FAILURE;
      }
 
  if(fitQuality->numberDoF() == 0) {
                                    ATH_MSG_WARNING("Number of DOF is zero !! Could not normalize chi2 ");
                    return StatusCode::FAILURE;
                                    }
  
  else m_chi2  = (float)fitQuality->chiSquared()/(float)fitQuality->numberDoF();
        
   return StatusCode::SUCCESS;
}

◆ fillHoleData()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillHoleData ( const Trk::TrackStateOnSurface & )

virtual

fill ntuple data of holes on track without writing the record

Implements Trk::IValidationNtupleTool.

Definition at line 382 of file BasicValTrkParticleNtupleTool.cxx.

383 {return StatusCode::SUCCESS;}

◆ fillMeasurementData()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillMeasurementData	(	const Trk::MeasurementBase *	,
		const Trk::TrackParameters *	)

virtual

fill ntuple data of given measurement and track parameters without writing the record

Implements Trk::IValidationNtupleTool.

Definition at line 371 of file BasicValTrkParticleNtupleTool.cxx.

373 {return StatusCode::SUCCESS;}

◆ fillOutlierData()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillOutlierData	(	const Trk::MeasurementBase *	,
		const Trk::TrackParameters *	,
		const Trk::FitQualityOnSurface *	)

virtual

fill ntuple data of an outlier measurement (without writing the record yet).

The TrackParameters and FQoS are optional, i.e. NULL pointers should be handled by the implementation

Implements Trk::IValidationNtupleTool.

Definition at line 376 of file BasicValTrkParticleNtupleTool.cxx.

379 {return StatusCode::SUCCESS;}

◆ fillTrackData()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillTrackData	(	const Trk::Track &	,
		const int	iterationIndex,
		const unsigned int	fitStatCode )

virtual

fill ntuple data of a given track without writing the record.

if this method is called twice without writing the ntuple inbetween the first data will be lost!

Implements Trk::IValidationNtupleTool.

Definition at line 359 of file BasicValTrkParticleNtupleTool.cxx.

362 {return StatusCode::SUCCESS;}

◆ fillTrackParameter()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillTrackParameter	(	const Trk::TrackParameters *	,
		const int	iterationIndex )

virtual

fill ntuple data of a given track without writing the record.

if this method is called twice without writing the ntuple inbetween the first data will be lost!

Implements Trk::IValidationNtupleTool.

Definition at line 366 of file BasicValTrkParticleNtupleTool.cxx.

368 {return StatusCode::SUCCESS;}

◆ fillTrackParticleData()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillTrackParticleData ( const Trk::TrackParticleBase & track )

virtual

fill ntuple data of a given track particle without writing the record.

fill track data into variables without actually writing the record

if this method is called twice without writing the ntuple inbetween the first data will be lost!

Implements Trk::IValidationNtupleTool.

Definition at line 190 of file BasicValTrkParticleNtupleTool.cxx.

                                       {
 
  // ---------------------------------------
    // detect new event, reset TrackParticle counter if new event
  SG::ReadHandle<xAOD::EventInfo> evt(m_evt);
    if (!evt.isValid()) {
      ATH_MSG_WARNING("Could not retrieve event info");
      m_runNumber   = (int)s_errorEntry;
      m_eventNumber = (int)s_errorEntry;
      return StatusCode::FAILURE;
      }
  
   if (m_lastEventNumber!=evt->eventNumber()) {
        // we have a new event, reset TrackParticleID:
        m_TrackIDcounter = 0;
        m_lastEventNumber = evt->eventNumber();
      }
    m_TrackIDcounter++;
    m_TrackID = (unsigned char)m_TrackIDcounter;
    m_eventNumber = evt->eventNumber();
    m_runNumber   = evt->runNumber(); 
   
    ATH_MSG_VERBOSE ("Event: " << m_eventNumber << ", Run: "<< m_runNumber  << " TrackID: " << m_TrackID);
            
    //----------------------------------------------
    // fill track parameters in ntuple
    const Trk::Perigee* perpars = track.perigee();
    if (perpars != nullptr && fillTrkParticlePerigee(perpars).isFailure())  ATH_MSG_WARNING("Perigee parameters could not be written to ntuple");
    
    const Trk::TrackSummary* summary = track.trackSummary();
    if((!summary) || fillTrkParticleSummary(summary).isFailure()) ATH_MSG_WARNING("Summary parameters could not be written to ntuple");
 
    const Trk::FitQuality* fitQuality = track.fitQuality();
    if((!fitQuality) || fillFitQualityData(fitQuality).isFailure() ) ATH_MSG_WARNING("Fit Quality data could not be written to ntuple");
 
    return StatusCode::SUCCESS;
}

◆ fillTrackTruthData()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillTrackTruthData	(	const TrackParameters *&	,
		const TrackTruth &	,
		const int	truthIndex = -1 )

virtual

fill ntuple data of holes on track without writing the record

Implements Trk::IValidationNtupleTool.

Definition at line 386 of file BasicValTrkParticleNtupleTool.cxx.

389 {return StatusCode::SUCCESS;}

◆ fillTrkParticlePerigee()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillTrkParticlePerigee ( const Trk::Perigee * perigee )

private

fill the perigee in ntuple

Definition at line 233 of file BasicValTrkParticleNtupleTool.cxx.

                                                                                             {
 
  ATH_MSG_VERBOSE ("in fillTrackPerigee");
 
    if (!perigee) {
        ATH_MSG_WARNING("Something is wrong - track has no perigee at all!");
        m_Rec_d0    = 0;
        m_Rec_z0    = 0;
        m_Rec_phi0  = 0;
        m_Rec_eta = 0;
        m_Rec_qOverP = 0;
        return StatusCode::FAILURE;
    }
 
    ATH_MSG_VERBOSE("Fill Perigee Parameters now");
    // track parameters
    m_Rec_d0    = perigee->parameters()[Trk::d0];
    m_Rec_z0    = perigee->parameters()[Trk::z0];
    m_Rec_phi0  = perigee->parameters()[Trk::phi0];
    m_Rec_eta   = perigee->eta();
    m_Rec_qOverP = perigee->parameters()[Trk::qOverP];
    
    if (perigee->covariance()) {
 
        // error of track parameters
        m_errord0       = Amg::error(*perigee->covariance(),Trk::d0);
        m_errorz0       = Amg::error(*perigee->covariance(),Trk::z0);
        m_errorphi0     = Amg::error(*perigee->covariance(),Trk::phi0);
        m_errortheta0   = Amg::error(*perigee->covariance(),Trk::theta);
        m_errorqoverp   = Amg::error(*perigee->covariance(),Trk::qOverP);
 
 
    } // end if measPerigee
    ATH_MSG_DEBUG ("Trackparameters: d0=" << m_Rec_d0 << ", z0=" << m_Rec_z0 << ", phi0=" << m_Rec_phi0 << ", eta=" << m_Rec_eta);
 
    return StatusCode::SUCCESS;
}

◆ fillTrkParticleSummary()

StatusCode Trk::BasicValTrkParticleNtupleTool::fillTrkParticleSummary ( const Trk::TrackSummary * summary )

private

fill the perigee in ntuple

Definition at line 274 of file BasicValTrkParticleNtupleTool.cxx.

                                                                                                  {
 
  ATH_MSG_VERBOSE ("in fillTrackSummary");
 
    if (!summary) {
        ATH_MSG_WARNING("Something is wrong - track has no summary at all!");
        m_numberOfPixelHits  = 0;
        m_numberOfSCTHits    = 0;
 
        return StatusCode::FAILURE;
    }
 
    ATH_MSG_VERBOSE("Fill Summary Parameters now");
    // summary of the track
    m_numberOfPixelHits            = (unsigned char)summary->get(Trk::numberOfPixelHits);
    m_numberOfContribPixelLayers   = (unsigned char)summary->get(Trk::numberOfContribPixelLayers);
    m_numberOfPixelHoles           = (unsigned char)summary->get(Trk::numberOfPixelHoles);
    m_numberOfPixelDeadSensors     = (unsigned char)summary->get(Trk::numberOfPixelDeadSensors);
    m_numberOfSCTHits              = (unsigned char)summary->get(Trk::numberOfSCTHits); 
    m_numberOfSCTHoles             = (unsigned char)summary->get(Trk::numberOfSCTHoles);
    m_numberOfSCTDeadSensors       = (unsigned char)summary->get(Trk::numberOfSCTDeadSensors);
    m_numberOfTRTHits              = (unsigned char)summary->get(Trk::numberOfTRTHits);
    m_numberOfTRTHoles             = (unsigned char)summary->get(Trk::numberOfTRTHoles);
    m_numberOfTRTDeadStraws        = (unsigned char)summary->get(Trk::numberOfTRTDeadStraws);
    m_numberOfTRTHighThresholdHits = (unsigned char)summary->get(Trk::numberOfTRTHighThresholdHits);
    
    std::bitset<Trk::numberOfDetectorTypes> hitPattern; // all position are set to 0
    for (int i=0; i<Trk::numberOfDetectorTypes; ++i) {
      if (summary->isHit(static_cast<Trk::DetectorType>(i))) hitPattern.set(i,1); // set position i to 1
                              }
    m_idHitPattern = hitPattern.to_ulong();
  
    ATH_MSG_DEBUG ("Track summary: number of Pixel hit=" << m_numberOfPixelHits << ", number of SCT hits=" << m_numberOfSCTHits );
 
    return StatusCode::SUCCESS;
}

◆ finalize()

StatusCode Trk::BasicValTrkParticleNtupleTool::finalize ( )

finalize

Definition at line 112 of file BasicValTrkParticleNtupleTool.cxx.

                                                      {
 
  ATH_MSG_DEBUG("start finalize() in " << name());
  if (m_nt) {
    delete m_nt;
    m_nt = nullptr;
  }
 
  return StatusCode::SUCCESS;
}

◆ initialize()

StatusCode Trk::BasicValTrkParticleNtupleTool::initialize ( )

Definition at line 81 of file BasicValTrkParticleNtupleTool.cxx.

                                                        {
 
  ATH_CHECK( m_evt.initialize() );
 
    // create ntuple tree
  if (m_bookNewNtuple) {
    // ---------------------------
    // retrieve pointer to THistSvc
    SmartIF<ITHistSvc> tHistSvc{service("THistSvc")};
    ATH_CHECK( tHistSvc.isValid() );
 
    // ---------------------------
    // create tree and register it to THistSvc
 
    m_nt = new TTree(TString(m_ntupleTreeName), "Track Particle Validation");
    // NB: we must not delete the tree, this is done by THistSvc
 
    std::string fullNtupleName =  "/"+m_ntupleFileName+"/"+m_ntupleDirName+"/"+m_ntupleTreeName;
    ATH_CHECK( tHistSvc->regTree(fullNtupleName, m_nt) );
 
    // add the ntuple branches (this function has to be called by the client of this tool, if m_bookNewNtuple is set to false...)
    ATH_CHECK(addNtupleItems(m_nt));
  }
   
    m_TrackIDcounter = 0;
    return StatusCode::SUCCESS;
}

◆ inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

const InterfaceID & Trk::IValidationNtupleTool::interfaceID ( )

inlinestaticinherited

Interface ID, declared here, and defined below.

Definition at line 105 of file IValidationNtupleTool.h.

                                                                {
    return IID_IValidationNtupleTool;
}

◆ msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

◆ outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

◆ resetVariables()

void Trk::BasicValTrkParticleNtupleTool::resetVariables ( )

private

Definition at line 337 of file BasicValTrkParticleNtupleTool.cxx.

                                                      {
   m_numberOfPixelHits=0;  
   m_numberOfContribPixelLayers=0;
   m_numberOfPixelHoles=0;
   m_numberOfPixelDeadSensors=0;
   m_numberOfSCTHits=0; 
   m_numberOfSCTHoles=0; 
   m_numberOfSCTDeadSensors=0;
   m_numberOfTRTHits=0;
   m_numberOfTRTHoles=0;
   m_numberOfTRTDeadStraws=0;
   m_numberOfTRTHighThresholdHits=0;
   m_idHitPattern=0;
}

◆ sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

overridevirtualinherited

Handle START transition.

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

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

◆ writeRecord()

StatusCode Trk::BasicValTrkParticleNtupleTool::writeRecord ( TTree * tree )

virtual

write the filled data into the ntuple

Implements Trk::IValidationNtupleTool.

Definition at line 329 of file BasicValTrkParticleNtupleTool.cxx.

                                                                    {
    if (!tree) return StatusCode::FAILURE;
    ATH_MSG_VERBOSE ("Writting Track Particles into the Tree");
    tree->Fill();
    resetVariables();
    return StatusCode::SUCCESS;
}

◆ writeTrackData()

StatusCode Trk::BasicValTrkParticleNtupleTool::writeTrackData	(	const Trk::Track &	,
		const int	iterationIndex,
		const unsigned int	fitStatCode )

virtual

fill AND write ntuple data of a given track

Implements Trk::IValidationNtupleTool.

Definition at line 352 of file BasicValTrkParticleNtupleTool.cxx.

355 {return StatusCode::SUCCESS;}

◆ writeTrackParticleData()

StatusCode Trk::BasicValTrkParticleNtupleTool::writeTrackParticleData ( const Trk::TrackParticleBase & track )

virtual

fill AND write ntuple data of a given track particle

Implements Trk::IValidationNtupleTool.

Definition at line 171 of file BasicValTrkParticleNtupleTool.cxx.

                                        {
    if (!m_nt) {
        ATH_MSG_ERROR("writeTrackParticleData(...) can only be used, if property BookNewNtuple is set to true" );
        return StatusCode::FAILURE;
    }
    StatusCode sc;
 
    ATH_MSG_DEBUG ("in writeTrackParticleData(trk, indx)");
 
    sc = fillTrackParticleData(track);
    if (sc.isFailure()) return sc;
    
    return writeRecord(m_nt);
}

Member Data Documentation

◆ m_bookNewNtuple

bool Trk::BasicValTrkParticleNtupleTool::m_bookNewNtuple

private

jobOption: book new ntuple?

Definition at line 114 of file BasicValTrkParticleNtupleTool.h.

◆ m_chi2

float Trk::BasicValTrkParticleNtupleTool::m_chi2

private

chi sqared of the track fit normalized to DOF

Definition at line 143 of file BasicValTrkParticleNtupleTool.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_errord0

float Trk::BasicValTrkParticleNtupleTool::m_errord0

private

error on local d0 at perigee UNIT:mm

Definition at line 136 of file BasicValTrkParticleNtupleTool.h.

◆ m_errorphi0

float Trk::BasicValTrkParticleNtupleTool::m_errorphi0

private

error on angle phi at perigee UNIT:rad

Definition at line 138 of file BasicValTrkParticleNtupleTool.h.

◆ m_errorqoverp

float Trk::BasicValTrkParticleNtupleTool::m_errorqoverp

private

error on curvature q/p UNIT:1/MeV

Definition at line 140 of file BasicValTrkParticleNtupleTool.h.

◆ m_errortheta0

float Trk::BasicValTrkParticleNtupleTool::m_errortheta0

private

error on angle theta at perigee UNIT:rad

Definition at line 139 of file BasicValTrkParticleNtupleTool.h.

◆ m_errorz0

float Trk::BasicValTrkParticleNtupleTool::m_errorz0

private

error on local z0 at perigee UNIT:mm

Definition at line 137 of file BasicValTrkParticleNtupleTool.h.

◆ m_eventNumber

int Trk::BasicValTrkParticleNtupleTool::m_eventNumber

private

event number the track belongs to

Definition at line 126 of file BasicValTrkParticleNtupleTool.h.

◆ m_evt

SG::ReadHandleKey<xAOD::EventInfo> Trk::BasicValTrkParticleNtupleTool::m_evt {this, "EvtInfo", "EventInfo", "EventInfo name"}

private

Definition at line 162 of file BasicValTrkParticleNtupleTool.h.

162{this, "EvtInfo", "EventInfo", "EventInfo name"};

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_idHitPattern

unsigned long Trk::BasicValTrkParticleNtupleTool::m_idHitPattern

private

bit word carrying information about hit layers

Definition at line 158 of file BasicValTrkParticleNtupleTool.h.

◆ m_lastEventNumber

unsigned int Trk::BasicValTrkParticleNtupleTool::m_lastEventNumber

private

help variable to link with event property tree

Definition at line 122 of file BasicValTrkParticleNtupleTool.h.

◆ m_nt

TTree* Trk::BasicValTrkParticleNtupleTool::m_nt

private

Pointer to the NTuple tree.

Definition at line 119 of file BasicValTrkParticleNtupleTool.h.

◆ m_ntupleDirName

std::string Trk::BasicValTrkParticleNtupleTool::m_ntupleDirName

private

jobOption: Ntuple directory name

Definition at line 116 of file BasicValTrkParticleNtupleTool.h.

◆ m_ntupleFileName

std::string Trk::BasicValTrkParticleNtupleTool::m_ntupleFileName

private

jobOption: Ntuple file name

Definition at line 115 of file BasicValTrkParticleNtupleTool.h.

◆ m_ntupleTreeName

std::string Trk::BasicValTrkParticleNtupleTool::m_ntupleTreeName

private

jobOption: Ntuple tree name

Definition at line 117 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfContribPixelLayers

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfContribPixelLayers

private

number of contributing layers of the pixel detector

Definition at line 147 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfPixelDeadSensors

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfPixelDeadSensors

private

number of dead pixel sensors crossed

Definition at line 149 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfPixelHits

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfPixelHits

private

number of measurements in the Pixels

Definition at line 146 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfPixelHoles

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfPixelHoles

private

number of pixel layers on track with absence of hits

Definition at line 148 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfSCTDeadSensors

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfSCTDeadSensors

private

number of dead SCT sensors crossed

Definition at line 152 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfSCTHits

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfSCTHits

private

number of measurements in the SCT

Definition at line 150 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfSCTHoles

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfSCTHoles

private

number of SCT holes

Definition at line 151 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfTRTDeadStraws

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfTRTDeadStraws

private

number of dead TRT straws crossed

Definition at line 155 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfTRTHighThresholdHits

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfTRTHighThresholdHits

private

number of TRT hits which pass the high threshold

Definition at line 156 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfTRTHits

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfTRTHits

private

number of TRT hits

Definition at line 153 of file BasicValTrkParticleNtupleTool.h.

◆ m_numberOfTRTHoles

unsigned char Trk::BasicValTrkParticleNtupleTool::m_numberOfTRTHoles

private

number of TRT holes

Definition at line 154 of file BasicValTrkParticleNtupleTool.h.

◆ m_Rec_d0

float Trk::BasicValTrkParticleNtupleTool::m_Rec_d0

private

reconstructed track params: local d0 at perigee UNIT:mm

Definition at line 130 of file BasicValTrkParticleNtupleTool.h.

◆ m_Rec_eta

float Trk::BasicValTrkParticleNtupleTool::m_Rec_eta

private

reconstructed track params: pseudorapidity UNIT:1

Definition at line 133 of file BasicValTrkParticleNtupleTool.h.

◆ m_Rec_phi0

float Trk::BasicValTrkParticleNtupleTool::m_Rec_phi0

private

reconstructed track params: angle phi at perigee UNIT:rad

Definition at line 132 of file BasicValTrkParticleNtupleTool.h.

◆ m_Rec_qOverP

float Trk::BasicValTrkParticleNtupleTool::m_Rec_qOverP

private

reconstructed track params: curvature q/p UNIT:1/MeV

Definition at line 134 of file BasicValTrkParticleNtupleTool.h.

◆ m_Rec_z0

float Trk::BasicValTrkParticleNtupleTool::m_Rec_z0

private

reconstructed track params: local z0 at perigee UNIT:mm

Definition at line 131 of file BasicValTrkParticleNtupleTool.h.

◆ m_runNumber

int Trk::BasicValTrkParticleNtupleTool::m_runNumber

private

run number the track belongs to

Definition at line 125 of file BasicValTrkParticleNtupleTool.h.

◆ m_TrackID

unsigned char Trk::BasicValTrkParticleNtupleTool::m_TrackID

private

number of the track within the current event

Definition at line 127 of file BasicValTrkParticleNtupleTool.h.

◆ m_TrackIDcounter

int Trk::BasicValTrkParticleNtupleTool::m_TrackIDcounter

private

help variable to link with event property tree

Definition at line 121 of file BasicValTrkParticleNtupleTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ s_errorEntry

const float Trk::BasicValTrkParticleNtupleTool::s_errorEntry = -99999.0

staticprivate

error entry costant

Definition at line 160 of file BasicValTrkParticleNtupleTool.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ BasicValTrkParticleNtupleTool()

◆ ~BasicValTrkParticleNtupleTool()

Member Function Documentation

◆ addNtupleItems()

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ extraDeps_update_handler()

◆ fillFitQualityData()

◆ fillHoleData()

◆ fillMeasurementData()

◆ fillOutlierData()

◆ fillTrackData()

◆ fillTrackParameter()

◆ fillTrackParticleData()

◆ fillTrackTruthData()

◆ fillTrkParticlePerigee()

◆ fillTrkParticleSummary()

◆ finalize()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ resetVariables()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

◆ writeRecord()

◆ writeTrackData()

◆ writeTrackParticleData()

Member Data Documentation

◆ m_bookNewNtuple

◆ m_chi2

◆ m_detStore

◆ m_errord0

◆ m_errorphi0

◆ m_errorqoverp

◆ m_errortheta0

◆ m_errorz0

◆ m_eventNumber

◆ m_evt

◆ m_evtStore

◆ m_idHitPattern

◆ m_lastEventNumber

◆ m_nt

◆ m_ntupleDirName

◆ m_ntupleFileName

◆ m_ntupleTreeName

◆ m_numberOfContribPixelLayers

◆ m_numberOfPixelDeadSensors

◆ m_numberOfPixelHits

◆ m_numberOfPixelHoles

◆ m_numberOfSCTDeadSensors

◆ m_numberOfSCTHits

◆ m_numberOfSCTHoles

◆ m_numberOfTRTDeadStraws

◆ m_numberOfTRTHighThresholdHits

◆ m_numberOfTRTHits

◆ m_numberOfTRTHoles

◆ m_Rec_d0

◆ m_Rec_eta

◆ m_Rec_phi0

◆ m_Rec_qOverP

◆ m_Rec_z0