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

Public Types

using ModuleType = EDAnalyzer
 

Public Member Functions

 TCShowerElectronLikelihood (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 getShowerProfile (detinfo::DetectorClocksData const &clockdata, detinfo::DetectorPropertiesData const &detProp, std::vector< art::Ptr< recob::Hit >> showerhits, TVector3 shwvtx, TVector3 shwdir)
 
void findEnergyBin ()
 
double getLongLikelihood ()
 
double getTranLikelihood ()
 
void resetProfiles ()
 

Private Attributes

std::string fTemplateFile
 
std::string fROOTfile
 
TH3F * longTemplate
 
TH3F * tranTemplate
 
TH3F * tranTemplate_1
 
TH3F * tranTemplate_2
 
TH3F * tranTemplate_3
 
TH3F * tranTemplate_4
 
TH3F * tranTemplate_5
 
TProfile2D * longTemplateProf2D
 
TProfile2D * tranTemplateProf2D
 
TProfile2D * tranTemplateProf2D_1
 
TProfile2D * tranTemplateProf2D_2
 
TProfile2D * tranTemplateProf2D_3
 
TProfile2D * tranTemplateProf2D_4
 
TProfile2D * tranTemplateProf2D_5
 
TH1F * energyDist
 
TH1F * longLikelihoodHist
 
TH1F * tranLikelihoodHist
 
TH1F * longProfHist
 
TH1F * tranProfHist_1
 
TH1F * tranProfHist_2
 
TH1F * tranProfHist_3
 
TH1F * tranProfHist_4
 
TH1F * tranProfHist_5
 
TH1F * longProfile
 
TH1F * tranProfile
 
TH1F * tranProfile_1
 
TH1F * tranProfile_2
 
TH1F * tranProfile_3
 
TH1F * tranProfile_4
 
TH1F * tranProfile_5
 
int energyGuess
 
double energyChi2
 
int maxt
 
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 double X0 = 14
 
std::string fHitModuleLabel
 
std::string fShowerModuleLabel
 
std::string fTemplateModuleLabel
 
std::string fGenieGenModuleLabel
 
std::string fDigitModuleLabel
 
calo::CalorimetryAlg fCalorimetryAlg
 

Detailed Description

Definition at line 37 of file TCShowerElectronLikelihood_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer
inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

Definition at line 122 of file TCShowerElectronLikelihood_module.cc.

References file, fROOTfile, fTemplateFile, fhicl::ParameterSet::get(), LBINS, LMAX, LMIN, longProfile, longTemplate, longTemplateProf2D, TBINS, TMAX, TMIN, tranProfile, tranProfile_1, tranProfile_2, tranProfile_3, tranProfile_4, tranProfile_5, tranTemplate, tranTemplate_1, tranTemplate_2, tranTemplate_3, tranTemplate_4, tranTemplate_5, tranTemplateProf2D, tranTemplateProf2D_1, tranTemplateProf2D_2, tranTemplateProf2D_3, tranTemplateProf2D_4, and tranTemplateProf2D_5.

123  : EDAnalyzer(pset)
124  , fHitModuleLabel(pset.get<std::string>("HitModuleLabel", "trajcluster"))
125  , fShowerModuleLabel(pset.get<std::string>("ShowerModuleLabel", "tcshower"))
126  , fGenieGenModuleLabel(pset.get<std::string>("GenieGenModuleLabel", "generator"))
127  , fCalorimetryAlg(pset.get<fhicl::ParameterSet>("CalorimetryAlg"))
128 {
129  fTemplateFile = pset.get<std::string>("TemplateFile");
130  cet::search_path sp("FW_SEARCH_PATH");
131  if (!sp.find_file(fTemplateFile, fROOTfile))
132  throw cet::exception("TCShowerElectronLikelihood")
133  << "cannot find the root template file: \n"
134  << fTemplateFile << "\n bail ungracefully.\n";
135 
136  std::unique_ptr<TFile> file{TFile::Open(fROOTfile.c_str())};
137 
138  longTemplate = file->Get<TH3F>("tcshowertemplate/fLongitudinal");
139  tranTemplate = file->Get<TH3F>("tcshowertemplate/fTransverse");
140  tranTemplate_1 = file->Get<TH3F>("tcshowertemplate/fTransverse_1");
141  tranTemplate_2 = file->Get<TH3F>("tcshowertemplate/fTransverse_2");
142  tranTemplate_3 = file->Get<TH3F>("tcshowertemplate/fTransverse_3");
143  tranTemplate_4 = file->Get<TH3F>("tcshowertemplate/fTransverse_4");
144  tranTemplate_5 = file->Get<TH3F>("tcshowertemplate/fTransverse_5");
145 
146  longTemplateProf2D = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoLong2D");
147  tranTemplateProf2D = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D");
148  tranTemplateProf2D_1 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_1");
149  tranTemplateProf2D_2 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_2");
150  tranTemplateProf2D_3 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_3");
151  tranTemplateProf2D_4 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_4");
152  tranTemplateProf2D_5 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_5");
153 
154  longProfile = new TH1F("longProfile", "longitudinal shower profile;t;Q", LBINS, LMIN, LMAX);
155  tranProfile = new TH1F("tranProfile", "transverse shower profile;dist (cm);Q", TBINS, TMIN, TMAX);
156 
157  tranProfile_1 = new TH1F(
158  "tranProfile_1", "transverse shower profile [0 <= t < 1];dist (cm);Q", TBINS, TMIN, TMAX);
159  tranProfile_2 = new TH1F(
160  "tranProfile_2", "transverse shower profile [1 <= t < 2];dist (cm);Q", TBINS, TMIN, TMAX);
161  tranProfile_3 = new TH1F(
162  "tranProfile_3", "transverse shower profile [2 <= t < 3];dist (cm);Q", TBINS, TMIN, TMAX);
163  tranProfile_4 = new TH1F(
164  "tranProfile_4", "transverse shower profile [3 <= t < 4];dist (cm);Q", TBINS, TMIN, TMAX);
165  tranProfile_5 = new TH1F(
166  "tranProfile_5", "transverse shower profile [4 <= t < 5];dist (cm);Q", TBINS, TMIN, TMAX);
167 } // TCShowerElectronLikelihood
art::EDAnalyzer::EDAnalyzer
EDAnalyzer(fhicl::ParameterSet const &pset)
Definition: EDAnalyzer.cc:6
file
TFile * file
Definition: comparison_ascii.C:5
fhicl::exception
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33

Member Function Documentation

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

Definition at line 198 of file TCShowerElectronLikelihood_module.cc.

References util::abs(), simb::MCNeutrino::CCNC(), simb::MCParticle::E(), energyDist, energyGuess, fGenieGenModuleLabel, fHitModuleLabel, art::fill_ptr_vector(), findEnergyBin(), fShowerModuleLabel, art::ProductRetriever::getByLabel(), getLongLikelihood(), simb::MCTruth::GetNeutrino(), getShowerProfile(), getTranLikelihood(), simb::MCNeutrino::Lepton(), longLikelihoodHist, longProfHist, longProfile, maxt, simb::MCTruth::NeutrinoSet(), simb::MCNeutrino::Nu(), simb::MCParticle::PdgCode(), resetProfiles(), tranLikelihoodHist, tranProfHist_1, tranProfHist_2, tranProfHist_3, tranProfHist_4, tranProfHist_5, tranProfile_1, tranProfile_2, tranProfile_3, tranProfile_4, tranProfile_5, X0, and Z.

199 {
200  resetProfiles();
201 
202  art::Handle<std::vector<recob::Hit>> hitListHandle;
203  std::vector<art::Ptr<recob::Hit>> hitlist;
204  if (evt.getByLabel(fHitModuleLabel, hitListHandle)) art::fill_ptr_vector(hitlist, hitListHandle);
205 
206  art::Handle<std::vector<recob::Shower>> showerListHandle;
207  std::vector<art::Ptr<recob::Shower>> showerlist;
208  if (evt.getByLabel(fShowerModuleLabel, showerListHandle))
209  art::fill_ptr_vector(showerlist, showerListHandle);
210 
211  art::Handle<std::vector<simb::MCTruth>> mctruthListHandle;
212  std::vector<art::Ptr<simb::MCTruth>> mclist;
213  if (evt.getByLabel(fGenieGenModuleLabel, mctruthListHandle))
214  art::fill_ptr_vector(mclist, mctruthListHandle);
215 
216  art::FindManyP<recob::Hit> shwfm(showerListHandle, evt, fShowerModuleLabel);
217 
218  if (showerlist.size()) {
219  auto const clock_data =
220  art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
221  auto const det_prop =
222  art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clock_data);
223  std::vector<art::Ptr<recob::Hit>> showerhits = shwfm.at(0);
224  getShowerProfile(
225  clock_data, det_prop, showerhits, showerlist[0]->ShowerStart(), showerlist[0]->Direction());
226 
227  maxt = std::ceil((90 - showerlist[0]->ShowerStart().Z()) / X0);
228 
229  findEnergyBin();
230 
231  longLikelihoodHist->Fill(getLongLikelihood());
232  tranLikelihoodHist->Fill(getTranLikelihood());
233 
234  // check true shower energy
235  if (mclist.size()) {
236  art::Ptr<simb::MCTruth> mctruth = mclist[0];
237  if (mctruth->NeutrinoSet()) {
238  if (std::abs(mctruth->GetNeutrino().Nu().PdgCode()) == 12 &&
239  mctruth->GetNeutrino().CCNC() == 0) {
240  double elep = mctruth->GetNeutrino().Lepton().E();
241  energyDist->Fill(elep - energyGuess);
242  } // cc nue
243  } // neutrinoset
244  } // mclist
245  } // showerlist
246 
247  longProfHist = longProfile;
248  tranProfHist_1 = tranProfile_1;
249  tranProfHist_2 = tranProfile_2;
250  tranProfHist_3 = tranProfile_3;
251  tranProfHist_4 = tranProfile_4;
252  tranProfHist_5 = tranProfile_5;
253 } // 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
shower::TCShowerElectronLikelihood::getShowerProfile
void getShowerProfile(detinfo::DetectorClocksData const &clockdata, detinfo::DetectorPropertiesData const &detProp, std::vector< art::Ptr< recob::Hit >> showerhits, TVector3 shwvtx, TVector3 shwdir)
Definition: TCShowerElectronLikelihood_module.cc:275
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
Z
Float_t Z
Definition: plot.C:37
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
void shower::TCShowerElectronLikelihood::beginJob ( )
overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 171 of file TCShowerElectronLikelihood_module.cc.

References energyDist, LBINS, LMAX, LMIN, longLikelihoodHist, longProfHist, TBINS, TMAX, TMIN, tranLikelihoodHist, tranProfHist_1, tranProfHist_2, tranProfHist_3, tranProfHist_4, and tranProfHist_5.

172 {
173 
174  art::ServiceHandle<art::TFileService const> tfs;
175 
176  energyDist = tfs->make<TH1F>("energyDist", "true energy - guess energy", 41, -20.5, 20.5);
177  longLikelihoodHist = tfs->make<TH1F>("longLikelihoodHist", "longitudinal likelihood", 20, 0, 2);
178  tranLikelihoodHist = tfs->make<TH1F>("tranLikelihoodHist", "transverse likelihood", 20, 0, 3);
179 
180  // just for printing purposes
181  longProfHist =
182  tfs->make<TH1F>("longProfHist", "longitudinal e- profile (reco);t;Q", LBINS, LMIN, LMAX);
183  tranProfHist_1 = tfs->make<TH1F>(
184  "tranProfHist", "transverse e- profile (reco) [0 <= t < 1];dist (cm);Q", TBINS, TMIN, TMAX);
185  tranProfHist_2 = tfs->make<TH1F>(
186  "tranProfHist", "transverse e- profile (reco) [1 <= t < 2];dist (cm);Q", TBINS, TMIN, TMAX);
187  tranProfHist_3 = tfs->make<TH1F>(
188  "tranProfHist", "transverse e- profile (reco) [2 <= t < 3];dist (cm);Q", TBINS, TMIN, TMAX);
189  tranProfHist_4 = tfs->make<TH1F>(
190  "tranProfHist", "transverse e- profile (reco) [3 <= t < 4];dist (cm);Q", TBINS, TMIN, TMAX);
191  tranProfHist_5 = tfs->make<TH1F>(
192  "tranProfHist", "transverse e- profile (reco) [4 <= t < 5];dist (cm);Q", TBINS, TMIN, TMAX);
193 
194 } // 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")
void shower::TCShowerElectronLikelihood::findEnergyBin ( )
private

