#include <DetailedIDNtupleTool.h>

Inheritance diagram for InDet::DetailedIDNtupleTool:

Collaboration diagram for InDet::DetailedIDNtupleTool:

Public Member Functions
	DetailedIDNtupleTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~DetailedIDNtupleTool ()
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual void	dumpTrack (int itrk, const Trk::AlignTrack *alignTrack) override
	fills track information to ntuple
virtual void	storeHitmap () override
	stores hitmap for writing to ntuple
virtual void	fillHitmap () override
	fills ntuple with hit information
virtual void	fillSummary () override
	fills ntuple with event and track summary information
virtual void	showStatistics () override
	write statistics out to log file
virtual void	setNtuple (TFile *ntuple) override
	sets the output stream for the logfile
virtual StatusCode	fillNtuple () override
	writes trees and histograms to ntuple
virtual void	fillChi2VChamberShift (const int, const int, int , double , double , double , double , double , double *)
	used for numerical derivatives to write chi2 vs.
virtual void	setLogStream (std::ostream *os)
	returns pointer to TFile containing 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 ()

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.

Protected Attributes
std::ostream *	m_logStream = nullptr
	logfile output stream

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
bool	retrieveTruthInfo (const Trk::AlignTrack *trk)
const Trk::TrackParameters *	perigeeParameter (const Trk::AlignTrack *track) const
void	initializeNtuple ()
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ToolHandle< Trk::ITrackSummaryTool >	m_trackSumTool
	Pointer to track summary tool.
ToolHandle< Trk::IAlignModuleTool >	m_alignModuleTool
	Pointer to AlignmModuleTool.
PublicToolHandle< Trk::ITruthToTrack >	m_truthToTrack {this,"TruthToTrack","Trk::TruthToTrack/InDetTruthToTrack",""}
	the truth to track Tool
PublicToolHandle< Trk::IExtrapolator >	m_extrapolator {this,"Extrapolator","Trk::Extrapolator/AtlasExtrapolator",""}
	track extrapolator
bool	m_storeNormalRefittedOnly
	only store the AlignTrack which is normally refitted
bool	m_storeConstrainedOnly
	onlystore the AlignTrack which is normally refitted
bool	m_storeTruth
	retrieve the truth information
double	m_matchProbability
	the probabililty cut in the truth matching
TFile *	m_file
TTree *	m_tree
std::string	m_filename
	name of ntuple file
std::string	m_filepath
	path to ntuple file
std::string	m_trackCollection
std::string	m_tracksTruthName
int	m_runNumber = 0
int	m_evtNumber = 0
double	m_original_xvtx = 0.0
double	m_original_yvtx = 0.0
double	m_original_zvtx = 0.0
double	m_original_d0 = 0.0
double	m_original_z0 = 0.0
double	m_original_phi0 = 0.0
double	m_original_theta = 0.0
double	m_original_qoverp = 0.0
double	m_original_pt = 0.0
double	m_original_eta = 0.0
double	m_original_chi2 = 0.0
int	m_original_ndof = 0
double	m_original_chi2prob = 0.0
double	m_original_err_d0 = 0.0
double	m_original_err_z0 = 0.0
double	m_original_err_phi0 = 0.0
double	m_original_err_theta = 0.0
double	m_original_err_qoverp = 0.0
double	m_truth_d0 = 0.0
double	m_truth_z0 = 0.0
double	m_truth_phi0 = 0.0
double	m_truth_theta = 0.0
double	m_truth_eta = 0.0
double	m_truth_pt = 0.0
double	m_truth_qoverp = 0.0
double	m_truth_qoverpt = 0.0
double	m_truth_charge = 0.0
double	m_truth_prod_x = 0.0
double	m_truth_prod_y = 0.0
double	m_truth_prod_z = 0.0
double	m_original_toRef_d0 = 0.0
double	m_original_toRef_z0 = 0.0
double	m_original_toRef_phi0 = 0.0
double	m_original_toRef_theta = 0.0
double	m_original_toRef_qoverp = 0.0
double	m_toRef_d0 = 0.0
double	m_toRef_z0 = 0.0
double	m_toRef_phi0 = 0.0
double	m_toRef_theta = 0.0
double	m_toRef_qoverp = 0.0
double	m_err_d0 = 0.0
double	m_err_z0 = 0.0
double	m_err_phi0 = 0.0
double	m_err_theta = 0.0
double	m_err_qoverp = 0.0
double	m_xvtx = 0.0
double	m_yvtx = 0.0
double	m_zvtx = 0.0
double	m_d0 = 0.0
double	m_z0 = 0.0
double	m_phi0 = 0.0
double	m_theta = 0.0
double	m_qoverp = 0.0
double	m_pt = 0.0
double	m_eta = 0.0
double	m_chi2 = 0.0
int	m_ndof = 0
double	m_chi2prob = 0.0
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 36 of file DetailedIDNtupleTool.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

◆ DetailedIDNtupleTool()

InDet::DetailedIDNtupleTool::DetailedIDNtupleTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Definition at line 38 of file DetailedIDNtupleTool.cxx.

    : AthAlgTool(type, name, parent)
    , m_trackSumTool("Trk::TrackSummaryTool/TrackSummaryTool", this)
    , m_alignModuleTool{}
    , m_storeNormalRefittedOnly{false}
    , m_storeConstrainedOnly{true}
    , m_storeTruth{true}
    , m_matchProbability{0.8}
    , m_file{}
    , m_tree{}
    , m_trackCollection{"Tracks"}
    , m_tracksTruthName{"TrackTruthCollection"} {
    declareInterface<Trk::IFillNtupleTool>(this);
 
    declareProperty("FileName", m_filename = "IDAlign.root");
    declareProperty("FilePath", m_filepath = "./");
    declareProperty("TrackCollection", m_trackCollection);
    declareProperty("TracksTruthName", m_tracksTruthName);
    declareProperty("TrackSummaryTool", m_trackSumTool);
    declareProperty("AlignModuleTool", m_alignModuleTool);
    declareProperty("MatchProbability", m_matchProbability);
    declareProperty("StoreTruth", m_storeTruth);
    declareProperty("MatchProbability", m_matchProbability);
    declareProperty("StoreNormalRefittedOnly", m_storeNormalRefittedOnly);
    declareProperty("StoreConstrainedOnly", m_storeConstrainedOnly);
  }

