iFatras::HitCreatorTRT Class Reference

RIOOnTrack creation, starting from intersection on an active surface. More...

#include <HitCreatorTRT.h>

Inheritance diagram for iFatras::HitCreatorTRT:

Collaboration diagram for iFatras::HitCreatorTRT:

Public Member Functions
	HitCreatorTRT (const std::string &, const std::string &, const IInterface *)
	Constructor.
virtual	~HitCreatorTRT ()=default
	Destructor.
StatusCode	initialize ()
	AlgTool initailize method.
void	handle (const Incident &inc)
	handle for incident service
void	createSimHit (const ISF::ISFParticle &isp, const Trk::TrackParameters &, double) const
	Record the hit.
const ParametersROT *	createHit (const ISF::ISFParticle &, const Trk::TrackParameters &) const
	Return the cluster on Track – the PrepRawData is contained in this one.
const std::vector< ParametersROT > *	createHits (const ISF::ISFParticle &, const ParametersLayer &) const
	Return the cluster on Track – the PrepRawData is contained in this one.

Private Attributes
ServiceHandle< IIncidentSvc >	m_incidentSvc {this, "IncidentService", "IncidentSvc"}
TRTUncompressedHitCollection *	m_hitColl {}
	the sim hit collection
StringProperty	m_collectionName {this, "CollectionName", "TRTUncompressedHits"}
	name of the collection on storegate
ServiceHandle< IAtRndmGenSvc >	m_randomSvc {this, "RandomNumberService", "AtDSFMTGenSvc"}
	Pointer to the random number generator service.
StringProperty	m_randomEngineName {this, "RandomStreamName", "FatrasRnd"}
	Name of the random number stream.
CLHEP::HepRandomEngine *	m_randomEngine {}
	Random Engine.
StringProperty	m_trtIdHelperName {this, "TRT_IdHelperName", "TRT_ID"}
	where to find the SCT helper
const TRT_ID *	m_trtIdHelper {}
	TRT ID helper.
ToolHandle< ITRT_StrawStatusSummaryTool >	m_trtStatusSummaryTool {this, "StrawStatusSummaryTool", "TRT_StrawStatusSummaryTool"}
	Handle to TRT conditions service.
bool	m_useConditionsSvc {false}

Detailed Description

RIOOnTrack creation, starting from intersection on an active surface.

Author

Andreas.Salzburger -at- cern.ch

Definition at line 44 of file HitCreatorTRT.h.

Constructor & Destructor Documentation

HitCreatorTRT()

iFatras::HitCreatorTRT::HitCreatorTRT	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p )

Constructor.

Definition at line 36 of file HitCreatorTRT.cxx.

                                                             :
  base_class(t,n,p)
{
}

~HitCreatorTRT()

virtual iFatras::HitCreatorTRT::~HitCreatorTRT ( )

virtualdefault

Destructor.

Member Function Documentation

createHit()

const ParametersROT * iFatras::HitCreatorTRT::createHit ( const ISF::ISFParticle & , const Trk::TrackParameters &  ) const

inline

Return the cluster on Track – the PrepRawData is contained in this one.

Definition at line 64 of file HitCreatorTRT.h.

{ return 0; };

createHits()

const std::vector< ParametersROT > * iFatras::HitCreatorTRT::createHits ( const ISF::ISFParticle & , const ParametersLayer &  ) const

inline

Return the cluster on Track – the PrepRawData is contained in this one.

Definition at line 67 of file HitCreatorTRT.h.

{ return 0; }

createSimHit()

void iFatras::HitCreatorTRT::createSimHit	(	const ISF::ISFParticle &	isp,
		const Trk::TrackParameters &	pars,
		double	time ) const

Record the hit.

<

Todo

: fill information, decide which one is pre/post step point

Definition at line 92 of file HitCreatorTRT.cxx.

