Public Member Functions | Protected Member Functions | Private Member Functions | List of all members
evd::RecoBaseDrawer Class Reference

Aid in the rendering of RecoBase objects. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-15/EventDisplay/RecoBaseDrawer.h"

Public Member Functions

 RecoBaseDrawer ()
 
 ~RecoBaseDrawer ()
 
void CellHit2D (art::Event const &evt, evdb::View2D *xzview, evdb::View2D *yzview)
 
void Vertex2D (art::Event const &evt, evdb::View2D *xzview, evdb::View2D *yzview)
 
void Cluster2D (art::Event const &evt, evdb::View2D *xzview, evdb::View2D *yzview)
 
void OfflineChans2D (const art::Event &evt, evdb::View2D *xzview, evdb::View2D *yzview)
 
void HoughResult2D (art::Event const &evt, evdb::View2D *xzview, evdb::View2D *yzview)
 
void Prong2D (art::Event const &evt, evdb::View2D *xzview, evdb::View2D *yzview)
 
void Track2D (const art::Event &evt, evdb::View2D *xzview, evdb::View2D *yzview)
 
void CellHit3D (art::Event const &evt, evdb::View3D *view)
 
void Cluster3D (art::Event const &evt, evdb::View3D *view)
 
void Prong3D (art::Event const &evt, evdb::View3D *view)
 
void FillTQHisto (art::Event const &evt, TH1F *thisto, TH1F *qhisto)
 

Protected Member Functions

void DrawHoughResults2D (std::vector< const rb::HoughResult * > &htr, evdb::View2D *xzview, evdb::View2D *yzview)
 
void DrawVertices2D (std::vector< const rb::Vertex * > &vert, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &vertexIndex)
 
void DrawCluster2D (std::vector< const rb::Cluster * > &clusts, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &clustIndex)
 
void DrawProng2D (std::vector< const rb::Prong * > &prong, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &prongIndex)
 
void DrawTrack2D (std::vector< const rb::Track * > &track, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &trackIndex)
 
void DrawProng3D (std::vector< const rb::Prong * > &prong, evdb::View3D *view, std::vector< int > &prongIndex)
 
void DrawTrack3D (std::vector< const rb::Track * > &track, evdb::View3D *view, std::vector< int > &trackIndex)
 

Private Member Functions

void GetClusterOutlines (const rb::Cluster *c, int view, std::vector< double > &vzVpts, std::vector< double > &vzZpts)
 
void GetClusterOutlinesConvexHull (const rb::Cluster *c, int viewint, std::vector< double > &vzVpts, std::vector< double > &vzZpts) const
 
int GetClusterColor (unsigned int idx) const
 
int DimOutOfTime (double t1, double t2, int c)
 
int DimFiltered (const rb::CellHit *chit, std::vector< const rb::CellHit * > unfiltered, int c)
 
double ScaleFactor (double, double)
 

Detailed Description

Aid in the rendering of RecoBase objects.

Definition at line 27 of file RecoBaseDrawer.h.

Constructor & Destructor Documentation

evd::RecoBaseDrawer::RecoBaseDrawer ( )

Definition at line 44 of file RecoBaseDrawer.cxx.

45  {
46  }
evd::RecoBaseDrawer::~RecoBaseDrawer ( )

Definition at line 50 of file RecoBaseDrawer.cxx.

51  {
52 
53  }

Member Function Documentation

void evd::RecoBaseDrawer::CellHit2D ( art::Event const &  evt,
evdb::View2D *  xzview,
evdb::View2D *  yzview 
)

Render CellHit objects on a 2D viewing canvas

Parameters
evt: Event handle to get data objects from
which: Which collection of object to render
xzview: Pointer to X-Z view to draw on
yzview: Pointer to Y-Z view to draw on

Definition at line 65 of file RecoBaseDrawer.cxx.

References evdb::View2D::AddBox(), plot_validation_datamc::c, rb::CellHit::Cell(), colors, DimFiltered(), evd::SliceNavigator::DimHits(), DimOutOfTime(), dz, evd::RecoDrawingOptions::fCellHitInstances, evd::RecoDrawingOptions::fCellHitModules, evd::RawDrawingOptions::fColor, evd::RawDrawingOptions::fWhichHits, evdb::ColorScale::GetColor(), evd::SliceNavigator::GetProducts(), make_syst_table_plots::h, MECModelEnuComparisons::i, evd::RawDrawingOptions::kCAL_HITS, evd::RawDrawingOptions::kCOLOR_BY_CHARGE, evd::RawDrawingOptions::kCOLOR_BY_TIME, geo::kX, geo::kY, getGoodRuns4SAM::n, nc, rb::CellHit::PE(), rb::CellHit::Plane(), geo::GeometryBase::Plane(), evdb::Colors::Scale(), ScaleFactor(), evd::SliceNavigator::SlicesToNavigate(), string, confusionMatrixTree::t, rb::CellHit::TNS(), registry_explorer::v, and rb::CellHit::View().

68  {
69  art::ServiceHandle<evd::RawDrawingOptions> rawopt;
70 
71  if (rawopt->fWhichHits!=evd::RawDrawingOptions::kCAL_HITS)
72  return;
73 
74  art::ServiceHandle<evd::SliceNavigator> nav;
75 
76  art::ServiceHandle<evd::RecoDrawingOptions> recoopt;
77  art::ServiceHandle<geo::Geometry> geo;
78 
79  std::vector<const rb::Cluster*> slicesToNavigate = nav->SlicesToNavigate(evt);
80  std::vector<const rb::CellHit*> unfiltHits;
81  for(int i=0, n = slicesToNavigate.size(); i < n; ++i ){
82  for( int c = 0, nc = slicesToNavigate[i]->NCell(); c < nc; ++c)
83  unfiltHits.push_back( slicesToNavigate[i]->Cell(c).get());
84  }
85 
86  //std::cout << recoopt->fCellHitModules.size() << " modules found with cell hits." << std::endl;
87 
88  for (unsigned int imod=0; imod<recoopt->fCellHitModules.size(); ++imod) {
89  const std::string& which = recoopt->fCellHitModules[imod].c_str();
90  const std::string& instance = recoopt->fCellHitInstances[imod];
91 
92  std::vector<const rb::CellHit*> cellhits;
93  nav->GetProducts(evt, which, instance, cellhits);
94 
95  // Display all hits on the two 2D views provided
96  for (unsigned int i=0; i<cellhits.size(); ++i) {
97  const rb::CellHit* h = cellhits[i];
98  evdb::View2D* v = 0;
99 
100  // Select which view to draw this cell hit on
101  if (h->View() == geo::kX) v = xzview;
102  if (h->View() == geo::kY) v = yzview;
103  if (v==0) continue;
104 
105  // Try to get the "best" charge measurement, ie. the one last in
106  // the calibration chain
107  float q = h->PE();
108  float t = (h->TNS())/1000.0; // thisto is binned in usec
109 
110  // Render cell hits as markers with size and color encoding time
111  // and charge
112  double dz = geo->Plane(h->Plane())->Cell(h->Cell())->HalfD();
113  double dv = geo->Plane(h->Plane())->Cell(h->Cell())->HalfW();
114  double zpos=0;
115  double vpos=0;
116 
117  double xyz[3];
118  geo->Plane(h->Plane())->Cell(h->Cell())->GetCenter(xyz);
119  zpos = xyz[2];
120  if (h->View() == geo::kX) vpos = xyz[0];
121  else vpos = xyz[1];
122 
123  double sf = this->ScaleFactor(t, q);
124  //double q0 = 400.0;
125  int c = 1;
126  art::ServiceHandle<evdb::Colors> colors;
127  switch (rawopt->fColor) {
128  case evd::RawDrawingOptions::kCOLOR_BY_CHARGE:
129  c = colors->Scale("CalQ").GetColor(q);
130  break;
131  case evd::RawDrawingOptions::kCOLOR_BY_TIME:
132  c = colors->Scale("CalT").GetColor(t);
133  //sf = sqrt(q/q0);
134  //if (sf>1.0) sf = 1.0;
135  break;
136  default:
137  c = 1;
138  }
139  c = this->DimOutOfTime(t, t, c);
140  if(nav->DimHits() == true) c = this->DimFiltered( h, unfiltHits, c );
141 
142  TBox& b1 = v->AddBox(zpos-sf*dz, vpos-sf*dv, zpos+sf*dz, vpos+sf*dv);
143  b1.SetFillStyle(1001);
144  b1.SetFillColor(c);
145  b1.SetBit(kCannotPick);
146  }
147  } // loop on imod folders
148  }
rb::CellHit::TNS
float TNS() const
Definition: CellHit.h:46
rb::CellHit::Plane
unsigned short Plane() const
Definition: CellHit.h:39
rb::CellHit::View
geo::View_t View() const
Definition: CellHit.h:41
evd::RecoDrawingOptions::fCellHitModules
std::vector< std::string > fCellHitModules
Module labels for CellHits.
Definition: RecoDrawingOptions.h:55
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
evd::RecoBaseDrawer::ScaleFactor
double ScaleFactor(double, double)
Definition: RecoBaseDrawer.cxx:248
evdb::ColorScale::GetColor
int GetColor(double x) const
Definition: ColorScale.cxx:126
geo::GeometryBase::Plane
const PlaneGeo * Plane(unsigned int i) const
Definition: GeometryBase.cxx:749
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
geo::kY
Horizontal planes which measure Y.
Definition: PlaneGeo.h:29
rb::CellHit::Cell
unsigned short Cell() const
Definition: CellHit.h:40
evd::RawDrawingOptions::kCOLOR_BY_TIME
static const int kCOLOR_BY_TIME
Definition: RawDrawingOptions.h:35
evdb::View2D::AddBox
TBox & AddBox(double x1, double y1, double x2, double y2)
Definition: View2D.cxx:263
colors
int colors[6]
Definition: tools.h:1
dz
double dz[NP][NC]
Definition: checkMeasured_Target.C:20
evd::RawDrawingOptions::kCAL_HITS
static const int kCAL_HITS
Definition: RawDrawingOptions.h:39
evd::SliceNavigator::DimHits
bool DimHits() const
Definition: SliceNavigator.h:71
evd::RecoDrawingOptions::fCellHitInstances
std::vector< std::string > fCellHitInstances
productInstance labels for CellHits
Definition: RecoDrawingOptions.h:56
make_syst_table_plots.h
h
Definition: make_syst_table_plots.py:240
evdb::Colors::Scale
ColorScale & Scale(const std::string &nm)
Definition: Colors_service.cc:90
evt
int evt
Definition: FillTrainingTrees.h:22
plot_validation_datamc.c
c
Definition: plot_validation_datamc.py:198
evd::RecoBaseDrawer::DimFiltered
int DimFiltered(const rb::CellHit *chit, std::vector< const rb::CellHit * > unfiltered, int c)
Definition: RecoBaseDrawer.cxx:1193
rb::CellHit::PE
float PE() const
Definition: CellHit.h:42
evd::RecoBaseDrawer::DimOutOfTime
int DimOutOfTime(double t1, double t2, int c)
Definition: RecoBaseDrawer.cxx:1178
evd::RawDrawingOptions::kCOLOR_BY_CHARGE
static const int kCOLOR_BY_CHARGE
Definition: RawDrawingOptions.h:34
rb::CellHit
A rawdata::RawDigit with channel information decoded.
Definition: CellHit.h:27
nc
enum BeamMode nc
registry_explorer.v
v
Definition: registry_explorer.py:22
geo
Helper for AttenCurve.
Definition: Path.h:10
evd::SliceNavigator::SlicesToNavigate
std::vector< const rb::Cluster * > SlicesToNavigate(const art::Event &evt)
Definition: SliceNavigator_service.cc:669
string
enum BeamMode string
void evd::RecoBaseDrawer::CellHit3D ( art::Event const &  evt,
evdb::View3D *  view 
)

