LArSoft  v06_85_00
Liquid Argon Software toolkit - http://larsoft.org/
caldata::CalWireAna Class Reference

Base class for creation of raw signals on wires. More...

Inheritance diagram for caldata::CalWireAna:
art::EDAnalyzer art::EventObserverBase art::Consumer art::EngineCreator

Public Types

using WorkerType = WorkerT< EDAnalyzer >
 
using ModuleType = EDAnalyzer
 

Public Member Functions

 CalWireAna (fhicl::ParameterSet const &pset)
 
virtual ~CalWireAna ()
 
void analyze (const art::Event &evt)
 read/write access to event More...
 
void beginJob ()
 
void endJob ()
 
std::string workerType () const
 
bool modifiesEvent () const
 
void registerProducts (MasterProductRegistry &, ProductDescriptions &, ModuleDescription const &)
 
std::string const & processName () const
 
bool wantAllEvents () const
 
bool wantEvent (Event const &e)
 
fhicl::ParameterSetID selectorConfig () const
 
art::Handle< art::TriggerResultsgetTriggerResults (Event const &e) const
 
template<typename T , BranchType = InEvent>
ProductToken< T > consumes (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< T > consumes (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< T > consumesView (InputTag const &it)
 
template<typename T , BranchType = InEvent>
ProductToken< T > mayConsume (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< T > mayConsume (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< T > mayConsumeView (InputTag const &it)
 
base_engine_tcreateEngine (seed_t seed)
 
base_engine_tcreateEngine (seed_t seed, std::string const &kind_of_engine_to_make)
 
base_engine_tcreateEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
 
seed_t get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)
 

Static Public Member Functions

static cet::exempt_ptr< Consumernon_module_context ()
 

Protected Member Functions

CurrentProcessingContext const * currentContext () const
 
detail::CachedProducts & cachedProducts ()
 
void validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
 
void prepareForJob (fhicl::ParameterSet const &pset)
 
void showMissingConsumes () const
 

Private Attributes

std::string fCalWireModuleLabel
 name of module that produced the wires More...
 
std::string fDetSimModuleLabel
 
TH1F * fDiffsR
 histogram of Raw tdc to tdc differences More...
 
TH1F * fDiffsW
 histogram of Wire (post-deconvoution) tdc to tdc differences More...
 
TH1F * fDiffsRW
 
TH1F * fDiffsRWgph
 
TH1F * fMin
 
TH1F * fMax
 
TH1F * fIR
 
TH1F * fCR
 
TH1F * fIW
 
TH1F * fCW
 
TH1F * fRawIndPeak
 
TH1F * fRawColPeak
 
TH1F * fCalIndPeak
 
TH1F * fCalColPeak
 
TH1F * fNoiseHist
 
TH1F * fNoiseRMS
 
TH2F * fWireSig
 
TH2F * fRawSig
 
TH2F * fRD_WireMeanDiff2D
 histogram of difference between original tdc value and compressesed value vs original value More...
 
TH2F * fRD_WireRMSDiff2D
 histogram of difference between original tdc value and compressesed value vs original value More...
 
TH2F * fDiffsRWvsR
 
TH2F * fDiffsRWvsRgph
 
TH2F * fWindow
 

Detailed Description

Base class for creation of raw signals on wires.

Definition at line 50 of file CalWireAna_module.cc.

Member Typedef Documentation

Definition at line 39 of file EDAnalyzer.h.

Definition at line 38 of file EDAnalyzer.h.

Constructor & Destructor Documentation

caldata::CalWireAna::CalWireAna ( fhicl::ParameterSet const &  pset)
explicit

Definition at line 100 of file CalWireAna_module.cc.

101  : EDAnalyzer (pset)
102  , fCalWireModuleLabel(pset.get< std::string >("CalWireModuleLabel"))
103  , fDetSimModuleLabel (pset.get< std::string >("DetSimModuleLabel"))
104  {
105 
106  }
EDAnalyzer(Table< Config > const &config)
Definition: EDAnalyzer.h:100
std::string fCalWireModuleLabel
name of module that produced the wires
std::string fDetSimModuleLabel
caldata::CalWireAna::~CalWireAna ( )
virtual

Definition at line 109 of file CalWireAna_module.cc.

110  {
111  }

Member Function Documentation

void caldata::CalWireAna::analyze ( const art::Event evt)

read/write access to event

loop over all the raw digits in a window around peak

Definition at line 156 of file CalWireAna_module.cc.

References util::LArFFT::AlignedSum(), DEFINE_ART_MODULE, util::LArFFT::DoFFT(), fCalColPeak, fCalIndPeak, fCalWireModuleLabel, fCR, fCW, fDetSimModuleLabel, fDiffsR, fDiffsRW, fDiffsRWgph, fDiffsRWvsR, fDiffsRWvsRgph, fDiffsW, util::LArFFT::FFTSize(), fIR, fIW, fMax, fMin, fNoiseHist, fNoiseRMS, fRawColPeak, fRawIndPeak, fRawSig, fRD_WireMeanDiff2D, fRD_WireRMSDiff2D, fWindow, fWireSig, art::DataViewImpl::getByLabel(), geo::kCollection, geo::kInduction, LOG_DEBUG, max, min, art::PtrVector< T >::push_back(), geo::GeometryCore::SignalType(), art::PtrVector< T >::size(), tmp, and w.

157  {
158 
159  // loop over the raw digits and get the adc vector for each, then compress it and uncompress it
160 
161  lariov::ChannelStatusProvider const& channelStatus
164  evt.getByLabel(fDetSimModuleLabel,rdHandle);
166  evt.getByLabel(fCalWireModuleLabel,wHandle);
167 
168  // get the raw::RawDigit associated by fCalWireModuleLabel to wires in wHandle;
169  // RawDigitsFromWire.at(index) will be a art::Ptr<raw::RawDigit>
170  art::FindOneP<raw::RawDigit> RawDigitsFromWire(wHandle, evt, fCalWireModuleLabel);
171 
173  for(unsigned int i = 0; i < wHandle->size(); ++i){
174  art::Ptr<recob::Wire> w(wHandle,i);
175  wvec.push_back(w);
176  }
178  for(unsigned int i = 0; i < rdHandle->size(); ++i){
179  art::Ptr<raw::RawDigit> r(rdHandle,i);
180  rdvec.push_back(r);
181  }
182 
185  double pedestal = rdvec[0]->GetPedestal();
186  double threshold = 9.0;
187  double signalSize = rdvec[0]->Samples();
188  uint32_t indChan0=64;
189  uint32_t indChan1=110;
190  uint32_t colChan0=312;
191  uint32_t colChan1=354;
192  std::vector<double> ir(fft->FFTSize()),iw(fft->FFTSize()),cr(fft->FFTSize()),cw(fft->FFTSize());
193  ir[signalSize]=iw[signalSize]=cr[signalSize]=cw[signalSize]=1.0;
195  for(unsigned int rd = 0; rd < rdvec.size(); ++rd){
196  // Find corresponding wire.
197  std::vector<double> signal(fft->FFTSize());
198  for(unsigned int wd = 0; wd < wvec.size(); ++wd){
199  // if (wvec[wd]->RawDigit() == rdvec[rd]){
200  if (RawDigitsFromWire.at(wd) == rdvec[rd]){
201  std::vector<float> wirSig = wvec[wd]->Signal();
202  if(wirSig.size() > signal.size()) {
203  LOG_DEBUG("CalWireAna")<<"Incompatible vector size "<<wirSig.size()
204  <<" "<<signal.size();
205  return;
206  }
207  for(unsigned int ii = 0; ii < wirSig.size(); ++ii) {
208  signal[ii] = wirSig[ii];
209  }
210  break;
211  }
212  if (wd == (wvec.size()-1) ){
213  LOG_DEBUG("CalWireAna") << "caldata::CalWireAna:Big problem! No matching Wire for RawDigit. Bailing." << rd;
214  return;
215  }
216  }
217 
218  std::vector<double> adc(fft->FFTSize());
219 
220  for(unsigned int t = 1; t < rdvec[rd]->Samples(); ++t){
221  fDiffsR->Fill(rdvec[rd]->ADC(t) - rdvec[rd]->ADC(t-1));
222  adc[t-1]=rdvec[rd]->ADC(t-1);
223  fRawSig->Fill(rd,rdvec[rd]->ADC(t));
224  }
225  //get the last one for the adc vector
226  adc[rdvec[rd]->Samples()-1] = rdvec[rd]->ADC(rdvec[rd]->Samples()-1);
227  if(!channelStatus.IsBad(rdvec[rd]->Channel()) &&
228  (*max_element(adc.begin(),adc.end()) < pedestal+threshold &&
229  *min_element(adc.begin(),adc.end()) >pedestal -threshold)) {
230  double sum=0;
231  for(int i = 0; i < signalSize; i++) sum+=pow(adc[i]-pedestal,2.0);
232  fNoiseRMS->Fill(TMath::Sqrt(sum/(double)signalSize));
233  std::vector<double> temp(fft->FFTSize());
234  std::vector<TComplex> fTemp(fft->FFTSize()/2+1);
235  for(int i = 0; i < signalSize; i++) temp[i]=(adc[i]-pedestal)*sin(TMath::Pi()*(double)i/signalSize);
236  fft->DoFFT(temp,fTemp);
237  for(int i = 0; i < fft->FFTSize()/2+1; i++) fNoiseHist->Fill(i,fTemp[i].Rho());
238  }
239  if(geom->SignalType(rdvec[rd]->Channel()) == geo::kInduction &&
240  rdvec[rd]->Channel() > indChan0 &&
241  rdvec[rd]->Channel() < indChan1){
242  fft->AlignedSum(ir,adc);
243  fft->AlignedSum(iw,signal);
244  }
245  if(geom->SignalType(rdvec[rd]->Channel()) == geo::kCollection &&
246  rdvec[rd]->Channel() > colChan0 &&
247  rdvec[rd]->Channel() < colChan1) {
248  fft->AlignedSum(cr,adc);
249  fft->AlignedSum(cw,signal);
250  }
251  if(geom->SignalType(rdvec[rd]->Channel()) == geo::kInduction) {
252  if(*max_element(adc.begin(),adc.end()) > pedestal+threshold)
253  fRawIndPeak->Fill(*max_element(adc.begin(),adc.end()));
254  if(*max_element(signal.begin(),signal.end()) > pedestal+threshold)
255  fCalIndPeak->Fill(*max_element(signal.begin(),signal.end()));
256  }
257  if(geom->SignalType(rdvec[rd]->Channel()) == geo::kCollection) {
258  if(*max_element(adc.begin(),adc.end()) > pedestal +threshold)
259  fRawColPeak->Fill(*max_element(adc.begin(),adc.end()));
260  if(*max_element(signal.begin(),signal.end()) > pedestal+threshold)
261  fCalColPeak->Fill(*max_element(signal.begin(),signal.end()));
262  }
263 
264  int window = 8;
265  static unsigned int pulseHeight = 5;
266  unsigned int tmin = 1;
267  unsigned int tmax = 1;
268  int indMax = TMath::LocMax(signalSize,&adc[0]);
269  double sigMin = 0.0;
270  double sigMax = TMath::MaxElement(signalSize,&adc[0]);
271  if(geom->SignalType(rdvec[rd]->Channel()) == geo::kInduction && sigMax>=pulseHeight) {
272  int indMin = TMath::LocMin(signalSize,&adc[0]);
273  sigMin = TMath::MinElement(signalSize,&adc[0]);
274  tmin = std::max(indMax-window,0);
275  tmax = std::min(indMin+window,(int)signalSize-1);
276  LOG_DEBUG("CalWireAna") << "Induction channel, indMin,tmin,tmax "
277  << rd<< " " << indMin<< " " << tmin << " " << tmax;
278  }
279  else if (sigMax>=pulseHeight){
280  tmin = std::max(indMax-window,0);
281  tmax = std::min(indMax+window,(int)signalSize-1);
282  LOG_DEBUG("CalWireAna") << "Collection channel, tmin,tmax "<< rd<< " " << tmin << " " << tmax;
283  }
284 
285  fMin->Fill(sigMin);
286  fMax->Fill(sigMax);
287  fWindow->Fill(indMax, tmax - tmin);
288 
289  std::vector<double> winDiffs;
290  int cnt=0;
291  // for(unsigned int t = tmin; t < tmax; ++t)
292  static unsigned int tRawLead = 0;
293  for(unsigned int t = 1; t < signalSize; ++t){
294  fDiffsW->Fill(signal[t]-signal[t-1]);
295  fWireSig->Fill(rd,signal[t]);
296 
297  if (t>=tmin && t<=tmax && tmax>=pulseHeight && (t+tRawLead)<signalSize){
298  cnt++;
299  winDiffs.push_back((adc[t+tRawLead]-signal[t])/adc[t+tRawLead]);
300  fDiffsRW->Fill(adc[t+tRawLead]-signal[t]);
301  fDiffsRWvsR->Fill(rd,adc[t+tRawLead]-signal[t]);
302  if (sigMax >= pulseHeight){
303  fDiffsRWgph->Fill(adc[t+tRawLead]-signal[t]);
304  fDiffsRWvsRgph->Fill(rd,adc[t+tRawLead]-signal[t]);
305  }
306  }
307  }
308 
309  LOG_DEBUG("CalWireAna") << "on channel " << rdvec[rd]->Channel();
310  // TMath::Mean with iterators doesn't work. EC,23-Sep-2010.
311  double tmp = TMath::Mean(winDiffs.size(),&winDiffs[0]);
312  double tmp2 = TMath::RMS(winDiffs.size(),&winDiffs[0]);
313  double tmp3=0;
314  for (int ii=0; ii<rdvec[rd]->Samples(); ii++) tmp3+=rdvec[rd]->ADC(ii);
315  for(int i = 0; i < fft->FFTSize(); i++) {
316  fIR->Fill(i,ir[i]);
317  fIW->Fill(i,iw[i]);
318  fCR->Fill(i,cr[i]);
319  fCW->Fill(i,cw[i]);
320  }
321  fRD_WireMeanDiff2D->Fill(rd,tmp);
322  fRD_WireRMSDiff2D->Fill(rd,tmp2);
323 
324  }//end loop over rawDigits
325 
326  return;
327  }//end analyze method
TH2F * fRD_WireMeanDiff2D
histogram of difference between original tdc value and compressesed value vs original value ...
Float_t tmp
Definition: plot.C:37
void DoFFT(std::vector< T > &input, std::vector< TComplex > &output)
Definition: LArFFT.h:97
SigType_t SignalType(geo::PlaneID const &pid) const
Returns the type of signal on the channels of specified TPC plane.
TH1F * fDiffsR
histogram of Raw tdc to tdc differences
int FFTSize() const
Definition: LArFFT.h:69
TH2F * fRD_WireRMSDiff2D
histogram of difference between original tdc value and compressesed value vs original value ...
Int_t max
Definition: plot.C:27
void push_back(Ptr< U > const &p)
Definition: PtrVector.h:441
Signal from induction planes.
Definition: geo_types.h:92
size_type size() const
Definition: PtrVector.h:308
std::string fCalWireModuleLabel
name of module that produced the wires
void AlignedSum(std::vector< T > &input, std::vector< T > &output, bool add=true)
Definition: LArFFT.h:243
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
Definition: DataViewImpl.h:344
Int_t min
Definition: plot.C:26
std::string fDetSimModuleLabel
TH1F * fDiffsW
histogram of Wire (post-deconvoution) tdc to tdc differences
#define LOG_DEBUG(id)
Float_t w
Definition: plot.C:23
Definition: fwd.h:25
Signal from collection planes.
Definition: geo_types.h:93
void caldata::CalWireAna::beginJob ( )
virtual

Reimplemented from art::EDAnalyzer.

Definition at line 114 of file CalWireAna_module.cc.

References fCalColPeak, fCalIndPeak, fCR, fCW, fDiffsR, fDiffsRW, fDiffsRWgph, fDiffsRWvsR, fDiffsRWvsRgph, fDiffsW, fIR, fIW, fMax, fMin, fNoiseHist, fNoiseRMS, fRawColPeak, fRawIndPeak, fRawSig, fRD_WireMeanDiff2D, fRD_WireRMSDiff2D, fWindow, fWireSig, and art::TFileDirectory::make().

115  {
116  // get access to the TFile service
118 
119  fDiffsR = tfs->make<TH1F>("One timestamp diffs in RawDigits", ";#Delta ADC;", 40, -19.5, 20.5);
120  fDiffsW = tfs->make<TH1F>("One timestamp diffs in Wires", ";#Delta ADC;", 20, -9.5, 10.5);
121  fDiffsRW = tfs->make<TH1F>("Same timestamp diffs in RD-Wires", ";#Delta ADC;", 20, -9.5, 10.5);
122  fDiffsRWvsR = tfs->make<TH2F>("Same timestamp diffs in RD-Wires vs R", ";#Delta ADC;", 481,-0.5,480.5, 20, -9.5, 10.5);
123  fDiffsRWgph = tfs->make<TH1F>("Same timestamp diffs in RD-Wires gph", ";#Delta ADC;", 20, -9.5, 10.5);
124  fDiffsRWvsRgph = tfs->make<TH2F>("Same timestamp diffs in RD-Wires vs R gph", ";#Delta ADC;", 481,-0.5,480.5, 20, -9.5, 10.5);
125  fRawSig = tfs->make<TH2F>("One event, one channel Raw", "timestamp", 481,-0.5,480.5, 21,-0.5, 20.5);
126  fWireSig = tfs->make<TH2F>("One event, one channel Wire", "timestamp", 481,-0.5,480.5, 21, -0.5, 20.5);
127 
128  fRD_WireMeanDiff2D = tfs->make<TH2F>("Mean (Raw-CALD)-over-Raw in Window gph","Wire number",481,-0.05,480.5, 40, -1., 1.);
129  fRD_WireRMSDiff2D = tfs->make<TH2F>("RMS (Raw-CALD)-over-Raw in Window gph","Wire number",481,-0.05,480.5, 10, 0., 2.);
130 
131  fWindow = tfs->make<TH2F>("tmax-tmin vs indMax", "ticks", 200, 0, 2000, 20, -2.5, 60.5);
132  fMin = tfs->make<TH1F>("Value of min", "ticks", 21, -20.5, 0.5);
133  fMax = tfs->make<TH1F>("Value of max", "ticks", 21, 0.5, 20.5);
134  fRawIndPeak = tfs->make<TH1F>("indPeakRaw", ";Induction Peaks Raw;",40,5,45);
135  fRawColPeak = tfs->make<TH1F>("colPeakRaw", ";Collection Peaks Raw;",40,5,45);
136  fCalIndPeak = tfs->make<TH1F>("indPeakCal", ";Induction Peaks Calibrated;",40,5,45);
137  fCalColPeak = tfs->make<TH1F>("colPeakCal", ";Collection Peaks Calibrated;",40,5,45);
138 
139  fIR = tfs->make<TH1F>("Raw Ind signal","time ticks",4096,0.0,4096.);
140  fCR = tfs->make<TH1F>("Raw Coll signal","time ticks",4096,0.0,4096.);
141  fIW = tfs->make<TH1F>("Wire Ind signal","time ticks",4096,0.0,4096.);
142  fCW = tfs->make<TH1F>("Wire Coll signal","time ticks",4096,0.0,4096.);
143  fNoiseHist = tfs->make<TH1F>("Noise Histogram","FFT Bins",2049,0,2049);
144  fNoiseRMS = tfs->make<TH1F>("Noise RMS","RMS",25,0,2.0);
145 
146  return;
147 
148  }
TH2F * fRD_WireMeanDiff2D
histogram of difference between original tdc value and compressesed value vs original value ...
TH1F * fDiffsR
histogram of Raw tdc to tdc differences
TH2F * fRD_WireRMSDiff2D
histogram of difference between original tdc value and compressesed value vs original value ...
T * make(ARGS...args) const
TH1F * fDiffsW
histogram of Wire (post-deconvoution) tdc to tdc differences
detail::CachedProducts& art::EventObserverBase::cachedProducts ( )
inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

Referenced by art::EDAnalyzer::doEvent(), and art::OutputModule::doWriteEvent().

80  {
81  return selectors_;
82  }
detail::CachedProducts selectors_
template<typename T , BranchType = InEvent>
ProductToken<T> art::Consumer::consumes ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it)
inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

148 {
149  if (!moduleContext_)
150  return ProductToken<T>::invalid();
151 
152  consumables_[BT].emplace_back(ConsumableType::Product,
153  TypeID{typeid(T)},
154  it.label(),
155  it.instance(),
156  it.process());
157  return ProductToken<T>{it};
158 }
static ProductToken< T > invalid()
Definition: ProductToken.h:47
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
template<typename T , art::BranchType BT>
void art::Consumer::consumesMany ( )
inherited

Definition at line 162 of file Consumer.h.

163 {
164  if (!moduleContext_)
165  return;
166 
167  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
168 }
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it)
inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

173 {
174  if (!moduleContext_)
175  return ViewToken<T>::invalid();
176 
177  consumables_[BT].emplace_back(ConsumableType::ViewElement,
178  TypeID{typeid(T)},
179  it.label(),
180  it.instance(),
181  it.process());
182  return ViewToken<T>{it};
183 }
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed,
std::string const &  kind_of_engine_to_make 
)
inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

