InDet::TRT_TrackExtensionToolCosmics Class Reference

#include <TRT_TrackExtensionToolCosmics.h>

Inheritance diagram for InDet::TRT_TrackExtensionToolCosmics:

Collaboration diagram for InDet::TRT_TrackExtensionToolCosmics:

Classes
class	EventData

Public Member Functions
	TRT_TrackExtensionToolCosmics (const std::string &, const std::string &, const IInterface *)
virtual	~TRT_TrackExtensionToolCosmics ()
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual std::vector< const Trk::MeasurementBase * > &	extendTrack (const EventContext &ctx, const Trk::Track &, InDet::ITRT_TrackExtensionTool::IEventData &virt_event_data, InDet::TRT_DetElementLink_xk::TRT_DetElemUsedMap &used) const override
virtual std::vector< const Trk::MeasurementBase * > &	extendTrack (const EventContext &ctx, const Trk::TrackParameters *, InDet::ITRT_TrackExtensionTool::IEventData &virt_event_data, InDet::TRT_DetElementLink_xk::TRT_DetElemUsedMap &used) const override
virtual Trk::TrackSegment *	findSegment (const EventContext &ctx, const Trk::TrackParameters *, InDet::ITRT_TrackExtensionTool::IEventData &virt_event_data, InDet::TRT_DetElementLink_xk::TRT_DetElemUsedMap &used) const override
virtual Trk::Track *	newTrack (const EventContext &ctx, const Trk::Track &, InDet::ITRT_TrackExtensionTool::IEventData &virt_event_data, InDet::TRT_DetElementLink_xk::TRT_DetElemUsedMap &used) const override
virtual std::unique_ptr< InDet::ITRT_TrackExtensionTool::IEventData >	newEvent (const EventContext &ctx) const override
virtual MsgStream &	dump (MsgStream &out) const override
virtual std::ostream &	dump (std::ostream &out) const override
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	analyze_tpars (const std::vector< const Trk::TrackParameters * > *tpars, InDet::TRT_TrackExtensionToolCosmics::EventData &event_data) const
Trk::Surface *	findBoundarySurface (const Trk::TrackParameters &par, Trk::PropDirection dir, InDet::TRT_TrackExtensionToolCosmics::EventData &event_data) const
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.

Static Protected Member Functions
static MsgStream &	dumpConditions (MsgStream &out)
static MsgStream &	dumpEvent (MsgStream &out)
static Amg::Vector3D	intersect (const Trk::Surface *surf, const Trk::Perigee *per)

Protected Attributes
ToolHandle< Trk::IExtrapolator >	m_extrapolator {this, "Extrapolator", "Trk::Extrapolator/InDetExtrapolator"}
ToolHandle< Trk::IRIO_OnTrackCreator >	m_riontrackN {this, "RIOonTrackToolNoDr", "InDet::TRT_DriftCircleOnTrackNoDriftTimeTool/TRT_DriftCircleOnTrackNoDriftTimeTool"}
int	m_outputlevel {}
int	m_nprint {}
DoubleProperty	m_roadwidth
DoubleProperty	m_roadwidth_locz
BooleanProperty	m_searchNeighbour
SG::ReadHandleKey< TRT_DriftCircleContainer >	m_trtname {this,"TRT_ClustersContainer","TRT_DriftCircles","RHK to retrieve TRT_DriftCircles"}
AtlasDetectorID *	m_idHelper {}
const TRT_ID *	m_trtid {}

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
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

Author

Chris.nosp@m.tian.nosp@m..Schm.nosp@m.itt@.nosp@m.cern..nosp@m.ch

Definition at line 49 of file TRT_TrackExtensionToolCosmics.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

TRT_TrackExtensionToolCosmics()

InDet::TRT_TrackExtensionToolCosmics::TRT_TrackExtensionToolCosmics	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p )

Definition at line 31 of file TRT_TrackExtensionToolCosmics.cxx.

  : AthAlgTool(t,n,p)
{
  declareInterface<ITRT_TrackExtensionTool>(this);
}

~TRT_TrackExtensionToolCosmics()

InDet::TRT_TrackExtensionToolCosmics::~TRT_TrackExtensionToolCosmics ( )

virtualdefault

Member Function Documentation

analyze_tpars()