Render CellHit objects on a 3D viewing canvas

Parameters
evt: Event handle to get data objects from
which: Which collection of object to render
view: Pointer to 3D viewer

Definition at line 159 of file RecoBaseDrawer.cxx.

References evdb::View3D::AddPolyLine3D(), plot_validation_datamc::c, rb::CellHit::Cell(), DimOutOfTime(), evd::RecoDrawingOptions::fCellHitInstances, evd::RecoDrawingOptions::fCellHitModules, evd::RawDrawingOptions::fWhichHits, evd::SliceNavigator::GetProducts(), make_syst_table_plots::h, MECModelEnuComparisons::i, evd::RawDrawingOptions::kCAL_HITS, geo::kX, rb::CellHit::PE(), rb::CellHit::Plane(), geo::GeometryBase::Plane(), std::sqrt(), string, confusionMatrixTree::t, rb::CellHit::TNS(), rb::CellHit::View(), x1, submit_syst::x2, and y1.

161  {
162  art::ServiceHandle<evd::RawDrawingOptions> rawopt;
163  if (rawopt->fWhichHits!=evd::RawDrawingOptions::kCAL_HITS) return;
164  art::ServiceHandle<evd::SliceNavigator> nav;
165 
166  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
167  art::ServiceHandle<geo::Geometry> geo;
168 
169  for (unsigned int imod=0; imod<drawopt->fCellHitModules.size(); ++imod) {
170  const char* which = drawopt->fCellHitModules[imod].c_str();
171  const std::string& instance = drawopt->fCellHitInstances[imod];
172 
173  std::vector<const rb::CellHit*> cellhits;
174  nav->GetProducts(evt, which, instance, cellhits);
175 
176  // Display all hits on the two 2D views provided
177  for (unsigned int i=0; i<cellhits.size(); ++i) {
178  const rb::CellHit* h = cellhits[i];
179 
180  // Try to get the "best" charge measurement, ie. the one last in
181  // the calibration chain
182  float q=h->PE();
183  // if (h->IsA()==rb::RecoHit::Class()) {
184  // q = ((rb::RecoHit*)h)->GeV(); // or MIP or PECorr()
185  // }
186 
187  float t=h->TNS();
188  // if (h->IsA()==rb::RecoHit::Class()) {
189  // t = ((rb::RecoHit*)h)->T();
190  // }
191 
192  // Render cell hits as markers with size and color encoding time
193  // and charge
194  static evdb::ColorScale tscale(0.0,30.0);
195  int c = tscale(t);
196  double q0 = 1.0;
197  double sz = 0.3*sqrt(q/q0);
198 
199  c = this->DimOutOfTime(t,t,c);
200 
201  TPolyLine3D& p = view->AddPolyLine3D(5,c,1,2);
202  double x1, x2, y1, z1, z2;
203  // double y2; // i get set but am not used
204 
205  double xyz[3];
206  geo->Plane(h->Plane())->Cell(h->Cell())->GetCenter(xyz);
207 
208  // if (h->IsA()==rb::RecoHit::Class()) {
209  // z1 = ((rb::RecoHit*)h)->Z();
210  // if (h->View() == geo::kX) {
211  // x1 = ((rb::RecoHit*)h)->V();
212  // y1 = xyz[1];
213  // }
214  // else {
215  // y1 = ((rb::RecoHit*)h)->V();
216  // x1 = xyz[1];
217  // }
218  // }
219  // else {
220  z1 = xyz[2];
221  if (h->View() == geo::kX) {
222  x1 = xyz[0];
223  y1 = xyz[1];
224  }
225  else {
226  x1 = xyz[1];
227  y1 = xyz[0];
228  }
229  // }
230 
231  x2 = x1+sz;
232  //y2 = y1+sz; // i am not used anywhere
233  z2 = z1+sz;
234  x1 -= sz;
235  y1 -= sz;
236  z1 -= sz;
237 
238  p.SetPoint(0,x1-sz,y1,z1);
239  p.SetPoint(1,x1+sz,y1,z1);
240  p.SetPoint(2,x2,y1,z2);
241  p.SetPoint(1,x1,y1,z2);
242  p.SetPoint(0,x1,y1,z1);
243 
244  } // loop on i cell hits
245  } // loop on imod folders
246  }
rb::CellHit::TNS
float TNS() const
Definition: CellHit.h:46
y1
Float_t y1[n_points_granero]
Definition: compare.C:5
rb::CellHit::Plane
unsigned short Plane() const
Definition: CellHit.h:39
x1
Float_t x1[n_points_granero]
Definition: compare.C:5
rb::CellHit::View
geo::View_t View() const
Definition: CellHit.h:41
p
const char * p
Definition: xmltok.h:285
std::sqrt
T sqrt(T number)
Definition: d0nt_math.hpp:156
evd::RecoDrawingOptions::fCellHitModules
std::vector< std::string > fCellHitModules
Module labels for CellHits.
Definition: RecoDrawingOptions.h:55
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
submit_syst.x2
x2
Definition: submit_syst.py:250
geo::GeometryBase::Plane
const PlaneGeo * Plane(unsigned int i) const
Definition: GeometryBase.cxx:749
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
rb::CellHit::Cell
unsigned short Cell() const
Definition: CellHit.h:40
evdb::View3D::AddPolyLine3D
TPolyLine3D & AddPolyLine3D(int n, int c, int w, int s)
Definition: View3D.cxx:105
evd::RawDrawingOptions::kCAL_HITS
static const int kCAL_HITS
Definition: RawDrawingOptions.h:39
evdb::ColorScale
Build an association between a numerical range and a ROOT color index for use in, eg...
Definition: ColorScale.h:44
evd::RecoDrawingOptions::fCellHitInstances
std::vector< std::string > fCellHitInstances
productInstance labels for CellHits
Definition: RecoDrawingOptions.h:56
make_syst_table_plots.h
h
Definition: make_syst_table_plots.py:240
evt
int evt
Definition: FillTrainingTrees.h:22
plot_validation_datamc.c
c
Definition: plot_validation_datamc.py:198
rb::CellHit::PE
float PE() const
Definition: CellHit.h:42
evd::RecoBaseDrawer::DimOutOfTime
int DimOutOfTime(double t1, double t2, int c)
Definition: RecoBaseDrawer.cxx:1178
rb::CellHit
A rawdata::RawDigit with channel information decoded.
Definition: CellHit.h:27
geo
Helper for AttenCurve.
Definition: Path.h:10
string
enum BeamMode string
void evd::RecoBaseDrawer::Cluster2D ( art::Event const &  evt,
evdb::View2D *  xzview,
evdb::View2D *  yzview 
)

Definition at line 396 of file RecoBaseDrawer.cxx.

References DrawCluster2D(), evd::RecoDrawingOptions::fClusterIndex, evd::RecoDrawingOptions::fClusterInstances, evd::RecoDrawingOptions::fClusterModules, evd::RecoDrawingOptions::fClusterStyle, evd::SliceNavigator::GetProducts(), and string.

399  {
400  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
401  if(drawopt->fClusterStyle == 0) return;
402 
403  art::ServiceHandle<evd::SliceNavigator> nav;
404 
405  for (unsigned int imod=0; imod<drawopt->fClusterModules.size(); ++imod) {
406  const std::string& which = drawopt->fClusterModules[imod];
407  const std::string& instance = drawopt->fClusterInstances[imod];
408 
409  std::vector<const rb::Cluster*> clust;
410  nav->GetProducts(evt, which, instance, clust);
411 
412  DrawCluster2D(clust, xzview, yzview, drawopt->fClusterIndex[imod]);
413 
414  } // loop on imod folders
415  }
evd::RecoDrawingOptions::fClusterModules
std::vector< std::string > fClusterModules
Module labels for Clusters.
Definition: RecoDrawingOptions.h:58
evd::RecoBaseDrawer::DrawCluster2D
void DrawCluster2D(std::vector< const rb::Cluster * > &clusts, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &clustIndex)
Definition: RecoBaseDrawer.cxx:262
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
evd::RecoDrawingOptions::fClusterInstances
std::vector< std::string > fClusterInstances
productInstance labels for Clusters
Definition: RecoDrawingOptions.h:59
evt
int evt
Definition: FillTrainingTrees.h:22
evd::RecoDrawingOptions::fClusterIndex
std::vector< std::vector< int > > fClusterIndex
Definition: RecoDrawingOptions.h:79
string
enum BeamMode string
void evd::RecoBaseDrawer::Cluster3D ( art::Event const &  evt,
evdb::View3D *  view 
)

Definition at line 419 of file RecoBaseDrawer.cxx.

421  {
422 
423  }
int evd::RecoBaseDrawer::DimFiltered ( const rb::CellHit *  chit,
std::vector< const rb::CellHit * >  unfiltered,
int  c 
)
private

Definition at line 1193 of file RecoBaseDrawer.cxx.

References plot_validation_datamc::c.

Referenced by CellHit2D().

1196  {
1197  if (std::find(unfilt.begin(),unfilt.end(),chit)==unfilt.end()) {
1198  return -1;
1199  };
1200  return c;
1201  }
plot_validation_datamc.c
c
Definition: plot_validation_datamc.py:198
int evd::RecoBaseDrawer::DimOutOfTime ( double  t1,
double  t2,
int  c 
)
private

Definition at line 1178 of file RecoBaseDrawer.cxx.

References colors, lem_server::cs, e, evdb::ColorScale::InBounds(), and evdb::Colors::Scale().

Referenced by CellHit2D(), CellHit3D(), DrawCluster2D(), DrawHoughResults2D(), DrawProng2D(), DrawProng3D(), DrawTrack2D(), DrawTrack3D(), and DrawVertices2D().

1179  {
1180  static const double ns2usec = 1e-3;
1181  t1 *= ns2usec;
1182  t2 *= ns2usec;
1183 
1184  art::ServiceHandle<evdb::Colors> colors;
1185  const evdb::ColorScale& cs = colors->Scale("RawT");
1186  if (cs.InBounds(t1) && cs.InBounds(t2)) return c;
1187 
1188  return -1;
1189  }
evdb::ColorScale::InBounds
bool InBounds(double x) const
Definition: ColorScale.cxx:113
colors
int colors[6]
Definition: tools.h:1
evdb::ColorScale
Build an association between a numerical range and a ROOT color index for use in, eg...
Definition: ColorScale.h:44
evdb::Colors::Scale
ColorScale & Scale(const std::string &nm)
Definition: Colors_service.cc:90
plot_validation_datamc.c
c
Definition: plot_validation_datamc.py:198
t2
double t2
Definition: cana_time_peak.C:50
lem_server.cs
cs
Definition: lem_server.py:582
e
Float_t e
Definition: plot.C:35
void evd::RecoBaseDrawer::DrawCluster2D ( std::vector< const rb::Cluster * > &  clusts,
evdb::View2D *  xzview,
evdb::View2D *  yzview,
std::vector< int > &  clustIndex 
)
protected

Definition at line 262 of file RecoBaseDrawer.cxx.