34 {
35  return rng()->createEngine(
36  placeholder_schedule_id(), seed, kind_of_engine_to_make);
37 }
long seed
Definition: chem4.cc:68
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed,
std::string const &  kind_of_engine_to_make,
label_t const &  engine_label 
)
inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

43 {
44  return rng()->createEngine(
45  placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
46 }
long seed
Definition: chem4.cc:68
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
CurrentProcessingContext const * art::EDAnalyzer::currentContext ( ) const
protectedinherited

Definition at line 114 of file EDAnalyzer.cc.

References art::EDAnalyzer::current_context_.

Referenced by art::EDAnalyzer::workerType().

115  {
116  return current_context_.get();
117  }
CPC_exempt_ptr current_context_
Definition: EDAnalyzer.h:179
void caldata::CalWireAna::endJob ( )
virtual

Reimplemented from art::EDAnalyzer.

Definition at line 151 of file CalWireAna_module.cc.

152  {
153  }
EngineCreator::seed_t EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset,
char const  key[] = "seed",
seed_t const  implicit_seed = -1 
)
inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

52 {
53  auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
54  return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
55 }
art::Handle<art::TriggerResults> art::EventObserverBase::getTriggerResults ( Event const &  e) const
inlineinherited

Definition at line 61 of file EventObserverBase.h.