{
  // get surface and DetElement base
  const Trk::Surface &hitSurface = pars.associatedSurface();
  const Trk::TrkDetElementBase* detElementBase = hitSurface.associatedDetectorElement();
  // the detector element cast -> needed
  const InDetDD::TRT_BaseElement* hitTrtDetElement = dynamic_cast<const InDetDD::TRT_BaseElement*>((detElementBase));
  //if cast fails, return 0
  if ( !hitTrtDetElement ) return;
  // checks for detector element and surface
  if (hitTrtDetElement) {
    // get the real identifier
    Identifier hitId = hitSurface.associatedDetectorElementIdentifier();
    // check conditions of the straw
    if (m_useConditionsSvc &&  !m_trtStatusSummaryTool->get_status(hitId, Gaudi::Hive::currentContext()) ){
      ATH_MSG_VERBOSE("[ trthit ] Straw is not active ( says ConditionsTool). Ignore. ");
      return;
    }
    // create entry/exit point
    const Amg::Vector3D& gpos1 = pars.position();
    Amg::Vector3D gpos2 = gpos1 + pars.momentum().unit();
    // get the rotation
    const Amg::Transform3D& sTransform = hitSurface.transform();
    // create the hep points
    Amg::Vector3D lpoint1(sTransform.inverse() * gpos1);
    Amg::Vector3D lpoint2(sTransform.inverse() * gpos2);
 
    // new direction in 3D frame of cylinder
    Amg::Vector3D direc((lpoint2 - lpoint1).unit());
 
    // get line and circle constants
    double k = (direc.y())/(direc.x());
    double d = (lpoint2.x()*lpoint1.y() - lpoint1.x()*lpoint2.y())/(direc.x());
    double R = hitSurface.bounds().r();
 
    // and solve the qaudratic equation
    Trk::RealQuadraticEquation pquad(1 + k*k, 2*k*d, d*d - R*R);
    if (pquad.solutions != Trk::none){
      double t1 = (pquad.first  - lpoint1.x())/direc.x();
      double t2 = (pquad.second - lpoint1.x())/direc.x();
      // the solutions in the 3D frame of the cylinder : local coordination
      Amg::Vector3D sol1raw( (lpoint1 + t1 * direc) );
      Amg::Vector3D sol2raw( (lpoint1 + t2 * direc) );
 
      ATH_MSG_VERBOSE("[ trthit ] Straw center : " << hitSurface.center());
      ATH_MSG_VERBOSE("[ trthit ] Entry / exit : " << sol1raw << ", " << sol2raw );

      const double mass               = isp.mass();
      const double mom2               = pars.momentum().mag2();
      const float  kineticEnergy      = sqrt(mom2+mass*mass)-mass;
      const float  energyDepositInKeV = 0.; // only used for TRT hits directly created by photons
 
      int barrel_endcap = 0;
      int ispos = 0;
      switch ( m_trtIdHelper->barrel_ec(hitId) ) {
      case -2:  barrel_endcap = 1; ispos = 1; break;
      case -1:  barrel_endcap = 0; ispos = 0; break;
      case  1:  barrel_endcap = 0; ispos = 1; break;
      case  2:  barrel_endcap = 1; ispos = 0; break;
      default:
        ATH_MSG_WARNING("[ trthit ] Problem estimating barrel_endcap and ispos" ); return;
      }
 
      const int ringwheel = m_trtIdHelper->layer_or_wheel(hitId);
      const int phisector = m_trtIdHelper->phi_module(hitId);
      const int layer     = m_trtIdHelper->straw_layer(hitId);
      const int istraw    = m_trtIdHelper->straw(hitId);
 
      const TRTHitIdHelper* hitid_helper = TRTHitIdHelper::GetHelper();
      const int hitID = hitid_helper->buildHitId( barrel_endcap, ispos, ringwheel, phisector,layer,istraw);
 
      ATH_MSG_VERBOSE("[ trthit ] Add a TRTUncompressedHit hit to the collection. ");
      // fill into the HitCollection
      HepMcParticleLink partLink(HepMC::uniqueID(isp), 0,
                                 HepMcParticleLink::IS_POSITION,
                                 HepMcParticleLink::IS_ID);
      m_hitColl->Emplace( hitID,
                          partLink,
                          isp.pdgCode(),
                          (float) kineticEnergy, (float) energyDepositInKeV,
                          (float) sol1raw.x(), (float) sol1raw.y(), sol1raw.z(),
                          (float) sol2raw.x(), (float) sol2raw.y(), sol2raw.z(),
                          (float) time );
    }
 
  }
}

handle()

void iFatras::HitCreatorTRT::handle

(

const Incident &

inc

)

handle for incident service

Definition at line 72 of file HitCreatorTRT.cxx.

                                                       {
  // check the incident type
  if ( inc.type() == IncidentType::BeginEvent ){
    // check if the hit collection already contains:
    // (a)     if yes ... try to retrieve it
    if ( evtStore()->contains<TRTUncompressedHitCollection>(m_collectionName) ){
      if ( (evtStore()->retrieve(m_hitColl , m_collectionName)).isFailure() )
        ATH_MSG_ERROR( "[ --- ] Unable to retrieve TRTUncompressedHitCollection " << m_collectionName);
      // (a)     if no ... try to create it
    } else {
      m_hitColl = new TRTUncompressedHitCollection( m_collectionName);
      if ( (evtStore()->record(m_hitColl, m_collectionName, true)).isFailure() ) {
        ATH_MSG_ERROR( "[ --- ] Unable to record TRTUncompressedHitCollection " << m_collectionName);
        delete m_hitColl; m_hitColl=0;
      }
    }
  }
  return;
}

