Inheritance diagram for lbne::PhotonCounterT0Matching:

Public Types
using	ModuleType = EDProducer

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Member Functions
	PhotonCounterT0Matching (fhicl::ParameterSet const &p)

	PhotonCounterT0Matching (PhotonCounterT0Matching const &)=delete

	PhotonCounterT0Matching (PhotonCounterT0Matching &&)=delete

PhotonCounterT0Matching &	operator= (PhotonCounterT0Matching const &)=delete

PhotonCounterT0Matching &	operator= (PhotonCounterT0Matching &&)=delete

void	produce (art::Event &e) override

void	beginJob () override

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)

void	fillProductDescriptions ()

void	registerProducts (ProductDescriptions &productsToRegister)

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
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	TrackProp (double TrackStart_X, double TrackEnd_X, double &TrackLength_X, double &TrackCentre_X, double TrackStart_Y, double TrackEnd_Y, double &TrackLength_Y, double &TrackCentre_Y, double TrackStart_Z, double TrackEnd_Z, double &TrackLength_Z, double &TrackCentre_Z, double trkTimeStart, double trkTimeEnd, double &trkTimeLengh, double &trkTimeCentre, double &TrackLength)

double	DistFromPoint (double StartY, double EndY, double StartZ, double EndZ, double PointY, double PointZ)

Private Attributes
std::string	fTrackModuleLabel

std::string	fShowerModuleLabel

std::string	fHitsModuleLabel

std::string	fFlashModuleLabel

std::string	fTruthT0ModuleLabel

double	fPredictedXConstant

double	fPredictedXPower = 1

double	fPredictedExpConstant

double	fPredictedExpGradient

double	fDriftWindowSize

double	fWeightOfDeltaYZ

double	fMatchCriteria

double	fPEThreshold

bool	fVerbosity

double	TrackLength_X

double	TrackCentre_X

double	BestTrackCentre_X

double	TrackLength_Y

double	TrackCentre_Y

double	TrackLength_Z

double	TrackCentre_Z

double	trkTimeStart

double	trkTimeEnd

double	trkTimeLengh

double	trkTimeCentre

double	BesttrkTimeCentre

double	TrackLength

double	BestTrackLength

double	PredictedX

double	BestPredictedX

double	TimeSepPredX

double	BestTimeSepPredX

double	DeltaPredX

double	BestDeltaPredX

double	minYZSep

double	BestminYZSep

double	FitParam

double	BestFitParam

double	FlashTime

double	BestFlashTime

double	TimeSep

double	BestTimeSep

int	BestFlash

int	FlashTriggerType = 1

double	YZSep

double	MCTruthT0

TTree *	fTree

TH2D *	hPredX_T

TH2D *	hPredX_PE

TH2D *	hPredX_T_PE

TH2D *	hdeltaX_deltaYZ

TH2D *	hdeltaYZ_Length

TH2D *	hFitParam_Length

TH2D *	hPhotonT0_MCT0

TH1D *	hT0_diff_full

TH1D *	hT0_diff_zoom

Detailed Description

Definition at line 65 of file PhotonCounterT0Matching_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 17 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 26 of file Producer.h.

Constructor & Destructor Documentation

lbne::PhotonCounterT0Matching::PhotonCounterT0Matching ( fhicl::ParameterSet const & p )

explicit

Definition at line 157 of file PhotonCounterT0Matching_module.cc.

References fDriftWindowSize, fFlashModuleLabel, fHitsModuleLabel, fMatchCriteria, fPEThreshold, fPredictedExpConstant, fPredictedExpGradient, fPredictedXConstant, fShowerModuleLabel, fTrackModuleLabel, fTruthT0ModuleLabel, fVerbosity, and fWeightOfDeltaYZ.

                                                                                : EDProducer{p}
 {
   // Call appropriate produces<>() functions here.
   produces<std::vector<anab::T0>>();
   produces<art::Assns<recob::Track, anab::T0>>();
   produces<art::Assns<recob::Shower, anab::T0>>();
 
   fTrackModuleLabel = (p.get<std::string>("TrackModuleLabel"));
   fShowerModuleLabel = (p.get<std::string>("ShowerModuleLabel"));
   fHitsModuleLabel = (p.get<std::string>("HitsModuleLabel"));
   fFlashModuleLabel = (p.get<std::string>("FlashModuleLabel"));
   fTruthT0ModuleLabel = (p.get<std::string>("TruthT0ModuleLabel"));
 
   fPredictedXConstant = (p.get<double>("PredictedXConstant"));
   fPredictedExpConstant = (p.get<double>("PredictedExpConstant"));
   fPredictedExpGradient = (p.get<double>("PredictedExpGradient"));
 
   fDriftWindowSize = (p.get<double>("DriftWindowSize"));
   fWeightOfDeltaYZ = (p.get<double>("WeightOfDeltaYZ"));
   fMatchCriteria = (p.get<double>("MatchCriteria"));
   fPEThreshold = (p.get<double>("PEThreshold"));
 
   fVerbosity = (p.get<bool>("Verbose", false));
 }

