cosmic::CosmicTrackTagger Class Reference

Inheritance diagram for cosmic::CosmicTrackTagger:

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

Public Member Functions
	CosmicTrackTagger (fhicl::ParameterSet const &p)
virtual	~CosmicTrackTagger ()
void	produce (art::Event &e) override
void	beginJob () override
void	reconfigure (fhicl::ParameterSet const &p)
void	endJob () override
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 Member Functions
double	length (art::Ptr< recob::Track > track)

Private Attributes
std::string	fTrackModuleLabel
int	fEndTickPadding
int	fDetectorWidthTicks
float	fTPCXBoundary
float	fTPCYBoundary
float	fTPCZBoundary
float	fDetHalfHeight
float	fDetWidth
float	fDetLength
int	fMinTickDrift
int	fMaxTickDrift

Detailed Description

Definition at line 60 of file CosmicTrackTagger_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

cosmic::CosmicTrackTagger::CosmicTrackTagger ( fhicl::ParameterSet const &  p )

explicit

Definition at line 96 of file CosmicTrackTagger_module.cc.

References reconfigure().

{
  this->reconfigure(p);

  // Call appropriate Produces<>() functions here.
  produces< std::vector<anab::CosmicTag> >();
  produces< art::Assns<recob::Track, anab::CosmicTag> >();
}

cosmic::CosmicTrackTagger::~CosmicTrackTagger ( )

virtual

Definition at line 107 of file CosmicTrackTagger_module.cc.

                                            {
  // Clean up dynamic memory and other resources here.
}

Member Function Documentation

void cosmic::CosmicTrackTagger::beginJob ( )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 372 of file CosmicTrackTagger_module.cc.

                                      {
}

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::consumes ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

{
  return rng()->createEngine(placeholder_schedule_id(), seed);
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make  )

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make);
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make, label_t const &  engine_label  )

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
}

CurrentProcessingContext const * art::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

  {
    return current_context_.get();
  }

void cosmic::CosmicTrackTagger::endJob ( )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 405 of file CosmicTrackTagger_module.cc.

References DEFINE_ART_MODULE.

                                     {
  // Implementation of optional member function here.
}

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

double cosmic::CosmicTrackTagger::length ( art::Ptr< recob::Track >  track )

private

Definition at line 367 of file CosmicTrackTagger_module.cc.

References recob::Track::Length().

Referenced by produce().

                                                              {
  return track->Length();
}

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 cosmic::CosmicTrackTagger::produce ( art::Event &  e )

overridevirtual

Implements art::EDProducer.

Definition at line 111 of file CosmicTrackTagger_module.cc.