References evdb::View2D::AddPolyLine(), evdb::View2D::AddPolyMarker(), plot_validation_datamc::c, rb::CellHit::Cell(), plot_validation_datamc_2018::color, DimOutOfTime(), evd::RecoDrawingOptions::fClusterStyle, GetClusterColor(), GetClusterOutlines(), GetClusterOutlinesConvexHull(), make_syst_table_plots::h, MECModelEnuComparisons::i, allTimeWatchdog::index, makeTrainCVSamples::int, rb::Cluster::Is2D(), rb::Cluster::IsNoise(), evd::RecoDrawingOptions::kCLUSTERS_AS_BOXES, evd::RecoDrawingOptions::kCLUSTERS_AS_HULLS, evd::RecoDrawingOptions::kCLUSTERS_AS_MARKERS, geo::kX, geo::kY, rb::Cluster::MaxTNS(), rb::Cluster::MinTNS(), rb::Cluster::NXCell(), rb::Cluster::NYCell(), plot_validation_datamc::p1, plot_validation_datamc::p2, rb::CellHit::Plane(), geo::GeometryBase::Plane(), rb::Cluster::XCell(), and rb::Cluster::YCell().

Referenced by Cluster2D(), Prong2D(), and Track2D().

266  {
267  art::ServiceHandle<geo::Geometry> geo;
268  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
269 
270  bool drawAsMarkers = false;
271  bool drawAsBoxes = false;
272  bool drawAsHulls = false;
273  if (drawopt->fClusterStyle & evd::RecoDrawingOptions::kCLUSTERS_AS_MARKERS) {
274  drawAsMarkers = true;
275  }
276  if (drawopt->fClusterStyle & evd::RecoDrawingOptions::kCLUSTERS_AS_BOXES) {
277  drawAsBoxes = true;
278  }
279  if (drawopt->fClusterStyle & evd::RecoDrawingOptions::kCLUSTERS_AS_HULLS) {
280  drawAsHulls = true;
281  }
282 
283 
284  for (unsigned int iClust=0; iClust<clusts.size(); ++iClust) {
285  const rb::Cluster* c = clusts[iClust];
286 
287  std::vector< int > clustIndex;
288  bool drawThisCluster = clustIndex.empty();
289 
290  // if the indices have been provided, check whether or
291  // not this cluster should be drawn
292  for(int index = 0; index < (int)clustIndex.size(); index++){
293  if( clustIndex[index] == -1 || clustIndex[index] == (int)iClust)
294  drawThisCluster = true;
295  }
296 
297  // Start color index further down for 2D prongs
298  unsigned int coloridx = 0;
299  if(c->Is2D()) coloridx = 5;
300 
301  // Don't render "noise" clusters or clusters you don't want to draw
302  if (c->IsNoise() || !drawThisCluster ) continue;
303 
304  if (drawAsMarkers) {
305  // draw cluster with unique marker
306  int color = this->GetClusterColor(iClust+coloridx);
307  color = this->DimOutOfTime(c->MinTNS(),c->MaxTNS(),color);
308 
309  if(xzview){
310  TPolyMarker& pm = xzview->AddPolyMarker(c->NXCell(),
311  color,
312  kFullCircle, 0.8);
313  for(unsigned int hitIdx = 0; hitIdx < c->NXCell(); ++hitIdx){
314  const art::Ptr<rb::CellHit>& h = c->XCell(hitIdx);
315  double xyz[3];
316  geo->Plane(h->Plane())->Cell(h->Cell())->GetCenter(xyz);
317  pm.SetPoint(hitIdx, xyz[2], xyz[0]);
318  }
319  } // end if xzview
320  if(yzview){
321  TPolyMarker& pm = yzview->AddPolyMarker(c->NYCell(),
322  color,
323  kFullCircle, 0.8);
324  for(unsigned int hitIdx = 0; hitIdx < c->NYCell(); ++hitIdx){
325  const art::Ptr<rb::CellHit>& h = c->YCell(hitIdx);
326  double xyz[3];
327  geo->Plane(h->Plane())->Cell(h->Cell())->GetCenter(xyz);
328  pm.SetPoint(hitIdx, xyz[2], xyz[1]);
329  }
330  } // end if yzview
331  }
332 
333  if(drawAsBoxes || drawAsHulls){
334  std::vector<double> xzZpts, xzXpts;
335  std::vector<double> yzZpts, yzYpts;
336 
337  if (drawAsBoxes) {
338  this->GetClusterOutlines(c, geo::kX, xzXpts, xzZpts);
339  this->GetClusterOutlines(c, geo::kY, yzYpts, yzZpts);
340  }
341  else if (drawAsHulls) {
342  this->GetClusterOutlinesConvexHull(c, geo::kX, xzXpts, xzZpts);
343  this->GetClusterOutlinesConvexHull(c, geo::kY, yzYpts, yzZpts);
344  }
345 
346  bool isX = !xzZpts.empty();
347  bool isY = !yzZpts.empty();
348  int lcolor = this->GetClusterColor(iClust);
349  int fcolor = (isX && isY) ? lcolor : 17; // fill color
350  int width = 1; // line width
351  int style = 1; // 1=solid line style
352 
353  lcolor = this->DimOutOfTime(c->MinTNS(), c->MaxTNS(), lcolor);
354  fcolor = this->DimOutOfTime(c->MinTNS(), c->MaxTNS(), fcolor);
355 
356  if (isX && xzview) {
357  TPolyLine& p1 = xzview->AddPolyLine(xzZpts.size(),
358  lcolor,
359  width,
360  style);
361  TPolyLine& p2 = xzview->AddPolyLine(xzZpts.size(),
362  lcolor,
363  width,
364  style);
365  p1.SetOption("f");
366  p1.SetFillStyle(3003);
367  p1.SetFillColor(fcolor);
368  for (unsigned int i=0; i<xzZpts.size(); ++i) {
369  p1.SetPoint(i, xzZpts[i], xzXpts[i]);
370  p2.SetPoint(i, xzZpts[i], xzXpts[i]);
371  } // loop on i points in ZX view
372  } // if we have a cluster in the ZX view
373  if (isY && yzview) {
374  TPolyLine& p1 = yzview->AddPolyLine(yzZpts.size(),
375  lcolor,
376  width,
377  style);
378  TPolyLine& p2 = yzview->AddPolyLine(yzZpts.size(),
379  lcolor,
380  width,
381  style);
382  p1.SetOption("f");
383  p1.SetFillStyle(3003);
384  p1.SetFillColor(fcolor);
385  for (unsigned int i=0; i<yzZpts.size(); ++i) {
386  p1.SetPoint(i, yzZpts[i], yzYpts[i]);
387  p2.SetPoint(i, yzZpts[i], yzYpts[i]);
388  } // loop on i points in ZY view
389  } // if we have a cluster in the ZY view
390  }
391  } // loop on iClust clusters
392  }
plot_validation_datamc.p2
p2
Definition: plot_validation_datamc.py:199
rb::Cluster::XCell
art::Ptr< rb::CellHit > XCell(unsigned int xIdx) const
Get the ith cell in the x-view.
Definition: Cluster.cxx:157
evdb::View2D::AddPolyLine
TPolyLine & AddPolyLine(int n, int c, int w, int s)
Definition: View2D.cxx:210
rb::Cluster::Is2D
bool Is2D() const
Definition: Cluster.h:96
rb::CellHit::Plane
unsigned short Plane() const
Definition: CellHit.h:39
plot_validation_datamc_2018.color
color
Definition: plot_validation_datamc_2018.py:516
evd::RecoDrawingOptions::kCLUSTERS_AS_HULLS
static const int kCLUSTERS_AS_HULLS
Definition: RecoDrawingOptions.h:36
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
rb::Cluster
A collection of associated CellHits.
Definition: Cluster.h:47
geo::GeometryBase::Plane
const PlaneGeo * Plane(unsigned int i) const
Definition: GeometryBase.cxx:749
geo::kY
Horizontal planes which measure Y.
Definition: PlaneGeo.h:29
evd::RecoDrawingOptions::kCLUSTERS_AS_BOXES
static const int kCLUSTERS_AS_BOXES
Definition: RecoDrawingOptions.h:34
rb::CellHit::Cell
unsigned short Cell() const
Definition: CellHit.h:40
rb::Cluster::MinTNS
double MinTNS() const
Definition: Cluster.cxx:482
evd::RecoDrawingOptions::kCLUSTERS_AS_MARKERS
static const int kCLUSTERS_AS_MARKERS
Definition: RecoDrawingOptions.h:35
rb::Cluster::YCell
art::Ptr< rb::CellHit > YCell(unsigned int yIdx) const
Get the ith cell in the y-view.
Definition: Cluster.cxx:165
make_syst_table_plots.h
h
Definition: make_syst_table_plots.py:240
evdb::View2D::AddPolyMarker
TPolyMarker & AddPolyMarker(int n, int c, int st, double sz)
Definition: View2D.cxx:157
plot_validation_datamc.c
c
Definition: plot_validation_datamc.py:198
evd::RecoBaseDrawer::DimOutOfTime
int DimOutOfTime(double t1, double t2, int c)
Definition: RecoBaseDrawer.cxx:1178
evd::RecoBaseDrawer::GetClusterOutlines
void GetClusterOutlines(const rb::Cluster *c, int view, std::vector< double > &vzVpts, std::vector< double > &vzZpts)
Definition: RecoBaseDrawer.cxx:474
rb::Cluster::MaxTNS
double MaxTNS() const
Definition: Cluster.cxx:528
rb::Cluster::NYCell
unsigned int NYCell() const
Number of cells in the y-view.
Definition: Cluster.h:108
evd::RecoBaseDrawer::GetClusterOutlinesConvexHull
void GetClusterOutlinesConvexHull(const rb::Cluster *c, int viewint, std::vector< double > &vzVpts, std::vector< double > &vzZpts) const
Definition: RecoBaseDrawer.cxx:545
rb::Cluster::NXCell
unsigned int NXCell() const
Number of cells in the x-view.
Definition: Cluster.h:106
rb::Cluster::IsNoise
bool IsNoise() const
Is the noise flag set?
Definition: Cluster.h:163
evd::RecoBaseDrawer::GetClusterColor
int GetClusterColor(unsigned int idx) const
Definition: RecoBaseDrawer.cxx:1152
geo
Helper for AttenCurve.
Definition: Path.h:10
plot_validation_datamc.p1
p1
Definition: plot_validation_datamc.py:199
void evd::RecoBaseDrawer::DrawHoughResults2D ( std::vector< const rb::HoughResult * > &  htr,
evdb::View2D *  xzview,
evdb::View2D *  yzview 
)
protected

Definition at line 644 of file RecoBaseDrawer.cxx.

References evdb::View2D::AddLine(), b, plot_validation_datamc::c, colors, geo::GeometryBase::DetHalfHeight(), geo::GeometryBase::DetHalfWidth(), geo::GeometryBase::DetLength(), DimOutOfTime(), evd::RecoDrawingOptions::fHoughOpt, geom(), GetClusterColor(), evdb::ColorScale::GetColor(), MECModelEnuComparisons::i, calib::j, evd::RecoDrawingOptions::kHOUGH_COLORS, geo::kX, geo::kY, submit_hadd::l, m, evdb::Colors::Scale(), x1, and submit_syst::x2.

Referenced by GetClusterOutlinesConvexHull(), and HoughResult2D().

