filter::MuonFilter Class Reference

Inheritance diagram for filter::MuonFilter:

Public Types
using	ModuleType = EDFilter
using	WorkerType = WorkerT< EDFilter >
template<typename UserConfig >
using	Table = ProducerBase::Table< UserConfig >

Public Member Functions
	MuonFilter (fhicl::ParameterSet const &)
virtual	~MuonFilter ()
bool	filter (art::Event &evt)
void	beginJob ()
void	endJob ()
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 ()

Static Public Attributes
static constexpr bool	Pass {true}
static constexpr bool	Fail {false}

Protected Member Functions
CurrentProcessingContext const *	currentContext () const
void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
void	prepareForJob (fhicl::ParameterSet const &pset)
void	showMissingConsumes () const

Protected Attributes
std::string	fClusterModuleLabel
std::string	fLineModuleLabel
std::vector< double >	fCuts
double	fDCenter
double	fDelay
double	fTolerance
double	fMaxIon
double	fIonFactor
int	fDeltaWire
	allowed differences in wire number between 2 planes More...

Detailed Description

Definition at line 53 of file MuonFilter_module.cc.

Member Typedef Documentation

using art::EDFilter::ModuleType = EDFilter

inherited

Definition at line 37 of file EDFilter.h.

template<typename UserConfig >

using art::EDFilter::Table = ProducerBase::Table<UserConfig>

inherited

Definition at line 46 of file EDFilter.h.

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

inherited

Definition at line 38 of file EDFilter.h.

Constructor & Destructor Documentation

filter::MuonFilter::MuonFilter ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 83 of file MuonFilter_module.cc.

References reconfigure().

  {
    this->reconfigure(pset);   
  }

filter::MuonFilter::~MuonFilter ( )

virtual

Definition at line 89 of file MuonFilter_module.cc.

  {
  }

Member Function Documentation

void filter::MuonFilter::beginJob ( )

virtual

Reimplemented from art::EDFilter.

Definition at line 108 of file MuonFilter_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::EDFilter::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDFilter.cc.

References art::EDFilter::current_context_.

  {
    return current_context_.get();
  }

void filter::MuonFilter::endJob ( )

virtual

Reimplemented from art::EDFilter.

Definition at line 113 of file MuonFilter_module.cc.

  {
  }

bool filter::MuonFilter::filter ( art::Event &  evt )

virtual

Todo:

: the delta wire numbers seem a bit magic,

Todo:

: should also change to using Geometry::ChannelsIntersect method

Implements art::EDFilter.

Definition at line 118 of file MuonFilter_module.cc.

