LArSoft  v10_04_05
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 42 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 156 of file TCShowerTemplateMaker_module.cc.

157  : EDAnalyzer(pset)
158  , fHitModuleLabel(pset.get<std::string>("HitModuleLabel", "trajcluster"))
159  , fShowerModuleLabel(pset.get<std::string>("ShowerModuleLabel", "tcshower"))
160  , fGenieGenModuleLabel(pset.get<std::string>("GenieGenModuleLabel", "generator"))
161  , fDigitModuleLabel(pset.get<std::string>("DigitModuleLabel", "largeant"))
162  , fCalorimetryAlg(pset.get<fhicl::ParameterSet>("CalorimetryAlg"))
163 {} // 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 582 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().

583 {
584 
585  art::Handle<std::vector<recob::Hit>> hitListHandle;
586  std::vector<art::Ptr<recob::Hit>> hitlist;
587  if (evt.getByLabel(fHitModuleLabel, hitListHandle)) art::fill_ptr_vector(hitlist, hitListHandle);
588 
589  art::Handle<std::vector<sim::SimChannel>> scListHandle;
590  std::vector<art::Ptr<sim::SimChannel>> simchanlist;
591  if (evt.getByLabel(fDigitModuleLabel, scListHandle))
592  art::fill_ptr_vector(simchanlist, scListHandle);
593 
594  art::Handle<std::vector<recob::Shower>> showerListHandle;
595  std::vector<art::Ptr<recob::Shower>> showerlist;
596  if (evt.getByLabel(fShowerModuleLabel, showerListHandle))
597  art::fill_ptr_vector(showerlist, showerListHandle);
598 
599  art::Handle<std::vector<simb::MCTruth>> mctruthListHandle;
600  std::vector<art::Ptr<simb::MCTruth>> mclist;
601  if (evt.getByLabel(fGenieGenModuleLabel, mctruthListHandle))
602  art::fill_ptr_vector(mclist, mctruthListHandle);
603 
604  art::FindManyP<recob::Hit> shwfm(showerListHandle, evt, fShowerModuleLabel);
605 
606  if (empty(mclist)) return;
607 
608  auto const clock_data = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
609  auto const det_prop =
610  art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clock_data);
611 
612  art::Ptr<simb::MCTruth> mctruth = mclist[0];
613  if (mctruth->NeutrinoSet()) {
614  if (std::abs(mctruth->GetNeutrino().Nu().PdgCode()) == 12 &&
615  mctruth->GetNeutrino().CCNC() == 0) {
616  double elep = mctruth->GetNeutrino().Lepton().E();
617  if (showerlist.size()) {
618  std::vector<art::Ptr<recob::Hit>> showerhits = shwfm.at(0);
619  showerProfile(clock_data,
620  det_prop,
621  showerhits,
622  showerlist[0]->ShowerStart(),
623  showerlist[0]->Direction(),
624  elep);
625  }
626  showerProfileTrue(clock_data, det_prop, hitlist, elep);
627  showerProfileTrue(simchanlist, mctruth->GetNeutrino().Lepton());
628  }
629  }
630 } // 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:634
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:744
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 167 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.

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

Referenced by analyze().

