Skip to content

create pT2 object in Segment Linking #138

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
YonsiG wants to merge 5 commits into
base: CMSSW_14_1_0_pre3_LST_X_LSTCore_realfiles_batch1_devel
Choose a base branch
from
Draft

create pT2 object in Segment Linking #138

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
cbfbd3c
first half of PT2 codes update
YonsiG Nov 20, 2024
5035716
successfully ran the LST PT2
YonsiG Dec 1, 2024
c5de71e
add T4 into performance plots
YonsiG Dec 1, 2024
f8d0fc0
Update RecoTracker/LSTCore/src/alpaka/PT2.h
YonsiG Dec 6, 2024
46cbc28
Update RecoTracker/LSTCore/src/alpaka/Segment.h
YonsiG Dec 6, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 4 additions & 0 deletions RecoTracker/LSTCore/interface/Constants.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -41,6 +41,7 @@ namespace lst {

constexpr unsigned int n_max_pixel_md_per_modules = 2 * n_max_pixel_segments_per_module;

constexpr unsigned int n_max_pixel_segments = 15000;
Copy link

@slava77 slava77 Dec 4, 2024

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

can be pt2 already

Suggested change
constexpr unsigned int n_max_pixel_segments = 15000;
constexpr unsigned int n_max_pt2s = 15000;

constexpr unsigned int n_max_pixel_triplets = 5000;
constexpr unsigned int n_max_pixel_quintuplets = 15000;

Expand All @@ -60,6 +61,9 @@ namespace lst {
struct Params_T3 {
static constexpr int kLayers = 3, kHits = 6;
};
struct Params_pT2 {
static constexpr int kLayers = 4, kHits = 8;
};
struct Params_pT3 {
static constexpr int kLayers = 5, kHits = 10;
};
Expand Down
224 changes: 223 additions & 1 deletion RecoTracker/LSTCore/src/alpaka/Event.dev.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -12,6 +12,7 @@ void lst::Event<Acc3D>::init(bool verbose) {
trackCandidatesInGPU = nullptr;
pixelTripletsInGPU = nullptr;
pixelQuintupletsInGPU = nullptr;
PT2sInGPU = nullptr;
Copy link

@slava77 slava77 Dec 4, 2024

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

variables and functions should start with lower case, upper case is for type names

rangesInGPU = nullptr;

hitsInCPU = nullptr;
Expand All @@ -24,6 +25,7 @@ void lst::Event<Acc3D>::init(bool verbose) {
quintupletsInCPU = nullptr;
pixelTripletsInCPU = nullptr;
pixelQuintupletsInCPU = nullptr;
PT2sInCPU = nullptr;

//reset the arrays
for (int i = 0; i < 6; i++) {
Expand Down Expand Up @@ -107,6 +109,11 @@ void lst::Event<Acc3D>::resetEvent() {
delete pixelQuintupletsBuffers;
pixelQuintupletsInGPU = nullptr;
}
if (PT2sInGPU) {
delete PT2sInGPU;
delete PT2sBuffers;
PT2sInGPU = nullptr;
}

if (hitsInCPU != nullptr) {
delete hitsInCPU;
Expand Down Expand Up @@ -140,6 +147,10 @@ void lst::Event<Acc3D>::resetEvent() {
delete pixelQuintupletsInCPU;
pixelQuintupletsInCPU = nullptr;
}
if (PT2sInCPU != nullptr) {
delete PT2sInCPU;
PT2sInCPU = nullptr;
}
if (trackCandidatesInCPU != nullptr) {
delete trackCandidatesInCPU;
trackCandidatesInCPU = nullptr;
Expand Down Expand Up @@ -730,6 +741,39 @@ void lst::Event<Acc3D>::createTrackCandidates(bool no_pls_dupclean, bool tc_pls_

alpaka::enqueue(queue, checkHitspLSTask);
}
/*
Vec3D const threadsPerBlock_crossCleanPT2{1, 16, 64};
Vec3D const blocksPerGrid_crossCleanPT2{1, 4, 20};
WorkDiv3D const crossClean2_workDiv =
createWorkDiv(blocksPerGrid_crossCleanPT2, threadsPerBlock_crossCleanPT2, elementsPerThread);

lst::crossCleanPT2 crossCleanPT2_kernel;
auto const crossCleanPT2Task(alpaka::createTaskKernel<Acc3D>(crossCleanPT2_workDiv,
crossCleanPT2_kernel,
*modulesBuffers_->data(),
*rangesInGPU,
*PT2sInGPU,
*segmentsInGPU,
*pixelTripletsInGPU,
*pixelQuintupletsInGPU));

alpaka::enqueue(queue, crossCleanPT2Task);
*/
Vec3D const threadsPerBlock_addPT2asTrackCandidatesInGPU{1, 1, 512};
Vec3D const blocksPerGrid_addPT2asTrackCandidatesInGPU{1, 1, 1};
WorkDiv3D const addPT2asTrackCandidatesInGPU_workDiv = createWorkDiv(
blocksPerGrid_addPT2asTrackCandidatesInGPU, threadsPerBlock_addPT2asTrackCandidatesInGPU, elementsPerThread);

lst::addPT2asTrackCandidatesInGPU addPT2asTrackCandidatesInGPU_kernel;
auto const addPT2asTrackCandidatesInGPUTask(alpaka::createTaskKernel<Acc3D>(addPT2asTrackCandidatesInGPU_workDiv,
addPT2asTrackCandidatesInGPU_kernel,
nLowerModules_,
*PT2sInGPU,
*trackCandidatesInGPU,
*segmentsInGPU,
*rangesInGPU));

alpaka::enqueue(queue, addPT2asTrackCandidatesInGPUTask);

Vec3D const threadsPerBlock_crossCleanpLS{1, 16, 32};
Vec3D const blocksPerGrid_crossCleanpLS{1, 4, 20};
Expand Down Expand Up @@ -768,19 +812,22 @@ void lst::Event<Acc3D>::createTrackCandidates(bool no_pls_dupclean, bool tc_pls_
// Check if either n_max_pixel_track_candidates or n_max_nonpixel_track_candidates was reached
auto nTrackCanpT5Host_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);
auto nTrackCanpT3Host_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);
auto nTrackCanPT2Host_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);
auto nTrackCanpLSHost_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);
auto nTrackCanT5Host_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);
alpaka::memcpy(queue, nTrackCanpT5Host_buf, trackCandidatesBuffers->nTrackCandidatespT5_buf);
alpaka::memcpy(queue, nTrackCanpT3Host_buf, trackCandidatesBuffers->nTrackCandidatespT3_buf);
alpaka::memcpy(queue, nTrackCanPT2Host_buf, trackCandidatesBuffers->nTrackCandidatesPT2_buf);
alpaka::memcpy(queue, nTrackCanpLSHost_buf, trackCandidatesBuffers->nTrackCandidatespLS_buf);
alpaka::memcpy(queue, nTrackCanT5Host_buf, trackCandidatesBuffers->nTrackCandidatesT5_buf);
alpaka::wait(queue);

int nTrackCandidatespT5 = *alpaka::getPtrNative(nTrackCanpT5Host_buf);
int nTrackCandidatespT3 = *alpaka::getPtrNative(nTrackCanpT3Host_buf);
int nTrackCandidatesPT2 = *alpaka::getPtrNative(nTrackCanPT2Host_buf);
int nTrackCandidatespLS = *alpaka::getPtrNative(nTrackCanpLSHost_buf);
int nTrackCandidatesT5 = *alpaka::getPtrNative(nTrackCanT5Host_buf);
if ((nTrackCandidatespT5 + nTrackCandidatespT3 + nTrackCandidatespLS == n_max_pixel_track_candidates) ||
if ((nTrackCandidatespT5 + nTrackCandidatespT3 + nTrackCandidatespLS + nTrackCandidatesPT2 == n_max_pixel_track_candidates) ||
(nTrackCandidatesT5 == n_max_nonpixel_track_candidates)) {
printf(
"****************************************************************************************************\n"
Expand Down Expand Up @@ -911,6 +958,125 @@ void lst::Event<Acc3D>::createPixelTriplets() {
alpaka::wait(queue);
}

void lst::Event<Acc3D>::createPT2s() {
if (PT2sInGPU == nullptr) {
PT2sInGPU = new lst::PT2s();
PT2sBuffers = new lst::PT2sBuffer<Device>(n_max_pixel_segments, devAcc, queue);
PT2sInGPU->setData(*PT2sBuffers);
}

unsigned int nInnerSegments;
Copy link

@slava77 slava77 Dec 4, 2024

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

why not call this nPLS or nPixelSegments to match the naming more consistently (this applies to pT3/5 code from where this is copied)

auto nInnerSegments_src_view = alpaka::createView(devHost, &nInnerSegments, (size_t)1u);

auto dev_view_nSegments = alpaka::createSubView(segmentsBuffers->nSegments_buf, (Idx)1u, (Idx)nLowerModules_);
Copy link

@slava77 slava77 Dec 4, 2024

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

this is reading nPLs, right?
It would be better to name accordingly
In some other code pixelModuleIndex is used (which I'm guessing is the same value as nLowerModules_
(this applies to pT3/5 code from where this is copied)


alpaka::memcpy(queue, nInnerSegments_src_view, dev_view_nSegments);
Copy link

@slava77 slava77 Dec 4, 2024

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

memcpy method is (queue, destination, source) please rename nInnerSegments_src_view to nInnerSegments_host_view to be less confusing
(this applies to pT3/5 code from where this is copied)

alpaka::wait(queue);

auto superbins_buf = allocBufWrapper<int>(devHost, n_max_pixel_segments_per_module, queue);
auto pixelTypes_buf = allocBufWrapper<int8_t>(devHost, n_max_pixel_segments_per_module, queue);

alpaka::memcpy(queue, superbins_buf, segmentsBuffers->superbin_buf);
alpaka::memcpy(queue, pixelTypes_buf, segmentsBuffers->pixelType_buf);
alpaka::wait(queue);

auto connectedPixelSize_host_buf = allocBufWrapper<unsigned int>(devHost, nInnerSegments, queue);
auto connectedPixelIndex_host_buf = allocBufWrapper<unsigned int>(devHost, nInnerSegments, queue);
auto connectedPixelSize_dev_buf = allocBufWrapper<unsigned int>(devAcc, nInnerSegments, queue);
auto connectedPixelIndex_dev_buf = allocBufWrapper<unsigned int>(devAcc, nInnerSegments, queue);

int* superbins = alpaka::getPtrNative(superbins_buf);
int8_t* pixelTypes = alpaka::getPtrNative(pixelTypes_buf);
unsigned int* connectedPixelSize_host = alpaka::getPtrNative(connectedPixelSize_host_buf);
unsigned int* connectedPixelIndex_host = alpaka::getPtrNative(connectedPixelIndex_host_buf);
alpaka::wait(queue);

int pixelIndexOffsetPos =
pixelMapping_->connectedPixelsIndex[size_superbins - 1] + pixelMapping_->connectedPixelsSizes[size_superbins - 1];
int pixelIndexOffsetNeg = pixelMapping_->connectedPixelsIndexPos[size_superbins - 1] +
pixelMapping_->connectedPixelsSizesPos[size_superbins - 1] + pixelIndexOffsetPos;

// TODO: check if a map/reduction to just eligible pLSs would speed up the kernel
// the current selection still leaves a significant fraction of unmatchable pLSs
for (unsigned int i = 0; i < nInnerSegments; i++) { // loop over # pLS
int8_t pixelType = pixelTypes[i]; // Get pixel type for this pLS
int superbin = superbins[i]; // Get superbin for this pixel
if ((superbin < 0) or (superbin >= (int)size_superbins) or (pixelType > 2) or (pixelType < 0)) {
connectedPixelSize_host[i] = 0;
connectedPixelIndex_host[i] = 0;
continue;
}

// Used pixel type to select correct size-index arrays
if (pixelType == 0) {
connectedPixelSize_host[i] =
pixelMapping_->connectedPixelsSizes[superbin]; // number of connected modules to this pixel
Comment on lines +1001 to +1013
Copy link

@slava77 slava77 Dec 4, 2024

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

it looks like the same operations on the host side are done for the third time now.
It would be better to carry the connectedPixel* device buffers from the ones filled in the pT5 creation

auto connectedIdxBase = pixelMapping_->connectedPixelsIndex[superbin];
connectedPixelIndex_host[i] =
connectedIdxBase; // index to get start of connected modules for this superbin in map
} else if (pixelType == 1) {
connectedPixelSize_host[i] =
pixelMapping_->connectedPixelsSizesPos[superbin]; // number of pixel connected modules
auto connectedIdxBase = pixelMapping_->connectedPixelsIndexPos[superbin] + pixelIndexOffsetPos;
connectedPixelIndex_host[i] = connectedIdxBase; // index to get start of connected pixel modules
} else if (pixelType == 2) {
connectedPixelSize_host[i] =
pixelMapping_->connectedPixelsSizesNeg[superbin]; // number of pixel connected modules
auto connectedIdxBase = pixelMapping_->connectedPixelsIndexNeg[superbin] + pixelIndexOffsetNeg;
connectedPixelIndex_host[i] = connectedIdxBase; // index to get start of connected pixel modules
}
}

alpaka::memcpy(queue, connectedPixelSize_dev_buf, connectedPixelSize_host_buf, nInnerSegments);
alpaka::memcpy(queue, connectedPixelIndex_dev_buf, connectedPixelIndex_host_buf, nInnerSegments);
alpaka::wait(queue);

Vec3D const threadsPerBlock{1, 4, 32};
Vec3D const blocksPerGrid{16 /* above median of connected modules*/, 4096, 1};
WorkDiv3D const createPT2sInGPUFromMapv2_workDiv =
createWorkDiv(blocksPerGrid, threadsPerBlock, elementsPerThread);

lst::createPT2sInGPUFromMapv2 createPT2sInGPUFromMapv2_kernel;
auto const createPT2sInGPUFromMapv2Task(
alpaka::createTaskKernel<Acc3D>(createPT2sInGPUFromMapv2_workDiv,
createPT2sInGPUFromMapv2_kernel,
*modulesBuffers_->data(),
*rangesInGPU,
*mdsInGPU,
*segmentsInGPU,
*PT2sInGPU,
alpaka::getPtrNative(connectedPixelSize_dev_buf),
alpaka::getPtrNative(connectedPixelIndex_dev_buf),
nInnerSegments,
ptCut));

alpaka::enqueue(queue, createPT2sInGPUFromMapv2Task);
alpaka::wait(queue);

#ifdef WARNINGS
auto nPT2s_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);

alpaka::memcpy(queue, nPT2s_buf, pT2sBuffers->nPT2s_buf);
alpaka::wait(queue);

std::cout << "number of pixel triplets = " << *alpaka::getPtrNative(nPT2s_buf) << std::endl;
#endif

//pT3s can be cleaned here because they're not used in making pT5s!
Vec3D const threadsPerBlockDupPixTrip{1, 16, 16};
//seems like more blocks lead to conflicting writes
Vec3D const blocksPerGridDupPixTrip{1, 40, 1};
WorkDiv3D const removeDupPT2sInGPUFromMap_workDiv =
createWorkDiv(blocksPerGridDupPixTrip, threadsPerBlockDupPixTrip, elementsPerThread);

lst::removeDupPT2sInGPUFromMap removeDupPT2sInGPUFromMap_kernel;
auto const removeDupPT2sInGPUFromMapTask(alpaka::createTaskKernel<Acc3D>(
removeDupPT2sInGPUFromMap_workDiv, removeDupPT2sInGPUFromMap_kernel, *PT2sInGPU));

alpaka::enqueue(queue, removeDupPT2sInGPUFromMapTask);
alpaka::wait(queue);
}

void lst::Event<Acc3D>::createQuintuplets() {
Vec3D const threadsPerBlockCreateQuints{1, 1, 1024};
Vec3D const blocksPerGridCreateQuints{1, 1, 1};
Expand Down Expand Up @@ -1382,6 +1548,17 @@ unsigned int lst::Event<Acc3D>::getNumberOfTripletsByLayerEndcap(unsigned int la
return n_triplets_by_layer_endcap_[layer];
}

int lst::Event<Acc3D>::getNumberOfPT2s() {
auto nPT2s_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);

alpaka::memcpy(queue, nPT2s_buf, PT2sBuffers->nPT2s_buf);
alpaka::wait(queue);

int nPT2s = *alpaka::getPtrNative(nPT2s_buf);

return nPT2s;
}

int lst::Event<Acc3D>::getNumberOfPixelTriplets() {
auto nPixelTriplets_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);

Expand Down Expand Up @@ -1464,6 +1641,17 @@ int lst::Event<Acc3D>::getNumberOfPT3TrackCandidates() {
return nTrackCandidatesPT3;
}

int lst::Event<Acc3D>::getNumberOfPT2TrackCandidates() {
auto nTrackCandidatesPT2_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);

alpaka::memcpy(queue, nTrackCandidatesPT2_buf, trackCandidatesBuffers->nTrackCandidatesPT2_buf);
alpaka::wait(queue);

int nTrackCandidatesPT2 = *alpaka::getPtrNative(nTrackCandidatesPT2_buf);

return nTrackCandidatesPT2;
}

int lst::Event<Acc3D>::getNumberOfPLSTrackCandidates() {
auto nTrackCandidatesPLS_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);

Expand Down Expand Up @@ -1572,6 +1760,10 @@ lst::MiniDoubletsBuffer<DevHost>* lst::Event<Acc3D>::getMiniDoublets() {
alpaka::memcpy(queue, mdsInCPU->dphichanges_buf, miniDoubletsBuffers->dphichanges_buf, nMemHost);
alpaka::memcpy(queue, mdsInCPU->nMDs_buf, miniDoubletsBuffers->nMDs_buf);
alpaka::memcpy(queue, mdsInCPU->totOccupancyMDs_buf, miniDoubletsBuffers->totOccupancyMDs_buf);
alpaka::memcpy(queue, mdsInCPU->anchorX_buf, miniDoubletsBuffers->anchorX_buf);
alpaka::memcpy(queue, mdsInCPU->anchorY_buf, miniDoubletsBuffers->anchorY_buf);
alpaka::memcpy(queue, mdsInCPU->anchorZ_buf, miniDoubletsBuffers->anchorZ_buf);
alpaka::memcpy(queue, mdsInCPU->anchorRt_buf, miniDoubletsBuffers->anchorRt_buf);
alpaka::wait(queue);
}
return mdsInCPU;
Expand Down Expand Up @@ -1687,6 +1879,36 @@ lst::QuintupletsBuffer<DevHost>* lst::Event<Acc3D>::getQuintuplets() {
return quintupletsInCPU;
}

lst::PT2sBuffer<DevHost>* lst::Event<Acc3D>::getPT2s() {
if (PT2sInCPU == nullptr) {
// Get nPT2s parameter to initialize host based PT2sInCPU
auto nPT2s_buf = allocBufWrapper<unsigned int>(devHost, 1, queue);
alpaka::memcpy(queue, nPT2s_buf, PT2sBuffers->nPT2s_buf);
alpaka::wait(queue);

unsigned int nPT2s = *alpaka::getPtrNative(nPT2s_buf);
PT2sInCPU = new lst::PT2sBuffer<DevHost>(nPT2s, devHost, queue);
Copy link

@slava77 slava77 Dec 4, 2024

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
PT2sInCPU = new lst::PT2sBuffer<DevHost>(nPT2s, devHost, queue);
pT2sInCPU = new lst::PT2sBuffer<DevHost>(nPT2s, devHost, queue);

PT2sInCPU->setData(*PT2sInCPU);

*alpaka::getPtrNative(PT2sInCPU->nPT2s_buf) = nPT2s;
alpaka::memcpy(
queue, PT2sInCPU->totOccupancyPT2s_buf, PT2sBuffers->totOccupancyPT2s_buf);
alpaka::memcpy(queue, PT2sInCPU->rzChiSquared_buf, PT2sBuffers->rzChiSquared_buf, nPT2s);
alpaka::memcpy(
queue, PT2sInCPU->segmentIndices_buf, PT2sBuffers->segmentIndices_buf, nPT2s);
alpaka::memcpy(queue,
PT2sInCPU->pixelSegmentIndices_buf,
PT2sBuffers->pixelSegmentIndices_buf,
nPT2s);
alpaka::memcpy(queue, PT2sInCPU->pixelRadius_buf, PT2sBuffers->pixelRadius_buf, nPT2s);
alpaka::memcpy(queue, PT2sInCPU->isDup_buf, PT2sBuffers->isDup_buf, nPT2s);
alpaka::memcpy(queue, PT2sInCPU->eta_buf, PT2sBuffers->eta_buf, nPT2s);
alpaka::memcpy(queue, PT2sInCPU->phi_buf, PT2sBuffers->phi_buf, nPT2s);
alpaka::wait(queue);
}
return PT2sInCPU;
}

lst::PixelTripletsBuffer<DevHost>* lst::Event<Acc3D>::getPixelTriplets() {
if (pixelTripletsInCPU == nullptr) {
// Get nPixelTriplets parameter to initialize host based quintupletsInCPU
Expand Down
8 changes: 8 additions & 0 deletions RecoTracker/LSTCore/src/alpaka/Event.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,6 +13,7 @@
#include "MiniDoublet.h"
#include "PixelQuintuplet.h"
#include "PixelTriplet.h"
#include "PT2.h"
#include "TrackCandidate.h"

#include "HeterogeneousCore/AlpakaInterface/interface/host.h"
Expand Down Expand Up @@ -62,6 +63,8 @@ namespace lst {
QuintupletsBuffer<Device>* quintupletsBuffers;
TrackCandidates* trackCandidatesInGPU;
TrackCandidatesBuffer<Device>* trackCandidatesBuffers;
PT2s* PT2sInGPU;
PT2sBuffer<Device>* PT2sBuffers;
PixelTriplets* pixelTripletsInGPU;
PixelTripletsBuffer<Device>* pixelTripletsBuffers;
PixelQuintuplets* pixelQuintupletsInGPU;
Expand All @@ -76,6 +79,7 @@ namespace lst {
TrackCandidatesBuffer<DevHost>* trackCandidatesInCPU;
ModulesBuffer<DevHost>* modulesInCPU;
QuintupletsBuffer<DevHost>* quintupletsInCPU;
PT2sBuffer<DevHost>* PT2sInCPU;
PixelTripletsBuffer<DevHost>* pixelTripletsInCPU;
PixelQuintupletsBuffer<DevHost>* pixelQuintupletsInCPU;

Expand Down Expand Up @@ -155,6 +159,7 @@ namespace lst {
void createTrackCandidates(bool no_pls_dupclean, bool tc_pls_triplets);
void createExtendedTracks();
void createQuintuplets();
void createPT2s();
void createPixelTriplets();
void createPixelQuintuplets();
void pixelLineSegmentCleaning(bool no_pls_dupclean);
Expand Down Expand Up @@ -183,6 +188,7 @@ namespace lst {
int getNumberOfPixelTrackCandidates();
int getNumberOfPT5TrackCandidates();
int getNumberOfPT3TrackCandidates();
int getNumberOfPT2TrackCandidates();
int getNumberOfT5TrackCandidates();
int getNumberOfPLSTrackCandidates();

Expand All @@ -191,6 +197,7 @@ namespace lst {
unsigned int getNumberOfQuintupletsByLayerBarrel(unsigned int layer);
unsigned int getNumberOfQuintupletsByLayerEndcap(unsigned int layer);

int getNumberOfPT2s();
int getNumberOfPixelTriplets();
int getNumberOfPixelQuintuplets();

Expand All @@ -203,6 +210,7 @@ namespace lst {
QuintupletsBuffer<DevHost>* getQuintuplets();
TrackCandidatesBuffer<DevHost>* getTrackCandidates();
TrackCandidatesBuffer<DevHost>* getTrackCandidatesInCMSSW();
PT2sBuffer<DevHost>* getPT2s();
PixelTripletsBuffer<DevHost>* getPixelTriplets();
PixelQuintupletsBuffer<DevHost>* getPixelQuintuplets();
ModulesBuffer<DevHost>* getModules(bool isFull = false);
Expand Down
Loading