void InDet::TRT_TrackExtensionToolCosmics::analyze_tpars ( const std::vector< const Trk::TrackParameters * > * tpars, InDet::TRT_TrackExtensionToolCosmics::EventData & event_data ) const

protected

Definition at line 230 of file TRT_TrackExtensionToolCosmics.cxx.

{
  msg(MSG::DEBUG)<<"Number of tpars: "<<tpars->size()<<endmsg;
 
  double lastz=-99999;
  std::vector< const Trk::TrackParameters* >::const_iterator parameterIter = tpars->begin();
  for ( ; parameterIter != tpars->end(); ++parameterIter) {
    msg(MSG::DEBUG)<< "par pos: " << (**parameterIter).position() <<endmsg;
      
    if ( (*parameterIter)->associatedSurface().associatedDetectorElementIdentifier()==0 ) {
      msg(MSG::DEBUG)<<"No DE identifier!!!"<<endmsg;
      continue; 
    }
 
    const Identifier& DCId = (*parameterIter)->associatedSurface().associatedDetectorElementIdentifier();
    if ( m_trtid->is_trt( DCId ) ) {
      
      IdentifierHash detElements[3];
      detElements[1]=(*parameterIter)->associatedSurface().associatedDetectorElement()->identifyHash();
      if(m_searchNeighbour){
    //fill entry 0 and 2 with the neighbours (different in phi but identical in layer)
    int bec=m_trtid->barrel_ec(DCId);
    int phi_mod=m_trtid->phi_module(DCId);
    int layer_or_wheel=m_trtid->layer_or_wheel(DCId);
    int slay=m_trtid->straw_layer(DCId);
    
    phi_mod-=1;
    if(phi_mod<0) phi_mod=31;
    
    Identifier temp=m_trtid->straw_id(bec,phi_mod,layer_or_wheel,slay,0);
    detElements[0]=m_trtid->straw_layer_hash(temp);
    
    phi_mod+=2;
    if(phi_mod>31) phi_mod=0;
    temp=m_trtid->straw_id(bec,phi_mod,layer_or_wheel,slay,0);
    detElements[2]=m_trtid->straw_layer_hash(temp);
      }
      double maxdist=m_roadwidth;
      const InDet::TRT_DriftCircle *circ=nullptr;
      
      for(int i=-1;i<2;i++) {
    if(m_searchNeighbour || i==0){      
      
      //check if this PRD exists
      // get the driftCircleCollection belonging to this id
      const InDet::TRT_DriftCircleCollection *container = event_data.m_trtcontainer->indexFindPtr(detElements[i+1]);
      
      if(container==nullptr) {
        msg(MSG::DEBUG)<<"for the current detectorElement no DriftCircleContainer seems to exist: "<<m_trtid->show_to_string(m_trtid->layer_id(detElements[i+1]))<<endmsg;
        continue;
      }
      
      msg(MSG::DEBUG)<< "There are "  << container->size() << " entries in the TRT_DriftCircleCollection "<<m_trtid->show_to_string(m_trtid->layer_id(detElements[i+1])) <<endmsg;
      
      //take the closest one in case it satisfies some default cuts
      InDet::TRT_DriftCircleCollection::const_iterator driftCircleIterator = container->begin();
      for (; driftCircleIterator != container->end(); ++driftCircleIterator) {
 
        //get the associated surface of the driftcircle
        const Trk::Surface &dc_surface=(*driftCircleIterator)->detectorElement()->surface((*driftCircleIterator)->identify());
 
        //get the local position of the track prediction in the frame of the driftcircle
      std::optional<Amg::Vector2D> lpos=dc_surface.globalToLocal((*parameterIter)->position());
 
        double distance=m_roadwidth+1;
        if(lpos){
          distance = std::abs(lpos->x());
          msg(MSG::DEBUG)<<"Hit "<<m_trtid->show_to_string((*driftCircleIterator)->identify())<<" has a distance of "<<distance<<endmsg;
 
          double dist_locz=std::abs(lpos->y());
          if(distance<m_roadwidth+1){
        if(!dc_surface.insideBounds(*lpos,m_roadwidth,m_roadwidth_locz)){
          msg(MSG::DEBUG)<<"Hit not inside surface bounds! "<<distance<<" , "<<dist_locz<<endmsg;
          msg(MSG::DEBUG)<<"\trejecting hit"<<endmsg;
          distance=m_roadwidth+1;
        }
          }
 
        }
 
        if(distance<maxdist){
          maxdist=distance;
          circ=(*driftCircleIterator);
        }
      }
    }
      }
      msg(MSG::DEBUG)<<"Maximal distance: "<<maxdist<<endmsg;
      if(circ){
    msg(MSG::DEBUG)<<"Found Driftcircle! Adding it to list ..."<<m_trtid->show_to_string(circ->identify())<<endmsg;
        if (lastz<-9999) lastz=(**parameterIter).position().z();
        if (std::abs(lastz-(**parameterIter).position().z())>500.) return;
        lastz=(**parameterIter).position().z();
        const Trk::StraightLineSurface *slsurf=dynamic_cast<const Trk::StraightLineSurface *>(&circ->detectorElement()->surface(circ->identify())); if(!slsurf) continue;
        Trk::AtaStraightLine atasl((**parameterIter).position(),(**parameterIter).parameters()[Trk::phi],(**parameterIter).parameters()[Trk::theta],(**parameterIter).parameters()[Trk::qOverP],*slsurf);
        const Trk::MeasurementBase *newmeas=m_riontrackN->correct(*circ,atasl,Gaudi::Hive::currentContext()); 
    event_data.m_measurement.push_back(newmeas);
 
      }   
    }
  }
}

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.

{ return m_detStore; }

dump() [1/2]

MsgStream & InDet::TRT_TrackExtensionToolCosmics::dump ( MsgStream & out ) const

overridevirtual

Implements InDet::ITRT_TrackExtensionTool.

Definition at line 101 of file TRT_TrackExtensionToolCosmics.cxx.

{
  out<<std::endl;
  if(m_nprint)  return dumpEvent(out);
  return dumpConditions(out);
}

dump() [2/2]

std::ostream & InDet::TRT_TrackExtensionToolCosmics::dump ( std::ostream & out ) const

overridevirtual

Implements InDet::ITRT_TrackExtensionTool.

Definition at line 132 of file TRT_TrackExtensionToolCosmics.cxx.

{
  return out;
}

dumpConditions()

MsgStream & InDet::TRT_TrackExtensionToolCosmics::dumpConditions ( MsgStream & out )

staticprotected

Definition at line 113 of file TRT_TrackExtensionToolCosmics.cxx.

{
 
  return out;
}

dumpEvent()

MsgStream & InDet::TRT_TrackExtensionToolCosmics::dumpEvent ( MsgStream & out )

staticprotected

Definition at line 123 of file TRT_TrackExtensionToolCosmics.cxx.

{
  return out;
}

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.

{ return m_evtStore; }

extendTrack() [1/2]

std::vector< const Trk::MeasurementBase * > & InDet::TRT_TrackExtensionToolCosmics::extendTrack ( const EventContext & ctx, const Trk::Track & Tr, InDet::ITRT_TrackExtensionTool::IEventData & virt_event_data, InDet::TRT_DetElementLink_xk::TRT_DetElemUsedMap & used ) const

overridevirtual

Implements InDet::ITRT_TrackExtensionTool.

Definition at line 192 of file TRT_TrackExtensionToolCosmics.cxx.

{ 
  InDet::TRT_TrackExtensionToolCosmics::EventData &
     event_data=InDet::TRT_TrackExtensionToolCosmics::EventData::getPrivateEventData(virt_event_data);
 
  event_data.m_measurement.clear();
 
  if(not event_data.m_trtcontainer) return event_data.m_measurement;
 
  const Trk::TrackStates* 
    tsos = Tr.trackStateOnSurfaces();
 
  const Trk::TrackParameters* 
    par = (*(tsos->rbegin()))->trackParameters(); if(!par ) return event_data.m_measurement;
  const Trk::TrackParameters* 
    parb = (*(tsos->begin()))->trackParameters();
 
  
  if(parb && par!=parb) {
 
    const Amg::Vector3D& g1 = par ->position();
    const Amg::Vector3D& g2 = parb->position();
    if((g2.x()*g2.x()+g2.y()*g2.y()) > (g1.x()*g1.x()+g1.y()*g1.y())) par=parb;
  }
 
  if(Tr.perigeeParameters()) {
    return extendTrack(ctx, Tr.perigeeParameters(),event_data, used);
  }
  event_data.m_measurement.clear();
  return event_data.m_measurement;
}

