Inheritance diagram for trkf::CosmicTracker:

Public Types
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::Table< UserConfig, KeysToIgnore >

Public Member Functions
	CosmicTracker (fhicl::ParameterSet const &pset)

void	reconfigure (fhicl::ParameterSet const &p)

void	produce (art::Event &evt)

void	beginJob ()

void	endJob ()

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
cluster::ClusterMatchTQ	fClusterMatch

trkf::CosmicTrackerAlg	fCTAlg

std::string	fClusterModuleLabel
	label for input cluster collection More...

std::string	fSortDir
	sort space points More...

bool	fStitchTracks
	Stitch tracks from different TPCs. More...

double	fDisCut
	Distance cut for track merging. More...

double	fAngCut
	Angle cut for track merging. More...

bool	fTrajOnly
	Only use trajectory points from TrackTrajectoryAlg for debugging. More...

Detailed Description

Definition at line 195 of file CosmicTracker_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

trkf::CosmicTracker::CosmicTracker ( fhicl::ParameterSet const & pset )

explicit

Definition at line 230 of file CosmicTracker_module.cc.

References reconfigure().

                                                             :
     fClusterMatch(pset.get< fhicl::ParameterSet >("ClusterMatch")),
     fCTAlg(pset.get< fhicl::ParameterSet >("CTAlg"))
   {
     this->reconfigure(pset);
     produces< std::vector<recob::Track>                        >();
     produces< std::vector<recob::SpacePoint>                   >();
     produces< art::Assns<recob::Track,      recob::Cluster>    >();
     produces< art::Assns<recob::Track,      recob::SpacePoint> >();
     produces< art::Assns<recob::SpacePoint, recob::Hit>        >();
     produces< art::Assns<recob::Track,      recob::Hit>        >();
 
   }

Member Function Documentation

void trkf::CosmicTracker::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 256 of file CosmicTracker_module.cc.

257 {

258 }

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

void trkf::CosmicTracker::endJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 261 of file CosmicTracker_module.cc.

262 {

263 }

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.

 {
   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 trkf::CosmicTracker::produce ( art::Event & evt )

virtual

Implements art::EDProducer.

Definition at line 266 of file CosmicTracker_module.cc.

References evd::details::begin(), cluster::ClusterMatchTQ::ClusterMatch(), recob::tracking::convertCollToPoint(), recob::tracking::convertCollToVector(), util::CreateAssn(), DEFINE_ART_MODULE, dir, trkPoint::dir, evd::details::end(), fAngCut, fClusterMatch, fClusterModuleLabel, fCTAlg, fDisCut, art::fill_ptr_vector(), fSortDir, fStitchTracks, fTrajOnly, art::DataViewImpl::getByLabel(), trkPoint::hit, hits(), util::kBogusD, cluster::ClusterMatchTQ::matchedclusters, MatchTrack(), trkPoint::pos, art::PtrVector< T >::push_back(), art::Event::put(), s, sp_sort_x0(), sp_sort_x1(), sp_sort_y0(), sp_sort_y1(), sp_sort_z0(), sp_sort_z1(), spt_sort_x0(), spt_sort_x1(), spt_sort_y0(), spt_sort_y1(), spt_sort_z0(), spt_sort_z1(), trkf::CosmicTrackerAlg::SPTReco(), tmp, trkf::CosmicTrackerAlg::trajHit, trkf::CosmicTrackerAlg::trajPos, trkf::CosmicTrackerAlg::trkDir, and trkf::CosmicTrackerAlg::trkPos.

                                           {
   
     // get services
     art::ServiceHandle<geo::Geometry> geom;
 
     std::unique_ptr<std::vector<recob::Track>      >              tcol (new std::vector<recob::Track>);           
     std::unique_ptr<std::vector<recob::SpacePoint> >                 spcol(new std::vector<recob::SpacePoint>);
     std::unique_ptr<art::Assns<recob::Track, recob::SpacePoint> > tspassn(new art::Assns<recob::Track, recob::SpacePoint>);
     std::unique_ptr<art::Assns<recob::Track, recob::Cluster> >    tcassn(new art::Assns<recob::Track, recob::Cluster>);
     std::unique_ptr<art::Assns<recob::Track, recob::Hit> >        thassn(new art::Assns<recob::Track, recob::Hit>);
     std::unique_ptr<art::Assns<recob::SpacePoint, recob::Hit> >   shassn(new art::Assns<recob::SpacePoint, recob::Hit>);
 
     //double timetick = detprop->SamplingRate()*1e-3;    //time sample in us
     //double presamplings = detprop->TriggerOffset(); // presamplings in ticks  
     //double plane_pitch = geom->PlanePitch(0,1);   //wire plane pitch in cm 
     //double wire_pitch = geom->WirePitch();    //wire pitch in cm
     //double Efield_drift = larprop->Efield(0);  // Electric Field in the drift region in kV/cm
     //double Temperature = detprop->Temperature();  // LAr Temperature in K
 
     //double driftvelocity = larprop->DriftVelocity(Efield_drift,Temperature);    //drift velocity in the drift region (cm/us)
     //double timepitch = driftvelocity*timetick;                         //time sample (cm) 
 
 
     // get input Cluster object(s).
     art::Handle< std::vector<recob::Cluster> > clusterListHandle;
     std::vector<art::Ptr<recob::Cluster> > clusterlist;
     if (evt.getByLabel(fClusterModuleLabel,clusterListHandle))
       art::fill_ptr_vector(clusterlist, clusterListHandle);
 
     art::FindManyP<recob::Hit> fm(clusterListHandle, evt, fClusterModuleLabel);
 
     // find matched clusters
     fClusterMatch.ClusterMatch(clusterlist,fm);
     std::vector<std::vector<unsigned int> > &matchedclusters = fClusterMatch.matchedclusters;
 
     // get track space points
     std::vector<std::vector<trkPoint>> trkpts(matchedclusters.size());
     for (size_t itrk = 0; itrk<matchedclusters.size(); ++itrk){//loop over tracks
 
       std::vector<art::Ptr<recob::Hit> > hitlist;
       for (size_t iclu = 0; iclu<matchedclusters[itrk].size(); ++iclu){//loop over clusters
 
         std::vector< art::Ptr<recob::Hit> > hits = fm.at(matchedclusters[itrk][iclu]);
         for (size_t ihit = 0; ihit<hits.size(); ++ihit){
           hitlist.push_back(hits[ihit]);
         }
       }
       //reconstruct space points and directions
       fCTAlg.SPTReco(hitlist);
       if (!fTrajOnly){
         if (!fCTAlg.trkPos.size()) continue;
         for (size_t i = 0; i<hitlist.size(); ++i){
           trkPoint trkpt;
           trkpt.pos = fCTAlg.trkPos[i];
           trkpt.dir = fCTAlg.trkDir[i];
           trkpt.hit = hitlist[i];
           trkpts[itrk].push_back(trkpt);
         }
         if (fSortDir=="+x") std::sort(trkpts[itrk].begin(),trkpts[itrk].end(),sp_sort_x0);
         if (fSortDir=="-x") std::sort(trkpts[itrk].begin(),trkpts[itrk].end(),sp_sort_x1);
         if (fSortDir=="+y") std::sort(trkpts[itrk].begin(),trkpts[itrk].end(),sp_sort_y0);
         if (fSortDir=="-y") std::sort(trkpts[itrk].begin(),trkpts[itrk].end(),sp_sort_y1);
         if (fSortDir=="+z") std::sort(trkpts[itrk].begin(),trkpts[itrk].end(),sp_sort_z0);
         if (fSortDir=="-z") std::sort(trkpts[itrk].begin(),trkpts[itrk].end(),sp_sort_z1);
       }
 
       if (fTrajOnly){//debug only
         if (!fCTAlg.trajPos.size()) continue;
         size_t spStart = spcol->size();
         std::vector<recob::SpacePoint> spacepoints;
         //      for (size_t ihit = 0; ihit<hitlist.size(); ++ihit){
         //      if (fCTAlg.usehit[ihit] == 1){
         for (size_t ipt = 0; ipt<fCTAlg.trajPos.size(); ++ipt){
           art::PtrVector<recob::Hit> sp_hits;
           //sp_hits.push_back(hitlist[ihit]);
           double hitcoord[3];
           hitcoord[0] = fCTAlg.trajPos[ipt].X();
           hitcoord[1] = fCTAlg.trajPos[ipt].Y();
           hitcoord[2] = fCTAlg.trajPos[ipt].Z();
 //        if (itrk==1){
 //          std::cout<<"hitcoord "<<hitcoord[0]<<" "<<hitcoord[1]<<" "<<hitcoord[2]<<std::endl;
 //        }
           double err[6] = {util::kBogusD};
           recob::SpacePoint mysp(hitcoord, 
                                  err, 
                                  util::kBogusD, 
                                  spStart + spacepoints.size());//3d point at end of track
           spacepoints.push_back(mysp);
           spcol->push_back(mysp);        
           util::CreateAssn(*this, evt, *spcol, fCTAlg.trajHit[ipt], *shassn);
           //}//
         }//ihit
         size_t spEnd = spcol->size();
         //sort in z direction
         //std::sort(spacepoints.begin(),spacepoints.end(),sp_sort_z0);
         //std::sort(spcol->begin()+spStart,spcol->begin()+spEnd,sp_sort_z0);
         if (fSortDir == "+x"){
           std::sort(spacepoints.begin(),spacepoints.end(),spt_sort_x0);
           std::sort(spcol->begin()+spStart,spcol->begin()+spEnd,spt_sort_x0);
         }
         if (fSortDir == "-x"){
           std::sort(spacepoints.begin(),spacepoints.end(),spt_sort_x1);
           std::sort(spcol->begin()+spStart,spcol->begin()+spEnd,spt_sort_x1);
         }
         if (fSortDir == "+y"){
           std::sort(spacepoints.begin(),spacepoints.end(),spt_sort_y0);
           std::sort(spcol->begin()+spStart,spcol->begin()+spEnd,spt_sort_y0);
         }
         if (fSortDir == "-y"){
           std::sort(spacepoints.begin(),spacepoints.end(),spt_sort_y1);
           std::sort(spcol->begin()+spStart,spcol->begin()+spEnd,spt_sort_y1);
         }
         if (fSortDir == "+z"){
           std::sort(spacepoints.begin(),spacepoints.end(),spt_sort_z0);
           std::sort(spcol->begin()+spStart,spcol->begin()+spEnd,spt_sort_z0);
         }
         if (fSortDir == "-z"){
           std::sort(spacepoints.begin(),spacepoints.end(),spt_sort_z1);
           std::sort(spcol->begin()+spStart,spcol->begin()+spEnd,spt_sort_z1);
         }
 
         if(spacepoints.size()>0){
           
           // make a vector of the trajectory points along the track
           std::vector<TVector3> xyz(spacepoints.size());
           for(size_t s = 0; s < spacepoints.size(); ++s){
             xyz[s] = TVector3(spacepoints[s].XYZ());
           }        
           //Calculate track direction cosines 
           TVector3 startpointVec,endpointVec, DirCos;
           startpointVec = xyz[0];
           endpointVec = xyz.back();
           DirCos = endpointVec - startpointVec;
           //SetMag casues a crash if the magnitude of the vector is zero
           try
             {
               DirCos.SetMag(1.0);//normalize vector
             }
           catch(...){std::cout<<"The Spacepoint is infinitely small"<<std::endl;
             continue;
           }
           //std::cout<<DirCos.x()<<" "<<DirCos.y()<<" "<<DirCos.z()<<std::endl;
           std::vector<TVector3> dircos(spacepoints.size(), DirCos);
           
           tcol->push_back(recob::Track(recob::TrackTrajectory(recob::tracking::convertCollToPoint(xyz),
                                                               recob::tracking::convertCollToVector(dircos),
                                                               recob::Track::Flags_t(xyz.size()), false),
                                        0, -1., 0, recob::tracking::SMatrixSym55(), recob::tracking::SMatrixSym55(), tcol->size()));
           
           // make associations between the track and space points
           util::CreateAssn(*this, evt, *tcol, *spcol, *tspassn, spStart, spEnd);
           
           // now the track and clusters
           //util::CreateAssn(*this, evt, *tcol, clustersPerTrack, *tcassn);
           
           // and the hits and track
           std::vector<art::Ptr<recob::Hit> > trkhits;       
           for (size_t ihit = 0; ihit<hitlist.size(); ++ihit){
             //if (fCTAlg.usehit[ihit] == 1){
             trkhits.push_back(hitlist[ihit]);
             //}
           }
           util::CreateAssn(*this, evt, *tcol, trkhits, *thassn);
         }
       }
     }//itrk
 
     // by default creates a spacepoint and direction for each hit
     if (!fTrajOnly){
       std::vector<std::vector<unsigned int>> trkidx;
       //stitch tracks from different TPCs
       if (fStitchTracks){
         for (size_t itrk1 = 0; itrk1<trkpts.size(); ++itrk1){
           for (size_t itrk2 = itrk1+1; itrk2<trkpts.size(); ++itrk2){
             if (MatchTrack(trkpts[itrk1], trkpts[itrk2], fDisCut, fAngCut)){
               int found1 = -1;
               int found2 = -1;
               for (size_t i = 0; i<trkidx.size(); ++i){
                 for (size_t j = 0; j<trkidx[i].size(); ++j){
                   if (trkidx[i][j]==itrk1) found1 = i;
                   if (trkidx[i][j]==itrk2) found2 = i;
                 }
               }
               //std::cout<<itrk1<<" "<<itrk2<<" "<<found1<<" "<<found2<<std::endl;
               if (found1==-1&&found2==-1){
                 std::vector<unsigned int> tmp;
                 tmp.push_back(itrk1);
                 tmp.push_back(itrk2);
                 trkidx.push_back(tmp);
               }
               else if(found1==-1&&found2!=-1){
                 trkidx[found2].push_back(itrk1);
               }
               else if(found1!=-1&&found2==-1){
                 trkidx[found1].push_back(itrk2);
               }
               else if (found1!=found2){//merge two vectors
                 trkidx[found1].insert(trkidx[found1].end(),
                                       trkidx[found2].begin(),
                                       trkidx[found2].end());
                 trkidx.erase(trkidx.begin()+found2);
               }
             }//found match
           }//itrk2
         }//itrk1
         for (size_t itrk = 0; itrk<trkpts.size(); ++itrk){
           bool found = false;
           for (size_t i = 0; i<trkidx.size(); ++i){
             for (size_t j = 0; j<trkidx[i].size(); ++j){
               if (trkidx[i][j]==itrk) found = true;
             }
           }
           if (!found) {
             std::vector<unsigned int> tmp;
             tmp.push_back(itrk);
             trkidx.push_back(tmp);
           }
         }
       }//stitch
       else{
         trkidx.resize(trkpts.size());
         for (size_t i = 0; i<trkpts.size(); ++i){
           trkidx[i].push_back(i);
         }
       }
       
       //make recob::track and associations
       for (size_t i = 0; i<trkidx.size(); ++i){
         //all track points
         std::vector<trkPoint> finaltrkpts;
         //all the clusters associated with the current track
         std::vector<art::Ptr<recob::Cluster>> clustersPerTrack;
         //all hits
         std::vector<art::Ptr<recob::Hit> > hitlist;
         for(size_t j = 0; j<trkidx[i].size(); ++j){
           for (size_t k = 0; k<trkpts[trkidx[i][j]].size(); ++k){
             finaltrkpts.push_back(trkpts[trkidx[i][j]][k]);
             hitlist.push_back(trkpts[trkidx[i][j]][k].hit);
           }//k
           for (size_t iclu = 0; iclu<matchedclusters[trkidx[i][j]].size(); ++iclu){
             art::Ptr <recob::Cluster> cluster(clusterListHandle,matchedclusters[trkidx[i][j]][iclu]);
             clustersPerTrack.push_back(cluster);
           }
           
         }//j
         if (fStitchTracks){
           if (fSortDir=="+x") std::sort(finaltrkpts.begin(),finaltrkpts.end(),sp_sort_x0);
           if (fSortDir=="-x") std::sort(finaltrkpts.begin(),finaltrkpts.end(),sp_sort_x1);
           if (fSortDir=="+y") std::sort(finaltrkpts.begin(),finaltrkpts.end(),sp_sort_y0);
           if (fSortDir=="-y") std::sort(finaltrkpts.begin(),finaltrkpts.end(),sp_sort_y1);
           if (fSortDir=="+z") std::sort(finaltrkpts.begin(),finaltrkpts.end(),sp_sort_z0);
           if (fSortDir=="-z") std::sort(finaltrkpts.begin(),finaltrkpts.end(),sp_sort_z1);
         }
         size_t spStart = spcol->size();
         std::vector<recob::SpacePoint> spacepoints;
         for (size_t ipt = 0; ipt<finaltrkpts.size(); ++ipt){
           art::PtrVector<recob::Hit> sp_hits;
           sp_hits.push_back(finaltrkpts[ipt].hit);
           double hitcoord[3];
           hitcoord[0] = finaltrkpts[ipt].pos.X();
           hitcoord[1] = finaltrkpts[ipt].pos.Y();
           hitcoord[2] = finaltrkpts[ipt].pos.Z();
           //std::cout<<"hitcoord "<<hitcoord[0]<<" "<<hitcoord[1]<<" "<<hitcoord[2]<<std::endl;
           double err[6] = {util::kBogusD};
           recob::SpacePoint mysp(hitcoord, 
                                  err, 
                                  util::kBogusD, 
                                  spStart + spacepoints.size());//3d point at end of track
           spacepoints.push_back(mysp);
           spcol->push_back(mysp);   
           util::CreateAssn(*this, evt, *spcol, sp_hits, *shassn);
         }//ipt
         size_t spEnd = spcol->size();
         if(spacepoints.size()>0){
           // make a vector of the trajectory points along the track
           std::vector<TVector3> xyz(spacepoints.size());
           std::vector<TVector3> dircos(spacepoints.size());
           for(size_t s = 0; s < spacepoints.size(); ++s){
             xyz[s] = TVector3(spacepoints[s].XYZ());
             dircos[s] = finaltrkpts[s].dir;
             //flip direction if needed.
             if (spacepoints.size()>1){
               if (s==0){
                 TVector3 xyz1 = TVector3(spacepoints[s+1].XYZ());
                 TVector3 dir = xyz1-xyz[s];
                 if (dir.Angle(dircos[s])>0.8*TMath::Pi()){
                   dircos[s] = -dircos[s];
                 }
               }
               else{
                 TVector3 dir = xyz[s]-xyz[s-1];
                 if (dir.Angle(dircos[s])>0.8*TMath::Pi()){
                   dircos[s] = -dircos[s];
                 }
               }
             }
             //std::cout<<s<<" "<<xyz[s].X()<<" "<<xyz[s].Y()<<" "<<xyz[s].Z()<<" "<<dircos[s].X()<<" "<<dircos[s].Y()<<" "<<dircos[s].Z()<<std::endl;
           }        
           tcol->push_back(recob::Track(recob::TrackTrajectory(recob::tracking::convertCollToPoint(xyz),
                                                               recob::tracking::convertCollToVector(dircos),
                                                               recob::Track::Flags_t(xyz.size()), false),
                                        0, -1., 0, recob::tracking::SMatrixSym55(), recob::tracking::SMatrixSym55(), tcol->size()));
           
           // make associations between the track and space points
           util::CreateAssn(*this, evt, *tcol, *spcol, *tspassn, spStart, spEnd);
           
           // now the track and clusters
           util::CreateAssn(*this, evt, *tcol, clustersPerTrack, *tcassn);
           
           // and the hits and track
           std::vector<art::Ptr<recob::Hit> > trkhits;       
           for (size_t ihit = 0; ihit<hitlist.size(); ++ihit){
             //if (fCTAlg.usehit[ihit] == 1){
             trkhits.push_back(hitlist[ihit]);
             //}
           }
           util::CreateAssn(*this, evt, *tcol, trkhits, *thassn);
         }
         
       }//i
     }
     mf::LogVerbatim("Summary") << std::setfill('-') 
                                << std::setw(175) 
                                << "-" 
                                << std::setfill(' ');
     mf::LogVerbatim("Summary") << "CosmicTracker Summary:";
     for(unsigned int i = 0; i<tcol->size(); ++i) mf::LogVerbatim("Summary") << tcol->at(i) ;
     mf::LogVerbatim("Summary") << "CosmicTracker Summary End:";
     
     evt.put(std::move(tcol));
     evt.put(std::move(spcol));
     evt.put(std::move(tspassn));
     evt.put(std::move(tcassn));
     evt.put(std::move(thassn));
     evt.put(std::move(shassn));
 
     return;
   }

void trkf::CosmicTracker::reconfigure ( fhicl::ParameterSet const & p )

Definition at line 245 of file CosmicTracker_module.cc.

References fAngCut, fClusterModuleLabel, fDisCut, fSortDir, fStitchTracks, fTrajOnly, and fhicl::ParameterSet::get().

Referenced by CosmicTracker().

   {
     fClusterModuleLabel     = pset.get< std::string >("ClusterModuleLabel");
     fSortDir                = pset.get< std::string >("SortDirection","+z");
     fStitchTracks           = pset.get< bool   >("StitchTracks");
     fDisCut                 = pset.get< double >("DisCut");
     fAngCut                 = pset.get< double >("AngCut");
     fTrajOnly               = pset.get< bool   >("TrajOnly");
   }

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

double trkf::CosmicTracker::fAngCut

private

Angle cut for track merging.

Definition at line 218 of file CosmicTracker_module.cc.

Referenced by produce(), and reconfigure().

cluster::ClusterMatchTQ trkf::CosmicTracker::fClusterMatch

private

Definition at line 209 of file CosmicTracker_module.cc.

Referenced by produce().

std::string trkf::CosmicTracker::fClusterModuleLabel

private

label for input cluster collection

Definition at line 212 of file CosmicTracker_module.cc.

Referenced by produce(), and reconfigure().

trkf::CosmicTrackerAlg trkf::CosmicTracker::fCTAlg

private

Definition at line 210 of file CosmicTracker_module.cc.

Referenced by produce().

double trkf::CosmicTracker::fDisCut

private

Distance cut for track merging.

Definition at line 217 of file CosmicTracker_module.cc.

Referenced by produce(), and reconfigure().

std::string trkf::CosmicTracker::fSortDir

private

sort space points

Definition at line 214 of file CosmicTracker_module.cc.

Referenced by produce(), and reconfigure().

bool trkf::CosmicTracker::fStitchTracks

private

Stitch tracks from different TPCs.

Definition at line 216 of file CosmicTracker_module.cc.

Referenced by produce(), and reconfigure().

bool trkf::CosmicTracker::fTrajOnly

private

Only use trajectory points from TrackTrajectoryAlg for debugging.

Definition at line 220 of file CosmicTracker_module.cc.

Referenced by produce(), and reconfigure().

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

/cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/TrackFinder/CosmicTracker_module.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation