Trk::PropResultRootWriterSvc Class Reference

#include <PropResultRootWriterSvc.h>

Inheritance diagram for Trk::PropResultRootWriterSvc:

Collaboration diagram for Trk::PropResultRootWriterSvc:

Public Member Functions
	PropResultRootWriterSvc (const std::string &name, ISvcLocator *svcloc)
	constructor More...
virtual	~PropResultRootWriterSvc ()
	destructor More...
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvInterface) override
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
template<typename T >
void	write (const T *initialPerigee, const T *fwdParameters=nullptr, double fwdtime=std::numeric_limits< float >::quiet_NaN(), const T *bkwParameters=nullptr, double bkwtime=std::numeric_limits< float >::quiet_NaN())
MsgStream &	msg () const
MsgStream &	msg (const MSG::Level lvl) const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()
	Retrieve interface ID. More...

Private Attributes
std::mutex	m_writeMutex
ServiceHandle< ITHistSvc >	m_thistSvc
TTree *	m_tree
Gaudi::Property< std::string >	m_ntupleDirName {this, "DirName", "/ExtrapolationStudies/", ""}
Gaudi::Property< std::string >	m_treeName {this, "TreeName", "ATLAS", ""}
int	m_eventNum = 0
float	m_start_d0 = 0.0F
float	m_start_z0 = 0.0F
float	m_start_phi = 0.0F
float	m_start_theta = 0.0F
float	m_start_qop = 0.0F
float	m_fwd_success = 0.0F
float	m_fwd_time = 0.0F
float	m_fwd_final_l0 = 0.0F
float	m_fwd_final_l1 = 0.0F
float	m_fwd_final_x = 0.0F
float	m_fwd_final_y = 0.0F
float	m_fwd_final_z = 0.0F
float	m_fwd_final_phi = 0.0F
float	m_fwd_final_theta = 0.0F
float	m_fwd_final_qop = 0.0F
float	m_fwd_final_sigma_l0 = 0.0F
float	m_fwd_final_sigma_l1 = 0.0F
float	m_fwd_final_sigma_phi = 0.0F
float	m_fwd_final_sigma_theta = 0.0F
float	m_fwd_final_sigma_qop = 0.0F
float	m_bkw_success = 0.0F
float	m_bkw_time = 0.0F
float	m_bkw_final_d0 = 0.0F
float	m_bkw_final_z0 = 0.0F
float	m_bkw_final_phi = 0.0F
float	m_bkw_final_theta = 0.0F
float	m_bkw_final_qop = 0.0F
float	m_bkw_final_sigma_d0 = 0.0F
float	m_bkw_final_sigma_z0 = 0.0F
float	m_bkw_final_sigma_phi = 0.0F
float	m_bkw_final_sigma_theta = 0.0F
float	m_bkw_final_sigma_qop = 0.0F

Detailed Description

Definition at line 22 of file PropResultRootWriterSvc.h.

Constructor & Destructor Documentation

PropResultRootWriterSvc()

Trk::PropResultRootWriterSvc::PropResultRootWriterSvc	(	const std::string &	name,
		ISvcLocator *	svcloc
	)

constructor

Definition at line 15 of file PropResultRootWriterSvc.cxx.

: AthService(name, svc),
  m_thistSvc("THistSvc", name),
  m_tree(nullptr) {
}

~PropResultRootWriterSvc()

Trk::PropResultRootWriterSvc::~PropResultRootWriterSvc ( )

virtualdefault

destructor

Member Function Documentation

finalize()

StatusCode Trk::PropResultRootWriterSvc::finalize ( )

overridevirtual

Definition at line 80 of file PropResultRootWriterSvc.cxx.

{
  return StatusCode::SUCCESS;
}

initialize()

StatusCode Trk::PropResultRootWriterSvc::initialize ( )

overridevirtual

Definition at line 23 of file PropResultRootWriterSvc.cxx.