◆ ~DetailedIDNtupleTool()

InDet::DetailedIDNtupleTool::~DetailedIDNtupleTool ( )

virtual

Definition at line 66 of file DetailedIDNtupleTool.cxx.

66 {

67 }

Member Function Documentation

◆ declareGaudiProperty()

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareProperty()

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

95{ return m_detStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ dumpTrack()

void InDet::DetailedIDNtupleTool::dumpTrack	(	int	itrk,
		const Trk::AlignTrack *	alignTrack )

overridevirtual

fills track information to ntuple

if match truth fail, return directly

Implements Trk::IFillNtupleTool.

Definition at line 106 of file DetailedIDNtupleTool.cxx.

                                                                             {
    ATH_MSG_DEBUG("In dumpTrack()");
    Trk::AlignTrack::AlignTrackType type = alignTrack->type();
    const EventContext& ctx = Gaudi::Hive::currentContext();
    // if hope to dump all the tracks, should set:
    // m_storeNormalRefittedOnly = false, m_storeConstrainedOnly=false
    if ((m_storeNormalRefittedOnly && (type != Trk::AlignTrack::NormalRefitted)) ||
        (m_storeConstrainedOnly && (type != Trk::AlignTrack::BeamspotConstrained))) {
      return;
    }
 
    // get run and event numbers
    ATH_MSG_DEBUG("Retrieving event info.");
    const xAOD::EventInfo* eventInfo;
    if (evtStore()->retrieve(eventInfo).isFailure()) {
      ATH_MSG_ERROR("Could not retrieve event info.");
    } else {
      m_runNumber = eventInfo->runNumber();
      m_evtNumber = eventInfo->eventNumber();
    }
 
    // initialize variables
    constexpr double invalidParameterValue{-999.};
    m_d0 = invalidParameterValue;
    m_z0 = invalidParameterValue;
    m_phi0 = invalidParameterValue;
    m_theta = invalidParameterValue;
    m_qoverp = invalidParameterValue;
    m_pt = invalidParameterValue;
    m_eta = invalidParameterValue;
    constexpr double invalidChiSq{-1e12};
    m_chi2 = invalidChiSq;
    m_chi2prob = invalidChiSq;
    constexpr int invalidDegreesOfFreedom{-999};
 
    m_ndof = invalidDegreesOfFreedom;
    m_xvtx = invalidParameterValue;
    m_yvtx = invalidParameterValue;
    m_zvtx = invalidParameterValue;
 
    m_err_d0 = invalidParameterValue;
    m_err_z0 = invalidParameterValue;
    m_err_phi0 = invalidParameterValue;
    m_err_theta = invalidParameterValue;
    m_err_qoverp = invalidParameterValue;
    m_toRef_d0 = invalidParameterValue;
    m_toRef_z0 = invalidParameterValue;
    m_toRef_phi0 = invalidParameterValue;
    m_toRef_theta = invalidParameterValue;
    m_toRef_qoverp = invalidParameterValue;
 
 
    m_original_d0 = invalidParameterValue;
    m_original_z0 = invalidParameterValue;
    m_original_phi0 = invalidParameterValue;
    m_original_theta = invalidParameterValue;
    m_original_qoverp = invalidParameterValue;
    m_original_pt =invalidParameterValue;
    m_original_eta = -999.;
    m_original_chi2 = invalidChiSq;
    m_original_chi2prob = invalidChiSq;
    m_original_ndof = invalidDegreesOfFreedom;
    m_original_xvtx = invalidParameterValue;
    m_original_yvtx = invalidParameterValue;
    m_original_zvtx = invalidParameterValue;
 
    m_original_err_d0 = invalidParameterValue;
    m_original_err_z0 = invalidParameterValue;
    m_original_err_phi0 = invalidParameterValue;
    m_original_err_theta = invalidParameterValue;
    m_original_err_qoverp = invalidParameterValue;;
    m_original_toRef_d0 = invalidParameterValue;
    m_original_toRef_z0 = invalidParameterValue;
    m_original_toRef_phi0 = invalidParameterValue;
    m_original_toRef_theta =invalidParameterValue;
    m_original_toRef_qoverp = invalidParameterValue;;
 
    m_truth_d0 = invalidParameterValue;
    m_truth_z0 = invalidParameterValue;
    m_truth_phi0 = invalidParameterValue;
    m_truth_theta = invalidParameterValue;
    m_truth_qoverp = invalidParameterValue;
    m_truth_pt = invalidParameterValue;
    m_truth_eta = invalidParameterValue;;
    m_truth_prod_x = invalidParameterValue;
    m_truth_prod_y = invalidParameterValue;
    m_truth_prod_z = invalidParameterValue;;
 
 
 
    const Trk::Perigee* aMeasPer = (alignTrack->perigeeParameters());
 
    if (not aMeasPer) {
      ATH_MSG_ERROR("Could not get Trk::Perigee of the alignTrack");
      return;
    } else {
      m_d0 = aMeasPer->parameters()[Trk::d0];
      m_z0 = aMeasPer->parameters()[Trk::z0];
      m_phi0 = aMeasPer->parameters()[Trk::phi0];
      m_theta = aMeasPer->parameters()[Trk::theta];
      m_qoverp = aMeasPer->parameters()[Trk::qOverP];
 
      const AmgSymMatrix(5) * locCov = aMeasPer->covariance();
      m_err_d0 = Amg::error(*locCov, Trk::d0);
      m_err_z0 = Amg::error(*locCov, Trk::z0);
      m_err_phi0 = Amg::error(*locCov, Trk::phi0);
      m_err_theta = Amg::error(*locCov, Trk::theta);
      m_err_qoverp = Amg::error(*locCov, Trk::qOverP);
 
      m_pt =
        std::sqrt((aMeasPer->momentum().x()) * (aMeasPer->momentum().x()) + (aMeasPer->momentum().y()) *
                  (aMeasPer->momentum().y()));
      m_eta = aMeasPer->eta();
 
      m_xvtx = aMeasPer->position().x();
      m_yvtx = aMeasPer->position().y();
      m_zvtx = aMeasPer->position().z();
 
 
 
      // get fit quality and chi2 probability of track
      const Trk::FitQuality* fitQual = alignTrack->fitQuality();
      if (not fitQual) ATH_MSG_ERROR("No fit quality assigned to the track");
      else {
        if (fitQual->chiSquared() > 0. && fitQual->numberDoF() > 0) {
          m_chi2 = fitQual->chiSquared();
          m_ndof = fitQual->numberDoF();
          m_chi2prob = TMath::Prob(m_chi2, m_ndof);
 
          ATH_MSG_DEBUG("  - chi2             : " << m_chi2);
          ATH_MSG_DEBUG("  - ndof             : " << m_ndof);
          ATH_MSG_DEBUG("  - chi2/ndof        : " << m_chi2 / (double) m_ndof);
          ATH_MSG_DEBUG("  - chi2 propability : " << m_chi2prob);
        }
      }
    }
 
    const Trk::MeasurementBase* measbase = *(alignTrack->measurementsOnTrack()->begin());
    Amg::Vector3D refPoint = measbase->associatedSurface().center();
    const Trk::VertexOnTrack* vot = dynamic_cast<const Trk::VertexOnTrack*> (measbase);
    if (!vot) {
      ATH_MSG_ERROR(" Seems the pseudo-measuremnt in the alignTrack not exist!");
      ATH_MSG_ERROR(" this pseudo-measurement has been rejected as outlier in the refitting!");
      return;
    } else {
      ATH_MSG_DEBUG(" VoT get from the alignTrack:" << *vot);
    }
 
    ATH_MSG_DEBUG(" the pseudo-measurement position: " << refPoint);
    const Trk::Track* originalTrack = alignTrack->originalTrack();
 
    const Trk::Perigee* originalMeasPer = originalTrack->perigeeParameters();
    if (!originalMeasPer) {
      ATH_MSG_ERROR("No original track!");
    } else {
      m_original_d0 = originalMeasPer->parameters()[Trk::d0];
      m_original_z0 = originalMeasPer->parameters()[Trk::z0];
      m_original_phi0 = originalMeasPer->parameters()[Trk::phi0];
      m_original_theta = originalMeasPer->parameters()[Trk::theta];
      m_original_qoverp = originalMeasPer->parameters()[Trk::qOverP];
 
      const AmgSymMatrix(5) * locError = originalMeasPer->covariance();
      m_original_err_d0 = Amg::error(*locError, Trk::d0);
      m_original_err_z0 = Amg::error(*locError, Trk::z0);
      m_original_err_phi0 = Amg::error(*locError, Trk::phi0);
      m_original_err_theta = Amg::error(*locError, Trk::theta);
      m_original_err_qoverp = Amg::error(*locError, Trk::qOverP);
 
      m_original_pt = (1.0 / m_original_qoverp) * std::sin(m_original_theta);
      m_original_eta = -std::log(std::tan(0.5 * m_original_theta));
 
      m_original_xvtx = originalMeasPer->position().x();
      m_original_yvtx = originalMeasPer->position().y();
      m_original_zvtx = originalMeasPer->position().z();
 
      // get fit quality and chi2 probability of original Track
      const Trk::FitQuality* originalFitQual = originalTrack->fitQuality();
      if (not originalFitQual) ATH_MSG_ERROR("No fit quality assigned to the track");
      else {
        if (originalFitQual->chiSquared() > 0. && originalFitQual->numberDoF() > 0) {
          m_original_chi2 = originalFitQual->chiSquared();
          m_original_ndof = originalFitQual->numberDoF();
          m_original_chi2prob = TMath::Prob(m_original_chi2, m_original_ndof);
 
          ATH_MSG_DEBUG("  - original chi2             : " << m_original_chi2);
          ATH_MSG_DEBUG("  - original ndof             : " << m_original_ndof);
          ATH_MSG_DEBUG("  - original chi2/ndof        : " << m_original_chi2 / (double) m_original_ndof);
          ATH_MSG_DEBUG("  - original chi2 propability : " << m_original_chi2prob);
        }
      }
 
 
      const Trk::PerigeeSurface persf(refPoint);
      const Trk::Perigee* originalPerigeeAtRef =
        dynamic_cast<const Trk::Perigee*>(m_extrapolator->extrapolateTrack(ctx, *originalTrack, persf).release());
      if (!originalPerigeeAtRef) {
        const Trk::Perigee* originalTrackPerigee = originalTrack->perigeeParameters();
        if (originalTrackPerigee && ((originalTrackPerigee->associatedSurface())) == persf) {
          ATH_MSG_DEBUG("Perigee of Track is already expressed to given vertex, a copy is returned.");
          originalPerigeeAtRef = originalTrackPerigee->clone();
        } else ATH_MSG_DEBUG("Extrapolation to Perigee failed, NULL pointer is returned.");
      }
 
      if (originalPerigeeAtRef) {
        std::unique_ptr<const Trk::Perigee > originalMeasPerAtRef(originalPerigeeAtRef);
        m_original_toRef_d0 = originalMeasPerAtRef->parameters()[Trk::d0];
        m_original_toRef_z0 = originalMeasPerAtRef->parameters()[Trk::z0];
        m_original_toRef_phi0 = originalMeasPerAtRef->parameters()[Trk::phi0];
        m_original_toRef_theta = originalMeasPerAtRef->parameters()[Trk::theta];
        m_original_toRef_qoverp = originalMeasPerAtRef->parameters()[Trk::qOverP];
      }
 
      const Trk::Perigee* PerigeeAtRef =
        dynamic_cast<const Trk::Perigee*>(m_extrapolator->extrapolateTrack(ctx, *alignTrack, persf).release());
      if (!PerigeeAtRef) {
        const Trk::Perigee* alignTrackPerigee = alignTrack->perigeeParameters();
        if (alignTrackPerigee && ((alignTrackPerigee->associatedSurface())) == persf) {
          ATH_MSG_DEBUG("Perigee of AlignTrack is already expressed to given vertex, a copy is returned.");
          PerigeeAtRef = alignTrackPerigee->clone();
        } else ATH_MSG_DEBUG("Extrapolation to Perigee failed, NULL pointer is returned.");
      }
 
      //post-eigen, can simply use the TrackParameters * returned by m_extrapolator->extrapolate?
      if (PerigeeAtRef) {
        std::unique_ptr<const Trk::Perigee > MeasPerAtRef((PerigeeAtRef));
        m_toRef_d0 = MeasPerAtRef->parameters()[Trk::d0];
        m_toRef_z0 = MeasPerAtRef->parameters()[Trk::z0];
        m_toRef_phi0 = MeasPerAtRef->parameters()[Trk::phi0];
        m_toRef_theta = MeasPerAtRef->parameters()[Trk::theta];
        m_toRef_qoverp = MeasPerAtRef->parameters()[Trk::qOverP];
      }
    }
    if (m_storeTruth) retrieveTruthInfo(alignTrack);
 
 
    m_file->cd();
    m_tree->Fill();
    ATH_MSG_DEBUG("tree filled");
 
    return;
  }

◆ 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

◆ fillChi2VChamberShift()

virtual void Trk::IFillNtupleTool::fillChi2VChamberShift	(	const int	,
		const int	,
		int *	,
		double **	,
		double **	,
		double **	,
		double **	,
		double **	,
		double **	)

inlinevirtualinherited

used for numerical derivatives to write chi2 vs.

chamber positions

Definition at line 59 of file IFillNtupleTool.h.

61 {}

◆ fillHitmap()

void InDet::DetailedIDNtupleTool::fillHitmap ( )

overridevirtual

fills ntuple with hit information

Implements Trk::IFillNtupleTool.

Definition at line 464 of file DetailedIDNtupleTool.cxx.

464 {

465 }

◆ fillNtuple()

StatusCode InDet::DetailedIDNtupleTool::fillNtuple ( )

overridevirtual

writes trees and histograms to ntuple

Implements Trk::IFillNtupleTool.

Definition at line 91 of file DetailedIDNtupleTool.cxx.

                                              {
    if (m_file && m_file->IsOpen()) {
      m_file->cd();
      if (m_tree) m_tree->Write();
    }
    return StatusCode::SUCCESS;
  }

◆ fillSummary()

void InDet::DetailedIDNtupleTool::fillSummary ( )

overridevirtual

fills ntuple with event and track summary information

Implements Trk::IFillNtupleTool.

Definition at line 467 of file DetailedIDNtupleTool.cxx.

467 {

468 }

◆ finalize()

StatusCode InDet::DetailedIDNtupleTool::finalize ( )

overridevirtual

Implements Trk::IFillNtupleTool.

Definition at line 100 of file DetailedIDNtupleTool.cxx.

                                            {
    ATH_MSG_DEBUG("finalize() of DetailedIDNtupleTool");
    return StatusCode::SUCCESS;
  }

◆ initialize()

StatusCode InDet::DetailedIDNtupleTool::initialize ( )

overridevirtual

Implements Trk::IFillNtupleTool.

Definition at line 70 of file DetailedIDNtupleTool.cxx.

                                              {
    ATH_MSG_DEBUG("initialize() of DetailedIDNtupleTool");
    ATH_MSG_DEBUG("creating file with name " << m_filepath << m_filename);
 
    if (!m_tree) initializeNtuple();
    ATH_CHECK(m_truthToTrack.retrieve());
 
    /* Retrieve extrapolator tool */
    ATH_CHECK(m_extrapolator.retrieve());
 
    // get AlignModuleTool
    if (m_alignModuleTool.empty()) {
      ATH_MSG_FATAL("m_alignModuleTool is empty in  DetailedIDNtupleTool::initialize()");
      return StatusCode::FAILURE;
    } else {
      ATH_CHECK(m_alignModuleTool.retrieve());
    }
    return StatusCode::SUCCESS;
  }

◆ initializeNtuple()

void InDet::DetailedIDNtupleTool::initializeNtuple ( )

private

Definition at line 475 of file DetailedIDNtupleTool.cxx.

                                              {
    //m_file = new TFile((m_filepath+m_filename).c_str(), "RECREATE");
    m_file->cd();
    m_tree = new TTree("IDAlign", "Inner Detector Alignment Ntuple");
 
    m_tree->Branch("run", &m_runNumber, "run/I");
    m_tree->Branch("evt", &m_evtNumber, "evt/I");
 
    m_tree->Branch("xvtx", &m_xvtx, "xvtx/D");
    m_tree->Branch("yvtx", &m_yvtx, "yvtx/D");
    m_tree->Branch("zvtx", &m_zvtx, "zvtx/D");
    m_tree->Branch("d0", &m_d0, "d0/D");
    m_tree->Branch("z0", &m_z0, "z0/D");
    m_tree->Branch("phi0", &m_phi0, "phi0/D");
    m_tree->Branch("theta", &m_theta, "theta/D");
    m_tree->Branch("qoverp", &m_qoverp, "qoverp/D");
    m_tree->Branch("pt", &m_pt, "pt/D");
    m_tree->Branch("eta", &m_eta, "eta/D");
    m_tree->Branch("chi2", &m_chi2, "chi2/D");
    m_tree->Branch("ndof", &m_ndof, "ndof/I");
    m_tree->Branch("chi2prob", &m_chi2prob, "chi2prob/D");
    m_tree->Branch("err_d0", &m_err_d0, "err_d0/D");
    m_tree->Branch("err_z0", &m_err_z0, "err_z0/D");
    m_tree->Branch("err_phi0", &m_err_phi0, "err_phi0/D");
    m_tree->Branch("err_theta", &m_err_theta, "err_theta/D");
    m_tree->Branch("err_qoverp", &m_err_qoverp, "err_qoverp/D");
 
    if (m_storeTruth) {
      m_tree->Branch("m_truth_prod_x", &m_truth_prod_x, "truth_prod_x/D");
      m_tree->Branch("m_truth_prod_y", &m_truth_prod_y, "truth_prod_y/D");
      m_tree->Branch("m_truth_prod_z", &m_truth_prod_z, "truth_prod_z/D");
      m_tree->Branch("truth_d0", &m_truth_d0, "truth_d0/D");
      m_tree->Branch("truth_z0", &m_truth_z0, "truth_z0/D");
      m_tree->Branch("truth_phi0", &m_truth_phi0, "truth_phi0/D");
      m_tree->Branch("truth_theta", &m_truth_theta, "truth_theta/D");
      m_tree->Branch("truth_qoverp", &m_truth_qoverp, "truth_qoverp/D");
      m_tree->Branch("truth_pt", &m_truth_pt, "truth_pt/D");
      m_tree->Branch("truth_eta", &m_truth_eta, "truth_eta/D");
    }
 
    m_tree->Branch("original_xvtx", &m_original_xvtx, "original_xvtx/D");
    m_tree->Branch("original_yvtx", &m_original_yvtx, "original_yvtx/D");
    m_tree->Branch("original_zvtx", &m_original_zvtx, "original_zvtx/D");
    m_tree->Branch("original_d0", &m_original_d0, "original_d0/D");
    m_tree->Branch("original_z0", &m_original_z0, "original_z0/D");
    m_tree->Branch("original_phi0", &m_original_phi0, "original_phi0/D");
    m_tree->Branch("original_theta", &m_original_theta, "original_theta/D");
    m_tree->Branch("original_qoverp", &m_original_qoverp, "original_qoverp/D");
    m_tree->Branch("original_pt", &m_original_pt, "original_pt/D");
    m_tree->Branch("original_eta", &m_original_eta, "original_eta/D");
    m_tree->Branch("original_chi2", &m_original_chi2, "original_chi2/D");
    m_tree->Branch("original_ndof", &m_original_ndof, "original_ndof/I");
    m_tree->Branch("original_chi2prob", &m_original_chi2prob, "original_chi2prob/D");
 
    m_tree->Branch("original_err_d0", &m_original_err_d0, "original_err_d0/D");
    m_tree->Branch("original_err_z0", &m_original_err_z0, "original_err_z0/D");
    m_tree->Branch("original_err_phi0", &m_original_err_phi0, "original_err_phi0/D");
    m_tree->Branch("original_err_theta", &m_original_err_theta, "original_err_theta/D");
    m_tree->Branch("original_err_qoverp", &m_original_err_qoverp, "original_err_qoverp/D");
 
    m_tree->Branch("toRef_d0", &m_toRef_d0, "toRef_d0/D");
    m_tree->Branch("toRef_z0", &m_toRef_z0, "toRef_z0/D");
    m_tree->Branch("toRef_phi0", &m_toRef_phi0, "toRef_phi0/D");
    m_tree->Branch("toRef_theta", &m_toRef_theta, "toRef_theta/D");
    m_tree->Branch("toRef_qoverp", &m_toRef_qoverp, "toRef_qoverp/D");
    m_tree->Branch("original_toRef_d0", &m_original_toRef_d0, "original_toRef_d0/D");
    m_tree->Branch("original_toRef_z0", &m_original_toRef_z0, "original_toRef_z0/D");
    m_tree->Branch("original_toRef_phi0", &m_original_toRef_phi0, "original_toRef_phi0/D");
    m_tree->Branch("original_toRef_theta", &m_original_toRef_theta, "original_toRef_theta/D");
    m_tree->Branch("original_toRef_qoverp", &m_original_toRef_qoverp, "original_toRef_qoverp/D");
 
    return;
  }

◆ 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 & IFillNtupleTool::interfaceID ( )

inlinestaticinherited

Definition at line 80 of file IFillNtupleTool.h.

81 { return IID_TRKALIGNINTERFACES_IFillNtupleTool; }

IID_TRKALIGNINTERFACES_IFillNtupleTool

static const InterfaceID IID_TRKALIGNINTERFACES_IFillNtupleTool("IFillNtupleTool", 1, 0)

◆ 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.

◆ perigeeParameter()

const Trk::TrackParameters * InDet::DetailedIDNtupleTool::perigeeParameter ( const Trk::AlignTrack * track ) const

private

Definition at line 445 of file DetailedIDNtupleTool.cxx.

                                                                                                        {
    const Trk::TrackParameters* aMeasPer {};
 
    for (const auto i:*(alignTrack->trackStateOnSurfaces())) {
      if (i->type(Trk::TrackStateOnSurface::Perigee)) {
        aMeasPer = (i->trackParameters());
        break;
      }
    }
    if (not aMeasPer) {
      ATH_MSG_ERROR("Could not get Trk::MeasuredPerigee of the alignTrack");
    }
    return aMeasPer;
  }

◆ 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();
  }