initialize()

StatusCode iFatras::HitCreatorTRT::initialize

(

)

AlgTool initailize method.

Definition at line 46 of file HitCreatorTRT.cxx.

{
  // Random number service
  ATH_CHECK ( m_randomSvc.retrieve() );
  //Get own engine with own seeds:
  m_randomEngine = m_randomSvc->GetEngine(m_randomEngineName);
  if (!m_randomEngine) {
    ATH_MSG_ERROR( "[ --- ] Could not get random engine '" << m_randomEngineName << "'" );
    return StatusCode::FAILURE;
  }
 
  // Tools & Services
  ATH_CHECK ( m_trtStatusSummaryTool.retrieve( DisableTool{ !m_useConditionsSvc } ) );
  // Get the TRT Identifier-helper:
  ATH_CHECK ( detStore()->retrieve(m_trtIdHelper, m_trtIdHelperName) );
 
  // Athena/Gaudi framework
  ATH_CHECK ( m_incidentSvc.retrieve() );
  // register to the incident service: BeginEvent for TrackCollection
  m_incidentSvc->addListener( this, IncidentType::BeginEvent );
 
  ATH_MSG_VERBOSE( "[ trthit ]  initialize() successful." );
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_collectionName

StringProperty iFatras::HitCreatorTRT::m_collectionName {this, "CollectionName", "TRTUncompressedHits"}

private

name of the collection on storegate

Definition at line 76 of file HitCreatorTRT.h.

{this, "CollectionName", "TRTUncompressedHits"};           

m_hitColl

TRTUncompressedHitCollection* iFatras::HitCreatorTRT::m_hitColl {}

private

the sim hit collection

Definition at line 75 of file HitCreatorTRT.h.

{};                  

m_incidentSvc

ServiceHandle<IIncidentSvc> iFatras::HitCreatorTRT::m_incidentSvc {this, "IncidentService", "IncidentSvc"}

private

Definition at line 72 of file HitCreatorTRT.h.

{this, "IncidentService", "IncidentSvc"};

m_randomEngine

CLHEP::HepRandomEngine* iFatras::HitCreatorTRT::m_randomEngine {}

private

Random Engine.

Definition at line 81 of file HitCreatorTRT.h.

{};                 

m_randomEngineName

StringProperty iFatras::HitCreatorTRT::m_randomEngineName {this, "RandomStreamName", "FatrasRnd"}

private

Name of the random number stream.

Definition at line 80 of file HitCreatorTRT.h.

{this, "RandomStreamName", "FatrasRnd"};             

m_randomSvc

ServiceHandle<IAtRndmGenSvc> iFatras::HitCreatorTRT::m_randomSvc {this, "RandomNumberService", "AtDSFMTGenSvc"}

private

Pointer to the random number generator service.

Random Svc

Definition at line 79 of file HitCreatorTRT.h.

{this, "RandomNumberService", "AtDSFMTGenSvc"};                    

m_trtIdHelper

const TRT_ID* iFatras::HitCreatorTRT::m_trtIdHelper {}

private

TRT ID helper.

Definition at line 84 of file HitCreatorTRT.h.

{};                  

m_trtIdHelperName

StringProperty iFatras::HitCreatorTRT::m_trtIdHelperName {this, "TRT_IdHelperName", "TRT_ID"}

private

where to find the SCT helper

Definition at line 83 of file HitCreatorTRT.h.

{this, "TRT_IdHelperName", "TRT_ID"};              

m_trtStatusSummaryTool

ToolHandle<ITRT_StrawStatusSummaryTool> iFatras::HitCreatorTRT::m_trtStatusSummaryTool {this, "StrawStatusSummaryTool", "TRT_StrawStatusSummaryTool"}

private

Handle to TRT conditions service.

Definition at line 86 of file HitCreatorTRT.h.

{this, "StrawStatusSummaryTool", "TRT_StrawStatusSummaryTool"};          

m_useConditionsSvc

bool iFatras::HitCreatorTRT::m_useConditionsSvc {false}

private

Definition at line 87 of file HitCreatorTRT.h.

{false};

---

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

• HitCreatorTRT.h
• HitCreatorTRT.cxx