LArSoft  v07_13_02
Liquid Argon Software toolkit - http://larsoft.org/
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
shower::TCShowerTemplateMaker Class Reference
Inheritance diagram for shower::TCShowerTemplateMaker:
art::EDAnalyzer art::EventObserverBase art::Consumer art::EngineCreator

Public Types

using WorkerType = WorkerT< EDAnalyzer >
 
using ModuleType = EDAnalyzer
 

Public Member Functions

 TCShowerTemplateMaker (fhicl::ParameterSet const &pset)
 
virtual ~TCShowerTemplateMaker ()
 
void reconfigure (fhicl::ParameterSet const &pset)
 
void beginJob ()
 
void analyze (const art::Event &evt)
 
void showerProfile (std::vector< art::Ptr< recob::Hit > > showerhits, TVector3 shwvtx, TVector3 shwdir, double elep)
 
void showerProfileTrue (std::vector< art::Ptr< recob::Hit > > allhits, double elep)
 
void showerProfileTrue (std::vector< art::Ptr< sim::SimChannel > > allchan, simb::MCParticle electron)
 
std::string workerType () const
 
bool modifiesEvent () const
 
void registerProducts (MasterProductRegistry &, ProductDescriptions &, ModuleDescription const &)
 
std::string const & processName () const
 
bool wantAllEvents () const
 
bool wantEvent (Event const &e)
 
fhicl::ParameterSetID selectorConfig () const
 
art::Handle< art::TriggerResults > getTriggerResults (Event const &e) const
 