647  {
648  art::ServiceHandle<evd::RawDrawingOptions> rawopt;
649  art::ServiceHandle<evd::RecoDrawingOptions> recoopt;
650  art::ServiceHandle<geo::Geometry> geom;
651  art::ServiceHandle<evdb::Colors> colors;
652 
653  bool colorByIndex =
654  (recoopt->fHoughOpt == evd::RecoDrawingOptions::kHOUGH_COLORS);
655 
656  for (unsigned int i=0; i<htr.size(); ++i) {
657  double tave = 0.5*(htr[i]->fTNSlo+htr[i]->fTNShi)/1000.0;
658  int c = colors->Scale("RawT").GetColor(tave);
659 
660  for (unsigned int j=0; j<htr[i]->fPeak.size(); ++j) {
661  if (colorByIndex) c = this->GetClusterColor(j);
662  c = this->DimOutOfTime(htr[i]->fTNSlo, htr[i]->fTNShi, c);
663 
664  //
665  // Extract line segment parameters and intersect with detector sides
666  //
667  double m, b;
668  htr[i]->SlopeIcept(j, &m, &b);
669 
670  double x1=0, x2=0;
671  if (htr[i]->fView == geo::kX) {
672  x1 = geom->DetHalfWidth();
673  x2 = -geom->DetHalfWidth();
674  }
675  else {
676  x1 = geom->DetHalfHeight();
677  x2 = -geom->DetHalfHeight();
678  }
679 
680  double z1, z2;
681  z1 = (x1-b)/m;
682  z2 = (x2-b)/m;
683 
684  if (z1<0.0) { z1 = 0.0; x1 = m*z1+b;}
685  if (z1>geom->DetLength()) { z1 = geom->DetLength(); x1 = m*z1+b;}
686  if (z2<0.0) { z2 = 0.0; x2 = m*z2+b;}
687  if (z2>geom->DetLength()) { z2 = geom->DetLength(); x2 = m*z2+b;}
688 
689  if (xzv && htr[i]->fView==geo::kX) {
690  TLine& l = xzv->AddLine(z1, x1, z2, x2);
691  l.SetLineColor(c);
692  l.SetLineStyle(1);
693  }
694 
695  if (yzv && htr[i]->fView==geo::kY) {
696  TLine& l = yzv->AddLine(z1, x1, z2, x2);
697  l.SetLineColor(c);
698  l.SetLineStyle(1);
699  }
700  }
701  }
702  }
x1
Float_t x1[n_points_granero]
Definition: compare.C:5
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
submit_syst.x2
x2
Definition: submit_syst.py:250
geo::GeometryBase::DetLength
double DetLength() const
Definition: GeometryBase.cxx:481
evdb::ColorScale::GetColor
int GetColor(double x) const
Definition: ColorScale.cxx:126
geo::kY
Horizontal planes which measure Y.
Definition: PlaneGeo.h:29
m
#define m
Definition: MakeGENIEReweightTable_module.cc:407
colors
int colors[6]
Definition: tools.h:1
evdb::Colors::Scale
ColorScale & Scale(const std::string &nm)
Definition: Colors_service.cc:90
plot_validation_datamc.c
c
Definition: plot_validation_datamc.py:198
calib::j
const double j
Definition: BetheBloch.cxx:29
evd::RecoBaseDrawer::DimOutOfTime
int DimOutOfTime(double t1, double t2, int c)
Definition: RecoBaseDrawer.cxx:1178
geo::GeometryBase::DetHalfHeight
double DetHalfHeight() const
Definition: GeometryBase.cxx:474
submit_hadd.l
l
Definition: submit_hadd.py:89
geo::GeometryBase::DetHalfWidth
double DetHalfWidth() const
Definition: GeometryBase.cxx:467
geom
void geom(int which=0)
Definition: geom.C:163
b
const hit & b
Definition: hits.cxx:21
evd::RecoDrawingOptions::kHOUGH_COLORS
static const int kHOUGH_COLORS
Definition: RecoDrawingOptions.h:53
evd::RecoBaseDrawer::GetClusterColor
int GetClusterColor(unsigned int idx) const
Definition: RecoBaseDrawer.cxx:1152
void evd::RecoBaseDrawer::DrawProng2D ( std::vector< const rb::Prong * > &  prong,
evdb::View2D *  xzview,
evdb::View2D *  yzview,
std::vector< int > &  prongIndex 
)
protected

Definition at line 746 of file RecoBaseDrawer.cxx.

References evdb::View2D::AddLine(), plot_validation_datamc_2018::color, DimOutOfTime(), dx, dy, dz, GetClusterColor(), MECModelEnuComparisons::i, allTimeWatchdog::index, makeTrainCVSamples::int, geo::kX, and geo::kY.

Referenced by Prong2D(), and Track2D().

750  {
751  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
752  art::ServiceHandle<geo::Geometry> geo;
753  art::ServiceHandle<evd::RawDrawingOptions> rawopt;
754 
755  // Render prongs at 5 meters by default
756  double ds_default = 500.0;
757 
758  for (unsigned int i=0; i<prong.size(); ++i) {
759 
760 
761  bool drawThisProng = false;
762  if( prongIndex.empty() )
763  drawThisProng = true;
764 
765  for(int index = 0; index < (int)prongIndex.size(); index++){
766  if( prongIndex[index] == -1 || prongIndex[index] == (int)i)
767  drawThisProng = true;
768  }
769 
770  if (xzview && prong[i]->View() != geo::kY && drawThisProng) {
771  const double x0 = prong[i]->Start().X();
772  const double z0 = prong[i]->Start().Z();
773  const double dx = prong[i]->Dir().X();
774  const double dz = prong[i]->Dir().Z();
775  double ds = prong[i]->TotalLength();
776  if (ds==0.0) ds = ds_default;
777 
778  int color = this->GetClusterColor(i);
779  color =
780  this->DimOutOfTime(prong[i]->MinTNS(), prong[i]->MaxTNS(),color);
781 
782  TLine& pzx = xzview->AddLine(z0, x0, z0+ds*dz, x0+ds*dx);
783  pzx.SetLineColor(color);
784  pzx.SetLineWidth(2);
785  pzx.SetLineStyle(2);
786  }
787  if (yzview && prong[i]->View() != geo::kX && drawThisProng) {
788  const double y0 = prong[i]->Start().Y();
789  const double z0 = prong[i]->Start().Z();
790  const double dy = prong[i]->Dir().Y();
791  const double dz = prong[i]->Dir().Z();
792  double ds = prong[i]->TotalLength();
793  if (ds==0.0) ds = ds_default;
794 
795  int color = this->GetClusterColor(i);
796  color = this->DimOutOfTime(prong[i]->MinTNS(),
797  prong[i]->MaxTNS(),
798  color);
799 
800  TLine& pzy = yzview->AddLine(z0, y0, z0+ds*dz, y0+ds*dy);
801  pzy.SetLineColor(color);
802  pzy.SetLineWidth(2);
803  pzy.SetLineStyle(2);
804  }
805  }
806  }
plot_validation_datamc_2018.color
color
Definition: plot_validation_datamc_2018.py:516
evdb::View2D::AddLine
TLine & AddLine(double x1, double y1, double x2, double y2)
Definition: View2D.cxx:187
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
geo::kY
Horizontal planes which measure Y.
Definition: PlaneGeo.h:29
dy
double dy[NP][NC]
Definition: checkMeasured_Target.C:19
dx
double dx[NP][NC]
Definition: checkMeasured_Target.C:18
dz
double dz[NP][NC]
Definition: checkMeasured_Target.C:20
evd::RecoBaseDrawer::DimOutOfTime
int DimOutOfTime(double t1, double t2, int c)
Definition: RecoBaseDrawer.cxx:1178
View
Definition: View.py:1
ds
A very simple service to remember what detector we&#39;re working in.
Definition: Detector_service.cc:13
evd::RecoBaseDrawer::GetClusterColor
int GetClusterColor(unsigned int idx) const
Definition: RecoBaseDrawer.cxx:1152
geo
Helper for AttenCurve.
Definition: Path.h:10
void evd::RecoBaseDrawer::DrawProng3D ( std::vector< const rb::Prong * > &  prong,
evdb::View3D *  view,
std::vector< int > &  prongIndex 
)
protected

Definition at line 879 of file RecoBaseDrawer.cxx.

References evdb::View3D::AddPolyLine3D(), plot_validation_datamc_2018::color, DimOutOfTime(), dx, dy, dz, GetClusterColor(), MECModelEnuComparisons::i, allTimeWatchdog::index, and makeTrainCVSamples::int.

Referenced by Prong3D().

882  {
883 
884  for (unsigned int i=0; i<prong.size(); ++i) {
885 
886 
887  bool drawThisProng = false;
888  if( prongIndex.empty() )
889  drawThisProng = true;
890 
891  for(int index = 0; index < (int)prongIndex.size(); index++){
892  if( prongIndex[index] == -1 || prongIndex[index] == (int)i)
893  drawThisProng = true;
894  }
895 
896  if(!drawThisProng )
897  continue;
898 
899  int color = this->GetClusterColor(i);
900  color = this->DimOutOfTime(prong[i]->MinTNS(),
901  prong[i]->MaxTNS(),
902  color);
903 
904  TPolyLine3D& p = view->AddPolyLine3D(2,color,1,1);
905 
906  const double x0 = prong[i]->Start().X();
907  const double y0 = prong[i]->Start().Y();
908  const double z0 = prong[i]->Start().Z();
909 
910  const double dx = prong[i]->Dir().X();
911  const double dy = prong[i]->Dir().Y();
912  const double dz = prong[i]->Dir().Z();
913 
914  p.SetPoint(0, x0, y0, z0);
915  // Prong line is always 5 m long
916  p.SetPoint(1, x0+500*dx, y0+500*dy, z0+500*dz);
917  }
918 
919  }
plot_validation_datamc_2018.color
color
Definition: plot_validation_datamc_2018.py:516
p
const char * p
Definition: xmltok.h:285
evdb::View3D::AddPolyLine3D
TPolyLine3D & AddPolyLine3D(int n, int c, int w, int s)
Definition: View3D.cxx:105
dy
double dy[NP][NC]
Definition: checkMeasured_Target.C:19
dx
double dx[NP][NC]
Definition: checkMeasured_Target.C:18
dz
double dz[NP][NC]
Definition: checkMeasured_Target.C:20
evd::RecoBaseDrawer::DimOutOfTime
int DimOutOfTime(double t1, double t2, int c)
Definition: RecoBaseDrawer.cxx:1178
evd::RecoBaseDrawer::GetClusterColor
int GetClusterColor(unsigned int idx) const
Definition: RecoBaseDrawer.cxx:1152
void evd::RecoBaseDrawer::DrawTrack2D ( std::vector< const rb::Track * > &  track,
evdb::View2D *  xzview,
evdb::View2D *  yzview,
std::vector< int > &  trackIndex 
)
protected

Definition at line 810 of file RecoBaseDrawer.cxx.

References evdb::View2D::AddMarker(), evdb::View2D::AddPolyLine(), rb::CellHit::Cell(), rb::Cluster::Cell(), cellMax, col, DimOutOfTime(), evd::RecoDrawingOptions::fTrackOpt, geom(), art::Ptr< T >::get(), GetClusterColor(), allTimeWatchdog::index, makeTrainCVSamples::int, evd::RecoDrawingOptions::kORTHO_HITS, geo::kX, geo::kY, rb::Cluster::MaxTNS(), rb::Cluster::MinTNS(), getGoodRuns4SAM::n, rb::Cluster::NCell(), rb::Track::NTrajectoryPoints(), rb::CellHit::Plane(), geo::GeometryBase::Plane(), rb::Track::TrajectoryPoint(), rb::CellHit::View(), rb::Cluster::View(), w, rb::Track::W(), submit_syst::x, submit_syst::y, and test::z.

Referenced by Track2D().

