Skip to content
This repository has been archived by the owner on Dec 9, 2024. It is now read-only.

Add Occupancy Branches Back #404

Merged
merged 3 commits into from
May 29, 2024
Merged

Add Occupancy Branches Back #404

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
f73eef7
add occ branches back
GNiendorf May 21, 2024
7749db7
small cleanup
GNiendorf May 21, 2024
5794ed1
add if statement back, fixed
GNiendorf May 28, 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
87 changes: 87 additions & 0 deletions code/core/write_sdl_ntuple.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -136,6 +136,23 @@ void createOptionalOutputBranches()
ana.tx->createBranch<vector<float>>("t5_rzChiSquared");
ana.tx->createBranch<vector<float>>("t5_nonAnchorChiSquared");

// Occupancy branches
ana.tx->createBranch<vector<int>>("module_layers");
ana.tx->createBranch<vector<int>>("module_subdets");
ana.tx->createBranch<vector<int>>("module_rings");
ana.tx->createBranch<vector<int>>("module_rods");
ana.tx->createBranch<vector<int>>("module_modules");
ana.tx->createBranch<vector<bool>>("module_isTilted");
ana.tx->createBranch<vector<float>>("module_eta");
ana.tx->createBranch<vector<float>>("module_r");
ana.tx->createBranch<vector<int>>("md_occupancies");
ana.tx->createBranch<vector<int>>("sg_occupancies");
ana.tx->createBranch<vector<int>>("t3_occupancies");
ana.tx->createBranch<int>("tc_occupancies");
ana.tx->createBranch<vector<int>>("t5_occupancies");
ana.tx->createBranch<int>("pT3_occupancies");
ana.tx->createBranch<int>("pT5_occupancies");

#endif
}

Expand Down Expand Up @@ -297,10 +314,80 @@ void setOptionalOutputBranches(SDL::Event<SDL::Acc>* event)
setPixelQuintupletOutputBranches(event);
setQuintupletOutputBranches(event);
setPixelTripletOutputBranches(event);
setOccupancyBranches(event);

#endif
}

//________________________________________________________________________________________________________________________________
void setOccupancyBranches(SDL::Event<SDL::Acc>* event)
{
SDL::trackCandidates& trackCandidatesInGPU = (*event->getTrackCandidates());
SDL::triplets& tripletsInGPU = (*event->getTriplets());
SDL::segments& segmentsInGPU = (*event->getSegments());
SDL::miniDoublets& mdsInGPU = (*event->getMiniDoublets());
//SDL::hits& hitsInGPU = (*event->getHits());
SDL::modules& modulesInGPU = (*event->getModules());
SDL::quintuplets& quintupletsInGPU = (*event->getQuintuplets());
SDL::pixelQuintuplets& pixelQuintupletsInGPU = (*event->getPixelQuintuplets());
SDL::pixelTriplets& pixelTripletsInGPU = (*event->getPixelTriplets());
//get the occupancies from these dudes
std::vector<int> moduleLayer;
std::vector<int> moduleSubdet;
std::vector<int> moduleRing;
std::vector<int> moduleRod;
std::vector<int> moduleModule;
std::vector<float> moduleEta;
std::vector<float> moduleR;
std::vector<bool> moduleIsTilted;
std::vector<int> trackCandidateOccupancy;
std::vector<int> tripletOccupancy;
std::vector<int> segmentOccupancy;
std::vector<int> mdOccupancy;
std::vector<int> quintupletOccupancy;

bool isTilted;

for(unsigned int idx = 0; idx <= *(modulesInGPU.nLowerModules); idx++)
{
//layer = 0, subdet = 0 => pixel module
//module, md and segment - need some gymnastics
unsigned int lowerIdx = idx;//modulesInGPU.lowerModuleIndices[idx];
moduleLayer.push_back(modulesInGPU.layers[lowerIdx]);
moduleSubdet.push_back(modulesInGPU.subdets[lowerIdx]);
moduleRing.push_back(modulesInGPU.rings[lowerIdx]);
moduleRod.push_back(modulesInGPU.rods[lowerIdx]);
moduleEta.push_back(modulesInGPU.eta[lowerIdx]);
moduleR.push_back(modulesInGPU.r[lowerIdx]);
isTilted = (modulesInGPU.subdets[lowerIdx] == 5 and modulesInGPU.sides[lowerIdx] != 3);
moduleIsTilted.push_back(isTilted);
moduleModule.push_back(modulesInGPU.modules[lowerIdx]);
segmentOccupancy.push_back(segmentsInGPU.totOccupancySegments[lowerIdx]);
mdOccupancy.push_back(mdsInGPU.totOccupancyMDs[lowerIdx]);

if(idx < *(modulesInGPU.nLowerModules))
GNiendorf marked this conversation as resolved.
Show resolved Hide resolved
{
quintupletOccupancy.push_back(quintupletsInGPU.totOccupancyQuintuplets[idx]);
tripletOccupancy.push_back(tripletsInGPU.totOccupancyTriplets[idx]);
}
}
ana.tx->setBranch<vector<int>>("module_layers",moduleLayer);
ana.tx->setBranch<vector<int>>("module_subdets",moduleSubdet);
ana.tx->setBranch<vector<int>>("module_rings",moduleRing);
ana.tx->setBranch<vector<int>>("module_rods",moduleRod);
ana.tx->setBranch<vector<int>>("module_modules",moduleModule);
ana.tx->setBranch<vector<bool>>("module_isTilted",moduleIsTilted);
ana.tx->setBranch<vector<float>>("module_eta",moduleEta);
ana.tx->setBranch<vector<float>>("module_r",moduleR);
ana.tx->setBranch<vector<int>>("md_occupancies",mdOccupancy);
ana.tx->setBranch<vector<int>>("sg_occupancies",segmentOccupancy);
ana.tx->setBranch<vector<int>>("t3_occupancies",tripletOccupancy);
ana.tx->setBranch<int>("tc_occupancies",*(trackCandidatesInGPU.nTrackCandidates));
ana.tx->setBranch<int>("pT3_occupancies", *(pixelTripletsInGPU.totOccupancyPixelTriplets));
ana.tx->setBranch<vector<int>>("t5_occupancies", quintupletOccupancy);
ana.tx->setBranch<int>("pT5_occupancies", *(pixelQuintupletsInGPU.totOccupancyPixelQuintuplets));
}

//________________________________________________________________________________________________________________________________
void setPixelQuintupletOutputBranches(SDL::Event<SDL::Acc>* event)
{
Expand Down
1 change: 1 addition & 0 deletions code/core/write_sdl_ntuple.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,7 @@ void createGnnNtupleBranches();
void fillOutputBranches(SDL::Event<SDL::Acc>* event);
void setOutputBranches(SDL::Event<SDL::Acc>* event);
void setOptionalOutputBranches(SDL::Event<SDL::Acc>* event);
void setOccupancyBranches(SDL::Event<SDL::Acc>* event);
void setPixelQuintupletOutputBranches(SDL::Event<SDL::Acc>* event);
void setQuintupletOutputBranches(SDL::Event<SDL::Acc>* event);
void setPixelTripletOutputBranches(SDL::Event<SDL::Acc> *event);
Expand Down