FPGAActsTrkConverter Class Reference

#include <FPGAActsTrkConverter.h>

Inheritance diagram for FPGAActsTrkConverter:

Collaboration diagram for FPGAActsTrkConverter:

Public Member Functions
	FPGAActsTrkConverter (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~FPGAActsTrkConverter ()=default
virtual StatusCode	initialize () override final
virtual StatusCode	findProtoTracks (const EventContext &ctx, const xAOD::PixelClusterContainer &pixelContainer, const xAOD::StripClusterContainer &stripContainer, std::vector< ActsTrk::ProtoTrack > &foundProtoTracks, const std::vector< std::vector< FPGATrackSimHit >> &hitsInRoads, const std::vector< FPGATrackSimRoad > &roads) const override final
virtual StatusCode	findProtoTracks (const EventContext &ctx, const xAOD::PixelClusterContainer &pixelContainer, const xAOD::StripClusterContainer &stripContainer, std::vector< ActsTrk::ProtoTrack > &foundProtoTracks, const std::vector< FPGATrackSimTrack > &tracks) const override final

Protected Member Functions
std::unique_ptr< Acts::BoundTrackParameters >	makeParams (const FPGATrackSimRoad &road) const
std::unique_ptr< Acts::BoundTrackParameters >	makeParams (const FPGATrackSimTrack &track) const

Private Attributes
const PixelID *	m_pixelId {nullptr}
const SCT_ID *	m_SCTId {nullptr}

Detailed Description

Definition at line 17 of file FPGAActsTrkConverter.h.

Constructor & Destructor Documentation

FPGAActsTrkConverter()

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

Definition at line 11 of file FPGAActsTrkConverter.cxx.

                                 : base_class(type,name,parent) { }

~FPGAActsTrkConverter()

virtual FPGAActsTrkConverter::~FPGAActsTrkConverter ( )

virtualdefault

Member Function Documentation

findProtoTracks() [1/2]

StatusCode FPGAActsTrkConverter::findProtoTracks ( const EventContext &  ctx, const xAOD::PixelClusterContainer &  pixelContainer, const xAOD::StripClusterContainer &  stripContainer, std::vector< ActsTrk::ProtoTrack > &  foundProtoTracks, const std::vector< FPGATrackSimTrack > &  tracks  ) const

finaloverridevirtual

Definition at line 76 of file FPGAActsTrkConverter.cxx.

                                                                    {
    
    ATH_MSG_INFO("Creating Acts proto-tracks from FPGA tracks...");
 
    if (tracks.size()>1){
      for(const FPGATrackSimTrack& track: tracks) { 
        std::vector<ActsTrk::ATLASUncalibSourceLink> points;  
        const std::vector <FPGATrackSimHit>& hits = track.getFPGATrackSimHits();
 
        if (hits.size()>0){
 
          // TODO: move from loops over cluster container to links 
 
          for(const FPGATrackSimHit& h : hits){
            IdentifierHash hash = h.getIdentifierHash();
            if (h.isReal()){
              if (h.isPixel()){
                Identifier wafer_id = m_pixelId->wafer_id(hash);
                Identifier id = m_pixelId->pixel_id(wafer_id, h.getPhiIndex(), h.getEtaIndex()); 
                for (const xAOD::PixelCluster *cl : pixelContainer){
                  if (id == cl->rdoList().front()) points.push_back(ActsTrk::makeATLASUncalibSourceLink(&pixelContainer,cl->index(),ctx));
                }
              }
              if (h.isStrip()) {
                int strip = static_cast<int>(h.getPhiCoord());
                Identifier wafer_id = m_SCTId->wafer_id(hash);
                Identifier id = m_SCTId->strip_id(wafer_id, strip); 
                for (const xAOD::StripCluster *cl : stripContainer){
                  if (id == cl->rdoList().front()) points.push_back(ActsTrk::makeATLASUncalibSourceLink(&stripContainer,cl->index(),ctx)); 
                }
              }
            }
          }
          ATH_MSG_INFO("\tMade a proto-track with " <<points.size()<<" clusters");
        
          // Make the intput perigee
          std::unique_ptr<Acts::BoundTrackParameters> inputPerigee = makeParams(track);
          foundProtoTracks.emplace_back(points,std::move(inputPerigee));
        }
      }
    }
 
    return StatusCode::SUCCESS;
}

findProtoTracks() [2/2]

StatusCode FPGAActsTrkConverter::findProtoTracks ( const EventContext &  ctx, const xAOD::PixelClusterContainer &  pixelContainer, const xAOD::StripClusterContainer &  stripContainer, std::vector< ActsTrk::ProtoTrack > &  foundProtoTracks, const std::vector< std::vector< FPGATrackSimHit >> &  hitsInRoads, const std::vector< FPGATrackSimRoad > &  roads  ) const

finaloverridevirtual

Definition at line 28 of file FPGAActsTrkConverter.cxx.

                                                                  {
 
    ATH_MSG_INFO("Creating Acts proto-tracks from FPGA roads...");
 
    if (hitsInRoads.size() > 0) {
 
      for(size_t roadIndex=0; roadIndex<=hitsInRoads.size()-1;roadIndex++) { 
        std::vector<ActsTrk::ATLASUncalibSourceLink> points;  
 
        if (hitsInRoads.at(roadIndex).size()>0) {
        // TODO: move from loops over cluster container to links 
          for(const FPGATrackSimHit& h : hitsInRoads.at(roadIndex)){
            IdentifierHash hash = h.getIdentifierHash();
            if (h.isReal()) {
              if (h.isPixel()){
                Identifier wafer_id = m_pixelId->wafer_id(hash);
                Identifier id = m_pixelId->pixel_id(wafer_id, h.getPhiIndex(), h.getEtaIndex()); 
                for (const xAOD::PixelCluster *cl : pixelContainer){
                  if (id == cl->rdoList().front()) points.push_back(ActsTrk::makeATLASUncalibSourceLink(&pixelContainer,cl->index(),ctx));
                }
              }
              if (h.isStrip()) {
                int strip = static_cast<int>(h.getPhiCoord());
                Identifier wafer_id = m_SCTId->wafer_id(hash);
                Identifier id = m_SCTId->strip_id(wafer_id, strip);
                for (const xAOD::StripCluster *cl : stripContainer){
                  if (id == cl->rdoList().front()) points.push_back(ActsTrk::makeATLASUncalibSourceLink(&stripContainer,cl->index(),ctx)); 
                }
              }
            }
          }
          ATH_MSG_INFO("\tMade a proto-track with " <<points.size()<<" clusters");
        
          // Make the input perigee
          std::unique_ptr<Acts::BoundTrackParameters> inputPerigee = makeParams(roads.at(roadIndex));
          foundProtoTracks.emplace_back(points,std::move(inputPerigee));   
        }
      }
    }
 
    return StatusCode::SUCCESS;
}

initialize()

StatusCode FPGAActsTrkConverter::initialize ( )

finaloverridevirtual

Definition at line 16 of file FPGAActsTrkConverter.cxx.

                                            {
 
  ATH_MSG_DEBUG("Initializing FPGAActsTrkConverter...");
 
  // Get SCT & pixel Identifier helpers
  ATH_CHECK(detStore()->retrieve(m_pixelId, "PixelID"));
  ATH_CHECK(detStore()->retrieve(m_SCTId, "SCT_ID"));
 
  return StatusCode::SUCCESS;
 
}

makeParams() [1/2]

std::unique_ptr< Acts::BoundTrackParameters > FPGAActsTrkConverter::makeParams ( const FPGATrackSimRoad &  road ) const

protected

Definition at line 125 of file FPGAActsTrkConverter.cxx.

                                                                                                           {
  using namespace Acts::UnitLiterals;
 
  std::shared_ptr<const Acts::Surface> actsSurface = Acts::Surface::makeShared<Acts::PerigeeSurface>(Acts::Vector3(0., 0., 0.));
  Acts::BoundVector params;
 
  constexpr double GeVToMeV = 1000;
  double d0=0.; //?
  double z0=0.; //?
  double phi=road.getX();
  double eta=0.2;
  double theta=2*std::atan(std::exp(-eta));
  double qop=road.getY()/GeVToMeV; //
  double t=0.; //?
  ATH_MSG_DEBUG("\tphi=" <<phi << " eta=" << eta << " qop=" << qop);
 
  params << d0, z0, phi, theta, qop, t; 
 
  // Covariance - TODO
  Acts::BoundSquareMatrix cov = Acts::BoundSquareMatrix::Identity();
  cov *= (GeVToMeV*GeVToMeV); 
 
  // some ACTS paperwork 
  Trk::ParticleHypothesis hypothesis = Trk::pion;
  float mass = Trk::ParticleMasses::mass[hypothesis] * Acts::UnitConstants::MeV;
  Acts::PdgParticle absPdg = Acts::makeAbsolutePdgParticle(Acts::ePionPlus);
  Acts::ParticleHypothesis actsHypothesis{
    absPdg, mass, Acts::AnyCharge{1.0f}};
 
  return std::make_unique<Acts::BoundTrackParameters>(actsSurface, params,
                                    cov, actsHypothesis);
 
}

makeParams() [2/2]

std::unique_ptr< Acts::BoundTrackParameters > FPGAActsTrkConverter::makeParams ( const FPGATrackSimTrack &  track ) const

protected

Definition at line 160 of file FPGAActsTrkConverter.cxx.

                                                                                                             {
 
  using namespace Acts::UnitLiterals;
  std::shared_ptr<const Acts::Surface> actsSurface = Acts::Surface::makeShared<Acts::PerigeeSurface>(Acts::Vector3(0., 0., 0.));
  Acts::BoundVector params;
 
  constexpr double GeVToMeV = 1000;
  double d0=track.getD0();
  double z0=track.getZ0();
  double phi=track.getPhi();
  double theta=track.getTheta();
  double qop=track.getQOverPt();
  double t=0.;
 
  params << d0, z0, phi, theta, qop, t;  
  ATH_MSG_DEBUG("\td0= " << d0 << " z0=" <<z0 << " phi=" <<phi << " theta=" << theta<< " qoverp=" << qop);
 
  // Covariance - let's be honest and say we have no clue ;-) 
  Acts::BoundSquareMatrix cov = Acts::BoundSquareMatrix::Identity();
  cov *= (GeVToMeV*GeVToMeV); 
 
  // some ACTS paperwork 
  Trk::ParticleHypothesis hypothesis = Trk::pion;
  float mass = Trk::ParticleMasses::mass[hypothesis] * Acts::UnitConstants::MeV;
  Acts::PdgParticle absPdg = Acts::makeAbsolutePdgParticle(Acts::ePionPlus);
  Acts::ParticleHypothesis actsHypothesis{
    absPdg, mass, Acts::AnyCharge{1.0f}};
 
  return std::make_unique<Acts::BoundTrackParameters>(actsSurface, params,
                                    cov, actsHypothesis);
 
}

Member Data Documentation

m_pixelId

const PixelID* FPGAActsTrkConverter::m_pixelId {nullptr}

private

Definition at line 41 of file FPGAActsTrkConverter.h.

m_SCTId

const SCT_ID* FPGAActsTrkConverter::m_SCTId {nullptr}

private

Definition at line 42 of file FPGAActsTrkConverter.h.

---

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

• FPGAActsTrkConverter.h
• FPGAActsTrkConverter.cxx