References art::detail::CachedProducts::getOneTriggerResults(), and art::EventObserverBase::selectors_.

Referenced by art::OutputModule::doWriteEvent().

62  {
64  }
detail::CachedProducts selectors_
art::Handle< art::TriggerResults > getOneTriggerResults(Event const &) const
Float_t e
Definition: plot.C:34
template<typename T , BranchType = InEvent>
ProductToken<T> art::Consumer::mayConsume ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ProductToken<T> art::Consumer::mayConsume ( InputTag const &  it)
inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

191 {
192  if (!moduleContext_)
193  return ProductToken<T>::invalid();
194 
195  consumables_[BT].emplace_back(ConsumableType::Product,
196  TypeID{typeid(T)},
197  it.label(),
198  it.instance(),
199  it.process());
200  return ProductToken<T>{it};
201 }
static ProductToken< T > invalid()
Definition: ProductToken.h:47
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
template<typename T , art::BranchType BT>
void art::Consumer::mayConsumeMany ( )
inherited

Definition at line 205 of file Consumer.h.

206 {
207  if (!moduleContext_)
208  return;
209 
210  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
211 }
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const &  it)
inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

216 {
217  if (!moduleContext_)
218  return ViewToken<T>::invalid();
219 
220  consumables_[BT].emplace_back(ConsumableType::ViewElement,
221  TypeID{typeid(T)},
222  it.label(),
223  it.instance(),
224  it.process());
225  return ViewToken<T>{it};
226 }
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
bool art::EventObserverBase::modifiesEvent ( ) const
inlineinherited