{
  ATH_CHECK( AthService::initialize() );
 
  CHECK(m_thistSvc.retrieve());
  
  const std::string treeName = m_treeName;  
  m_tree = new TTree(treeName.data(), "MyStudies");
  
  std::string fullNtupleName = static_cast<std::string>(m_ntupleDirName) + m_treeName;
  CHECK(m_thistSvc->regTree(fullNtupleName, m_tree));
 
  if(m_tree == nullptr) {
    ATH_MSG_ERROR("Unable to create TTree");
    return StatusCode::FAILURE;
  }
 
  m_tree->Branch("event_nr"   , &m_eventNum   );
  
  m_tree->Branch("start_d0"   , &m_start_d0   );
  m_tree->Branch("start_z0"   , &m_start_z0   );
  m_tree->Branch("start_phi"  , &m_start_phi  );
  m_tree->Branch("start_theta", &m_start_theta);
  m_tree->Branch("start_qop"  , &m_start_qop  );
  
  m_tree->Branch("fwd_success"    , &m_fwd_success    );
  m_tree->Branch("fwd_time"       , &m_fwd_time       );
  m_tree->Branch("fwd_final_l0"   , &m_fwd_final_l0   );
  m_tree->Branch("fwd_final_l1"   , &m_fwd_final_l1   );
  m_tree->Branch("fwd_final_x"    , &m_fwd_final_x    );
  m_tree->Branch("fwd_final_y"    , &m_fwd_final_y    );
  m_tree->Branch("fwd_final_z"    , &m_fwd_final_z    );
  m_tree->Branch("fwd_final_phi"  , &m_fwd_final_phi  );
  m_tree->Branch("fwd_final_theta", &m_fwd_final_theta);
  m_tree->Branch("fwd_final_qop"  , &m_fwd_final_qop  );
  m_tree->Branch("fwd_final_sigma_l0"   , &m_fwd_final_sigma_l0   );
  m_tree->Branch("fwd_final_sigma_l1"   , &m_fwd_final_sigma_l1   );
  m_tree->Branch("fwd_final_sigma_phi"  , &m_fwd_final_sigma_phi  );
  m_tree->Branch("fwd_final_sigma_theta", &m_fwd_final_sigma_theta);
  m_tree->Branch("fwd_final_sigma_qop"  , &m_fwd_final_sigma_qop  );
  
  m_tree->Branch("bkw_success"    , &m_bkw_success    );
  m_tree->Branch("bkw_time"       , &m_bkw_time       );
  m_tree->Branch("bkw_final_d0"   , &m_bkw_final_d0   );
  m_tree->Branch("bkw_final_z0"   , &m_bkw_final_z0   );
  m_tree->Branch("bkw_final_phi"  , &m_bkw_final_phi  );
  m_tree->Branch("bkw_final_theta", &m_bkw_final_theta);
  m_tree->Branch("bkw_final_qop"  , &m_bkw_final_qop  );
  m_tree->Branch("bkw_final_sigma_d0"   , &m_bkw_final_sigma_d0   );
  m_tree->Branch("bkw_final_sigma_z0"   , &m_bkw_final_sigma_z0   );
  m_tree->Branch("bkw_final_sigma_phi"  , &m_bkw_final_sigma_phi  );
  m_tree->Branch("bkw_final_sigma_theta", &m_bkw_final_sigma_theta);
  m_tree->Branch("bkw_final_sigma_qop"  , &m_bkw_final_sigma_qop  );
  
  return StatusCode::SUCCESS;
}

interfaceID()

static const InterfaceID& Trk::PropResultRootWriterSvc::interfaceID ( )

inlinestatic

Retrieve interface ID.

Declaration of the interface ID ( interface id, major version, minor version)

Definition at line 30 of file PropResultRootWriterSvc.h.

                                            {
      static const InterfaceID IID_PropResultRootWriterSvc("PropResultRootWriterSvc", 1, 0);
      return IID_PropResultRootWriterSvc;    
    }

msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

MsgStream& AthCommonMsg< Service >::msg ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

queryInterface()

virtual StatusCode Trk::PropResultRootWriterSvc::queryInterface ( const InterfaceID &  riid, void **  ppvInterface  )

inlineoverridevirtual

Definition at line 36 of file PropResultRootWriterSvc.h.

                                                                                             {
      ATH_MSG_DEBUG("in queryInterface()");    
      if (PropResultRootWriterSvc::interfaceID().versionMatch(riid)) {
        ATH_MSG_DEBUG("matched PropResultRootWriterSvc");
        *ppvInterface=(PropResultRootWriterSvc*)this;
      } else {
        ATH_MSG_DEBUG("Default to Service interface");
        return AthService::queryInterface(riid,ppvInterface);
      }
      return StatusCode::SUCCESS;
    }

