LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
Public Types | Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | List of all members
shower::TCShowerTemplateMaker Class Reference
Inheritance diagram for shower::TCShowerTemplateMaker:
art::EDAnalyzer art::detail::Analyzer art::detail::LegacyModule art::Observer art::ModuleBase

Public Types

using ModuleType = EDAnalyzer
 

Public Member Functions

 TCShowerTemplateMaker (fhicl::ParameterSet const &pset)
 
void doBeginJob (SharedResources const &resources)
 
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)
 
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)
 
std::unique_ptr< Worker > makeWorker (WorkerParams const &wp)
 
template<typename T , BranchType BT>
ViewToken< T > consumesView (InputTag const &tag)
 
template<typename T , BranchType BT>
ViewToken< T > mayConsumeView (InputTag const &tag)
 

Protected Member Functions

std::string const & processName () const
 
bool wantAllEvents () const noexcept
 
bool wantEvent (ScheduleID id, Event const &e) const
 
Handle< TriggerResults > getTriggerResults (Event const &e) const
 
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 ()
 

Private Member Functions

void beginJob () override
 
void analyze (const art::Event &evt) override
 
void showerProfile (detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, std::vector< art::Ptr< recob::Hit >> showerhits, TVector3 shwvtx, TVector3 shwdir, double elep)
 
void showerProfileTrue (detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, std::vector< art::Ptr< recob::Hit >> allhits, double elep)
 
void showerProfileTrue (std::vector< art::Ptr< sim::SimChannel >> allchan, simb::MCParticle electron)
 

Private Attributes

TProfile * fShowerProfileSimLong
 
TProfile * fShowerProfileHitLong
 
TProfile * fShowerProfileRecoLong
 
TProfile2D * fShowerProfileSimLong2D
 
TProfile2D * fShowerProfileHitLong2D
 
TProfile2D * fShowerProfileRecoLong2D
 
TProfile * fShowerProfileSimTrans
 
TProfile * fShowerProfileHitTrans
 
TProfile * fShowerProfileRecoTrans
 
TProfile2D * fShowerProfileSimTrans2D
 
TProfile2D * fShowerProfileHitTrans2D
 
TProfile2D * fShowerProfileRecoTrans2D
 
TProfile2D * fShowerProfileSimTrans2D_1
 
TProfile2D * fShowerProfileHitTrans2D_1
 
TProfile2D * fShowerProfileRecoTrans2D_1
 
TProfile2D * fShowerProfileSimTrans2D_2
 
TProfile2D * fShowerProfileHitTrans2D_2
 
TProfile2D * fShowerProfileRecoTrans2D_2
 
TProfile2D * fShowerProfileSimTrans2D_3
 
TProfile2D * fShowerProfileHitTrans2D_3
 
TProfile2D * fShowerProfileRecoTrans2D_3
 
TProfile2D * fShowerProfileSimTrans2D_4
 
TProfile2D * fShowerProfileHitTrans2D_4
 
TProfile2D * fShowerProfileRecoTrans2D_4
 
TProfile2D * fShowerProfileSimTrans2D_5
 
TProfile2D * fShowerProfileHitTrans2D_5
 
TProfile2D * fShowerProfileRecoTrans2D_5
 
TH3F * fLongitudinal
 
TH3F * fTransverse
 
TH3F * fTransverse_1
 
TH3F * fTransverse_2
 
TH3F * fTransverse_3
 
TH3F * fTransverse_4
 
TH3F * fTransverse_5
 
TProfile * fLongitudinal_electron
 
TProfile * fTransverse1_electron
 
TProfile * fTransverse2_electron
 
TProfile * fTransverse3_electron
 
TProfile * fTransverse4_electron
 
TProfile * fTransverse5_electron
 
TProfile * fLongitudinal_photon
 
TProfile * fTransverse1_photon
 
TProfile * fTransverse2_photon
 
TProfile * fTransverse3_photon
 
TProfile * fTransverse4_photon
 
TProfile * fTransverse5_photon
 
TProfile * fLongitudinal_other
 
TProfile * fTransverse1_other
 
TProfile * fTransverse2_other
 
TProfile * fTransverse3_other
 
TProfile * fTransverse4_other
 
TProfile * fTransverse5_other
 
const int LBINS = 20
 
const int LMIN = 0
 
const int LMAX = 5
 
const int TBINS = 20
 
const int TMIN = -5
 
const int TMAX = 5
 
const int EBINS = 20
 
const double EMIN = 0.5
 
const double EMAX = 20.5
 
const double X0 = 14
 
std::string fHitModuleLabel
 
std::string fShowerModuleLabel
 
std::string fGenieGenModuleLabel
 
std::string fDigitModuleLabel
 
calo::CalorimetryAlg fCalorimetryAlg
 

Detailed Description

Definition at line 37 of file TCShowerTemplateMaker_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer
inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

shower::TCShowerTemplateMaker::TCShowerTemplateMaker ( fhicl::ParameterSet const &  pset)
explicit

Definition at line 151 of file TCShowerTemplateMaker_module.cc.

152  : EDAnalyzer(pset)
153  , fHitModuleLabel(pset.get<std::string>("HitModuleLabel", "trajcluster"))
154  , fShowerModuleLabel(pset.get<std::string>("ShowerModuleLabel", "tcshower"))
155  , fGenieGenModuleLabel(pset.get<std::string>("GenieGenModuleLabel", "generator"))
156  , fDigitModuleLabel(pset.get<std::string>("DigitModuleLabel", "largeant"))
157  , fCalorimetryAlg(pset.get<fhicl::ParameterSet>("CalorimetryAlg"))
158 {} // TCShowerTemplateMaker
art::EDAnalyzer::EDAnalyzer
EDAnalyzer(fhicl::ParameterSet const &pset)
Definition: EDAnalyzer.cc:6

Member Function Documentation

void shower::TCShowerTemplateMaker::analyze ( const art::Event &  evt)
overrideprivate

Definition at line 577 of file TCShowerTemplateMaker_module.cc.

References util::abs(), simb::MCNeutrino::CCNC(), simb::MCParticle::E(), util::empty(), fDigitModuleLabel, fGenieGenModuleLabel, fHitModuleLabel, art::fill_ptr_vector(), fShowerModuleLabel, art::ProductRetriever::getByLabel(), simb::MCTruth::GetNeutrino(), simb::MCNeutrino::Lepton(), simb::MCTruth::NeutrinoSet(), simb::MCNeutrino::Nu(), simb::MCParticle::PdgCode(), showerProfile(), and showerProfileTrue().

578 {
579 
580  art::Handle<std::vector<recob::Hit>> hitListHandle;
581  std::vector<art::Ptr<recob::Hit>> hitlist;
582  if (evt.getByLabel(fHitModuleLabel, hitListHandle)) art::fill_ptr_vector(hitlist, hitListHandle);
583 
584  art::Handle<std::vector<sim::SimChannel>> scListHandle;
585  std::vector<art::Ptr<sim::SimChannel>> simchanlist;
586  if (evt.getByLabel(fDigitModuleLabel, scListHandle))
587  art::fill_ptr_vector(simchanlist, scListHandle);
588 
589  art::Handle<std::vector<recob::Shower>> showerListHandle;
590  std::vector<art::Ptr<recob::Shower>> showerlist;
591  if (evt.getByLabel(fShowerModuleLabel, showerListHandle))
592  art::fill_ptr_vector(showerlist, showerListHandle);
593 
594  art::Handle<std::vector<simb::MCTruth>> mctruthListHandle;
595  std::vector<art::Ptr<simb::MCTruth>> mclist;
596  if (evt.getByLabel(fGenieGenModuleLabel, mctruthListHandle))
597  art::fill_ptr_vector(mclist, mctruthListHandle);
598 
599  art::FindManyP<recob::Hit> shwfm(showerListHandle, evt, fShowerModuleLabel);
600 
601  if (empty(mclist)) return;
602 
603  auto const clock_data = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
604  auto const det_prop =
605  art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clock_data);
606 
607  art::Ptr<simb::MCTruth> mctruth = mclist[0];
608  if (mctruth->NeutrinoSet()) {
609  if (std::abs(mctruth->GetNeutrino().Nu().PdgCode()) == 12 &&
610  mctruth->GetNeutrino().CCNC() == 0) {
611  double elep = mctruth->GetNeutrino().Lepton().E();
612  if (showerlist.size()) {
613  std::vector<art::Ptr<recob::Hit>> showerhits = shwfm.at(0);
614  showerProfile(clock_data,
615  det_prop,
616  showerhits,
617  showerlist[0]->ShowerStart(),
618  showerlist[0]->Direction(),
619  elep);
620  }
621  showerProfileTrue(clock_data, det_prop, hitlist, elep);
622  showerProfileTrue(simchanlist, mctruth->GetNeutrino().Lepton());
623  }
624  }
625 } // analyze
simb::MCParticle::E
double E(const int i=0) const
Definition: MCParticle.h:234
simb::MCParticle::PdgCode
int PdgCode() const
Definition: MCParticle.h:213
simb::MCNeutrino::CCNC
int CCNC() const
Definition: MCNeutrino.h:148
simb::MCTruth::GetNeutrino
const simb::MCNeutrino & GetNeutrino() const
Definition: MCTruth.h:77
simb::MCNeutrino::Nu
const simb::MCParticle & Nu() const
Definition: MCNeutrino.h:146
util::abs
constexpr auto abs(T v)
Returns the absolute value of the argument.
Definition: constexpr_math.h:40
shower::TCShowerTemplateMaker::showerProfile
void showerProfile(detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, std::vector< art::Ptr< recob::Hit >> showerhits, TVector3 shwvtx, TVector3 shwdir, double elep)
Definition: TCShowerTemplateMaker_module.cc:629
shower::TCShowerTemplateMaker::showerProfileTrue
void showerProfileTrue(detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, std::vector< art::Ptr< recob::Hit >> allhits, double elep)
Definition: TCShowerTemplateMaker_module.cc:740
simb::MCNeutrino::Lepton
const simb::MCParticle & Lepton() const
Definition: MCNeutrino.h:147
art::ProductRetriever::getByLabel
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
Definition: ProductRetriever.h:444
art::Direction
Direction
Definition: types.h:12
simb::MCTruth::NeutrinoSet
bool NeutrinoSet() const
Definition: MCTruth.h:78
art::fill_ptr_vector
void fill_ptr_vector(std::vector< Ptr< T >> &ptrs, H const &h)
Definition: Ptr.h:306
util::empty
decltype(auto) constexpr empty(T &&obj)
ADL-aware version of std::empty.
Definition: StdUtils.h:109
void shower::TCShowerTemplateMaker::beginJob ( )
overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 162 of file TCShowerTemplateMaker_module.cc.

