opreco::OpticalRecoAna Class Reference

#include "OpticalRecoAna.h"

Inheritance diagram for opreco::OpticalRecoAna:

Public Types
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer

Public Member Functions
	OpticalRecoAna (const fhicl::ParameterSet &)
virtual	~OpticalRecoAna ()
void	beginJob ()
void	analyze (const art::Event &)
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

Private Member Functions
void	get_MC_particle_list (sim::ParticleList const &, std::vector< simb::MCParticle > &)
simb::Origin_t	get_MC_particle_origin (simb::MCParticle const &)
float	update_MC_particle_time (simb::MCParticle const &, bool &, geo::Geometry const &)
void	match_flashes_to_tracks (std::vector< recob::OpFlash > const &, std::vector< recob::Track > const &)
void	compare_track_and_flash (recob::Track const &, recob::OpFlash const &, bool &)
void	match_flashes_to_particles (std::vector< recob::OpFlash > const &, std::vector< simb::MCParticle > const &, float const &, geo::Geometry const &)
void	compare_particle_and_flash (simb::MCParticle const &, recob::OpFlash const &, bool &, float const &, geo::Geometry const &)
void	match_flashes_to_particles (std::vector< recob::Track > const &, std::vector< simb::MCParticle > const &, geo::Geometry const &)
void	compare_particle_and_track (simb::MCParticle const &, recob::Track const &, bool &, geo::Geometry const &)
void	FillMatchTree_PF (std::vector< recob::OpFlash > const &, std::vector< simb::MCParticle > const &, float const &, geo::Geometry const &)

Private Attributes
std::string	fFlashModuleLabel
std::string	fTrackModuleLabel
float	fKineticEnergyMin
float	fPEMin
float	fTimeMatchMax
std::vector< flash_match >	fFlash_match_vector
std::vector< track_match >	fTrack_match_vector
std::vector< particle_match >	fParticle_match_vector
TH1F *	fTimeDiff
TH1F *	fTimeDiff_fine
TH1I *	fNBeamFlashes
TTree *	fMatchTree_PF
TTree *	fMatchTree_PF_NotNu
int	fRun
int	fEvent
int	fParticleID
int	fParticleMother
int	fParticleTrackID
float	fParticleTime
float	fParticleVx
float	fParticleVy
float	fParticleVz
int	fFlashID
float	fFlashTime
float	fFlashTimeWidth
float	fFlashY
float	fFlashZ
float	fFlashYWidth
float	fFlashZWidth
float	fFlashPE
int	fFlashOnBeamTime
int	fMatchIndex
int	fMatchIndex_NotNu

Detailed Description

Definition at line 42 of file OpticalRecoAna.h.

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

opreco::OpticalRecoAna::OpticalRecoAna

(

const fhicl::ParameterSet &

)

Definition at line 28 of file OpticalRecoAna_module.cc.

References fFlash_match_vector, fFlashModuleLabel, fKineticEnergyMin, fParticle_match_vector, fPEMin, fTimeMatchMax, fTrack_match_vector, fTrackModuleLabel, and fhicl::ParameterSet::get().

                                                                 : EDAnalyzer(pset)
{
  
  // Indicate that the Input Module comes from .fcl
  fFlashModuleLabel  = pset.get<std::string>("FlashModule");
  fTrackModuleLabel  = pset.get<std::string>("TrackModule");
  fKineticEnergyMin  = pset.get<float>("KineticEnergyMin");
  fPEMin             = pset.get<float>("PEMin");
  fTimeMatchMax      = pset.get<float>("TimeMatchMax");
  
  fFlash_match_vector.clear();
  fTrack_match_vector.clear();
  fParticle_match_vector.clear();
}

opreco::OpticalRecoAna::~OpticalRecoAna ( )

virtual

Definition at line 46 of file OpticalRecoAna_module.cc.

{}

Member Function Documentation

void opreco::OpticalRecoAna::analyze