814  {
815  art::ServiceHandle<geo::Geometry> geom;
816  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
817 
818  bool orthoHits =
819  (drawopt->fTrackOpt & evd::RecoDrawingOptions::kORTHO_HITS);
820 
821  const size_t trkMax = trks.size();
822  for(size_t trkIdx = 0; trkIdx < trkMax; ++trkIdx){
823  const rb::Track* trk = trks[trkIdx];
824 
825  bool drawThisTrack = false;
826  if( trackIndex.empty() )
827  drawThisTrack = true;
828 
829  for(int index = 0; index < (int)trackIndex.size(); index++){
830  if( trackIndex[index] == -1 || trackIndex[index] == (int)trkIdx)
831  drawThisTrack = true;
832  }
833 
834  if(trk->NTrajectoryPoints() < 2 || !drawThisTrack ) continue;
835 
836  int col = GetClusterColor(trkIdx);
837  col = this->DimOutOfTime(trk->MinTNS(), trk->MaxTNS(), col);
838 
839  TPolyLine* px = 0;
840  TPolyLine* py = 0;
841  if(xzview && trk->View() != geo::kY)
842  px = &xzview->AddPolyLine(trk->NTrajectoryPoints(), col, 2, 2);
843  if(yzview && trk->View() != geo::kX)
844  py = &yzview->AddPolyLine(trk->NTrajectoryPoints(), col, 2, 2);
845 
846  for(size_t n = 0; n < trk->NTrajectoryPoints(); ++n){
847  const double x = trk->TrajectoryPoint(n).X();
848  const double y = trk->TrajectoryPoint(n).Y();
849  const double z = trk->TrajectoryPoint(n).Z();
850 
851  if(px) px->SetPoint(n, z, x);
852  if(py) py->SetPoint(n, z, y);
853  } // end for n
854 
855  if (orthoHits) {
856  const unsigned int cellMax = trk->NCell();
857  for(unsigned int cellIdx = 0; cellIdx < cellMax; ++cellIdx){
858  art::Ptr<rb::CellHit> chit = trk->Cell(cellIdx);
859 
860  const double w = trk->W(chit.get());
861 
862  double xyz[3];
863  geom->Plane(chit->Plane())->Cell(chit->Cell())->GetCenter(xyz);
864  const double z = xyz[2];
865 
866  if(chit->View() == geo::kX && yzview){
867  yzview->AddMarker(z, w, col, kFullCircle, .75);
868  }
869  if(chit->View() == geo::kY && xzview){
870  xzview->AddMarker(z, w, col, kFullCircle, .75);
871  }
872  } // end for cellIdx
873  } // end if
874  } // end for trkIdx
875  }
rb::Cluster::View
virtual geo::View_t View() const
kXorY for 3D clusters.
Definition: Cluster.h:99
rb::Track::NTrajectoryPoints
size_t NTrajectoryPoints() const
Definition: Track.h:83
rb::Cluster::NCell
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
Definition: Cluster.cxx:134
rb::Track::TrajectoryPoint
TVector3 TrajectoryPoint(unsigned int i) const
The ith point on the trajectory, a 3-vector in cm.
Definition: Track.cxx:158
evd::RecoDrawingOptions::kORTHO_HITS
static const int kORTHO_HITS
Definition: RecoDrawingOptions.h:44
rb::Track::W
virtual double W(const rb::CellHit *chit) const
Estimate the unmeasured coordinate of chit.
Definition: Track.cxx:281
evdb::View2D::AddPolyLine
TPolyLine & AddPolyLine(int n, int c, int w, int s)
Definition: View2D.cxx:210
rb::CellHit::Plane
unsigned short Plane() const
Definition: CellHit.h:39
rb::CellHit::View
geo::View_t View() const
Definition: CellHit.h:41
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
rb::Track
A rb::Prong with full reconstructed trajectory.
Definition: Track.h:20
geo::GeometryBase::Plane
const PlaneGeo * Plane(unsigned int i) const
Definition: GeometryBase.cxx:749
geo::kY
Horizontal planes which measure Y.
Definition: PlaneGeo.h:29
submit_syst.y
y
Definition: submit_syst.py:249
rb::CellHit::Cell
unsigned short Cell() const
Definition: CellHit.h:40
trk
Track finder for cosmic rays.
Definition: CosmicTrackAna_module.cc:182
rb::Cluster::MinTNS
double MinTNS() const
Definition: Cluster.cxx:482
submit_syst.x
x
Definition: submit_syst.py:249
col
Int_t col[ntarg]
Definition: Style.C:29
evd::RecoBaseDrawer::DimOutOfTime
int DimOutOfTime(double t1, double t2, int c)
Definition: RecoBaseDrawer.cxx:1178
test.z
z
Definition: test.py:28
rb::Cluster::MaxTNS
double MaxTNS() const
Definition: Cluster.cxx:528
rb::Cluster::Cell
art::Ptr< rb::CellHit > Cell(geo::View_t view, unsigned int viewIdx) const
Get the ith cell from view view.
Definition: Cluster.cxx:145
art::Ptr::get
T const * get() const
Definition: Ptr.h:321
cellMax
const int cellMax
Definition: fit_thresh_corrs.C:30
geom
void geom(int which=0)
Definition: geom.C:163
evd::RecoBaseDrawer::GetClusterColor
int GetClusterColor(unsigned int idx) const
Definition: RecoBaseDrawer.cxx:1152
evdb::View2D::AddMarker
TMarker & AddMarker(double x, double y, int c, int st, double sz)
Definition: View2D.cxx:124
w
Float_t w
Definition: plot.C:20
void evd::RecoBaseDrawer::DrawTrack3D ( std::vector< const rb::Track * > &  track,
evdb::View3D *  view,
std::vector< int > &  trackIndex 
)
protected

Definition at line 923 of file RecoBaseDrawer.cxx.

References evdb::View3D::AddPolyLine3D(), plot_validation_datamc_2018::color, DimOutOfTime(), GetClusterColor(), MECModelEnuComparisons::i, allTimeWatchdog::index, makeTrainCVSamples::int, getGoodRuns4SAM::n, submit_syst::x, submit_syst::y, and test::z.

Referenced by Prong3D().

926  {
927 
928  for(unsigned int i = 0; i < track.size(); ++i){
929 
930  bool drawThisTrack = false;
931  if( trackIndex.empty() )
932  drawThisTrack = true;
933 
934  for(int index = 0; index < (int)trackIndex.size(); index++){
935  if( trackIndex[index] == -1 || trackIndex[index] == (int)i)
936  drawThisTrack = true;
937  }
938 
939  if(track[i]->NTrajectoryPoints() < 2 || !drawThisTrack) continue;
940 
941  int color = this->GetClusterColor(i);
942  color = this->DimOutOfTime(track[i]->MinTNS(),
943  track[i]->MaxTNS(),
944  color);
945  TPolyLine3D& p = view->AddPolyLine3D(2,color,1,1);
946  for(unsigned int n = 0; n < track[i]->NTrajectoryPoints(); ++n){
947  const double x = track[i]->TrajectoryPoint(n).X();
948  const double y = track[i]->TrajectoryPoint(n).Y();
949  const double z = track[i]->TrajectoryPoint(n).Z();
950 
951  p.SetPoint(n, x , y , z);
952  }
953  }
954  }
plot_validation_datamc_2018.color
color
Definition: plot_validation_datamc_2018.py:516
p
const char * p
Definition: xmltok.h:285
track
Definition: event.h:19
submit_syst.y
y
Definition: submit_syst.py:249
evdb::View3D::AddPolyLine3D
TPolyLine3D & AddPolyLine3D(int n, int c, int w, int s)
Definition: View3D.cxx:105
submit_syst.x
x
Definition: submit_syst.py:249
evd::RecoBaseDrawer::DimOutOfTime
int DimOutOfTime(double t1, double t2, int c)
Definition: RecoBaseDrawer.cxx:1178
test.z
z
Definition: test.py:28
evd::RecoBaseDrawer::GetClusterColor
int GetClusterColor(unsigned int idx) const
Definition: RecoBaseDrawer.cxx:1152
void evd::RecoBaseDrawer::DrawVertices2D ( std::vector< const rb::Vertex * > &  vert,
evdb::View2D *  xzview,
evdb::View2D *  yzview,
std::vector< int > &  vertexIndex 
)
protected

Definition at line 706 of file RecoBaseDrawer.cxx.

References evdb::View2D::AddMarker(), col, DimOutOfTime(), GetClusterColor(), MECModelEnuComparisons::i, allTimeWatchdog::index, makeTrainCVSamples::int, confusionMatrixTree::t, submit_syst::x, submit_syst::y, and test::z.

710  {
711 
712  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
713  art::ServiceHandle<geo::Geometry> geo;
714 
715  const int mstyle = 28; //set marker style (cross)
716  const double msize = 2.0; //set marker size
717 
718  for (unsigned int i=0; i<vert.size(); ++i) {
719 
720  bool drawThisVertex = false;
721  if( vertexIndex.empty() )
722  drawThisVertex = true;
723 
724  for(int index = 0; index < (int)vertexIndex.size(); index++){
725  if( vertexIndex[index] == -1 || vertexIndex[index] == (int)i)
726  drawThisVertex = true;
727  }
728 
729  if(!drawThisVertex) continue;
730 
731  const double x = vert[i]->GetX();
732  const double y = vert[i]->GetY();
733  const double z = vert[i]->GetZ();
734  const double t = vert[i]->GetT();
735 
736  int col = GetClusterColor(i+6);// Start this color at i=6
737  col = this->DimOutOfTime(t, t, col);
738 
739  if (xzview) xzview->AddMarker(z,x,col,mstyle,msize);
740  if (yzview) yzview->AddMarker(z,y,col,mstyle,msize);
741  }
742  }
submit_syst.y
y
Definition: submit_syst.py:249
submit_syst.x
x
Definition: submit_syst.py:249
col
Int_t col[ntarg]
Definition: Style.C:29
evd::RecoBaseDrawer::DimOutOfTime
int DimOutOfTime(double t1, double t2, int c)
Definition: RecoBaseDrawer.cxx:1178
test.z
z
Definition: test.py:28
evd::RecoBaseDrawer::GetClusterColor
int GetClusterColor(unsigned int idx) const
Definition: RecoBaseDrawer.cxx:1152
geo
Helper for AttenCurve.
Definition: Path.h:10
evdb::View2D::AddMarker
TMarker & AddMarker(double x, double y, int c, int st, double sz)
Definition: View2D.cxx:124
void evd::RecoBaseDrawer::FillTQHisto ( art::Event const &  evt,
TH1F *  thisto,
TH1F *  qhisto 
)

Definition at line 1205 of file RecoBaseDrawer.cxx.

References evd::RecoDrawingOptions::fCellHitInstances, evd::RecoDrawingOptions::fCellHitModules, evd::SliceNavigator::GetProducts(), make_syst_table_plots::h, MECModelEnuComparisons::i, rb::CellHit::PE(), string, confusionMatrixTree::t, and rb::CellHit::TNS().

Referenced by evd::TQPad::Draw().

1208  {
1209  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
1210  art::ServiceHandle<evd::SliceNavigator> nav;
1211 
1212  for (unsigned int imod=0; imod<drawopt->fCellHitModules.size(); ++imod) {
1213  std::string const& which = drawopt->fCellHitModules[imod].c_str();
1214  std::string const& instance = drawopt->fCellHitInstances[imod].c_str();
1215 
1216  std::vector<const rb::CellHit*> cellhits;
1217  nav->GetProducts(evt, which, instance, cellhits);
1218 
1219  for (unsigned int i=0; i<cellhits.size(); ++i) {
1220  const rb::CellHit* h = cellhits[i];
1221 
1222  // best available charge, time
1223  float q=h->PE();
1224  // if (h->IsA()==rb::RecoHit::Class()) {
1225  // q = ((rb::RecoHit*)h)->PECorr(); // or GeV or MIP?
1226  // }
1227  float t=(h->TNS())/1000.0; // thisto is binned in usec
1228  // if (h->IsA()==rb::RecoHit::Class()) {
1229  // t = ((rb::RecoHit*)h)->T();
1230  // }
1231 
1232  if (thisto!=0) thisto->Fill(t);
1233  if (qhisto!=0) qhisto->Fill(q);
1234  }
1235  }
1236 
1237  }
rb::CellHit::TNS
float TNS() const
Definition: CellHit.h:46
evd::RecoDrawingOptions::fCellHitModules
std::vector< std::string > fCellHitModules
Module labels for CellHits.
Definition: RecoDrawingOptions.h:55
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
evd::RecoDrawingOptions::fCellHitInstances
std::vector< std::string > fCellHitInstances
productInstance labels for CellHits
Definition: RecoDrawingOptions.h:56
make_syst_table_plots.h
h
Definition: make_syst_table_plots.py:240
evt
int evt
Definition: FillTrainingTrees.h:22
rb::CellHit::PE
float PE() const
Definition: CellHit.h:42
rb::CellHit
A rawdata::RawDigit with channel information decoded.
Definition: CellHit.h:27
string
enum BeamMode string
int evd::RecoBaseDrawer::GetClusterColor ( unsigned int  idx) const
private

Definition at line 1152 of file RecoBaseDrawer.cxx.

References colors, evdb::Colors::Foreground(), kBlue, kGreen, kOrange, kRed, and kViolet.

Referenced by DrawCluster2D(), DrawHoughResults2D(), DrawProng2D(), DrawProng3D(), DrawTrack2D(), DrawTrack3D(), DrawVertices2D(), and OfflineChans2D().

1153  {
1154  art::ServiceHandle<evdb::Colors> colors;
1155  bool isBonW = (colors->Foreground(0) == kBlack);
1156 
1157  // Background dependent shades
1158  const int kNCol = 12;
1159  const int kCols[kNCol] ={
1160  (isBonW) ? kRed-4 : kAzure+6,
1161  (isBonW) ? kGreen-3 : kPink+1,
1162  (isBonW) ? kBlue-7 : kGreen-9,
1163  (isBonW) ? kOrange+1 : kOrange+1,
1164  (isBonW) ? kViolet-1 : kViolet-5,
1165  (isBonW) ? kTeal+1 : kTeal+1,
1166  (isBonW) ? kOrange-3 : kYellow-9,
1167  (isBonW) ? kCyan+1 : kCyan+1,
1168  (isBonW) ? kMagenta+1 : kRed-7,
1169  (isBonW) ? kSpring+1 : kSpring+10,
1170  (isBonW) ? kAzure+1 : kAzure+6,
1171  (isBonW) ? kPink+1 : kPink+1
1172  };
1173  return kCols[idx%kNCol];
1174  }
kOrange
enum BeamMode kOrange
kRed
enum BeamMode kRed
compare_h5_caf.idx
int idx
Definition: compare_h5_caf.py:69
colors
int colors[6]
Definition: tools.h:1
evdb::Colors::Foreground
int Foreground(int i=0)
Definition: Colors_service.cc:127
kViolet
enum BeamMode kViolet
kGreen
enum BeamMode kGreen
kBlue
enum BeamMode kBlue
void evd::RecoBaseDrawer::GetClusterOutlines ( const rb::Cluster *  c,
int  viewint,
std::vector< double > &  vzVpts,
std::vector< double > &  vzZpts 
)
private

Make a set of points which outline a cluster

Parameters
c: Reco base cluster to outline
view: Detector view to consider
vzXpts: X or Y coordinates of points in XZ or YZ view that outline cluster
vzZpts: Z coordinates of points in XZ ir YZ view that outline cluster

Definition at line 474 of file RecoBaseDrawer.cxx.

References getBrightness::cell, rb::CellHit::Cell(), rb::Cluster::Cell(), geo::GeometryBase::CellInfo(), GetClusterOutlinesConvexHull(), calib::j, rb::Cluster::NCell(), rb::CellHit::Plane(), NDAPDHVSetting::plane, elec2geo::pos, and POTSpillRate::view.

Referenced by DrawCluster2D().

478  {
479  art::ServiceHandle<geo::Geometry> geo;
480 
481  geo::View_t view = geo::View_t(viewint);
482 
483  const int ncell = c->NCell(view);
484  if(ncell == 0) return;
485 
486  // Map plane numbers to highest and lowest cells in the plane
487  std::map<int,int> celllo, cellhi;
488  // On first pass, initialize
489  for (int j=0; j<ncell; ++j) {
490  const int plane = c->Cell(view, j)->Plane();
491  const int cell = c->Cell(view, j)->Cell();
492  celllo[plane] = cell;
493  cellhi[plane] = cell;
494  }
495  // Finalize on second pass
496  for (int j=0; j<ncell; ++j) {
497  const int plane = c->Cell(view, j)->Plane();
498  const int cell = c->Cell(view, j)->Cell();
499  if (cell < celllo[plane]) celllo[plane] = cell;
500  if (cell > cellhi[plane]) cellhi[plane] = cell;
501  }
502 
503  // Loop over planes and low cells to make lines along bottom
504  // edge. Work from upstream edge to downstream edge
505  std::map<int,int>::iterator itr (celllo.begin());
506  std::map<int,int>::iterator itrEnd(celllo.end());
507  for (; itr!=itrEnd; ++itr) {
508  const int plane = itr->first;
509  const int cell = itr->second;
510  double pos[3], dpos[3];
511  geo->CellInfo(plane, cell, &view, pos, dpos);
512 
513  vzZpts.push_back(pos[2]-dpos[2]);
514  vzVpts.push_back(pos[view]-dpos[view]);
515 
516  vzZpts.push_back(pos[2]+dpos[2]);
517  vzVpts.push_back(pos[view]-dpos[view]);
518  }
519 
520  // Loop over planes and high cells to make lines along top
521  // edge. Work from downstream edge toward upstream edge
522  std::map<int,int>::reverse_iterator ritr (cellhi.rbegin());
523  std::map<int,int>::reverse_iterator ritrEnd(cellhi.rend());
524  for (; ritr!=ritrEnd; ++ritr) {
525  const int plane = ritr->first;
526  const int cell = ritr->second;
527  double pos[3], dpos[3];
528  geo->CellInfo(plane, cell, &view, pos, dpos);
529 
530  vzZpts.push_back(pos[2]+dpos[2]);
531  vzVpts.push_back(pos[view]+dpos[view]);
532 
533  vzZpts.push_back(pos[2]-dpos[2]);
534  vzVpts.push_back(pos[view]+dpos[view]);
535  }
536 
537  // Add link to starting point to close the box
538  vzZpts.push_back(vzZpts[0]);
539  vzVpts.push_back(vzVpts[0]);
540  }
rb::Cluster::NCell
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
Definition: Cluster.cxx:134
geo::View_t
enum geo::_plane_proj View_t
Enumerate the possible plane projections.
rb::CellHit::Plane
unsigned short Plane() const
Definition: CellHit.h:39
geo::GeometryBase::CellInfo
void CellInfo(unsigned int ip, unsigned int ic, View_t *view=0, double *pos=0, double *dpos=0) const
Definition: GeometryBase.cxx:813
rb::CellHit::Cell
unsigned short Cell() const
Definition: CellHit.h:40
calib::j
const double j
Definition: BetheBloch.cxx:29
rb::Cluster::Cell
art::Ptr< rb::CellHit > Cell(geo::View_t view, unsigned int viewIdx) const
Get the ith cell from view view.
Definition: Cluster.cxx:145
geo
Helper for AttenCurve.
Definition: Path.h:10
void evd::RecoBaseDrawer::GetClusterOutlinesConvexHull ( const rb::Cluster *  c,
int  viewint,
std::vector< double > &  vzVpts,
std::vector< double > &  vzZpts 
) const
private

Definition at line 545 of file RecoBaseDrawer.cxx.

References ana::assert(), rb::CellHit::Cell(), rb::Cluster::Cell(), stan::math::dot(), DrawHoughResults2D(), dz, e, geom(), geo::CellGeo::GetCenter(), geo::CellGeo::HalfD(), geo::CellGeo::HalfW(), getGoodRuns4SAM::n, rb::Cluster::NCell(), rb::CellHit::Plane(), geo::GeometryBase::Plane(), util::pythag(), and POTSpillRate::view.

Referenced by DrawCluster2D(), and GetClusterOutlines().

549  {
550  // This is the "wrapping-paper" algorithm
551 
552  const geo::View_t view = geo::View_t(viewint);
553 
554  std::vector<double> vs;
555  std::vector<double> zs;
556 
557  art::ServiceHandle<geo::Geometry> geom;
558 
559  for(unsigned int cellIdx = 0; cellIdx < c->NCell(view); ++cellIdx){
560  const art::Ptr<rb::CellHit>& chit = c->Cell(view, cellIdx);
561  const geo::CellGeo* geocell = geom->Plane(chit->Plane())->Cell(chit->Cell());
562  double xyz[3];
563  geocell->GetCenter(xyz);
564  const double v0 = xyz[view];
565  const double z0 = xyz[2];
566  const double dv = geocell->HalfW();
567  const double dz = geocell->HalfD();
568  vs.push_back(v0-dv); zs.push_back(z0-dz);
569  vs.push_back(v0-dv); zs.push_back(z0+dz);
570  vs.push_back(v0+dv); zs.push_back(z0+dz);
571  vs.push_back(v0+dv); zs.push_back(z0-dz);
572  } // end for cellIdx
573 
574  const unsigned int N = vs.size();
575  if(N == 0) return;
576 
577  unsigned int minIdx = 0;
578 
579  double minz = +1e10;
580  for(unsigned int n = 0; n < N; ++n){
581  if(zs[n] < minz){
582  minz = zs[n];
583  minIdx = n;
584  }
585  } // end for n
586 
587  vzVpts.push_back(vs[minIdx]);
588  vzZpts.push_back(zs[minIdx]);
589 
590  // The way the "wrapping-paper" is pointing. We're working clockwise from the far left.
591  double dv = 1;
592  double dz = 0;
593  unsigned int curIdx = minIdx;
594 
595  while(true){
596  const double curV = vs[curIdx];
597  const double curZ = zs[curIdx];
598 
599  double bestDot = -2;
600  int bestIdx = -1;
601  for(unsigned int n = 0; n < N; ++n){
602  if(n == curIdx) continue;
603 
604  const double ourdV = vs[n]-curV;
605  const double ourdZ = zs[n]-curZ;
606  const double dot = (ourdV*dv+ourdZ*dz)/util::pythag(ourdV, ourdZ);;
607  // Need to give the starting point a little bit of an advantage,
608  // otherwise we can miss it out and end up in an infinite loop.
609  if(dot > bestDot+((bestIdx == int(minIdx)) ? 1e-6 : 0) ||
610  (dot > bestDot-1e-3 && n == minIdx)){
611  bestDot = dot;
612  bestIdx = n;
613  }
614  } // end for n
615 
616  assert(bestIdx >= 0);
617 
618  vzVpts.push_back(vs[bestIdx]);
619  vzZpts.push_back(zs[bestIdx]);
620 
621  dv = vs[bestIdx]-vs[curIdx];
622  dz = zs[bestIdx]-zs[curIdx];
623  const double mag = util::pythag(dv, dz);
624  dv /= mag;
625  dz /= mag;
626  curIdx = bestIdx;
627 
628  if(curIdx == minIdx) break;
629  }
630 
631  /*
632  double area = 0;
633  for(unsigned int n = 1; n < vzZpts.size()-2; ++n){
634  const double xprod = (vzZpts[n+1]-vzZpts[0])*(vzVpts[n]-vzVpts[0])-(vzVpts[n+1]-vzVpts[0])*(vzZpts[n]-vzZpts[0]);
635  assert(xprod > -1e-6);
636  area += xprod/2;
637  }
638  */
639  }
geo::CellGeo::HalfD
double HalfD() const
Definition: CellGeo.cxx:205
rb::Cluster::NCell
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
Definition: Cluster.cxx:134
geo::CellGeo::GetCenter
void GetCenter(double *xyz, double localz=0.0) const
Definition: CellGeo.cxx:159
geo::View_t
enum geo::_plane_proj View_t
Enumerate the possible plane projections.
rb::CellHit::Plane
unsigned short Plane() const
Definition: CellHit.h:39
geo::CellGeo::HalfW
double HalfW() const
Definition: CellGeo.cxx:191
geo::GeometryBase::Plane
const PlaneGeo * Plane(unsigned int i) const
Definition: GeometryBase.cxx:749
rb::CellHit::Cell
unsigned short Cell() const
Definition: CellHit.h:40
dz
double dz[NP][NC]
Definition: checkMeasured_Target.C:20
rb::Cluster::Cell
art::Ptr< rb::CellHit > Cell(geo::View_t view, unsigned int viewIdx) const
Get the ith cell from view view.
Definition: Cluster.cxx:145
stan::math::dot
double dot(const std::vector< double > &x, const std::vector< double > &y)
Definition: dot.hpp:10
util::pythag
double pythag(double x, double y)
2D Euclidean distance
Definition: MathUtil.h:29
geom
void geom(int which=0)
Definition: geom.C:163
ana::assert
assert(nhit_max >=nhit_nbins)
e
Float_t e
Definition: plot.C:35
geo::CellGeo
Encapsulate the cell geometry.
Definition: CellGeo.h:25
void evd::RecoBaseDrawer::HoughResult2D ( art::Event const &  evt,
evdb::View2D *  xzview,
evdb::View2D *  yzview 
)