References EBINS, EMAX, EMIN, fLongitudinal, fLongitudinal_electron, fLongitudinal_other, fLongitudinal_photon, fShowerProfileHitLong, fShowerProfileHitLong2D, fShowerProfileHitTrans, fShowerProfileHitTrans2D, fShowerProfileHitTrans2D_1, fShowerProfileHitTrans2D_2, fShowerProfileHitTrans2D_3, fShowerProfileHitTrans2D_4, fShowerProfileHitTrans2D_5, fShowerProfileRecoLong, fShowerProfileRecoLong2D, fShowerProfileRecoTrans, fShowerProfileRecoTrans2D, fShowerProfileRecoTrans2D_1, fShowerProfileRecoTrans2D_2, fShowerProfileRecoTrans2D_3, fShowerProfileRecoTrans2D_4, fShowerProfileRecoTrans2D_5, fShowerProfileSimLong, fShowerProfileSimLong2D, fShowerProfileSimTrans, fShowerProfileSimTrans2D, fShowerProfileSimTrans2D_1, fShowerProfileSimTrans2D_2, fShowerProfileSimTrans2D_3, fShowerProfileSimTrans2D_4, fShowerProfileSimTrans2D_5, fTransverse, fTransverse1_electron, fTransverse1_other, fTransverse1_photon, fTransverse2_electron, fTransverse2_other, fTransverse2_photon, fTransverse3_electron, fTransverse3_other, fTransverse3_photon, fTransverse4_electron, fTransverse4_other, fTransverse4_photon, fTransverse5_electron, fTransverse5_other, fTransverse5_photon, fTransverse_1, fTransverse_2, fTransverse_3, fTransverse_4, fTransverse_5, LBINS, LMAX, LMIN, TBINS, TMAX, and TMIN.