(

const art::Event &

evt

)

Definition at line 106 of file OpticalRecoAna_module.cc.

References art::EventID::event(), fEvent, fFlash_match_vector, fFlashModuleLabel, fNBeamFlashes, fParticle_match_vector, fRun, get_MC_particle_list(), art::DataViewImpl::getByLabel(), art::Event::id(), art::Event::isRealData(), LOG_INFO, match_flashes_to_particles(), cheat::ParticleInventoryService::ParticleList(), and art::EventID::run().

{
  
  //get flag for MC
  const bool is_MC= !(evt.isRealData());

  fRun = evt.id().run();
  fEvent = evt.id().event();
  
  //first get out track, flash, and particle list handles
  art::Handle< std::vector<recob::OpFlash> > flash_handle;
  evt.getByLabel(fFlashModuleLabel, flash_handle);
  std::vector<recob::OpFlash> const& flash_vector(*flash_handle);
  fFlash_match_vector.resize(flash_vector.size());
  
  LOG_INFO ("OpticalRecoAna")  
    << "Number of flashes is " << flash_vector.size() << std::flush;
  
  //art::Handle< std::vector<recob::Track> > track_handle;
  //evt.getByLabel(fTrackModuleLabel, track_handle);
  //std::vector<recob::Track> const& track_vector(*track_handle);
  //fTrack_match_vector.resize(track_vector.size());
  
  //LOG_INFO ("OpticalRecoAna")  
  //<< "Number of tracks is " << track_vector.size() << std::flush;
  
  //match_flashes_to_tracks(flash_vector, track_vector);

  //all this for the MC matching
  if(is_MC){
    
    art::ServiceHandle<cheat::ParticleInventoryService> pi_serv;
    std::vector<simb::MCParticle> particle_list;
    get_MC_particle_list(pi_serv->ParticleList(),particle_list);
    fParticle_match_vector.resize(particle_list.size());
    
    mf::LogInfo("OpticalRecoAna")  
      << "Number of MC particles is " << particle_list.size();
    

    art::ServiceHandle<opdet::OpDigiProperties> odp;
    const float ns_per_PMT_tick = 1e3;// already converted to microseconds//( 1e3 / odp->SampleFreq()) ; //SampleFreq is in MHz
    art::ServiceHandle<geo::Geometry> geometry_handle;

    match_flashes_to_particles(flash_vector,
                               particle_list,
                               ns_per_PMT_tick,
                               *geometry_handle);
    /*
    FillMatchTree_PF(flash_vector,
                     particle_list,
                     ns_per_PMT_tick,
                     *geometry_handle);
                     */
    int n_flashes = 0;
    for(auto const& my_flash : flash_vector){
      if(my_flash.OnBeamTime()==1) n_flashes++;
    }
    fNBeamFlashes->Fill(n_flashes);

    //check_flash_matches();

    //match_tracks_to_particles(track_handle,particle_list);
    
  }
  
}

void opreco::OpticalRecoAna::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 52 of file OpticalRecoAna_module.cc.

References fEvent, fFlashID, fFlashOnBeamTime, fFlashPE, fFlashTime, fFlashTimeWidth, fFlashY, fFlashYWidth, fFlashZ, fMatchIndex, fMatchIndex_NotNu, fMatchTree_PF, fMatchTree_PF_NotNu, fNBeamFlashes, fParticleID, fParticleMother, fParticleTime, fParticleTrackID, fParticleVx, fParticleVy, fParticleVz, fRun, fTimeDiff, fTimeDiff_fine, and art::TFileDirectory::make().

