lbne::PhotonCounterT0Matching Class Reference

Inheritance diagram for lbne::PhotonCounterT0Matching:

Public Types
using	ModuleType = EDProducer
using	WorkerType = WorkerT< EDProducer >
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::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	reconfigure (fhicl::ParameterSet const &p)
template<typename PROD , BranchType B = InEvent>
ProductID	getProductID (std::string const &instanceName={}) const
template<typename PROD , BranchType B>
ProductID	getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const
bool	modifiesEvent () 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
void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
void	prepareForJob (fhicl::ParameterSet const &pset)
void	showMissingConsumes () const

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 86 of file PhotonCounterT0Matching_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 43 of file EDProducer.h.

using art::EDProducer::WorkerType = WorkerT<EDProducer>

inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 164 of file PhotonCounterT0Matching_module.cc.

References reconfigure().

{
  // Call appropriate produces<>() functions here.
  produces< std::vector<anab::T0>               >();
  produces< art::Assns<recob::Track , anab::T0> >();
  produces< art::Assns<recob::Shower, anab::T0> > ();
  reconfigure(p);
}

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 196 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, art::TFileDirectory::make(), and MCTruthT0.

{
  // Implementation of optional member function here.
  art::ServiceHandle<art::TFileService> 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 = 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().

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

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

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

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

{
  return rng()->createEngine(placeholder_schedule_id(), seed);
}

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

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make);
}

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

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
}

CurrentProcessingContext const * art::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

  {
    return current_context_.get();
  }

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 477 of file PhotonCounterT0Matching_module.cc.

References DEFINE_ART_MODULE.

Referenced by produce().

                                                                                                                                         {
  double Length = hypot ( fabs( EndY - StartY), fabs ( EndZ - StartZ ) );
  //  double distance = (double)((point.x - line_start.x) * (line_end.y - line_start.y) - (point.y - line_start.y) * (line_end.x - line_start.x)) / normalLength;
  double distance = ( (PointZ - StartZ) * (EndY - StartY) - (PointY - StartY) * (EndZ - StartZ) ) / Length;
  return fabs(distance);
}

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

{
  auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
  return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
}

template<typename PROD , BranchType B>

ProductID art::EDProducer::getProductID ( std::string const &  instanceName = {} ) const

inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

  {
    return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                               instanceName);
  }

template<typename PROD , BranchType B>

ProductID art::ProducerBase::getProductID ( ModuleDescription const &  moduleDescription, std::string const &  instanceName  ) const

inherited

Definition at line 56 of file ProducerBase.h.

References B, and art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

  {
    auto const& pd =
      get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
    return pd.productID();
  }

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

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

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

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

    {
      return true;
    }

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

{
  static Consumer invalid{InvalidTag{}};
  return &invalid;
}

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

delete

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

delete

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

{
  if (!moduleContext_)
    return;

  pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
  for (auto& consumablesPerBranch : consumables_) {
    cet::sort_all(consumablesPerBranch);
  }
}

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

overridevirtual

Implements art::EDProducer.

Definition at line 228 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, Fill(), art::fill_ptr_vector(), FitParam, FlashTime, FlashTriggerType, fMatchCriteria, fPEThreshold, fPredictedExpConstant, fPredictedExpGradient, fPredictedXConstant, fPredictedXPower, fShowerModuleLabel, fTrackModuleLabel, fTree, fTruthT0ModuleLabel, fVerbosity, fWeightOfDeltaYZ, art::DataViewImpl::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> geom;
  auto const* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
  auto const* timeservice = lar::providerFrom<detinfo::DetectorClocksService>();
  //Commenting this out since it doesn't appear to actually be used.
//  art::ServiceHandle<cheat::BackTrackerService> bt_serv;

  //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() / timeservice->TPCClock().Frequency(); //Got in ticks, now in us!
      trkTimeEnd   = allHits[0]->PeakTime() / timeservice->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/timeservice->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------
        //double YFitRegion = (-1 * DeltaPredX ) + 80;
        //if ( minYZSep > YFitRegion ) continue;
        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 );
            //std::cout << "Size " << T0s.size() << " " << MCTruthT0 << " " << BestFlashTime << std::endl;      
          }
        }
        // -- Fill TTree --
        fTree->Fill();
        //Make Association
        T0col->push_back(anab::T0(BestFlashTime * 1e3,
                                  FlashTriggerType,
                                  (int)BestFlash,
                                  (*T0col).size(),
                                  BestFitParam
                                  ));
        util::CreateAssn(*this, evt, *T0col, tracklist[iTrk], *Trackassn);
      } // Valid Track
    } // Loop over tracks   
  }
  /* // ------------------------------------------------- SHOWER STUFF -------------------------------------------------
  if (showerListHandle.isValid()){
    art::FindManyP<recob::Hit> fmsht(showerListHandle,evt, fShowerModuleLabel);
    // Now Loop over showers....
    size_t NShowers = showerlist.size();
    for (size_t Shower = 0; Shower < NShowers; ++Shower) {
      ShowerMatchID     = 0;
      ShowerID          = 0;
      ShowerT0          = 0;
      std::vector< art::Ptr<recob::Hit> > allHits = fmsht.at(Shower);
      
      std::map<int,double> showeride;
      for(size_t h = 0; h < allHits.size(); ++h){
        art::Ptr<recob::Hit> hit = allHits[h];
        std::vector<sim::IDE> ides;
        std::vector<sim::TrackIDE> TrackIDs = bt->HitToTrackID(hit);
        
        for(size_t e = 0; e < TrackIDs.size(); ++e){
          showeride[TrackIDs[e].trackID] += TrackIDs[e].energy;
        }
      }
      // Work out which IDE despoited the most charge in the hit if there was more than one.
      double maxe = -1;
      double tote = 0;
      for (std::map<int,double>::iterator ii = showeride.begin(); ii!=showeride.end(); ++ii){
        tote += ii->second;
        if ((ii->second)>maxe){
          maxe = ii->second;
          ShowerID = ii->first;
        }
      }
      // Now have MCParticle trackID corresponding to shower, so get PdG code and T0 etc.
      const simb::MCParticle *particle = bt->TrackIDToParticle(ShowerID);
      ShowerT0 = particle->T();
      ShowerID = particle->TrackId();
      ShowerTriggerType = 1; // Using PhotonCounter as trigger, so tigger type is 1.
      
      T0col->push_back(anab::T0(ShowerT0,
                                ShowerTriggerType,
                                ShowerID,
                                (*T0col).size()
                                ));
      util::CreateAssn(*this, evt, *T0col, showerlist[Shower], *Showerassn);    
    }// Loop over showers
  }
  */ 
  evt.put(std::move(T0col));
  evt.put(std::move(Trackassn));
  evt.put(std::move(Showerassn));
  
} // Produce

void lbne::PhotonCounterT0Matching::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 175 of file PhotonCounterT0Matching_module.cc.

References fDriftWindowSize, fFlashModuleLabel, fHitsModuleLabel, fMatchCriteria, fPEThreshold, fPredictedExpConstant, fPredictedExpGradient, fPredictedXConstant, fShowerModuleLabel, fTrackModuleLabel, fTruthT0ModuleLabel, fVerbosity, fWeightOfDeltaYZ, and fhicl::ParameterSet::get().

Referenced by PhotonCounterT0Matching().

{
  // Implementation of optional member function here.
  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) );
}

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

{
  if (!moduleContext_)
    return;

  // If none of the branches have missing consumes statements, exit early.
  if (std::all_of(cbegin(missingConsumes_),
                  cend(missingConsumes_),
                  [](auto const& perBranch) { return perBranch.empty(); }))
    return;

  constexpr cet::HorizontalRule rule{60};
  mf::LogPrint log{"MTdiagnostics"};
  log << '\n'
      << rule('=') << '\n'
      << "The following consumes (or mayConsume) statements are missing from\n"
      << module_context(moduleDescription_) << '\n'
      << rule('-') << '\n';

  cet::for_all_with_index(
    missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
      for (auto const& pi : perBranch) {
        log << "  "
            << assemble_consumes_statement(static_cast<BranchType>(i), pi)
            << '\n';
      }
    });
  log << rule('=');
}

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 451 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;
}

void art::Consumer::validateConsumedProduct ( BranchType const  bt, ProductInfo const &  pi  )

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

{
  // Early exits if consumes tracking has been disabled or if the
  // consumed product is an allowed consumable.
  if (!moduleContext_)
    return;

  if (cet::binary_search_all(consumables_[bt], pi))
    return;

  if (requireConsumes_) {
    throw Exception(errors::ProductRegistrationFailure,
                    "Consumer: an error occurred during validation of a "
                    "retrieved product\n\n")
      << "The following consumes (or mayConsume) statement is missing from\n"
      << module_context(moduleDescription_) << ":\n\n"
      << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
  }

  missingConsumes_[bt].insert(pi);
}

Member Data Documentation

double lbne::PhotonCounterT0Matching::BestDeltaPredX

private

Definition at line 141 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::BestFitParam

private

Definition at line 143 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

int lbne::PhotonCounterT0Matching::BestFlash

private

Definition at line 146 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::BestFlashTime

private

Definition at line 144 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestminYZSep

private

Definition at line 142 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestPredictedX

private

Definition at line 139 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestTimeSep

private

Definition at line 145 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestTimeSepPredX

private

Definition at line 140 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestTrackCentre_X

private

Definition at line 132 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::BestTrackLength

private

Definition at line 138 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::BesttrkTimeCentre

private

Definition at line 137 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::DeltaPredX

private

Definition at line 141 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::fDriftWindowSize

private

Definition at line 125 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

std::string lbne::PhotonCounterT0Matching::fFlashModuleLabel

private

Definition at line 119 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

std::string lbne::PhotonCounterT0Matching::fHitsModuleLabel

private

Definition at line 118 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

double lbne::PhotonCounterT0Matching::FitParam

private

Definition at line 143 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::FlashTime

private

Definition at line 144 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

int lbne::PhotonCounterT0Matching::FlashTriggerType =1

private

Definition at line 147 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::fMatchCriteria

private

Definition at line 127 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

double lbne::PhotonCounterT0Matching::fPEThreshold

private

Definition at line 128 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

double lbne::PhotonCounterT0Matching::fPredictedExpConstant

private

Definition at line 123 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

double lbne::PhotonCounterT0Matching::fPredictedExpGradient

private

Definition at line 124 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

double lbne::PhotonCounterT0Matching::fPredictedXConstant

private

Definition at line 121 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

double lbne::PhotonCounterT0Matching::fPredictedXPower = 1

private

Definition at line 122 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

std::string lbne::PhotonCounterT0Matching::fShowerModuleLabel

private

Definition at line 117 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

std::string lbne::PhotonCounterT0Matching::fTrackModuleLabel

private

Definition at line 116 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

TTree* lbne::PhotonCounterT0Matching::fTree

private

Definition at line 151 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

std::string lbne::PhotonCounterT0Matching::fTruthT0ModuleLabel

private

Definition at line 120 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

bool lbne::PhotonCounterT0Matching::fVerbosity

private

Definition at line 129 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

double lbne::PhotonCounterT0Matching::fWeightOfDeltaYZ

private

Definition at line 126 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and reconfigure().

TH2D* lbne::PhotonCounterT0Matching::hdeltaX_deltaYZ

private

Definition at line 155 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hdeltaYZ_Length

private

Definition at line 156 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hFitParam_Length

private

Definition at line 157 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hPhotonT0_MCT0

private

Definition at line 158 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hPredX_PE

private

Definition at line 153 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hPredX_T

private

Definition at line 152 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH2D* lbne::PhotonCounterT0Matching::hPredX_T_PE

private

Definition at line 154 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH1D* lbne::PhotonCounterT0Matching::hT0_diff_full

private

Definition at line 159 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

TH1D* lbne::PhotonCounterT0Matching::hT0_diff_zoom

private

Definition at line 160 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::MCTruthT0

private

Definition at line 149 of file PhotonCounterT0Matching_module.cc.

Referenced by beginJob(), and produce().

double lbne::PhotonCounterT0Matching::minYZSep

private

Definition at line 142 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::PredictedX

private

Definition at line 139 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TimeSep

private

Definition at line 145 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TimeSepPredX

private

Definition at line 140 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackCentre_X

private

Definition at line 132 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackCentre_Y

private

Definition at line 133 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackCentre_Z

private

Definition at line 134 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackLength

private

Definition at line 138 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::TrackLength_X

private

Definition at line 132 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::TrackLength_Y

private

Definition at line 133 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::TrackLength_Z

private

Definition at line 134 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::trkTimeCentre

private

Definition at line 137 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::trkTimeEnd

private

Definition at line 135 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

double lbne::PhotonCounterT0Matching::trkTimeLengh

private

Definition at line 135 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::trkTimeStart

private

Definition at line 135 of file PhotonCounterT0Matching_module.cc.

Referenced by produce(), and TrackProp().

double lbne::PhotonCounterT0Matching::YZSep

private

Definition at line 149 of file PhotonCounterT0Matching_module.cc.

Referenced by produce().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larana/v07_04_02/source/larana/T0Finder/PhotonCounterT0Matching_module.cc