Definition at line 352 of file TCShowerElectronLikelihood_module.cc.

References energyChi2, energyGuess, longProfile, longTemplate, longTemplateProf2D, tranProfile, tranProfile_1, tranProfile_2, tranProfile_3, tranProfile_4, tranProfile_5, tranTemplate, tranTemplateProf2D_1, tranTemplateProf2D_2, tranTemplateProf2D_3, tranTemplateProf2D_4, and tranTemplateProf2D_5.

Referenced by analyze().

353 {
354 
355  if (longProfile->GetNbinsX() != longTemplate->GetNbinsX())
356  throw cet::exception("TCShowerElectronLikelihood")
357  << "Bin mismatch in longitudinal profile template \n";
358 
359  if (tranProfile->GetNbinsX() != tranTemplate->GetNbinsX())
360  throw cet::exception("TCShowerElectronLikelihood")
361  << "Bin mismatch in transverse profile template \n";
362 
363  double chi2min = 999999;
364  double bestbin = -1;
365 
366  int ebins = longTemplate->GetNbinsY();
367  int lbins = longTemplate->GetNbinsX();
368  int tbins = tranTemplate->GetNbinsX();
369 
370  for (int i = 0; i < ebins; ++i) {
371  double thischi2 = 0;
372 
373  TProfile* ltemp = longTemplateProf2D->ProfileX("_x", i + 1, i + 1);
374  TProfile* ttemp_1 = tranTemplateProf2D_1->ProfileX("_x_1", i + 1, i + 1);
375  TProfile* ttemp_2 = tranTemplateProf2D_2->ProfileX("_x_2", i + 1, i + 1);
376  TProfile* ttemp_3 = tranTemplateProf2D_3->ProfileX("_x_3", i + 1, i + 1);
377  TProfile* ttemp_4 = tranTemplateProf2D_4->ProfileX("_x_4", i + 1, i + 1);
378  TProfile* ttemp_5 = tranTemplateProf2D_5->ProfileX("_x_5", i + 1, i + 1);
379 
380  int nlbins = 0;
381  int ntbins = 0;
382 
383  for (int j = 0; j < lbins; ++j) {
384  double obs = longProfile->GetBinContent(j + 1);
385  double exp = ltemp->GetBinContent(j + 1);
386  if (obs != 0) {
387  thischi2 += cet::square(obs - exp) / exp;
388  ++nlbins;
389  }
390  } // loop through longitudinal bins
391 
392  for (int j = 0; j < tbins; ++j) {
393  double obs = tranProfile_1->GetBinContent(j + 1);
394  double exp = ttemp_1->GetBinContent(j + 1);
395  if (obs != 0) {
396  thischi2 += cet::square(obs - exp) / exp;
397  ++ntbins;
398  }
399 
400  obs = tranProfile_2->GetBinContent(j + 1);
401  exp = ttemp_2->GetBinContent(j + 1);
402  if (obs != 0) {
403  thischi2 += cet::square(obs - exp) / exp;
404  ++ntbins;
405  }
406 
407  obs = tranProfile_3->GetBinContent(j + 1);
408  exp = ttemp_3->GetBinContent(j + 1);
409  if (obs != 0) {
410  thischi2 += cet::square(obs - exp) / exp;
411  ++ntbins;
412  }
413 
414  obs = tranProfile_4->GetBinContent(j + 1);
415  exp = ttemp_4->GetBinContent(j + 1);
416  if (obs != 0) {
417  thischi2 += cet::square(obs - exp) / exp;
418  ++ntbins;
419  }
420 
421  obs = tranProfile_5->GetBinContent(j + 1);
422  exp = ttemp_5->GetBinContent(j + 1);
423  if (obs != 0) {
424  thischi2 += cet::square(obs - exp) / exp;
425  ++ntbins;
426  }
427  } // loop through longitudinal bins
428 
429  thischi2 /= (nlbins + ntbins);
430 
431  if (thischi2 < chi2min) {
432  chi2min = thischi2;
433  bestbin = i;
434  }
435 
436  } // loop through energy bins
437 
438  energyChi2 = chi2min;
439  energyGuess = bestbin + 1;
440 } // findEnergyBin
fhicl::exception
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
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
double shower::TCShowerElectronLikelihood::getLongLikelihood ( )
private