lbne::PhotonCounterT0Matching::PhotonCounterT0Matching ( PhotonCounterT0Matching const & )

delete

lbne::PhotonCounterT0Matching::PhotonCounterT0Matching ( PhotonCounterT0Matching && )

delete

Member Function Documentation

void lbne::PhotonCounterT0Matching::beginJob ( )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 182 of file PhotonCounterT0Matching_module.cc.

References BestFitParam, BestFlashTime, BestminYZSep, BestPredictedX, BestTimeSep, BestTimeSepPredX, BestTrackCentre_X, BesttrkTimeCentre, fTree, hdeltaX_deltaYZ, hdeltaYZ_Length, hFitParam_Length, hPhotonT0_MCT0, hPredX_PE, hPredX_T, hPredX_T_PE, hT0_diff_full, hT0_diff_zoom, and MCTruthT0.

 {
   // Implementation of optional member function here.
   art::ServiceHandle<art::TFileService const> tfs;
   fTree = tfs->make<TTree>("PhotonCounterT0Matching", "PhotonCounterT0");
   fTree->Branch("TrackCentre_X", &BestTrackCentre_X, "TrackCentre_X/D");
   fTree->Branch("PredictedX", &BestPredictedX, "PredictedX/D");
   fTree->Branch("TrackTimeCent", &BesttrkTimeCentre, "TimeSepPredX/D");
   fTree->Branch("FlashTime", &BestFlashTime, "FlashTime/D");
   fTree->Branch("TimeSep", &BestTimeSep, "TimeSep/D");
   fTree->Branch("TimeSepPredX", &BestTimeSepPredX, "TimeSepPredX/D");
   fTree->Branch("minYZSep", &BestminYZSep, "minYZSep/D");
   fTree->Branch("FitParam", &BestFitParam, "FitParam/D");
   fTree->Branch("MCTruthT0", &MCTruthT0, "MCTruthT0/D");
 
   hPredX_T = tfs->make<TH2D>("hPredX_T",
                              "Predicted X from timing information against reconstructed X; "
                              "Reconstructed X (cm); Predicted X (cm)",
                              30,
                              0,
                              300,
                              30,
                              0,
                              300);
   hPredX_PE = tfs->make<TH2D>("hPredX_PE",
                               "Predicted X from PE information against reconstructed X; "
                               "Reconstructed X (cm); Predicted X (cm)",
                               30,
                               0,
                               300,
                               30,
                               0,
                               300);
   hPredX_T_PE = tfs->make<TH2D>("hPredX_T_PE",
                                 "Predicted X position from time and PE information; Predicted X "
                                 "from timing information (cm); Predicted X from PE information",
                                 60,
                                 0,
                                 300,
                                 60,
                                 0,
                                 300);
   hdeltaX_deltaYZ = tfs->make<TH2D>(
     "hdeltaX_deltaYZ",
     "Difference between X predicted from PE's and T agaisnt distance of flash from track in YZ; "
     "Difference in X predicted from PE's and T (cm); Distance of flash from track in YZ (cm)",
     40,
     0,
     200,
     40,
     0,
     100);
   hdeltaYZ_Length = tfs->make<TH2D>("hdeltaYZ_Length",
                                     "Distance of flash from track against track length; Distance "
                                     "from flash to track (cm); Track length (cm)",
                                     20,
                                     0,
                                     100,
                                     60,
                                     0,
                                     300);
   hFitParam_Length =
     tfs->make<TH2D>("hFitParam_Length",
                     "How fit correlates with track length; Fit correlation; Track Length (cm)",
                     50,
                     0,
                     250,
                     30,
                     0,
                     300);
   hPhotonT0_MCT0 = tfs->make<TH2D>("hPhotonT0_MCT0",
                                    "Comparing Photon Counter reconstructed T0 against MCTruth T0; "
                                    "Photon Counter T0 (us); MCTruthT0 T0 (us)",
                                    1760,
                                    -1600,
                                    16000,
                                    1760,
                                    -1600,
                                    16000);
   hT0_diff_full = tfs->make<TH1D>(
     "hT0_diff_full",
     "Difference between MCTruth T0 and photon detector T0; Time difference (us); Number",
     500,
     -2500,
     2500);
   hT0_diff_zoom = tfs->make<TH1D>(
     "hT0_diff_zoom",
     "Difference between MCTruth T0 and photon detector T0; Time difference (us); Number",
     320,
     -1.6,
     1.6);
 }

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().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

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().

   {
     collector_.consumesMany<T, BT>();
   }

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().

   {
     return collector_.consumesView<T, BT>(tag);
   }

double lbne::PhotonCounterT0Matching::DistFromPoint	(	double	StartY,
		double	EndY,
		double	StartZ,
		double	EndZ,
		double	PointY,
		double	PointZ
	)

private

Calculate the distance between the centre of the flash and the centre of a line connecting two adjacent space points.

Definition at line 522 of file PhotonCounterT0Matching_module.cc.

References DEFINE_ART_MODULE, and tca::Length().

Referenced by produce().

 {
   double Length = hypot(fabs(EndY - StartY), fabs(EndZ - StartZ));
   double distance =
     ((PointZ - StartZ) * (EndY - StartY) - (PointY - StartY) * (EndZ - StartZ)) / Length;
   return fabs(distance);
 }

void art::detail::Producer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Producer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

   {
     auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

   {
     auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

void art::detail::Producer::doEndJob ( )

inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Producer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto r = rp.makeRun(mc, rp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto sr = srp.makeSubRun(mc, srp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

bool art::detail::Producer::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 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

   {
     auto e = ep.makeEvent(mc);
     ++counts_run;
     ProcessingFrame const frame{mc.scheduleID()};
     produceWithFrame(e, frame);
     e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
     ++counts_passed;
     return true;
   }

void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void art::Modifier::fillProductDescriptions ( )

inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

   {
     ProductRegistryHelper::fillDescriptions(moduleDescription());
   }

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().

   {
     return collector_.getConsumables();
   }

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.

   {
     return doMakeWorker(wp);
   }

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().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

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().

   {
     collector_.mayConsumeMany<T, BT>();
   }

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().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

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().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

PhotonCounterT0Matching& lbne::PhotonCounterT0Matching::operator= ( PhotonCounterT0Matching const & )

delete

PhotonCounterT0Matching& lbne::PhotonCounterT0Matching::operator= ( PhotonCounterT0Matching && )

delete

void lbne::PhotonCounterT0Matching::produce ( art::Event & e )

overridevirtual

Implements art::EDProducer.

Definition at line 275 of file PhotonCounterT0Matching_module.cc.

References BestDeltaPredX, BestFitParam, BestFlash, BestFlashTime, BestminYZSep, BestPredictedX, BestTimeSep, BestTimeSepPredX, BestTrackCentre_X, BestTrackLength, BesttrkTimeCentre, util::CreateAssn(), DeltaPredX, DistFromPoint(), fDriftWindowSize, fFlashModuleLabel, fHitsModuleLabel, art::fill_ptr_vector(), FitParam, FlashTime, FlashTriggerType, fMatchCriteria, fPEThreshold, fPredictedExpConstant, fPredictedExpGradient, fPredictedXConstant, fPredictedXPower, fShowerModuleLabel, fTrackModuleLabel, fTree, fTruthT0ModuleLabel, fVerbosity, fWeightOfDeltaYZ, art::ProductRetriever::getByLabel(), hdeltaX_deltaYZ, hdeltaYZ_Length, hFitParam_Length, hPhotonT0_MCT0, hPredX_PE, hPredX_T, hPredX_T_PE, hT0_diff_full, hT0_diff_zoom, art::Handle< T >::isValid(), MCTruthT0, minYZSep, PredictedX, art::Event::put(), TimeSep, TimeSepPredX, TrackCentre_X, TrackCentre_Y, TrackCentre_Z, TrackLength, TrackLength_X, TrackLength_Y, TrackLength_Z, TrackProp(), trkTimeCentre, trkTimeEnd, trkTimeLengh, trkTimeStart, and YZSep.

 {
   // Access art services...
   art::ServiceHandle<geo::Geometry const> geom;
   auto const clock_data = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
   auto const detprop =
     art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clock_data);
 
   //TrackList handle
   art::Handle<std::vector<recob::Track>> trackListHandle;
   std::vector<art::Ptr<recob::Track>> tracklist;
   if (evt.getByLabel(fTrackModuleLabel, trackListHandle))
     art::fill_ptr_vector(tracklist, trackListHandle);
 
   //ShowerList handle
   art::Handle<std::vector<recob::Shower>> showerListHandle;
   std::vector<art::Ptr<recob::Shower>> showerlist;
   if (evt.getByLabel(fShowerModuleLabel, showerListHandle))
     art::fill_ptr_vector(showerlist, showerListHandle);
 
   //HitList Handle
   art::Handle<std::vector<recob::Hit>> hitListHandle;
   std::vector<art::Ptr<recob::Hit>> hitlist;
   if (evt.getByLabel(fHitsModuleLabel, hitListHandle)) art::fill_ptr_vector(hitlist, hitListHandle);
 
   //FlashList Handle
   art::Handle<std::vector<recob::OpFlash>> flashListHandle;
   std::vector<art::Ptr<recob::OpFlash>> flashlist;
   if (evt.getByLabel(fFlashModuleLabel, flashListHandle))
     art::fill_ptr_vector(flashlist, flashListHandle);
 
   // Create anab::T0 objects and make association with recob::Track
 
   std::unique_ptr<std::vector<anab::T0>> T0col(new std::vector<anab::T0>);
   std::unique_ptr<art::Assns<recob::Track, anab::T0>> Trackassn(
     new art::Assns<recob::Track, anab::T0>);
   std::unique_ptr<art::Assns<recob::Shower, anab::T0>> Showerassn(
     new art::Assns<recob::Shower, anab::T0>);
 
   if (trackListHandle.isValid() && flashListHandle.isValid()) {
     //Access tracks and hits
     art::FindManyP<recob::Hit> fmtht(trackListHandle, evt, fTrackModuleLabel);
     art::FindMany<anab::T0> fmtruth(trackListHandle, evt, fTruthT0ModuleLabel);
 
     size_t NTracks = tracklist.size();
     size_t NFlashes = flashlist.size();
 
     if (fVerbosity)
       std::cout << "There were " << NTracks << " tracks and " << NFlashes
                 << " flashes in this event." << std::endl;
 
     // Now to access PhotonCounter for each track...
     for (size_t iTrk = 0; iTrk < NTracks; ++iTrk) {
       if (fVerbosity) std::cout << "\n New Track " << (int)iTrk << std::endl;
       // Reset Variables.
       BestFlashTime = BestFitParam = BestTrackCentre_X = BestTrackLength = 9999;
       BestTimeSepPredX = BestPredictedX = BestDeltaPredX = BestminYZSep = MCTruthT0 = 9999;
       bool ValidTrack = false;
 
       // Work out Properties of the track.
       recob::Track::Point_t trackStart, trackEnd;
       std::tie(trackStart, trackEnd) = tracklist[iTrk]->Extent();
       std::vector<art::Ptr<recob::Hit>> allHits = fmtht.at(iTrk);
       size_t nHits = allHits.size();
       trkTimeStart = allHits[nHits - 1]->PeakTime() /
                      clock_data.TPCClock().Frequency(); // Got in ticks, now in us!
       trkTimeEnd =
         allHits[0]->PeakTime() / clock_data.TPCClock().Frequency(); // Got in ticks, now in us!
       TrackProp(trackStart.X(),
                 trackEnd.X(),
                 TrackLength_X,
                 TrackCentre_X,
                 trackStart.Y(),
                 trackEnd.Y(),
                 TrackLength_Y,
                 TrackCentre_Y,
                 trackStart.Z(),
                 trackEnd.Z(),
                 TrackLength_Z,
                 TrackCentre_Z,
                 trkTimeStart,
                 trkTimeEnd,
                 trkTimeLengh,
                 trkTimeCentre, // times in us!
                 TrackLength);
 
       // Some cout statement about track properties.
       if (fVerbosity) {
         std::cout << trackStart.X() << " " << trackEnd.X() << " " << TrackLength_X << " "
                   << TrackCentre_X << "\n"
                   << trackStart.Y() << " " << trackEnd.Y() << " " << TrackLength_Y << " "
                   << TrackCentre_Y << "\n"
                   << trackStart.Z() << " " << trackEnd.Z() << " " << TrackLength_Z << " "
                   << TrackCentre_Z << "\n"
                   << trkTimeStart << " " << trkTimeEnd << " " << trkTimeLengh << " "
                   << trkTimeCentre << std::endl;
       }
       // ----- Loop over flashes ------
       for (size_t iFlash = 0; iFlash < NFlashes; ++iFlash) {
         //Reset some flash specific quantities
         YZSep = minYZSep = 9999;
         FlashTime = TimeSep = 9999;
         PredictedX = TimeSepPredX = DeltaPredX = FitParam = 9999;
         // Check flash could be caused by track...
         FlashTime = flashlist[iFlash]->Time(); // Got in us!
         TimeSep = trkTimeCentre - FlashTime;   // Time in us!
         if (TimeSep < 0 || TimeSep > (fDriftWindowSize / clock_data.TPCClock().Frequency()))
           continue; // Times compared in us!
 
         // Check flash has enough PE's to satisfy our threshold
         if (flashlist[iFlash]->TotalPE() < fPEThreshold) continue;
 
         // Work out some quantities for this flash...
         // PredictedX = ( A / x^n ) + exp ( B + Cx )
         PredictedX =
           (fPredictedXConstant / pow(flashlist[iFlash]->TotalPE(), fPredictedXPower)) +
           (exp(fPredictedExpConstant + (fPredictedExpGradient * flashlist[iFlash]->TotalPE())));
         TimeSepPredX = TimeSep * detprop.DriftVelocity(); // us * cm/us = cm!
         DeltaPredX = fabs(TimeSepPredX - PredictedX);
         // Dependant on each point...
         for (size_t Point = 1; Point < tracklist[iTrk]->NumberTrajectoryPoints(); ++Point) {
           auto NewPoint = tracklist[iTrk]->LocationAtPoint(Point);
           auto PrevPoint = tracklist[iTrk]->LocationAtPoint(Point - 1);
           YZSep = DistFromPoint(NewPoint.Y(),
                                 PrevPoint.Y(),
                                 NewPoint.Z(),
                                 PrevPoint.Z(),
                                 flashlist[iFlash]->YCenter(),
                                 flashlist[iFlash]->ZCenter());
           if (Point == 1) minYZSep = YZSep;
           if (YZSep < minYZSep) minYZSep = YZSep;
         }
 
         // Determine how well matched this track is......
         if (fMatchCriteria == 0)
           FitParam =
             pow(((DeltaPredX * DeltaPredX) + (minYZSep * minYZSep * fWeightOfDeltaYZ)), 0.5);
         else if (fMatchCriteria == 1)
           FitParam = minYZSep;
         else if (fMatchCriteria == 2)
           FitParam = DeltaPredX;
 
         //----FLASH INFO-----
         if (fVerbosity) {
           std::cout << "\nFlash " << (int)iFlash << " " << TrackCentre_X << ", " << TimeSepPredX
                     << " - " << PredictedX << " = " << DeltaPredX << ", " << minYZSep << " -> "
                     << FitParam << std::endl;
         }
         //----Select best flash------
         if (FitParam < BestFitParam) {
           ValidTrack = true;
           BestFlash = (int)iFlash;
           BestFitParam = FitParam;
           BestTrackCentre_X = TrackCentre_X;
           BestTrackLength = TrackLength;
           BesttrkTimeCentre = trkTimeCentre;
           BestTimeSepPredX = TimeSepPredX;
           BestPredictedX = PredictedX;
           BestDeltaPredX = DeltaPredX;
           BestminYZSep = minYZSep;
           BestFlashTime = FlashTime;
           BestTimeSep = TimeSep;
         } // Find best Flash
       }   // Loop over Flashes
 
       // ---- Now Make association and fill TTree/Histos with the best matched flash.....
       if (ValidTrack) {
 
         // -- Fill Histos --
         hPredX_T->Fill(BestTrackCentre_X, BestTimeSepPredX);
         hPredX_PE->Fill(BestTrackCentre_X, BestPredictedX);
         hPredX_T_PE->Fill(BestTimeSepPredX, BestPredictedX);
         hdeltaX_deltaYZ->Fill(BestDeltaPredX, BestminYZSep);
         hdeltaYZ_Length->Fill(BestminYZSep, BestTrackLength);
         hFitParam_Length->Fill(BestFitParam, BestTrackLength);
         // ------ Compare Photon Matched to MCTruth Matched -------
         if (fmtruth.isValid()) {
           std::vector<const anab::T0*> T0s = fmtruth.at((int)iTrk);
           for (size_t i = 0; i < T0s.size(); ++i) {
             MCTruthT0 = T0s[i]->Time() / 1e3; // Got in ns, now in us!!
             hPhotonT0_MCT0->Fill(BestFlashTime, MCTruthT0);
             hT0_diff_full->Fill(MCTruthT0 - BestFlashTime);
             hT0_diff_zoom->Fill(MCTruthT0 - BestFlashTime);
           }
         }
         // -- Fill TTree --
         fTree->Fill();
         //Make Association
         T0col->push_back(anab::T0(
           BestFlashTime * 1e3, FlashTriggerType, (int)BestFlash, (*T0col).size(), BestFitParam));
         util::CreateAssn(evt, *T0col, tracklist[iTrk], *Trackassn);
       } // Valid Track
     }   // Loop over tracks
   }
   evt.put(std::move(T0col));
   evt.put(std::move(Trackassn));
   evt.put(std::move(Showerassn));
 
 } // Produce

void art::Modifier::registerProducts ( ProductDescriptions & productsToRegister )

inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

   {
     ProductRegistryHelper::registerProducts(productsToRegister,
                                             moduleDescription());
   }

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

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().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }

void lbne::PhotonCounterT0Matching::TrackProp	(	double	TrackStart_X,
		double	TrackEnd_X,
		double &	TrackLength_X,
		double &	TrackCentre_X,
		double	TrackStart_Y,
		double	TrackEnd_Y,
		double &	TrackLength_Y,
		double &	TrackCentre_Y,
		double	TrackStart_Z,
		double	TrackEnd_Z,
		double &	TrackLength_Z,
		double &	TrackCentre_Z,
		double	trkTimeStart,
		double	trkTimeEnd,
		double &	trkTimeLengh,
		double &	trkTimeCentre,
		double &	TrackLength
	)

private

Calculate central values for track X, Y, Z and time, as well as lengths and overall track length.

Definition at line 475 of file PhotonCounterT0Matching_module.cc.

References TrackLength_X, TrackLength_Y, TrackLength_Z, trkTimeLengh, and trkTimeStart.

Referenced by produce().

 {
   TrackLength_X = fabs(TrackEnd_X - TrackStart_X);
   if (TrackStart_X < TrackEnd_X)
     TrackCentre_X = TrackStart_X + 0.5 * TrackLength_X;
   else
     TrackCentre_X = TrackStart_X - 0.5 * TrackLength_X;
 
   TrackLength_Y = fabs(TrackEnd_Y - TrackStart_Y);
   if (TrackStart_Y < TrackEnd_Y)
     TrackCentre_Y = TrackStart_Y + 0.5 * TrackLength_Y;
   else
     TrackCentre_Y = TrackStart_Y - 0.5 * TrackLength_Y;
 
   TrackLength_Z = fabs(TrackEnd_Z - TrackStart_Z);
   if (TrackStart_Z < TrackEnd_Z)
     TrackCentre_Z = TrackStart_Z + 0.5 * TrackLength_Z;
   else
     TrackCentre_Z = TrackStart_Z - 0.5 * TrackLength_Z;
 
   trkTimeLengh = trkTimeEnd - trkTimeStart;
   trkTimeCentre = trkTimeStart + 0.5 * trkTimeLengh;
 
   TrackLength = pow(pow((TrackEnd_X - TrackStart_X), 2) + pow((TrackEnd_Y - TrackStart_Y), 2) +
                       pow((TrackEnd_Z - TrackStart_Z), 2),
                     0.5);
 
   return;
 }

Member Data Documentation

double lbne::PhotonCounterT0Matching::BestDeltaPredX

private

Definition at line 135 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::BestFitParam

private

Definition at line 137 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

int lbne::PhotonCounterT0Matching::BestFlash

private

Definition at line 140 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::BestFlashTime

private

Definition at line 138 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestminYZSep

private

Definition at line 136 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestPredictedX

private

Definition at line 133 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestTimeSep

private

Definition at line 139 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestTimeSepPredX

private

Definition at line 134 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestTrackCentre_X

private

Definition at line 126 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestTrackLength

private

Definition at line 132 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::BesttrkTimeCentre

private

Definition at line 131 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::DeltaPredX

private

Definition at line 135 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::fDriftWindowSize

private

Definition at line 119 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

std::string lbne::PhotonCounterT0Matching::fFlashModuleLabel

private

Definition at line 113 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

std::string lbne::PhotonCounterT0Matching::fHitsModuleLabel

private

Definition at line 112 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

double lbne::PhotonCounterT0Matching::FitParam

private

Definition at line 137 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::FlashTime

private

Definition at line 138 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

int lbne::PhotonCounterT0Matching::FlashTriggerType = 1

private

Definition at line 141 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::fMatchCriteria

private

Definition at line 121 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

double lbne::PhotonCounterT0Matching::fPEThreshold

private

Definition at line 122 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

double lbne::PhotonCounterT0Matching::fPredictedExpConstant

private

Definition at line 117 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

double lbne::PhotonCounterT0Matching::fPredictedExpGradient

private

Definition at line 118 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

double lbne::PhotonCounterT0Matching::fPredictedXConstant

private

Definition at line 115 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

double lbne::PhotonCounterT0Matching::fPredictedXPower = 1

private

Definition at line 116 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

std::string lbne::PhotonCounterT0Matching::fShowerModuleLabel

private

Definition at line 111 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

std::string lbne::PhotonCounterT0Matching::fTrackModuleLabel

private

Definition at line 110 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

TTree* lbne::PhotonCounterT0Matching::fTree

private

Definition at line 145 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

std::string lbne::PhotonCounterT0Matching::fTruthT0ModuleLabel

private

Definition at line 114 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

bool lbne::PhotonCounterT0Matching::fVerbosity

private

Definition at line 123 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

double lbne::PhotonCounterT0Matching::fWeightOfDeltaYZ

private

Definition at line 120 of file PhotonCounterT0Matching_module.cc.

Referenced by PhotonCounterT0Matching(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hdeltaX_deltaYZ

private

Definition at line 149 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hdeltaYZ_Length

private

Definition at line 150 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hFitParam_Length

private

Definition at line 151 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hPhotonT0_MCT0

private

Definition at line 152 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hPredX_PE

private

Definition at line 147 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hPredX_T

private

Definition at line 146 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hPredX_T_PE

private

Definition at line 148 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH1D* lbne::PhotonCounterT0Matching::hT0_diff_full

private

Definition at line 153 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH1D* lbne::PhotonCounterT0Matching::hT0_diff_zoom

private

Definition at line 154 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::MCTruthT0

private

Definition at line 143 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::minYZSep

private

Definition at line 136 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::PredictedX

private

Definition at line 133 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TimeSep

private

Definition at line 139 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TimeSepPredX

private

Definition at line 134 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackCentre_X

private

Definition at line 126 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackCentre_Y

private

Definition at line 127 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackCentre_Z

private

Definition at line 128 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackLength

private

Definition at line 132 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackLength_X

private

Definition at line 126 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::TrackLength_Y

private

Definition at line 127 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::TrackLength_Z

private

Definition at line 128 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::trkTimeCentre

private

Definition at line 131 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::trkTimeEnd

private

Definition at line 129 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::trkTimeLengh

private

Definition at line 129 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::trkTimeStart

private

Definition at line 129 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::YZSep

private

Definition at line 143 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

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

larana/v09_15_05/source/larana/T0Finder/PhotonCounterT0Matching_module.cc

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation