17 if (tjList.empty())
return;
21 CTP_t inCTP = slc.
tjs[tjList[0][0] - 1].CTP;
22 for (
unsigned short it1 = 0; it1 < slc.
tjs.size(); ++it1) {
24 if (tj1.
CTP != inCTP)
continue;
30 bool inShower =
false;
32 for (
size_t l1 = 0; l1 < tjList.size(); ++l1) {
34 for (
size_t l2 = 0; l2 < tjList[l1].size(); ++l2) {
36 if (trajID == tjList[l1][l2]) {
47 for (
int i = 0; i < 8; i++) {
62 if (ss.
ID == 0)
return;
66 for (
unsigned short it1 = 0; it1 < slc.
tjs.size(); ++it1) {
77 bool isShowerTj =
false;
78 for (
size_t i = 0; i <
stv.
TjID.size(); ++i) {
90 if (isShowerTj)
continue;
93 for (
size_t i = 0; i < ss.
TjIDs.size(); ++i) {
94 if (trajID == ss.
TjIDs[i]) {
127 for (
int i = 0; i < 8; i++) {
182 for (
unsigned short cotIndex = 0; cotIndex < slc.
cots.size(); ++cotIndex) {
183 auto& ss = slc.
cots[cotIndex];
184 if (ss.CTP != inCTP)
continue;
185 if (ss.ID == 0)
continue;
193 for (
unsigned short cotIndex = 0; cotIndex < slc.
cots.size(); ++cotIndex) {
194 auto& ss = slc.
cots[cotIndex];
195 if (ss.ID == 0)
continue;
203 bool existingStage =
false;
204 for (
unsigned short i = 0; i < stv.
StageName.size(); ++i) {
206 existingStage =
true;
211 if (!existingStage) {
std::vector< int > EnvStage
std::vector< int > IsShowerParent
std::vector< Trajectory > tjs
vector of all trajectories in each plane
std::vector< float > EndWir
std::vector< float > EndAng
std::vector< ShowerStruct > cots
std::vector< float > BeginTim
std::vector< Point2_t > Envelope
std::vector< float > BeginAng
The data type to uniquely identify a Plane.
void ClearShowerTree(ShowerTreeVars &stv)
short MCSMom
Normalized RMS using ALL hits. Assume it is 50% to start.
void SaveAllCots(TCSlice &slc, const CTP_t &inCTP, std::string someText)
std::vector< float > EndTim
std::vector< int > ShowerID
auto vector(Vector const &v)
Returns a manipulator which will print the specified array.
std::vector< std::string > StageName
float unitsPerTick
scale factor from Tick to WSE equivalent units
std::vector< short > BeginVtx
CTP_t CTP
Cryostat, TPC, Plane code.
std::vector< TrajPoint > Pts
Trajectory points.
std::vector< short > EndVtx
void SaveTjInfo(TCSlice &slc, std::vector< std::vector< int >> &tjList, std::string stageName)
std::array< unsigned short, 2 > EndPt
First and last point in the trajectory that has charge.
std::vector< float > Envelope
PlaneID_t Plane
Index of the plane within its TPC.
Definition of data types for geometry description.
std::vector< float > BeginChg
int ID
ID that is local to one slice.
std::array< unsigned short, 2 > VtxID
ID of 2D vertex.
std::bitset< 16 > modes
number of points to find AveChg
std::vector< int > EnvPlane
std::vector< short > MCSMom
std::vector< int > StageNum
void SaveTjInfoStuff(TCSlice &slc, Trajectory &tj, int stageNum)
std::bitset< 128 > AlgMod
Bit set if algorithm AlgBit_t modifed the trajectory.
std::vector< float > BeginWir
geo::PlaneID DecodeCTP(CTP_t CTP)
std::vector< float > EndChg
int GetStageNum(ShowerTreeVars &stv, std::string stageName)
std::vector< int > EnvShowerID
std::vector< int > IsShowerTj
std::vector< short > PlaneNum