Definition at line 444 of file TCShowerElectronLikelihood_module.cc.

References energyGuess, LBINS, longProfile, and longTemplate.

Referenced by analyze().

445 {
446  if (energyGuess < 0) return -9999.;
447  int energyBin = energyGuess;
448 
449  double longLikelihood = 0.;
450  int nbins = 0;
451 
452  for (int i = 0; i < LBINS; ++i) {
453  double qval = longProfile->GetBinContent(i + 1);
454  int qbin = longTemplate->GetZaxis()->FindBin(qval);
455  int binentries = longTemplate->GetBinContent(i + 1, energyBin, qbin);
456  int totentries = longTemplate->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
457  if (qval > 0) {
458  ++nbins;
459  double prob = (double)binentries / totentries * 100;
460  if (binentries > 0) longLikelihood += log(prob);
461  }
462  } // loop through
463 
464  longLikelihood /= nbins;
465 
466  std::cout << longLikelihood << std::endl;
467  return longLikelihood;
468 } // getLongLikelihood
void shower::TCShowerElectronLikelihood::getShowerProfile ( detinfo::DetectorClocksData const &  clockdata,
detinfo::DetectorPropertiesData const &  detProp,
std::vector< art::Ptr< recob::Hit >>  showerhits,
TVector3  shwvtx,
TVector3  shwdir 
)
private

Definition at line 275 of file TCShowerElectronLikelihood_module.cc.

References util::abs(), detinfo::DetectorPropertiesData::ConvertXToTicks(), detinfo::DetectorPropertiesData::DriftVelocity(), detinfo::DetectorPropertiesData::Efield(), fCalorimetryAlg, calo::CalorimetryAlg::LifetimeCorrection(), longProfile, detinfo::sampling_rate(), detinfo::DetectorPropertiesData::Temperature(), geo::vect::toPoint(), geo::vect::toVector(), tranProfile, tranProfile_1, tranProfile_2, tranProfile_3, tranProfile_4, tranProfile_5, geo::GeometryCore::WireCoordinate(), geo::GeometryCore::WirePitch(), and X0.

Referenced by analyze().

281 {
282  art::ServiceHandle<geo::Geometry const> geom;
283 
284  auto collectionPlane = geo::PlaneID(0, 0, 1);
285 
286  using namespace geo::vect;
287  auto const shwvtx_p = toPoint(shwvtx);
288  auto const shwvtx_p2 = shwvtx_p + toVector(shwdir);
289 
290  double shwVtxTime = detProp.ConvertXToTicks(shwvtx_p.X(), collectionPlane);
291  double shwVtxWire = geom->WireCoordinate(shwvtx_p, collectionPlane);
292 
293  double shwTwoTime = detProp.ConvertXToTicks(shwvtx_p2.X(), collectionPlane);
294  double shwTwoWire = geom->WireCoordinate(shwvtx_p2, collectionPlane);
295 
296  for (size_t i = 0; i < showerhits.size(); ++i) {
297  if (showerhits[i]->WireID().Plane != collectionPlane.Plane) continue;
298 
299  double wirePitch = geom->WirePitch(showerhits[i]->WireID());
300  double tickToDist = detProp.DriftVelocity(detProp.Efield(), detProp.Temperature());
301  tickToDist *= 1.e-3 * sampling_rate(clockdata); // 1e-3 is conversion of 1/us to 1/ns
302 
303  double xvtx = shwVtxTime * tickToDist;
304  double yvtx = shwVtxWire * wirePitch;
305 
306  double xtwo = shwTwoTime * tickToDist;
307  double ytwo = shwTwoWire * wirePitch;
308 
309  double xtwoorth = (ytwo - yvtx) + xvtx;
310  double ytwoorth = -(xtwo - xvtx) + yvtx;
311 
312  double xhit = showerhits[i]->PeakTime() * tickToDist;
313  double yhit = showerhits[i]->WireID().Wire * wirePitch;
314 
315  double ldist = std::abs((ytwoorth - yvtx) * xhit - (xtwoorth - xvtx) * yhit + xtwoorth * yvtx -
316  ytwoorth * xvtx) /
317  std::hypot(ytwoorth - yvtx, xtwoorth - xvtx);
318  double tdist = ((ytwo - yvtx) * xhit - (xtwo - xvtx) * yhit + xtwo * yvtx - ytwo * xvtx) /
319  std::hypot(ytwo - yvtx, xtwo - xvtx);
320 
321  double to3D = 1. / std::hypot(xvtx - xtwo,
322  yvtx - ytwo); // distance between two points in 3D space is one
323  ldist *= to3D;
324  tdist *= to3D;
325  double t = ldist / X0;
326 
327  double Q = showerhits[i]->Integral() *
328  fCalorimetryAlg.LifetimeCorrection(clockdata, detProp, showerhits[i]->PeakTime());
329 
330  longProfile->Fill(t, Q);
331  tranProfile->Fill(tdist, Q);
332 
333  if (t < 1)
334  tranProfile_1->Fill(tdist, Q);
335  else if (t < 2)
336  tranProfile_2->Fill(tdist, Q);
337  else if (t < 3)
338  tranProfile_3->Fill(tdist, Q);
339  else if (t < 4)
340  tranProfile_4->Fill(tdist, Q);
341  else if (t < 5)
342  tranProfile_5->Fill(tdist, Q);
343 
344  } // loop through showerhits
345 
346  return;
347 
348 } // getShowerProfile
geo::GeometryCore::WireCoordinate
Length_t WireCoordinate(Point_t const &pos, PlaneID const &planeid) const
Returns the index of the nearest wire to the specified position.
Definition: GeometryCore.cxx:1093
geo::vect::toPoint
::geo::Point_t toPoint(Point const &p)
Convert the specified point into a geo::Point_t.
Definition: geo_vectors_utils.h:1584
util::abs
constexpr auto abs(T v)
Returns the absolute value of the argument.
Definition: constexpr_math.h:40
geo::vect::toVector
::geo::Vector_t toVector(Vector const &v)
Convert the specified vector into a geo::Vector_t.
Definition: geo_vectors_utils.h:1591
geo::fhicl::PlaneID
IDparameter< geo::PlaneID > PlaneID
Member type of validated geo::PlaneID parameter.
Definition: geo_types_fhicl.h:355
geo::vect
Utilities to manipulate geometry vectors.The utilities include generic vector interface facilities al...
Definition: geo_vectors_utils.h:55
geo::fhicl::WireID
IDparameter< geo::WireID > WireID
Member type of validated geo::WireID parameter.
Definition: geo_types_fhicl.h:383
calo::CalorimetryAlg::LifetimeCorrection
double LifetimeCorrection(detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, double time, double T0=0) const
Definition: CalorimetryAlg.cxx:221
detinfo::sampling_rate
double sampling_rate(DetectorClocksData const &data)
Returns the period of the TPC readout electronics clock.
Definition: DetectorClocksData.h:575
geo::GeometryCore::WirePitch
Length_t WirePitch(PlaneID const &planeid=plane_zero) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:763
double shower::TCShowerElectronLikelihood::getTranLikelihood ( )
private

Definition at line 472 of file TCShowerElectronLikelihood_module.cc.

References DEFINE_ART_MODULE, energyGuess, TBINS, tranProfile_1, tranProfile_2, tranProfile_3, tranProfile_4, tranProfile_5, tranTemplate_1, tranTemplate_2, tranTemplate_3, tranTemplate_4, and tranTemplate_5.

Referenced by analyze().

473 {
474  if (energyGuess < 0) return -9999.;
475  int energyBin = energyGuess;
476 
477  double tranLikelihood_1 = 0;
478  double tranLikelihood_2 = 0;
479  double tranLikelihood_3 = 0;
480  double tranLikelihood_4 = 0;
481  double tranLikelihood_5 = 0;
482 
483  double qval;
484  int qbin, binentries, totentries;
485 
486  int nbins = 0;
487 
488  for (int i = 0; i < TBINS; ++i) {
489  qval = tranProfile_1->GetBinContent(i + 1);
490  qbin = tranTemplate_1->GetZaxis()->FindBin(qval);
491  binentries = tranTemplate_1->GetBinContent(i + 1, energyBin, qbin);
492  totentries = tranTemplate_1->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
493  if (qval > 0) {
494  ++nbins;
495  double prob = (double)binentries / totentries * 100;
496  if (binentries > 0) tranLikelihood_1 += log(prob);
497  }
498 
499  qval = tranProfile_2->GetBinContent(i + 1);
500  qbin = tranTemplate_2->GetZaxis()->FindBin(qval);
501  binentries = tranTemplate_2->GetBinContent(i + 1, energyBin, qbin);
502  totentries = tranTemplate_2->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
503  if (qval > 0) {
504  ++nbins;
505  double prob = (double)binentries / totentries * 100;
506  if (binentries > 0) tranLikelihood_2 += log(prob);
507  }
508 
509  qval = tranProfile_3->GetBinContent(i + 1);
510  qbin = tranTemplate_3->GetZaxis()->FindBin(qval);
511  binentries = tranTemplate_3->GetBinContent(i + 1, energyBin, qbin);
512  totentries = tranTemplate_3->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
513  if (qval > 0) {
514  ++nbins;
515  double prob = (double)binentries / totentries * 100;
516  if (binentries > 0) tranLikelihood_3 += log(prob);
517  }
518 
519  qval = tranProfile_4->GetBinContent(i + 1);
520  qbin = tranTemplate_4->GetZaxis()->FindBin(qval);
521  binentries = tranTemplate_4->GetBinContent(i + 1, energyBin, qbin);
522  totentries = tranTemplate_4->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
523  if (qval > 0) {
524  ++nbins;
525  double prob = (double)binentries / totentries * 100;
526  if (binentries > 0) tranLikelihood_4 += log(prob);
527  }
528 
529  qval = tranProfile_5->GetBinContent(i + 1);
530  qbin = tranTemplate_5->GetZaxis()->FindBin(qval);
531  binentries = tranTemplate_5->GetBinContent(i + 1, energyBin, qbin);
532  totentries = tranTemplate_5->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
533  if (qval > 0) {
534  ++nbins;
535  double prob = (double)binentries / totentries * 100;
536  if (binentries > 0) tranLikelihood_5 += log(prob);
537  }
538 
539  } // loop through
540 
541  double const tranLikelihood =
542  (tranLikelihood_1 + tranLikelihood_2 + tranLikelihood_3 + tranLikelihood_4 + tranLikelihood_5) /
543  nbins;
544  std::cout << tranLikelihood << std::endl;
545  return tranLikelihood;
546 } // getTranLikelihood
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 shower::TCShowerElectronLikelihood::resetProfiles ( )
private

Definition at line 257 of file TCShowerElectronLikelihood_module.cc.

References energyChi2, energyGuess, longProfile, maxt, tranProfile, tranProfile_1, tranProfile_2, tranProfile_3, tranProfile_4, and tranProfile_5.

Referenced by analyze().

258 {
259 
260  longProfile->Reset();
261  tranProfile->Reset();
262  tranProfile_1->Reset();
263  tranProfile_2->Reset();
264  tranProfile_3->Reset();
265  tranProfile_4->Reset();
266  tranProfile_5->Reset();
267 
268  energyGuess = -9999;
269  energyChi2 = -9999;
270  maxt = -9999;
271 } // resetProfiles
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 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

double shower::TCShowerElectronLikelihood::energyChi2
private

Definition at line 95 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), and resetProfiles().

TH1F* shower::TCShowerElectronLikelihood::energyDist
private

Definition at line 75 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and beginJob().

int shower::TCShowerElectronLikelihood::energyGuess
private

Definition at line 94 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), findEnergyBin(), getLongLikelihood(), getTranLikelihood(), and resetProfiles().

calo::CalorimetryAlg shower::TCShowerElectronLikelihood::fCalorimetryAlg
private

Definition at line 114 of file TCShowerElectronLikelihood_module.cc.

Referenced by getShowerProfile().

std::string shower::TCShowerElectronLikelihood::fDigitModuleLabel
private

Definition at line 112 of file TCShowerElectronLikelihood_module.cc.

std::string shower::TCShowerElectronLikelihood::fGenieGenModuleLabel
private

Definition at line 111 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze().

std::string shower::TCShowerElectronLikelihood::fHitModuleLabel
private

Definition at line 108 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze().

std::string shower::TCShowerElectronLikelihood::fROOTfile
private

Definition at line 57 of file TCShowerElectronLikelihood_module.cc.

Referenced by TCShowerElectronLikelihood().

std::string shower::TCShowerElectronLikelihood::fShowerModuleLabel
private

Definition at line 109 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze().

std::string shower::TCShowerElectronLikelihood::fTemplateFile
private

Definition at line 56 of file TCShowerElectronLikelihood_module.cc.

Referenced by TCShowerElectronLikelihood().

std::string shower::TCShowerElectronLikelihood::fTemplateModuleLabel
private

Definition at line 110 of file TCShowerElectronLikelihood_module.cc.

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

Definition at line 98 of file TCShowerElectronLikelihood_module.cc.

Referenced by beginJob(), getLongLikelihood(), and TCShowerElectronLikelihood().

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

Definition at line 100 of file TCShowerElectronLikelihood_module.cc.

Referenced by beginJob(), and TCShowerElectronLikelihood().

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

Definition at line 99 of file TCShowerElectronLikelihood_module.cc.

Referenced by beginJob(), and TCShowerElectronLikelihood().

TH1F* shower::TCShowerElectronLikelihood::longLikelihoodHist
private

Definition at line 76 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and beginJob().

TH1F* shower::TCShowerElectronLikelihood::longProfHist
private

Definition at line 79 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and beginJob().

TH1F* shower::TCShowerElectronLikelihood::longProfile
private

Definition at line 86 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), findEnergyBin(), getLongLikelihood(), getShowerProfile(), resetProfiles(), and TCShowerElectronLikelihood().