Definition at line 984 of file RecoBaseDrawer.cxx.

References DrawHoughResults2D(), evd::RecoDrawingOptions::fHoughOpt, evd::RecoDrawingOptions::fHoughResultInstances, evd::RecoDrawingOptions::fHoughResultModules, evd::SliceNavigator::GetProducts(), evd::RecoDrawingOptions::kOFF, and string.

987  {
988  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
989  art::ServiceHandle<evd::SliceNavigator> nav;
990  unsigned int imod;
991 
992  if(drawopt->fHoughOpt == evd::RecoDrawingOptions::kOFF) return;
993 
994  for (imod=0; imod<drawopt->fHoughResultModules.size(); ++imod) {
995  std::string const& which = drawopt->fHoughResultModules[imod];
996  std::string const& instance = drawopt->fHoughResultInstances[imod];
997 
998  std::vector<const rb::HoughResult*> htr;
999  nav->GetProducts(evt, which, instance, htr);
1000  this->DrawHoughResults2D(htr, xzview, yzview);
1001  }
1002  }
evd::RecoBaseDrawer::DrawHoughResults2D
void DrawHoughResults2D(std::vector< const rb::HoughResult * > &htr, evdb::View2D *xzview, evdb::View2D *yzview)
Definition: RecoBaseDrawer.cxx:644
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
evd::RecoDrawingOptions::fHoughResultModules
std::vector< std::string > fHoughResultModules
Module labels for HoughResults.
Definition: RecoDrawingOptions.h:64
evt
int evt
Definition: FillTrainingTrees.h:22
evd::RecoDrawingOptions::fHoughResultInstances
std::vector< std::string > fHoughResultInstances
productInstance labels for HoughResults
Definition: RecoDrawingOptions.h:65
string
enum BeamMode string
void evd::RecoBaseDrawer::OfflineChans2D ( const art::Event &  evt,
evdb::View2D *  xzview,
evdb::View2D *  yzview 
)

Definition at line 427 of file RecoBaseDrawer.cxx.

References evdb::View2D::AddMarker(), plot_validation_datamc::c, geo::OfflineChan::Cell(), plot_validation_datamc_2018::color, evd::RecoDrawingOptions::fOfflineChanInstances, evd::RecoDrawingOptions::fOfflineChanModules, evd::RecoDrawingOptions::fOfflineChanOpt, geom(), GetClusterColor(), evd::SliceNavigator::GetProducts(), evd::RecoDrawingOptions::kOFF, geo::kX, geo::OfflineChan::Plane(), geo::GeometryBase::Plane(), string, and POTSpillRate::view.

430  {
431  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
432  if(drawopt->fOfflineChanOpt == evd::RecoDrawingOptions::kOFF) return;
433 
434  art::ServiceHandle<geo::Geometry> geom;
435  art::ServiceHandle<evd::SliceNavigator> nav;
436 
437  for(unsigned int imod = 0; imod < drawopt->fOfflineChanModules.size(); ++imod){
438  const char* which = drawopt->fOfflineChanModules[imod].c_str();
439  const std::string& instance = drawopt->fOfflineChanInstances[imod];
440 
441  std::vector<const geo::OfflineChan*> chans;
442  nav->GetProducts(evt, which, instance, chans);
443  for(unsigned int chanIdx = 0; chanIdx < chans.size(); ++chanIdx){
444  const geo::OfflineChan* c = chans[chanIdx];
445 
446  const geo::View_t view = geom->Plane(c->Plane())->View();
447 
448  evdb::View2D* viewobj = (view == geo::kX) ? xzview : yzview;
449  if(!viewobj) continue;
450 
451  double xyz[3];
452  geom->Plane(c->Plane())->Cell(c->Cell())->GetCenter(xyz);
453 
454  int color = GetClusterColor(imod);
455  // color = this->DimOutOfTime(c->MinTNS(), c->MaxTNS(), color);
456 
457  viewobj->AddMarker(xyz[2], xyz[view],
458  color,
459  kFullCircle, 0.8);
460  } // end for chanIdx
461  } // end for imod
462  }
geo::View_t
enum geo::_plane_proj View_t
Enumerate the possible plane projections.
plot_validation_datamc_2018.color
color
Definition: plot_validation_datamc_2018.py:516
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
geo::GeometryBase::Plane
const PlaneGeo * Plane(unsigned int i) const
Definition: GeometryBase.cxx:749
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
plot_validation_datamc.c
c
Definition: plot_validation_datamc.py:198
evd::RecoDrawingOptions::fOfflineChanInstances
std::vector< std::string > fOfflineChanInstances
productInstance labels for OfflineChans
Definition: RecoDrawingOptions.h:62
View
Definition: View.py:1
geo::OfflineChan::Plane
unsigned short Plane() const
Definition: OfflineChan.h:31
geo::OfflineChan::Cell
unsigned short Cell() const
Definition: OfflineChan.h:32
evd::RecoDrawingOptions::fOfflineChanModules
std::vector< std::string > fOfflineChanModules
Module labels for OfflineChans.
Definition: RecoDrawingOptions.h:61
geo::OfflineChan
A (plane, cell) pair.
Definition: OfflineChan.h:17
geom
void geom(int which=0)
Definition: geom.C:163
evd::RecoBaseDrawer::GetClusterColor
int GetClusterColor(unsigned int idx) const
Definition: RecoBaseDrawer.cxx:1152
evdb::View2D::AddMarker
TMarker & AddMarker(double x, double y, int c, int st, double sz)
Definition: View2D.cxx:124
string
enum BeamMode string
void evd::RecoBaseDrawer::Prong2D ( art::Event const &  evt,
evdb::View2D *  xzview,
evdb::View2D *  yzview 
)

Definition at line 1006 of file RecoBaseDrawer.cxx.

References DrawCluster2D(), DrawProng2D(), evd::RecoDrawingOptions::fProngIndex, evd::RecoDrawingOptions::fProngInstances, evd::RecoDrawingOptions::fProngModules, evd::RecoDrawingOptions::fProngOpt, evd::RecoDrawingOptions::fShowerOpt, evd::SliceNavigator::GetProducts(), evd::RecoDrawingOptions::kAS_CLUSTER, evd::RecoDrawingOptions::kAS_PRONG, evd::RecoDrawingOptions::kOFF, and string.

1009  {
1010  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
1011  art::ServiceHandle<evd::SliceNavigator> nav;
1012 
1013  bool drawProng =(drawopt->fProngOpt != 0);
1014  bool drawShowers=(drawopt->fShowerOpt != evd::RecoDrawingOptions::kOFF);
1015 
1016  if (drawProng || drawShowers) {
1017  for (unsigned int imod=0; imod<drawopt->fProngModules.size(); ++imod) {
1018  std::string const& which = drawopt->fProngModules[imod];
1019  std::string const& instance = drawopt->fProngInstances[imod];
1020 
1021  std::vector<const rb::Prong*> prongs;
1022  nav->GetProducts(evt, which, instance, prongs);
1023  std::vector<int> prongIndex;
1024  if( drawopt->fProngIndex.size() > imod )
1025  prongIndex = drawopt->fProngIndex[imod];
1026 
1027  if(drawopt->fProngOpt & evd::RecoDrawingOptions::kAS_CLUSTER){
1028  std::vector<const rb::Cluster*> clusts(prongs.begin(), prongs.end());
1029  DrawCluster2D(clusts, xzview, yzview, prongIndex);
1030  }
1031 
1032  if(drawopt->fProngOpt & evd::RecoDrawingOptions::kAS_PRONG){
1033  DrawProng2D(prongs, xzview, yzview, prongIndex);
1034  }
1035  }
1036  }
1037  }
evd::RecoBaseDrawer::DrawCluster2D
void DrawCluster2D(std::vector< const rb::Cluster * > &clusts, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &clustIndex)
Definition: RecoBaseDrawer.cxx:262
evd::RecoDrawingOptions::kAS_CLUSTER
static const int kAS_CLUSTER
Definition: RecoDrawingOptions.h:40
evd::RecoDrawingOptions::kAS_PRONG
static const int kAS_PRONG
Definition: RecoDrawingOptions.h:41
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
evd::RecoBaseDrawer::DrawProng2D
void DrawProng2D(std::vector< const rb::Prong * > &prong, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &prongIndex)
Definition: RecoBaseDrawer.cxx:746
evt
int evt
Definition: FillTrainingTrees.h:22
evd::RecoDrawingOptions::fProngInstances
std::vector< std::string > fProngInstances
productInstance labels for Prongs
Definition: RecoDrawingOptions.h:71
evd::RecoDrawingOptions::fProngModules
std::vector< std::string > fProngModules
Module labels for Prongs.
Definition: RecoDrawingOptions.h:70
evd::RecoDrawingOptions::fProngIndex
std::vector< std::vector< int > > fProngIndex
Definition: RecoDrawingOptions.h:82
string
enum BeamMode string
void evd::RecoBaseDrawer::Prong3D ( art::Event const &  evt,
evdb::View3D *  view 
)

Definition at line 1091 of file RecoBaseDrawer.cxx.

References DrawProng3D(), DrawTrack3D(), evd::RecoDrawingOptions::fClusterModules, evd::RecoDrawingOptions::fProngIndex, evd::RecoDrawingOptions::fProngInstances, evd::RecoDrawingOptions::fProngModules, evd::RecoDrawingOptions::fProngOpt, evd::RecoDrawingOptions::fShowerOpt, evd::RecoDrawingOptions::fTrackIndex, evd::RecoDrawingOptions::fTrackInstances, evd::RecoDrawingOptions::fTrackModules, evd::RecoDrawingOptions::fTrackOpt, evd::SliceNavigator::GetProducts(), makeTrainCVSamples::int, evd::RecoDrawingOptions::kOFF, string, and track.