write()

template<typename T >

void Trk::PropResultRootWriterSvc::write	(	const T *	initialPerigee,
		const T *	fwdParameters = nullptr,
		double	fwdtime = std::numeric_limits<float>::quiet_NaN(),
		const T *	bkwParameters = nullptr,
		double	bkwtime = std::numeric_limits<float>::quiet_NaN()
	)

Definition at line 105 of file PropResultRootWriterSvc.h.

                                                                               {
  
    const auto& ctx = Gaudi::Hive::currentContext();
  
    m_eventNum = ctx.eventID().event_number();
    
    if (initialPerigee) {
      m_start_d0   = initialPerigee->parameters()[0];
      m_start_z0   = initialPerigee->parameters()[1];
      m_start_phi  = initialPerigee->parameters()[2];
      m_start_theta= initialPerigee->parameters()[3];
      m_start_qop  = initialPerigee->parameters()[4];
    } else {
      m_start_d0   = std::numeric_limits<float>::quiet_NaN();
      m_start_z0   = std::numeric_limits<float>::quiet_NaN();
      m_start_phi  = std::numeric_limits<float>::quiet_NaN();
      m_start_theta= std::numeric_limits<float>::quiet_NaN();
      m_start_qop  = std::numeric_limits<float>::quiet_NaN();
    }
 
    if (fwdParameters) {
      m_fwd_success     = 1;   
      m_fwd_time        = fwdtime;
      m_fwd_final_l0    = fwdParameters->parameters()[0];
      m_fwd_final_l1    = fwdParameters->parameters()[1];
      m_fwd_final_phi   = fwdParameters->parameters()[2];
      m_fwd_final_theta = fwdParameters->parameters()[3];
      m_fwd_final_qop   = fwdParameters->parameters()[4];
      m_fwd_final_x     = fwdParameters->position()[0];
      m_fwd_final_y     = fwdParameters->position()[1];
      m_fwd_final_z     = fwdParameters->position()[2];
      if (fwdParameters->covariance()) {
        m_fwd_final_sigma_l0    = (*fwdParameters->covariance())(0,0);
        m_fwd_final_sigma_l1    = (*fwdParameters->covariance())(1,1);
        m_fwd_final_sigma_phi   = (*fwdParameters->covariance())(2,2);
        m_fwd_final_sigma_theta = (*fwdParameters->covariance())(3,3);
        m_fwd_final_sigma_qop   = (*fwdParameters->covariance())(4,4);
      } else {        
        m_fwd_final_sigma_l0    = std::numeric_limits<float>::quiet_NaN();
        m_fwd_final_sigma_l1    = std::numeric_limits<float>::quiet_NaN();
        m_fwd_final_sigma_phi   = std::numeric_limits<float>::quiet_NaN();
        m_fwd_final_sigma_theta = std::numeric_limits<float>::quiet_NaN();
        m_fwd_final_sigma_qop   = std::numeric_limits<float>::quiet_NaN();
      }
    } else {
      m_fwd_success           = std::numeric_limits<float>::quiet_NaN();
      m_fwd_time              = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_l0          = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_l1          = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_phi         = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_theta       = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_qop         = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_x           = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_y           = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_z           = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_sigma_l0    = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_sigma_l1    = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_sigma_phi   = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_sigma_theta = std::numeric_limits<float>::quiet_NaN();
      m_fwd_final_sigma_qop   = std::numeric_limits<float>::quiet_NaN();
    }
    
    if (bkwParameters) {
      m_bkw_success     = 1;   
      m_bkw_time        = bkwtime;
      m_bkw_final_d0    = bkwParameters->parameters()[0];
      m_bkw_final_z0    = bkwParameters->parameters()[1];
      m_bkw_final_phi   = bkwParameters->parameters()[2];
      m_bkw_final_theta = bkwParameters->parameters()[3];
      m_bkw_final_qop   = bkwParameters->parameters()[4];
      if (bkwParameters->covariance()) {
        m_bkw_final_sigma_d0    = (*bkwParameters->covariance())(0,0);
        m_bkw_final_sigma_z0    = (*bkwParameters->covariance())(1,1);
        m_bkw_final_sigma_phi   = (*bkwParameters->covariance())(2,2);
        m_bkw_final_sigma_theta = (*bkwParameters->covariance())(3,3);
        m_bkw_final_sigma_qop   = (*bkwParameters->covariance())(4,4);
      } else {        
        m_bkw_final_sigma_d0    = std::numeric_limits<float>::quiet_NaN();
        m_bkw_final_sigma_z0    = std::numeric_limits<float>::quiet_NaN();
        m_bkw_final_sigma_phi   = std::numeric_limits<float>::quiet_NaN();
        m_bkw_final_sigma_theta = std::numeric_limits<float>::quiet_NaN();
        m_bkw_final_sigma_qop   = std::numeric_limits<float>::quiet_NaN();
      }
    } else {
      m_bkw_success            = std::numeric_limits<float>::quiet_NaN();
      m_bkw_time               = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_d0           = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_z0           = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_phi          = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_theta        = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_qop          = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_sigma_d0     = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_sigma_z0     = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_sigma_phi    = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_sigma_theta  = std::numeric_limits<float>::quiet_NaN();
      m_bkw_final_sigma_qop    = std::numeric_limits<float>::quiet_NaN();
    }
    m_tree->Fill();
  }