extendTrack() [2/2]

std::vector< const Trk::MeasurementBase * > & InDet::TRT_TrackExtensionToolCosmics::extendTrack ( const EventContext & ctx, const Trk::TrackParameters * par, InDet::ITRT_TrackExtensionTool::IEventData & virt_event_data, InDet::TRT_DetElementLink_xk::TRT_DetElemUsedMap & used ) const

overridevirtual

Implements InDet::ITRT_TrackExtensionTool.

Definition at line 352 of file TRT_TrackExtensionToolCosmics.cxx.

{
  InDet::TRT_TrackExtensionToolCosmics::EventData &
     event_data=InDet::TRT_TrackExtensionToolCosmics::EventData::getPrivateEventData(virt_event_data);
 
  event_data.m_measurement.clear();
 
  std::vector<Identifier> vecID;
  std::vector<const Trk::TrackParameters*> vecTP;
 
 
 
  std::vector<const Trk::TrackParameters* >* tpars_down=nullptr;
  std::vector<const Trk::TrackParameters* >* tpars_up=nullptr;
  const Trk::Perigee *per=dynamic_cast<const Trk::Perigee *>(par);
  if (!per) {
    msg(MSG::FATAL)<<"Track perigee not found!"<<endmsg;
    return event_data.m_measurement;
  }
 
  if (!event_data.m_trtcontainer) {
    return event_data.m_measurement;
  }
 
  InDet::TRT_DriftCircleContainer::const_iterator
      w = event_data.m_trtcontainer->begin(),we = event_data.m_trtcontainer->end();
  for(; w!=we; ++w) {
     if ((**w).empty()) continue; 
     const Trk::Surface &surf=(**(**w).begin()).detectorElement()->surface();
     Amg::Vector3D pos=intersect(&surf,per);
     Amg::Vector3D locintersec = (surf.transform().inverse())*pos;
     Amg::Vector2D locpos(locintersec.x(), locintersec.y());
     if (pos.perp()<500. || !surf.insideBounds(locpos,50.,50.)) continue;
 
     Amg::Vector3D pos2=surf.transform()*Amg::Vector3D(locintersec.x(),locintersec.y(),0);
     
     Trk::TrackParameters *newpar=nullptr;
     if (auto st = surf.type(); st == Trk::SurfaceType::Plane){
       const Trk::PlaneSurface *plsurf=static_cast<const Trk::PlaneSurface *>(&surf);
       newpar=new Trk::AtaPlane(pos2,per->parameters()[Trk::phi],per->parameters()[Trk::theta],per->parameters()[Trk::qOverP],*plsurf);
     } else if (st == Trk::SurfaceType::Disc){
       const Trk::DiscSurface *discsurf=static_cast<const Trk::DiscSurface *>(&surf);
       newpar=new Trk::AtaDisc(pos2,per->parameters()[Trk::phi],per->parameters()[Trk::theta],per->parameters()[Trk::qOverP],*discsurf);
     }
     vecTP.push_back(newpar);
   }
  tpars_down=new std::vector<const Trk::TrackParameters* >;
  tpars_up=new std::vector<const Trk::TrackParameters* >;
 
  if(!tpars_down ||  !tpars_up) return event_data.m_measurement;
 
  tp_sort_cosmics sorter(per->parameters()[Trk::theta]);
  std::sort(vecTP.begin(),vecTP.end(),sorter);  
  for (const auto *tmppar : vecTP){ 
    if ((per->parameters()[Trk::theta]>M_PI/2 && per->position().z()<tmppar->position().z()) || (per->parameters()[Trk::theta]<M_PI/2 && per->position().z()>tmppar->position().z())) tpars_up->push_back(tmppar); 
    else tpars_down->push_back(tmppar);
 
  }
  if (!tpars_up->empty()) std::reverse(tpars_up->begin(),tpars_up->end());
  
 
  if(tpars_down){
    analyze_tpars(tpars_down,event_data);
  }
  if(tpars_up){
    analyze_tpars(tpars_up, event_data);
  }
 
  //clean up
  if(tpars_up){
    std::vector< const Trk::TrackParameters* >::const_iterator parameterIter= tpars_up->begin();
    for ( ; parameterIter != tpars_up->end(); ++parameterIter) {
      delete *parameterIter;
    }
    delete tpars_up;
  }
  
  if(tpars_down){
    std::vector< const Trk::TrackParameters* >::const_iterator parameterIter = tpars_down->begin();
    for ( ; parameterIter != tpars_down->end(); ++parameterIter) {
      delete *parameterIter;
    }
    delete tpars_down;
  }
 
  msg(MSG::DEBUG)<<"Found "<<event_data.m_measurement.size()<<" driftcircles"<<endmsg;
 
  return event_data.m_measurement;
}

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

