66 typedef std::map< unsigned int, FebPlotGroup >
FebPlotSet;
76 unsigned int&,
unsigned int&,
unsigned int&,
unsigned int&);
127 group,diblock,ldcm,lfeb);
150 tfs->
make<TH2F>(Form(
"HBookedChanCP_%s_%d",viewName.c_str(),
group),
151 Form(
"Booked Channel, %s, group %d",viewName.c_str(),
group),
160 for (
unsigned int ipix=0;
162 const unsigned int chan = cmap->
Map()->
encodeDChan(det,diblock,ldcm,lfeb,ipix);
164 const unsigned int plannum = cmap->
Map()->
getPlane(logchan);
166 const unsigned int view=(
unsigned int)plane->
View();
176 std::cout<<
"$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ geo::kXorY = " <<
"Here" <<
std::endl;
191 tfs->
make<TH2F>(Form(
"HRawDigits_%s_%d",viewName.c_str(),
cut),
192 Form(
"Raw Digits, %s, %s",viewName.c_str(),cutname.c_str()),
204 tfs->
make<TH1F>(Form(
"HExtentPlaneTrk_%d",
cut),
205 Form(
"Track Extent, Planes, %s",cutname.c_str()),
208 tfs->
make<TH1F>(Form(
"HNCellTrk_%d",
cut),
209 Form(
"Track Number of Cells, %s",cutname.c_str()),
212 tfs->
make<TH1F>(Form(
"HNCellViewTrk_X_%d",
cut),
213 Form(
"Track Number of Cells, X, %s",cutname.c_str()),
216 tfs->
make<TH1F>(Form(
"HNCellViewTrk_Y_%d",
cut),
217 Form(
"Track Number of Cells, Y, %s",cutname.c_str()),
228 unsigned int hrfeb,
int det,
229 unsigned int &
group,
unsigned int &
diblock,
unsigned int &ldcm,
234 group= (hrfeb/1000000)%100;
235 diblock=(hrfeb/ 10000)%100;
236 ldcm= (hrfeb/ 100)%100;
237 lfeb= (hrfeb/ 1)%100;
239 std::string label(Form(
"g%d_%d_%2.2d_%2.2d_%2.2d",group,det,diblock,ldcm,lfeb));
263 for (
unsigned int id=0;
id<digs->size();
id++)
272 const unsigned int plannum = cmap->
Map()->
getPlane(logchan);
274 const unsigned int view = (
unsigned int)plane->
View();
295 if (!tracks->size()) {
306 mf::LogError(
"CPLTWarn") <<
"\nFound Tracks but no CellHits!\n";
310 for(
unsigned int tr = 0;
tr < tracks->size();
tr++)
318 for (
unsigned int ic=0; ic<track->
NCell(); ic++)
323 const unsigned int plannum = cmap->
Map()->
getPlane(logchan);
351 unsigned int lookupchan=chan &
355 for (FebPlotSet::iterator iter=plots.begin(); iter!=plots.end(); iter++)
359 if (iter->second.find(lookupchan) != iter->second.end()) {
360 iter->second[lookupchan]->AddHit(dig->
ADC(),
393 FebPlotSet::iterator siter=
fFebRawAll.begin();
395 FebPlotGroup::iterator giter=(siter->second).begin();
396 for (;giter!=(siter->second).
end();giter++)
397 (giter->second)->FinishEvent();
402 FebPlotGroup::iterator giter=(siter->second).begin();
403 for (;giter!=(siter->second).
end();giter++)
404 (giter->second)->FinishEvent();
410 FebPlotSet::iterator siter=
fFebRawAll.begin();
412 FebPlotGroup::iterator giter=(siter->second).begin();
413 for (;giter!=(siter->second).
end();giter++)
414 (giter->second)->Reset();
419 FebPlotGroup::iterator giter=(siter->second).begin();
420 for (;giter!=(siter->second).
end();giter++)
421 (giter->second)->Reset();
void beginRun(const art::Run &run)
const std::string fRawDigitLabel
diblock
print "ROW IS " print row
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
std::map< unsigned int, FebPlotGroup > FebPlotSet
std::map< unsigned int, RawFebPlots * > FebPlotGroup
APD Pixel Number (8bit), valid range 0-31.
enum geo::_plane_proj View_t
Enumerate the possible plane projections.
FebPlotSet fFebRawOnTrack
Vertical planes which measure X.
unsigned int Ncells() const
Number of cells in this plane.
TH1F * fHNCellTrk[kNumLevel]
const daqchannelmap::DAQChannelMap * Map() const
const unsigned int fMinCellTrack
const PlaneGeo * Plane(unsigned int i) const
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
DEFINE_ART_MODULE(TestTMapFile)
cell_t getCell(lchan logicalchan) const
Decode the cell number from an lchan.
lchan encodeLChan(int detId, plane_t plane, cell_t cell) const
Horizontal planes which measure Y.
const std::vector< unsigned int > fHumanReadableFEBs
std::string HRFebToStringLabel(unsigned int, int, unsigned int &, unsigned int &, unsigned int &, unsigned int &)
View_t View() const
Which coordinate does this plane measure.
uint32_t DaqChannel() const
const std::string fTrackLabel
correl_yv Fill(-(dy[iP-1][iC-1]), hyv->GetBinContent(iP, iC))
const std::vector< Plot > plots
novadaq::cnv::DetId DetId() const
Prefer ds::DetectorService::DetId() instead.
Geometry information for a single readout plane.
EDAnalyzer(Table< Config > const &config)
enum dprf::ChannelPlots::_cutlevel Cuts_t
unsigned int NYCell() const
Number of cells in the y-view.
virtual void analyze(art::Event const &)
art::Ptr< rb::CellHit > Cell(geo::View_t view, unsigned int viewIdx) const
Get the ith cell from view view.
pixel_t getPixel(dchan daqchan) const
Decode the pixel id from a dchan.
TH2F * fHRawDigits[geo::kXorY][kNumLevel]
std::map< unsigned int, ViewH > fHBookedChanCP
T * make(ARGS...args) const
const unsigned int fMinCellTrackView
int16_t ADC(uint32_t i) const
void FillFebRawPlots(FebPlotSet &, art::Ptr< rawdata::RawDigit >)
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
unsigned int NXCell() const
Number of cells in the x-view.
const unsigned int fMinPlaneTrack
const std::set< unsigned int > & GetPlanesByView(View_t v=kXorY) const
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
unsigned int ExtentPlane(geo::View_t view=geo::kXorY) const
const fhicl::ParameterSet & fPSet
TH1F * fHNCellViewTrk[geo::kXorY][kNumLevel]
plane_t getPlane(lchan logicalchan) const
Decode the plane number from an lchan.
assert(nhit_max >=nhit_nbins)
dchan encodeDChan(int detID, diblock_t diblock, dcm_id_t dcm, feb_t feb, pixel_t pixel) const
ChannelPlots(fhicl::ParameterSet const &p)
unsigned int NPlanes() const
const std::string fCellHitLabel
uint32_t dchan
< DAQ Channel Map Package
bool TrackOK(art::Ptr< rb::Track >)
Encapsulate the geometry of one entire detector (near, far, ndos)
TH1F * fHExtentPlaneTrk[kNumLevel]