163 {
164 
165  art::ServiceHandle<art::TFileService const> tfs;
166  fShowerProfileSimLong =
167  tfs->make<TProfile>("fShowerProfileSimLong",
168  "longitudinal e- profile (true, simchannel);t;E (MeV)",
169  LBINS,
170  LMIN,
171  LMAX);
172  fShowerProfileHitLong = tfs->make<TProfile>(
173  "fShowerProfileHitLong", "longitudinal e- profile (true, hit);t;E (MeV)", LBINS, LMIN, LMAX);
174  fShowerProfileRecoLong = tfs->make<TProfile>(
175  "fShowerProfileRecoLong", "longitudinal e- profile (reco);t;Q", LBINS, LMIN, LMAX);
176 
177  fShowerProfileSimLong2D = tfs->make<TProfile2D>(
178  "fShowerProfileSimLong2D",
179  "longitudinal e- profile (true, simchannel);t;electron energy (MeV);E (MeV)",
180  LBINS,
181  LMIN,
182  LMAX,
183  EBINS,
184  EMIN,
185  EMAX);
186  fShowerProfileHitLong2D =
187  tfs->make<TProfile2D>("fShowerProfileHitLong2D",
188  "longitudinal e- profile (true, hit);t;electron energy (MeV);E (MeV)",
189  LBINS,
190  LMIN,
191  LMAX,
192  EBINS,
193  EMIN,
194  EMAX);
195  fShowerProfileRecoLong2D =
196  tfs->make<TProfile2D>("fShowerProfileRecoLong2D",
197  "longitudinal e- profile (reco);t;electron energy (MeV);Q",
198  LBINS,
199  LMIN,
200  LMAX,
201  EBINS,
202  EMIN,
203  EMAX);
204 
205  fShowerProfileSimTrans =
206  tfs->make<TProfile>("fShowerProfileSimTrans",
207  "transverse e- profile (true, simchannel);dist (cm);E (MeV)",
208  TBINS,
209  TMIN,
210  TMAX);
211  fShowerProfileHitTrans =
212  tfs->make<TProfile>("fShowerProfileHitTrans",
213  "transverse e- profile (true, hit);dist (cm);E (MeV)",
214  TBINS,
215  TMIN,
216  TMAX);
217  fShowerProfileRecoTrans = tfs->make<TProfile>(
218  "fShowerProfileRecoTrans", "transverse e- profile (reco);dist (cm);Q", TBINS, TMIN, TMAX);
219 
220  fShowerProfileSimTrans2D = tfs->make<TProfile2D>(
221  "fShowerProfileSimTrans2D",
222  "transverse e- profile (true, simchannel);t;electron energy (MeV);E (MeV)",
223  TBINS,
224  TMIN,
225  TMAX,
226  EBINS,
227  EMIN,
228  EMAX);
229  fShowerProfileHitTrans2D =
230  tfs->make<TProfile2D>("fShowerProfileHitTrans2D",
231  "transverse e- profile (true, hit);t;electron energy (MeV);E (MeV)",
232  TBINS,
233  TMIN,
234  TMAX,
235  EBINS,
236  EMIN,
237  EMAX);
238  fShowerProfileRecoTrans2D =
239  tfs->make<TProfile2D>("fShowerProfileRecoTrans2D",
240  "transverse e- profile (reco);t;electron energy (MeV);Q",
241  TBINS,
242  TMIN,
243  TMAX,
244  EBINS,
245  EMIN,
246  EMAX);
247 
248  fShowerProfileSimTrans2D_1 = tfs->make<TProfile2D>(
249  "fShowerProfileSimTrans2D_1",
250  "transverse e- profile [0 <= t < 1] (true, simchannel);t;electron energy (MeV);E (MeV)",
251  TBINS,
252  TMIN,
253  TMAX,
254  EBINS,
255  EMIN,
256  EMAX);
257  fShowerProfileHitTrans2D_1 = tfs->make<TProfile2D>(
258  "fShowerProfileHitTrans2D_1",
259  "transverse e- profile [0 <= t < 1] (true, hit);t;electron energy (MeV);E (MeV)",
260  TBINS,
261  TMIN,
262  TMAX,
263  EBINS,
264  EMIN,
265  EMAX);
266  fShowerProfileRecoTrans2D_1 =
267  tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_1",
268  "transverse e- profile [0 <= t < 1] (reco);t;electron energy (MeV);Q",
269  TBINS,
270  TMIN,
271  TMAX,
272  EBINS,
273  EMIN,
274  EMAX);
275 
276  fShowerProfileSimTrans2D_2 = tfs->make<TProfile2D>(
277  "fShowerProfileSimTrans2D_2",
278  "transverse e- profile [1 <= t < 2] (true, simchannel);t;electron energy (MeV);E (MeV)",
279  TBINS,
280  TMIN,
281  TMAX,
282  EBINS,
283  EMIN,
284  EMAX);
285  fShowerProfileHitTrans2D_2 = tfs->make<TProfile2D>(
286  "fShowerProfileHitTrans2D_2",
287  "transverse e- profile [1 <= t < 2] (true, hit);t;electron energy (MeV);E (MeV)",
288  TBINS,
289  TMIN,
290  TMAX,
291  EBINS,
292  EMIN,
293  EMAX);
294  fShowerProfileRecoTrans2D_2 =
295  tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_2",
296  "transverse e- profile [1 <= t < 2] (reco);t;electron energy (MeV);Q",
297  TBINS,
298  TMIN,
299  TMAX,
300  EBINS,
301  EMIN,
302  EMAX);
303 
304  fShowerProfileSimTrans2D_3 = tfs->make<TProfile2D>(
305  "fShowerProfileSimTrans2D_3",
306  "transverse e- profile [2 <= t < 3] (true, simchannel);t;electron energy (MeV);E (MeV)",
307  TBINS,
308  TMIN,
309  TMAX,
310  EBINS,
311  EMIN,
312  EMAX);
313  fShowerProfileHitTrans2D_3 = tfs->make<TProfile2D>(
314  "fShowerProfileHitTrans2D_3",
315  "transverse e- profile [2 <= t < 3] (true, hit);t;electron energy (MeV);E (MeV)",
316  TBINS,
317  TMIN,
318  TMAX,
319  EBINS,
320  EMIN,
321  EMAX);
322  fShowerProfileRecoTrans2D_3 =
323  tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_3",
324  "transverse e- profile [2 <= t < 3] (reco);t;electron energy (MeV);Q",
325  TBINS,
326  TMIN,
327  TMAX,
328  EBINS,
329  EMIN,
330  EMAX);
331 
332  fShowerProfileSimTrans2D_4 = tfs->make<TProfile2D>(
333  "fShowerProfileSimTrans2D_4",
334  "transverse e- profile [3 <= t < 4] (true, simchannel);t;electron energy (MeV);E (MeV)",
335  TBINS,
336  TMIN,
337  TMAX,
338  EBINS,
339  EMIN,
340  EMAX);
341  fShowerProfileHitTrans2D_4 = tfs->make<TProfile2D>(
342  "fShowerProfileHitTrans2D_4",
343  "transverse e- profile [3 <= t < 4] (true, hit);t;electron energy (MeV);E (MeV)",
344  TBINS,
345  TMIN,
346  TMAX,
347  EBINS,
348  EMIN,
349  EMAX);
350  fShowerProfileRecoTrans2D_4 =
351  tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_4",
352  "transverse e- profile [3 <= t < 4] (reco);t;electron energy (MeV);Q",
353  TBINS,
354  TMIN,
355  TMAX,
356  EBINS,
357  EMIN,
358  EMAX);
359 
360  fShowerProfileSimTrans2D_5 = tfs->make<TProfile2D>(
361  "fShowerProfileSimTrans2D_5",
362  "transverse e- profile [4 <= t < 5] (true, simchannel);t;electron energy (MeV);E (MeV)",
363  TBINS,
364  TMIN,
365  TMAX,
366  EBINS,
367  EMIN,
368  EMAX);
369  fShowerProfileHitTrans2D_5 = tfs->make<TProfile2D>(
370  "fShowerProfileHitTrans2D_5",
371  "transverse e- profile [4 <= t < 5] (true, hit);t;electron energy (MeV);E (MeV)",
372  TBINS,
373  TMIN,
374  TMAX,
375  EBINS,
376  EMIN,
377  EMAX);
378  fShowerProfileRecoTrans2D_5 =
379  tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_5",
380  "transverse e- profile [4 <= t < 5] (reco);t;electron energy (MeV);Q",
381  TBINS,
382  TMIN,
383  TMAX,
384  EBINS,
385  EMIN,
386  EMAX);
387 
388  fLongitudinal = tfs->make<TH3F>("fLongitudinal",
389  "longitudinal e- profile;t;electron energy (MeV);Q",
390  LBINS,
391  LMIN,
392  LMAX,
393  EBINS,
394  EMIN,
395  EMAX,
396  100,
397  0,
398  150000);
399  fTransverse = tfs->make<TH3F>("fTransverse",
400  "transverse e- profile;dist (cm);electron energy (MeV);Q",
401  TBINS,
402  TMIN,
403  TMAX,
404  EBINS,
405  EMIN,
406  EMAX,
407  100,
408  0,
409  150000);
410  fTransverse_1 =
411  tfs->make<TH3F>("fTransverse_1",
412  "transverse e- profile [0 <= t < 1];dist (cm);electron energy (MeV);Q",
413  TBINS,
414  TMIN,
415  TMAX,
416  EBINS,
417  EMIN,
418  EMAX,
419  100,
420  0,
421  100000);
422  fTransverse_2 =
423  tfs->make<TH3F>("fTransverse_2",
424  "transverse e- profile [1 <= t < 2];dist (cm);electron energy (MeV);Q",
425  TBINS,
426  TMIN,
427  TMAX,
428  EBINS,
429  EMIN,
430  EMAX,
431  100,
432  0,
433  100000);
434  fTransverse_3 =
435  tfs->make<TH3F>("fTransverse_3",
436  "transverse e- profile [2 <= t < 3];dist (cm);electron energy (MeV);Q",
437  TBINS,
438  TMIN,
439  TMAX,
440  EBINS,
441  EMIN,
442  EMAX,
443  100,
444  0,
445  100000);
446  fTransverse_4 =
447  tfs->make<TH3F>("fTransverse_4",
448  "transverse e- profile [3 <= t < 4];dist (cm);electron energy (MeV);Q",
449  TBINS,
450  TMIN,
451  TMAX,
452  EBINS,
453  EMIN,
454  EMAX,
455  100,
456  0,
457  100000);
458  fTransverse_5 =
459  tfs->make<TH3F>("fTransverse_5",
460  "transverse e- profile [4 <= t < 5];dist (cm);electron energy (MeV);Q",
461  TBINS,
462  TMIN,
463  TMAX,
464  EBINS,
465  EMIN,
466  EMAX,
467  100,
468  0,
469  100000);
470 
471  // electrons
472  fLongitudinal_electron = tfs->make<TProfile>(
473  "fLongitudinal_electron", "longitudinal e- profile (reco);t;Q", LBINS, LMIN, LMAX);
474  fTransverse1_electron =
475  tfs->make<TProfile>("fTransverse1_electron",
476  "transverse e- profile [0 <= t < 1] (reco);dist (cm);Q",
477  TBINS,
478  TMIN,
479  TMAX);
480  fTransverse2_electron =
481  tfs->make<TProfile>("fTransverse2_electron",
482  "transverse e- profile [1 <= t < 2] (reco);dist (cm);Q",
483  TBINS,
484  TMIN,
485  TMAX);
486  fTransverse3_electron =
487  tfs->make<TProfile>("fTransverse3_electron",
488  "transverse e- profile [2 <= t < 3] (reco);dist (cm);Q",
489  TBINS,
490  TMIN,
491  TMAX);
492  fTransverse4_electron =
493  tfs->make<TProfile>("fTransverse4_electron",
494  "transverse e- profile [3 <= t < 4] (reco);dist (cm);Q",
495  TBINS,
496  TMIN,
497  TMAX);
498  fTransverse5_electron =
499  tfs->make<TProfile>("fTransverse5_electron",
500  "transverse e- profile [4 <= t < 5] (reco);dist (cm);Q",
501  TBINS,
502  TMIN,
503  TMAX);
504 
505  // photons
506  fLongitudinal_photon = tfs->make<TProfile>(
507  "fLongitudinal_photon", "longitudinal photon profile (reco);t;Q", LBINS, LMIN, LMAX);
508  fTransverse1_photon =
509  tfs->make<TProfile>("fTransverse1_photon",
510  "transverse photon profile [0 <= t < 1] (reco);dist (cm);Q",
511  TBINS,
512  TMIN,
513  TMAX);
514  fTransverse2_photon =
515  tfs->make<TProfile>("fTransverse2_photon",
516  "transverse photon profile [1 <= t < 2] (reco);dist (cm);Q",
517  TBINS,
518  TMIN,
519  TMAX);
520  fTransverse3_photon =
521  tfs->make<TProfile>("fTransverse3_photon",
522  "transverse photon profile [2 <= t < 3] (reco);dist (cm);Q",
523  TBINS,
524  TMIN,
525  TMAX);
526  fTransverse4_photon =
527  tfs->make<TProfile>("fTransverse4_photon",
528  "transverse photon profile [3 <= t < 4] (reco);dist (cm);Q",
529  TBINS,
530  TMIN,
531  TMAX);
532  fTransverse5_photon =
533  tfs->make<TProfile>("fTransverse5_photon",
534  "transverse photon profile [4 <= t < 5] (reco);dist (cm);Q",
535  TBINS,
536  TMIN,
537  TMAX);
538 
539  // other
540  fLongitudinal_other = tfs->make<TProfile>(
541  "fLongitudinal_other", "longitudinal other profile (reco);t;Q", LBINS, LMIN, LMAX);
542  fTransverse1_other =
543  tfs->make<TProfile>("fTransverse1_other",
544  "transverse other profile [0 <= t < 1] (reco);dist (cm);Q",
545  TBINS,
546  TMIN,
547  TMAX);
548  fTransverse2_other =
549  tfs->make<TProfile>("fTransverse2_other",
550  "transverse other profile [1 <= t < 2] (reco);dist (cm);Q",
551  TBINS,
552  TMIN,
553  TMAX);
554  fTransverse3_other =
555  tfs->make<TProfile>("fTransverse3_other",
556  "transverse other profile [2 <= t < 3] (reco);dist (cm);Q",
557  TBINS,
558  TMIN,
559  TMAX);
560  fTransverse4_other =
561  tfs->make<TProfile>("fTransverse4_other",
562  "transverse other profile [3 <= t < 4] (reco);dist (cm);Q",
563  TBINS,
564  TMIN,
565  TMAX);
566  fTransverse5_other =
567  tfs->make<TProfile>("fTransverse5_other",
568  "transverse other profile [4 <= t < 5] (reco);dist (cm);Q",
569  TBINS,
570  TMIN,
571  TMAX);
572 
573 } // beginJob
template<typename T , BranchType BT>
ProductToken< T > art::ModuleBase::consumes ( InputTag const &  tag)
protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