template<typename T , BranchType = InEvent>
ProductToken< T > consumes (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< T > consumes (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< T > consumesView (InputTag const &it)
 
template<typename T , BranchType = InEvent>
ProductToken< T > mayConsume (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< T > mayConsume (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< T > mayConsumeView (InputTag const &it)
 
base_engine_t & createEngine (seed_t seed)
 
base_engine_t & createEngine (seed_t seed, std::string const &kind_of_engine_to_make)
 
base_engine_t & createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
 
seed_t get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)
 

Static Public Member Functions

static cet::exempt_ptr< Consumer > non_module_context ()
 

Protected Member Functions

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

Protected Attributes

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

Detailed Description

Definition at line 49 of file TCShowerTemplateMaker_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer
inherited

Definition at line 39 of file EDAnalyzer.h.

using art::EDAnalyzer::WorkerType = WorkerT<EDAnalyzer>
inherited

Definition at line 38 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

Definition at line 166 of file TCShowerTemplateMaker_module.cc.

References reconfigure().

166  :
167  EDAnalyzer(pset),
168  fHitModuleLabel (pset.get< std::string >("HitModuleLabel", "trajcluster" ) ),
169  fShowerModuleLabel (pset.get< std::string >("ShowerModuleLabel", "tcshower" ) ),
170  fGenieGenModuleLabel (pset.get< std::string >("GenieGenModuleLabel", "generator") ),
171  fDigitModuleLabel (pset.get< std::string >("DigitModuleLabel", "largeant") ),
172  fCalorimetryAlg (pset.get< fhicl::ParameterSet >("CalorimetryAlg") ) {
173  this->reconfigure(pset);
174 } // TCShowerTemplateMaker
art::EDAnalyzer::EDAnalyzer
EDAnalyzer(Table< Config > const &config)
Definition: EDAnalyzer.h:100
shower::TCShowerTemplateMaker::reconfigure
void reconfigure(fhicl::ParameterSet const &pset)
Definition: TCShowerTemplateMaker_module.cc:183
shower::TCShowerTemplateMaker::~TCShowerTemplateMaker ( )
virtual

Definition at line 178 of file TCShowerTemplateMaker_module.cc.

178  {
179 } // ~TCShowerTemplateMaker

Member Function Documentation

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

Definition at line 267 of file TCShowerTemplateMaker_module.cc.

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

267  {
268 
269  art::Handle< std::vector<recob::Hit> > hitListHandle;
270  std::vector<art::Ptr<recob::Hit> > hitlist;
271  if (evt.getByLabel(fHitModuleLabel,hitListHandle))
272  art::fill_ptr_vector(hitlist, hitListHandle);
273 
274  art::Handle< std::vector<sim::SimChannel> > scListHandle;
275  std::vector<art::Ptr<sim::SimChannel> > simchanlist;
276  if (evt.getByLabel(fDigitModuleLabel,scListHandle))
277  art::fill_ptr_vector(simchanlist, scListHandle);
278 
279  art::Handle< std::vector<recob::Shower> > showerListHandle;
280  std::vector<art::Ptr<recob::Shower> > showerlist;
281  if (evt.getByLabel(fShowerModuleLabel,showerListHandle))
282  art::fill_ptr_vector(showerlist, showerListHandle);
283 
284  art::Handle< std::vector<simb::MCTruth> > mctruthListHandle;
285  std::vector<art::Ptr<simb::MCTruth> > mclist;
286  if (evt.getByLabel(fGenieGenModuleLabel,mctruthListHandle))
287  art::fill_ptr_vector(mclist, mctruthListHandle);
288 
289  art::FindManyP<recob::Hit> shwfm(showerListHandle, evt, fShowerModuleLabel);
290 
291  if (mclist.size()) {
292  art::Ptr<simb::MCTruth> mctruth = mclist[0];
293  if (mctruth->NeutrinoSet()) {
294  if (std::abs(mctruth->GetNeutrino().Nu().PdgCode()) == 12 && mctruth->GetNeutrino().CCNC() == 0) {
295  double elep = mctruth->GetNeutrino().Lepton().E();
296  if (showerlist.size()) {
297  std::vector< art::Ptr<recob::Hit> > showerhits = shwfm.at(0);
298  showerProfile(showerhits, showerlist[0]->ShowerStart(), showerlist[0]->Direction(), elep);
299  }
300  showerProfileTrue(hitlist, elep);
301  showerProfileTrue(simchanlist, mctruth->GetNeutrino().Lepton());
302  }
303  // use this to look at background
304  /*
305  if (!(std::abs(mctruth->GetNeutrino().Nu().PdgCode()) == 12 && mctruth->GetNeutrino().CCNC() == 0)) {
306  if (showerlist.size()) {
307  std::vector< art::Ptr<recob::Hit> > showerhits = shwfm.at(0);
308  showerProfile(showerhits, showerlist[0]->ShowerStart(), showerlist[0]->Direction(), 0);
309  }
310  }
311  */
312 
313  }
314  }
315 
316  //fTree->Fill();
317 
318 } // analyze
simb::MCParticle::E
double E(const int i=0) const
Definition: MCParticle.h:237
simb::MCParticle::PdgCode
int PdgCode() const
Definition: MCParticle.h:216
simb::MCNeutrino::CCNC
int CCNC() const
Definition: MCNeutrino.h:152
simb::MCTruth::GetNeutrino
const simb::MCNeutrino & GetNeutrino() const
Definition: MCTruth.h:74
simb::MCNeutrino::Nu
const simb::MCParticle & Nu() const
Definition: MCNeutrino.h:150
simb::MCNeutrino::Lepton
const simb::MCParticle & Lepton() const
Definition: MCNeutrino.h:151
shower::TCShowerTemplateMaker::showerProfile
void showerProfile(std::vector< art::Ptr< recob::Hit > > showerhits, TVector3 shwvtx, TVector3 shwdir, double elep)
Definition: TCShowerTemplateMaker_module.cc:322
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
Definition: DataViewImpl.h:344
art::Direction
Direction
Definition: AssnsIter.h:24
simb::MCTruth::NeutrinoSet
bool NeutrinoSet() const
Definition: MCTruth.h:75
art::fill_ptr_vector
void fill_ptr_vector(std::vector< Ptr< T >> &ptrs, H const &h)
Definition: Ptr.h:464
shower::TCShowerTemplateMaker::showerProfileTrue
void showerProfileTrue(std::vector< art::Ptr< recob::Hit > > allhits, double elep)
Definition: TCShowerTemplateMaker_module.cc:417
void shower::TCShowerTemplateMaker::beginJob ( )
virtual

Reimplemented from art::EDAnalyzer.

Definition at line 188 of file TCShowerTemplateMaker_module.cc.

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

188  {
189 
190  art::ServiceHandle<art::TFileService> tfs;
191 
192  // fTree = tfs->make<TTree>("tcshowerana", "tcshowerana");
193 
194  fShowerProfileSimLong = tfs->make<TProfile>("fShowerProfileSimLong", "longitudinal e- profile (true, simchannel);t;E (MeV)", LBINS, LMIN, LMAX);
195  fShowerProfileHitLong = tfs->make<TProfile>("fShowerProfileHitLong", "longitudinal e- profile (true, hit);t;E (MeV)", LBINS, LMIN, LMAX);
196  fShowerProfileRecoLong = tfs->make<TProfile>("fShowerProfileRecoLong", "longitudinal e- profile (reco);t;Q", LBINS, LMIN, LMAX);
197 
198  fShowerProfileSimLong2D = tfs->make<TProfile2D>("fShowerProfileSimLong2D", "longitudinal e- profile (true, simchannel);t;electron energy (MeV);E (MeV)", LBINS, LMIN, LMAX, EBINS, EMIN, EMAX);
199  fShowerProfileHitLong2D = tfs->make<TProfile2D>("fShowerProfileHitLong2D", "longitudinal e- profile (true, hit);t;electron energy (MeV);E (MeV)", LBINS, LMIN, LMAX, EBINS, EMIN, EMAX);
200  fShowerProfileRecoLong2D = tfs->make<TProfile2D>("fShowerProfileRecoLong2D", "longitudinal e- profile (reco);t;electron energy (MeV);Q", LBINS, LMIN, LMAX, EBINS, EMIN, EMAX);
201 
202  fShowerProfileSimTrans = tfs->make<TProfile>("fShowerProfileSimTrans", "transverse e- profile (true, simchannel);dist (cm);E (MeV)", TBINS, TMIN, TMAX);
203  fShowerProfileHitTrans = tfs->make<TProfile>("fShowerProfileHitTrans", "transverse e- profile (true, hit);dist (cm);E (MeV)", TBINS, TMIN, TMAX);
204  fShowerProfileRecoTrans = tfs->make<TProfile>("fShowerProfileRecoTrans", "transverse e- profile (reco);dist (cm);Q", TBINS, TMIN, TMAX);
205 
206  fShowerProfileSimTrans2D = tfs->make<TProfile2D>("fShowerProfileSimTrans2D", "transverse e- profile (true, simchannel);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
207  fShowerProfileHitTrans2D = tfs->make<TProfile2D>("fShowerProfileHitTrans2D", "transverse e- profile (true, hit);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
208  fShowerProfileRecoTrans2D = tfs->make<TProfile2D>("fShowerProfileRecoTrans2D", "transverse e- profile (reco);t;electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
209 
210  fShowerProfileSimTrans2D_1 = tfs->make<TProfile2D>("fShowerProfileSimTrans2D_1", "transverse e- profile [0 <= t < 1] (true, simchannel);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
211  fShowerProfileHitTrans2D_1 = tfs->make<TProfile2D>("fShowerProfileHitTrans2D_1", "transverse e- profile [0 <= t < 1] (true, hit);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
212  fShowerProfileRecoTrans2D_1 = tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_1", "transverse e- profile [0 <= t < 1] (reco);t;electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
213 
214  fShowerProfileSimTrans2D_2 = tfs->make<TProfile2D>("fShowerProfileSimTrans2D_2", "transverse e- profile [1 <= t < 2] (true, simchannel);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
215  fShowerProfileHitTrans2D_2 = tfs->make<TProfile2D>("fShowerProfileHitTrans2D_2", "transverse e- profile [1 <= t < 2] (true, hit);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
216  fShowerProfileRecoTrans2D_2 = tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_2", "transverse e- profile [1 <= t < 2] (reco);t;electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
217 
218  fShowerProfileSimTrans2D_3 = tfs->make<TProfile2D>("fShowerProfileSimTrans2D_3", "transverse e- profile [2 <= t < 3] (true, simchannel);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
219  fShowerProfileHitTrans2D_3 = tfs->make<TProfile2D>("fShowerProfileHitTrans2D_3", "transverse e- profile [2 <= t < 3] (true, hit);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
220  fShowerProfileRecoTrans2D_3 = tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_3", "transverse e- profile [2 <= t < 3] (reco);t;electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
221 
222  fShowerProfileSimTrans2D_4 = tfs->make<TProfile2D>("fShowerProfileSimTrans2D_4", "transverse e- profile [3 <= t < 4] (true, simchannel);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
223  fShowerProfileHitTrans2D_4 = tfs->make<TProfile2D>("fShowerProfileHitTrans2D_4", "transverse e- profile [3 <= t < 4] (true, hit);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
224  fShowerProfileRecoTrans2D_4 = tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_4", "transverse e- profile [3 <= t < 4] (reco);t;electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
225 
226  fShowerProfileSimTrans2D_5 = tfs->make<TProfile2D>("fShowerProfileSimTrans2D_5", "transverse e- profile [4 <= t < 5] (true, simchannel);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
227  fShowerProfileHitTrans2D_5 = tfs->make<TProfile2D>("fShowerProfileHitTrans2D_5", "transverse e- profile [4 <= t < 5] (true, hit);t;electron energy (MeV);E (MeV)", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
228  fShowerProfileRecoTrans2D_5 = tfs->make<TProfile2D>("fShowerProfileRecoTrans2D_5", "transverse e- profile [4 <= t < 5] (reco);t;electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX);
229 
230  fLongitudinal = tfs->make<TH3F>("fLongitudinal", "longitudinal e- profile;t;electron energy (MeV);Q", LBINS, LMIN, LMAX, EBINS, EMIN, EMAX, 100, 0, 150000);
231  fTransverse = tfs->make<TH3F>("fTransverse", "transverse e- profile;dist (cm);electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX, 100, 0, 150000);
232  fTransverse_1 = tfs->make<TH3F>("fTransverse_1", "transverse e- profile [0 <= t < 1];dist (cm);electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX, 100, 0, 100000);
233  fTransverse_2 = tfs->make<TH3F>("fTransverse_2", "transverse e- profile [1 <= t < 2];dist (cm);electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX, 100, 0, 100000);
234  fTransverse_3 = tfs->make<TH3F>("fTransverse_3", "transverse e- profile [2 <= t < 3];dist (cm);electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX, 100, 0, 100000);
235  fTransverse_4 = tfs->make<TH3F>("fTransverse_4", "transverse e- profile [3 <= t < 4];dist (cm);electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX, 100, 0, 100000);
236  fTransverse_5 = tfs->make<TH3F>("fTransverse_5", "transverse e- profile [4 <= t < 5];dist (cm);electron energy (MeV);Q", TBINS, TMIN, TMAX, EBINS, EMIN, EMAX, 100, 0, 100000);
237 
238  // electrons
239  fLongitudinal_electron = tfs->make<TProfile>("fLongitudinal_electron", "longitudinal e- profile (reco);t;Q", LBINS, LMIN, LMAX);
240  fTransverse1_electron = tfs->make<TProfile>("fTransverse1_electron", "transverse e- profile [0 <= t < 1] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
241  fTransverse2_electron = tfs->make<TProfile>("fTransverse2_electron", "transverse e- profile [1 <= t < 2] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
242  fTransverse3_electron = tfs->make<TProfile>("fTransverse3_electron", "transverse e- profile [2 <= t < 3] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
243  fTransverse4_electron = tfs->make<TProfile>("fTransverse4_electron", "transverse e- profile [3 <= t < 4] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
244  fTransverse5_electron = tfs->make<TProfile>("fTransverse5_electron", "transverse e- profile [4 <= t < 5] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
245 
246  // photons
247  fLongitudinal_photon = tfs->make<TProfile>("fLongitudinal_photon", "longitudinal photon profile (reco);t;Q", LBINS, LMIN, LMAX);
248  fTransverse1_photon = tfs->make<TProfile>("fTransverse1_photon", "transverse photon profile [0 <= t < 1] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
249  fTransverse2_photon = tfs->make<TProfile>("fTransverse2_photon", "transverse photon profile [1 <= t < 2] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
250  fTransverse3_photon = tfs->make<TProfile>("fTransverse3_photon", "transverse photon profile [2 <= t < 3] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
251  fTransverse4_photon = tfs->make<TProfile>("fTransverse4_photon", "transverse photon profile [3 <= t < 4] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
252  fTransverse5_photon = tfs->make<TProfile>("fTransverse5_photon", "transverse photon profile [4 <= t < 5] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
253 
254  // other
255  fLongitudinal_other = tfs->make<TProfile>("fLongitudinal_other", "longitudinal other profile (reco);t;Q", LBINS, LMIN, LMAX);
256  fTransverse1_other = tfs->make<TProfile>("fTransverse1_other", "transverse other profile [0 <= t < 1] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
257  fTransverse2_other = tfs->make<TProfile>("fTransverse2_other", "transverse other profile [1 <= t < 2] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
258  fTransverse3_other = tfs->make<TProfile>("fTransverse3_other", "transverse other profile [2 <= t < 3] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
259  fTransverse4_other = tfs->make<TProfile>("fTransverse4_other", "transverse other profile [3 <= t < 4] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
260  fTransverse5_other = tfs->make<TProfile>("fTransverse5_other", "transverse other profile [4 <= t < 5] (reco);dist (cm);Q", TBINS, TMIN, TMAX);
261 
262 
263 } // beginJob
art::TFileDirectory::make
T * make(ARGS...args) const
Definition: TFileDirectory.h:82
detail::CachedProducts& art::EventObserverBase::cachedProducts ( )
inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

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

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

Definition at line 147 of file Consumer.h.

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

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

Definition at line 162 of file Consumer.h.

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

Definition at line 172 of file Consumer.h.

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

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

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

27 {
28  return rng()->createEngine(placeholder_schedule_id(), seed);
29 }
seed
long seed
Definition: chem4.cc:68
art::EngineCreator::rng
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
Definition: EngineCreator.cc:58
EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed,
std::string const &  kind_of_engine_to_make 
)
inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

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

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

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

Definition at line 114 of file EDAnalyzer.cc.

References art::EDAnalyzer::current_context_.

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

115  {
116  return current_context_.get();
117  }
art::EDAnalyzer::current_context_
CPC_exempt_ptr current_context_
Definition: EDAnalyzer.h:179
EngineCreator::seed_t EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset,
char const  key[] = "seed",
seed_t const  implicit_seed = -1 
)
inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

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

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

Definition at line 61 of file EventObserverBase.h.

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

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

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

Definition at line 190 of file Consumer.h.

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

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

Definition at line 205 of file Consumer.h.

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

Definition at line 215 of file Consumer.h.

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

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

Definition at line 25 of file EventObserverBase.h.

26  {
27  return false;
28  }
cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )
staticinherited

Definition at line 76 of file Consumer.cc.

Referenced by art::RootOutput::beginSubRun(), art::OutputModule::doBeginRun(), art::OutputModule::doBeginSubRun(), art::OutputModule::doEndRun(), art::OutputModule::doEndSubRun(), art::ProducingService::doPostReadEvent(), art::ProducingService::doPostReadRun(), art::ProducingService::doPostReadSubRun(), art::OutputModule::doWriteEvent(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EventProcessor::readEvent(), art::EventProcessor::readRun(), art::EmptyEvent::readRun_(), art::EventProcessor::readSubRun(), and art::EmptyEvent::readSubRun_().

77 {
78  static Consumer invalid{InvalidTag{}};
79  return &invalid;
80 }
void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset)
protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

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

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

Definition at line 41 of file EventObserverBase.h.

References art::EventObserverBase::process_name_.

Referenced by art::FileDumperOutput::printPrincipal(), and art::RootOutput::RootOutput().

42  {
43  return process_name_;
44  }
void shower::TCShowerTemplateMaker::reconfigure ( fhicl::ParameterSet const &  pset)

Definition at line 183 of file TCShowerTemplateMaker_module.cc.

Referenced by TCShowerTemplateMaker().

183  {
184 } // reconfigure
void art::EventObserverBase::registerProducts ( MasterProductRegistry &  ,
ProductDescriptions &  ,
ModuleDescription const &   
)
inlineinherited

Definition at line 33 of file EventObserverBase.h.

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

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

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

57  {
58  return selector_config_id_;
59  }
art::EventObserverBase::selector_config_id_
fhicl::ParameterSetID selector_config_id_
Definition: EventObserverBase.h:93
void shower::TCShowerTemplateMaker::showerProfile ( std::vector< art::Ptr< recob::Hit > >  showerhits,
TVector3  shwvtx,
TVector3  shwdir,
double  elep 
)

Definition at line 322 of file TCShowerTemplateMaker_module.cc.

References fCalorimetryAlg, fLongitudinal, fShowerProfileRecoLong, fShowerProfileRecoLong2D, fShowerProfileRecoTrans, fShowerProfileRecoTrans2D, fShowerProfileRecoTrans2D_1, fShowerProfileRecoTrans2D_2, fShowerProfileRecoTrans2D_3, fShowerProfileRecoTrans2D_4, fShowerProfileRecoTrans2D_5, fTransverse, fTransverse_1, fTransverse_2, fTransverse_3, fTransverse_4, fTransverse_5, LBINS, calo::CalorimetryAlg::LifetimeCorrection(), LMAX, LMIN, TBINS, TMAX, TMIN, and X0.

Referenced by analyze().

322  {
323 
324  auto const* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
325  art::ServiceHandle<geo::Geometry> geom;
326 
327  auto collectionPlane = geo::PlaneID(0, 0, 1);
328 
329  double shwVtxTime = detprop->ConvertXToTicks(shwvtx[0], collectionPlane);
330  double shwVtxWire = geom->WireCoordinate(shwvtx[1], shwvtx[2], collectionPlane);
331 
332  double shwTwoTime = detprop->ConvertXToTicks(shwvtx[0]+shwdir[0], collectionPlane);
333  double shwTwoWire = geom->WireCoordinate(shwvtx[1]+shwdir[1], shwvtx[2]+shwdir[2], collectionPlane);
334 
335  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
336  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
337 
338  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
339  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
340  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
341  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
342  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
343 
344  for (size_t i = 0; i < showerhits.size(); ++i) {
345  if (showerhits[i]->WireID().Plane != collectionPlane.Plane) continue;
346 
347  double wirePitch = geom->WirePitch(showerhits[i]->WireID());
348  double tickToDist = detprop->DriftVelocity(detprop->Efield(),detprop->Temperature());
349  tickToDist *= 1.e-3 * detprop->SamplingRate(); // 1e-3 is conversion of 1/us to 1/ns
350 
351  double xvtx = shwVtxTime * tickToDist;
352  double yvtx = shwVtxWire * wirePitch;
353 
354  double xtwo = shwTwoTime * tickToDist;
355  double ytwo = shwTwoWire * wirePitch;
356 
357  double xtwoorth = (ytwo - yvtx) + xvtx;
358  double ytwoorth = -(xtwo - xvtx) + yvtx;
359 
360  double xhit = showerhits[i]->PeakTime() * tickToDist;
361  double yhit = showerhits[i]->WireID().Wire * wirePitch;
362 
363  double ldist = std::abs((ytwoorth-yvtx)*xhit - (xtwoorth-xvtx)*yhit + xtwoorth*yvtx - ytwoorth*xvtx)/std::sqrt( pow((ytwoorth-yvtx), 2) + pow((xtwoorth-xvtx), 2) );
364  double tdist = ((ytwo-yvtx)*xhit - (xtwo-xvtx)*yhit + xtwo*yvtx - ytwo*xvtx)/std::sqrt( pow((ytwo-yvtx), 2) + pow((xtwo-xvtx), 2) );
365 
366  double to3D = 1. / sqrt( pow(xvtx-xtwo,2) + pow(yvtx-ytwo,2) ) ; // distance between two points in 3D space is one
367  ldist *= to3D;
368  tdist *= to3D;
369  double t = ldist/X0;
370 
371  double Q = showerhits[i]->Integral() * fCalorimetryAlg.LifetimeCorrection(showerhits[i]->PeakTime());
372 
373  ltemp->Fill(t, Q);
374  ttemp->Fill(tdist, Q);
375 
376  if (t < 1) ttemp_1->Fill(tdist, Q);
377  else if (t < 2) ttemp_2->Fill(tdist, Q);
378  else if (t < 3) ttemp_3->Fill(tdist, Q);
379  else if (t < 4) ttemp_4->Fill(tdist, Q);
380  else if (t < 5) ttemp_5->Fill(tdist, Q);
381 
382  } // loop through showerhits
383 
384  for (int i = 0; i < LBINS; ++i) {
385  if (ltemp->GetBinContent(i+1) == 0) continue;
386  fShowerProfileRecoLong->Fill(ltemp->GetBinCenter(i+1), ltemp->GetBinContent(i+1));
387  fShowerProfileRecoLong2D->Fill(ltemp->GetBinCenter(i+1), elep, ltemp->GetBinContent(i+1));
388  fLongitudinal->Fill(ltemp->GetBinCenter(i+1), elep, ltemp->GetBinContent(i+1));
389  }
390 
391  for (int i = 0; i < TBINS; ++i) {
392  if (ttemp->GetBinContent(i+1) == 0) continue;
393  fShowerProfileRecoTrans->Fill(ttemp->GetBinCenter(i+1), ttemp->GetBinContent(i+1));
394  fShowerProfileRecoTrans2D->Fill(ttemp->GetBinCenter(i+1), elep, ttemp->GetBinContent(i+1));
395  fTransverse->Fill(ttemp->GetBinCenter(i+1), elep, ttemp->GetBinContent(i+1));
396 
397  fShowerProfileRecoTrans2D_1->Fill(ttemp_1->GetBinCenter(i+1), elep, ttemp_1->GetBinContent(i+1));
398  fShowerProfileRecoTrans2D_2->Fill(ttemp_2->GetBinCenter(i+1), elep, ttemp_2->GetBinContent(i+1));
399  fShowerProfileRecoTrans2D_3->Fill(ttemp_3->GetBinCenter(i+1), elep, ttemp_3->GetBinContent(i+1));
400  fShowerProfileRecoTrans2D_4->Fill(ttemp_4->GetBinCenter(i+1), elep, ttemp_4->GetBinContent(i+1));
401  fShowerProfileRecoTrans2D_5->Fill(ttemp_5->GetBinCenter(i+1), elep, ttemp_5->GetBinContent(i+1));
402 
403  fTransverse_1->Fill(ttemp_1->GetBinCenter(i+1), elep, ttemp_1->GetBinContent(i+1));
404  fTransverse_2->Fill(ttemp_2->GetBinCenter(i+1), elep, ttemp_2->GetBinContent(i+1));
405  fTransverse_3->Fill(ttemp_3->GetBinCenter(i+1), elep, ttemp_3->GetBinContent(i+1));
406  fTransverse_4->Fill(ttemp_4->GetBinCenter(i+1), elep, ttemp_4->GetBinContent(i+1));
407  fTransverse_5->Fill(ttemp_5->GetBinCenter(i+1), elep, ttemp_5->GetBinContent(i+1));
408 
409  }
410 
411  return;
412 
413 } // showerProfile
geo::PlaneID
The data type to uniquely identify a Plane.
Definition: geo_types.h:250
calo::CalorimetryAlg::LifetimeCorrection
double LifetimeCorrection(double time, double T0=0) const
Definition: CalorimetryAlg.cxx:133
void shower::TCShowerTemplateMaker::showerProfileTrue ( std::vector< art::Ptr< recob::Hit > >  allhits,
double  elep 
)

Definition at line 417 of file TCShowerTemplateMaker_module.cc.

References energy, fShowerProfileHitLong, fShowerProfileHitLong2D, fShowerProfileHitTrans, fShowerProfileHitTrans2D, fShowerProfileHitTrans2D_1, fShowerProfileHitTrans2D_2, fShowerProfileHitTrans2D_3, fShowerProfileHitTrans2D_4, fShowerProfileHitTrans2D_5, cheat::BackTrackerService::HitToEveTrackIDEs(), LBINS, LMAX, LMIN, recob::Hit::PeakTime(), TBINS, TMAX, TMIN, cheat::ParticleInventoryService::TrackIdToParticle_P(), geo::WireID::Wire, recob::Hit::WireID(), and X0.

Referenced by analyze().

417  {
418 
419  auto const* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
420  art::ServiceHandle<geo::Geometry> geom;
421  auto collectionPlane = geo::PlaneID(0, 0, 1);
422  art::ServiceHandle<cheat::BackTrackerService> btserv;
423  art::ServiceHandle<cheat::ParticleInventoryService> piserv;
424  std::map<int,double> trkID_E;
425 
426  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
427  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
428 
429  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
430  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
431  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
432  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
433  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
434 
435  double xvtx = -999;
436  double yvtx = -999;
437  double zvtx = -999;
438  double xtwo = -999;
439  double ytwo = -999;
440  double ztwo = -999;
441  double shwvtxT = -999;
442  double shwvtxW = -999;
443  double shwtwoT = -999;
444  double shwtwoW = -999;
445 
446  double shwvtxx = -999;
447  double shwvtxy = -999;
448  double shwtwox = -999;
449  double shwtwoy = -999;
450  double xtwoorth = -999;
451  double ytwoorth = -999;
452 
453  double wirePitch = -999;
454  double tickToDist = -999;
455 
456  bool foundParent = false;
457 
458  for (size_t i = 0; i < allhits.size(); ++i) {
459  if (allhits[i]->WireID().Plane != collectionPlane.Plane) continue;
460 
461  art::Ptr<recob::Hit> hit = allhits[i];
462  std::vector<sim::TrackIDE> trackIDs = btserv->HitToEveTrackIDEs(hit);
463 
464  for (size_t j = 0; j < trackIDs.size(); ++j) {
465  // only want energy associated with the electron and electron must have neutrino mother
466  if ( std::abs((piserv->TrackIdToParticle_P(trackIDs[j].trackID))->PdgCode()) != 11) continue;
467  if ( std::abs((piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Mother()) != 0) continue;
468 
469  if (!foundParent) {
470  xvtx = (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Vx();
471  yvtx = (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Vy();
472  zvtx = (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Vz();
473 
474  xtwo = xvtx + (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Px();
475  ytwo = yvtx + (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Py();
476  ztwo = zvtx + (piserv->TrackIdToParticle_P(trackIDs[j].trackID))->Pz();
477 
478  shwvtxT = detprop->ConvertXToTicks(xvtx, collectionPlane);
479  shwvtxW = geom->WireCoordinate(yvtx, zvtx, collectionPlane);
480 
481  shwtwoT = detprop->ConvertXToTicks(xtwo, collectionPlane);
482  shwtwoW = geom->WireCoordinate(ytwo, ztwo, collectionPlane);
483 
484  wirePitch = geom->WirePitch(hit->WireID());
485  tickToDist = detprop->DriftVelocity(detprop->Efield(),detprop->Temperature());
486  tickToDist *= 1.e-3 * detprop->SamplingRate(); // 1e-3 is conversion of 1/us to 1/ns
487 
488  shwvtxx = shwvtxT * tickToDist;
489  shwvtxy = shwvtxW * wirePitch;
490 
491  shwtwox = shwtwoT * tickToDist;
492  shwtwoy = shwtwoW * wirePitch;
493 
494  xtwoorth = (shwtwoy - shwvtxy) + shwvtxx;
495  ytwoorth = -(shwtwox - shwvtxx) + shwvtxy;
496 
497  foundParent = true;
498  }
499  double xhit = hit->PeakTime() * tickToDist;
500  double yhit = hit->WireID().Wire * wirePitch;
501 
502  double ldist = abs((ytwoorth-shwvtxy)*xhit - (xtwoorth-shwvtxx)*yhit + xtwoorth*shwvtxy - ytwoorth*shwvtxx)/sqrt( pow((ytwoorth-shwvtxy), 2) + pow((xtwoorth-shwvtxx), 2) );
503  double tdist = ((shwtwoy-shwvtxy)*xhit - (shwtwox-shwvtxx)*yhit + shwtwox*shwvtxy - shwtwoy*shwvtxx)/sqrt( pow((shwtwoy-shwvtxy), 2) + pow((shwtwox-shwvtxx), 2) );
504 
505  double to3D = sqrt( pow(xvtx-xtwo , 2) + pow(yvtx-ytwo , 2) + pow(zvtx-ztwo , 2) ) / sqrt( pow(shwvtxx-shwtwox,2) + pow(shwvtxy-shwtwoy,2) ) ; // distance between two points in 3D space is one
506  ldist *= to3D;
507  tdist *= to3D;
508  double t = ldist/X0;
509 
510  double energy = trackIDs[j].energy;
511 
512  ltemp->Fill(t, energy);
513  ttemp->Fill(tdist, energy);
514 
515  if (t < 1) ttemp_1->Fill(tdist, energy);
516  else if (t < 2) ttemp_2->Fill(tdist, energy);
517  else if (t < 3) ttemp_3->Fill(tdist, energy);
518  else if (t < 4) ttemp_4->Fill(tdist, energy);
519  else if (t < 5) ttemp_5->Fill(tdist, energy);
520 
521  } // loop through track IDE
522 
523  } // loop through all hits
524 
525  for (int i = 0; i < LBINS; ++i) {
526  if (ltemp->GetBinContent(i+1) == 0) continue;
527  fShowerProfileHitLong->Fill(ltemp->GetBinCenter(i+1), ltemp->GetBinContent(i+1));
528  fShowerProfileHitLong2D->Fill(ltemp->GetBinCenter(i+1), elep, ltemp->GetBinContent(i+1));
529  }
530 
531  for (int i = 0; i < TBINS; ++i) {
532  if (ttemp->GetBinContent(i+1) == 0) continue;
533  fShowerProfileHitTrans->Fill(ttemp->GetBinCenter(i+1), ttemp->GetBinContent(i+1));
534  fShowerProfileHitTrans2D->Fill(ttemp->GetBinCenter(i+1), elep, ttemp->GetBinContent(i+1));
535 
536  fShowerProfileHitTrans2D_1->Fill(ttemp_1->GetBinCenter(i+1), elep, ttemp_1->GetBinContent(i+1));
537  fShowerProfileHitTrans2D_2->Fill(ttemp_2->GetBinCenter(i+1), elep, ttemp_2->GetBinContent(i+1));
538  fShowerProfileHitTrans2D_3->Fill(ttemp_3->GetBinCenter(i+1), elep, ttemp_3->GetBinContent(i+1));
539  fShowerProfileHitTrans2D_4->Fill(ttemp_4->GetBinCenter(i+1), elep, ttemp_4->GetBinContent(i+1));
540  fShowerProfileHitTrans2D_5->Fill(ttemp_5->GetBinCenter(i+1), elep, ttemp_5->GetBinContent(i+1));
541  }
542 
543  return;
544 
545 } // showerProfileTrue
cheat::ParticleInventoryService::TrackIdToParticle_P
const simb::MCParticle * TrackIdToParticle_P(int const &id)
Definition: ParticleInventoryService_service.cc:134
recob::Hit::WireID
geo::WireID WireID() const
Initial tdc tick for hit.
Definition: Hit.h:234
geo::PlaneID
The data type to uniquely identify a Plane.
Definition: geo_types.h:250
geo::WireID::Wire
WireID_t Wire
Index of the wire within its plane.
Definition: geo_types.h:313
energy
double energy
Definition: plottest35.C:25
hit
Detector simulation of raw signals on wires.
Definition: DumpHits_module.cc:25
recob::Hit::PeakTime
float PeakTime() const
Time of the signal peak, in tick units.
Definition: Hit.h:219
cheat::BackTrackerService::HitToEveTrackIDEs
const std::vector< sim::TrackIDE > HitToEveTrackIDEs(recob::Hit const &hit)
Definition: BackTrackerService_service.cc:166
void shower::TCShowerTemplateMaker::showerProfileTrue ( std::vector< art::Ptr< sim::SimChannel > >  allchan,
simb::MCParticle  electron 
)

Definition at line 549 of file TCShowerTemplateMaker_module.cc.

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

549  {
550  art::ServiceHandle<cheat::ParticleInventoryService> piserv;
551  art::ServiceHandle<geo::Geometry> geom;
552 
553  std::vector<sim::MCEnDep> alledep;
554 
555  TH1F* ltemp = new TH1F("ltemp", "ltemp", LBINS, LMIN, LMAX);
556  TH1F* ttemp = new TH1F("ttemp", "ttemp", TBINS, TMIN, TMAX);
557 
558  TH1F* ttemp_1 = new TH1F("ttemp_1", "ttemp_1", TBINS, TMIN, TMAX);
559  TH1F* ttemp_2 = new TH1F("ttemp_2", "ttemp_2", TBINS, TMIN, TMAX);
560  TH1F* ttemp_3 = new TH1F("ttemp_3", "ttemp_3", TBINS, TMIN, TMAX);
561  TH1F* ttemp_4 = new TH1F("ttemp_4", "ttemp_4", TBINS, TMIN, TMAX);
562  TH1F* ttemp_5 = new TH1F("ttemp_5", "ttemp_5", TBINS, TMIN, TMAX);
563 
564  // get electron energy depositions
565  for (size_t i = 0; i < allchan.size(); ++i) {
566  art::Ptr<sim::SimChannel> simchan = allchan[i];
567  if ( geom->View(simchan->Channel()) != geo::kV ) continue;
568  auto tdc_ide_map = simchan->TDCIDEMap();
569 
570  for(auto const& tdc_ide_pair : tdc_ide_map) {
571  for (auto const& ide : tdc_ide_pair.second) {
572  if (piserv->TrackIdToMotherParticle_P(ide.trackID) == NULL) continue;
573  if ( std::abs(piserv->TrackIdToMotherParticle_P(ide.trackID)->PdgCode()) != 11 ) continue;
574 
575  sim::MCEnDep edep;
576  edep.SetVertex(ide.x,ide.y,ide.z);
577  edep.SetEnergy(ide.energy);
578  edep.SetTrackId(ide.trackID);
579 
580  alledep.push_back(edep);
581 
582  } // loop through ide
583  } // loop through tdc_ide
584 
585  } // loop through channels
586 
587  double x0 = electron.Vx();
588  double y0 = electron.Vy();
589  double z0 = electron.Vz();
590 
591  double x2 = electron.Px();
592  double y2 = electron.Py();
593  double z2 = electron.Pz();
594 
595  TVector3 v0(x2, y2, z2);
596  v0 = v0.Unit();
597 
598  for (size_t i = 0; i < alledep.size(); ++i) {
599  double x = (double)alledep[i].Vertex()[0];
600  double y = (double)alledep[i].Vertex()[1];
601  double z = (double)alledep[i].Vertex()[2];
602 
603  TVector3 v1(x-x0, y-y0, z-z0);
604 
605  double ldist = v0.Dot(v1);
606  double t = ldist/X0;
607  double tdist = (v0.Orthogonal()).Dot(v1);
608 
609  double energy = alledep[i].Energy();
610 
611  ltemp->Fill(t, energy);
612  ttemp->Fill(tdist, energy);
613 
614  if (t < 1) ttemp_1->Fill(tdist, energy);
615  else if (t < 2) ttemp_2->Fill(tdist, energy);
616  else if (t < 3) ttemp_3->Fill(tdist, energy);
617  else if (t < 4) ttemp_4->Fill(tdist, energy);
618  else if (t < 5) ttemp_5->Fill(tdist, energy);
619 
620  }
621 
622  for (int i = 0; i < LBINS; ++i) {
623  if (ltemp->GetBinContent(i+1) == 0) continue;
624  fShowerProfileSimLong->Fill(ltemp->GetBinCenter(i+1), ltemp->GetBinContent(i+1));
625  fShowerProfileSimLong2D->Fill(ltemp->GetBinCenter(i+1), electron.E(), ltemp->GetBinContent(i+1));
626  }
627 
628  for (int i = 0; i < TBINS; ++i) {
629  if (ttemp->GetBinContent(i+1) == 0) continue;
630  fShowerProfileSimTrans->Fill(ttemp->GetBinCenter(i+1), ttemp->GetBinContent(i+1));
631  fShowerProfileSimTrans2D->Fill(ttemp->GetBinCenter(i+1), electron.E(), ttemp->GetBinContent(i+1));
632 
633  fShowerProfileSimTrans2D_1->Fill(ttemp_1->GetBinCenter(i+1), electron.E(), ttemp_1->GetBinContent(i+1));
634  fShowerProfileSimTrans2D_2->Fill(ttemp_2->GetBinCenter(i+1), electron.E(), ttemp_2->GetBinContent(i+1));
635  fShowerProfileSimTrans2D_3->Fill(ttemp_3->GetBinCenter(i+1), electron.E(), ttemp_3->GetBinContent(i+1));
636  fShowerProfileSimTrans2D_4->Fill(ttemp_4->GetBinCenter(i+1), electron.E(), ttemp_4->GetBinContent(i+1));
637  fShowerProfileSimTrans2D_5->Fill(ttemp_5->GetBinCenter(i+1), electron.E(), ttemp_5->GetBinContent(i+1));
638  }
639 
640  return;
641 
642 } // showerProfileTrue
x
Float_t x
Definition: compare.C:6
simb::MCParticle::E
double E(const int i=0) const
Definition: MCParticle.h:237
simb::MCParticle::PdgCode
int PdgCode() const
Definition: MCParticle.h:216
simb::MCParticle::Py
double Py(const int i=0) const
Definition: MCParticle.h:235
sim::MCEnDep::SetEnergy
void SetEnergy(float e)
Definition: MCDataHolder.h:39
geo::kV
Planes which measure V.
Definition: geo_types.h:77
y
Float_t y
Definition: compare.C:6
z
Double_t z
Definition: plot.C:279
simb::MCParticle::Px
double Px(const int i=0) const
Definition: MCParticle.h:234
sim::MCEnDep::SetTrackId
void SetTrackId(unsigned int id)
Definition: MCDataHolder.h:41
y2
Float_t y2[n_points_geant4]
Definition: compare.C:26
energy
double energy
Definition: plottest35.C:25
sim::MCEnDep::SetVertex
void SetVertex(float x, float y, float z)
Definition: MCDataHolder.h:32
edep
Double_t edep
Definition: macro.C:13
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:560
sim::SimChannel::Channel
raw::ChannelID_t Channel() const
Returns the readout channel this object describes.
Definition: SimChannel.h:332
simb::MCParticle::Vx
double Vx(const int i=0) const
Definition: MCParticle.h:225
simb::MCParticle::Pz
double Pz(const int i=0) const
Definition: MCParticle.h:236
simb::MCParticle::Vz
double Vz(const int i=0) const
Definition: MCParticle.h:227
sim::SimChannel::TDCIDEMap
TDCIDEs_t const & TDCIDEMap() const
Returns all the deposited energy information as stored.
Definition: SimChannel.h:331
cheat::ParticleInventoryService::TrackIdToMotherParticle_P
const simb::MCParticle * TrackIdToMotherParticle_P(int const &id)
Definition: ParticleInventoryService_service.cc:141
x2
Float_t x2[n_points_geant4]
Definition: compare.C:26
simb::MCParticle::Vy
double Vy(const int i=0) const
Definition: MCParticle.h:226
void art::Consumer::showMissingConsumes ( ) const
protectedinherited

Definition at line 125 of file Consumer.cc.

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

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

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

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

Definition at line 46 of file EventObserverBase.h.

References art::EventObserverBase::wantAllEvents_.

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

47  {
48  return wantAllEvents_;
49  }
bool art::EventObserverBase::wantEvent ( Event const &  e)
inlineinherited

Definition at line 51 of file EventObserverBase.h.

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

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

52  {
53  return selectors_.wantEvent(e);
54  }
art::EventObserverBase::selectors_
detail::CachedProducts selectors_
Definition: EventObserverBase.h:89
e
Float_t e
Definition: plot.C:34
art::detail::CachedProducts::wantEvent
bool wantEvent(Event const &)
Definition: CachedProducts.cc:52
std::string art::EDAnalyzer::workerType ( ) const
inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

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

Member Data Documentation

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

Definition at line 147 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

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

Definition at line 149 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

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

Definition at line 148 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

calo::CalorimetryAlg shower::TCShowerTemplateMaker::fCalorimetryAlg
protected

Definition at line 158 of file TCShowerTemplateMaker_module.cc.

Referenced by showerProfile().

std::string shower::TCShowerTemplateMaker::fDigitModuleLabel
protected

Definition at line 156 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

std::string shower::TCShowerTemplateMaker::fGenieGenModuleLabel
protected

Definition at line 155 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

std::string shower::TCShowerTemplateMaker::fHitModuleLabel
protected

Definition at line 153 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

TH3F* shower::TCShowerTemplateMaker::fLongitudinal
protected

Definition at line 103 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_electron
protected

Definition at line 112 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_other
protected

Definition at line 126 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fLongitudinal_photon
protected

Definition at line 119 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

std::string shower::TCShowerTemplateMaker::fShowerModuleLabel
protected

Definition at line 154 of file TCShowerTemplateMaker_module.cc.

Referenced by analyze().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileHitLong
protected

Definition at line 68 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitLong2D
protected

Definition at line 72 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileHitTrans
protected

Definition at line 76 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D
protected

Definition at line 80 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_1
protected

Definition at line 84 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_2
protected

Definition at line 88 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_3
protected

Definition at line 92 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_4
protected

Definition at line 96 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileHitTrans2D_5
protected

Definition at line 100 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileRecoLong
protected

Definition at line 69 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoLong2D
protected

Definition at line 73 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans
protected

Definition at line 77 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D
protected

Definition at line 81 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_1
protected

Definition at line 85 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_2
protected

Definition at line 89 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_3
protected

Definition at line 93 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_4
protected

Definition at line 97 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileRecoTrans2D_5
protected

Definition at line 101 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileSimLong
protected

Definition at line 67 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimLong2D
protected

Definition at line 71 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile* shower::TCShowerTemplateMaker::fShowerProfileSimTrans
protected

Definition at line 75 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D
protected

Definition at line 79 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_1
protected

Definition at line 83 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_2
protected

Definition at line 87 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_3
protected

Definition at line 91 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_4
protected

Definition at line 95 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TProfile2D* shower::TCShowerTemplateMaker::fShowerProfileSimTrans2D_5
protected

Definition at line 99 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfileTrue().

TH3F* shower::TCShowerTemplateMaker::fTransverse
protected

Definition at line 104 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TProfile* shower::TCShowerTemplateMaker::fTransverse1_electron
protected

Definition at line 113 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse1_other
protected

Definition at line 127 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse1_photon
protected

Definition at line 120 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse2_electron
protected

Definition at line 114 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse2_other
protected

Definition at line 128 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse2_photon
protected

Definition at line 121 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse3_electron
protected

Definition at line 115 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse3_other
protected

Definition at line 129 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse3_photon
protected

Definition at line 122 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse4_electron
protected

Definition at line 116 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse4_other
protected

Definition at line 130 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse4_photon
protected

Definition at line 123 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse5_electron
protected

Definition at line 117 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse5_other
protected

Definition at line 131 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TProfile* shower::TCShowerTemplateMaker::fTransverse5_photon
protected

Definition at line 124 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob().

TH3F* shower::TCShowerTemplateMaker::fTransverse_1
protected

Definition at line 105 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_2
protected

Definition at line 106 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_3
protected

Definition at line 107 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_4
protected

Definition at line 108 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

TH3F* shower::TCShowerTemplateMaker::fTransverse_5
protected

Definition at line 109 of file TCShowerTemplateMaker_module.cc.

Referenced by beginJob(), and showerProfile().

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

Definition at line 133 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 135 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 134 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 137 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 139 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 138 of file TCShowerTemplateMaker_module.cc.

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

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

Definition at line 151 of file TCShowerTemplateMaker_module.cc.

Referenced by showerProfile(), and showerProfileTrue().

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