{
  art::ServiceHandle<art::TFileService> tfs;
  fTimeDiff = tfs->make<TH1F>("htdiff","Time difference between particles and flashes; t_diff (ns); flash/particle pairs",1e3,-5e6,5e6);
  fTimeDiff_fine = tfs->make<TH1F>("htdiff_fine","Time difference between particles and flashes; t_diff (ns); flash/particle pairs",100,-1000,1000);
  fNBeamFlashes = tfs->make<TH1I>("hNBeamFlashes","Number of flashes OnBeamTime per event; N_{Flashes}; Events",5,0,5);

  fMatchTree_PF = tfs->make<TTree>("MatchTree_PF","MatchTree_PF");
  fMatchTree_PF->Branch("Run",&fRun,"Run/I");
  fMatchTree_PF->Branch("Event",&fEvent,"Event/I");
  fMatchTree_PF->Branch("p_id",&fParticleID,"p_id/I");
  fMatchTree_PF->Branch("p_time",&fParticleTime,"p_time/F");
  fMatchTree_PF->Branch("p_vx",&fParticleVx,"p_vx/F");
  fMatchTree_PF->Branch("p_vy",&fParticleVy,"p_vy/F");
  fMatchTree_PF->Branch("p_vz",&fParticleVz,"p_vz/F");
  fMatchTree_PF->Branch("p_mid",&fParticleMother,"p_mid/I");
  fMatchTree_PF->Branch("p_trackid",&fParticleTrackID,"p_trackid/I");
  fMatchTree_PF->Branch("FlashID",&fFlashID,"flash_id/I");
  fMatchTree_PF->Branch("f_time",&fFlashTime,"f_time/F");
  fMatchTree_PF->Branch("f_timewidth",&fFlashTimeWidth,"f_timewidth/F");
  fMatchTree_PF->Branch("f_y",&fFlashY,"f_y/F");
  fMatchTree_PF->Branch("f_z",&fFlashZ,"f_z/F");
  fMatchTree_PF->Branch("f_ywidth",&fFlashYWidth,"f_ywidth/F");
  fMatchTree_PF->Branch("f_zwidth",&fFlashYWidth,"f_zwidth/F");
  fMatchTree_PF->Branch("f_onbeamtime",&fFlashOnBeamTime,"f_onbeamtime/I");
  fMatchTree_PF->Branch("f_pe",&fFlashPE,"f_pe/F");
  fMatchTree_PF->Branch("matchIndex",&fMatchIndex,"matchIndex/I");
  
  fMatchTree_PF_NotNu = tfs->make<TTree>("MatchTree_PF_NotNu","MatchTree_PF_NotNu");
  fMatchTree_PF_NotNu->Branch("Run",&fRun,"Run/I");
  fMatchTree_PF_NotNu->Branch("Event",&fEvent,"Event/I");
  fMatchTree_PF_NotNu->Branch("p_id",&fParticleID,"p_id/I");
  fMatchTree_PF_NotNu->Branch("p_time",&fParticleTime,"p_time/F");
  fMatchTree_PF_NotNu->Branch("p_vx",&fParticleVx,"p_vx/F");
  fMatchTree_PF_NotNu->Branch("p_vy",&fParticleVy,"p_vy/F");
  fMatchTree_PF_NotNu->Branch("p_vz",&fParticleVz,"p_vz/F");
  fMatchTree_PF_NotNu->Branch("p_mid",&fParticleMother,"p_mid/I");
  fMatchTree_PF_NotNu->Branch("p_trackid",&fParticleTrackID,"p_trackid/I");
  fMatchTree_PF_NotNu->Branch("FlashID",&fFlashID,"flash_id/I");
  fMatchTree_PF_NotNu->Branch("f_time",&fFlashTime,"f_time/F");
  fMatchTree_PF_NotNu->Branch("f_timewidth",&fFlashTimeWidth,"f_timewidth/F");
  fMatchTree_PF_NotNu->Branch("f_y",&fFlashY,"f_y/F");
  fMatchTree_PF_NotNu->Branch("f_z",&fFlashZ,"f_z/F");
  fMatchTree_PF_NotNu->Branch("f_ywidth",&fFlashYWidth,"f_ywidth/F");
  fMatchTree_PF_NotNu->Branch("f_zwidth",&fFlashYWidth,"f_zwidth/F");
  fMatchTree_PF_NotNu->Branch("f_onbeamtime",&fFlashOnBeamTime,"f_onbeamtime/I");
  fMatchTree_PF_NotNu->Branch("f_pe",&fFlashPE,"f_pe/F");
  fMatchTree_PF_NotNu->Branch("matchIndex",&fMatchIndex_NotNu,"matchIndex/I");

}

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

  {
    return selectors_;
  }