TH3F* shower::TCShowerElectronLikelihood::longTemplate
private

Definition at line 59 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), getLongLikelihood(), and TCShowerElectronLikelihood().

TProfile2D* shower::TCShowerElectronLikelihood::longTemplateProf2D
private

Definition at line 66 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), and TCShowerElectronLikelihood().

int shower::TCShowerElectronLikelihood::maxt
private

Definition at line 96 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and resetProfiles().

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

Definition at line 102 of file TCShowerElectronLikelihood_module.cc.

Referenced by beginJob(), getTranLikelihood(), and TCShowerElectronLikelihood().

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

Definition at line 104 of file TCShowerElectronLikelihood_module.cc.

Referenced by beginJob(), and TCShowerElectronLikelihood().

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

Definition at line 103 of file TCShowerElectronLikelihood_module.cc.

Referenced by beginJob(), and TCShowerElectronLikelihood().

TH1F* shower::TCShowerElectronLikelihood::tranLikelihoodHist
private

Definition at line 77 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and beginJob().

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_1
private

Definition at line 80 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and beginJob().

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_2
private

Definition at line 81 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and beginJob().

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_3
private

Definition at line 82 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and beginJob().

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_4
private

Definition at line 83 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and beginJob().

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_5
private

Definition at line 84 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and beginJob().

TH1F* shower::TCShowerElectronLikelihood::tranProfile
private

Definition at line 87 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), getShowerProfile(), resetProfiles(), and TCShowerElectronLikelihood().

TH1F* shower::TCShowerElectronLikelihood::tranProfile_1
private

Definition at line 88 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), findEnergyBin(), getShowerProfile(), getTranLikelihood(), resetProfiles(), and TCShowerElectronLikelihood().

TH1F* shower::TCShowerElectronLikelihood::tranProfile_2
private

Definition at line 89 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), findEnergyBin(), getShowerProfile(), getTranLikelihood(), resetProfiles(), and TCShowerElectronLikelihood().

TH1F* shower::TCShowerElectronLikelihood::tranProfile_3
private

Definition at line 90 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), findEnergyBin(), getShowerProfile(), getTranLikelihood(), resetProfiles(), and TCShowerElectronLikelihood().

TH1F* shower::TCShowerElectronLikelihood::tranProfile_4
private

Definition at line 91 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), findEnergyBin(), getShowerProfile(), getTranLikelihood(), resetProfiles(), and TCShowerElectronLikelihood().

TH1F* shower::TCShowerElectronLikelihood::tranProfile_5
private

Definition at line 92 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), findEnergyBin(), getShowerProfile(), getTranLikelihood(), resetProfiles(), and TCShowerElectronLikelihood().

TH3F* shower::TCShowerElectronLikelihood::tranTemplate
private

Definition at line 60 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), and TCShowerElectronLikelihood().

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_1
private

Definition at line 61 of file TCShowerElectronLikelihood_module.cc.

Referenced by getTranLikelihood(), and TCShowerElectronLikelihood().

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_2
private

Definition at line 62 of file TCShowerElectronLikelihood_module.cc.

Referenced by getTranLikelihood(), and TCShowerElectronLikelihood().

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_3
private

Definition at line 63 of file TCShowerElectronLikelihood_module.cc.

Referenced by getTranLikelihood(), and TCShowerElectronLikelihood().

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_4
private

Definition at line 64 of file TCShowerElectronLikelihood_module.cc.

Referenced by getTranLikelihood(), and TCShowerElectronLikelihood().

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_5
private

Definition at line 65 of file TCShowerElectronLikelihood_module.cc.

Referenced by getTranLikelihood(), and TCShowerElectronLikelihood().

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D
private

Definition at line 67 of file TCShowerElectronLikelihood_module.cc.

Referenced by TCShowerElectronLikelihood().

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_1
private

Definition at line 68 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), and TCShowerElectronLikelihood().

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_2
private

Definition at line 69 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), and TCShowerElectronLikelihood().

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_3
private

Definition at line 70 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), and TCShowerElectronLikelihood().

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_4
private

Definition at line 71 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), and TCShowerElectronLikelihood().

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_5
private

Definition at line 72 of file TCShowerElectronLikelihood_module.cc.

Referenced by findEnergyBin(), and TCShowerElectronLikelihood().

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

Definition at line 106 of file TCShowerElectronLikelihood_module.cc.

Referenced by analyze(), and getShowerProfile().

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