Definition at line 25 of file EventObserverBase.h.

26  {
27  return false;
28  }
void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset)
protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

90 {
91  if (!moduleContext_)
92  return;
93 
94  pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
95  for (auto& consumablesPerBranch : consumables_) {
96  cet::sort_all(consumablesPerBranch);
97  }
98 }
bool requireConsumes_
Definition: Consumer.h:137
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
std::string const& art::EventObserverBase::processName ( ) const
inlineinherited
void art::EventObserverBase::registerProducts ( MasterProductRegistry ,
ProductDescriptions ,
ModuleDescription const &   
)
inlineinherited

Definition at line 33 of file EventObserverBase.h.

36  {}
fhicl::ParameterSetID art::EventObserverBase::selectorConfig ( ) const
inlineinherited

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

Referenced by art::RootOutputFile::writeOne().

57  {
58  return selector_config_id_;
59  }
fhicl::ParameterSetID selector_config_id_
void art::Consumer::showMissingConsumes ( ) const
protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

126 {
127  if (!moduleContext_)
128  return;
129 
130  // If none of the branches have missing consumes statements, exit early.
131  if (std::all_of(cbegin(missingConsumes_),
132  cend(missingConsumes_),
133  [](auto const& perBranch) { return perBranch.empty(); }))
134  return;
135 
136  constexpr cet::HorizontalRule rule{60};
137  mf::LogPrint log{"MTdiagnostics"};
138  log << '\n'
139  << rule('=') << '\n'
140  << "The following consumes (or mayConsume) statements are missing from\n"
141  << module_context(moduleDescription_) << '\n'
142  << rule('-') << '\n';
143 
144  cet::for_all_with_index(
145  missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
146  for (auto const& pi : perBranch) {
147  log << " "
148  << assemble_consumes_statement(static_cast<BranchType>(i), pi)
149  << '\n';
150  }
151  });
152  log << rule('=');
153 }
cet::exempt_ptr< ModuleDescription const > moduleDescription_
Definition: Consumer.h:140
constexpr T pi()
Returns the constant pi (up to 35 decimal digits of precision)
bool moduleContext_
Definition: Consumer.h:136
ConsumableProductSets missingConsumes_
Definition: Consumer.h:139
void art::Consumer::validateConsumedProduct ( BranchType const  bt,
ProductInfo const &  pi 
)
protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

