Public Member Functions | Private Attributes | List of all members
AnalysisExample::DbigEVD4apaFD Class Reference
Inheritance diagram for AnalysisExample::DbigEVD4apaFD:
art::EDAnalyzer art::detail::Analyzer art::detail::LegacyModule art::Observer art::ModuleBase

Public Member Functions

 DbigEVD4apaFD (fhicl::ParameterSet const &pset)
 
virtual ~DbigEVD4apaFD ()
 
void beginJob ()
 
void beginRun (const art::Run &run)
 
void reconfigure (fhicl::ParameterSet const &pset)
 
void analyze (const art::Event &evt)
 
- Public Member Functions inherited from art::EDAnalyzer
 EDAnalyzer (fhicl::ParameterSet const &pset)
 
template<typename Config >
 EDAnalyzer (Table< Config > const &config)
 
std::string workerType () const
 
- Public Member Functions inherited from art::detail::Analyzer
virtual ~Analyzer () noexcept
 
 Analyzer (fhicl::ParameterSet const &pset)
 
template<typename Config >
 Analyzer (Table< Config > const &config)
 
void doBeginJob ()
 
void doEndJob ()
 
void doRespondToOpenInputFile (FileBlock const &fb)
 
void doRespondToCloseInputFile (FileBlock const &fb)
 
void doRespondToOpenOutputFiles (FileBlock const &fb)
 
void doRespondToCloseOutputFiles (FileBlock const &fb)
 
bool doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
 
bool doEndRun (RunPrincipal &rp, ModuleContext const &mc)
 
bool doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
 
bool doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
 
bool doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
 
- Public Member Functions inherited from art::Observer
 ~Observer () noexcept
 
 Observer (Observer const &)=delete
 
 Observer (Observer &&)=delete
 
Observer & operator= (Observer const &)=delete
 
Observer & operator= (Observer &&)=delete
 
void registerProducts (ProductDescriptions &, ModuleDescription const &)
 
void fillDescriptions (ModuleDescription const &)
 
std::string const & processName () const
 
bool wantAllEvents () const
 
bool wantEvent (Event const &e)
 
fhicl::ParameterSetID selectorConfig () const
 
Handle< TriggerResults > getTriggerResults (Event const &e) const
 
- Public Member Functions inherited from art::ModuleBase
virtual ~ModuleBase () noexcept
 
 ModuleBase ()
 
ModuleDescription const & moduleDescription () const
 
void setModuleDescription (ModuleDescription const &)
 
std::array< std::vector< ProductInfo >, NumBranchTypes > const & getConsumables () const
 
void sortConsumables (std::string const &current_process_name)
 
template<typename T , BranchType BT>
ViewToken< T > consumesView (InputTag const &tag)
 
template<typename T , BranchType BT>
ViewToken< T > mayConsumeView (InputTag const &tag)
 

Private Attributes

std::string fChanHitLabel
 
std::string fWidHitLabel
 
std::string fHitCheatLabel
 
std::string fClusLabel
 
unsigned int fAPA
 
unsigned int fEvent
 
std::map< unsigned int, unsigned int > fAPAToNumHits
 
std::vector< std::vector< TH2I * > > fAmbigWireU0
 
std::vector< std::vector< TH2I * > > fAmbigWireV0
 
std::vector< std::vector< TH2I * > > fAmbigWireU1
 
std::vector< std::vector< TH2I * > > fAmbigWireV1
 
std::vector< std::vector< TH2I * > > fDisambigWireU0
 
std::vector< std::vector< TH2I * > > fDisambigWireV0
 
std::vector< std::vector< TH2I * > > fDisambigWireU1
 
std::vector< std::vector< TH2I * > > fDisambigWireV1
 
std::vector< std::vector< TH2I * > > fCheatWireU0
 
std::vector< std::vector< TH2I * > > fCheatWireV0
 
std::vector< std::vector< TH2I * > > fCheatWireU1
 
std::vector< std::vector< TH2I * > > fCheatWireV1
 
std::vector< std::vector< TH2I * > > fTimeChanZ0
 
std::vector< std::vector< TH2I * > > fTimeChanZ1
 
unsigned int fNticks
 
unsigned int fUChMin
 
unsigned int fUChMax
 
unsigned int fVChMin
 
unsigned int fVChMax
 
unsigned int fZ0ChMin
 
unsigned int fZ0ChMax
 
unsigned int fZ1ChMin
 
unsigned int fZ1ChMax
 
unsigned int fnAPAs
 
unsigned int fnEvts
 
unsigned int fChansPerAPA
 
unsigned int fUWireMax
 
unsigned int fVWireMax
 
unsigned int fZ0WireMax
 
unsigned int fZ1WireMax
 
unsigned int fMinT
 
unsigned int fMaxT
 
unsigned int fMaxTimeRange
 
bool fTimeReverseZ0 = true
 
bool fGotFirstEvtNum = false
 
int fFirstEvt
 
art::ServiceHandle< geo::Geometry > fGeom
 

Additional Inherited Members

- Public Types inherited from art::EDAnalyzer
using WorkerType = WorkerT< EDAnalyzer >
 
using ModuleType = EDAnalyzer
 
- Protected Member Functions inherited from art::Observer
 Observer (fhicl::ParameterSet const &config)
 
 Observer (std::vector< std::string > const &paths, fhicl::ParameterSet const &config)
 
detail::ProcessAndEventSelectors & processAndEventSelectors ()
 
- Protected Member Functions inherited from art::ModuleBase
ConsumesCollector & consumesCollector ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > consumes (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > mayConsume (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 

Detailed Description

Definition at line 45 of file DbigEVD4apaFD_module.cc.

Constructor & Destructor Documentation

AnalysisExample::DbigEVD4apaFD::DbigEVD4apaFD ( fhicl::ParameterSet const &  pset)
explicit

Definition at line 134 of file DbigEVD4apaFD_module.cc.

135  : EDAnalyzer( parameterSet )
136  , fEvent(-1)
137  {
138  this->reconfigure(parameterSet);
139  }
art::EDAnalyzer::EDAnalyzer
EDAnalyzer(fhicl::ParameterSet const &pset)
Definition: EDAnalyzer.h:27
AnalysisExample::DbigEVD4apaFD::reconfigure
void reconfigure(fhicl::ParameterSet const &pset)
Definition: DbigEVD4apaFD_module.cc:144
AnalysisExample::DbigEVD4apaFD::~DbigEVD4apaFD ( )
virtual

Definition at line 162 of file DbigEVD4apaFD_module.cc.

162  {
163  }

Member Function Documentation

void AnalysisExample::DbigEVD4apaFD::analyze ( const art::Event &  evt)

Definition at line 313 of file DbigEVD4apaFD_module.cc.

313  {
314 
315  // Get the index for Hist vectors - not 1 if reading a file that
316  // started at a custom firstEvent in the fcl file
317  if(!fGotFirstEvtNum){ fFirstEvt = event.id().event(); fGotFirstEvtNum = true; }
318  // This works if all evts are contiguously numbered, which if produced by a filter
319  // may not be the case.
320  // unsigned int evtIndex = event.id().event() - fFirstEvt;
321  unsigned int evtIndex = ++fEvent;
322 
323  art::Handle< std::vector<recob::Hit> > ChHits;
324  event.getByLabel(fChanHitLabel, ChHits);
325  art::Handle< std::vector<recob::Hit> > WidHits;
326  event.getByLabel(fWidHitLabel, WidHits);
327  art::Handle< std::vector<recob::Hit> > CheatedHits;
328  event.getByLabel(fHitCheatLabel, CheatedHits);
329 
330 
331  // plot hits on a given side of an APA vs all
332  // wireIDs corresponding to the it channel
333  for( auto const& hit : (*ChHits) ){
334  uint32_t chan = hit.Channel();
335  unsigned int apa = std::floor( chan/fChansPerAPA );
336  double startT = hit.PeakTimeMinusRMS();
337  double endT = hit.PeakTimePlusRMS();
338  double charge = hit.PeakAmplitude();
339  std::vector< geo::WireID > wids = fGeom->ChannelToWire(chan);
340 
341  //if( fAPAToNumHits.count(apa) == 0 ) fAPAToNumHits[apa] = 0;
342  //fAPAToNumHits[apa]++;
343 
344  for(size_t i=0; i<wids.size(); i++){
345 
346  unsigned int wire = wids[i].Wire;
347 
348  if( wids[i].TPC % 2 == 0 ){
349 
350  if( fGeom->View(chan) == geo::kU ){
351  for( double t = startT; t <= endT; t++ ){
352  fAmbigWireU0[evtIndex][apa]->Fill(wire, t, charge);
353  }
354  }
355  else if( fGeom->View(chan) == geo::kV ){
356  for( double t = startT; t <= endT; t++ ){
357  fAmbigWireV0[evtIndex][apa]->Fill(wire, t, charge);
358  }
359  }
360  else if ( fGeom->View(chan) == geo::kZ ){
361  for( double t = startT; t <= endT; t++ ){
362  double Z0fillT = t;
363  if(fTimeReverseZ0) Z0fillT = -t; // time reverse for geometrically logical display
364  fTimeChanZ0[evtIndex][apa]->Fill(chan, Z0fillT, charge);
365  }
366  }
367 
368  } else if( wids[i].TPC % 2 == 1 ){
369 
370  if( fGeom->View(chan) == geo::kU ){
371  for( double t = startT; t <= endT; t++ ){
372  fAmbigWireU1[evtIndex][apa]->Fill(wire, t, charge);
373  }
374  }
375  else if( fGeom->View(chan) == geo::kV ){
376  for( double t = startT; t <= endT; t++ ){
377  fAmbigWireV1[evtIndex][apa]->Fill(wire, t, charge);
378  }
379  }
380  else if ( hit.View() == geo::kZ ){
381  for( double t = startT; t <= endT; t++ ){
382  fTimeChanZ1[evtIndex][apa]->Fill(chan, t, charge);
383  }
384  }
385 
386  }
387 
388  } // end wid loop
389  } // end channel hit per wireID loop
390 
391 
392 
393  // for now, let user know where other activity is from here
394 
395  //std::map<unsigned int, unsigned int>::iterator apa_it;
396  //for( apa_it=fAPAToNumHits.begin(); apa_it!=fAPAToNumHits.end(); apa_it++ )
397  // mf::LogVerbatim("DbigEVD4apaFD") << "APA " << apa_it->first
398  // << " has " << apa_it->second << " hits.";
399 
400 
401 
402  // plot hits on a given side of an APA vs. disambiguated wireID
403  for( auto const& hit : (*WidHits) ){
404  uint32_t chan = hit.Channel();
405  unsigned int apa = std::floor( chan/fChansPerAPA );
406  unsigned int wire = hit.WireID().Wire;
407  double startT = hit.PeakTimeMinusRMS();
408  double endT = hit.PeakTimePlusRMS();
409  double charge = hit.PeakAmplitude();
410 
411  if(hit.WireID().TPC % 2 == 0){
412 
413  if( fGeom->View(chan) == geo::kU ){
414  for( double t = startT; t <= endT; t++ ){
415  fDisambigWireU0[evtIndex][apa]->Fill(wire, t, charge);
416  }
417  }
418  if( fGeom->View(chan) == geo::kV ){
419  for( double t = startT; t <= endT; t++ ){
420  fDisambigWireV0[evtIndex][apa]->Fill(wire, t, charge);
421  }
422  }
423 
424  } else if(hit.WireID().TPC % 2 == 1){
425 
426  if( fGeom->View(chan) == geo::kU ){
427  for( double t = startT; t <= endT; t++ ){
428  fDisambigWireU1[evtIndex][apa]->Fill(wire, t, charge);
429  }
430  }
431  if( fGeom->View(chan) == geo::kV ){
432  for( double t = startT; t <= endT; t++ ){
433  fDisambigWireV1[evtIndex][apa]->Fill(wire, t, charge);
434  }
435  }
436 
437  }
438 
439  } // end wireID hit per wireID loop
440 
441 
442  // plot cheated hits on a given side of an APA vs. cheated wireID
443  for( auto const& hit : (*CheatedHits) ){
444  uint32_t chan = hit.Channel();
445  unsigned int apa = std::floor( chan/fChansPerAPA );
446  unsigned int wire = hit.WireID().Wire;
447  double startT = hit.PeakTimeMinusRMS();
448  double endT = hit.PeakTimePlusRMS();
449  double charge = hit.PeakAmplitude();
450 
451  if(hit.WireID().TPC % 2 == 0){
452 
453  if( fGeom->View(chan) == geo::kU ){
454  for( double t = startT; t <= endT; t++ ){
455  fCheatWireU0[evtIndex][apa]->Fill(wire, t, charge);
456  }
457  }
458  if( fGeom->View(chan) == geo::kV ){
459  for( double t = startT; t <= endT; t++ ){
460  fCheatWireV0[evtIndex][apa]->Fill(wire, t, charge);
461  }
462  }
463 
464  } else if(hit.WireID().TPC % 2 == 1){
465 
466  if( fGeom->View(chan) == geo::kU ){
467  for( double t = startT; t <= endT; t++ ){
468  fCheatWireU1[evtIndex][apa]->Fill(wire, t, charge);
469  }
470  }
471  if( fGeom->View(chan) == geo::kV ){
472  for( double t = startT; t <= endT; t++ ){
473  fCheatWireV1[evtIndex][apa]->Fill(wire, t, charge);
474  }
475  }
476 
477  }
478 
479  } // end cheated hit loop
480 
481 
482  return;
483  }
AnalysisExample::DbigEVD4apaFD::fAmbigWireU0
std::vector< std::vector< TH2I * > > fAmbigWireU0
Definition: DbigEVD4apaFD_module.cc:74
generate_CCQE_events.t
t
Definition: generate_CCQE_events.py:68
geo::kV
Planes which measure V.
Definition: geo_types.h:126
AnalysisExample::DbigEVD4apaFD::fDisambigWireU1
std::vector< std::vector< TH2I * > > fDisambigWireU1
Definition: DbigEVD4apaFD_module.cc:80
geo::GeometryCore::ChannelToWire
std::vector< geo::WireID > ChannelToWire(raw::ChannelID_t const channel) const
Returns a list of wires connected to the specified TPC channel.
Definition: GeometryCore.cxx:1121
geo::kZ
Planes which measure Z direction.
Definition: geo_types.h:128
AnalysisExample::DbigEVD4apaFD::fCheatWireV1
std::vector< std::vector< TH2I * > > fCheatWireV1
Definition: DbigEVD4apaFD_module.cc:85
AnalysisExample::DbigEVD4apaFD::fAmbigWireU1
std::vector< std::vector< TH2I * > > fAmbigWireU1
Definition: DbigEVD4apaFD_module.cc:76
geo::kU
Planes which measure U.
Definition: geo_types.h:125
AnalysisExample::DbigEVD4apaFD::fTimeChanZ0
std::vector< std::vector< TH2I * > > fTimeChanZ0
Definition: DbigEVD4apaFD_module.cc:86
AnalysisExample::DbigEVD4apaFD::fDisambigWireU0
std::vector< std::vector< TH2I * > > fDisambigWireU0
Definition: DbigEVD4apaFD_module.cc:78
AnalysisExample::DbigEVD4apaFD::fCheatWireV0
std::vector< std::vector< TH2I * > > fCheatWireV0
Definition: DbigEVD4apaFD_module.cc:83
uint32_t
unsigned int uint32_t
Definition: stdint.h:126
AnalysisExample::DbigEVD4apaFD::fAmbigWireV0
std::vector< std::vector< TH2I * > > fAmbigWireV0
Definition: DbigEVD4apaFD_module.cc:75
AnalysisExample::DbigEVD4apaFD::fGeom
art::ServiceHandle< geo::Geometry > fGeom
Definition: DbigEVD4apaFD_module.cc:127
geo::GeometryCore::View
View_t View(geo::PlaneID const &pid) const
Returns the view (wire orientation) on the channels of specified TPC plane.
Definition: GeometryCore.cxx:603
AnalysisExample::DbigEVD4apaFD::fCheatWireU1
std::vector< std::vector< TH2I * > > fCheatWireU1
Definition: DbigEVD4apaFD_module.cc:84
hit
Detector simulation of raw signals on wires.
Definition: RobustMCAna_module.cc:49
AnalysisExample::DbigEVD4apaFD::fDisambigWireV1
std::vector< std::vector< TH2I * > > fDisambigWireV1
Definition: DbigEVD4apaFD_module.cc:81
AnalysisExample::DbigEVD4apaFD::fAmbigWireV1
std::vector< std::vector< TH2I * > > fAmbigWireV1
Definition: DbigEVD4apaFD_module.cc:77
AnalysisExample::DbigEVD4apaFD::fDisambigWireV0
std::vector< std::vector< TH2I * > > fDisambigWireV0
Definition: DbigEVD4apaFD_module.cc:79
AnalysisExample::DbigEVD4apaFD::fTimeChanZ1
std::vector< std::vector< TH2I * > > fTimeChanZ1
Definition: DbigEVD4apaFD_module.cc:87
AnalysisExample::DbigEVD4apaFD::fCheatWireU0
std::vector< std::vector< TH2I * > > fCheatWireU0
Definition: DbigEVD4apaFD_module.cc:82
void AnalysisExample::DbigEVD4apaFD::beginJob ( )
virtual

Reimplemented from art::EDAnalyzer.

Definition at line 167 of file DbigEVD4apaFD_module.cc.

167  {
168 
169  // loop through channels in the first APA to find the channel boundaries for each view
170  // will adjust for desired APA after
171  fUChMin = 0;
172  fChansPerAPA = ( fGeom->Nchannels() )/ (fGeom->NTPC()/2 );
173  fZ1ChMax = ( fGeom->Nchannels() )/ (fGeom->NTPC()/2 ) - 1;
174  for ( unsigned int c = fUChMin + 1; c < fZ1ChMax; c++ ){
175  if ( fGeom->View(c) == geo::kV && fGeom->View(c-1) == geo::kU ){
176  fVChMin = c;
177  fUChMax = c - 1;
178  }
179  if ( fGeom->View(c) == geo::kZ && fGeom->View(c-1) == geo::kV ){
180  fZ0ChMin = c;
181  fVChMax = c-1;
182  }
183  if ( fGeom->View(c) == geo::kZ && fGeom->ChannelToWire(c)[0].TPC == fGeom->ChannelToWire(c-1)[0].TPC + 1 ){
184  fZ1ChMin = c;
185  fZ0ChMax = c-1;
186  }
187  }
188 
189  //fUWireMax = fGeom->Nwires(0, std::floor( fAPA/2 ), 0);
190  //fVWireMax = fGeom->Nwires(1, std::floor( fAPA/2 ), 0);
191  //fZ0WireMax = fGeom->Nwires(2, std::floor( fAPA/2 ), 0);
192  fUWireMax = fGeom->Nwires(0, 0, 0);
193  fVWireMax = fGeom->Nwires(1, 0, 0);
194  fZ0WireMax = fGeom->Nwires(2, 0, 0);
195  fZ1WireMax = fZ0WireMax;
196 
197 
198  unsigned int minT = 0;
199  unsigned int maxT = fNticks;
200 
201  float ThumbFactor = 0.85;
202  unsigned int binT = (maxT-minT)*ThumbFactor;
203 
204  // Access ART's TFileService, which will handle creating and writing
205  // histograms and n-tuples for us.
206  art::ServiceHandle<art::TFileService> tfs;
207 
208  int minZ0T = minT;
209  int maxZ0T = maxT;
210  // makes for more logical display in the root display script
211  if(fTimeReverseZ0){ minZ0T = -fNticks;
212  maxZ0T = 0; }
213 
214 
215  std::stringstream hname;
216  fAmbigWireU0.resize(fnEvts); fAmbigWireU1.resize(fnEvts);
217  fCheatWireU0.resize(fnEvts); fCheatWireU1.resize(fnEvts);
218  fDisambigWireU0.resize(fnEvts); fDisambigWireU1.resize(fnEvts);
219  fAmbigWireV0.resize(fnEvts); fAmbigWireV1.resize(fnEvts);
220  fCheatWireV0.resize(fnEvts); fCheatWireV1.resize(fnEvts);
221  fDisambigWireV0.resize(fnEvts); fDisambigWireV1.resize(fnEvts);
222  fTimeChanZ0.resize(fnEvts); fTimeChanZ1.resize(fnEvts);
223 
224  // art starts events at number 1
225  for(unsigned int evt=1; evt<=fnEvts; evt++){
226  // index with evt-1
227  fAmbigWireU0[evt-1].resize(fnAPAs); fAmbigWireU1[evt-1].resize(fnAPAs);
228  fCheatWireU0[evt-1].resize(fnAPAs); fCheatWireU1[evt-1].resize(fnAPAs);
229  fDisambigWireU0[evt-1].resize(fnAPAs); fDisambigWireU1[evt-1].resize(fnAPAs);
230  fAmbigWireV0[evt-1].resize(fnAPAs); fAmbigWireV1[evt-1].resize(fnAPAs);
231  fCheatWireV0[evt-1].resize(fnAPAs); fCheatWireV1[evt-1].resize(fnAPAs);
232  fDisambigWireV0[evt-1].resize(fnAPAs); fDisambigWireV1[evt-1].resize(fnAPAs);
233  fTimeChanZ0[evt-1].resize(fnAPAs); fTimeChanZ1[evt-1].resize(fnAPAs);
234 
235 
236  for(unsigned int apa=0; apa<fnAPAs; apa++){
237  //unsigned int fUChanMin = fUChMin + apa*fChansPerAPA;
238  //unsigned int fUChanMax = fUChMax + apa*fChansPerAPA;
239  //unsigned int fVChanMin = fVChMin + apa*fChansPerAPA;
240  //unsigned int fVChanMax = fVChMax + apa*fChansPerAPA;
241  unsigned int fZ0ChanMin = fZ0ChMin + apa*fChansPerAPA;
242  unsigned int fZ0ChanMax = fZ0ChMax + apa*fChansPerAPA;
243  unsigned int fZ1ChanMin = fZ1ChMin + apa*fChansPerAPA;
244  unsigned int fZ1ChanMax = fZ1ChMax + apa*fChansPerAPA;
245 
246  hname.str(""); hname << "fAmbigWireU0_" << apa << "_" << evt;
247  fAmbigWireU0[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "U Ambiguous Hits (side 0)", fUWireMax, 0, fUWireMax, binT, minT, maxT);
248  hname.str(""); hname << "fCheatWireU0_" << apa << "_" << evt;
249  fCheatWireU0[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "U Cheated Hits (side 0)", fUWireMax, 0, fUWireMax, binT, minT, maxT);
250  hname.str(""); hname << "fDisambigWireU0_" << apa << "_" << evt;
251  fDisambigWireU0[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "U Disambiguated Hits (side 0)", fUWireMax, 0, fUWireMax, binT, minT, maxT);
252 
253  hname.str(""); hname << "fAmbigWireV0_" << apa << "_" << evt;
254  fAmbigWireV0[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "V Ambiguous Hits (side 0)", fVWireMax, 0, fVWireMax, binT, minT, maxT);
255  hname.str(""); hname << "fCheatWireV0_" << apa << "_" << evt;
256  fCheatWireV0[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "V Cheated Hits (side 0)", fVWireMax, 0, fVWireMax, binT, minT, maxT);
257  hname.str(""); hname << "fDisambigWireV0_" << apa << "_" << evt;
258  fDisambigWireV0[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "V Disambiguated Hits (side 0)", fVWireMax, 0, fVWireMax, binT, minT, maxT);
259 
260  hname.str(""); hname << "fAmbigWireU1_" << apa << "_" << evt;
261  fAmbigWireU1[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "U Ambiguous Hits (side 1)", fUWireMax, 0, fUWireMax, binT, minT, maxT);
262  hname.str(""); hname << "fCheatWireU1_" << apa << "_" << evt;
263  fCheatWireU1[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "U Cheated Hits (side 1)", fUWireMax, 0, fUWireMax, binT, minT, maxT);
264  hname.str(""); hname << "fDisambigWireU1_" << apa << "_" << evt;
265  fDisambigWireU1[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "U Disambiguated Hits (side 1)", fUWireMax, 0, fUWireMax, binT, minT, maxT);
266 
267  hname.str(""); hname << "fAmbigWireV1_" << apa << "_" << evt;
268  fAmbigWireV1[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "V Ambiguous Hits (side 1)", fVWireMax, 0, fVWireMax, binT, minT, maxT);
269  hname.str(""); hname << "fCheatWireV1_" << apa << "_" << evt;
270  fCheatWireV1[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "V Cheated Hits (side 1)", fVWireMax, 0, fVWireMax, binT, minT, maxT);
271  hname.str(""); hname << "fDisambigWireV1_" << apa << "_" << evt;
272  fDisambigWireV1[evt-1][apa] = tfs->make<TH2I>(hname.str().c_str(), "V Disambiguated Hits (side 1)", fVWireMax, 0, fVWireMax, binT, minT, maxT);
273 
274 
275  fAmbigWireU0[evt-1][apa]->SetStats(0); fAmbigWireV0[evt-1][apa]->SetStats(0);
276  fAmbigWireU1[evt-1][apa]->SetStats(0); fAmbigWireV1[evt-1][apa]->SetStats(0);
277  fDisambigWireU0[evt-1][apa]->SetStats(0); fDisambigWireV0[evt-1][apa]->SetStats(0);
278  fDisambigWireU1[evt-1][apa]->SetStats(0); fDisambigWireV1[evt-1][apa]->SetStats(0);
279  fCheatWireU0[evt-1][apa]->SetStats(0); fCheatWireV0[evt-1][apa]->SetStats(0);
280  fCheatWireU1[evt-1][apa]->SetStats(0); fCheatWireV1[evt-1][apa]->SetStats(0);
281 
282  fDisambigWireU0[evt-1][apa]->SetMarkerColor(2); fDisambigWireV0[evt-1][apa]->SetMarkerColor(2);
283  fDisambigWireU1[evt-1][apa]->SetMarkerColor(2); fDisambigWireV1[evt-1][apa]->SetMarkerColor(2);
284 
285  hname.str(""); hname << "fTimeChanZ0_" << apa << "_" << evt;
286  fTimeChanZ0[evt-1][apa] = tfs->make<TH2I>( hname.str().c_str(), "",
287  (fZ0ChanMax - fZ0ChanMin + 1)*ThumbFactor,
288  fZ0ChanMin, fZ0ChanMax,
289  binT, minZ0T, maxZ0T); // NOTE: Time reversed with a fcl trigger
290 
291  hname.str(""); hname << "fTimeChanZ1_" << apa << "_" << evt;
292  fTimeChanZ1[evt-1][apa] = tfs->make<TH2I>( hname.str().c_str(), "",
293  (fZ1ChanMax - fZ1ChanMin + 1)*ThumbFactor,
294  fZ1ChanMin, fZ1ChanMax, binT, minT, maxT);
295 
296  fTimeChanZ0[evt-1][apa]->SetStats(0); fTimeChanZ1[evt-1][apa]->SetStats(0);
297 
298  } // apas
299  } // evts
300 
301 
302  }
AnalysisExample::DbigEVD4apaFD::fAmbigWireU0
std::vector< std::vector< TH2I * > > fAmbigWireU0
Definition: DbigEVD4apaFD_module.cc:74
geo::kV
Planes which measure V.
Definition: geo_types.h:126
AnalysisExample::DbigEVD4apaFD::fDisambigWireU1
std::vector< std::vector< TH2I * > > fDisambigWireU1
Definition: DbigEVD4apaFD_module.cc:80
geo::GeometryCore::ChannelToWire
std::vector< geo::WireID > ChannelToWire(raw::ChannelID_t const channel) const
Returns a list of wires connected to the specified TPC channel.
Definition: GeometryCore.cxx:1121
geo::kZ
Planes which measure Z direction.
Definition: geo_types.h:128
AnalysisExample::DbigEVD4apaFD::fCheatWireV1
std::vector< std::vector< TH2I * > > fCheatWireV1
Definition: DbigEVD4apaFD_module.cc:85
geo::GeometryCore::Nwires
unsigned int Nwires(unsigned int p, unsigned int tpc=0, unsigned int cstat=0) const
Returns the total number of wires in the specified plane.
Definition: GeometryCore.h:3347
AnalysisExample::DbigEVD4apaFD::fAmbigWireU1
std::vector< std::vector< TH2I * > > fAmbigWireU1
Definition: DbigEVD4apaFD_module.cc:76
geo::GeometryCore::Nchannels
unsigned int Nchannels() const
Returns the number of TPC readout channels in the detector.
Definition: GeometryCore.cxx:210
geo::kU
Planes which measure U.
Definition: geo_types.h:125
AnalysisExample::DbigEVD4apaFD::fTimeChanZ0
std::vector< std::vector< TH2I * > > fTimeChanZ0
Definition: DbigEVD4apaFD_module.cc:86
AnalysisExample::DbigEVD4apaFD::fDisambigWireU0
std::vector< std::vector< TH2I * > > fDisambigWireU0
Definition: DbigEVD4apaFD_module.cc:78
AnalysisExample::DbigEVD4apaFD::fCheatWireV0
std::vector< std::vector< TH2I * > > fCheatWireV0
Definition: DbigEVD4apaFD_module.cc:83
AnalysisExample::DbigEVD4apaFD::fAmbigWireV0
std::vector< std::vector< TH2I * > > fAmbigWireV0
Definition: DbigEVD4apaFD_module.cc:75
AnalysisExample::DbigEVD4apaFD::fGeom
art::ServiceHandle< geo::Geometry > fGeom
Definition: DbigEVD4apaFD_module.cc:127
geo::GeometryCore::View
View_t View(geo::PlaneID const &pid) const
Returns the view (wire orientation) on the channels of specified TPC plane.
Definition: GeometryCore.cxx:603
AnalysisExample::DbigEVD4apaFD::fCheatWireU1
std::vector< std::vector< TH2I * > > fCheatWireU1
Definition: DbigEVD4apaFD_module.cc:84
geo::GeometryCore::NTPC
unsigned int NTPC(unsigned int cstat=0) const
Returns the total number of TPCs in the specified cryostat.
Definition: GeometryCore.h:2321
AnalysisExample::DbigEVD4apaFD::fDisambigWireV1
std::vector< std::vector< TH2I * > > fDisambigWireV1
Definition: DbigEVD4apaFD_module.cc:81
AnalysisExample::DbigEVD4apaFD::fAmbigWireV1
std::vector< std::vector< TH2I * > > fAmbigWireV1
Definition: DbigEVD4apaFD_module.cc:77
tca::evt
TCEvent evt
Definition: DataStructs.cxx:7
AnalysisExample::DbigEVD4apaFD::fDisambigWireV0
std::vector< std::vector< TH2I * > > fDisambigWireV0
Definition: DbigEVD4apaFD_module.cc:79
AnalysisExample::DbigEVD4apaFD::fTimeChanZ1
std::vector< std::vector< TH2I * > > fTimeChanZ1
Definition: DbigEVD4apaFD_module.cc:87
AnalysisExample::DbigEVD4apaFD::fCheatWireU0
std::vector< std::vector< TH2I * > > fCheatWireU0
Definition: DbigEVD4apaFD_module.cc:82
void AnalysisExample::DbigEVD4apaFD::beginRun ( const art::Run &  run)

Definition at line 306 of file DbigEVD4apaFD_module.cc.

306  {
307 
308  }
void AnalysisExample::DbigEVD4apaFD::reconfigure ( fhicl::ParameterSet const &  pset)

Definition at line 144 of file DbigEVD4apaFD_module.cc.

144  {
145  fChanHitLabel = p.get< std::string >("ChanHitLabel");
146  fWidHitLabel = p.get< std::string >("WidHitLabel");
147  fHitCheatLabel = p.get< std::string >("HitCheatLabel");
148  fClusLabel = p.get< std::string >("ClusLabel");
149  fnEvts = p.get< unsigned int >("nEvents");
150 
151  fTimeReverseZ0 = p.get< bool >("TimeReverseZ0");
152 
153  auto const *fDetProp = lar::providerFrom<detinfo::DetectorPropertiesService>();
154  fNticks = fDetProp->NumberTimeSamples();
155  fnAPAs = fGeom->NTPC()/2;
156  return;
157  }
string
std::string string
Definition: nybbler.cc:12
AnalysisExample::DbigEVD4apaFD::fGeom
art::ServiceHandle< geo::Geometry > fGeom
Definition: DbigEVD4apaFD_module.cc:127
geo::GeometryCore::NTPC
unsigned int NTPC(unsigned int cstat=0) const
Returns the total number of TPCs in the specified cryostat.
Definition: GeometryCore.h:2321
test.p
p
Definition: test.py:223

Member Data Documentation

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fAmbigWireU0
private

Definition at line 74 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fAmbigWireU1
private

Definition at line 76 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fAmbigWireV0
private

Definition at line 75 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fAmbigWireV1
private

Definition at line 77 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fAPA
private

Definition at line 66 of file DbigEVD4apaFD_module.cc.

std::map< unsigned int, unsigned int > AnalysisExample::DbigEVD4apaFD::fAPAToNumHits
private

Definition at line 71 of file DbigEVD4apaFD_module.cc.

std::string AnalysisExample::DbigEVD4apaFD::fChanHitLabel
private

Definition at line 62 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fChansPerAPA
private

Definition at line 114 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fCheatWireU0
private

Definition at line 82 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fCheatWireU1
private

Definition at line 84 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fCheatWireV0
private

Definition at line 83 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fCheatWireV1
private

Definition at line 85 of file DbigEVD4apaFD_module.cc.

std::string AnalysisExample::DbigEVD4apaFD::fClusLabel
private

Definition at line 65 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fDisambigWireU0
private

Definition at line 78 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fDisambigWireU1
private

Definition at line 80 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fDisambigWireV0
private

Definition at line 79 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fDisambigWireV1
private

Definition at line 81 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fEvent
private

Definition at line 67 of file DbigEVD4apaFD_module.cc.

int AnalysisExample::DbigEVD4apaFD::fFirstEvt
private

Definition at line 125 of file DbigEVD4apaFD_module.cc.

art::ServiceHandle<geo::Geometry> AnalysisExample::DbigEVD4apaFD::fGeom
private

Definition at line 127 of file DbigEVD4apaFD_module.cc.

bool AnalysisExample::DbigEVD4apaFD::fGotFirstEvtNum = false
private

Definition at line 124 of file DbigEVD4apaFD_module.cc.

std::string AnalysisExample::DbigEVD4apaFD::fHitCheatLabel
private

Definition at line 64 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fMaxT
private

Definition at line 121 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fMaxTimeRange
private

Definition at line 121 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fMinT
private

Definition at line 121 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fnAPAs
private

Definition at line 112 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fnEvts
private

Definition at line 113 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fNticks
private

Definition at line 100 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fTimeChanZ0
private

Definition at line 86 of file DbigEVD4apaFD_module.cc.

std::vector<std::vector<TH2I*> > AnalysisExample::DbigEVD4apaFD::fTimeChanZ1
private

Definition at line 87 of file DbigEVD4apaFD_module.cc.

bool AnalysisExample::DbigEVD4apaFD::fTimeReverseZ0 = true
private

Definition at line 122 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fUChMax
private

Definition at line 104 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fUChMin
private

Definition at line 103 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fUWireMax
private

Definition at line 116 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fVChMax
private

Definition at line 106 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fVChMin
private

Definition at line 105 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fVWireMax
private

Definition at line 117 of file DbigEVD4apaFD_module.cc.

std::string AnalysisExample::DbigEVD4apaFD::fWidHitLabel
private

Definition at line 63 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fZ0ChMax
private

Definition at line 108 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fZ0ChMin
private

Definition at line 107 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fZ0WireMax
private

Definition at line 118 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fZ1ChMax
private

Definition at line 110 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fZ1ChMin
private

Definition at line 109 of file DbigEVD4apaFD_module.cc.

unsigned int AnalysisExample::DbigEVD4apaFD::fZ1WireMax
private

Definition at line 119 of file DbigEVD4apaFD_module.cc.

The documentation for this class was generated from the following file:
Generated by   doxygen 1.8.11