Public 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 Member Functions

 TCShowerTemplateMaker (fhicl::ParameterSet const &pset)
 
- Public Member Functions inherited from art::EDAnalyzer
 EDAnalyzer (fhicl::ParameterSet const &pset)
 
template<typename Config >
 EDAnalyzer (Table< Config > const &config)
 
std::string workerType () const
 
- Public Member Functions inherited from art::detail::Analyzer
virtual ~Analyzer () noexcept
 
 Analyzer (fhicl::ParameterSet const &pset)
 
template<typename Config >
 Analyzer (Table< Config > const &config)
 
void doBeginJob (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)
 
- Public Member Functions inherited from art::Observer
 ~Observer () noexcept
 
 Observer (Observer const &)=delete
 
 Observer (Observer &&)=delete
 
Observer & operator= (Observer const &)=delete
 
Observer & operator= (Observer &&)=delete
 
void registerProducts (ProductDescriptions &, ModuleDescription const &)
 
void fillDescriptions (ModuleDescription const &)
 
fhicl::ParameterSetID selectorConfig () const
 
- Public Member Functions inherited from art::ModuleBase
virtual ~ModuleBase () noexcept
 
 ModuleBase ()
 
ModuleDescription const & moduleDescription () const
 
void setModuleDescription (ModuleDescription const &)
 
std::array< std::vector< ProductInfo >, NumBranchTypes > const & getConsumables () const
 
void sortConsumables (std::string const &current_process_name)
 
template<typename T , BranchType BT>
ViewToken< T > consumesView (InputTag const &tag)
 
template<typename T , BranchType BT>
ViewToken< T > mayConsumeView (InputTag const &tag)
 

Private 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
 

Additional Inherited Members

- Public Types inherited from art::EDAnalyzer
using WorkerType = WorkerT< EDAnalyzer >
 
using ModuleType = EDAnalyzer
 
- Protected Member Functions inherited from art::Observer
std::string const & processName () const
 
bool wantAllEvents () const noexcept
 
bool wantEvent (ScheduleID id, Event const &e) const
 
Handle< TriggerResults > getTriggerResults (Event const &e) const
 
 Observer (fhicl::ParameterSet const &config)
 
 Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
 