640 {
641  auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
642  auto const& collectionPlane = wireReadoutGeom.Plane(collectionPlaneID);
643 
644  double shwVtxTime = detProp.ConvertXToTicks(shwvtx[0], collectionPlaneID);
645  using geo::vect::toPoint;
646  double shwVtxWire = collectionPlane.WireCoordinate(toPoint(shwvtx));
647 
648  double shwTwoTime = detProp.ConvertXToTicks(shwvtx[0] + shwdir[0], collectionPlaneID);
649  double shwTwoWire = collectionPlane.WireCoordinate(toPoint(shwvtx + shwdir));
650 
651  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
652  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
653 
654  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
655  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
656  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
657  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
658  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
659 
660  for (size_t i = 0; i < showerhits.size(); ++i) {
661  if (showerhits[i]->WireID().Plane != collectionPlaneID.Plane) continue;
662 
663  double wirePitch = wireReadoutGeom.Plane(showerhits[i]->WireID()).WirePitch();
664  double tickToDist = detProp.DriftVelocity(detProp.Efield(), detProp.Temperature());
665  tickToDist *= 1.e-3 * sampling_rate(clockData); // 1e-3 is conversion of 1/us to 1/ns
666 
667  double xvtx = shwVtxTime * tickToDist;
668  double yvtx = shwVtxWire * wirePitch;
669 
670  double xtwo = shwTwoTime * tickToDist;
671  double ytwo = shwTwoWire * wirePitch;
672 
673  double xtwoorth = (ytwo - yvtx) + xvtx;
674  double ytwoorth = -(xtwo - xvtx) + yvtx;
675 
676  double xhit = showerhits[i]->PeakTime() * tickToDist;
677  double yhit = showerhits[i]->WireID().Wire * wirePitch;
678 
679  double ldist = std::abs((ytwoorth - yvtx) * xhit - (xtwoorth - xvtx) * yhit + xtwoorth * yvtx -
680  ytwoorth * xvtx) /
681  std::hypot(ytwoorth - yvtx, xtwoorth - xvtx);
682  double tdist = ((ytwo - yvtx) * xhit - (xtwo - xvtx) * yhit + xtwo * yvtx - ytwo * xvtx) /
683  std::hypot(ytwo - yvtx, xtwo - xvtx);
684 
685  double to3D = 1. / std::hypot(xvtx - xtwo,
686  yvtx - ytwo); // distance between two points in 3D space is one
687  ldist *= to3D;
688  tdist *= to3D;
689  double t = ldist / X0;
690 
691  double Q = showerhits[i]->Integral() *
692  fCalorimetryAlg.LifetimeCorrection(clockData, detProp, showerhits[i]->PeakTime());
693 
694  ltemp->Fill(t, Q);
695  ttemp->Fill(tdist, Q);
696 
697  if (t < 1)
698  ttemp_1->Fill(tdist, Q);
699  else if (t < 2)
700  ttemp_2->Fill(tdist, Q);
701  else if (t < 3)
702  ttemp_3->Fill(tdist, Q);
703  else if (t < 4)
704  ttemp_4->Fill(tdist, Q);
705  else if (t < 5)
706  ttemp_5->Fill(tdist, Q);
707 
708  } // loop through showerhits
709 
710  for (int i = 0; i < LBINS; ++i) {
711  if (ltemp->GetBinContent(i + 1) == 0) continue;
712  fShowerProfileRecoLong->Fill(ltemp->GetBinCenter(i + 1), ltemp->GetBinContent(i + 1));
713  fShowerProfileRecoLong2D->Fill(ltemp->GetBinCenter(i + 1), elep, ltemp->GetBinContent(i + 1));
714  fLongitudinal->Fill(ltemp->GetBinCenter(i + 1), elep, ltemp->GetBinContent(i + 1));
715  }
716 
717  for (int i = 0; i < TBINS; ++i) {
718  if (ttemp->GetBinContent(i + 1) == 0) continue;
719  fShowerProfileRecoTrans->Fill(ttemp->GetBinCenter(i + 1), ttemp->GetBinContent(i + 1));
720  fShowerProfileRecoTrans2D->Fill(ttemp->GetBinCenter(i + 1), elep, ttemp->GetBinContent(i + 1));
721  fTransverse->Fill(ttemp->GetBinCenter(i + 1), elep, ttemp->GetBinContent(i + 1));
722 
723  fShowerProfileRecoTrans2D_1->Fill(
724  ttemp_1->GetBinCenter(i + 1), elep, ttemp_1->GetBinContent(i + 1));
725  fShowerProfileRecoTrans2D_2->Fill(
726  ttemp_2->GetBinCenter(i + 1), elep, ttemp_2->GetBinContent(i + 1));
727  fShowerProfileRecoTrans2D_3->Fill(
728  ttemp_3->GetBinCenter(i + 1), elep, ttemp_3->GetBinContent(i + 1));
729  fShowerProfileRecoTrans2D_4->Fill(
730  ttemp_4->GetBinCenter(i + 1), elep, ttemp_4->GetBinContent(i + 1));
731  fShowerProfileRecoTrans2D_5->Fill(
732  ttemp_5->GetBinCenter(i + 1), elep, ttemp_5->GetBinContent(i + 1));
733 
734  fTransverse_1->Fill(ttemp_1->GetBinCenter(i + 1), elep, ttemp_1->GetBinContent(i + 1));
735  fTransverse_2->Fill(ttemp_2->GetBinCenter(i + 1), elep, ttemp_2->GetBinContent(i + 1));
736  fTransverse_3->Fill(ttemp_3->GetBinCenter(i + 1), elep, ttemp_3->GetBinContent(i + 1));
737  fTransverse_4->Fill(ttemp_4->GetBinCenter(i + 1), elep, ttemp_4->GetBinContent(i + 1));
738  fTransverse_5->Fill(ttemp_5->GetBinCenter(i + 1), elep, ttemp_5->GetBinContent(i + 1));
739  }
740 } // showerProfile
util::abs
constexpr auto abs(T v)
Returns the absolute value of the argument.
Definition: constexpr_math.h:40
Get
cout<< "Opened file "<< fin<< " ixs= "<< ixs<< endl;if(ixs==0) hhh=(TH1F *) fff-> Get("h1")
Definition: AddMC.C:8
geo::fhicl::WireID
IDparameter< geo::WireID > WireID
Member type of validated geo::WireID parameter.
Definition: geo_types_fhicl.h:383
geo::vect::toPoint
Point_t toPoint(Point const &p)
Convert the specified point into a geo::Point_t.
Definition: geo_vectors_utils.h:1562
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
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 744 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, Get, cheat::BackTrackerService::HitToEveTrackIDEs(), LBINS, LMAX, LMIN, detinfo::sampling_rate(), TBINS, detinfo::DetectorPropertiesData::Temperature(), TMAX, TMIN, cheat::ParticleInventoryService::TrackIdToParticle_P(), and X0.

Referenced by analyze().

749 {
750  auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
751  auto const& collectionPlane = wireReadoutGeom.Plane(collectionPlaneID);
752  art::ServiceHandle<cheat::BackTrackerService const> btserv;
753  art::ServiceHandle<cheat::ParticleInventoryService const> piserv;
754  std::map<int, double> trkID_E;
755 
756  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
757  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
758 
759  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
760  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
761  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
762  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
763  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
764 
765  double xvtx = -999;
766  double yvtx = -999;
767  double zvtx = -999;
768  double xtwo = -999;
769  double ytwo = -999;
770  double ztwo = -999;
771  double shwvtxT = -999;
772  double shwvtxW = -999;
773  double shwtwoT = -999;
774  double shwtwoW = -999;
775 
776  double shwvtxx = -999;
777  double shwvtxy = -999;
778  double shwtwox = -999;
779  double shwtwoy = -999;
780  double xtwoorth = -999;
781  double ytwoorth = -999;
782 
783  double wirePitch = -999;
784  double tickToDist = -999;
785 
786  bool foundParent = false;
787 
788  for (auto const& hit : allhits) {
789  if (hit->WireID().Plane != collectionPlaneID.Plane) continue;
790 
791  // art::Ptr<recob::Hit> hit = allhits[i];
792  std::vector<sim::TrackIDE> trackIDs = btserv->HitToEveTrackIDEs(clockData, hit);
793 
794  for (size_t j = 0; j < trackIDs.size(); ++j) {
795  // only want energy associated with the electron and electron must have neutrino mother
796  if (std::abs((piserv->TrackIdToParticle_P(trackIDs[j].trackID))->PdgCode()) != 11) continue;
797  if (std::abs((piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Mother()) != 0) continue;
798 
799  if (!foundParent) {
800  xvtx = (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Vx();
801  yvtx = (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Vy();
802  zvtx = (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Vz();
803 
804  xtwo = xvtx + (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Px();
805  ytwo = yvtx + (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Py();
806  ztwo = zvtx + (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Pz();
807 
808  shwvtxT = detProp.ConvertXToTicks(xvtx, collectionPlaneID);
809  geo::Point_t const vtx{xvtx, yvtx, zvtx};
810  shwvtxW = collectionPlane.WireCoordinate(vtx);
811 
812  shwtwoT = detProp.ConvertXToTicks(xtwo, collectionPlaneID);
813  geo::Point_t const vtwo{xtwo, ytwo, ztwo};
814  shwtwoW = collectionPlane.WireCoordinate(vtwo);
815 
816  wirePitch = wireReadoutGeom.Plane(hit->WireID()).WirePitch();
817  tickToDist = detProp.DriftVelocity(detProp.Efield(), detProp.Temperature());
818  tickToDist *= 1.e-3 * sampling_rate(clockData); // 1e-3 is conversion of 1/us to 1/ns
819 
820  shwvtxx = shwvtxT * tickToDist;
821  shwvtxy = shwvtxW * wirePitch;
822 
823  shwtwox = shwtwoT * tickToDist;
824  shwtwoy = shwtwoW * wirePitch;
825 
826  xtwoorth = (shwtwoy - shwvtxy) + shwvtxx;
827  ytwoorth = -(shwtwox - shwvtxx) + shwvtxy;
828 
829  foundParent = true;
830  }
831  double xhit = hit->PeakTime() * tickToDist;
832  double yhit = hit->WireID().Wire * wirePitch;
833 
834  double ldist = abs((ytwoorth - shwvtxy) * xhit - (xtwoorth - shwvtxx) * yhit +
835  xtwoorth * shwvtxy - ytwoorth * shwvtxx) /
836  std::hypot(ytwoorth - shwvtxy, xtwoorth - shwvtxx);
837  double tdist = ((shwtwoy - shwvtxy) * xhit - (shwtwox - shwvtxx) * yhit + shwtwox * shwvtxy -
838  shwtwoy * shwvtxx) /
839  std::hypot(shwtwoy - shwvtxy, shwtwox - shwvtxx);
840 
841  double to3D = std::hypot(xvtx - xtwo, yvtx - ytwo, zvtx - ztwo) /
842  std::hypot(shwvtxx - shwtwox,
843  shwvtxy - shwtwoy); // distance between two points in 3D space is one
844  ldist *= to3D;
845  tdist *= to3D;
846  double t = ldist / X0;
847 
848  double energy = trackIDs[j].energy;
849 
850  ltemp->Fill(t, energy);
851  ttemp->Fill(tdist, energy);
852 
853  if (t < 1)
854  ttemp_1->Fill(tdist, energy);
855  else if (t < 2)
856  ttemp_2->Fill(tdist, energy);
857  else if (t < 3)
858  ttemp_3->Fill(tdist, energy);
859  else if (t < 4)
860  ttemp_4->Fill(tdist, energy);
861  else if (t < 5)
862  ttemp_5->Fill(tdist, energy);
863 
864  } // loop through track IDE
865 
866  } // loop through all hits
867 
868  for (int i = 0; i < LBINS; ++i) {
869  if (ltemp->GetBinContent(i + 1) == 0) continue;
870  fShowerProfileHitLong->Fill(ltemp->GetBinCenter(i + 1), ltemp->GetBinContent(i + 1));
871  fShowerProfileHitLong2D->Fill(ltemp->GetBinCenter(i + 1), elep, ltemp->GetBinContent(i + 1));
872  }
873 
874  for (int i = 0; i < TBINS; ++i) {
875  if (ttemp->GetBinContent(i + 1) == 0) continue;
876  fShowerProfileHitTrans->Fill(ttemp->GetBinCenter(i + 1), ttemp->GetBinContent(i + 1));
877  fShowerProfileHitTrans2D->Fill(ttemp->GetBinCenter(i + 1), elep, ttemp->GetBinContent(i + 1));
878 
879  fShowerProfileHitTrans2D_1->Fill(
880  ttemp_1->GetBinCenter(i + 1), elep, ttemp_1->GetBinContent(i + 1));
881  fShowerProfileHitTrans2D_2->Fill(
882  ttemp_2->GetBinCenter(i + 1), elep, ttemp_2->GetBinContent(i + 1));
883  fShowerProfileHitTrans2D_3->Fill(
884  ttemp_3->GetBinCenter(i + 1), elep, ttemp_3->GetBinContent(i + 1));
885  fShowerProfileHitTrans2D_4->Fill(
886  ttemp_4->GetBinCenter(i + 1), elep, ttemp_4->GetBinContent(i + 1));
887  fShowerProfileHitTrans2D_5->Fill(
888  ttemp_5->GetBinCenter(i + 1), elep, ttemp_5->GetBinContent(i + 1));
889  }
890 } // showerProfileTrue
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
Get
cout<< "Opened file "<< fin<< " ixs= "<< ixs<< endl;if(ixs==0) hhh=(TH1F *) fff-> Get("h1")
Definition: AddMC.C:8
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:162
detinfo::sampling_rate
double sampling_rate(DetectorClocksData const &data)
Returns the period of the TPC readout electronics clock.
Definition: DetectorClocksData.h:575
void shower::TCShowerTemplateMaker::showerProfileTrue ( std::vector< art::Ptr< sim::SimChannel >>  allchan,
simb::MCParticle  electron 
)
private

Definition at line 894 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, Get, 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(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), simb::MCParticle::Vz(), x, X0, x2, y, y2, and z.

897 {
898  art::ServiceHandle<cheat::ParticleInventoryService const> piserv;
899  auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
900 
901  std::vector<sim::MCEnDep> alledep;
902 
903  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
904  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
905 
906  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
907  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
908  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
909  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
910  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
911 
912  // get electron energy depositions
913  for (size_t i = 0; i < allchan.size(); ++i) {
914  art::Ptr<sim::SimChannel> simchan = allchan[i];
915  if (wireReadoutGeom.View(simchan->Channel()) != geo::kV) continue;
916  auto tdc_ide_map = simchan->TDCIDEMap();
917 
918  for (auto const& tdc_ide_pair : tdc_ide_map) {
919  for (auto const& ide : tdc_ide_pair.second) {
920  if (piserv->TrackIdToMotherParticle_P(ide.trackID) == nullptr) continue;
921  if (std::abs(piserv->TrackIdToMotherParticle_P(ide.trackID)->PdgCode()) != 11) continue;
922 
923  sim::MCEnDep edep;
924  edep.SetVertex(ide.x, ide.y, ide.z);
925  edep.SetEnergy(ide.energy);
926  edep.SetTrackId(ide.trackID);
927 
928  alledep.push_back(edep);
929 
930  } // loop through ide
931  } // loop through tdc_ide
932 
933  } // loop through channels
934 
935  double x0 = electron.Vx();
936  double y0 = electron.Vy();
937  double z0 = electron.Vz();
938 
939  double x2 = electron.Px();
940  double y2 = electron.Py();
941  double z2 = electron.Pz();
942 
943  TVector3 v0(x2, y2, z2);
944  v0 = v0.Unit();
945 
946  for (size_t i = 0; i < alledep.size(); ++i) {
947  double x = (double)alledep[i].Vertex()[0];
948  double y = (double)alledep[i].Vertex()[1];
949  double z = (double)alledep[i].Vertex()[2];
950 
951  TVector3 v1(x - x0, y - y0, z - z0);
952 
953  double ldist = v0.Dot(v1);
954  double t = ldist / X0;
955  double tdist = (v0.Orthogonal()).Dot(v1);
956 
957  double energy = alledep[i].Energy();
958 
959  ltemp->Fill(t, energy);
960  ttemp->Fill(tdist, energy);
961 
962  if (t < 1)
963  ttemp_1->Fill(tdist, energy);
964  else if (t < 2)
965  ttemp_2->Fill(tdist, energy);
966  else if (t < 3)
967  ttemp_3->Fill(tdist, energy);
968  else if (t < 4)
969  ttemp_4->Fill(tdist, energy);
970  else if (t < 5)
971  ttemp_5->Fill(tdist, energy);
972  }
973 
974  for (int i = 0; i < LBINS; ++i) {
975  if (ltemp->GetBinContent(i + 1) == 0) continue;
976  fShowerProfileSimLong->Fill(ltemp->GetBinCenter(i + 1), ltemp->GetBinContent(i + 1));
977  fShowerProfileSimLong2D->Fill(
978  ltemp->GetBinCenter(i + 1), electron.E(), ltemp->GetBinContent(i + 1));
979  }
980 
981  for (int i = 0; i < TBINS; ++i) {
982  if (ttemp->GetBinContent(i + 1) == 0) continue;
983  fShowerProfileSimTrans->Fill(ttemp->GetBinCenter(i + 1), ttemp->GetBinContent(i + 1));
984  fShowerProfileSimTrans2D->Fill(
985  ttemp->GetBinCenter(i + 1), electron.E(), ttemp->GetBinContent(i + 1));
986 
987  fShowerProfileSimTrans2D_1->Fill(
988  ttemp_1->GetBinCenter(i + 1), electron.E(), ttemp_1->GetBinContent(i + 1));
989  fShowerProfileSimTrans2D_2->Fill(
990  ttemp_2->GetBinCenter(i + 1), electron.E(), ttemp_2->GetBinContent(i + 1));
991  fShowerProfileSimTrans2D_3->Fill(
992  ttemp_3->GetBinCenter(i + 1), electron.E(), ttemp_3->GetBinContent(i + 1));
993  fShowerProfileSimTrans2D_4->Fill(
994  ttemp_4->GetBinCenter(i + 1), electron.E(), ttemp_4->GetBinContent(i + 1));
995  fShowerProfileSimTrans2D_5->Fill(
996  ttemp_5->GetBinCenter(i + 1), electron.E(), ttemp_5->GetBinContent(i + 1));
997  }
998 } // 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:132
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
Get
cout<< "Opened file "<< fin<< " ixs= "<< ixs<< endl;if(ixs==0) hhh=(TH1F *) fff-> Get("h1")
Definition: AddMC.C:8
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
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 137 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

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

Definition at line 139 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

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

Definition at line 138 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

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

Definition at line 148 of file TCShowerTemplateMaker_module.cc.

Referenced by showerProfile().

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

Definition at line 146 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

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

Definition at line 145 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

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

Definition at line 143 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

TH3F* shower::TCShowerTemplateMaker::fLongitudinal
private

Definition at line 99 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_electron
private

Definition at line 108 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_other
private

Definition at line 122 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_photon
private

Definition at line 115 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

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

Definition at line 144 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileHitLong
private

Definition at line 64 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitLong2D
private

Definition at line 68 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileHitTrans
private

Definition at line 72 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D
private

Definition at line 76 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_1
private

Definition at line 80 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_2
private

Definition at line 84 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_3
private

Definition at line 88 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_4
private

Definition at line 92 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_5
private

Definition at line 96 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileRecoLong
private

Definition at line 65 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoLong2D
private

Definition at line 69 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans
private

Definition at line 73 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D
private

Definition at line 77 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_1
private

Definition at line 81 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_2
private

Definition at line 85 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_3
private

Definition at line 89 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_4
private

Definition at line 93 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_5
private

Definition at line 97 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileSimLong
private

Definition at line 63 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimLong2D
private

Definition at line 67 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileSimTrans
private

Definition at line 71 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D
private

Definition at line 75 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_1
private

Definition at line 79 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_2
private

Definition at line 83 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_3
private

Definition at line 87 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_4
private

Definition at line 91 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_5
private

Definition at line 95 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TH3F* shower::TCShowerTemplateMaker::fTransverse
private

Definition at line 100 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fTransverse1_electron
private

Definition at line 109 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse1_other
private

Definition at line 123 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse1_photon
private

Definition at line 116 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse2_electron
private

Definition at line 110 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse2_other
private

Definition at line 124 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse2_photon
private

Definition at line 117 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse3_electron
private

Definition at line 111 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse3_other
private

Definition at line 125 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse3_photon
private

Definition at line 118 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse4_electron
private

Definition at line 112 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse4_other
private

Definition at line 126 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse4_photon
private

Definition at line 119 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse5_electron
private

Definition at line 113 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse5_other
private

Definition at line 127 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse5_photon
private

Definition at line 120 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TH3F* shower::TCShowerTemplateMaker::fTransverse_1
private

Definition at line 101 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_2
private

Definition at line 102 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_3
private

Definition at line 103 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_4
private

Definition at line 104 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_5
private

Definition at line 105 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

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

Definition at line 129 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 131 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 130 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 133 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 135 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 134 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 141 of file TCShowerTemplateMaker_module.cc.

Referenced by showerProfile(), and showerProfileTrue().

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