103 {
104  // Early exits if consumes tracking has been disabled or if the
105  // consumed product is an allowed consumable.
106  if (!moduleContext_)
107  return;
108 
109  if (cet::binary_search_all(consumables_[bt], pi))
110  return;
111 
112  if (requireConsumes_) {
114  "Consumer: an error occurred during validation of a "
115  "retrieved product\n\n")
116  << "The following consumes (or mayConsume) statement is missing from\n"
117  << module_context(moduleDescription_) << ":\n\n"
118  << " " << assemble_consumes_statement(bt, pi) << "\n\n";
119  }
120 
121  missingConsumes_[bt].insert(pi);
122 }
cet::exempt_ptr< ModuleDescription const > moduleDescription_
Definition: Consumer.h:140
bool requireConsumes_
Definition: Consumer.h:137
cet::coded_exception< errors::ErrorCodes, ExceptionDetail::translate > Exception
Definition: Exception.h:66
constexpr T pi()
Returns the constant pi (up to 35 decimal digits of precision)
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
ConsumableProductSets missingConsumes_
Definition: Consumer.h:139
bool art::EventObserverBase::wantAllEvents ( ) const
inlineinherited
bool art::EventObserverBase::wantEvent ( Event const &  e)
inlineinherited
std::string art::EDAnalyzer::workerType ( ) const
inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

110  {
111  return "WorkerT<EDAnalyzer>";
112  }

Member Data Documentation

TH1F* caldata::CalWireAna::fCalColPeak
private

Definition at line 80 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fCalIndPeak
private

Definition at line 79 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

std::string caldata::CalWireAna::fCalWireModuleLabel
private

name of module that produced the wires

Definition at line 64 of file CalWireAna_module.cc.

Referenced by analyze().

TH1F* caldata::CalWireAna::fCR
private

Definition at line 74 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fCW
private

Definition at line 76 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

std::string caldata::CalWireAna::fDetSimModuleLabel
private

Definition at line 65 of file CalWireAna_module.cc.

Referenced by analyze().

TH1F* caldata::CalWireAna::fDiffsR
private

histogram of Raw tdc to tdc differences

Definition at line 67 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fDiffsRW
private

Definition at line 69 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fDiffsRWgph
private

Definition at line 70 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH2F* caldata::CalWireAna::fDiffsRWvsR
private

Definition at line 87 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH2F* caldata::CalWireAna::fDiffsRWvsRgph
private

Definition at line 88 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fDiffsW
private

histogram of Wire (post-deconvoution) tdc to tdc differences

Definition at line 68 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fIR
private

Definition at line 73 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fIW
private

Definition at line 75 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fMax
private

Definition at line 72 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fMin
private

Definition at line 71 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fNoiseHist
private

Definition at line 81 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fNoiseRMS
private

Definition at line 82 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fRawColPeak
private

Definition at line 78 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* caldata::CalWireAna::fRawIndPeak
private

Definition at line 77 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH2F* caldata::CalWireAna::fRawSig
private

Definition at line 84 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH2F* caldata::CalWireAna::fRD_WireMeanDiff2D
private

histogram of difference between original tdc value and compressesed value vs original value

Definition at line 85 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH2F* caldata::CalWireAna::fRD_WireRMSDiff2D
private

histogram of difference between original tdc value and compressesed value vs original value

Definition at line 86 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH2F* caldata::CalWireAna::fWindow
private

Definition at line 89 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().

TH2F* caldata::CalWireAna::fWireSig
private

Definition at line 83 of file CalWireAna_module.cc.

Referenced by analyze(), and beginJob().


The documentation for this class was generated from the following file: