ShowerRecoTools::ShowerTrackColinearTrajPointDirection Class Reference

Inheritance diagram for ShowerRecoTools::ShowerTrackColinearTrajPointDirection:

Public Member Functions
	ShowerTrackColinearTrajPointDirection (const fhicl::ParameterSet &pset)
int	CalculateElement (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder) override

Private Member Functions
int	RunShowerTool (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder, std::string evd_display_name_append="")
virtual void	InitialiseProducers ()
void	SetPtr (art::ProducesCollector *collector)
void	InitaliseProducerPtr (reco::shower::ShowerProducedPtrsHolder &uniqueproducerPtrs)
virtual int	AddAssociations (const art::Ptr< recob::PFParticle > &pfpPtr, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder)
const shower::LArPandoraShowerAlg &	GetLArPandoraShowerAlg () const
template<class T >
art::Ptr< T >	GetProducedElementPtr (std::string Name, reco::shower::ShowerElementHolder &ShowerEleHolder, int iter=-1)
template<class T >
void	InitialiseProduct (std::string Name, std::string InstanceName="")
template<class T , class A , class B >
void	AddSingle (A &a, B &b, std::string Name)
int	GetVectorPtrSize (std::string Name)
void	PrintPtrs ()
void	PrintPtr (std::string Name)

Private Attributes
int	fVerbose
bool	fUsePandoraVertex
bool	fAllowDynamicSliding
bool	fUsePositionInfo
bool	fUseStartPos
float	fAngleCut
std::string	fInitialTrackInputLabel
std::string	fShowerStartPositionInputLabel
std::string	fShowerDirectionOutputLabel

Detailed Description

Definition at line 24 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Constructor & Destructor Documentation

ShowerRecoTools::ShowerTrackColinearTrajPointDirection::ShowerTrackColinearTrajPointDirection

(

const fhicl::ParameterSet &

pset

)

Definition at line 55 of file ShowerTrackColinearTrajPointDirection_tool.cc.

    : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
    , fVerbose(pset.get<int>("Verbose"))
    , fUsePandoraVertex(pset.get<bool>("UsePandoraVertex"))
    , fAllowDynamicSliding(pset.get<bool>("AllowDynamicSliding"))
    , fUsePositionInfo(pset.get<bool>("UsePositionInfo"))
    , fUseStartPos(pset.get<bool>("UseStartPos"))
    , fAngleCut(pset.get<float>("AngleCut"))
    , fInitialTrackInputLabel(pset.get<std::string>("InitialTrackInputLabel"))
    , fShowerStartPositionInputLabel(pset.get<std::string>("ShowerStartPositionInputLabel"))
    , fShowerDirectionOutputLabel(pset.get<std::string>("ShowerDirectionOutputLabel"))

  {}

Member Function Documentation

virtual int ShowerRecoTools::IShowerTool::AddAssociations ( const art::Ptr< recob::PFParticle > &  pfpPtr, art::Event &  Event, reco::shower::ShowerElementHolder &  ShowerEleHolder  )

inlinevirtualinherited

Reimplemented in ShowerRecoTools::ShowerPCADirection, ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder, ShowerRecoTools::ShowerExampleTool, and ShowerRecoTools::ShowerSkeletonTool.

Definition at line 74 of file IShowerTool.h.

    {
      return 0;
    }

template<class T , class A , class B >

void ShowerRecoTools::IShowerTool::AddSingle ( A &  a, B &  b, std::string  Name  )

inlineprotectedinherited

Definition at line 152 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::AddSingle().

    {
      UniquePtrs->AddSingle<T>(a, b, Name);
    }

int ShowerRecoTools::ShowerTrackColinearTrajPointDirection::CalculateElement ( const art::Ptr< recob::PFParticle > &  pfparticle, art::Event &  Event, reco::shower::ShowerElementHolder &  ShowerEleHolder  )

overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 70 of file ShowerTrackColinearTrajPointDirection_tool.cc.