finalize()

StatusCode InDet::TRT_TrackExtensionToolCosmics::finalize ( )

overridevirtual

Definition at line 90 of file TRT_TrackExtensionToolCosmics.cxx.

{
   StatusCode sc = AlgTool::finalize(); 
 
   return sc;
}

findBoundarySurface()

Trk::Surface * InDet::TRT_TrackExtensionToolCosmics::findBoundarySurface ( const Trk::TrackParameters & par, Trk::PropDirection dir, InDet::TRT_TrackExtensionToolCosmics::EventData & event_data ) const

protected

Definition at line 460 of file TRT_TrackExtensionToolCosmics.cxx.

{
 
  const EventContext& ctx = Gaudi::Hive::currentContext();
  const Trk::TrackParameters* test=m_extrapolator->extrapolateDirectly(ctx,
                                                                       par,
                                                                       *event_data.m_trtcylinder,
                                                                       dir,true,Trk::muon).release();
  if(test){
    delete test;
    return event_data.m_trtcylinder.get();
  }
 
  test=m_extrapolator->extrapolateDirectly(ctx,
                                           par,
                                           *event_data.m_trtdiscA,dir,true,Trk::muon).release();
  if(test){
    delete test;
    return event_data.m_trtdiscA.get();
  }
 
  test=m_extrapolator->extrapolateDirectly(ctx,
                                           par,
                                           *event_data.m_trtdiscC,dir,true,Trk::muon).release();
  if(test){
    delete test;
    return event_data.m_trtdiscC.get();
  }
 
  return nullptr;
}

findSegment()

Trk::TrackSegment * InDet::TRT_TrackExtensionToolCosmics::findSegment ( const EventContext & ctx, const Trk::TrackParameters * , InDet::ITRT_TrackExtensionTool::IEventData & virt_event_data, InDet::TRT_DetElementLink_xk::TRT_DetElemUsedMap & used ) const

overridevirtual

Implements InDet::ITRT_TrackExtensionTool.

Definition at line 450 of file TRT_TrackExtensionToolCosmics.cxx.

{
  return nullptr;
}

initialize()

StatusCode InDet::TRT_TrackExtensionToolCosmics::initialize ( )

overridevirtual

Definition at line 49 of file TRT_TrackExtensionToolCosmics.cxx.

{
  StatusCode sc = AlgTool::initialize(); 
 
  // Get RIO_OnTrack creator without drift time information
  //
  if(m_riontrackN.retrieve().isFailure()) {
    ATH_MSG_FATAL("Failed to retrieve tool "<< m_riontrackN);
    return StatusCode::FAILURE;
  }
  ATH_MSG_DEBUG("Retrieved tool " << m_riontrackN);
 
  ATH_CHECK( m_extrapolator.retrieve() );
 
  if ((detStore()->retrieve(m_trtid)).isFailure()) {
    ATH_MSG_FATAL("Problem retrieving TRTID helper");
    return StatusCode::FAILURE;
  }
 
  // Get output print level
  //
  m_outputlevel = msg().level()-MSG::DEBUG;
  if(m_outputlevel<=0) {
    m_nprint=0; msg(MSG::DEBUG)<<(*this)<<endmsg;
  }
 
 
  //Initialize container
  ATH_CHECK(m_trtname.initialize());
 
  return sc;
 
 
  
 
}

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 & InDet::ITRT_TrackExtensionTool::interfaceID ( )

inlinestaticinherited

Definition at line 116 of file ITRT_TrackExtensionTool.h.

    {
      return IID_ITRT_TrackExtensionTool;
    }

intersect()

Amg::Vector3D InDet::TRT_TrackExtensionToolCosmics::intersect ( const Trk::Surface * surf, const Trk::Perigee * per )

staticprotected