References util::CreateAssn(), recob::Track::End(), fDetHalfHeight, fDetLength, fDetWidth, art::fill_ptr_vector(), fMaxTickDrift, fTPCXBoundary, fTPCYBoundary, fTPCZBoundary, fTrackModuleLabel, art::DataViewImpl::getByLabel(), art::detail::getType(), recob::Track::ID(), art::Ptr< T >::key(), anab::kGeometry_XX, anab::kGeometry_Y, anab::kGeometry_YY, anab::kGeometry_YZ, anab::kGeometry_Z, anab::kGeometry_ZZ, anab::kNotTagged, anab::kOutsideDrift_Partial, length(), recob::Track::Length(), recob::Track::NumberTrajectoryPoints(), art::Event::put(), and recob::Track::Vertex().

                                                  {
  // Implementation of required member function here.

  std::unique_ptr< std::vector< anab::CosmicTag > >              cosmicTagTrackVector( new std::vector<anab::CosmicTag> );
  std::unique_ptr< art::Assns<recob::Track, anab::CosmicTag > >  assnOutCosmicTagTrack( new art::Assns<recob::Track, anab::CosmicTag>);

  TStopwatch ts;

  art::Handle<std::vector<recob::Track> > Trk_h;
  e.getByLabel( fTrackModuleLabel, Trk_h );
  std::vector<art::Ptr<recob::Track> > TrkVec;
  art::fill_ptr_vector(TrkVec, Trk_h);

  // LOOPING OVER INSPILL TRACKS
  
    art::FindManyP<recob::Hit>        hitsSpill   (Trk_h, e, fTrackModuleLabel);
    //    art::FindManyP<recob::Cluster>    ClusterSpill(Trk_h, e, fTrackModuleLabel);

    for( unsigned int iTrack=0; iTrack<Trk_h->size(); iTrack++ ) {

        int isCosmic    =  0;
        anab::CosmicTagID_t tag_id = anab::CosmicTagID_t::kNotTagged;

        art::Ptr<recob::Track>              tTrack  = TrkVec.at(iTrack);
        std::vector<art::Ptr<recob::Hit> >  HitVec  = hitsSpill.at(iTrack);
        
        if (iTrack != tTrack.key())
        {
            std::cout << "Mismatch in track index/key" << std::endl;
        }

        // A BETTER WAY OF FINDING END POINTS:
        auto tVector1 = tTrack->Vertex();
        auto tVector2 = tTrack->End();

        float trackEndPt1_X = tVector1.X();
        float trackEndPt1_Y = tVector1.Y();
        float trackEndPt1_Z = tVector1.Z();
        float trackEndPt2_X = tVector2.X();
        float trackEndPt2_Y = tVector2.Y();
        float trackEndPt2_Z = tVector2.Z();


        if( trackEndPt1_X != trackEndPt1_X ||
            trackEndPt1_Y != trackEndPt1_Y ||
            trackEndPt1_Z != trackEndPt1_Z ||
            trackEndPt2_X != trackEndPt2_X ||
            trackEndPt2_Y != trackEndPt2_Y ||
            trackEndPt2_Z != trackEndPt2_Z ) {
                std::cerr << "!!! FOUND A PROBLEM... the length is: " << tTrack->Length() <<
                    " np: " << tTrack->NumberTrajectoryPoints() << " id: " << tTrack->ID() << " " << tTrack << std::endl;
                //for( size_t hh=0; hh<tTrack->NumberTrajectoryPoints(); hh++) {
                //  std::cerr << hh << " " << tTrack->LocationAtPoint(hh)[0] << ", " <<
                //    tTrack->LocationAtPoint(hh)[1] << ", " <<
                //    tTrack->LocationAtPoint(hh)[2] << std::endl;
                //}
                std::vector<float> tempPt1, tempPt2;
                tempPt1.push_back(-999);
                tempPt1.push_back(-999);
                tempPt1.push_back(-999);
                tempPt2.push_back(-999);
                tempPt2.push_back(-999);
                tempPt2.push_back(-999);
                cosmicTagTrackVector->emplace_back( tempPt1, tempPt2, -999, tag_id);
                util::CreateAssn(*this, e, *cosmicTagTrackVector, tTrack, *assnOutCosmicTagTrack );
                continue; // I don't want to deal with these "tracks"
        }



        // Getting first and last ticks
        float tick1 =  9999;
        float tick2 = -9999;
        bool dumpMe(false);
      
        for ( unsigned int p = 0; p < HitVec.size(); p++) {
            if (dumpMe)
            {
                std::cout << "###>> Hit key: " << HitVec[p].key() << ", peak - RMS: " << HitVec[p]->PeakTimeMinusRMS() << ", peak + RMS: " << HitVec[p]->PeakTimePlusRMS() << std::endl;
            }
            if( HitVec[p]->PeakTimeMinusRMS() < tick1 ) tick1 =  HitVec[p]->PeakTimeMinusRMS();
            if( HitVec[p]->PeakTimePlusRMS()  > tick2 ) tick2 =  HitVec[p]->PeakTimePlusRMS();
        }
      

        // Are any of the ticks outside of the ReadOutWindow ?
        //if(tick1 < fDetectorWidthTicks || tick2 > 2*fDetectorWidthTicks ) {
        if(tick1 < fMinTickDrift || tick2 > fMaxTickDrift ) {
            isCosmic = 1;
            tag_id = anab::CosmicTagID_t::kOutsideDrift_Partial;
        }



        // Now check Y & Z boundaries:
        int nBdY = 0, nBdZ=0;
        if( isCosmic == 0 ) {

            // Checking lower side of TPC
            if( fabs( fDetHalfHeight + trackEndPt1_Y ) < fTPCYBoundary ||
               fabs( fDetHalfHeight + trackEndPt2_Y ) < fTPCYBoundary ||
               trackEndPt1_Y < -fDetHalfHeight ||
               trackEndPt2_Y < -fDetHalfHeight ) nBdY++;
        
            // Checking upper side of TPC
            if( fabs( fDetHalfHeight - trackEndPt1_Y ) < fTPCYBoundary ||
               fabs( fDetHalfHeight - trackEndPt2_Y ) < fTPCYBoundary ||
               trackEndPt1_Y > fDetHalfHeight ||
               trackEndPt2_Y > fDetHalfHeight ) nBdY++;
        
            if(fabs(trackEndPt1_Z - fDetLength)<fTPCZBoundary || fabs(trackEndPt2_Z - fDetLength) < fTPCZBoundary ) nBdZ++;
            if(fabs(trackEndPt1_Z )< fTPCZBoundary || fabs(trackEndPt2_Z )< fTPCZBoundary ) nBdZ++;
            if( (nBdY+nBdZ)>1 ) {
                isCosmic = 2;
                if(nBdY>1) tag_id = anab::CosmicTagID_t::kGeometry_YY;
                else if(nBdZ>1) tag_id = anab::CosmicTagID_t::kGeometry_ZZ;
                else tag_id = anab::CosmicTagID_t::kGeometry_YZ;
            }
            else if( (nBdY+nBdZ)==1) {
                isCosmic = 3 ;
                if(nBdY==1) tag_id = anab::CosmicTagID_t::kGeometry_Y;
                else if(nBdZ==1) tag_id = anab::CosmicTagID_t::kGeometry_Z;
            }
        }

        std::vector<float> endPt1;
        std::vector<float> endPt2;
        endPt1.push_back( trackEndPt1_X );
        endPt1.push_back( trackEndPt1_Y );
        endPt1.push_back( trackEndPt1_Z );
        endPt2.push_back( trackEndPt2_X );
        endPt2.push_back( trackEndPt2_Y );
        endPt2.push_back( trackEndPt2_Z );

        float cosmicScore = isCosmic > 0 ? 1 : 0;
        if( isCosmic==3 ) cosmicScore = 0.5;


        // Doing a very basic check on X boundaries
        // this gets the types of tracks that go through both X boundaries of the detector
        if( fabs( trackEndPt1_X - trackEndPt2_X ) > fDetWidth-fTPCXBoundary ) {
            cosmicScore = 1 ;
            isCosmic = 4;
            tag_id = anab::CosmicTagID_t::kGeometry_XX;
        }

        /*
         // Now check for X boundary
         if( isCosmic==0 ) {
         //anab::CosmicTag cctt = anab::CosmicTag(endPt1, endPt2, tag_id, isCosmic );
         int nBdX =0;
         float xBnd1 = -9999; //cctt.getXInteraction(endPt1[0], 2.0*fDetWidth, fReadOutWindowSize, trackTime, std::floor(tick1) );
         float xBnd2 = -9999; //cctt.getXInteraction(endPt1[0], 2.0*fDetWidth, fReadOutWindowSize, trackTime, std::floor(tick2) );
         if(xBnd1 < fTPCXBoundary || xBnd2 < fTPCXBoundary) nBdX++;
         if( ( fDetWidth - xBnd1 < fTPCXBoundary ) || ( fDetWidth - xBnd1 < fTPCXBoundary ) ) nBdX++;
         if(  nBdX+nBdY+nBdZ>1 && 0 ) isCosmic = 3; // THIS ISN'T SETUP YET -- NEED A HANDLE TO TIME INFO
         if( nBd >0 ) {isCosmic=3; cosmicScore = 0.5;}
         }
         */

        cosmicTagTrackVector->emplace_back( endPt1,
                                          endPt2,
                                          cosmicScore,
                                          tag_id);

//      std::cerr << "The IsCosmic value is "<< isCosmic << " end pts " 
//              << trackEndPt1_X<<","<< trackEndPt1_Y << "," << trackEndPt1_Z<< " | | " 
//              << trackEndPt2_X<< ","<< trackEndPt2_Y <<"," << trackEndPt2_Z << std::endl;
      
      //outTracksForTags->push_back( *tTrack );

        util::CreateAssn(*this, e, *cosmicTagTrackVector, tTrack, *assnOutCosmicTagTrack );
        //util::CreateAssn(*this, e, *cosmicTagTrackVector, HitVec, *assnOutCosmicTagHit);
    }
    // END OF LOOPING OVER INSPILL TRACKS
    

    float         dE=0, dS=0, temp=0, IScore=0;
    unsigned int  IndexE = 0, iTrk1=0, iTrk=0;
    anab::CosmicTagID_t IType=anab::CosmicTagID_t::kNotTagged;
        
    for(iTrk=0; iTrk<Trk_h->size(); iTrk++ ){
        art::Ptr<recob::Track> tTrk  = TrkVec.at(iTrk);
        if ((*cosmicTagTrackVector)[iTrk].CosmicScore()==0){
            auto tStart = tTrk->Vertex();
            auto tEnd   = tTrk->End();
            unsigned int l=0;
            for(iTrk1=0; iTrk1<Trk_h->size(); iTrk1++ ){
                art::Ptr<recob::Track> tTrk1  = TrkVec.at(iTrk1);
                float getScore = (*cosmicTagTrackVector)[iTrk1].CosmicScore();
                if (getScore == 1 || getScore == 0.5){
                    anab::CosmicTagID_t getType = (*cosmicTagTrackVector)[iTrk1].CosmicType();
                    auto tStart1 = tTrk1->Vertex();
                    auto tEnd1   = tTrk1->End();
                    auto NumE    = (tEnd-tStart1).Cross(tEnd-tEnd1);
                    auto DenE    = tEnd1-tStart1;
                    dE = NumE.R()/DenE.R();
                    if (l==0){
                        temp = dE;
                        IndexE = iTrk1;
                        IScore = getScore;
                        IType  = getType;
                    }
                    if (dE<temp){
                        temp = dE;
                        IndexE = iTrk1;
                        IScore = getScore;
                        IType  = getType;
                    }
                    l++;
                }
            }//End Trk1 loop
            art::Ptr<recob::Track> tTrkI = TrkVec.at(IndexE);
            auto tStartI = tTrkI->Vertex();
            auto tEndI   = tTrkI->End();
            auto NumS    = (tStart-tStartI).Cross(tStart-tEndI);
            auto DenS    = tEndI-tStartI;
            dS = NumS.R()/DenS.R();
            if (((dS<5 && temp<5) || (dS<temp && dS<5)) && (length(tTrk)<60)){
                (*cosmicTagTrackVector)[iTrk].CosmicScore() = IScore-0.05;
                (*cosmicTagTrackVector)[iTrk].CosmicType()  = IType;
          //util::CreateAssn(*this, e, *cosmicTagTrackVector, tTrk, *assnOutCosmicTagTrack, iTrk);
            }
        }//end cosmicScore==0 loop
    }//end iTrk loop      
 
 /*std::cout<<"\n"<<Trk_h->size()<<"\t"<<(*cosmicTagTrackVector).size();
   for(unsigned int f=0;f<Trk_h->size();f++){
        std::cout<<"\n\t"<<f<<"\t"<<(*cosmicTagTrackVector)[f].CosmicScore()<<"\t"<<(*cosmicTagTrackVector)[f].CosmicType();
   }*/
 
  // e.put( std::move(outTracksForTags) );
  e.put( std::move(cosmicTagTrackVector) );
  e.put( std::move(assnOutCosmicTagTrack) );

  TrkVec.clear();

} // end of produce

void cosmic::CosmicTrackTagger::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 376 of file CosmicTrackTagger_module.cc.

References fDetectorWidthTicks, fDetHalfHeight, fDetLength, fDetWidth, fEndTickPadding, fMaxTickDrift, fMinTickDrift, fTPCXBoundary, fTPCYBoundary, fTPCZBoundary, fTrackModuleLabel, and fhicl::ParameterSet::get().

Referenced by CosmicTrackTagger().

                                                                     {
  // Implementation of optional member function here.
  
  auto const* detp = lar::providerFrom<detinfo::DetectorPropertiesService>();
  auto const* geo = lar::providerFrom<geo::Geometry>();
  auto const* ts = lar::providerFrom<detinfo::DetectorClocksService>();

  fDetHalfHeight = geo->DetHalfHeight();
  fDetWidth      = 2.*geo->DetHalfWidth();
  fDetLength     = geo->DetLength();

  float fSamplingRate = detp->SamplingRate();

  fTrackModuleLabel = p.get< std::string >("TrackModuleLabel", "track");
  fEndTickPadding   = p.get<    int      >("EndTickPadding",   50);

  fTPCXBoundary = p.get< float >("TPCXBoundary", 5);
  fTPCYBoundary = p.get< float >("TPCYBoundary", 5);
  fTPCZBoundary = p.get< float >("TPCZBoundary", 5);

  const double driftVelocity = detp->DriftVelocity( detp->Efield(), detp->Temperature() ); // cm/us

  //std::cerr << "Drift velocity is " << driftVelocity << " cm/us.  Sampling rate is: "<< fSamplingRate << " detector width: " <<  2*geo->DetHalfWidth() << std::endl;
  fDetectorWidthTicks = 2*geo->DetHalfWidth()/(driftVelocity*fSamplingRate/1000); // ~3200 for uB
  fMinTickDrift = ts->Time2Tick(ts->TriggerTime());
  fMaxTickDrift = fMinTickDrift + fDetectorWidthTicks + fEndTickPadding;
}

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

int cosmic::CosmicTrackTagger::fDetectorWidthTicks

private

Definition at line 89 of file CosmicTrackTagger_module.cc.

Referenced by reconfigure().

float cosmic::CosmicTrackTagger::fDetHalfHeight

private

Definition at line 91 of file CosmicTrackTagger_module.cc.

Referenced by produce(), and reconfigure().

float cosmic::CosmicTrackTagger::fDetLength

private

Definition at line 91 of file CosmicTrackTagger_module.cc.

Referenced by produce(), and reconfigure().

float cosmic::CosmicTrackTagger::fDetWidth

private

Definition at line 91 of file CosmicTrackTagger_module.cc.

Referenced by produce(), and reconfigure().

int cosmic::CosmicTrackTagger::fEndTickPadding

private

Definition at line 88 of file CosmicTrackTagger_module.cc.

Referenced by reconfigure().

int cosmic::CosmicTrackTagger::fMaxTickDrift

private

Definition at line 92 of file CosmicTrackTagger_module.cc.

Referenced by produce(), and reconfigure().

int cosmic::CosmicTrackTagger::fMinTickDrift

private

Definition at line 92 of file CosmicTrackTagger_module.cc.

Referenced by reconfigure().

float cosmic::CosmicTrackTagger::fTPCXBoundary

private

Definition at line 90 of file CosmicTrackTagger_module.cc.

Referenced by produce(), and reconfigure().

float cosmic::CosmicTrackTagger::fTPCYBoundary

private

Definition at line 90 of file CosmicTrackTagger_module.cc.

Referenced by produce(), and reconfigure().

float cosmic::CosmicTrackTagger::fTPCZBoundary

private

Definition at line 90 of file CosmicTrackTagger_module.cc.

Referenced by produce(), and reconfigure().

std::string cosmic::CosmicTrackTagger::fTrackModuleLabel

private

Definition at line 82 of file CosmicTrackTagger_module.cc.

Referenced by produce(), and reconfigure().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larana/v07_04_02/source/larana/CosmicRemoval/CosmicTrackTagger_module.cc