vertex::VertexFinder2D Class Reference

Inheritance diagram for vertex::VertexFinder2D:

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

Public Member Functions
	VertexFinder2D (fhicl::ParameterSet const &pset)
virtual	~VertexFinder2D ()
void	beginJob ()
void	reconfigure (fhicl::ParameterSet const &p)
void	produce (art::Event &evt)
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
TH1D *	dtIC
std::string	fClusterModuleLabel

Detailed Description

Definition at line 76 of file VertexFinder2D_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

vertex::VertexFinder2D::VertexFinder2D ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 101 of file VertexFinder2D_module.cc.

  {  
    this->reconfigure(pset);    
    produces< std::vector<recob::Vertex> >();
    produces< std::vector<recob::EndPoint2D> >();
    produces< art::Assns<recob::EndPoint2D, recob::Hit> >();
    produces< art::Assns<recob::Vertex, recob::Hit> >();
    produces< art::Assns<recob::Vertex, recob::Shower> >();
    produces< art::Assns<recob::Vertex, recob::Track> >();
  }

vertex::VertexFinder2D::~VertexFinder2D ( )

virtual

Definition at line 112 of file VertexFinder2D_module.cc.

  {
  }

Member Function Documentation

void vertex::VertexFinder2D::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 123 of file VertexFinder2D_module.cc.

References art::TFileDirectory::make().

                               {
    // get access to the TFile service
    art::ServiceHandle<art::TFileService> tfs;
    dtIC = tfs->make<TH1D>("dtIC","It0-Ct0",100,-5,5);
    dtIC->Sumw2();
  }

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

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 vertex::VertexFinder2D::produce ( art::Event &  evt )

virtual

Todo:

should really get the actual vector of hits corresponding to end point

Todo:

for now will get all hits from the current cluster

Todo:

need to actually make tracks and showers to go into 3D vertex

Todo:

currently just passing empty collections to the ctor

Todo:

uncomment following lines when the shower vector actually contains showers from the art::Event

Todo:

uncomment following lines when the track vector actually contains tracks from the art::Event

Implements art::EDProducer.

Definition at line 131 of file VertexFinder2D_module.cc.

References evd::details::begin(), c1, c2, geo::GeometryCore::ChannelsIntersect(), detinfo::DetectorProperties::ConvertTicksToX(), util::CreateAssn(), DEFINE_ART_MODULE, geo::GeometryCore::DetHalfHeight(), detinfo::DetectorProperties::DriftVelocity(), detinfo::DetectorProperties::Efield(), evd::details::end(), art::DataViewImpl::getByLabel(), geo::GeometryCore::GetLArTPCVolumeName(), hits(), CluLen::index, geo::kU, geo::kV, geo::kZ, CluLen::length, LOG_VERBATIM, min, myfunction(), n, geo::GeometryCore::Plane(), geo::GeometryCore::PlaneWireToChannel(), art::PtrVector< T >::push_back(), art::Event::put(), detinfo::DetectorProperties::SamplingRate(), art::PtrVector< T >::size(), SortByWire(), t1, t2, detinfo::DetectorProperties::Temperature(), geo::WireGeo::ThetaZ(), detinfo::DetectorProperties::TriggerOffset(), geo::GeometryCore::Views(), geo::PlaneGeo::Wire(), geo::GeometryCore::WirePitch(), and zz.

  {

    
    art::ServiceHandle<geo::Geometry> geom;
    const detinfo::DetectorProperties* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
    // define TPC parameters
    TString tpcName = geom->GetLArTPCVolumeName();
    
    double YC =  (geom->DetHalfHeight())*2.;

    // wire angle with respect to the vertical direction
    double Angle = geom->Plane(1).Wire(0).ThetaZ(false)-TMath::Pi()/2.; 
    
    // Parameters temporary defined here, but possibly to be retrieved somewhere in the code
    double timetick = detprop->SamplingRate()*1.e-3; //time sample in us
    double presamplings = detprop->TriggerOffset(); //trigger offset

    double wire_pitch   = geom->WirePitch(); //wire pitch in cm
    double Efield_drift = detprop->Efield();      // Electric Field in the drift region in kV/cm
    double Temperature  = detprop->Temperature(); // LAr Temperature in K
    
    //drift velocity in the drift region (cm/us)
    double driftvelocity = detprop->DriftVelocity(Efield_drift,Temperature); 

    //time sample (cm) 
    double timepitch = driftvelocity*timetick; 
    
    art::Handle< std::vector<recob::Cluster> > clusterListHandle;
    evt.getByLabel(fClusterModuleLabel,clusterListHandle);

    art::PtrVector<recob::Cluster> clusters;
    for (unsigned int ii = 0; ii <  clusterListHandle->size(); ++ii){
      art::Ptr<recob::Cluster> clusterHolder(clusterListHandle,ii);
      clusters.push_back(clusterHolder);
    }

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

    //Point to a collection of vertices to output.
    std::unique_ptr<std::vector<recob::Vertex> >                 vcol(new std::vector<recob::Vertex>);          //3D vertex
    std::unique_ptr<std::vector<recob::EndPoint2D> >             epcol(new std::vector<recob::EndPoint2D>);  //2D vertex
    std::unique_ptr< art::Assns<recob::EndPoint2D, recob::Hit> > assnep(new art::Assns<recob::EndPoint2D, recob::Hit>);
    std::unique_ptr< art::Assns<recob::Vertex, recob::Shower> >  assnsh(new art::Assns<recob::Vertex, recob::Shower>);
    std::unique_ptr< art::Assns<recob::Vertex, recob::Track> >   assntr(new art::Assns<recob::Vertex, recob::Track>);
    std::unique_ptr< art::Assns<recob::Vertex, recob::Hit> >     assnh(new art::Assns<recob::Vertex, recob::Hit>);

    // nplanes here is really being used as a proxy for the 
    // number of views in the detector
    int nplanes = geom->Views().size();
        
    std::vector< std::vector<int> > Cls(nplanes); //index to clusters in each view
    std::vector< std::vector<CluLen> > clulens(nplanes);

    std::vector<double> dtdwstart;
      
    //loop over clusters
    for(size_t iclu = 0; iclu < clusters.size(); ++iclu){
      
      float w0 = clusters[iclu]->StartWire();
      float w1 = clusters[iclu]->EndWire();
      float t0 = clusters[iclu]->StartTick();
      float t1 = clusters[iclu]->EndTick();
//      t0 -= detprop->GetXTicksOffset(clusters[iclu]->View(),0,0);
//      t1 -= detprop->GetXTicksOffset(clusters[iclu]->View(),0,0);

      CluLen clulen;
      clulen.index = iclu;
      clulen.length = sqrt(pow((w0-w1)*wire_pitch,2)+pow(detprop->ConvertTicksToX(t0,clusters[iclu]->View(),0,0)-detprop->ConvertTicksToX(t1,clusters[iclu]->View(),0,0),2));

      switch(clusters[iclu]->View()){
        
      case geo::kU :
        clulens[0].push_back(clulen);
        break;
      case geo::kV :
        clulens[1].push_back(clulen);
        break;
      case geo::kZ :
        clulens[2].push_back(clulen);
        break;
      default :
        break;
      }

      std::vector<double> wires;
      std::vector<double> times;
      
      std::vector< art::Ptr<recob::Hit> > hit = fmh.at(iclu);
      std::sort(hit.begin(), hit.end(), SortByWire());
      int n = 0;
      for(size_t i = 0; i < hit.size(); ++i){
        wires.push_back(hit[i]->WireID().Wire);
        times.push_back(hit[i]->PeakTime());
        ++n;
      }
      if(n>=2){
        TGraph *the2Dtrack = new TGraph(std::min(10,n),&wires[0],&times[0]);           
        try{
          the2Dtrack->Fit("pol1","Q");
          TF1 *pol1=(TF1*) the2Dtrack->GetFunction("pol1");
          double par[2];
          pol1->GetParameters(par);
          //std::cout<<iclu<<" "<<par[1]<<" "<<clusters[iclu]->dTdW()<<std::endl;
          dtdwstart.push_back(par[1]);
        }
        catch(...){
          mf::LogWarning("VertexFinder2D") << "Fitter failed";
          delete the2Dtrack;
          dtdwstart.push_back(std::tan(clusters[iclu]->StartAngle()));
          continue;
        }
        delete the2Dtrack;
      }
      else dtdwstart.push_back(std::tan(clusters[iclu]->StartAngle()));
    
    }
    
    //sort clusters based on 2D length
    for (size_t i = 0; i<clulens.size(); ++i){
      std::sort (clulens[i].begin(),clulens[i].end(), myfunction);
      for (size_t j = 0; j<clulens[i].size(); ++j){
        Cls[i].push_back(clulens[i][j].index);
      }
    }

    std::vector< std::vector<int> > cluvtx(nplanes);
    std::vector<double> vtx_w;
    std::vector<double> vtx_t;
    
    for (int i = 0; i < nplanes; ++i){
      if (Cls[i].size() >= 1){
        //at least one cluster
        //find the longest two clusters
        int c1 = -1;
        int c2 = -1;
        double ww0 = -999;
        double wb1 = -999;
        double we1 = -999;
        double wb2 = -999;
        double we2 = -999;
        double tt1 = -999;
        double tt2 = -999;
        double dtdw1 = -999;
        double dtdw2 = -999;
        double lclu1 = -999;
        double lclu2 = -999;
        for (unsigned j = 0; j<Cls[i].size(); ++j){
          double lclu = std::sqrt(pow((clusters[Cls[i][j]]->StartWire()-clusters[Cls[i][j]]->EndWire())*13.5,2)
                                  +pow(clusters[Cls[i][j]]->StartTick()-clusters[Cls[i][j]]->EndTick(),2));
          bool rev = false;
          bool deltaraylike = false;
          bool enoughhits = false;
          if (c1 != -1){
            double wb = clusters[Cls[i][j]]->StartWire();
            double we = clusters[Cls[i][j]]->EndWire();
            double tt = clusters[Cls[i][j]]->StartTick();
            double dtdw = dtdwstart[Cls[i][j]];
            int nhits = fmh.at(Cls[i][j]).size();
            ww0 = (tt-tt1+dtdw1*wb1-dtdw*wb)/(dtdw1-dtdw);
            if (std::abs(wb1-ww0) > std::abs(we1-ww0)) rev = true;//reverse cluster dir
            if ((!rev && ww0 > wb1+15)||(rev && ww0 < we1-15)) deltaraylike = true;
            if (((!rev && ww0 > wb1+10)||(rev && ww0 < we1-10)) && nhits < 5) deltaraylike = true;
            if (wb > wb1+20 && nhits < 20) deltaraylike = true;
            if (wb > wb1+50 && nhits < 20) deltaraylike = true;
            if (wb > wb1+8 && TMath::Abs(dtdw1-dtdw) < 0.15) deltaraylike = true;
            if (std::abs(wb-wb1) > 30 && std::abs(we-we1) > 30) deltaraylike = true;
            if (std::abs(tt-tt1) > 100) deltaraylike = true; //not really deltaray, but isolated cluster
            //make sure there are enough hits in the cluster
            //at leaset 2 hits if goes horizentally, at leaset 4 hits if goes vertically
            double alpha = std::atan(dtdw);
            if (nhits >= int(2+3*(1-std::abs(std::cos(alpha))))) enoughhits = true;
            if (nhits < 5 && (ww0 < wb1-20 || ww0 > we1+20)) enoughhits = false;
            
          }
          //do not replace the second cluster if the 3rd cluster is not consistent with the existing 2
          bool replace = true;
          if (c1 != -1 && c2 != -1){
            double wb = clusters[Cls[i][j]]->StartWire();
            double we = clusters[Cls[i][j]]->EndWire();
            ww0 = (tt2-tt1+dtdw1*wb1-dtdw2*wb2)/(dtdw1-dtdw2);
            if ((std::abs(ww0-wb1) < 10 || std::abs(ww0-we1) < 10) &&
                (std::abs(ww0-wb2) < 10 || std::abs(ww0-we2) < 10)){
              if (std::abs(ww0-wb) > 15 && std::abs(ww0-we) > 15) replace = false;
            }
            //std::cout<<c1<<" "<<c2<<" "<<ww0<<" "<<wb1<<" "<<wb2<<" "<<wb<<" "<<we<<std::endl;
          }
          if (lclu1 < lclu){
            if (c1 != -1 && !deltaraylike && enoughhits){
              lclu2 = lclu1;
              c2 = c1;
              wb2 = wb1;
              we2 = we1;
              tt2 = tt1;
              dtdw2 = dtdw1;
            }
            lclu1 = lclu;
            c1 = Cls[i][j];
            wb1 = clusters[Cls[i][j]]->StartWire();
            we1 = clusters[Cls[i][j]]->EndWire();
            tt1 = clusters[Cls[i][j]]->StartTick();
            if (wb1>we1){
              wb1 = clusters[Cls[i][j]]->EndWire();
              we1 = clusters[Cls[i][j]]->StartWire();
              tt1 = clusters[Cls[i][j]]->EndTick();
            }
            dtdw1 = dtdwstart[Cls[i][j]];
          }
          else if (lclu2 < lclu){
            if (!deltaraylike && enoughhits && replace){
              lclu2 = lclu;
              c2 = Cls[i][j];
              wb2 = clusters[Cls[i][j]]->StartWire();
              we2 = clusters[Cls[i][j]]->EndWire();
              tt2 = clusters[Cls[i][j]]->StartTick();
              dtdw2 = dtdwstart[Cls[i][j]];
            }
          }
        }
        if (c1 != -1 && c2 != -1){
          cluvtx[i].push_back(c1);
          cluvtx[i].push_back(c2);
          
          double w1 = clusters[c1]->StartWire();
          double t1 = clusters[c1]->StartTick();
          if (clusters[c1]->StartWire()>clusters[c1]->EndWire()){
            w1 = clusters[c1]->EndWire();
            t1 = clusters[c1]->EndTick();
          }
          double k1 = dtdwstart[c1];
          double w2 = clusters[c2]->StartWire();
          double t2 = clusters[c2]->StartTick();
          if (clusters[c2]->StartWire()>clusters[c2]->EndWire()){
            w1 = clusters[c2]->EndWire();
            t1 = clusters[c2]->EndTick();
          }
          double k2 = dtdwstart[c2];
//        std::cout<<c1<<" "<<w1<<" "<<t1<<" "<<k1<<" "<<std::endl;
//        std::cout<<c2<<" "<<w2<<" "<<t2<<" "<<k2<<" "<<std::endl;
          //calculate the vertex
          if (std::abs(k1-k2) < 0.5){
            vtx_w.push_back(w1);
            vtx_t.push_back(t1);
          }
          else{
            double t0 = (k1*k2*(w1-w2)+k1*t2-k2*t1)/(k1-k2);
            double w0 = (t2-t1+k1*w1-k2*w2)/(k1-k2);
            vtx_w.push_back(w0);
            vtx_t.push_back(t0);
          }
        }
        else if (Cls[i].size() >= 1){
          if (c1 != -1){
            cluvtx[i].push_back(c1);
            vtx_w.push_back(wb1);
            vtx_t.push_back(tt1);
          }
          else{
            cluvtx[i].push_back(Cls[i][0]);
            vtx_w.push_back(clusters[Cls[i][0]]->StartWire());
            vtx_t.push_back(clusters[Cls[i][0]]->StartTick());
          }
        }
        //save 2D vertex
        // make an empty art::PtrVector of hits
        std::vector< art::Ptr<recob::Hit> > hits = fmh.at(Cls[i][0]);
        double totalQ = 0.;
        for(size_t h = 0; h < hits.size(); ++h) totalQ += hits[h]->Integral();
        
        geo::WireID wireID(hits[0]->WireID().Cryostat,
                           hits[0]->WireID().TPC,
                           hits[0]->WireID().Plane,
                           (unsigned int)vtx_w.back());  //for update to EndPoint2D ... WK 4/22/13
        
        recob::EndPoint2D vertex(vtx_t.back(),
                                 wireID, //for update to EndPoint2D ... WK 4/22/13
                                 1,
                                 epcol->size(),
                                 clusters[Cls[i][0]]->View(),
                                 totalQ);
        epcol->push_back(vertex);
        
        util::CreateAssn(*this, evt, *epcol, hits, *assnep);
        
      }
      else{
        //no cluster found
        vtx_w.push_back(-1);
        vtx_t.push_back(-1);
      }
    }
    //std::cout<<vtx_w[0]<<" "<<vtx_t[0]<<" "<<vtx_w[1]<<" "<<vtx_t[1]<<std::endl;
    
    Double_t vtxcoord[3];
    if (Cls[0].size()>0&&Cls[1].size()>0){//ignore w view
      double Iw0 = (vtx_w[0]+3.95)*wire_pitch;
      double Cw0 = (vtx_w[1]+1.84)*wire_pitch;
      
      double It0 = vtx_t[0] - presamplings;
      It0 *= timepitch;
      double Ct0 = vtx_t[1] - presamplings ;
      Ct0 *= timepitch;
      vtxcoord[0] = detprop->ConvertTicksToX(vtx_t[1],1,0,0);
      vtxcoord[1] = (Cw0-Iw0)/(2.*std::sin(Angle));
      vtxcoord[2] = (Cw0+Iw0)/(2.*std::cos(Angle))-YC/2.*std::tan(Angle);
      
      double yy,zz;       
      if(vtx_w[0]>=0&&vtx_w[0]<=239&&vtx_w[1]>=0&&vtx_w[1]<=239){
        if(geom->ChannelsIntersect(geom->PlaneWireToChannel(0,(int)((Iw0/wire_pitch)-3.95)),      
                                   geom->PlaneWireToChannel(1,(int)((Cw0/wire_pitch)-1.84)),
                                   yy,zz)){
          //channelsintersect provides a slightly more accurate set of y and z coordinates. 
          // use channelsintersect in case the wires in question do cross.
          vtxcoord[1] = yy;
          vtxcoord[2] = zz;
        }
        else{
          vtxcoord[0] = -99999;
          vtxcoord[1] = -99999;
          vtxcoord[2] = -99999;
        }
      } 
      dtIC->Fill(It0-Ct0);
    }
    else{
      vtxcoord[0] = -99999;
      vtxcoord[1] = -99999;
      vtxcoord[2] = -99999;
    }
    
    art::PtrVector<recob::Track> vTracks_vec;
    art::PtrVector<recob::Shower> vShowers_vec;
    
    recob::Vertex the3Dvertex(vtxcoord, vcol->size());
    vcol->push_back(the3Dvertex);
    
    if(vShowers_vec.size() > 0){
      util::CreateAssn(*this, evt, *vcol, vShowers_vec, *assnsh);
      // get the hits associated with each track and associate those with the vertex
      //          for(size_t p = 0; p < vShowers_vec.size(); ++p){
      //            std::vector< art::Ptr<recob::Hit> > hits = fms.at(p);
      //            util::CreateAssn(*this, evt, *vcol, hits, *assnh);
      //          }
    }
    
    if(vTracks_vec.size() > 0){
      util::CreateAssn(*this, evt, *vcol, vTracks_vec, *assntr);
      // get the hits associated with each track and associate those with the vertex
      //          for(size_t p = 0; p < vTracks_vec.size(); ++p){
      //            std::vector< art::Ptr<recob::Hit> > hits = fmt.at(p);
      //            util::CreateAssn(*this, evt, *vcol, hits, *assnh);
      //          }
    }
    
    LOG_VERBATIM("Summary") << std::setfill('-') << std::setw(175) << "-" << std::setfill(' ');
    LOG_VERBATIM("Summary") << "VertexFinder2D Summary:";
    for(size_t i = 0; i<epcol->size(); ++i) LOG_VERBATIM("Summary") << epcol->at(i) ;
    for(size_t i = 0; i<vcol->size(); ++i) LOG_VERBATIM("Summary") << vcol->at(i) ;
    
    evt.put(std::move(epcol));
    evt.put(std::move(vcol));
    evt.put(std::move(assnep));
    evt.put(std::move(assntr));
    evt.put(std::move(assnsh));
    evt.put(std::move(assnh));

  } // end of produce

void vertex::VertexFinder2D::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 117 of file VertexFinder2D_module.cc.

References fhicl::ParameterSet::get().

  {
    fClusterModuleLabel  = p.get< std::string >("ClusterModuleLabel");
    return;
  }

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

TH1D* vertex::VertexFinder2D::dtIC

private

Definition at line 90 of file VertexFinder2D_module.cc.

std::string vertex::VertexFinder2D::fClusterModuleLabel

private

Definition at line 92 of file VertexFinder2D_module.cc.

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/VertexFinder/VertexFinder2D_module.cc