References a1, a2, DEFINE_ART_MODULE, recob::Cluster::EndTick(), recob::Cluster::EndWire(), fClusterModuleLabel, fCuts, fDCenter, fDelay, fDeltaWire, fIonFactor, fLineModuleLabel, fMaxIon, fTolerance, art::DataViewImpl::getByLabel(), hits(), geo::GeometryCore::IntersectionPoint(), geo::GeometryCore::Nplanes(), recob::Cluster::Plane(), geo::GeometryCore::Plane(), art::PtrVector< T >::push_back(), art::PtrVector< T >::reserve(), art::PtrVector< T >::size(), recob::Cluster::StartTick(), recob::Cluster::StartWire(), geo::PlaneGeo::View(), x1, x2, y1, and y2.

  { 
    art::ServiceHandle<geo::Geometry> geom;
    //    art::ServiceHandle<detinfo::LArPropertiesService> larprop_s;
    //    art::ServiceHandle<detinfo::DetectorPropertiesService> detprop_s;
    auto const * detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
    
    //Drift Velocity in cm/us Sampling rate in ns
    double drift = detprop->DriftVelocity(detprop->Efield(), detprop->Temperature())*detprop->SamplingRate()/1000.0; 

    //This code only works comparing 2 planes so for now these are the 
    // last induction plane and collection plane
    int vPlane = geom->Nplanes() - 1;
    geo::View_t vView = geom->Plane(vPlane).View();
    int uPlane = vPlane-1;
    geo::View_t uView = geom->Plane(uPlane).View();
    art::Handle< std::vector< recob::Cluster > > clustHandle;
    evt.getByLabel(fClusterModuleLabel,clustHandle);

    art::FindManyP<recob::Hit> fmh(clustHandle, evt, fClusterModuleLabel);

    art::PtrVector<recob::Cluster> clusters;
    clusters.reserve(clustHandle->size());
    for(unsigned int i = 0; i < clustHandle->size(); ++i) {
      clusters.push_back(art::Ptr<recob::Cluster>(clustHandle,i));
    }
    double indIon(0),colIon(0);
    std::map<int,int> indMap;
    std::map<int,int> colMap;
    std::vector<std::pair<int,int> > rLook;
    int matchNum=1;
    std::vector<std::vector<double> > tGoing;
    std::vector<std::vector<double> > matched;
    std::vector<double> pointTemp(6);
    std::pair<int,int> pairTemp;
    double ionSum(0.0);
    for(size_t cluster = 0; cluster < clusters.size(); ++cluster) {
      ionSum=0.0;
      std::vector< art::Ptr<recob::Hit> > hits = fmh.at(cluster);
      for(unsigned int hit = 0; hit < hits.size(); hit++) {
        ionSum+=hits[hit]->PeakAmplitude();
      }
      if(clusters[cluster]->View() == uView) indIon+=ionSum;
      else if(clusters[cluster]->View() == vView) colIon+=ionSum;
    }
    mf::LogInfo("MuonFilter") << "Ionizations: " << indIon << " " << colIon ;
    art::Handle<std::vector<recob::Cluster > > lines;
    art::PtrVector<recob::Cluster> inductionSegments, collectionSegments;
    evt.getByLabel(fLineModuleLabel,lines);
    art::PtrVector<recob::Cluster> lineVec;
    lineVec.reserve(lines->size());
    for(unsigned int i = 0; i < lines->size(); ++i) {
      lineVec.push_back(art::Ptr<recob::Cluster>(lines,i));
    }

    for(size_t cl = 0;cl < clusters.size(); cl++) {
      std::vector< art::Ptr<recob::Hit> > hits = fmh.at(cl);
      if (hits.size()>0 && clusters[cl]->View()==uView)
        inductionSegments.push_back(clusters[cl]);
      else if(hits.size()>0 && clusters[cl]->View() == vView) collectionSegments.push_back(clusters[cl]);
    } 
    
    art::FindManyP<recob::Hit> fmhi(inductionSegments,  evt, fClusterModuleLabel);
    art::FindManyP<recob::Hit> fmhc(collectionSegments, evt, fClusterModuleLabel);

    if(inductionSegments.size() == 0 || collectionSegments.size() == 0) { 
      mf::LogInfo("MuonFilter") << "At least one plane with no track";
    }
    else {  

      for(unsigned int i = 0; i < inductionSegments.size(); i++) { 
        if(indMap[i]) continue;
        for(unsigned int j = 0; j < collectionSegments.size(); j++) {
          if(colMap[j]) continue;        

          art::Ptr<recob::Cluster>  indSeg = inductionSegments[i];
          art::Ptr<recob::Cluster>  colSeg = collectionSegments[j];

          std::vector< art::Ptr<recob::Hit> > indHits = fmhi.at(i);
          
          std::vector< art::Ptr<recob::Hit> > colHits = fmhc.at(j);
          
          double trk1Start = indSeg->StartTick()+fDelay;
          double trk1End = indSeg->EndTick()+fDelay;
          double trk2Start =colSeg->StartTick();
          double trk2End =colSeg->EndTick();
          
          int uPos1 = indSeg->StartWire();
          int uPos2 = indSeg->EndWire(); 
          int const vPos1 = colSeg->StartWire();
          int const vPos2 = colSeg->EndWire();
          mf::LogInfo("MuonFilter") << "I J " << i <<" " << j ;
          mf::LogInfo("MuonFilter") << "Start/end " << indSeg->StartWire() 
                                    <<" "<< colSeg->StartWire() 
                                    <<" "<< indSeg->EndWire() 
                                    <<" "<< colSeg->EndWire() ;
          mf::LogInfo("MuonFilter")<<"U's "<< uPos1 <<" " << uPos2 
                                   <<"V's "<< vPos1 <<" " << vPos2 
                                   << " times " << trk1End <<" "<< trk2End 
                                   <<" "<< trk1Start <<" "<< trk2Start ;
          //need to have the corresponding endpoints matched
          //check if they match in this order else switch
          //really should use the crossing function and then have limits
          //on distance outide tpc, or some other way of dealing with
          //imperfect matches
          if((TMath::Abs(uPos1-vPos1)>fDeltaWire || TMath::Abs(uPos2-vPos2)>fDeltaWire) &&
             (TMath::Abs(uPos1-vPos2)<=fDeltaWire && TMath::Abs(uPos2-vPos1)<=fDeltaWire)) {
            mf::LogInfo("MuonFilter") << "Swapped1" ;
            std::swap(uPos1,uPos2);
          }
          //check for time tolerance 
          if((TMath::Abs(trk1Start-trk2Start) > fTolerance && TMath::Abs(trk1End-trk2End) > fTolerance) &&  
             (TMath::Abs(trk1Start-trk2End) < fTolerance && TMath::Abs(trk1End-trk2Start) < fTolerance)) {
            std::swap(trk1Start,trk1End);
            std::swap(uPos1,uPos2);
            mf::LogInfo("MuonFilter") << "Swapped2" ;
          }
          mf::LogInfo("MuonFilter") << "Times: " << trk1Start <<" "<< trk2Start <<" "<<trk1End <<" "<<trk2End;
          //again needs to be fixed
          if((TMath::Abs(trk1Start-trk2Start) < fTolerance && TMath::Abs(trk1End-trk2End) < fTolerance)    && 
             (TMath::Abs(uPos1-vPos1) <=fDeltaWire+2 && TMath::Abs(uPos2-vPos2) <= fDeltaWire+2))
          {
            geo::WireID u_wID1(indSeg->Plane(), uPos1);
            geo::WireID u_wID2(indSeg->Plane(), uPos2);
            geo::WireID v_wID1(colSeg->Plane(), vPos1);
            geo::WireID v_wID2(colSeg->Plane(), vPos2);
            
            double y1, y2, z1, z2;
            geom->IntersectionPoint(u_wID1, v_wID1, y1, z1);
            geom->IntersectionPoint(u_wID2, v_wID2, y2, z2);
            
            double const x1 = (trk1Start+trk2Start)/2.0*drift-fDCenter;
            double const x2 = (trk1End+trk2End)/2.0*drift-fDCenter;
            mf::LogInfo("MuonFilter") <<"Match " << matchNum 
                                      <<" " << x1 << " " << y1 << " " << z1 
                                      <<" " << x2 << " " << y2 << " " << z2;
            bool x1edge,x2edge,y1edge, y2edge,z1edge,z2edge;
            indMap[i]=matchNum;
            colMap[j]=matchNum;
            matchNum++;
            pointTemp[0]=x1;
            pointTemp[1]=y1;
            pointTemp[2]=z1;
            pointTemp[3]=x2;
            pointTemp[4]=y2;
            pointTemp[5]=z2;        
            x1edge =(TMath::Abs(x1) -fCuts[0] > 0);
            x2edge =(TMath::Abs(x2) -fCuts[0] > 0);
            y1edge =(TMath::Abs(y1) -fCuts[1] > 0);
            y2edge =(TMath::Abs(y2) -fCuts[1] > 0);
            z1edge =(TMath::Abs(z1) -fCuts[2] > 0);
            z2edge =(TMath::Abs(z2) -fCuts[2] > 0);
            if((x1edge||y1edge||z1edge) && (x2edge||y2edge||z2edge)) {
              tGoing.push_back(pointTemp);
              mf::LogInfo("MuonFilter") << "outside   Removed induction ion: ";          

              for(size_t h = 0; h < indHits.size(); h++){
                mf::LogInfo("MuonFilter") << indHits[h]->PeakAmplitude() << " ";
                indIon -= indHits[h]->PeakAmplitude();
              }
              mf::LogInfo("MuonFilter")  <<"Removed collection ion: ";

              for(size_t h = 0; h < colHits.size(); h++){ 
                    mf::LogInfo("MuonFilter") << colHits[h]->PeakAmplitude() << " ";
                    colIon -= colHits[h]->PeakAmplitude();
              }
              mf::LogInfo("MuonFilter")<<"Ionization outside track I/C: " << indIon << " "<<colIon;         
            }
            else if((x1edge || y1edge || z1edge)  && 
                    !(x2edge || y2edge || z2edge) && 
                    (z2-z1)>1.2){
              tGoing.push_back(pointTemp);
              mf::LogInfo("MuonFilter") << "stopping   Removed induction ion: ";          
              for(size_t h = 0; h < indHits.size(); h++){
                mf::LogInfo("MuonFilter") <<indHits[h]->PeakAmplitude() << " ";
                indIon -= indHits[h]->PeakAmplitude();
              }
              mf::LogInfo("MuonFilter")  <<"Removed collection ion: ";
              for(size_t h = 0; h < colHits.size(); h++){ 
                mf::LogInfo("MuonFilter") << colHits[h]->PeakAmplitude()<< " ";
                colIon -= colHits[h]->PeakAmplitude();
              }
              mf::LogInfo("MuonFilter") <<"Ionization outside track I/C: " << indIon << " "<<colIon;
            }
            else {
              pairTemp = std::make_pair(i,j);
              mf::LogInfo("MuonFilter") << "rLook matchnum " <<matchNum << " "<<i << " "<<j ; 
              rLook.push_back(pairTemp);
              matched.push_back(pointTemp);  
            }
            break;  //advances i, makes j=0;                    
          }
        }
      }
    }
    //after all matches are made, remove deltas
    double distance=0;
    for(unsigned int i = 0; i < tGoing.size(); i++) 
      for(unsigned int j = 0; j < matched.size();j++){
        mf::LogInfo("MuonFilter") << tGoing.size() <<" " << matched.size() << " " << i << " " << j;
        //test if one is contained within the other in the z-direction 
        if((tGoing[i][2] <= matched[j][2]) && 
           (tGoing[i][5] >= matched[j][5])) { 
          TVector3 a1(&tGoing[i][0]);
          TVector3 a2(&tGoing[i][3]);
          TVector3 b1(&matched[j][0]);
          distance = TMath::Abs((((a1-a2).Cross((a1-a2).Cross(a1-b1))).Unit()).Dot(a1-b1));
          mf::LogInfo("MuonFilter") <<"distance "<< distance;
          if(distance < 6){ 
            mf::LogInfo("MuonFilter") <<"Removing delta ion "<<rLook.size()<<" "<<rLook[j].first<<" "<<matchNum;
            std::vector< art::Ptr<recob::Hit> > temp = fmhi.at(rLook[j].first);
            for(unsigned int h = 0; h < temp.size();h++)         
              indIon -= temp[h]->PeakAmplitude();
            temp = fmhc.at(rLook[j].second);
            for(unsigned int h = 0; h < temp.size();h++)           
              colIon -= temp[h]->PeakAmplitude();   
          }
        }
      } 
    mf::LogInfo("MuonFilter") <<"indIon "<<indIon*fIonFactor <<" colIon " << colIon ;
    if((indIon*fIonFactor > fMaxIon) && (colIon > fMaxIon)) 
      return true;
    else return false;
  }

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::EDFilter::getProductID ( std::string const &  instanceName = {} ) const

