Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
Public Member Functions | Protected Member Functions | Private Attributes | List of all members
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
 
std::vector< Identifier > getRdoIdList (const FPGATrackSimHit &hit) const
 
template<typename XAOD_CLUSTER >
StatusCode matchTrackMeasurements (const EventContext &ctx, const DataVector< XAOD_CLUSTER > &clusterContainer, const std::vector< Identifier > &rdoIDs, std::vector< ActsTrk::ATLASUncalibSourceLink > &measurements) const
 
StatusCode findPrototrackMeasurements (const EventContext &ctx, const xAOD::PixelClusterContainer &pixelClusterContainer, const xAOD::StripClusterContainer &stripClusterContainer, std::vector< ActsTrk::ATLASUncalibSourceLink > &measurements, const std::vector< FPGATrackSimHit > &hits) 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 13 of file FPGAActsTrkConverter.cxx.

15  : base_class(type,name,parent) { }

◆ ~FPGAActsTrkConverter()

virtual FPGAActsTrkConverter::~FPGAActsTrkConverter ( )
virtualdefault

Member Function Documentation

◆ findPrototrackMeasurements()

StatusCode FPGAActsTrkConverter::findPrototrackMeasurements ( const EventContext &  ctx,
const xAOD::PixelClusterContainer &  pixelClusterContainer,
const xAOD::StripClusterContainer &  stripClusterContainer,
std::vector< ActsTrk::ATLASUncalibSourceLink > &  measurements,
const std::vector< FPGATrackSimHit > &  hits 
) const
protected

Definition at line 77 of file FPGAActsTrkConverter.cxx.

81  {
82  if (hits.empty()) {
83  ATH_MSG_ERROR("Found FPGATrack without hits");
84  return StatusCode::FAILURE;
85  }
86  // TODO: move from loops over cluster container to links
87  for (const FPGATrackSimHit& h : hits) {
88  if (h.isReal()) {
89  std::vector<Identifier> ids = getRdoIdList(h);
90  if (h.isPixel()) {
91  ATH_MSG_DEBUG("Looking for Pixel cluster to match");
92  ATH_CHECK(matchTrackMeasurements<xAOD::PixelCluster>(ctx, pixelClusterContainer, ids, measurements));
93  }
94  else if (h.isStrip()) {
95  ATH_MSG_DEBUG("Looking for Strip cluster to match");
96  ATH_CHECK(matchTrackMeasurements<xAOD::StripCluster>(ctx, stripClusterContainer, ids, measurements));
97  }
98  else {
99  ATH_MSG_ERROR("FPGA hit not classified as pixel or strip");
100  return StatusCode::FAILURE;
101  }
102  }
103  else
104  ATH_MSG_DEBUG("Skipping hit as non-Real");
105  }
106  return StatusCode::SUCCESS;
107 }

◆ 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 55 of file FPGAActsTrkConverter.cxx.

59  {
60 
61  ATH_MSG_INFO("Creating Acts proto-tracks from FPGA tracks...");
62 
63  for (const FPGATrackSimTrack& track : tracks) {
64  if (not track.passedOR()) continue;
65  std::vector<ActsTrk::ATLASUncalibSourceLink> points;
66  const std::vector <FPGATrackSimHit>& hits = track.getFPGATrackSimHits();
67  ATH_CHECK(findPrototrackMeasurements(ctx, pixelContainer, stripContainer, points, hits));
68  if (points.size()){
69  ATH_MSG_INFO("\tMaking a proto-track with " << points.size() << " clusters");
70  std::unique_ptr<Acts::BoundTrackParameters> inputPerigee = makeParams(track);
71  foundProtoTracks.emplace_back(points, std::move(inputPerigee));
72  }
73  }
74  return StatusCode::SUCCESS;
75 }

◆ 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 30 of file FPGAActsTrkConverter.cxx.

35  {
36 
37  ATH_MSG_INFO("Creating Acts proto-tracks from FPGA roads...");
38 
39  if (hitsInRoads.size() > 0) {
40 
41  for(size_t roadIndex=0; roadIndex<=hitsInRoads.size()-1;roadIndex++) {
42  std::vector<ActsTrk::ATLASUncalibSourceLink> points;
43  ATH_CHECK(findPrototrackMeasurements(ctx, pixelContainer, stripContainer, points, hitsInRoads.at(roadIndex)));
44  if (points.size()) {
45  std::unique_ptr<Acts::BoundTrackParameters> inputPerigee = makeParams(roads.at(roadIndex));
46  foundProtoTracks.emplace_back(points, std::move(inputPerigee));
47  ATH_MSG_INFO("Made a prototrack with " << points.size() << " measurements");
48  }
49  }
50  }
51 
52  return StatusCode::SUCCESS;
53 }

◆ getRdoIdList()

std::vector< Identifier > FPGAActsTrkConverter::getRdoIdList ( const FPGATrackSimHit &  hit) const
protected

Definition at line 109 of file FPGAActsTrkConverter.cxx.

110 {
111  std::vector<Identifier> ids;
112  const auto& idHashVec = hit.getIDHashVec();
113  const auto& phiIndexVec = hit.getPhiIndexVec();
114  const auto& etaIndexVec = hit.getEtaIndexVec();
115 
116  ids.reserve(idHashVec.size());
117 
118  if (hit.isPixel())
119  {
120  for (size_t i = 0; i < idHashVec.size(); ++i)
121  {
122  ids.emplace_back(m_pixelId->pixel_id(
123  m_pixelId->wafer_id(idHashVec[i]),
124  phiIndexVec[i],
125  etaIndexVec[i]
126  ));
127  }
128  }
129  else if (hit.isStrip())
130  {
131  for (size_t i = 0; i < idHashVec.size(); ++i)
132  {
133  ids.emplace_back(m_SCTId->strip_id(
134  m_SCTId->wafer_id(idHashVec[i]),
135  phiIndexVec[i]
136  ));
137  }
138  }
139  return ids;
140 }

◆ initialize()

StatusCode FPGAActsTrkConverter::initialize ( )
finaloverridevirtual

Definition at line 18 of file FPGAActsTrkConverter.cxx.

18  {
19 
20  ATH_MSG_DEBUG("Initializing FPGAActsTrkConverter...");
21 
22  // Get SCT & pixel Identifier helpers
23  ATH_CHECK(detStore()->retrieve(m_pixelId, "PixelID"));
24  ATH_CHECK(detStore()->retrieve(m_SCTId, "SCT_ID"));
25 
26  return StatusCode::SUCCESS;
27 
28 }

◆ makeParams() [1/2]

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

Definition at line 174 of file FPGAActsTrkConverter.cxx.

174  {
175  using namespace Acts::UnitLiterals;
176 
177  std::shared_ptr<const Acts::Surface> actsSurface = Acts::Surface::makeShared<Acts::PerigeeSurface>(Acts::Vector3(0., 0., 0.));
178  Acts::BoundVector params;
179 
180  constexpr double GeVToMeV = 1000;
181  double d0=0.; //?
182  double z0=0.; //?
183  double phi=road.getX();
184  double eta=0.2;
185  double theta=2*std::atan(std::exp(-eta));
186  double qop = (std::abs(road.getY()) > 1E-9) ? road.getY()/GeVToMeV : 1E-12;
187  double t=0.; //?
188  ATH_MSG_DEBUG("\tphi=" <<phi << " eta=" << eta << " qop=" << qop);
189 
190  params << d0, z0, phi, theta, qop, t;
191 
192  // Covariance - TODO
193  Acts::BoundSquareMatrix cov = Acts::BoundSquareMatrix::Identity();
194  cov *= (GeVToMeV*GeVToMeV);
195 
196  // some ACTS paperwork
197  Trk::ParticleHypothesis hypothesis = Trk::pion;
198  float mass = Trk::ParticleMasses::mass[hypothesis] * Acts::UnitConstants::MeV;
199  Acts::PdgParticle absPdg = Acts::makeAbsolutePdgParticle(Acts::ePionPlus);
200  Acts::ParticleHypothesis actsHypothesis{
201  absPdg, mass, Acts::AnyCharge{1.0f}};
202 
203  return std::make_unique<Acts::BoundTrackParameters>(actsSurface, params,
204  cov, actsHypothesis);
205 
206 }

◆ makeParams() [2/2]

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

Definition at line 209 of file FPGAActsTrkConverter.cxx.

209  {
210 
211  using namespace Acts::UnitLiterals;
212  std::shared_ptr<const Acts::Surface> actsSurface = Acts::Surface::makeShared<Acts::PerigeeSurface>(Acts::Vector3(0., 0., 0.));
213  Acts::BoundVector params;
214 
215  constexpr double GeVToMeV = 1000;
216  double d0=track.getD0();
217  double z0=track.getZ0();
218  double phi=track.getPhi();
219  double theta=track.getTheta();
220  double qop=track.getQOverPt();
221  double t=0.;
222 
223  params << d0, z0, phi, theta, qop, t;
224  ATH_MSG_DEBUG("\td0= " << d0 << " z0=" <<z0 << " phi=" <<phi << " theta=" << theta<< " qoverp=" << qop);
225 
226  // Covariance - let's be honest and say we have no clue ;-)
227  Acts::BoundSquareMatrix cov = Acts::BoundSquareMatrix::Identity();
228  cov *= (GeVToMeV*GeVToMeV);
229 
230  // some ACTS paperwork
231  Trk::ParticleHypothesis hypothesis = Trk::pion;
232  float mass = Trk::ParticleMasses::mass[hypothesis] * Acts::UnitConstants::MeV;
233  Acts::PdgParticle absPdg = Acts::makeAbsolutePdgParticle(Acts::ePionPlus);
234  Acts::ParticleHypothesis actsHypothesis{
235  absPdg, mass, Acts::AnyCharge{1.0f}};
236 
237  return std::make_unique<Acts::BoundTrackParameters>(actsSurface, params,
238  cov, actsHypothesis);
239 
240 }

◆ matchTrackMeasurements()

template<typename XAOD_CLUSTER >
StatusCode FPGAActsTrkConverter::matchTrackMeasurements ( const EventContext &  ctx,
const DataVector< XAOD_CLUSTER > &  clusterContainer,
const std::vector< Identifier > &  rdoIDs,
std::vector< ActsTrk::ATLASUncalibSourceLink > &  measurements 
) const
protected

Definition at line 144 of file FPGAActsTrkConverter.cxx.

148 {
149  for (const XAOD_CLUSTER* cl : clusterContainer) {
150  const auto& rdoList = cl->rdoList();
151  if(rdoIDs.size() != rdoList.size()) continue;
152  size_t matchedCounter=0;
153  for (const Identifier& id: rdoIDs){
154  if(std::find(rdoList.begin(), rdoList.end(),id) != rdoList.end()) matchedCounter++;
155  }
156  if(matchedCounter == rdoList.size())
157  measurements.emplace_back(ActsTrk::makeATLASUncalibSourceLink(&clusterContainer, cl->index(), ctx));
158  else if (matchedCounter>0) {
159  std::stringstream ss;
160  ss << "List of xAOD cluster rdoIDs:\n";
161  for (const Identifier& id: rdoList) ss << id << "\n";
162 
163  ss << "\nList of FPGA cluster rdoIDs:\n";
164  for (const Identifier& id: rdoIDs) ss << id << "\n";
165 
166  ATH_MSG_ERROR(std::format("Noticed rdoID mismatch: commonHits = {} | xAODClusterHits = {} | FPGAClusterHits = {}\n{}",
167  matchedCounter, rdoList.size(), rdoIDs.size(), ss.str()));
168  return StatusCode::FAILURE;
169  }
170  }
171  return StatusCode::SUCCESS;
172 }

Member Data Documentation

◆ m_pixelId

const PixelID* FPGAActsTrkConverter::m_pixelId {nullptr}
private

Definition at line 53 of file FPGAActsTrkConverter.h.

◆ m_SCTId

const SCT_ID* FPGAActsTrkConverter::m_SCTId {nullptr}
private

Definition at line 54 of file FPGAActsTrkConverter.h.

The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
FPGATrackSimHit::getEtaIndexVec
const std::vector< int > & getEtaIndexVec() const
Definition: FPGATrackSimHit.h:189
TRTCalib_Extractor.hits
hits
Definition: TRTCalib_Extractor.py:35
FPGATrackSimHit::isStrip
bool isStrip() const
Definition: FPGATrackSimHit.h:65
PowhegControl_ttHplus_NLO.ss
ss
Definition: PowhegControl_ttHplus_NLO.py:83
phi
Scalar phi() const
phi method
Definition: AmgMatrixBasePlugin.h:67
vtune_athena.format
format
Definition: vtune_athena.py:14
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
find
std::string find(const std::string &s)
return a remapped string
Definition: hcg.cxx:135
FPGAActsTrkConverter::getRdoIdList
std::vector< Identifier > getRdoIdList(const FPGATrackSimHit &hit) const
Definition: FPGAActsTrkConverter.cxx:109
Base_Fragment.mass
mass
Definition: Sherpa_i/share/common/Base_Fragment.py:59
FPGATrackSimTrack
Definition: FPGATrackSimTrack.h:18
FPGAActsTrkConverter::m_SCTId
const SCT_ID * m_SCTId
Definition: FPGAActsTrkConverter.h:54
eta
Scalar eta() const
pseudorapidity method
Definition: AmgMatrixBasePlugin.h:83
theta
Scalar theta() const
theta method
Definition: AmgMatrixBasePlugin.h:75
FPGATrackSimHit::getPhiIndexVec
const std::vector< int > & getPhiIndexVec() const
Definition: FPGATrackSimHit.h:188
python.SystemOfUnits.MeV
int MeV
Definition: SystemOfUnits.py:154
plotBeamSpotVxVal.cov
cov
Definition: plotBeamSpotVxVal.py:201
FPGATrackSimHit::getIDHashVec
const std::vector< unsigned > & getIDHashVec() const
Definition: FPGATrackSimHit.h:190
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
drawFromPickle.exp
exp
Definition: drawFromPickle.py:36
DataPrepToActsConfig.stripContainer
stripContainer
Definition: DataPrepToActsConfig.py:10
FPGATrackSimRoad::getX
float getX() const
Definition: FPGATrackSimRoad.h:87
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
PixelID::wafer_id
Identifier wafer_id(int barrel_ec, int layer_disk, int phi_module, int eta_module) const
For a single crystal.
Definition: PixelID.h:364
drawFromPickle.atan
atan
Definition: drawFromPickle.py:36
FPGATrackSimHit
Definition: FPGATrackSimHit.h:41
Trk::ParticleHypothesis
ParticleHypothesis
Definition: ParticleHypothesis.h:25
ActsTrk::makeATLASUncalibSourceLink
ATLASUncalibSourceLink makeATLASUncalibSourceLink(const xAOD::UncalibratedMeasurementContainer *container, std::size_t index, [[maybe_unused]] const EventContext &ctx)
Definition: ATLASSourceLink.h:30
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
lumiFormat.i
int i
Definition: lumiFormat.py:85
h
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
Trk::pion
@ pion
Definition: ParticleHypothesis.h:29
TRT::Track::d0
@ d0
Definition: InnerDetector/InDetCalibEvent/TRT_CalibData/TRT_CalibData/TrackInfo.h:62
test_pyathena.parent
parent
Definition: test_pyathena.py:15
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
TRT::Track::z0
@ z0
Definition: InnerDetector/InDetCalibEvent/TRT_CalibData/TRT_CalibData/TrackInfo.h:63
FPGATrackSimHit::isPixel
bool isPixel() const
Definition: FPGATrackSimHit.h:64
FPGAActsTrkConverter::m_pixelId
const PixelID * m_pixelId
Definition: FPGAActsTrkConverter.h:53
Trk::ParticleMasses::mass
constexpr double mass[PARTICLEHYPOTHESES]
the array of masses
Definition: ParticleHypothesis.h:53
python.PyKernel.detStore
detStore
Definition: PyKernel.py:41
FPGAActsTrkConverter::findPrototrackMeasurements
StatusCode findPrototrackMeasurements(const EventContext &ctx, const xAOD::PixelClusterContainer &pixelClusterContainer, const xAOD::StripClusterContainer &stripClusterContainer, std::vector< ActsTrk::ATLASUncalibSourceLink > &measurements, const std::vector< FPGATrackSimHit > &hits) const
Definition: FPGAActsTrkConverter.cxx:77
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
python.subdetectors.mmg.ids
ids
Definition: mmg.py:8
VP1PartSpect::E
@ E
Definition: VP1PartSpectFlags.h:21
xAOD::ParticleHypothesis
ParticleHypothesis
Definition: TrackingPrimitives.h:193
FPGATrackSimRoad::getY
float getY() const
Definition: FPGATrackSimRoad.h:88
PixelID::pixel_id
Identifier pixel_id(int barrel_ec, int layer_disk, int phi_module, int eta_module, int phi_index, int eta_index) const
For an individual pixel.
Definition: PixelID.h:432
FPGAActsTrkConverter::makeParams
std::unique_ptr< Acts::BoundTrackParameters > makeParams(const FPGATrackSimRoad &road) const
Definition: FPGAActsTrkConverter.cxx:174
xAOD::track
@ track
Definition: TrackingPrimitives.h:513
PowhegControl_ttFCNC_NLO.params
params
Definition: PowhegControl_ttFCNC_NLO.py:226
SCT_ID::wafer_id
Identifier wafer_id(int barrel_ec, int layer_disk, int phi_module, int eta_module, int side) const
For a single side of module.
Definition: SCT_ID.h:464
DataPrepToActsConfig.pixelContainer
pixelContainer
Definition: DataPrepToActsConfig.py:9
dq_make_web_display.cl
cl
print [x.__class__ for x in toList(dqregion.getSubRegions()) ]
Definition: dq_make_web_display.py:26
jobOptions.points
points
Definition: jobOptions.GenevaPy8_Zmumu.py:97
SCT_ID::strip_id
Identifier strip_id(int barrel_ec, int layer_disk, int phi_module, int eta_module, int side, int strip) const
For an individual strip.
Definition: SCT_ID.h:535
Identifier
Definition: IdentifierFieldParser.cxx:14
Generated on Fri Apr 18 2025 21:34:30 for ATLAS Offline Software by   doxygen 1.8.18