◆ retrieveTruthInfo()

bool InDet::DetailedIDNtupleTool::retrieveTruthInfo ( const Trk::AlignTrack * trk )

private

Definition at line 350 of file DetailedIDNtupleTool.cxx.

                                                                              {
    // although we select tracks using the TrackSelectionTool, we still need to get a complete TrackCollection
    // from StoreGate for use in the track-truth map, otherwise the track-truth matching is screwed up
 
    const TrackCollection* RecCollection = nullptr;
 
    if (evtStore()->retrieve(RecCollection, m_trackCollection).isFailure()) {
      ATH_MSG_WARNING("Track collection \"" << m_trackCollection << "\" not found.");
      return false;
    }
    if (RecCollection) {
      ATH_MSG_DEBUG(
        "Retrieved " << m_trackCollection << " with size " << RecCollection->size() <<
      " reconstructed tracks from storegate");
    } else {
      ATH_MSG_WARNING("RecCollection is null pointer in DetailedIDNtupleTool::retrieveTruthInfo");
      return false;
    }
 
    // get TrackTruthCollection
    const TrackTruthCollection* TruthMap = nullptr;
    if (StatusCode::SUCCESS != evtStore()->retrieve(TruthMap, m_tracksTruthName)) {
      ATH_MSG_DEBUG("Cannot find " << m_tracksTruthName);
      return false;
    }
 
    ATH_MSG_DEBUG(
      "Track Truth Collection with name " << m_tracksTruthName << " with size " << TruthMap->size() <<
    " found in StoreGate");
 
    bool flag = false;
    // get fit quality and chi2 probability of track
    const Trk::Perigee* startPerigee = alignTrack->perigeeParameters();
    const Trk::Perigee* measPer = startPerigee;
 
    if (measPer == nullptr) {
      ATH_MSG_DEBUG("No measured perigee parameters assigned to the track");
      return false;
    }
 
    ElementLink<TrackCollection> tracklink;
    tracklink.setElement(const_cast<Trk::Track*>(alignTrack->originalTrack()));
    tracklink.setStorableObject(*RecCollection);
    const ElementLink<TrackCollection> tracklink2 = tracklink;
 
    TrackTruthCollection::const_iterator found = TruthMap->find(tracklink2);
    if ((found != TruthMap->end()) && (found->second.probability() > m_matchProbability)) {
      TrackTruth trtruth = found->second;
      const HepMcParticleLink& HMPL = trtruth.particleLink();
 
      if (HMPL.isValid()) {
#ifdef HEPMC3
        HepMC::ConstGenParticlePtr genparptr = HMPL.scptr();
#else
        const HepMC::GenParticle* genparptr = HMPL.cptr();
#endif
 
        if (genparptr) {
          if (genparptr->production_vertex()) {
            if (genparptr->pdg_id() == 0) {
              ATH_MSG_INFO("PDG ID is zero in DetailedIDNtupleTool::retrieveTruthInfo");
            } else {
              const Trk::TrackParameters* generatedTrackPerigee = m_truthToTrack->makePerigeeParameters(genparptr);
              if (!generatedTrackPerigee) ATH_MSG_WARNING("Unable to extrapolate genparticle to perigee!");
              else {
                flag = true;
                m_truth_qoverpt = generatedTrackPerigee->parameters()[Trk::qOverP] / std::sin(
                  generatedTrackPerigee->parameters()[Trk::theta]);
                m_truth_qoverp = generatedTrackPerigee->parameters()[Trk::qOverP];
                m_truth_phi0 = generatedTrackPerigee->parameters()[Trk::phi0];
                m_truth_d0 = generatedTrackPerigee->parameters()[Trk::d0];
                m_truth_z0 = generatedTrackPerigee->parameters()[Trk::z0];
                m_truth_theta = generatedTrackPerigee->parameters()[Trk::theta];
                m_truth_eta = generatedTrackPerigee->eta();
                m_truth_prod_x = genparptr->production_vertex()->position().x();
                m_truth_prod_y = genparptr->production_vertex()->position().y();
                m_truth_prod_z = genparptr->production_vertex()->position().z();
 
                delete  generatedTrackPerigee;
                m_truth_pt = 1. / std::abs(m_truth_qoverpt);
                m_truth_charge = 1;
                if (m_truth_qoverpt < 0) m_truth_charge = -1;
                if (m_truth_phi0 < 0) m_truth_phi0 += 2 * M_PI;
                ATH_MSG_DEBUG("Found matched truth track with phi, PT = " << m_truth_phi0 << ", " << m_truth_pt);
              }
            }
          }
        }
      }
    }
 
    return flag;
  }