Definition at line 494 of file TRT_TrackExtensionToolCosmics.cxx.

                                                                                                      {
  // Calculate intersection of helix with silicon module. Assume barrel modules parallel to z-axis, endcap modules perpendicular to z-axis
 
  double sinTheta = std::sin(per->parameters()[3]);
  double r= (std::abs(per->parameters()[Trk::qOverP]) > 1e-10) ? -sinTheta/(per->parameters()[Trk::qOverP]*0.6) : 1e6;
  double xc=per->position().x()-r*std::sin(per->parameters()[2]);
  double yc=per->position().y()+r*std::cos(per->parameters()[2]);
  double phi0=std::atan2(per->position().y()-yc,per->position().x()-xc);
  double theta=per->parameters()[Trk::theta];
  double z0=per->position().z();
 
  if (std::abs(surf->normal().z())>0.5){ // endcap module
    double delta_s=(surf->center().z()-z0)/cos(theta);
    double delta_phi=delta_s*std::sin(theta)/r;
    double x=xc+std::abs(r)*std::cos(phi0+delta_phi);
    double y=yc+std::abs(r)*std::sin(phi0+delta_phi);
    return Amg::Vector3D(x,y,surf->center().z());
 
 
 
  }
  else { // barrel module
    double x1=surf->center().x();
    double y1=surf->center().y();
    double x2=x1+surf->transform()(0,0);
    double y2=y1+surf->transform()(1,0);
    double a=(x2-x1)*(x2-x1)+(y2-y1)*(y2-y1);
    double b=2*( (x2 - x1)*(x1 - xc) + (y2 - y1)*(y1 - yc));
    double c=xc*xc + yc*yc + x1*x1 + y1*y1 - 2*(xc*x1 + yc*y1) - r*r;
    double discr=b*b-4*a*c;
    if (discr<0) return Amg::Vector3D(0,0,0);
    double u1=(-b-std::sqrt(discr))/(2*a);
    double u2=(-b+std::sqrt(discr))/(2*a);
    double u=(std::abs(u1)<std::abs(u2)) ? u1 : u2;
    double x=x1+u*(x2-x1);
    double y=y1+u*(y2-y1);
    double phi=std::atan2(y-yc,x-xc);
    double delta_phi=phi-phi0;
    if (std::abs(std::abs(delta_phi)-2*M_PI)<std::abs(delta_phi)){
      if (delta_phi<0) delta_phi+=2*M_PI;
      else delta_phi-=2*M_PI;
 
    }
    double delta_z=r*delta_phi/std::tan(theta);
 
    double z=z0+delta_z;
    return Amg::Vector3D(x,y,z);
  }
}

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

newEvent()

std::unique_ptr< InDet::ITRT_TrackExtensionTool::IEventData > InDet::TRT_TrackExtensionToolCosmics::newEvent ( const EventContext & ctx ) const

overridevirtual

Implements InDet::ITRT_TrackExtensionTool.

Definition at line 162 of file TRT_TrackExtensionToolCosmics.cxx.

{
  //create the boundary surfaces
  //
  SG::ReadHandle<TRT_DriftCircleContainer> trtcontainer(m_trtname, ctx);
  if(not trtcontainer.isValid() && m_outputlevel<=0) {
    std::stringstream msg;
    msg << "Missing TRT_DriftCircleContainer " << m_trtname.key();
    throw std::runtime_error( msg.str() );
  }
 
  std::unique_ptr<EventData> event_data(new EventData(trtcontainer.cptr()));
 
  Amg::RotationMatrix3D r; r.setIdentity();
  Amg::Transform3D t = Amg::Transform3D(r * Amg::Translation3D(Amg::Vector3D::Zero()));
  event_data->m_trtcylinder = std::make_unique<Trk::CylinderSurface>(t,1150.,3000.);
  Amg::Transform3D transf = Amg::Transform3D(r * Amg::Translation3D(Amg::Vector3D(0.,0.,3000)));
  event_data->m_trtdiscA   = std::make_unique<Trk::DiscSurface>    (transf,1.,1200.);
  transf = Amg::Transform3D(r * Amg::Translation3D(Amg::Vector3D(0.,0.,-3000)));
  event_data->m_trtdiscC   = std::make_unique<Trk::DiscSurface>    (transf,1.,1200.);
 
  return std::unique_ptr<InDet::ITRT_TrackExtensionTool::IEventData>(event_data.release());
 
}

