32 unsigned int vtxPlane,
unsigned int vtxXCell,
unsigned int vtxYCell,
57 for(
const auto& pe:
fPE){
68 unsigned int view = plane % 2;
97 if(localplaneid>=-30&&localplaneid<=120&&localcellid>=-70&&localcellid<=70){
98 int index = 141*(localplaneid+30) + (localcellid+70);
99 if(index>=0&&index<(
int)
fPEX.size()){
109 if(localplaneid>=-30&&localplaneid<=120&&localcellid>=-70&&localcellid<=70){
110 int index = 141*(localplaneid+30) + (localcellid+70);
111 if(index>=0&&index<(
int)
fPEY.size()){
122 const unsigned int&
cell)
124 unsigned int view = plane % 2;
137 const unsigned int&
cell)
const 146 const unsigned int&
cell)
148 unsigned int view = plane % 2;
168 for(
unsigned int iCell = 0; iCell <
fNCell; ++iCell)
170 for(
unsigned int iPlane = 0; iPlane <
fNPlane; iPlane += 2)
186 for(
unsigned int iCell = 0; iCell <
fNCell; ++iCell)
188 for(
unsigned int iPlane = 1; iPlane <
fNPlane; iPlane += 2)
215 for(
unsigned int iPlane = 0; iPlane <
fNPlane; ++iPlane)
217 for(
unsigned int iCell = 0; iCell <
fNCell; ++iCell)
220 hist->SetBinContent(iPlane+1, iCell + fNCell*(iPlane%2) + 1,
235 for(
unsigned int iPlane = 0; iPlane <
fNPlane; ++iPlane)
237 for(
unsigned int iCell = 0; iCell <
fNCell; ++iCell)
243 if ( pe > 0 && lab < 1 ) lab = 0.1;
246 hist->SetBinContent(iPlane+1, iCell + fNCell*(iPlane%2) + 1, lab);
259 for(
unsigned int iPlane = 0; iPlane < ((
fNPlane/2)+(
fNPlane%2)); ++iPlane)
261 for(
unsigned int iCell = 0; iCell <
fNCell; ++iCell)
267 hist->SetBinContent(iPlane+1, iCell+1, pe);
280 for(
unsigned int iPlane = 0; iPlane < ((
fNPlane/2)+(
fNPlane%2)); ++iPlane)
282 for(
unsigned int iCell = 0; iCell <
fNCell; ++iCell)
288 if ( pe > 0 && lab < 1 ) lab = 0.1;
291 hist->SetBinContent(iPlane+1, iCell+1, lab);
300 os <<
"RegPixelMap with " << m.
NPixel() <<
" pixels, " 301 << m.
NPlane() <<
" planes" 302 <<
" by " << m.
NCell() <<
" cells" ;
std::vector< HitType > fLabY
Vector of Y Truth labels for pixels.
int FirstCell(const unsigned int &view) const
std::vector< double > fPurX
Vector of X purity for pixels.
void Add(const unsigned int &plane, const unsigned int &cell, const float &pe, const HitType label, const double purity)
unsigned int GlobalToIndex(const unsigned int &plane, const unsigned int &cell)
Take global plane, cell (detector) and return index in fPE vector.
unsigned int fNPlane
Number of planes, length of pixel map.
unsigned int NPixel() const
Total number of pixels in map.
std::ostream & operator<<(std::ostream &os, const PixelMapProducer &p)
void FillInputVector(float *input) const
TH2F * ToTH2(TString name) const
Return the pixel map as a 2D histogram for visualization.
std::vector< float > fPEX
Vector of X PE measurements for pixels.
Defines an enumeration for prong classification.
std::vector< HitType > fLabX
Vector of X Truth labels for pixels.
std::vector< double > fPur
Vector of purity for pixels.
unsigned int fNCell
Number of cells, width of pixel map.
TH2F * SingleViewToLabTH2(const unsigned int &view, TString name) const
unsigned int GlobalToIndexSingle(const unsigned int &plane, const unsigned int &cell)
Take global plane, cell (detector) and return index in fPE vector.
std::vector< HitType > fLab
Vector of Truth labels for pixels.
unsigned int LocalToIndex(const unsigned int &plane, const unsigned int &cell) const
Take local plane, cell (within map) and return index in fPE vector.
unsigned int NCell() const
Width in cells.
std::vector< float > fPEY
Vector of Y PE measurements for pixels.
unsigned int NPlane() const
Length in planes.
RegPixelMap, basic input to CVN neural net.
std::vector< float > fPE
Vector of PE measurements for pixels.
assert(nhit_max >=nhit_nbins)
TH2F * SingleViewToTH2(const unsigned int &view, TString name) const
RegPixelMap for Regression modified from Dominick's PixelMap.
std::vector< double > fPurY
Vector of Y purity for pixels.
TH2F * ToLabTH2(TString name) const
bool IsWithin(const unsigned int &plane, const unsigned int &cell, const unsigned int &view=0)