opdet::TrackTimeAssoc Class Reference

Inheritance diagram for opdet::TrackTimeAssoc:

Public Types
using	ModuleType = EDProducer
using	WorkerType = WorkerT< EDProducer >
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::Table< UserConfig, KeysToIgnore >

Public Member Functions
	TrackTimeAssoc (const fhicl::ParameterSet &)
virtual	~TrackTimeAssoc ()
void	produce (art::Event &)
void	reconfigure (fhicl::ParameterSet const &p)
std::vector< double >	GetMIPHypotheses (trkf::BezierTrack *BTrack, double XOffset=0)
std::vector< std::vector< double > >	ScanMIPHypotheses (trkf::BezierTrack *Btrack)
void	PrintHypotheses (std::vector< std::vector< double > > TrackHypotheses)
double	GetChi2 (std::vector< double > signal, std::vector< double > hypothesis, double UpperLim=0)
double	GetMinChi2 (std::vector< std::vector< double > > ScannedHypotheses, std::vector< double > FlashShape)
void	StoreFlashMatches (std::vector< art::Ptr< recob::Track > > &Tracks, std::vector< art::Ptr< recob::OpFlash > > &Flashes, std::vector< anab::FlashMatch > &Matches, art::Event &evt)
void	beginJob ()
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 Attributes
std::string	fTrackModuleLabel
std::string	fFlashModuleLabel
int	fBezierResolution
int	fPairingMode
double	fLengthCut
double	fPECut

Detailed Description

Definition at line 41 of file TrackTimeAssoc_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

opdet::TrackTimeAssoc::TrackTimeAssoc

(

const fhicl::ParameterSet &

)

Definition at line 146 of file TrackTimeAssoc_module.cc.

  {

    produces< std::vector<anab::FlashMatch> >();
    produces< art::Assns<recob::Track, anab::FlashMatch> >();
    produces< art::Assns<recob::OpFlash, anab::FlashMatch> >();

    this->reconfigure(pset);
   }

opdet::TrackTimeAssoc::~TrackTimeAssoc ( )

virtual

Definition at line 182 of file TrackTimeAssoc_module.cc.

  {
  }

Member Function Documentation

void opdet::TrackTimeAssoc::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 174 of file TrackTimeAssoc_module.cc.

  {
  }

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

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

double opdet::TrackTimeAssoc::GetChi2	(	std::vector< double >	signal,
		std::vector< double >	hypothesis,
		double	UpperLim = 0
	)

Definition at line 502 of file TrackTimeAssoc_module.cc.

  {
       
    double SignalIntegral = 0;
    double HypoIntegral = 0;

    for(size_t i=0; i!=signal.size(); ++i)
      {
        SignalIntegral+=signal.at(i);
        HypoIntegral+=hypothesis.at(i);
      }

    // If light yield indicates >1 MIP light yield, 
    //   Normalize the hypothesis to the signal size.
    //   We do not allow <1 MIP hypotheses.

    double NormFactor = SignalIntegral/HypoIntegral;
    if(NormFactor > 1) 
      {
        for(size_t i=0; i!=hypothesis.size(); ++i)
          {
            hypothesis.at(i) *= NormFactor;
          }
      }
    
    double Chi2=0;
    
    for(size_t i=0; i!=signal.size(); ++i)
      {    
        // We assume width = sqrt(hypothesis)

        if(hypothesis.at(i)!=0)
          Chi2 += pow(hypothesis.at(i) - signal.at(i),2)/(hypothesis.at(i));    
        
        if(Chi2>UpperLim)
          break;
      }
    return Chi2;
  }

double opdet::TrackTimeAssoc::GetMinChi2	(	std::vector< std::vector< double > >	ScannedHypotheses,
		std::vector< double >	FlashShape
	)

Definition at line 258 of file TrackTimeAssoc_module.cc.

  {
    double MinChi2  = 10000;
    if(FlashShape.size()==0) return MinChi2;
    
    for(size_t i=0; i!=ScannedHypotheses.size(); ++i)
      {
        if(ScannedHypotheses.at(i).size()>0)
          {
            double Chi2 = GetChi2(FlashShape, ScannedHypotheses.at(i), MinChi2);
            if(Chi2 < MinChi2)
              {
                MinChi2 = Chi2;
              }
          }
      }
    return MinChi2;
  }

std::vector< double > opdet::TrackTimeAssoc::GetMIPHypotheses	(	trkf::BezierTrack *	BTrack,
		double	XOffset = 0
	)

Definition at line 281 of file TrackTimeAssoc_module.cc.

References trkf::BezierTrack::GetLength(), trkf::BezierTrack::GetTrackPoint(), geo::GeometryCore::NOpDets(), and s.

  {
    art::ServiceHandle<geo::Geometry> geom;
    std::vector<double> ReturnVector(geom->NOpDets(),0);
    
    art::ServiceHandle<phot::PhotonVisibilityService> pvs;

    float TrackLength = Btrack->GetLength();

    double xyz[3];
    for (int b=0; b!=fBezierResolution; b++)
      {
        float s               = float(b) / float(fBezierResolution);
        float dQdx            = 2.1;    // Assume MIP value
        
        Btrack->GetTrackPoint(s,xyz);
        xyz[0]+=XOffset;
        const float* PointVisibility = pvs->GetAllVisibilities(xyz);
        if (!PointVisibility) continue; // point not covered by the service
        float LightAmount = dQdx*TrackLength/float(fBezierResolution);
        
        for(size_t OpDet =0; OpDet!=pvs->NOpChannels();  OpDet++)
          {
            ReturnVector.at(OpDet)+= PointVisibility[OpDet] * LightAmount;
          }
      }
    return ReturnVector;
  }

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

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 opdet::TrackTimeAssoc::PrintHypotheses

(

std::vector< std::vector< double > >

TrackHypotheses

)

Definition at line 545 of file TrackTimeAssoc_module.cc.

  {
    // List the light hypotheses per track, per PMT
    for (size_t i=0; i!=TrackHypotheses.size(); ++i)
      {
        mf::LogVerbatim("TrackTimeAssoc")<< "Visbility for track " << i <<std::endl;

        for(size_t j=0; j!=TrackHypotheses.at(i).size(); ++j)
          {
            mf::LogVerbatim("TrackTimeAssoc") << "Signal at PMT " << j << ", "  << TrackHypotheses.at(i).at(j)<<std::endl;
          }
      }
    
  }

void opdet::TrackTimeAssoc::produce ( art::Event &  evt )

virtual

Implements art::EDProducer.

Definition at line 313 of file TrackTimeAssoc_module.cc.