void opreco::OpticalRecoAna::compare_particle_and_flash ( simb::MCParticle const &  particle, recob::OpFlash const &  flash, bool &  matching, float const &  ns_per_PMT_tick, geo::Geometry const &  geometry  )

private

Definition at line 399 of file OpticalRecoAna_module.cc.

                                                                                    {
} 

void opreco::OpticalRecoAna::compare_particle_and_track ( simb::MCParticle const &  particle, recob::Track const &  track, bool &  matching, geo::Geometry const &  geometry  )

private

Definition at line 432 of file OpticalRecoAna_module.cc.

Referenced by match_flashes_to_particles().

                                                                                    {
}

void opreco::OpticalRecoAna::compare_track_and_flash ( recob::Track const &  track, recob::OpFlash const &  flash, bool &  matching  )

private

Definition at line 271 of file OpticalRecoAna_module.cc.

Referenced by match_flashes_to_tracks().

                                                                     {
}

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::EDAnalyzer::currentContext ( ) const

protectedinherited

Definition at line 114 of file EDAnalyzer.cc.

References art::EDAnalyzer::current_context_.

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

  {
    return current_context_.get();
  }

void opreco::OpticalRecoAna::FillMatchTree_PF ( std::vector< recob::OpFlash > const &  flash_vector, std::vector< simb::MCParticle > const &  particle_vector, float const &  ns_per_PMT_tick, geo::Geometry const &  geometry  )

private

Definition at line 357 of file OpticalRecoAna_module.cc.

References fFlash_match_vector, fFlashID, fFlashTime, fFlashY, fFlashYWidth, fFlashZ, fFlashZWidth, fMatchTree_PF, fParticleID, fParticleTime, fParticleVx, fParticleVy, fParticleVz, recob::OpFlash::Time(), and update_MC_particle_time().

                                                                          {

    for(size_t i_flash=0; i_flash < flash_vector.size(); i_flash++){

      recob::OpFlash const& my_flash( flash_vector.at(i_flash) );
      const float flash_time = my_flash.Time()*ns_per_PMT_tick;

      for(auto i_particle : fFlash_match_vector.at(i_flash).particle_indices){

        simb::MCParticle const& my_particle( particle_vector.at(i_particle) );
        bool pass_check = false;
        const float particle_time = update_MC_particle_time(my_particle,pass_check,geometry);
        
        if(my_particle.Mother()==0 && my_particle.TrackId()>1e4){

          fParticleID = i_particle;
          fParticleTime = particle_time;
          fParticleVx = my_particle.Vx(0);
          fParticleVy = my_particle.Vy(0);
          fParticleVz = my_particle.Vz(0);
          fFlashID = i_flash;
          fFlashTime = flash_time;
          fFlashY = my_flash.YCenter();
          fFlashZ = my_flash.ZCenter();
          fFlashYWidth = my_flash.YWidth();
          fFlashZWidth = my_flash.ZWidth();

          fMatchTree_PF->Fill();

          return;
        }

      }

  }

}

void opreco::OpticalRecoAna::get_MC_particle_list ( sim::ParticleList const &  plist, std::vector< simb::MCParticle > &  particle_vector  )

private

Definition at line 179 of file OpticalRecoAna_module.cc.

References sim::ParticleList::begin(), simb::MCParticle::E(), sim::ParticleList::end(), fKineticEnergyMin, ipart, simb::MCParticle::Mass(), and simb::MCParticle::PdgCode().

