Inheritance diagram for sbnd::CRTT0Matching:

Public Member Functions
	CRTT0Matching (fhicl::ParameterSet const &p)

	CRTT0Matching (CRTT0Matching const &)=delete

	CRTT0Matching (CRTT0Matching &&)=delete

CRTT0Matching &	operator= (CRTT0Matching const &)=delete

CRTT0Matching &	operator= (CRTT0Matching &&)=delete

void	produce (art::Event &e) override

void	beginJob () override

void	endJob () override

void	reconfigure (fhicl::ParameterSet const &p)

Private Attributes
art::InputTag	fTpcTrackModuleLabel
	name of track producer More...

art::InputTag	fCrtHitModuleLabel
	name of crt producer More...

CRTT0MatchAlg	t0Alg

Detailed Description

Definition at line 61 of file sbndcode/sbndcode/CRT/CRTTools/CRTT0Matching_module.cc.

Constructor & Destructor Documentation

sbnd::CRTT0Matching::CRTT0Matching ( fhicl::ParameterSet const & p )

explicit

Definition at line 96 of file sbndcode/sbndcode/CRT/CRTTools/CRTT0Matching_module.cc.

     : EDProducer(p), t0Alg(p.get<fhicl::ParameterSet>("T0Alg"))
   // Initialize member data here, if know don't want to reconfigure on the fly
   {
 
     // Call appropriate produces<>() functions here.
     produces< std::vector<anab::T0>                  >();
     produces< art::Assns<recob::Track , anab::T0>    >();
     produces< art::Assns<sbn::crt::CRTHit, anab::T0> >();
     
     reconfigure(p);
 
   } // CRTT0Matching()

sbnd::CRTT0Matching::CRTT0Matching ( CRTT0Matching const & )

delete

sbnd::CRTT0Matching::CRTT0Matching ( CRTT0Matching && )

delete

Member Function Documentation

void sbnd::CRTT0Matching::beginJob ( )

override

Definition at line 120 of file sbndcode/sbndcode/CRT/CRTTools/CRTT0Matching_module.cc.

   {
 
   } // CRTT0Matching::beginJob()

void sbnd::CRTT0Matching::endJob ( )

override

Definition at line 192 of file sbndcode/sbndcode/CRT/CRTTools/CRTT0Matching_module.cc.

   {
 
   } // CRTT0Matching::endJob()

CRTT0Matching& sbnd::CRTT0Matching::operator= ( CRTT0Matching const & )

delete

CRTT0Matching& sbnd::CRTT0Matching::operator= ( CRTT0Matching && )

delete

void sbnd::CRTT0Matching::produce ( art::Event & e )

override

Definition at line 125 of file sbndcode/sbndcode/CRT/CRTTools/CRTT0Matching_module.cc.

   {
 
     // Create anab::T0 objects and make association with recob::Track
     std::unique_ptr< std::vector<anab::T0> > T0col( new std::vector<anab::T0>);
     std::unique_ptr< art::Assns<recob::Track, anab::T0> > Trackassn( new art::Assns<recob::Track, anab::T0>);
     std::unique_ptr< art::Assns <sbn::crt::CRTHit, anab::T0> > t0_crthit_assn( new art::Assns<sbn::crt::CRTHit, anab::T0> );
 
     // Retrieve CRT hit list
     art::Handle<std::vector<sbn::crt::CRTHit>> crtListHandle;
     std::vector<art::Ptr<sbn::crt::CRTHit>> crtList;
     if(event.getByLabel(fCrtHitModuleLabel, crtListHandle))
       art::fill_ptr_vector(crtList, crtListHandle);
 
     std::vector<sbn::crt::CRTHit> crtHits;
     for (auto const& crtHit : crtList){
       crtHits.push_back(*crtHit);
     }
 
     // Retrieve track list
     art::Handle< std::vector<recob::Track> > trackListHandle;
     std::vector<art::Ptr<recob::Track> > trackList;
     if (event.getByLabel(fTpcTrackModuleLabel,trackListHandle))
       art::fill_ptr_vector(trackList, trackListHandle);   
 
     mf::LogInfo("CRTT0Matching")
       <<"Number of reconstructed tracks = "<<trackList.size()<<"\n"
       <<"Number of CRT hits = "<<crtList.size();
    
     if (trackListHandle.isValid() && crtListHandle.isValid() ){
 
       auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(event);
       // Loop over all the reconstructed tracks 
       for(size_t track_i = 0; track_i < trackList.size(); track_i++) {
 
         // std::pair<double, double> matchedTime = t0Alg.T0AndDCAFromCRTHits(detProp, *trackList[track_i], crtHits, event);
         matchCand closest = t0Alg.GetClosestCRTHit(detProp, *trackList[track_i], crtHits, event);
 
         if(closest.dca >=0 ){
           mf::LogInfo("CRTT0Matching")
             <<"Matched time = "<<closest.t0<<" [us] to track "<<trackList[track_i]->ID()<<" with DCA = "<<closest.dca;
           T0col->push_back(anab::T0(closest.t0*1e3, trackList[track_i]->ID(),  closest.thishit.plane, (int)closest.extrapLen, closest.dca));
           util::CreateAssn(*this, event, *T0col, trackList[track_i], *Trackassn);
 
           //find this CRThit in the collection
           // note this does not work (both the loop and the assoc !!)
           unsigned CRThitIndex = std::numeric_limits<unsigned>::max();
           for (int ic=0; ic<(int)crtList.size(); ++ic){
             if (crtList[ic]->ts0_ns==closest.thishit.ts0_ns && crtList[ic]->z_pos==closest.thishit.z_pos && crtList[ic]->peshit==closest.thishit.peshit)
               CRThitIndex=ic;
           }
           if (CRThitIndex != std::numeric_limits<unsigned>::max())
             util::CreateAssn(*this, event, *T0col, crtList[CRThitIndex], *t0_crthit_assn);
 
         } // DCA check
 
       } // Loop over tracks  
 
     } // Validity check
 
     event.put(std::move(T0col));
     event.put(std::move(Trackassn));
     event.put(std::move(t0_crthit_assn));
 
   } // CRTT0Matching::produce()

void sbnd::CRTT0Matching::reconfigure ( fhicl::ParameterSet const & p )

Definition at line 111 of file sbndcode/sbndcode/CRT/CRTTools/CRTT0Matching_module.cc.

   {
 
     fTpcTrackModuleLabel = (p.get<art::InputTag> ("TpcTrackModuleLabel"));
     fCrtHitModuleLabel   = (p.get<art::InputTag> ("CrtHitModuleLabel")); 
 
   } // CRTT0Matching::reconfigure()