- Protected Member Functions inherited from art::ModuleBase
ConsumesCollector & consumesCollector ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > consumes (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > mayConsume (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 

Detailed Description

Definition at line 36 of file TCShowerTemplateMaker_module.cc.

Constructor & Destructor Documentation

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

Definition at line 150 of file TCShowerTemplateMaker_module.cc.

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

Member Function Documentation

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

Definition at line 578 of file TCShowerTemplateMaker_module.cc.

579 {
580 
581  art::Handle<std::vector<recob::Hit>> hitListHandle;
582  std::vector<art::Ptr<recob::Hit>> hitlist;
583  if (evt.getByLabel(fHitModuleLabel, hitListHandle)) art::fill_ptr_vector(hitlist, hitListHandle);
584 
585  art::Handle<std::vector<sim::SimChannel>> scListHandle;
586  std::vector<art::Ptr<sim::SimChannel>> simchanlist;
587  if (evt.getByLabel(fDigitModuleLabel, scListHandle))
588  art::fill_ptr_vector(simchanlist, scListHandle);
589 
590  art::Handle<std::vector<recob::Shower>> showerListHandle;
591  std::vector<art::Ptr<recob::Shower>> showerlist;
592  if (evt.getByLabel(fShowerModuleLabel, showerListHandle))
593  art::fill_ptr_vector(showerlist, showerListHandle);
594 
595  art::Handle<std::vector<simb::MCTruth>> mctruthListHandle;
596  std::vector<art::Ptr<simb::MCTruth>> mclist;
597  if (evt.getByLabel(fGenieGenModuleLabel, mctruthListHandle))
598  art::fill_ptr_vector(mclist, mctruthListHandle);
599 
600  art::FindManyP<recob::Hit> shwfm(showerListHandle, evt, fShowerModuleLabel);
601 
602  if (empty(mclist)) return;
603 
604  auto const clock_data = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
605  auto const det_prop =
606  art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clock_data);
607 
608  art::Ptr<simb::MCTruth> mctruth = mclist[0];
609  if (mctruth->NeutrinoSet()) {
610  if (std::abs(mctruth->GetNeutrino().Nu().PdgCode()) == 12 &&
611  mctruth->GetNeutrino().CCNC() == 0) {
612  double elep = mctruth->GetNeutrino().Lepton().E();
613  if (showerlist.size()) {
614  std::vector<art::Ptr<recob::Hit>> showerhits = shwfm.at(0);
615  showerProfile(clock_data,
616  det_prop,
617  showerhits,
618  showerlist[0]->ShowerStart(),
619  showerlist[0]->Direction(),
620  elep);
621  }
622  showerProfileTrue(clock_data, det_prop, hitlist, elep);
623  showerProfileTrue(simchanlist, mctruth->GetNeutrino().Lepton());
624  }
625  }
626 } // analyze
simb::MCParticle::E
double E(const int i=0) const
Definition: MCParticle.h:233
simb::MCParticle::PdgCode
int PdgCode() const
Definition: MCParticle.h:212
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
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:631
abs
T abs(T value)
Definition: sparse_vector_test.cc:30
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
Definition: DataViewImpl.h:633
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:743
simb::MCNeutrino::Lepton
const simb::MCParticle & Lepton() const
Definition: MCNeutrino.h:147
art::Direction
Direction
Definition: AssnsIter.h:13
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:297
util::empty
decltype(auto) constexpr empty(T &&obj)
ADL-aware version of std::empty.
Definition: StdUtils.h:97
void shower::TCShowerTemplateMaker::beginJob ( )
overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 162 of file TCShowerTemplateMaker_module.cc.

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
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 631 of file TCShowerTemplateMaker_module.cc.

637 {
638  art::ServiceHandle<geo::Geometry const> geom;
639 
640  auto collectionPlane = geo::PlaneID(0, 0, 1);
641 
642  double shwVtxTime = detProp.ConvertXToTicks(shwvtx[0], collectionPlane);
643  double shwVtxWire = geom->WireCoordinate(shwvtx[1], shwvtx[2], collectionPlane);
644 
645  double shwTwoTime = detProp.ConvertXToTicks(shwvtx[0] + shwdir[0], collectionPlane);
646  double shwTwoWire =
647  geom->WireCoordinate(shwvtx[1] + shwdir[1], shwvtx[2] + shwdir[2], collectionPlane);
648 
649  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
650  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
651 
652  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
653  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
654  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
655  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
656  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
657 
658  for (size_t i = 0; i < showerhits.size(); ++i) {
659  if (showerhits[i]->WireID().Plane != collectionPlane.Plane) continue;
660 
661  double wirePitch = geom->WirePitch(showerhits[i]->WireID());
662  double tickToDist = detProp.DriftVelocity(detProp.Efield(), detProp.Temperature());
663  tickToDist *= 1.e-3 * sampling_rate(clockData); // 1e-3 is conversion of 1/us to 1/ns
664 
665  double xvtx = shwVtxTime * tickToDist;
666  double yvtx = shwVtxWire * wirePitch;
667 
668  double xtwo = shwTwoTime * tickToDist;
669  double ytwo = shwTwoWire * wirePitch;
670 
671  double xtwoorth = (ytwo - yvtx) + xvtx;
672  double ytwoorth = -(xtwo - xvtx) + yvtx;
673 
674  double xhit = showerhits[i]->PeakTime() * tickToDist;
675  double yhit = showerhits[i]->WireID().Wire * wirePitch;
676 
677  double ldist = std::abs((ytwoorth - yvtx) * xhit - (xtwoorth - xvtx) * yhit + xtwoorth * yvtx -
678  ytwoorth * xvtx) /
679  std::hypot(ytwoorth - yvtx, xtwoorth - xvtx);
680  double tdist = ((ytwo - yvtx) * xhit - (xtwo - xvtx) * yhit + xtwo * yvtx - ytwo * xvtx) /
681  std::hypot(ytwo - yvtx, xtwo - xvtx);
682 
683  double to3D = 1. / std::hypot(xvtx - xtwo,
684  yvtx - ytwo); // distance between two points in 3D space is one
685  ldist *= to3D;
686  tdist *= to3D;
687  double t = ldist / X0;
688 
689  double Q = showerhits[i]->Integral() *
690  fCalorimetryAlg.LifetimeCorrection(clockData, detProp, showerhits[i]->PeakTime());
691 
692  ltemp->Fill(t, Q);
693  ttemp->Fill(tdist, Q);
694 
695  if (t < 1)
696  ttemp_1->Fill(tdist, Q);
697  else if (t < 2)
698  ttemp_2->Fill(tdist, Q);
699  else if (t < 3)
700  ttemp_3->Fill(tdist, Q);
701  else if (t < 4)
702  ttemp_4->Fill(tdist, Q);
703  else if (t < 5)
704  ttemp_5->Fill(tdist, Q);
705 
706  } // loop through showerhits
707 
708  for (int i = 0; i < LBINS; ++i) {
709  if (ltemp->GetBinContent(i + 1) == 0) continue;
710  fShowerProfileRecoLong->Fill(ltemp->GetBinCenter(i + 1), ltemp->GetBinContent(i + 1));
711  fShowerProfileRecoLong2D->Fill(ltemp->GetBinCenter(i + 1), elep, ltemp->GetBinContent(i + 1));
712  fLongitudinal->Fill(ltemp->GetBinCenter(i + 1), elep, ltemp->GetBinContent(i + 1));
713  }
714 
715  for (int i = 0; i < TBINS; ++i) {
716  if (ttemp->GetBinContent(i + 1) == 0) continue;
717  fShowerProfileRecoTrans->Fill(ttemp->GetBinCenter(i + 1), ttemp->GetBinContent(i + 1));
718  fShowerProfileRecoTrans2D->Fill(ttemp->GetBinCenter(i + 1), elep, ttemp->GetBinContent(i + 1));
719  fTransverse->Fill(ttemp->GetBinCenter(i + 1), elep, ttemp->GetBinContent(i + 1));
720 
721  fShowerProfileRecoTrans2D_1->Fill(
722  ttemp_1->GetBinCenter(i + 1), elep, ttemp_1->GetBinContent(i + 1));
723  fShowerProfileRecoTrans2D_2->Fill(
724  ttemp_2->GetBinCenter(i + 1), elep, ttemp_2->GetBinContent(i + 1));
725  fShowerProfileRecoTrans2D_3->Fill(
726  ttemp_3->GetBinCenter(i + 1), elep, ttemp_3->GetBinContent(i + 1));
727  fShowerProfileRecoTrans2D_4->Fill(
728  ttemp_4->GetBinCenter(i + 1), elep, ttemp_4->GetBinContent(i + 1));
729  fShowerProfileRecoTrans2D_5->Fill(
730  ttemp_5->GetBinCenter(i + 1), elep, ttemp_5->GetBinContent(i + 1));
731 
732  fTransverse_1->Fill(ttemp_1->GetBinCenter(i + 1), elep, ttemp_1->GetBinContent(i + 1));
733  fTransverse_2->Fill(ttemp_2->GetBinCenter(i + 1), elep, ttemp_2->GetBinContent(i + 1));
734  fTransverse_3->Fill(ttemp_3->GetBinCenter(i + 1), elep, ttemp_3->GetBinContent(i + 1));
735  fTransverse_4->Fill(ttemp_4->GetBinCenter(i + 1), elep, ttemp_4->GetBinContent(i + 1));
736  fTransverse_5->Fill(ttemp_5->GetBinCenter(i + 1), elep, ttemp_5->GetBinContent(i + 1));
737  }
738 } // showerProfile
geo::GeometryCore::WireCoordinate
geo::Length_t WireCoordinate(double YPos, double ZPos, geo::PlaneID const &planeid) const
Returns the index of the nearest wire to the specified position.
Definition: GeometryCore.cxx:1128
geo::GeometryCore::WirePitch
geo::Length_t WirePitch(geo::PlaneID const &planeid) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:824
cet::hypot
std::enable_if_t< std::is_arithmetic_v< T >, T > hypot(T x, T y)
Definition: hypot.h:60
abs
T abs(T value)
Definition: sparse_vector_test.cc:30
geo::fhicl::WireID
IDparameter< geo::WireID > WireID
Member type of validated geo::WireID parameter.
Definition: geo_types_fhicl.h:407
reco_momentum_tuples.t
t
Definition: reco_momentum_tuples.py:25
CRT::PlaneID
detail::Node< FrameID, bool > PlaneID
Definition: CRTID.h:125
calo::CalorimetryAlg::LifetimeCorrection
double LifetimeCorrection(detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, double time, double T0=0) const
Definition: CalorimetryAlg.cxx:246
detinfo::sampling_rate
double sampling_rate(DetectorClocksData const &data)
Returns the period of the TPC readout electronics clock.
Definition: DetectorClocksData.h:643
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 743 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 891 of file TCShowerTemplateMaker_module.cc.

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

Member Data Documentation

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

Definition at line 131 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 133 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 132 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 142 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 140 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 139 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 137 of file TCShowerTemplateMaker_module.cc.

TH3F* shower::TCShowerTemplateMaker::fLongitudinal
private

Definition at line 93 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_electron
private

Definition at line 102 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_other
private

Definition at line 116 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_photon
private

Definition at line 109 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 138 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fShowerProfileHitLong
private

Definition at line 58 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitLong2D
private

Definition at line 62 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fShowerProfileHitTrans
private

Definition at line 66 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D
private

Definition at line 70 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_1
private

Definition at line 74 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_2
private

Definition at line 78 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_3
private

Definition at line 82 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_4
private

Definition at line 86 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_5
private

Definition at line 90 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fShowerProfileRecoLong
private

Definition at line 59 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoLong2D
private

Definition at line 63 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans
private

Definition at line 67 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D
private

Definition at line 71 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_1
private

Definition at line 75 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_2
private

Definition at line 79 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_3
private

Definition at line 83 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_4
private

Definition at line 87 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_5
private

Definition at line 91 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fShowerProfileSimLong
private

Definition at line 57 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimLong2D
private

Definition at line 61 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fShowerProfileSimTrans
private

Definition at line 65 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D
private

Definition at line 69 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_1
private

Definition at line 73 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_2
private

Definition at line 77 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_3
private

Definition at line 81 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_4
private

Definition at line 85 of file TCShowerTemplateMaker_module.cc.

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_5
private

Definition at line 89 of file TCShowerTemplateMaker_module.cc.

TH3F* shower::TCShowerTemplateMaker::fTransverse
private

Definition at line 94 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse1_electron
private

Definition at line 103 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse1_other
private

Definition at line 117 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse1_photon
private

Definition at line 110 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse2_electron
private

Definition at line 104 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse2_other
private

Definition at line 118 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse2_photon
private

Definition at line 111 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse3_electron
private

Definition at line 105 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse3_other
private

Definition at line 119 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse3_photon
private

Definition at line 112 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse4_electron
private

Definition at line 106 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse4_other
private

Definition at line 120 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse4_photon
private

Definition at line 113 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse5_electron
private

Definition at line 107 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse5_other
private

Definition at line 121 of file TCShowerTemplateMaker_module.cc.

TProfile* shower::TCShowerTemplateMaker::fTransverse5_photon
private

Definition at line 114 of file TCShowerTemplateMaker_module.cc.

TH3F* shower::TCShowerTemplateMaker::fTransverse_1
private

Definition at line 95 of file TCShowerTemplateMaker_module.cc.

TH3F* shower::TCShowerTemplateMaker::fTransverse_2
private

Definition at line 96 of file TCShowerTemplateMaker_module.cc.

TH3F* shower::TCShowerTemplateMaker::fTransverse_3
private

Definition at line 97 of file TCShowerTemplateMaker_module.cc.

TH3F* shower::TCShowerTemplateMaker::fTransverse_4
private

Definition at line 98 of file TCShowerTemplateMaker_module.cc.

TH3F* shower::TCShowerTemplateMaker::fTransverse_5
private

Definition at line 99 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 123 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 125 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 124 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 127 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 129 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 128 of file TCShowerTemplateMaker_module.cc.

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

Definition at line 135 of file TCShowerTemplateMaker_module.cc.

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