inlineinherited

Definition at line 131 of file EDFilter.h.

References art::EDFilter::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 filter::MuonFilter::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 94 of file MuonFilter_module.cc.

References fClusterModuleLabel, fCuts, fDCenter, fDelay, fDeltaWire, fIonFactor, fLineModuleLabel, fMaxIon, fTolerance, and fhicl::ParameterSet::get().

Referenced by MuonFilter().

  {
    fClusterModuleLabel = p.get< std::string         >("ClusterModuleLabel"); 
    fLineModuleLabel    = p.get< std::string               >("LineModuleLabel");    
    fTolerance          = p.get< double                    >("Tolerance");               
    fDelay              = p.get< double                    >("Delay");               
    fDCenter            = p.get< double                    >("DCenter");               
    fMaxIon             = p.get< double                    >("MaxIon");               
    fIonFactor          = p.get< double                    >("IonFactor");          
    fCuts               = p.get< std::vector<double> >("Cuts");
    fDeltaWire          = p.get< int                 >("DeltaWire");
  }

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

constexpr bool art::EDFilter::Fail {false}

staticinherited

Definition at line 33 of file EDFilter.h.

std::string filter::MuonFilter::fClusterModuleLabel

protected

Definition at line 70 of file MuonFilter_module.cc.

Referenced by filter(), and reconfigure().