References reco::shower::ShowerElementHolder::CheckElement(), DEFINE_ART_CLASS_TOOL, recob::Track::DirectionAtPoint(), fAllowDynamicSliding, fAngleCut, fInitialTrackInputLabel, recob::Track::FlagsAtPoint(), fShowerDirectionOutputLabel, fShowerStartPositionInputLabel, fUsePandoraVertex, fUsePositionInfo, fUseStartPos, fVerbose, reco::shower::ShowerElementHolder::GetElement(), recob::Track::LocationAtPoint(), recob::TrajectoryPointFlagTraits::NoPoint, recob::Track::NumberTrajectoryPoints(), reco::shower::ShowerElementHolder::SetElement(), and recob::Track::Start().

  {

    //Check the Track has been defined
    if (!ShowerEleHolder.CheckElement(fInitialTrackInputLabel)) {
      if (fVerbose)
        mf::LogError("ShowerTrackColinearTrajPointDirection")
          << "Initial track not set" << std::endl;
      return 1;
    }
    recob::Track InitialTrack;
    ShowerEleHolder.GetElement(fInitialTrackInputLabel, InitialTrack);

    //Smartly choose the which trajectory point to look at by ignoring the smush of hits at the vertex.
    if (InitialTrack.NumberTrajectoryPoints() == 1) {
      if (fVerbose)
        mf::LogError("ShowerTrackColinearTrajPointDirection")
          << "Not Enough trajectory points." << std::endl;
      return 1;
    }

    //Trajectory point which the direction is calcualted for.
    int trajpoint = 0;
    geo::Vector_t Direction;

    if (fUsePositionInfo) {
      //Get the start position.
      geo::Point_t StartPosition;

      if (fUsePandoraVertex) {
        //Check the Track has been defined
        if (!ShowerEleHolder.CheckElement(fShowerStartPositionInputLabel)) {
          if (fVerbose)
            mf::LogError("ShowerTrackColinearTrajPointDirection")
              << "Shower start position not set" << std::endl;
          return 1;
        }
        ShowerEleHolder.GetElement(fShowerStartPositionInputLabel, StartPosition);
      }
      else {
        StartPosition = InitialTrack.Start();
      }

      //Loop over the trajectory points and find two corresponding trajectory points where the angle between themselves (or themsleves and the start position) is less the fMinAngle.
      for (unsigned int traj = 0; traj < InitialTrack.NumberTrajectoryPoints() - 2; ++traj) {
        ++trajpoint;

        //ignore bogus info.
        auto trajflags = InitialTrack.FlagsAtPoint(trajpoint);
        if (trajflags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) { continue; }

        bool bail = false;

        //ignore bogus info.
        auto flags = InitialTrack.FlagsAtPoint(traj);
        if (flags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) { continue; }

        //find the next non bogus traj point.
        int nexttraj = traj + 1;
        auto nextflags = InitialTrack.FlagsAtPoint(nexttraj);
        while (nextflags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) {
          if (nexttraj == (int)InitialTrack.NumberTrajectoryPoints() - 2) {
            bail = true;
            break;
          }
          ++nexttraj;
          nextflags = InitialTrack.FlagsAtPoint(nexttraj);
        }

        //find the next next non bogus traj point.
        int nextnexttraj = nexttraj + 1;
        auto nextnextflags = InitialTrack.FlagsAtPoint(nextnexttraj);
        while (nextnextflags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) {
          if (nexttraj == (int)InitialTrack.NumberTrajectoryPoints() - 1) {
            bail = true;
            break;
          }
          ++nextnexttraj;
          nextnextflags = InitialTrack.FlagsAtPoint(nextnexttraj);
        }

        if (bail) {
          if (fVerbose)
            mf::LogError("ShowerTrackColinearTrajPointDirection")
              << "Trajectory point not set as rest of the traj points are bogus." << std::endl;
          break;
        }

        //Get the directions.
        geo::Vector_t TrajPosition = InitialTrack.LocationAtPoint(traj) - StartPosition;
        geo::Vector_t NextTrajPosition;
        geo::Vector_t NextNextTrajPosition;
        if (fUseStartPos) {
          NextTrajPosition = InitialTrack.LocationAtPoint(nexttraj) - StartPosition;
          NextNextTrajPosition = InitialTrack.LocationAtPoint(nextnexttraj) - StartPosition;
        }
        else {
          NextTrajPosition =
            InitialTrack.LocationAtPoint(nexttraj) - InitialTrack.LocationAtPoint(traj);
          NextNextTrajPosition =
            InitialTrack.LocationAtPoint(nextnexttraj) - InitialTrack.LocationAtPoint(traj + 1);
        }

        //Might still be bogus and we can't use the start point
        if (TrajPosition.R() == 0) { continue; }
        if (NextTrajPosition.R() == 0) { continue; }
        if (NextNextTrajPosition.R() == 0) { continue; }

        //Check to see if the angle between the directions is small enough.
        if (Angle(TrajPosition, NextTrajPosition) < fAngleCut &&
            Angle(TrajPosition, NextNextTrajPosition) < fAngleCut) {
          break;
        }

        //Move the start position onwords.
        if (fAllowDynamicSliding) { StartPosition = InitialTrack.LocationAtPoint(traj); }
      }

      geo::Point_t TrajPosition = InitialTrack.LocationAtPoint(trajpoint);
      Direction = (TrajPosition - StartPosition).Unit();
    }
    else {
      //Loop over the trajectory points and find two corresponding trajectory points where the angle between themselves (or themsleves and the start position) is less the fMinAngle.
      for (unsigned int traj = 0; traj < InitialTrack.NumberTrajectoryPoints() - 2; ++traj) {
        ++trajpoint;

        //ignore bogus info.
        auto trajflags = InitialTrack.FlagsAtPoint(trajpoint);
        if (trajflags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) { continue; }

        //ignore bogus info.
        auto flags = InitialTrack.FlagsAtPoint(traj);
        if (flags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) { continue; }

        bool bail = false;

        geo::Vector_t TrajDirection;

        //Get the next non bogus trajectory points
        if (fUseStartPos) {
          int prevtraj = 0;
          auto prevflags = InitialTrack.FlagsAtPoint(prevtraj);
          while (prevflags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) {
            if (prevtraj == (int)InitialTrack.NumberTrajectoryPoints() - 2) {
              bail = true;
              break;
            }
            ++prevtraj;
            prevflags = InitialTrack.FlagsAtPoint(prevtraj);
          }
          TrajDirection = InitialTrack.DirectionAtPoint(prevtraj);
        }
        else if (fAllowDynamicSliding && traj != 0) {
          int prevtraj = traj - 1;
          auto prevflags = InitialTrack.FlagsAtPoint(prevtraj);
          while (prevflags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) {
            if (prevtraj == 0) {
              bail = true;
              break;
            }
            --prevtraj;
            prevflags = InitialTrack.FlagsAtPoint(prevtraj);
          }
          TrajDirection = InitialTrack.DirectionAtPoint(prevtraj);
        }
        else {
          TrajDirection = InitialTrack.DirectionAtPoint(traj);
        }

        //find the next non bogus traj point.
        int nexttraj = traj + 1;
        auto nextflags = InitialTrack.FlagsAtPoint(nexttraj);
        while (nextflags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) {
          if (nexttraj == (int)InitialTrack.NumberTrajectoryPoints() - 2) {
            bail = true;
            break;
          }
          ++nexttraj;
          nextflags = InitialTrack.FlagsAtPoint(nexttraj);
        }

        //find the next next non bogus traj point.
        int nextnexttraj = nexttraj + 1;
        auto nextnextflags = InitialTrack.FlagsAtPoint(nextnexttraj);
        while (nextnextflags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) {
          if (nexttraj == (int)InitialTrack.NumberTrajectoryPoints() - 1) {
            bail = true;
            break;
          }
          ++nextnexttraj;
          nextnextflags = InitialTrack.FlagsAtPoint(nextnexttraj);
        }

        if (bail) {
          if (fVerbose)
            mf::LogError("ShowerTrackColinearTrajPointDirection")
              << "Trajectory point not set as rest of the traj points are bogus." << std::endl;
          break;
        }

        //Get the directions.
        geo::Vector_t NextTrajDirection = InitialTrack.DirectionAtPoint(nexttraj);
        geo::Vector_t NextNextTrajDirection = InitialTrack.DirectionAtPoint(nextnexttraj);

        //Might still be bogus and we can't use the start point
        if (TrajDirection.R() == 0) { continue; }
        if (NextTrajDirection.R() == 0) { continue; }
        if (NextNextTrajDirection.R() == 0) { continue; }

        //See if the angle is small enough.
        if (Angle(TrajDirection, NextTrajDirection) < fAngleCut &&
            Angle(TrajDirection, NextNextTrajDirection) < fAngleCut) {
          break;
        }
      }
      Direction = InitialTrack.DirectionAtPoint(trajpoint).Unit();
    }

    if (trajpoint == (int)InitialTrack.NumberTrajectoryPoints() - 3) {
      if (fVerbose)
        mf::LogError("ShowerSmartTrackTrajectoryPointDirectio")
          << "Trajectory point not set." << std::endl;
      return 1;
    }

    //Set the direction.
    geo::Vector_t DirectionErr = {-999, -999, -999};
    ShowerEleHolder.SetElement(Direction, DirectionErr, fShowerDirectionOutputLabel);
    return 0;
  }