Member Data Documentation

m_bkw_final_d0

float Trk::PropResultRootWriterSvc::m_bkw_final_d0 = 0.0F

private

Definition at line 91 of file PropResultRootWriterSvc.h.

m_bkw_final_phi

float Trk::PropResultRootWriterSvc::m_bkw_final_phi = 0.0F

private

Definition at line 93 of file PropResultRootWriterSvc.h.

m_bkw_final_qop

float Trk::PropResultRootWriterSvc::m_bkw_final_qop = 0.0F

private

Definition at line 95 of file PropResultRootWriterSvc.h.

m_bkw_final_sigma_d0

float Trk::PropResultRootWriterSvc::m_bkw_final_sigma_d0 = 0.0F

private

Definition at line 96 of file PropResultRootWriterSvc.h.

m_bkw_final_sigma_phi

float Trk::PropResultRootWriterSvc::m_bkw_final_sigma_phi = 0.0F

private

Definition at line 98 of file PropResultRootWriterSvc.h.

m_bkw_final_sigma_qop

float Trk::PropResultRootWriterSvc::m_bkw_final_sigma_qop = 0.0F

private

Definition at line 100 of file PropResultRootWriterSvc.h.

m_bkw_final_sigma_theta

float Trk::PropResultRootWriterSvc::m_bkw_final_sigma_theta = 0.0F

private

Definition at line 99 of file PropResultRootWriterSvc.h.

m_bkw_final_sigma_z0

float Trk::PropResultRootWriterSvc::m_bkw_final_sigma_z0 = 0.0F

private

Definition at line 97 of file PropResultRootWriterSvc.h.

m_bkw_final_theta

float Trk::PropResultRootWriterSvc::m_bkw_final_theta = 0.0F

private

Definition at line 94 of file PropResultRootWriterSvc.h.

m_bkw_final_z0

float Trk::PropResultRootWriterSvc::m_bkw_final_z0 = 0.0F

private

Definition at line 92 of file PropResultRootWriterSvc.h.

m_bkw_success

float Trk::PropResultRootWriterSvc::m_bkw_success = 0.0F

private

Definition at line 89 of file PropResultRootWriterSvc.h.

m_bkw_time

float Trk::PropResultRootWriterSvc::m_bkw_time = 0.0F

private

Definition at line 90 of file PropResultRootWriterSvc.h.

m_eventNum

int Trk::PropResultRootWriterSvc::m_eventNum = 0

private

Definition at line 65 of file PropResultRootWriterSvc.h.

m_fwd_final_l0

float Trk::PropResultRootWriterSvc::m_fwd_final_l0 = 0.0F

private

Definition at line 75 of file PropResultRootWriterSvc.h.

m_fwd_final_l1

float Trk::PropResultRootWriterSvc::m_fwd_final_l1 = 0.0F

private

Definition at line 76 of file PropResultRootWriterSvc.h.

m_fwd_final_phi

float Trk::PropResultRootWriterSvc::m_fwd_final_phi = 0.0F

private