std::vector<double> filter::MuonFilter::fCuts

protected

Definition at line 72 of file MuonFilter_module.cc.

Referenced by filter(), and reconfigure().

double filter::MuonFilter::fDCenter

protected

Definition at line 73 of file MuonFilter_module.cc.

Referenced by filter(), and reconfigure().

double filter::MuonFilter::fDelay

protected

Definition at line 74 of file MuonFilter_module.cc.

Referenced by filter(), and reconfigure().

int filter::MuonFilter::fDeltaWire

protected

allowed differences in wire number between 2 planes

Definition at line 78 of file MuonFilter_module.cc.

Referenced by filter(), and reconfigure().

double filter::MuonFilter::fIonFactor

protected

Definition at line 77 of file MuonFilter_module.cc.

Referenced by filter(), and reconfigure().

std::string filter::MuonFilter::fLineModuleLabel

protected

Definition at line 71 of file MuonFilter_module.cc.

Referenced by filter(), and reconfigure().

double filter::MuonFilter::fMaxIon

protected

Definition at line 76 of file MuonFilter_module.cc.

Referenced by filter(), and reconfigure().

double filter::MuonFilter::fTolerance

protected

Definition at line 75 of file MuonFilter_module.cc.

Referenced by filter(), and reconfigure().

constexpr bool art::EDFilter::Pass {true}

staticinherited

Definition at line 32 of file EDFilter.h.

---

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

• larevt/v06_22_04/source/larevt/Filters/MuonFilter_module.cc