Referenced by analyze().

                                                                                                                        {
  
  for(sim::ParticleList::const_iterator ipart = plist.begin(); ipart != plist.end(); ++ipart) {
    
    const simb::MCParticle* particle = (*ipart).second;
    /*
    std::cout << "Particle " << particle_vector.size() << " (" << (*ipart).first << ") is " << particle->PdgCode() << " (mother=" << particle->Mother() 
                << ") with K_energy " << particle->E()-0.001*particle->Mass() << std::endl;
      std::cout << "\t(X,Y,Z,T)=(" << particle->Vx() << "," << particle->Vy() << "," << particle->Vz() << "," << particle->T() << ")" << std::endl;
    */
    //do not store if it's below our energy cut
    if( particle->E() < 0.001*particle->Mass() +  fKineticEnergyMin ) continue;
    
    //check to see if it's a charged particle we expect to leave ionization
    const int pdg = std::abs(particle->PdgCode());
    if( pdg==11 // electron
        || pdg==13 //muon
        || pdg==15 //tau
        || pdg==211 //charged pion
        || pdg==321 //charged kaon
        || pdg==2212 //proton
        || pdg==2214 //delta
        || pdg==2224 //delta
        || pdg==3222 //sigma
        || pdg==3112 //sigma
        || pdg==3312 //xi
        || pdg==3334 //omega
        ) {
      particle_vector.emplace_back(*particle);
      //std::cout << "Particle " << particle_vector.size() << " (" << (*ipart).first << ") is " << pdg << " (status=" << particle->StatusCode() 
      //        << ") with K_energy " << particle->E()-0.001*particle->Mass() << std::endl;
      //std::cout << "\t(X,Y,Z,T)=(" << particle->Vx() << "," << particle->Vy() << "," << particle->Vz() << "," << particle->T() << ")" << std::endl;
    }
  }
  
}

simb::Origin_t opreco::OpticalRecoAna::get_MC_particle_origin ( simb::MCParticle const &  particle )

private

Definition at line 174 of file OpticalRecoAna_module.cc.

References simb::MCParticle::TrackId(), and cheat::ParticleInventoryService::TrackIdToMCTruth_P().

Referenced by match_flashes_to_particles().

                                                                                       {
  art::ServiceHandle<cheat::ParticleInventoryService> pi_serv;
  return (pi_serv->TrackIdToMCTruth_P(particle.TrackId()))->Origin();
}

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();
}

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

  {
    return selectors_.getOneTriggerResults(e);
  }

void opreco::OpticalRecoAna::match_flashes_to_particles ( std::vector< recob::OpFlash > const &  flash_vector, std::vector< simb::MCParticle > const &  particle_vector, float const &  ns_per_PMT_tick, geo::Geometry const &  geometry  )

private

Definition at line 276 of file OpticalRecoAna_module.cc.

References fFlash_match_vector, fFlashID, fFlashOnBeamTime, fFlashPE, fFlashTime, fFlashTimeWidth, fFlashY, fFlashYWidth, fFlashZ, fFlashZWidth, fMatchIndex, fMatchIndex_NotNu, fMatchTree_PF, fMatchTree_PF_NotNu, fParticle_match_vector, fParticleID, fParticleMother, fParticleTime, fParticleTrackID, fParticleVx, fParticleVy, fParticleVz, fPEMin, fTimeDiff, fTimeDiff_fine, fTimeMatchMax, get_MC_particle_origin(), simb::kBeamNeutrino, recob::OpFlash::Time(), and update_MC_particle_time().