newTrack()

Trk::Track * InDet::TRT_TrackExtensionToolCosmics::newTrack ( const EventContext & ctx, const Trk::Track & , InDet::ITRT_TrackExtensionTool::IEventData & virt_event_data, InDet::TRT_DetElementLink_xk::TRT_DetElemUsedMap & used ) const

overridevirtual

Implements InDet::ITRT_TrackExtensionTool.

Definition at line 550 of file TRT_TrackExtensionToolCosmics.cxx.

{ 
  return nullptr;
}

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

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 >, and AthCheckedComponent<::AthAlgTool >.

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_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_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

ToolHandle<Trk::IExtrapolator> InDet::TRT_TrackExtensionToolCosmics::m_extrapolator {this, "Extrapolator", "Trk::Extrapolator/InDetExtrapolator"}

protected

Definition at line 112 of file TRT_TrackExtensionToolCosmics.h.

{this, "Extrapolator", "Trk::Extrapolator/InDetExtrapolator"};

m_idHelper

AtlasDetectorID* InDet::TRT_TrackExtensionToolCosmics::m_idHelper {}

protected

Definition at line 148 of file TRT_TrackExtensionToolCosmics.h.

{}; //<! Detector ID helper

m_nprint

int InDet::TRT_TrackExtensionToolCosmics::m_nprint {}

protected

Definition at line 118 of file TRT_TrackExtensionToolCosmics.h.

{}         ; // Kind of print

m_outputlevel

int InDet::TRT_TrackExtensionToolCosmics::m_outputlevel {}

protected

Definition at line 117 of file TRT_TrackExtensionToolCosmics.h.

{}    ; // Print level

m_riontrackN

ToolHandle<Trk::IRIO_OnTrackCreator> InDet::TRT_TrackExtensionToolCosmics::m_riontrackN {this, "RIOonTrackToolNoDr", "InDet::TRT_DriftCircleOnTrackNoDriftTimeTool/TRT_DriftCircleOnTrackNoDriftTimeTool"}

protected

Definition at line 114 of file TRT_TrackExtensionToolCosmics.h.

{this, "RIOonTrackToolNoDr", "InDet::TRT_DriftCircleOnTrackNoDriftTimeTool/TRT_DriftCircleOnTrackNoDriftTimeTool"}; //

m_roadwidth

DoubleProperty InDet::TRT_TrackExtensionToolCosmics::m_roadwidth

protected

Initial value:

{this, "RoadWidth", 10.,
    "Max width of the road"}

Definition at line 119 of file TRT_TrackExtensionToolCosmics.h.

                                {this, "RoadWidth", 10.,
    "Max width of the road"};

m_roadwidth_locz

DoubleProperty InDet::TRT_TrackExtensionToolCosmics::m_roadwidth_locz

protected

Initial value:

{this, "BoundaryLocZTolerance", 10.,
    "Max width of the road along the straw"}

Definition at line 121 of file TRT_TrackExtensionToolCosmics.h.

                                     {this, "BoundaryLocZTolerance", 10.,
    "Max width of the road along the straw"};

m_searchNeighbour

BooleanProperty InDet::TRT_TrackExtensionToolCosmics::m_searchNeighbour

protected

Initial value:

{this, "SearchNeighbour", false,
    "Also search neighbouring detector elements?"}

Definition at line 123 of file TRT_TrackExtensionToolCosmics.h.

                                       {this, "SearchNeighbour", false,
    "Also search neighbouring detector elements?"};

m_trtid

const TRT_ID* InDet::TRT_TrackExtensionToolCosmics::m_trtid {}

protected

Definition at line 150 of file TRT_TrackExtensionToolCosmics.h.

{} ;

m_trtname

SG::ReadHandleKey<TRT_DriftCircleContainer> InDet::TRT_TrackExtensionToolCosmics::m_trtname {this,"TRT_ClustersContainer","TRT_DriftCircles","RHK to retrieve TRT_DriftCircles"}

protected

Definition at line 126 of file TRT_TrackExtensionToolCosmics.h.

{this,"TRT_ClustersContainer","TRT_DriftCircles","RHK to retrieve TRT_DriftCircles"};

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.

---

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

• TRT_TrackExtensionToolCosmics.h
• TRT_TrackExtensionToolCosmics.cxx