const shower::LArPandoraShowerAlg& ShowerRecoTools::IShowerTool::GetLArPandoraShowerAlg ( ) const

inlineprotectedinherited

Definition at line 82 of file IShowerTool.h.

Referenced by ShowerRecoTools::ShowerLengthPercentile::CalculateElement(), ShowerRecoTools::ShowerUnidirectiondEdx::CalculateElement(), ShowerRecoTools::Shower3DCylinderTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPCAPropergationStartPosition::CalculateElement(), ShowerRecoTools::Shower2DLinearRegressionTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPFPVertexStartPosition::CalculateElement(), ShowerRecoTools::ShowerDirectionCheater::CalculateElement(), ShowerRecoTools::ShowerPCADirection::CalculateElement(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerTrajPointdEdx::CalculateElement(), ShowerRecoTools::ShowerPCADirection::CalculateShowerPCA(), ShowerRecoTools::Shower2DLinearRegressionTrackHitFinder::FindInitialTrackHits(), ShowerRecoTools::Shower3DCylinderTrackHitFinder::FindTrackSpacePoints(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::FitSegmentAndCalculateResidual(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::PruneFrontOfSPSPool(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::PruneTrack(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::RunTestOfIncrementalSpacePointFinder(), ShowerRecoTools::ShowerTrackPCADirection::ShowerPCAVector(), and ShowerRecoTools::ShowerIncrementalTrackHitFinder::ShowerPCAVector().

    {
      return fLArPandoraShowerAlg;
    };

template<class T >

art::Ptr<T> ShowerRecoTools::IShowerTool::GetProducedElementPtr ( std::string  Name, reco::shower::ShowerElementHolder &  ShowerEleHolder, int  iter = -1  )

inlineprotectedinherited

Definition at line 102 of file IShowerTool.h.

References reco::shower::ShowerElementHolder::CheckElement(), reco::shower::ShowerProducedPtrsHolder::CheckUniqueProduerPtr(), reco::shower::ShowerProducedPtrsHolder::GetArtPtr(), and reco::shower::ShowerElementHolder::GetShowerNumber().

    {

      //Check the element has been set
      bool check_element = ShowerEleHolder.CheckElement(Name);
      if (!check_element) {
        throw cet::exception("IShowerTool") << "tried to get a element that does not exist. Failed "
                                               "at making the art ptr for Element: "
                                            << Name << std::endl;
      }

      //Check the unique ptr has been set.
      bool check_ptr = UniquePtrs->CheckUniqueProduerPtr(Name);
      if (!check_ptr) {
        throw cet::exception("IShowerTool")
          << "tried to get a ptr that does not exist. Failed at making the art ptr for Element"
          << Name;
      }

      //Check if the user has defined an index if not just use the current shower index/
      int index;
      if (iter != -1) { index = iter; }
      else {
        index = ShowerEleHolder.GetShowerNumber();
      }

      //Make the ptr
      return UniquePtrs->GetArtPtr<T>(Name, index);
    }

int ShowerRecoTools::IShowerTool::GetVectorPtrSize ( std::string  Name )

inlineprotectedinherited

Definition at line 158 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::GetVectorPtrSize().

Referenced by ShowerRecoTools::ShowerExampleTool::AddAssociations(), ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::AddAssociations(), and ShowerRecoTools::ShowerPCADirection::AddAssociations().

{ return UniquePtrs->GetVectorPtrSize(Name); }

void ShowerRecoTools::IShowerTool::InitaliseProducerPtr ( reco::shower::ShowerProducedPtrsHolder &  uniqueproducerPtrs )

inlineinherited

Definition at line 68 of file IShowerTool.h.

    {
      UniquePtrs = &uniqueproducerPtrs;
    }

virtual void ShowerRecoTools::IShowerTool::InitialiseProducers ( )

inlinevirtualinherited

Reimplemented in ShowerRecoTools::ShowerPCADirection, ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder, and ShowerRecoTools::ShowerExampleTool.

Definition at line 62 of file IShowerTool.h.

{}

template<class T >

void ShowerRecoTools::IShowerTool::InitialiseProduct ( std::string  Name, std::string  InstanceName = ""  )

inlineprotectedinherited

Definition at line 137 of file IShowerTool.h.

References art::ProducesCollector::produces(), and reco::shower::ShowerProducedPtrsHolder::SetShowerUniqueProduerPtr().

    {

      if (collectorPtr == nullptr) {
        mf::LogWarning("IShowerTool") << "The art::ProducesCollector ptr has not been set";
        return;
      }

      collectorPtr->produces<T>(InstanceName);
      UniquePtrs->SetShowerUniqueProduerPtr(type<T>(), Name, InstanceName);
    }

void ShowerRecoTools::IShowerTool::PrintPtr ( std::string  Name )

inlineprotectedinherited

Definition at line 162 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::PrintPtr().

Referenced by ShowerRecoTools::ShowerExampleTool::CalculateElement().

{ UniquePtrs->PrintPtr(Name); }

void ShowerRecoTools::IShowerTool::PrintPtrs ( )

inlineprotectedinherited

Definition at line 160 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::PrintPtrs().

Referenced by ShowerRecoTools::ShowerExampleTool::CalculateElement().

{ UniquePtrs->PrintPtrs(); }

int ShowerRecoTools::IShowerTool::RunShowerTool ( const art::Ptr< recob::PFParticle > &  pfparticle, art::Event &  Event, reco::shower::ShowerElementHolder &  ShowerEleHolder, std::string  evd_display_name_append = ""  )

inlineinherited

Definition at line 46 of file IShowerTool.h.

    {

      int calculation_status = CalculateElement(pfparticle, Event, ShowerEleHolder);
      if (calculation_status != 0) return calculation_status;
      if (fRunEventDisplay) {
        IShowerTool::GetLArPandoraShowerAlg().DebugEVD(
          pfparticle, Event, ShowerEleHolder, evd_display_name_append);
      }
      return calculation_status;
    }

void ShowerRecoTools::IShowerTool::SetPtr ( art::ProducesCollector *  collector )

inlineinherited

Definition at line 65 of file IShowerTool.h.

{ collectorPtr = collector; }

Member Data Documentation

bool ShowerRecoTools::ShowerTrackColinearTrajPointDirection::fAllowDynamicSliding

private

Definition at line 40 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Referenced by CalculateElement().

float ShowerRecoTools::ShowerTrackColinearTrajPointDirection::fAngleCut

private

Definition at line 48 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackColinearTrajPointDirection::fInitialTrackInputLabel

private

Definition at line 50 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackColinearTrajPointDirection::fShowerDirectionOutputLabel

private

Definition at line 52 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackColinearTrajPointDirection::fShowerStartPositionInputLabel

private

Definition at line 51 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Referenced by CalculateElement().

bool ShowerRecoTools::ShowerTrackColinearTrajPointDirection::fUsePandoraVertex

private

Definition at line 37 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Referenced by CalculateElement().

bool ShowerRecoTools::ShowerTrackColinearTrajPointDirection::fUsePositionInfo

private

Definition at line 42 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Referenced by CalculateElement().

bool ShowerRecoTools::ShowerTrackColinearTrajPointDirection::fUseStartPos

private

Definition at line 45 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Referenced by CalculateElement().

int ShowerRecoTools::ShowerTrackColinearTrajPointDirection::fVerbose

private

Definition at line 36 of file ShowerTrackColinearTrajPointDirection_tool.cc.

Referenced by CalculateElement().

---

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

• larpandora/v09_22_05/source/larpandora/LArPandoraEventBuilding/LArPandoraShower/Tools/ShowerTrackColinearTrajPointDirection_tool.cc