Definition at line 80 of file PropResultRootWriterSvc.h.

m_fwd_final_qop

float Trk::PropResultRootWriterSvc::m_fwd_final_qop = 0.0F

private

Definition at line 82 of file PropResultRootWriterSvc.h.

m_fwd_final_sigma_l0

float Trk::PropResultRootWriterSvc::m_fwd_final_sigma_l0 = 0.0F

private

Definition at line 83 of file PropResultRootWriterSvc.h.

m_fwd_final_sigma_l1

float Trk::PropResultRootWriterSvc::m_fwd_final_sigma_l1 = 0.0F

private

Definition at line 84 of file PropResultRootWriterSvc.h.

m_fwd_final_sigma_phi

float Trk::PropResultRootWriterSvc::m_fwd_final_sigma_phi = 0.0F

private

Definition at line 85 of file PropResultRootWriterSvc.h.

m_fwd_final_sigma_qop

float Trk::PropResultRootWriterSvc::m_fwd_final_sigma_qop = 0.0F

private

Definition at line 87 of file PropResultRootWriterSvc.h.

m_fwd_final_sigma_theta

float Trk::PropResultRootWriterSvc::m_fwd_final_sigma_theta = 0.0F

private

Definition at line 86 of file PropResultRootWriterSvc.h.

m_fwd_final_theta

float Trk::PropResultRootWriterSvc::m_fwd_final_theta = 0.0F

private

Definition at line 81 of file PropResultRootWriterSvc.h.

m_fwd_final_x

float Trk::PropResultRootWriterSvc::m_fwd_final_x = 0.0F

private

Definition at line 77 of file PropResultRootWriterSvc.h.

m_fwd_final_y

float Trk::PropResultRootWriterSvc::m_fwd_final_y = 0.0F

private

Definition at line 78 of file PropResultRootWriterSvc.h.

m_fwd_final_z

float Trk::PropResultRootWriterSvc::m_fwd_final_z = 0.0F

private

Definition at line 79 of file PropResultRootWriterSvc.h.

m_fwd_success

float Trk::PropResultRootWriterSvc::m_fwd_success = 0.0F

private

Definition at line 73 of file PropResultRootWriterSvc.h.

m_fwd_time

float Trk::PropResultRootWriterSvc::m_fwd_time = 0.0F

private

Definition at line 74 of file PropResultRootWriterSvc.h.

m_ntupleDirName

Gaudi::Property<std::string> Trk::PropResultRootWriterSvc::m_ntupleDirName {this, "DirName", "/ExtrapolationStudies/", ""}

private

Definition at line 62 of file PropResultRootWriterSvc.h.

m_start_d0

float Trk::PropResultRootWriterSvc::m_start_d0 = 0.0F

private

Definition at line 67 of file PropResultRootWriterSvc.h.

m_start_phi

float Trk::PropResultRootWriterSvc::m_start_phi = 0.0F

private

Definition at line 69 of file PropResultRootWriterSvc.h.

m_start_qop

float Trk::PropResultRootWriterSvc::m_start_qop = 0.0F

private

Definition at line 71 of file PropResultRootWriterSvc.h.

m_start_theta

float Trk::PropResultRootWriterSvc::m_start_theta = 0.0F

private

Definition at line 70 of file PropResultRootWriterSvc.h.

m_start_z0

float Trk::PropResultRootWriterSvc::m_start_z0 = 0.0F

private

Definition at line 68 of file PropResultRootWriterSvc.h.

m_thistSvc

ServiceHandle<ITHistSvc> Trk::PropResultRootWriterSvc::m_thistSvc

private

Definition at line 60 of file PropResultRootWriterSvc.h.

m_tree

TTree* Trk::PropResultRootWriterSvc::m_tree

private

Definition at line 61 of file PropResultRootWriterSvc.h.

m_treeName

Gaudi::Property<std::string> Trk::PropResultRootWriterSvc::m_treeName {this, "TreeName", "ATLAS", ""}

private

Definition at line 63 of file PropResultRootWriterSvc.h.

m_writeMutex

std::mutex Trk::PropResultRootWriterSvc::m_writeMutex

private

Definition at line 57 of file PropResultRootWriterSvc.h.

---

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

• PropResultRootWriterSvc.h
• PropResultRootWriterSvc.cxx