Referenced by analyze().

                                                                                    {

  fMatchIndex=0;
  fMatchIndex_NotNu=0;
  int lastFlashID = -1;

  for(size_t i_particle=0; i_particle < particle_vector.size(); i_particle++){
    
    simb::MCParticle const& my_particle( particle_vector.at(i_particle) );
    bool pass_check = false;
    const float particle_time = update_MC_particle_time(my_particle,pass_check,geometry);
    if(!pass_check) continue;
    
    //std::cout << "Particle " << i_particle << " (" << my_particle.PdgCode() << ")" << particle_time << std::endl;
    
    for(size_t i_flash=0; i_flash < flash_vector.size(); i_flash++){
      
      recob::OpFlash const& my_flash( flash_vector.at(i_flash) );
      if(my_flash.TotalPE() < fPEMin) continue;
      
      const float flash_time = my_flash.Time()*ns_per_PMT_tick;

      //if(my_flash.OnBeamTime()) 
        //std::cout << "\tFlash " << i_flash << " time is " << flash_time << std::endl;

      fTimeDiff->Fill(particle_time-flash_time);
      fTimeDiff_fine->Fill(particle_time-flash_time);
      
      fParticleID = i_particle;
      fParticleTime = particle_time;
      fParticleVx = my_particle.Vx(0);
      fParticleVy = my_particle.Vy(0);
      fParticleVz = my_particle.Vz(0);
      fParticleMother = my_particle.Mother();
      fParticleTrackID = my_particle.TrackId();
      fFlashID = i_flash;
      fFlashTime = flash_time;
      fFlashTimeWidth = my_flash.TimeWidth()*ns_per_PMT_tick;
      fFlashY = my_flash.YCenter();
      fFlashZ = my_flash.ZCenter();
      fFlashYWidth = my_flash.YWidth();
      fFlashZWidth = my_flash.ZWidth();
      fFlashPE = my_flash.TotalPE();
      fFlashOnBeamTime = my_flash.OnBeamTime();

      bool beam_match = ((std::abs(particle_time - flash_time )/fFlashTimeWidth)<fTimeMatchMax &&
                         get_MC_particle_origin(my_particle)==simb::kBeamNeutrino);
      bool nonbeam_match = ((std::abs(particle_time - flash_time )/fFlashTimeWidth)<fTimeMatchMax && 
                            get_MC_particle_origin(my_particle)!=simb::kBeamNeutrino);

      if( beam_match ){ 
        fMatchTree_PF->Fill();
        fMatchIndex++;

        //std::cout << "\t\tParticle (x,y,z)=(" << my_particle.Vx() << "," << my_particle.Vy() << "," << my_particle.Vz() << std::endl;
        //std::cout << "\t\tParticle PDG = " << my_particle.PdgCode() << " mother=" << my_particle.Mother() << std::endl;
        //std::cout << "\t\tFlash (y,z) = (" 
        //        << my_flash.YCenter() << " +/- " << my_flash.YWidth() << ","
        //        << my_flash.ZCenter() << " +/- " << my_flash.ZWidth() << std::endl;

        fFlash_match_vector.at(i_flash).particle_indices.push_back(i_particle);
        fParticle_match_vector.at(i_particle).flash_indices.push_back(i_flash);
      }
      else if( nonbeam_match ){ 
        if(lastFlashID!=fFlashID){
          fMatchTree_PF_NotNu->Fill();
          fMatchIndex_NotNu++;
          lastFlashID = fFlashID;
        }
      }


    }//end inner loop over particles
    
  }//end loop over flashes
  
}//end match_flashes_to_particles

void opreco::OpticalRecoAna::match_flashes_to_particles ( std::vector< recob::Track > const &  track_vector, std::vector< simb::MCParticle > const &  particle_vector, geo::Geometry const &  geometry  )

private

Definition at line 406 of file OpticalRecoAna_module.cc.