References evd::details::begin(), f, art::DataViewImpl::getByLabel(), recob::Track::Length(), geo::GeometryCore::NOpChannels(), geo::GeometryCore::NOpDets(), geo::GeometryCore::OpDetFromOpChannel(), recob::OpFlash::TotalPE(), track, and opdet::TrackTimeAssoc_tracksort().

  {
    
    //    int EventID = evt.id().event();
    
    
    // Read in flashes from the event
    art::Handle< std::vector<recob::OpFlash> > flashh;
    evt.getByLabel(fFlashModuleLabel, flashh);
    std::vector<art::Ptr<recob::OpFlash> > Flashes;
    for(unsigned int i=0; i < flashh->size(); ++i)
      {
        art::Ptr<recob::OpFlash> flash(flashh,i);
        if(flash->TotalPE()>fPECut) Flashes.push_back(flash);
      }

    // Read in tracks from the event
    art::Handle< std::vector<recob::Track> > trackh;
    evt.getByLabel(fTrackModuleLabel, trackh);
    std::vector<art::Ptr<recob::Track> >  Tracks;
    for(unsigned int i=0; i < trackh->size(); ++i)
      {
        art::Ptr<recob::Track> track(trackh,i);
        if(track->Length()>fLengthCut) Tracks.push_back(track);
      }

    std::sort(Tracks.begin(), Tracks.end(), TrackTimeAssoc_tracksort);

    // Use these to produce Bezier tracks
    std::vector<bool> TracksToCut(Tracks.size(),false);
    std::vector<trkf::BezierTrack*> BTracks;
    BTracks.clear();
    for(size_t i=0; i!=Tracks.size(); i++)
      BTracks.push_back(new trkf::BezierTrack(*Tracks.at(i)));
        
    art::ServiceHandle<geo::Geometry> geom;
    size_t NOpDets = geom->NOpDets();
    
    std::map<int, bool> OnBeamFlashes;
    
    std::vector<std::vector<std::vector<double> > > TrackHypotheses;
    std::vector<std::vector<double> > FlashShapes;
       

    // For each track
    for (size_t i=0; i!=BTracks.size(); ++i)
      {
        TrackHypotheses.push_back(ScanMIPHypotheses(BTracks.at(i)));
      }
    

    for(size_t f=0; f!=Flashes.size(); ++f)
      {

        std::vector<double> ThisFlashShape(NOpDets,0);
            
        //      if(Flashes.at(f)->InBeamFrame())
        //        {
        for (unsigned int c = 0; c < geom->NOpChannels(); c++){
          unsigned int o = geom->OpDetFromOpChannel(c);
          ThisFlashShape[o]+=Flashes.at(f)->PE(c);
        }
        if(Flashes.at(f)->OnBeamTime()) OnBeamFlashes[f]=true;
        //        }
        FlashShapes.push_back(ThisFlashShape);

      }

    // This map stores the Chi2 for every match:
    //    Chi2Map[Track][Flash] = chi2
    std::map<int, std::map<int, double> > Chi2Map;

    // This map sorts the preferences of each track for each flash
    //    SortedPrefs[TrackID][Chi2] = {flashid1, flashid2, ...}
    std::vector<std::map<double, std::vector<int> > > SortedPrefs(Tracks.size());
    

    for(size_t i=0; i!=TrackHypotheses.size(); ++i)
      {
        for(size_t j=0; j!=FlashShapes.size(); ++j)         
          {
            double Chi2 = GetMinChi2(TrackHypotheses.at(i), FlashShapes.at(j));
            
            Chi2Map[i][j]=Chi2;
            SortedPrefs[i][Chi2].push_back(j);

          }
      }


    // This will hold the list of matches
    std::vector<anab::FlashMatch> Matches;
        


    if(fPairingMode==0)
      {
        // In pairing mode 0, store all combinatoric matches and let the user
        //   deal with Chi2s later
        for(size_t i=0; i!=BTracks.size(); ++i)
          for(size_t j=0; j!=Flashes.size(); ++j)
            Matches.push_back( anab::FlashMatch(Chi2Map[i][j], j, i, (Flashes.at(j)->OnBeamTime()>0)));
      }
    
    else if(fPairingMode==1)
      {
        // In pairing mode 1, use the stable marriage algorithm to make a guess
        //   at good 1<->1 pairings
        
        bool StillPairing =true;
        std::vector<int> FlashesPaired(Flashes.size(),-1);
        std::vector<int> TracksPaired(BTracks.size(),-1);
        
        // If we made a new match in the last round, don't stop
        while(StillPairing)
          {
            StillPairing=false;
            for(size_t i=0; i!=BTracks.size(); ++i)
              {
                // If this track still to be paired
                if(TracksPaired[i]<0)
                  {
                    // Find the flash with best remaining chi2 
                    bool MadeMatch  = false;
                    for(auto itPref = SortedPrefs[i].begin(); itPref!=SortedPrefs[i].end(); ++itPref)
                      {
                        for(size_t iflash =0; iflash!=itPref->second.size(); ++iflash)
                          {
                            int FlashID           = itPref->second.at(iflash);
                            int FlashExistingPref = FlashesPaired[FlashID];
                            if(FlashExistingPref < 0)
                              {
                                // This flash is available - make pairing
                                TracksPaired[i]         = FlashID;
                                
                                // if the flash is on beam, claim to be 
                                //  satisfied, but don't occupy flash
                                // if flash is cosmic, claim it.
                                if(!OnBeamFlashes[FlashID])
                                  FlashesPaired[FlashID]  = i;
                                
                                StillPairing = true;
                                MadeMatch    = true;
                              }
                            else
                              {
                                // This flash taken - flash gets to vote
                                if(Chi2Map[i][FlashID] < Chi2Map[FlashExistingPref][FlashID])
                                  {
                                    // If the flash prefers the new guy, switch
                                    FlashesPaired[FlashID]          = i;
                                    TracksPaired[i]                 = FlashID;
                                    TracksPaired[FlashExistingPref] = -1;
                                    MadeMatch    = true;
                                    StillPairing = true;
                                    break;
                                  }
                                // or else just roll on...                          
                              }
                          }
                        if(MadeMatch)  break;
                      } // end loop over chi2s
                  } // end if unpaired track
              } // end loop over tracks
          } // end loop until no more pairing
        
                
        for(size_t i=0; i!=BTracks.size(); ++i)
          {
            if(TracksPaired[i]>0)
              {
                int TrackID = i;
                int FlashID = TracksPaired[i];
                
                Matches.push_back( anab::FlashMatch(Chi2Map[TrackID][FlashID], FlashID, TrackID, Flashes.at(FlashID)->OnBeamTime() ));
              }
          }
      }


    StoreFlashMatches(Tracks, Flashes, Matches, evt);
    
  }

void opdet::TrackTimeAssoc::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 159 of file TrackTimeAssoc_module.cc.

References fhicl::ParameterSet::get().

  {
    fTrackModuleLabel = pset.get<std::string>("TrackModuleLabel");   
    fFlashModuleLabel = pset.get<std::string>("FlashModuleLabel");
    fBezierResolution = pset.get<int>("BezierResolution");
    fLengthCut        = pset.get<double>("LengthCut");
    fPECut            = pset.get<double>("PECut");
    fPairingMode      = pset.get<int>("PairingMode");
  

  }

std::vector< std::vector< double > > opdet::TrackTimeAssoc::ScanMIPHypotheses

(

trkf::BezierTrack *

Btrack

)

Definition at line 188 of file TrackTimeAssoc_module.cc.

References phot::PhotonVisibilityService::GetAllVisibilities(), trkf::BezierTrack::GetTrackPoint(), trkf::BezierTrack::GetTrajectory(), recob::Trajectory::Length(), phot::PhotonVisibilityService::NOpChannels(), geo::GeometryCore::NOpDets(), s, and recob::Trajectory::Start().

  {

    double MinX   = 0;    //temporary
    double MaxX   = 250;  //temporary
    size_t XSteps = 50;  //temporary

    art::ServiceHandle<geo::Geometry> geom;

    std::vector<std::vector<double> > ReturnVector(XSteps);
    for(size_t i=0; i!=XSteps; ++i)
      ReturnVector[i].resize(geom->NOpDets());
    
    art::ServiceHandle<phot::PhotonVisibilityService> pvs;

    float TrackLength = Btrack->GetTrajectory().Length();
    float OldVertex   = Btrack->GetTrajectory().Start().X();
    
    std::vector<bool> ValidTrajectory(XSteps, true);
    
    double xyz[3];
    for (int b=0; b!=fBezierResolution; b++)
      {
        auto const* larp = lar::providerFrom<detinfo::LArPropertiesService>();

        double MIPYield   = larp->ScintYield();
        double QE         = 0.01;
        double MIPdQdx    = 2.1;
        double PromptFrac = 0.25;
        double PromptMIPScintYield = MIPYield * QE * MIPdQdx * PromptFrac;
        
        //      std::cout<<"check : " << PromptMIPScintYield<<std::endl;
        
        float s               = float(b) / float(fBezierResolution);
        float LightAmount     = PromptMIPScintYield * TrackLength/float(fBezierResolution);
        
        
        for(size_t i=0; i!=XSteps; ++i)
          {
            if(ValidTrajectory[i])
              {
                float NewVertex = MinX + float(i)/float(XSteps)*(MaxX-MinX);    
                Btrack->GetTrackPoint(s,xyz);
                xyz[0] += NewVertex-OldVertex;
                
                if((xyz[0] > MaxX) || (xyz[0] < MinX) ) ValidTrajectory[i]=false; 
                
                const float* PointVisibility = pvs->GetAllVisibilities(xyz);
                if (!PointVisibility) continue; // point not covered by visibility service
                
                for(size_t OpDet =0; OpDet!=pvs->NOpChannels();  OpDet++)
                  {
                    ReturnVector.at(i).at(OpDet) += PointVisibility[OpDet] * LightAmount;
                  }
              }
          }
      }
    for(size_t i=0; i!=ReturnVector.size(); ++i)
      if(!ValidTrajectory[i]) ReturnVector.at(i).clear();
    return ReturnVector;
  }

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 opdet::TrackTimeAssoc::StoreFlashMatches	(	std::vector< art::Ptr< recob::Track > > &	Tracks,
		std::vector< art::Ptr< recob::OpFlash > > &	Flashes,
		std::vector< anab::FlashMatch > &	Matches,
		art::Event &	evt
	)

Definition at line 564 of file TrackTimeAssoc_module.cc.

References util::CreateAssn(), and art::Event::put().

  {
    std::unique_ptr< std::vector<anab::FlashMatch> > flash_matches ( new std::vector<anab::FlashMatch>);
    std::unique_ptr< art::Assns<recob::Track, anab::FlashMatch > > assn_track( new art::Assns<recob::Track, anab::FlashMatch>);
    std::unique_ptr< art::Assns<recob::OpFlash, anab::FlashMatch > > assn_flash( new art::Assns<recob::OpFlash, anab::FlashMatch>);

    for(size_t i=0; i!=Matches.size(); ++i)
      {
        flash_matches->push_back(Matches.at(i));
        
        util::CreateAssn(*this, evt, *(flash_matches.get()), Tracks.at(Matches.at(i).SubjectID()), *(assn_track.get()), i); 

        util::CreateAssn(*this, evt, *(flash_matches.get()), Flashes.at(Matches.at(i).FlashID()), *(assn_flash.get()), i); 
      }
    

    evt.put(std::move(flash_matches));
    evt.put(std::move(assn_track));
    evt.put(std::move(assn_flash));
  }

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

int opdet::TrackTimeAssoc::fBezierResolution

private

Definition at line 66 of file TrackTimeAssoc_module.cc.

std::string opdet::TrackTimeAssoc::fFlashModuleLabel

private

Definition at line 65 of file TrackTimeAssoc_module.cc.

double opdet::TrackTimeAssoc::fLengthCut

private

Definition at line 68 of file TrackTimeAssoc_module.cc.

int opdet::TrackTimeAssoc::fPairingMode

private

Definition at line 67 of file TrackTimeAssoc_module.cc.

double opdet::TrackTimeAssoc::fPECut

private

Definition at line 69 of file TrackTimeAssoc_module.cc.

std::string opdet::TrackTimeAssoc::fTrackModuleLabel

private

Definition at line 64 of file TrackTimeAssoc_module.cc.

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larana/v07_04_02/source/larana/OpticalDetector/TrackTimeAssoc_module.cc