62  {
63  return collector_.consumes<T, BT>(tag);
64  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::consumes
ProductToken< T > consumes(InputTag const &)
Definition: ConsumesCollector.h:50
ConsumesCollector & art::ModuleBase::consumesCollector ( )
protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

58  {
59  return collector_;
60  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
template<typename T , BranchType BT>
void art::ModuleBase::consumesMany ( )
protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

76  {
77  collector_.consumesMany<T, BT>();
78  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::ModuleBase::consumesView ( InputTag const &  )
protectedinherited
template<typename T , BranchType BT>
ViewToken<T> art::ModuleBase::consumesView ( InputTag const &  tag)
inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

69  {
70  return collector_.consumesView<T, BT>(tag);
71  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::consumesView
ViewToken< Element > consumesView(InputTag const &)
void art::detail::Analyzer::doBeginJob ( SharedResources const &  resources)
inherited

Definition at line 25 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

26  {
27  setupQueues(resources);
28  ProcessingFrame const frame{ScheduleID{}};
29  beginJobWithFrame(frame);
30  }
art::detail::Analyzer::beginJobWithFrame
virtual void beginJobWithFrame(ProcessingFrame const &)=0
art::detail::Analyzer::setupQueues
virtual void setupQueues(SharedResources const &)=0
bool art::detail::Analyzer::doBeginRun ( RunPrincipal &  rp,
ModuleContext const &  mc 
)
inherited

Definition at line 68 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

69  {
70  ProcessingFrame const frame{mc.scheduleID()};
71  beginRunWithFrame(std::as_const(rp).makeRun(mc), frame);
72  return true;
73  }
art::detail::Analyzer::beginRunWithFrame
virtual void beginRunWithFrame(Run const &, ProcessingFrame const &)=0
bool art::detail::Analyzer::doBeginSubRun ( SubRunPrincipal &  srp,
ModuleContext const &  mc 
)
inherited

Definition at line 84 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

85  {
86  ProcessingFrame const frame{mc.scheduleID()};
87  beginSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
88  return true;
89  }
art::detail::Analyzer::beginSubRunWithFrame
virtual void beginSubRunWithFrame(SubRun const &, ProcessingFrame const &)=0
void art::detail::Analyzer::doEndJob ( )
inherited

Definition at line 33 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

34  {
35  ProcessingFrame const frame{ScheduleID{}};
36  endJobWithFrame(frame);
37  }
art::detail::Analyzer::endJobWithFrame
virtual void endJobWithFrame(ProcessingFrame const &)=0
bool art::detail::Analyzer::doEndRun ( RunPrincipal &  rp,
ModuleContext const &  mc 
)
inherited

Definition at line 76 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

77  {
78  ProcessingFrame const frame{mc.scheduleID()};
79  endRunWithFrame(std::as_const(rp).makeRun(mc), frame);
80  return true;
81  }
art::detail::Analyzer::endRunWithFrame
virtual void endRunWithFrame(Run const &, ProcessingFrame const &)=0
bool art::detail::Analyzer::doEndSubRun ( SubRunPrincipal &  srp,
ModuleContext const &  mc 
)
inherited

Definition at line 92 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

93  {
94  ProcessingFrame const frame{mc.scheduleID()};
95  endSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
96  return true;
97  }
art::detail::Analyzer::endSubRunWithFrame
virtual void endSubRunWithFrame(SubRun const &, ProcessingFrame const &)=0
bool art::detail::Analyzer::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 
)
inherited

Definition at line 100 of file Analyzer.cc.

References e, and art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

105  {
106  auto const e = std::as_const(ep).makeEvent(mc);
107  if (wantEvent(mc.scheduleID(), e)) {
108  ++counts_run;
109  ProcessingFrame const frame{mc.scheduleID()};
110  analyzeWithFrame(e, frame);
111  ++counts_passed;
112  }
113  return true;
114  }
art::Observer::wantEvent
bool wantEvent(ScheduleID id, Event const &e) const
Definition: Observer.cc:63
e
Float_t e
Definition: plot.C:35
art::detail::Analyzer::analyzeWithFrame
virtual void analyzeWithFrame(Event const &, ProcessingFrame const &)=0
void art::detail::Analyzer::doRespondToCloseInputFile ( FileBlock const &  fb)
inherited

Definition at line 47 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

48  {
49  ProcessingFrame const frame{ScheduleID{}};
50  respondToCloseInputFileWithFrame(fb, frame);
51  }
fb
TFile fb("Li6.root")
art::detail::Analyzer::respondToCloseInputFileWithFrame
virtual void respondToCloseInputFileWithFrame(FileBlock const &, ProcessingFrame const &)=0
void art::detail::Analyzer::doRespondToCloseOutputFiles ( FileBlock const &  fb)
inherited

Definition at line 61 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

62  {
63  ProcessingFrame const frame{ScheduleID{}};
64  respondToCloseOutputFilesWithFrame(fb, frame);
65  }
art::detail::Analyzer::respondToCloseOutputFilesWithFrame
virtual void respondToCloseOutputFilesWithFrame(FileBlock const &, ProcessingFrame const &)=0
fb
TFile fb("Li6.root")
void art::detail::Analyzer::doRespondToOpenInputFile ( FileBlock const &  fb)
inherited

Definition at line 40 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

41  {
42  ProcessingFrame const frame{ScheduleID{}};
43  respondToOpenInputFileWithFrame(fb, frame);
44  }
fb
TFile fb("Li6.root")
art::detail::Analyzer::respondToOpenInputFileWithFrame
virtual void respondToOpenInputFileWithFrame(FileBlock const &, ProcessingFrame const &)=0
void art::detail::Analyzer::doRespondToOpenOutputFiles ( FileBlock const &  fb)
inherited

Definition at line 54 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

55  {
56  ProcessingFrame const frame{ScheduleID{}};
57  respondToOpenOutputFilesWithFrame(fb, frame);
58  }
art::detail::Analyzer::respondToOpenOutputFilesWithFrame
virtual void respondToOpenOutputFilesWithFrame(FileBlock const &, ProcessingFrame const &)=0
fb
TFile fb("Li6.root")
std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const
inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

44  {
45  return collector_.getConsumables();
46  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::getConsumables
std::array< std::vector< ProductInfo >, NumBranchTypes > const & getConsumables() const
Definition: ConsumesCollector.cc:13
Handle< TriggerResults > art::Observer::getTriggerResults ( Event const &  e) const
protectedinherited

Definition at line 75 of file Observer.cc.

References art::ProductRetriever::get(), and art::Observer::selectors_.

Referenced by art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

76  {
77  if (selectors_) {
78  return selectors_->getOneTriggerResults(e);
79  }
80 
81  // The following applies for cases where no SelectEvents entries
82  // exist.
83  Handle<TriggerResults> h;
84  if (e.get(empty_process_name, h)) {
85  return h;
86  }
87  return Handle<TriggerResults>{};
88  }
e
Float_t e
Definition: plot.C:35
art::Observer::selectors_
std::optional< detail::ProcessAndEventSelectors > selectors_
Definition: Observer.h:79
std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const &  wp)
inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

38  {
39  return doMakeWorker(wp);
40  }
art::ModuleBase::doMakeWorker
virtual std::unique_ptr< Worker > doMakeWorker(WorkerParams const &wp)=0
template<typename T , BranchType BT>
ProductToken< T > art::ModuleBase::mayConsume ( InputTag const &  tag)
protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

83  {
84  return collector_.mayConsume<T, BT>(tag);
85  }
art::ConsumesCollector::mayConsume
ProductToken< T > mayConsume(InputTag const &)
Definition: ConsumesCollector.h:75
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
template<typename T , BranchType BT>
void art::ModuleBase::mayConsumeMany ( )
protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

97  {
98  collector_.mayConsumeMany<T, BT>();
99  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const &  )
protectedinherited
template<typename T , BranchType BT>
ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const &  tag)
inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

90  {
91  return collector_.mayConsumeView<T, BT>(tag);
92  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::mayConsumeView
ViewToken< Element > mayConsumeView(InputTag const &)
ModuleDescription const & art::ModuleBase::moduleDescription ( ) const
inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

14  {
15  if (md_.has_value()) {
16  return *md_;
17  }
18 
19  throw Exception{errors::LogicError,
20  "There was an error while calling moduleDescription().\n"}
21  << "The moduleDescription() base-class member function cannot be called\n"
22  "during module construction. To determine which module is "
23  "responsible\n"
24  "for calling it, find the '<module type>:<module "
25  "label>@Construction'\n"
26  "tag in the message prefix above. Please contact artists@fnal.gov\n"
27  "for guidance.\n";
28  }
art::Exception
cet::coded_exception< errors::ErrorCodes, ExceptionDetail::translate > Exception
Definition: Exception.h:66
art::ModuleBase::md_
std::optional< ModuleDescription > md_
Definition: ModuleBase.h:55
string const & art::Observer::processName ( ) const
protectedinherited

Definition at line 57 of file Observer.cc.

References art::Observer::process_name_.

Referenced by art::FileDumperOutput::printPrincipal().

58  {
59  return process_name_;
60  }
art::Observer::process_name_
std::string process_name_
Definition: Observer.h:76
void art::ModuleBase::setModuleDescription ( ModuleDescription const &  md)
inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

32  {
33  md_ = md;
34  }
art::ModuleBase::md_
std::optional< ModuleDescription > md_
Definition: ModuleBase.h:55
void shower::TCShowerTemplateMaker::showerProfile ( detinfo::DetectorClocksData const &  clockData,
detinfo::DetectorPropertiesData const &  detProp,
std::vector< art::Ptr< recob::Hit >>  showerhits,
TVector3  shwvtx,
TVector3  shwdir,
double  elep 
)
private

Definition at line 629 of file TCShowerTemplateMaker_module.cc.

References util::abs(), detinfo::DetectorPropertiesData::ConvertXToTicks(), detinfo::DetectorPropertiesData::DriftVelocity(), detinfo::DetectorPropertiesData::Efield(), fCalorimetryAlg, fLongitudinal, fShowerProfileRecoLong, fShowerProfileRecoLong2D, fShowerProfileRecoTrans, fShowerProfileRecoTrans2D, fShowerProfileRecoTrans2D_1, fShowerProfileRecoTrans2D_2, fShowerProfileRecoTrans2D_3, fShowerProfileRecoTrans2D_4, fShowerProfileRecoTrans2D_5, fTransverse, fTransverse_1, fTransverse_2, fTransverse_3, fTransverse_4, fTransverse_5, LBINS, calo::CalorimetryAlg::LifetimeCorrection(), LMAX, LMIN, detinfo::sampling_rate(), TBINS, detinfo::DetectorPropertiesData::Temperature(), TMAX, TMIN, geo::vect::toPoint(), geo::GeometryCore::WireCoordinate(), geo::GeometryCore::WirePitch(), and X0.

Referenced by analyze().

635 {
636  art::ServiceHandle<geo::Geometry const> geom;
637 
638  auto collectionPlane = geo::PlaneID(0, 0, 1);
639 
640  double shwVtxTime = detProp.ConvertXToTicks(shwvtx[0], collectionPlane);
641  using geo::vect::toPoint;
642  double shwVtxWire = geom->WireCoordinate(toPoint(shwvtx), collectionPlane);
643 
644  double shwTwoTime = detProp.ConvertXToTicks(shwvtx[0] + shwdir[0], collectionPlane);
645  double shwTwoWire = geom->WireCoordinate(toPoint(shwvtx + shwdir), collectionPlane);
646 
647  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
648  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
649 
650  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
651  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
652  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
653  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
654  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
655 
656  for (size_t i = 0; i < showerhits.size(); ++i) {
657  if (showerhits[i]->WireID().Plane != collectionPlane.Plane) continue;
658 
659  double wirePitch = geom->WirePitch(showerhits[i]->WireID());
660  double tickToDist = detProp.DriftVelocity(detProp.Efield(), detProp.Temperature());
661  tickToDist *= 1.e-3 * sampling_rate(clockData); // 1e-3 is conversion of 1/us to 1/ns
662 
663  double xvtx = shwVtxTime * tickToDist;
664  double yvtx = shwVtxWire * wirePitch;
665 
666  double xtwo = shwTwoTime * tickToDist;
667  double ytwo = shwTwoWire * wirePitch;
668 
669  double xtwoorth = (ytwo - yvtx) + xvtx;
670  double ytwoorth = -(xtwo - xvtx) + yvtx;
671 
672  double xhit = showerhits[i]->PeakTime() * tickToDist;
673  double yhit = showerhits[i]->WireID().Wire * wirePitch;
674 
675  double ldist = std::abs((ytwoorth - yvtx) * xhit - (xtwoorth - xvtx) * yhit + xtwoorth * yvtx -
676  ytwoorth * xvtx) /
677  std::hypot(ytwoorth - yvtx, xtwoorth - xvtx);
678  double tdist = ((ytwo - yvtx) * xhit - (xtwo - xvtx) * yhit + xtwo * yvtx - ytwo * xvtx) /
679  std::hypot(ytwo - yvtx, xtwo - xvtx);
680 
681  double to3D = 1. / std::hypot(xvtx - xtwo,
682  yvtx - ytwo); // distance between two points in 3D space is one
683  ldist *= to3D;
684  tdist *= to3D;
685  double t = ldist / X0;
686 
687  double Q = showerhits[i]->Integral() *
688  fCalorimetryAlg.LifetimeCorrection(clockData, detProp, showerhits[i]->PeakTime());
689 
690  ltemp->Fill(t, Q);
691  ttemp->Fill(tdist, Q);
692 
693  if (t < 1)
694  ttemp_1->Fill(tdist, Q);
695  else if (t < 2)
696  ttemp_2->Fill(tdist, Q);
697  else if (t < 3)
698  ttemp_3->Fill(tdist, Q);
699  else if (t < 4)
700  ttemp_4->Fill(tdist, Q);
701  else if (t < 5)
702  ttemp_5->Fill(tdist, Q);
703 
704  } // loop through showerhits
705 
706  for (int i = 0; i < LBINS; ++i) {
707  if (ltemp->GetBinContent(i + 1) == 0) continue;
708  fShowerProfileRecoLong->Fill(ltemp->GetBinCenter(i + 1), ltemp->GetBinContent(i + 1));
709  fShowerProfileRecoLong2D->Fill(ltemp->GetBinCenter(i + 1), elep, ltemp->GetBinContent(i + 1));
710  fLongitudinal->Fill(ltemp->GetBinCenter(i + 1), elep, ltemp->GetBinContent(i + 1));
711  }
712 
713  for (int i = 0; i < TBINS; ++i) {
714  if (ttemp->GetBinContent(i + 1) == 0) continue;
715  fShowerProfileRecoTrans->Fill(ttemp->GetBinCenter(i + 1), ttemp->GetBinContent(i + 1));
716  fShowerProfileRecoTrans2D->Fill(ttemp->GetBinCenter(i + 1), elep, ttemp->GetBinContent(i + 1));
717  fTransverse->Fill(ttemp->GetBinCenter(i + 1), elep, ttemp->GetBinContent(i + 1));
718 
719  fShowerProfileRecoTrans2D_1->Fill(
720  ttemp_1->GetBinCenter(i + 1), elep, ttemp_1->GetBinContent(i + 1));
721  fShowerProfileRecoTrans2D_2->Fill(
722  ttemp_2->GetBinCenter(i + 1), elep, ttemp_2->GetBinContent(i + 1));
723  fShowerProfileRecoTrans2D_3->Fill(
724  ttemp_3->GetBinCenter(i + 1), elep, ttemp_3->GetBinContent(i + 1));
725  fShowerProfileRecoTrans2D_4->Fill(
726  ttemp_4->GetBinCenter(i + 1), elep, ttemp_4->GetBinContent(i + 1));
727  fShowerProfileRecoTrans2D_5->Fill(
728  ttemp_5->GetBinCenter(i + 1), elep, ttemp_5->GetBinContent(i + 1));
729 
730  fTransverse_1->Fill(ttemp_1->GetBinCenter(i + 1), elep, ttemp_1->GetBinContent(i + 1));
731  fTransverse_2->Fill(ttemp_2->GetBinCenter(i + 1), elep, ttemp_2->GetBinContent(i + 1));
732  fTransverse_3->Fill(ttemp_3->GetBinCenter(i + 1), elep, ttemp_3->GetBinContent(i + 1));
733  fTransverse_4->Fill(ttemp_4->GetBinCenter(i + 1), elep, ttemp_4->GetBinContent(i + 1));
734  fTransverse_5->Fill(ttemp_5->GetBinCenter(i + 1), elep, ttemp_5->GetBinContent(i + 1));
735  }
736 } // showerProfile
geo::GeometryCore::WireCoordinate
Length_t WireCoordinate(Point_t const &pos, PlaneID const &planeid) const
Returns the index of the nearest wire to the specified position.
Definition: GeometryCore.cxx:1093
geo::vect::toPoint
::geo::Point_t toPoint(Point const &p)
Convert the specified point into a geo::Point_t.
Definition: geo_vectors_utils.h:1584
util::abs
constexpr auto abs(T v)
Returns the absolute value of the argument.
Definition: constexpr_math.h:40
geo::fhicl::PlaneID
IDparameter< geo::PlaneID > PlaneID
Member type of validated geo::PlaneID parameter.
Definition: geo_types_fhicl.h:355
geo::fhicl::WireID
IDparameter< geo::WireID > WireID
Member type of validated geo::WireID parameter.
Definition: geo_types_fhicl.h:383
calo::CalorimetryAlg::LifetimeCorrection
double LifetimeCorrection(detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, double time, double T0=0) const
Definition: CalorimetryAlg.cxx:221
detinfo::sampling_rate
double sampling_rate(DetectorClocksData const &data)
Returns the period of the TPC readout electronics clock.
Definition: DetectorClocksData.h:575
geo::GeometryCore::WirePitch
Length_t WirePitch(PlaneID const &planeid=plane_zero) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:763
void shower::TCShowerTemplateMaker::showerProfileTrue ( detinfo::DetectorClocksData const &  clockData,
detinfo::DetectorPropertiesData const &  detProp,
std::vector< art::Ptr< recob::Hit >>  allhits,
double  elep 
)
private

Definition at line 740 of file TCShowerTemplateMaker_module.cc.

References util::abs(), detinfo::DetectorPropertiesData::ConvertXToTicks(), detinfo::DetectorPropertiesData::DriftVelocity(), detinfo::DetectorPropertiesData::Efield(), energy, fShowerProfileHitLong, fShowerProfileHitLong2D, fShowerProfileHitTrans, fShowerProfileHitTrans2D, fShowerProfileHitTrans2D_1, fShowerProfileHitTrans2D_2, fShowerProfileHitTrans2D_3, fShowerProfileHitTrans2D_4, fShowerProfileHitTrans2D_5, cheat::BackTrackerService::HitToEveTrackIDEs(), LBINS, LMAX, LMIN, detinfo::sampling_rate(), TBINS, detinfo::DetectorPropertiesData::Temperature(), TMAX, TMIN, cheat::ParticleInventoryService::TrackIdToParticle_P(), geo::GeometryCore::WireCoordinate(), geo::GeometryCore::WirePitch(), and X0.

Referenced by analyze().

745 {
746  art::ServiceHandle<geo::Geometry const> geom;
747  auto collectionPlane = geo::PlaneID(0, 0, 1);
748  art::ServiceHandle<cheat::BackTrackerService const> btserv;
749  art::ServiceHandle<cheat::ParticleInventoryService const> piserv;
750  std::map<int, double> trkID_E;
751 
752  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
753  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
754 
755  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
756  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
757  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
758  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
759  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
760 
761  double xvtx = -999;
762  double yvtx = -999;
763  double zvtx = -999;
764  double xtwo = -999;
765  double ytwo = -999;
766  double ztwo = -999;
767  double shwvtxT = -999;
768  double shwvtxW = -999;
769  double shwtwoT = -999;
770  double shwtwoW = -999;
771 
772  double shwvtxx = -999;
773  double shwvtxy = -999;
774  double shwtwox = -999;
775  double shwtwoy = -999;
776  double xtwoorth = -999;
777  double ytwoorth = -999;
778 
779  double wirePitch = -999;
780  double tickToDist = -999;
781 
782  bool foundParent = false;
783 
784  for (auto const& hit : allhits) {
785  if (hit->WireID().Plane != collectionPlane.Plane) continue;
786 
787  // art::Ptr<recob::Hit> hit = allhits[i];
788  std::vector<sim::TrackIDE> trackIDs = btserv->HitToEveTrackIDEs(clockData, hit);
789 
790  for (size_t j = 0; j < trackIDs.size(); ++j) {
791  // only want energy associated with the electron and electron must have neutrino mother
792  if (std::abs((piserv->TrackIdToParticle_P(trackIDs[j].trackID))->PdgCode()) != 11) continue;
793  if (std::abs((piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Mother()) != 0) continue;
794 
795  if (!foundParent) {
796  xvtx = (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Vx();
797  yvtx = (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Vy();
798  zvtx = (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Vz();
799 
800  xtwo = xvtx + (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Px();
801  ytwo = yvtx + (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Py();
802  ztwo = zvtx + (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Pz();
803 
804  shwvtxT = detProp.ConvertXToTicks(xvtx, collectionPlane);
805  geo::Point_t const vtx{xvtx, yvtx, zvtx};
806  shwvtxW = geom->WireCoordinate(vtx, collectionPlane);
807 
808  shwtwoT = detProp.ConvertXToTicks(xtwo, collectionPlane);
809  shwtwoW = geom->WireCoordinate(vtx, collectionPlane);
810 
811  wirePitch = geom->WirePitch(hit->WireID());
812  tickToDist = detProp.DriftVelocity(detProp.Efield(), detProp.Temperature());
813  tickToDist *= 1.e-3 * sampling_rate(clockData); // 1e-3 is conversion of 1/us to 1/ns
814 
815  shwvtxx = shwvtxT * tickToDist;
816  shwvtxy = shwvtxW * wirePitch;
817 
818  shwtwox = shwtwoT * tickToDist;
819  shwtwoy = shwtwoW * wirePitch;
820 
821  xtwoorth = (shwtwoy - shwvtxy) + shwvtxx;
822  ytwoorth = -(shwtwox - shwvtxx) + shwvtxy;
823 
824  foundParent = true;
825  }
826  double xhit = hit->PeakTime() * tickToDist;
827  double yhit = hit->WireID().Wire * wirePitch;
828 
829  double ldist = abs((ytwoorth - shwvtxy) * xhit - (xtwoorth - shwvtxx) * yhit +
830  xtwoorth * shwvtxy - ytwoorth * shwvtxx) /
831  std::hypot(ytwoorth - shwvtxy, xtwoorth - shwvtxx);
832  double tdist = ((shwtwoy - shwvtxy) * xhit - (shwtwox - shwvtxx) * yhit + shwtwox * shwvtxy -
833  shwtwoy * shwvtxx) /
834  std::hypot(shwtwoy - shwvtxy, shwtwox - shwvtxx);
835 
836  double to3D = std::hypot(xvtx - xtwo, yvtx - ytwo, zvtx - ztwo) /
837  std::hypot(shwvtxx - shwtwox,
838  shwvtxy - shwtwoy); // distance between two points in 3D space is one
839  ldist *= to3D;
840  tdist *= to3D;
841  double t = ldist / X0;
842 
843  double energy = trackIDs[j].energy;
844 
845  ltemp->Fill(t, energy);
846  ttemp->Fill(tdist, energy);
847 
848  if (t < 1)
849  ttemp_1->Fill(tdist, energy);
850  else if (t < 2)
851  ttemp_2->Fill(tdist, energy);
852  else if (t < 3)
853  ttemp_3->Fill(tdist, energy);
854  else if (t < 4)
855  ttemp_4->Fill(tdist, energy);
856  else if (t < 5)
857  ttemp_5->Fill(tdist, energy);
858 
859  } // loop through track IDE
860 
861  } // loop through all hits
862 
863  for (int i = 0; i < LBINS; ++i) {
864  if (ltemp->GetBinContent(i + 1) == 0) continue;
865  fShowerProfileHitLong->Fill(ltemp->GetBinCenter(i + 1), ltemp->GetBinContent(i + 1));
866  fShowerProfileHitLong2D->Fill(ltemp->GetBinCenter(i + 1), elep, ltemp->GetBinContent(i + 1));
867  }
868 
869  for (int i = 0; i < TBINS; ++i) {
870  if (ttemp->GetBinContent(i + 1) == 0) continue;
871  fShowerProfileHitTrans->Fill(ttemp->GetBinCenter(i + 1), ttemp->GetBinContent(i + 1));
872  fShowerProfileHitTrans2D->Fill(ttemp->GetBinCenter(i + 1), elep, ttemp->GetBinContent(i + 1));
873 
874  fShowerProfileHitTrans2D_1->Fill(
875  ttemp_1->GetBinCenter(i + 1), elep, ttemp_1->GetBinContent(i + 1));
876  fShowerProfileHitTrans2D_2->Fill(
877  ttemp_2->GetBinCenter(i + 1), elep, ttemp_2->GetBinContent(i + 1));
878  fShowerProfileHitTrans2D_3->Fill(
879  ttemp_3->GetBinCenter(i + 1), elep, ttemp_3->GetBinContent(i + 1));
880  fShowerProfileHitTrans2D_4->Fill(
881  ttemp_4->GetBinCenter(i + 1), elep, ttemp_4->GetBinContent(i + 1));
882  fShowerProfileHitTrans2D_5->Fill(
883  ttemp_5->GetBinCenter(i + 1), elep, ttemp_5->GetBinContent(i + 1));
884  }
885 } // showerProfileTrue
geo::GeometryCore::WireCoordinate
Length_t WireCoordinate(Point_t const &pos, PlaneID const &planeid) const
Returns the index of the nearest wire to the specified position.
Definition: GeometryCore.cxx:1093
cheat::ParticleInventoryService::TrackIdToParticle_P
const simb::MCParticle * TrackIdToParticle_P(int id) const
Definition: ParticleInventoryService.cc:154
util::abs
constexpr auto abs(T v)
Returns the absolute value of the argument.
Definition: constexpr_math.h:40
geo::fhicl::PlaneID
IDparameter< geo::PlaneID > PlaneID
Member type of validated geo::PlaneID parameter.
Definition: geo_types_fhicl.h:355
energy
double energy
Definition: plottest35.C:25
hit
Detector simulation of raw signals on wires.
Definition: DumpHits_module.cc:25
geo::Point_t
ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Point_t
Type for representation of position in physical 3D space.
Definition: geo_vectors.h:180
cheat::BackTrackerService::HitToEveTrackIDEs
std::vector< sim::TrackIDE > HitToEveTrackIDEs(detinfo::DetectorClocksData const &clockData, recob::Hit const &hit) const
Definition: BackTrackerService.cc:192
detinfo::sampling_rate
double sampling_rate(DetectorClocksData const &data)
Returns the period of the TPC readout electronics clock.
Definition: DetectorClocksData.h:575
geo::GeometryCore::WirePitch
Length_t WirePitch(PlaneID const &planeid=plane_zero) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:763
void shower::TCShowerTemplateMaker::showerProfileTrue ( std::vector< art::Ptr< sim::SimChannel >>  allchan,
simb::MCParticle  electron 
)
private

Definition at line 889 of file TCShowerTemplateMaker_module.cc.

References util::abs(), sim::SimChannel::Channel(), DEFINE_ART_MODULE, simb::MCParticle::E(), edep, energy, fShowerProfileSimLong, fShowerProfileSimLong2D, fShowerProfileSimTrans, fShowerProfileSimTrans2D, fShowerProfileSimTrans2D_1, fShowerProfileSimTrans2D_2, fShowerProfileSimTrans2D_3, fShowerProfileSimTrans2D_4, fShowerProfileSimTrans2D_5, geo::kV, LBINS, LMAX, LMIN, simb::MCParticle::PdgCode(), simb::MCParticle::Px(), simb::MCParticle::Py(), simb::MCParticle::Pz(), sim::MCEnDep::SetEnergy(), sim::MCEnDep::SetTrackId(), sim::MCEnDep::SetVertex(), TBINS, sim::SimChannel::TDCIDEMap(), TMAX, TMIN, cheat::ParticleInventoryService::TrackIdToMotherParticle_P(), geo::GeometryCore::View(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), simb::MCParticle::Vz(), x, X0, x2, y, y2, and z.

892 {
893  art::ServiceHandle<cheat::ParticleInventoryService const> piserv;
894  art::ServiceHandle<geo::Geometry const> geom;
895 
896  std::vector<sim::MCEnDep> alledep;
897 
898  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
899  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
900 
901  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
902  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
903  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
904  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
905  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
906 
907  // get electron energy depositions
908  for (size_t i = 0; i < allchan.size(); ++i) {
909  art::Ptr<sim::SimChannel> simchan = allchan[i];
910  if (geom->View(simchan->Channel()) != geo::kV) continue;
911  auto tdc_ide_map = simchan->TDCIDEMap();
912 
913  for (auto const& tdc_ide_pair : tdc_ide_map) {
914  for (auto const& ide : tdc_ide_pair.second) {
915  if (piserv->TrackIdToMotherParticle_P(ide.trackID) == NULL) continue;
916  if (std::abs(piserv->TrackIdToMotherParticle_P(ide.trackID)->PdgCode()) != 11) continue;
917 
918  sim::MCEnDep edep;
919  edep.SetVertex(ide.x, ide.y, ide.z);
920  edep.SetEnergy(ide.energy);
921  edep.SetTrackId(ide.trackID);
922 
923  alledep.push_back(edep);
924 
925  } // loop through ide
926  } // loop through tdc_ide
927 
928  } // loop through channels
929 
930  double x0 = electron.Vx();
931  double y0 = electron.Vy();
932  double z0 = electron.Vz();
933 
934  double x2 = electron.Px();
935  double y2 = electron.Py();
936  double z2 = electron.Pz();
937 
938  TVector3 v0(x2, y2, z2);
939  v0 = v0.Unit();
940 
941  for (size_t i = 0; i < alledep.size(); ++i) {
942  double x = (double)alledep[i].Vertex()[0];
943  double y = (double)alledep[i].Vertex()[1];
944  double z = (double)alledep[i].Vertex()[2];
945 
946  TVector3 v1(x - x0, y - y0, z - z0);
947 
948  double ldist = v0.Dot(v1);
949  double t = ldist / X0;
950  double tdist = (v0.Orthogonal()).Dot(v1);
951 
952  double energy = alledep[i].Energy();
953 
954  ltemp->Fill(t, energy);
955  ttemp->Fill(tdist, energy);
956 
957  if (t < 1)
958  ttemp_1->Fill(tdist, energy);
959  else if (t < 2)
960  ttemp_2->Fill(tdist, energy);
961  else if (t < 3)
962  ttemp_3->Fill(tdist, energy);
963  else if (t < 4)
964  ttemp_4->Fill(tdist, energy);
965  else if (t < 5)
966  ttemp_5->Fill(tdist, energy);
967  }
968 
969  for (int i = 0; i < LBINS; ++i) {
970  if (ltemp->GetBinContent(i + 1) == 0) continue;
971  fShowerProfileSimLong->Fill(ltemp->GetBinCenter(i + 1), ltemp->GetBinContent(i + 1));
972  fShowerProfileSimLong2D->Fill(
973  ltemp->GetBinCenter(i + 1), electron.E(), ltemp->GetBinContent(i + 1));
974  }
975 
976  for (int i = 0; i < TBINS; ++i) {
977  if (ttemp->GetBinContent(i + 1) == 0) continue;
978  fShowerProfileSimTrans->Fill(ttemp->GetBinCenter(i + 1), ttemp->GetBinContent(i + 1));
979  fShowerProfileSimTrans2D->Fill(
980  ttemp->GetBinCenter(i + 1), electron.E(), ttemp->GetBinContent(i + 1));
981 
982  fShowerProfileSimTrans2D_1->Fill(
983  ttemp_1->GetBinCenter(i + 1), electron.E(), ttemp_1->GetBinContent(i + 1));
984  fShowerProfileSimTrans2D_2->Fill(
985  ttemp_2->GetBinCenter(i + 1), electron.E(), ttemp_2->GetBinContent(i + 1));
986  fShowerProfileSimTrans2D_3->Fill(
987  ttemp_3->GetBinCenter(i + 1), electron.E(), ttemp_3->GetBinContent(i + 1));
988  fShowerProfileSimTrans2D_4->Fill(
989  ttemp_4->GetBinCenter(i + 1), electron.E(), ttemp_4->GetBinContent(i + 1));
990  fShowerProfileSimTrans2D_5->Fill(
991  ttemp_5->GetBinCenter(i + 1), electron.E(), ttemp_5->GetBinContent(i + 1));
992  }
993 } // showerProfileTrue
x
Float_t x
Definition: compare.C:6
simb::MCParticle::E
double E(const int i=0) const
Definition: MCParticle.h:234
simb::MCParticle::PdgCode
int PdgCode() const
Definition: MCParticle.h:213
simb::MCParticle::Py
double Py(const int i=0) const
Definition: MCParticle.h:232
sim::MCEnDep::SetEnergy
void SetEnergy(float e)
Definition: MCDataHolder.h:36
geo::kV
Planes which measure V.
Definition: geo_types.h:136
y
Float_t y
Definition: compare.C:6
z
Double_t z
Definition: plot.C:276
simb::MCParticle::Px
double Px(const int i=0) const
Definition: MCParticle.h:231
util::abs
constexpr auto abs(T v)
Returns the absolute value of the argument.
Definition: constexpr_math.h:40
sim::MCEnDep::SetTrackId
void SetTrackId(unsigned int id)
Definition: MCDataHolder.h:38
y2
Float_t y2[n_points_geant4]
Definition: compare.C:26
energy
double energy
Definition: plottest35.C:25
sim::MCEnDep::SetVertex
void SetVertex(float x, float y, float z)
Definition: MCDataHolder.h:29
edep
Double_t edep
Definition: macro.C:13
sim::SimChannel::Channel
raw::ChannelID_t Channel() const
Returns the readout channel this object describes.
Definition: SimChannel.h:323
simb::MCParticle::Vx
double Vx(const int i=0) const
Definition: MCParticle.h:222
geo::GeometryCore::View
View_t View(PlaneID const &pid) const
Returns the view (wire orientation) on the channels of specified TPC plane.
Definition: GeometryCore.cxx:489
simb::MCParticle::Pz
double Pz(const int i=0) const
Definition: MCParticle.h:233
simb::MCParticle::Vz
double Vz(const int i=0) const
Definition: MCParticle.h:224
sim::SimChannel::TDCIDEMap
TDCIDEs_t const & TDCIDEMap() const
Returns all the deposited energy information as stored.
Definition: SimChannel.h:319
art::detail::Vertex
boost::graph_traits< ModuleGraph >::vertex_descriptor Vertex
Definition: ModuleGraph.h:25
cheat::ParticleInventoryService::TrackIdToMotherParticle_P
const simb::MCParticle * TrackIdToMotherParticle_P(int id) const
Definition: ParticleInventoryService.cc:161
x2
Float_t x2[n_points_geant4]
Definition: compare.C:26
simb::MCParticle::Vy
double Vy(const int i=0) const
Definition: MCParticle.h:223
void art::ModuleBase::sortConsumables ( std::string const &  current_process_name)
inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

50  {
51  // Now that we know we have seen all the consumes declarations,
52  // sort the results for fast lookup later.
53  collector_.sortConsumables(current_process_name);
54  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::sortConsumables
void sortConsumables(std::string const &current_process_name)
Definition: ConsumesCollector.cc:19
bool art::Observer::wantAllEvents ( ) const
inlineprotectednoexceptinherited

Definition at line 31 of file Observer.h.

References e, art::Observer::getTriggerResults(), art::Observer::wantAllEvents_, and art::Observer::wantEvent().

32  {
33  return wantAllEvents_;
34  }
art::Observer::wantAllEvents_
bool wantAllEvents_
Definition: Observer.h:75
bool art::Observer::wantEvent ( ScheduleID  id,
Event const &  e 
) const
protectedinherited

Definition at line 63 of file Observer.cc.

References art::Observer::rejectors_, art::Observer::selectors_, and art::Observer::wantAllEvents_.

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

64  {
65  if (wantAllEvents_) {
66  return true;
67  }
68  bool const select_event = selectors_ ? selectors_->matchEvent(id, e) : true;
69  bool const reject_event =
70  rejectors_ ? rejectors_->matchEvent(id, e) : false;
71  return select_event and not reject_event;
72  }
art::Observer::wantAllEvents_
bool wantAllEvents_
Definition: Observer.h:75
art::Observer::rejectors_
std::optional< detail::ProcessAndEventSelectors > rejectors_
Definition: Observer.h:80
e
Float_t e
Definition: plot.C:35
art::Observer::selectors_
std::optional< detail::ProcessAndEventSelectors > selectors_
Definition: Observer.h:79

Member Data Documentation

const int shower::TCShowerTemplateMaker::EBINS = 20
private

Definition at line 132 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

const double shower::TCShowerTemplateMaker::EMAX = 20.5
private

Definition at line 134 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

const double shower::TCShowerTemplateMaker::EMIN = 0.5
private

Definition at line 133 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

calo::CalorimetryAlg shower::TCShowerTemplateMaker::fCalorimetryAlg
private

Definition at line 143 of file TCShowerTemplateMaker_module.cc.

Referenced by showerProfile().

std::string shower::TCShowerTemplateMaker::fDigitModuleLabel
private

Definition at line 141 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

std::string shower::TCShowerTemplateMaker::fGenieGenModuleLabel
private

Definition at line 140 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

std::string shower::TCShowerTemplateMaker::fHitModuleLabel
private

Definition at line 138 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

TH3F* shower::TCShowerTemplateMaker::fLongitudinal
private

Definition at line 94 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_electron
private

Definition at line 103 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_other
private

Definition at line 117 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_photon
private

Definition at line 110 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

std::string shower::TCShowerTemplateMaker::fShowerModuleLabel
private

Definition at line 139 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileHitLong
private

Definition at line 59 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitLong2D
private

Definition at line 63 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileHitTrans
private

Definition at line 67 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D
private

Definition at line 71 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_1
private

Definition at line 75 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_2
private

Definition at line 79 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_3
private

Definition at line 83 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_4
private

Definition at line 87 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_5
private

Definition at line 91 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileRecoLong
private

Definition at line 60 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoLong2D
private

Definition at line 64 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans
private

Definition at line 68 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D
private

Definition at line 72 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_1
private

Definition at line 76 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_2
private

Definition at line 80 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_3
private

Definition at line 84 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_4
private

Definition at line 88 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_5
private

Definition at line 92 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileSimLong
private

Definition at line 58 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimLong2D
private

Definition at line 62 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileSimTrans
private

Definition at line 66 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D
private

Definition at line 70 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_1
private

Definition at line 74 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_2
private

Definition at line 78 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_3
private

Definition at line 82 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_4
private

Definition at line 86 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_5
private

Definition at line 90 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TH3F* shower::TCShowerTemplateMaker::fTransverse
private

Definition at line 95 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fTransverse1_electron
private

Definition at line 104 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse1_other
private

Definition at line 118 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse1_photon
private

Definition at line 111 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse2_electron
private

Definition at line 105 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse2_other
private

Definition at line 119 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse2_photon
private

Definition at line 112 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse3_electron
private

Definition at line 106 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse3_other
private

Definition at line 120 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse3_photon
private

Definition at line 113 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse4_electron
private

Definition at line 107 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse4_other
private

Definition at line 121 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse4_photon
private

Definition at line 114 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse5_electron
private

Definition at line 108 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse5_other
private

Definition at line 122 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse5_photon
private

Definition at line 115 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TH3F* shower::TCShowerTemplateMaker::fTransverse_1
private

Definition at line 96 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_2
private

Definition at line 97 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_3
private

Definition at line 98 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_4
private

Definition at line 99 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_5
private

Definition at line 100 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

const int shower::TCShowerTemplateMaker::LBINS = 20
private

Definition at line 124 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), showerProfile(), and showerProfileTrue().

const int shower::TCShowerTemplateMaker::LMAX = 5
private

Definition at line 126 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), showerProfile(), and showerProfileTrue().

const int shower::TCShowerTemplateMaker::LMIN = 0
private

Definition at line 125 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), showerProfile(), and showerProfileTrue().

const int shower::TCShowerTemplateMaker::TBINS = 20
private

Definition at line 128 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), showerProfile(), and showerProfileTrue().

const int shower::TCShowerTemplateMaker::TMAX = 5
private

Definition at line 130 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), showerProfile(), and showerProfileTrue().

const int shower::TCShowerTemplateMaker::TMIN = -5
private

Definition at line 129 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), showerProfile(), and showerProfileTrue().

const double shower::TCShowerTemplateMaker::X0 = 14
private

Definition at line 136 of file TCShowerTemplateMaker_module.cc.

Referenced by showerProfile(), and showerProfileTrue().

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