References compare_particle_and_track(), fParticle_match_vector, and fTrack_match_vector.

                                                                                    {
  bool matching=false;
  
  for(size_t i_track=0; i_track < track_vector.size(); i_track++){

    recob::Track const& my_track( track_vector.at(i_track) );

    for(size_t i_particle=0; i_particle < particle_vector.size(); i_particle++){
      
      simb::MCParticle const& my_particle( particle_vector.at(i_particle) );
      
      matching=false;
      compare_particle_and_track(my_particle,my_track,matching,geometry);
      if(matching){
        fTrack_match_vector.at(i_track).particle_indices.push_back(i_track);
        fParticle_match_vector.at(i_particle).track_indices.push_back(i_track);
      }

    }//end inner loop over particles
    
  }//end loop over tracks
  
}//end match_tracks_to_particles

void opreco::OpticalRecoAna::match_flashes_to_tracks ( std::vector< recob::OpFlash > const &  flash_vector, std::vector< recob::Track > const &  track_vector  )

private

Definition at line 245 of file OpticalRecoAna_module.cc.

References compare_track_and_flash(), fFlash_match_vector, fPEMin, and fTrack_match_vector.

                                                                                               {
  bool matching=false;
  
  for(size_t i_flash=0; i_flash < flash_vector.size(); i_flash++){

    recob::OpFlash const& my_flash( flash_vector.at(i_flash) );
    if(my_flash.TotalPE() < fPEMin) continue;

    for(size_t i_track=0; i_track < track_vector.size(); i_track++){
      
      recob::Track const& my_track( track_vector.at(i_track) );
      
      matching=false;
      compare_track_and_flash(my_track,my_flash,matching);
      if(matching){
        fFlash_match_vector.at(i_flash).track_indices.push_back(i_track);
        fTrack_match_vector.at(i_track).flash_indices.push_back(i_flash);
      }

    }//end inner loop over tracks
    
  }//end loop over flashes
  
}//end match_flashes_to_tracks

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::EventObserverBase::modifiesEvent ( ) const

inlineinherited

Definition at line 25 of file EventObserverBase.h.

  {
    return false;
  }

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

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

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

  {
    return process_name_;
  }

void art::EventObserverBase::registerProducts ( MasterProductRegistry &  , ProductDescriptions &  , ModuleDescription const &    )

inlineinherited

Definition at line 33 of file EventObserverBase.h.

  {}

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

  {
    return selector_config_id_;
  }

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('=');
}

float opreco::OpticalRecoAna::update_MC_particle_time ( simb::MCParticle const &  particle, bool &  pass_check, geo::Geometry const &  geometry  )

private

Definition at line 216 of file OpticalRecoAna_module.cc.

References e, simb::MCParticle::NumberTrajectoryPoints(), geo::GeometryCore::PositionToTPC(), simb::MCParticle::T(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), and simb::MCParticle::Vz().

Referenced by FillMatchTree_PF(), and match_flashes_to_particles().

                                                                                                                                 {

  pass_check = false;

  unsigned int tpc   = 0;
  unsigned int cstat = 0;

  const size_t numtraj = particle.NumberTrajectoryPoints();
  size_t t_iter = 0;
  while(t_iter < numtraj){
    try{
      // check if the particle is inside a TPC                                                                                         
      double pos[3] = {particle.Vx(t_iter), particle.Vy(t_iter), particle.Vz(t_iter)};
      geometry.PositionToTPC(pos, tpc, cstat);
    }
    catch(cet::exception &e){
      t_iter++;
      continue;
    }
    break;
  }

  if(t_iter == numtraj)
    return 0;

  pass_check = true;
  return particle.T(t_iter);
}

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

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

  {
    return wantAllEvents_;
  }

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

  {
    return selectors_.wantEvent(e);
  }

std::string art::EDAnalyzer::workerType ( ) const

inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

    {
      return "WorkerT<EDAnalyzer>";
    }

Member Data Documentation

int opreco::OpticalRecoAna::fEvent

private

Definition at line 106 of file OpticalRecoAna.h.

Referenced by analyze(), and beginJob().

std::vector<flash_match> opreco::OpticalRecoAna::fFlash_match_vector

private

Definition at line 60 of file OpticalRecoAna.h.

Referenced by analyze(), FillMatchTree_PF(), match_flashes_to_particles(), match_flashes_to_tracks(), and OpticalRecoAna().

int opreco::OpticalRecoAna::fFlashID

private

Definition at line 114 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

std::string opreco::OpticalRecoAna::fFlashModuleLabel

private

Definition at line 54 of file OpticalRecoAna.h.

Referenced by analyze(), and OpticalRecoAna().

int opreco::OpticalRecoAna::fFlashOnBeamTime

private

Definition at line 122 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fFlashPE

private

Definition at line 121 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fFlashTime

private

Definition at line 115 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fFlashTimeWidth

private

Definition at line 116 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fFlashY

private

Definition at line 117 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fFlashYWidth

private

Definition at line 119 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fFlashZ

private

Definition at line 118 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fFlashZWidth

private

Definition at line 120 of file OpticalRecoAna.h.

Referenced by FillMatchTree_PF(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fKineticEnergyMin

private

Definition at line 56 of file OpticalRecoAna.h.

Referenced by get_MC_particle_list(), and OpticalRecoAna().

int opreco::OpticalRecoAna::fMatchIndex

private

Definition at line 123 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

int opreco::OpticalRecoAna::fMatchIndex_NotNu

private

Definition at line 124 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

TTree* opreco::OpticalRecoAna::fMatchTree_PF

private

Definition at line 103 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

TTree* opreco::OpticalRecoAna::fMatchTree_PF_NotNu

private

Definition at line 104 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

TH1I* opreco::OpticalRecoAna::fNBeamFlashes

private

Definition at line 101 of file OpticalRecoAna.h.

Referenced by analyze(), and beginJob().

std::vector<particle_match> opreco::OpticalRecoAna::fParticle_match_vector

private

Definition at line 62 of file OpticalRecoAna.h.

Referenced by analyze(), match_flashes_to_particles(), and OpticalRecoAna().

int opreco::OpticalRecoAna::fParticleID

private

Definition at line 107 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

int opreco::OpticalRecoAna::fParticleMother

private

Definition at line 108 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fParticleTime

private

Definition at line 110 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

int opreco::OpticalRecoAna::fParticleTrackID

private

Definition at line 109 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fParticleVx

private

Definition at line 111 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fParticleVy

private

Definition at line 112 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fParticleVz

private

Definition at line 113 of file OpticalRecoAna.h.

Referenced by beginJob(), FillMatchTree_PF(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fPEMin

private

Definition at line 57 of file OpticalRecoAna.h.

Referenced by match_flashes_to_particles(), match_flashes_to_tracks(), and OpticalRecoAna().

int opreco::OpticalRecoAna::fRun

private

Definition at line 105 of file OpticalRecoAna.h.

Referenced by analyze(), and beginJob().

TH1F* opreco::OpticalRecoAna::fTimeDiff

private

Definition at line 99 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

TH1F* opreco::OpticalRecoAna::fTimeDiff_fine

private

Definition at line 100 of file OpticalRecoAna.h.

Referenced by beginJob(), and match_flashes_to_particles().

float opreco::OpticalRecoAna::fTimeMatchMax

private

Definition at line 58 of file OpticalRecoAna.h.

Referenced by match_flashes_to_particles(), and OpticalRecoAna().

std::vector<track_match> opreco::OpticalRecoAna::fTrack_match_vector

private

Definition at line 61 of file OpticalRecoAna.h.

Referenced by match_flashes_to_particles(), match_flashes_to_tracks(), and OpticalRecoAna().

std::string opreco::OpticalRecoAna::fTrackModuleLabel

private

Definition at line 55 of file OpticalRecoAna.h.

Referenced by OpticalRecoAna().

---

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

• larana/v06_21_01/source/larana/OpticalDetector/OpticalRecoAna.h
• larana/v06_21_01/source/larana/OpticalDetector/OpticalRecoAna_module.cc