◆ setLogStream()

virtual void Trk::IFillNtupleTool::setLogStream ( std::ostream * os )

inlinevirtualinherited

returns pointer to TFile containing ntuple

sets the output stream for the logfile

Definition at line 67 of file IFillNtupleTool.h.

67{ m_logStream = os; }

Trk::IFillNtupleTool::m_logStream

std::ostream * m_logStream

logfile output stream

Definition IFillNtupleTool.h:76

ReadFromCoolCompare.os

int os

Definition ReadFromCoolCompare.py:226

◆ setNtuple()

virtual void InDet::DetailedIDNtupleTool::setNtuple ( TFile * ntuple )

inlineoverridevirtual

sets the output stream for the logfile

Implements Trk::IFillNtupleTool.

Definition at line 54 of file DetailedIDNtupleTool.h.

54{ m_file=ntuple; }

◆ showStatistics()

void InDet::DetailedIDNtupleTool::showStatistics ( )

overridevirtual

write statistics out to log file

Implements Trk::IFillNtupleTool.

Definition at line 471 of file DetailedIDNtupleTool.cxx.

471 {

472 }

◆ storeHitmap()

void InDet::DetailedIDNtupleTool::storeHitmap ( )

overridevirtual

stores hitmap for writing to ntuple

Implements Trk::IFillNtupleTool.

Definition at line 460 of file DetailedIDNtupleTool.cxx.

460 {

461 }

◆ sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

overridevirtualinherited

Handle START transition.

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

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_alignModuleTool

ToolHandle<Trk::IAlignModuleTool> InDet::DetailedIDNtupleTool::m_alignModuleTool

private

Pointer to AlignmModuleTool.

Definition at line 63 of file DetailedIDNtupleTool.h.

◆ m_chi2

double InDet::DetailedIDNtupleTool::m_chi2 = 0.0

private

Definition at line 155 of file DetailedIDNtupleTool.h.

◆ m_chi2prob

double InDet::DetailedIDNtupleTool::m_chi2prob = 0.0

private

Definition at line 157 of file DetailedIDNtupleTool.h.

◆ m_d0

double InDet::DetailedIDNtupleTool::m_d0 = 0.0

private

Definition at line 148 of file DetailedIDNtupleTool.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_err_d0

double InDet::DetailedIDNtupleTool::m_err_d0 = 0.0

private

Definition at line 139 of file DetailedIDNtupleTool.h.

◆ m_err_phi0

double InDet::DetailedIDNtupleTool::m_err_phi0 = 0.0

private

Definition at line 141 of file DetailedIDNtupleTool.h.

◆ m_err_qoverp

double InDet::DetailedIDNtupleTool::m_err_qoverp = 0.0

private

Definition at line 143 of file DetailedIDNtupleTool.h.

◆ m_err_theta

double InDet::DetailedIDNtupleTool::m_err_theta = 0.0

private

Definition at line 142 of file DetailedIDNtupleTool.h.

◆ m_err_z0

double InDet::DetailedIDNtupleTool::m_err_z0 = 0.0

private

Definition at line 140 of file DetailedIDNtupleTool.h.

◆ m_eta

double InDet::DetailedIDNtupleTool::m_eta = 0.0

private

Definition at line 154 of file DetailedIDNtupleTool.h.

◆ m_evtNumber

int InDet::DetailedIDNtupleTool::m_evtNumber = 0

private

Definition at line 85 of file DetailedIDNtupleTool.h.

◆ 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_extrapolator

PublicToolHandle<Trk::IExtrapolator> InDet::DetailedIDNtupleTool::m_extrapolator {this,"Extrapolator","Trk::Extrapolator/AtlasExtrapolator",""}

private

track extrapolator

Definition at line 66 of file DetailedIDNtupleTool.h.

67{this,"Extrapolator","Trk::Extrapolator/AtlasExtrapolator",""};

◆ m_file

TFile* InDet::DetailedIDNtupleTool::m_file

private

Definition at line 74 of file DetailedIDNtupleTool.h.

◆ m_filename

std::string InDet::DetailedIDNtupleTool::m_filename

private

name of ntuple file

Definition at line 77 of file DetailedIDNtupleTool.h.

◆ m_filepath

std::string InDet::DetailedIDNtupleTool::m_filepath

private

path to ntuple file

Definition at line 78 of file DetailedIDNtupleTool.h.

◆ m_logStream

std::ostream* Trk::IFillNtupleTool::m_logStream = nullptr

protectedinherited

logfile output stream

Definition at line 76 of file IFillNtupleTool.h.

◆ m_matchProbability

double InDet::DetailedIDNtupleTool::m_matchProbability

private

the probabililty cut in the truth matching

Definition at line 72 of file DetailedIDNtupleTool.h.

◆ m_ndof

int InDet::DetailedIDNtupleTool::m_ndof = 0

private

Definition at line 156 of file DetailedIDNtupleTool.h.

◆ m_original_chi2

double InDet::DetailedIDNtupleTool::m_original_chi2 = 0.0

private

Definition at line 100 of file DetailedIDNtupleTool.h.

◆ m_original_chi2prob

double InDet::DetailedIDNtupleTool::m_original_chi2prob = 0.0

private

Definition at line 102 of file DetailedIDNtupleTool.h.

◆ m_original_d0

double InDet::DetailedIDNtupleTool::m_original_d0 = 0.0

private

Definition at line 92 of file DetailedIDNtupleTool.h.

◆ m_original_err_d0

double InDet::DetailedIDNtupleTool::m_original_err_d0 = 0.0

private

Definition at line 104 of file DetailedIDNtupleTool.h.

◆ m_original_err_phi0

double InDet::DetailedIDNtupleTool::m_original_err_phi0 = 0.0

private

Definition at line 106 of file DetailedIDNtupleTool.h.

◆ m_original_err_qoverp

double InDet::DetailedIDNtupleTool::m_original_err_qoverp = 0.0

private

Definition at line 108 of file DetailedIDNtupleTool.h.

◆ m_original_err_theta

double InDet::DetailedIDNtupleTool::m_original_err_theta = 0.0

private

Definition at line 107 of file DetailedIDNtupleTool.h.

◆ m_original_err_z0

double InDet::DetailedIDNtupleTool::m_original_err_z0 = 0.0

private

Definition at line 105 of file DetailedIDNtupleTool.h.

◆ m_original_eta

double InDet::DetailedIDNtupleTool::m_original_eta = 0.0

private

Definition at line 99 of file DetailedIDNtupleTool.h.

◆ m_original_ndof

int InDet::DetailedIDNtupleTool::m_original_ndof = 0

private

Definition at line 101 of file DetailedIDNtupleTool.h.

◆ m_original_phi0

double InDet::DetailedIDNtupleTool::m_original_phi0 = 0.0

private

Definition at line 94 of file DetailedIDNtupleTool.h.

◆ m_original_pt

double InDet::DetailedIDNtupleTool::m_original_pt = 0.0

private

Definition at line 97 of file DetailedIDNtupleTool.h.

◆ m_original_qoverp

double InDet::DetailedIDNtupleTool::m_original_qoverp = 0.0

private

Definition at line 96 of file DetailedIDNtupleTool.h.

◆ m_original_theta

double InDet::DetailedIDNtupleTool::m_original_theta = 0.0

private

Definition at line 95 of file DetailedIDNtupleTool.h.

◆ m_original_toRef_d0

double InDet::DetailedIDNtupleTool::m_original_toRef_d0 = 0.0

private

Definition at line 125 of file DetailedIDNtupleTool.h.

◆ m_original_toRef_phi0

double InDet::DetailedIDNtupleTool::m_original_toRef_phi0 = 0.0

private

Definition at line 127 of file DetailedIDNtupleTool.h.

◆ m_original_toRef_qoverp

double InDet::DetailedIDNtupleTool::m_original_toRef_qoverp = 0.0

private

Definition at line 129 of file DetailedIDNtupleTool.h.

◆ m_original_toRef_theta

double InDet::DetailedIDNtupleTool::m_original_toRef_theta = 0.0

private

Definition at line 128 of file DetailedIDNtupleTool.h.

◆ m_original_toRef_z0

double InDet::DetailedIDNtupleTool::m_original_toRef_z0 = 0.0

private

Definition at line 126 of file DetailedIDNtupleTool.h.

◆ m_original_xvtx

double InDet::DetailedIDNtupleTool::m_original_xvtx = 0.0

private

Definition at line 88 of file DetailedIDNtupleTool.h.

◆ m_original_yvtx

double InDet::DetailedIDNtupleTool::m_original_yvtx = 0.0

private

Definition at line 89 of file DetailedIDNtupleTool.h.

◆ m_original_z0

double InDet::DetailedIDNtupleTool::m_original_z0 = 0.0

private

Definition at line 93 of file DetailedIDNtupleTool.h.

◆ m_original_zvtx

double InDet::DetailedIDNtupleTool::m_original_zvtx = 0.0

private

Definition at line 90 of file DetailedIDNtupleTool.h.

◆ m_phi0

double InDet::DetailedIDNtupleTool::m_phi0 = 0.0

private

Definition at line 150 of file DetailedIDNtupleTool.h.

◆ m_pt

double InDet::DetailedIDNtupleTool::m_pt = 0.0

private

Definition at line 153 of file DetailedIDNtupleTool.h.

◆ m_qoverp

double InDet::DetailedIDNtupleTool::m_qoverp = 0.0

private

Definition at line 152 of file DetailedIDNtupleTool.h.

◆ m_runNumber

int InDet::DetailedIDNtupleTool::m_runNumber = 0

private

Definition at line 84 of file DetailedIDNtupleTool.h.

◆ m_storeConstrainedOnly

bool InDet::DetailedIDNtupleTool::m_storeConstrainedOnly

private

onlystore the AlignTrack which is normally refitted

Definition at line 70 of file DetailedIDNtupleTool.h.

◆ m_storeNormalRefittedOnly

bool InDet::DetailedIDNtupleTool::m_storeNormalRefittedOnly

private

only store the AlignTrack which is normally refitted

Definition at line 69 of file DetailedIDNtupleTool.h.

◆ m_storeTruth

bool InDet::DetailedIDNtupleTool::m_storeTruth

private

retrieve the truth information

Definition at line 71 of file DetailedIDNtupleTool.h.

◆ m_theta

double InDet::DetailedIDNtupleTool::m_theta = 0.0

private

Definition at line 151 of file DetailedIDNtupleTool.h.

◆ m_toRef_d0

double InDet::DetailedIDNtupleTool::m_toRef_d0 = 0.0

private

Definition at line 132 of file DetailedIDNtupleTool.h.

◆ m_toRef_phi0

double InDet::DetailedIDNtupleTool::m_toRef_phi0 = 0.0

private

Definition at line 134 of file DetailedIDNtupleTool.h.

◆ m_toRef_qoverp

double InDet::DetailedIDNtupleTool::m_toRef_qoverp = 0.0

private

Definition at line 136 of file DetailedIDNtupleTool.h.

◆ m_toRef_theta

double InDet::DetailedIDNtupleTool::m_toRef_theta = 0.0

private

Definition at line 135 of file DetailedIDNtupleTool.h.

◆ m_toRef_z0

double InDet::DetailedIDNtupleTool::m_toRef_z0 = 0.0

private

Definition at line 133 of file DetailedIDNtupleTool.h.

◆ m_trackCollection

std::string InDet::DetailedIDNtupleTool::m_trackCollection

private

Definition at line 80 of file DetailedIDNtupleTool.h.

◆ m_tracksTruthName

std::string InDet::DetailedIDNtupleTool::m_tracksTruthName

private

Definition at line 81 of file DetailedIDNtupleTool.h.

◆ m_trackSumTool

ToolHandle<Trk::ITrackSummaryTool> InDet::DetailedIDNtupleTool::m_trackSumTool

private

Pointer to track summary tool.

Definition at line 62 of file DetailedIDNtupleTool.h.

◆ m_tree

TTree* InDet::DetailedIDNtupleTool::m_tree

private

Definition at line 75 of file DetailedIDNtupleTool.h.

◆ m_truth_charge

double InDet::DetailedIDNtupleTool::m_truth_charge = 0.0

private

Definition at line 119 of file DetailedIDNtupleTool.h.

◆ m_truth_d0

double InDet::DetailedIDNtupleTool::m_truth_d0 = 0.0

private

Definition at line 111 of file DetailedIDNtupleTool.h.

◆ m_truth_eta

double InDet::DetailedIDNtupleTool::m_truth_eta = 0.0

private

Definition at line 115 of file DetailedIDNtupleTool.h.

◆ m_truth_phi0

double InDet::DetailedIDNtupleTool::m_truth_phi0 = 0.0

private

Definition at line 113 of file DetailedIDNtupleTool.h.

◆ m_truth_prod_x

double InDet::DetailedIDNtupleTool::m_truth_prod_x = 0.0

private

Definition at line 120 of file DetailedIDNtupleTool.h.

◆ m_truth_prod_y

double InDet::DetailedIDNtupleTool::m_truth_prod_y = 0.0

private

Definition at line 121 of file DetailedIDNtupleTool.h.

◆ m_truth_prod_z

double InDet::DetailedIDNtupleTool::m_truth_prod_z = 0.0

private

Definition at line 122 of file DetailedIDNtupleTool.h.

◆ m_truth_pt

double InDet::DetailedIDNtupleTool::m_truth_pt = 0.0

private

Definition at line 116 of file DetailedIDNtupleTool.h.

◆ m_truth_qoverp

double InDet::DetailedIDNtupleTool::m_truth_qoverp = 0.0

private

Definition at line 117 of file DetailedIDNtupleTool.h.

◆ m_truth_qoverpt

double InDet::DetailedIDNtupleTool::m_truth_qoverpt = 0.0

private

Definition at line 118 of file DetailedIDNtupleTool.h.

◆ m_truth_theta

double InDet::DetailedIDNtupleTool::m_truth_theta = 0.0

private

Definition at line 114 of file DetailedIDNtupleTool.h.

◆ m_truth_z0

double InDet::DetailedIDNtupleTool::m_truth_z0 = 0.0

private

Definition at line 112 of file DetailedIDNtupleTool.h.

◆ m_truthToTrack

PublicToolHandle<Trk::ITruthToTrack> InDet::DetailedIDNtupleTool::m_truthToTrack {this,"TruthToTrack","Trk::TruthToTrack/InDetTruthToTrack",""}

private

the truth to track Tool

Definition at line 64 of file DetailedIDNtupleTool.h.

65{this,"TruthToTrack","Trk::TruthToTrack/InDetTruthToTrack",""};

◆ 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.

◆ m_xvtx

double InDet::DetailedIDNtupleTool::m_xvtx = 0.0

private

Definition at line 145 of file DetailedIDNtupleTool.h.

◆ m_yvtx

double InDet::DetailedIDNtupleTool::m_yvtx = 0.0

private

Definition at line 146 of file DetailedIDNtupleTool.h.

◆ m_z0

double InDet::DetailedIDNtupleTool::m_z0 = 0.0

private

Definition at line 149 of file DetailedIDNtupleTool.h.

◆ m_zvtx

double InDet::DetailedIDNtupleTool::m_zvtx = 0.0

private

Definition at line 147 of file DetailedIDNtupleTool.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ DetailedIDNtupleTool()

◆ ~DetailedIDNtupleTool()

Member Function Documentation

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ dumpTrack()

◆ evtStore()

◆ extraDeps_update_handler()

◆ fillChi2VChamberShift()

◆ fillHitmap()

◆ fillNtuple()

◆ fillSummary()

◆ finalize()

◆ initialize()

◆ initializeNtuple()

◆ inputHandles()

◆ interfaceID()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ perigeeParameter()

◆ renounce()

◆ renounceArray()

◆ retrieveTruthInfo()

◆ setLogStream()

◆ setNtuple()

◆ showStatistics()

◆ storeHitmap()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_alignModuleTool

◆ m_chi2

◆ m_chi2prob

◆ m_d0

◆ m_detStore

◆ m_err_d0

◆ m_err_phi0

◆ m_err_qoverp

◆ m_err_theta

◆ m_err_z0

◆ m_eta

◆ m_evtNumber

◆ m_evtStore

◆ m_extrapolator

◆ m_file

◆ m_filename

◆ m_filepath

◆ m_logStream

◆ m_matchProbability

◆ m_ndof

◆ m_original_chi2

◆ m_original_chi2prob

◆ m_original_d0

◆ m_original_err_d0

◆ m_original_err_phi0

◆ m_original_err_qoverp

◆ m_original_err_theta

◆ m_original_err_z0

◆ m_original_eta

◆ m_original_ndof

◆ m_original_phi0

◆ m_original_pt

◆ m_original_qoverp

◆ m_original_theta

◆ m_original_toRef_d0

◆ m_original_toRef_phi0