1093  {
1094  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
1095  art::ServiceHandle<evd::SliceNavigator> nav;
1096 
1097  bool drawProng =(drawopt->fProngOpt != 0);
1098  bool drawShowers=(drawopt->fShowerOpt != evd::RecoDrawingOptions::kOFF);
1099  bool drawTracks =(drawopt->fTrackOpt != 0);
1100 
1101  if (drawProng || drawShowers) {
1102  for (unsigned int imod=0; imod<drawopt->fProngModules.size(); ++imod) {
1103  std::string const& which = drawopt->fProngModules[imod];
1104  std::string const& instance = drawopt->fProngInstances[imod];
1105 
1106  std::vector<const rb::Prong*> prong;
1107  nav->GetProducts(evt, which, instance, prong);
1108 
1109  std::vector<int> prongIndex;
1110  if( drawopt->fProngIndex.size() > imod )
1111  prongIndex = drawopt->fProngIndex[imod];
1112 
1113  DrawProng3D(prong, view, prongIndex);
1114  }
1115  }
1116 
1117  if (drawTracks) {
1118  for (unsigned int imod=0; imod<drawopt->fTrackModules.size(); ++imod) {
1119  std::string const& whichTrack = drawopt->fTrackModules[imod];
1120  std::string whichMichel = "";
1121  std::string const& instance = drawopt->fTrackInstances[imod];
1122 
1123  for (int jmod=0; jmod< (int)drawopt->fClusterModules.size(); ++jmod){
1124 
1125  std::string const& clustLabel = drawopt->fClusterModules[jmod];
1126  std::size_t foundMichel = clustLabel.find("michel");
1127  std::size_t foundTrack = clustLabel.find(whichTrack);
1128 
1129  if( foundMichel!=std::string::npos &&
1130  foundTrack!=std::string::npos)
1131  whichMichel = clustLabel;
1132  }
1133 
1134 
1135  std::vector<const rb::Track*> track;
1136  if(whichMichel=="")
1137  nav->GetProducts(evt, whichTrack, instance, track);
1138  else
1139  nav->GetProducts(evt, whichTrack, whichMichel, instance, track);
1140 
1141  std::vector<int> trackIndex;
1142  if( drawopt->fTrackIndex.size() > imod )
1143  trackIndex = drawopt->fTrackIndex[imod];
1144 
1145  DrawTrack3D(track, view, trackIndex);
1146  }
1147  }
1148  }
evd::RecoDrawingOptions::fClusterModules
std::vector< std::string > fClusterModules
Module labels for Clusters.
Definition: RecoDrawingOptions.h:58
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
evd::RecoDrawingOptions::fTrackModules
std::vector< std::string > fTrackModules
Module labels for Tracks.
Definition: RecoDrawingOptions.h:73
evt
int evt
Definition: FillTrainingTrees.h:22
evd::RecoDrawingOptions::fProngInstances
std::vector< std::string > fProngInstances
productInstance labels for Prongs
Definition: RecoDrawingOptions.h:71
evd::RecoBaseDrawer::DrawProng3D
void DrawProng3D(std::vector< const rb::Prong * > &prong, evdb::View3D *view, std::vector< int > &prongIndex)
Definition: RecoBaseDrawer.cxx:879
evd::RecoBaseDrawer::DrawTrack3D
void DrawTrack3D(std::vector< const rb::Track * > &track, evdb::View3D *view, std::vector< int > &trackIndex)
Definition: RecoBaseDrawer.cxx:923
evd::RecoDrawingOptions::fProngModules
std::vector< std::string > fProngModules
Module labels for Prongs.
Definition: RecoDrawingOptions.h:70
evd::RecoDrawingOptions::fTrackInstances
std::vector< std::string > fTrackInstances
productInstance labels for Tracks
Definition: RecoDrawingOptions.h:74
evd::RecoDrawingOptions::fTrackIndex
std::vector< std::vector< int > > fTrackIndex
Definition: RecoDrawingOptions.h:80
track
Float_t track
Definition: plot.C:35
evd::RecoDrawingOptions::fProngIndex
std::vector< std::vector< int > > fProngIndex
Definition: RecoDrawingOptions.h:82
string
enum BeamMode string
double evd::RecoBaseDrawer::ScaleFactor ( double  ,
double  q 
)
private

Definition at line 248 of file RecoBaseDrawer.cxx.

References evd::RawDrawingOptions::fRawDrawingOpt, evd::RawDrawingOptions::fScaleFactor, evd::RawDrawingOptions::kSCALE_HITS_BY_CHARGE, and std::sqrt().

Referenced by CellHit2D().

249 {
250  art::ServiceHandle<evd::RawDrawingOptions> drawopt;
251 
252  double sf = drawopt->fScaleFactor;
253  double q0 = 3.0;
254  if (drawopt->fRawDrawingOpt & evd::RawDrawingOptions::kSCALE_HITS_BY_CHARGE)
255  {
256  sf *= 0.3*sqrt(q/q0);
257  }
258  return sf;
259 }
std::sqrt
T sqrt(T number)
Definition: d0nt_math.hpp:156
evd::RawDrawingOptions::kSCALE_HITS_BY_CHARGE
static const int kSCALE_HITS_BY_CHARGE
Definition: RawDrawingOptions.h:47
void evd::RecoBaseDrawer::Track2D ( const art::Event &  evt,
evdb::View2D *  xzview,
evdb::View2D *  yzview 
)

Definition at line 1041 of file RecoBaseDrawer.cxx.

References DrawCluster2D(), DrawProng2D(), DrawTrack2D(), evd::RecoDrawingOptions::fClusterModules, evd::RecoDrawingOptions::fTrackIndex, evd::RecoDrawingOptions::fTrackInstances, evd::RecoDrawingOptions::fTrackModules, evd::RecoDrawingOptions::fTrackOpt, evd::SliceNavigator::GetProducts(), evd::RecoDrawingOptions::kAS_CLUSTER, evd::RecoDrawingOptions::kAS_PRONG, evd::RecoDrawingOptions::kAS_TRACK, and string.

1044  {
1045  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
1046  art::ServiceHandle<evd::SliceNavigator> nav;
1047 
1048  if(drawopt->fTrackOpt == 0) return;
1049 
1050  for(unsigned int imod = 0; imod < drawopt->fTrackModules.size(); ++imod){
1051  const std::string& which = drawopt->fTrackModules[imod];
1052  const std::string& instance = drawopt->fTrackInstances[imod];
1053 
1054  std::string whichMichel;
1055  for(unsigned int jmod = 0; jmod < drawopt->fClusterModules.size(); ++jmod){
1056  const std::string& clustLabel = drawopt->fClusterModules[jmod];
1057  bool foundMichel = clustLabel.find("michel") != std::string::npos;
1058  bool foundTrack = clustLabel.find(which) != std::string::npos;
1059 
1060  if(foundMichel && foundTrack) whichMichel = clustLabel;
1061  }
1062 
1063  std::vector<const rb::Track*> tracks;
1064  if(whichMichel.empty())
1065  nav->GetProducts(evt, which, instance, tracks);
1066  else
1067  nav->GetProducts(evt, which, whichMichel, instance, tracks);
1068 
1069  std::vector<int> trackIndex;
1070  if(drawopt->fTrackIndex.size() > imod)
1071  trackIndex = drawopt->fTrackIndex[imod];
1072 
1073  if(drawopt->fTrackOpt & evd::RecoDrawingOptions::kAS_CLUSTER){
1074  std::vector<const rb::Cluster*> clusts(tracks.begin(), tracks.end());
1075  DrawCluster2D(clusts, xzview, yzview, trackIndex);
1076  }
1077 
1078  if(drawopt->fTrackOpt & evd::RecoDrawingOptions::kAS_PRONG){
1079  std::vector<const rb::Prong*> prongs(tracks.begin(), tracks.end());
1080  DrawProng2D(prongs, xzview, yzview, trackIndex);
1081  }
1082 
1083  if(drawopt->fTrackOpt & evd::RecoDrawingOptions::kAS_TRACK){
1084  DrawTrack2D(tracks, xzview, yzview, trackIndex);
1085  }
1086  }
1087  }
evd::RecoDrawingOptions::fClusterModules
std::vector< std::string > fClusterModules
Module labels for Clusters.
Definition: RecoDrawingOptions.h:58
evd::RecoBaseDrawer::DrawCluster2D
void DrawCluster2D(std::vector< const rb::Cluster * > &clusts, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &clustIndex)
Definition: RecoBaseDrawer.cxx:262
evd::RecoDrawingOptions::kAS_CLUSTER
static const int kAS_CLUSTER
Definition: RecoDrawingOptions.h:40
evd::RecoDrawingOptions::kAS_PRONG
static const int kAS_PRONG
Definition: RecoDrawingOptions.h:41
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
evd::RecoBaseDrawer::DrawProng2D
void DrawProng2D(std::vector< const rb::Prong * > &prong, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &prongIndex)
Definition: RecoBaseDrawer.cxx:746
evd::RecoBaseDrawer::DrawTrack2D
void DrawTrack2D(std::vector< const rb::Track * > &track, evdb::View2D *xzview, evdb::View2D *yzview, std::vector< int > &trackIndex)
Definition: RecoBaseDrawer.cxx:810
evd::RecoDrawingOptions::fTrackModules
std::vector< std::string > fTrackModules
Module labels for Tracks.
Definition: RecoDrawingOptions.h:73
evd::RecoDrawingOptions::fTrackInstances
std::vector< std::string > fTrackInstances
productInstance labels for Tracks
Definition: RecoDrawingOptions.h:74
evd::RecoDrawingOptions::kAS_TRACK
static const int kAS_TRACK
Definition: RecoDrawingOptions.h:42
evd::RecoDrawingOptions::fTrackIndex
std::vector< std::vector< int > > fTrackIndex
Definition: RecoDrawingOptions.h:80
string
enum BeamMode string
void evd::RecoBaseDrawer::Vertex2D ( art::Event const &  evt,
evdb::View2D *  xzview,
evdb::View2D *  yzview 
)

Definition at line 958 of file RecoBaseDrawer.cxx.

References evd::RecoDrawingOptions::fVertexIndex, evd::RecoDrawingOptions::fVertexInstances, evd::RecoDrawingOptions::fVertexModules, evd::RecoDrawingOptions::fVertexOpt, evd::SliceNavigator::GetProducts(), evd::RecoDrawingOptions::kOFF, and string.

961  {
962  art::ServiceHandle<evd::RecoDrawingOptions> drawopt;
963 
964  if(drawopt->fVertexOpt == evd::RecoDrawingOptions::kOFF) return;
965 
966  art::ServiceHandle<evd::SliceNavigator> nav;
967 
968  unsigned int imod;
969  for (imod=0; imod<drawopt->fVertexModules.size(); ++imod) {
970  std::string const& which = drawopt->fVertexModules[imod];
971  std::string const& instance = drawopt->fVertexInstances[imod];
972 
973  std::vector<const rb::Vertex*> vert;
974  nav->GetProducts(evt, which, instance, vert);
975  std::vector<int> vertexIndex;
976  if( drawopt->fVertexIndex.size() > imod )
977  vertexIndex = drawopt->fVertexIndex[imod];
978 
979  this->DrawVertices2D(vert, xzview, yzview, vertexIndex);
980  }
981  }
evd::RecoDrawingOptions::fVertexInstances
std::vector< std::string > fVertexInstances
productInstance labels for Vertexes
Definition: RecoDrawingOptions.h:68
evd::SliceNavigator::GetProducts
size_t GetProducts(const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
Definition: SliceNavigator_service.cc:206
evd::RecoDrawingOptions::fVertexModules
std::vector< std::string > fVertexModules
Module labels for Vertices.
Definition: RecoDrawingOptions.h:67
evt
int evt
Definition: FillTrainingTrees.h:22
evd::RecoDrawingOptions::fVertexIndex
std::vector< std::vector< int > > fVertexIndex
Definition: RecoDrawingOptions.h:81
string
enum BeamMode string
The documentation for this class was generated from the following files:
Generated on Fri Jan 15 2021 14:17:48 for NOvA by   doxygen 1.8.11