Inheritance diagram for ShowerRecoTools::ShowerLengthPercentile:

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

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

Public Member Functions inherited from ShowerRecoTools::IShowerTool
	IShowerTool (const fhicl::ParameterSet &pset)

virtual	~IShowerTool () noexcept=default

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)

Private Attributes
float	fPercentile

art::InputTag	fPFParticleLabel

int	fVerbose

std::string	fShowerStartPositionInputLabel

std::string	fShowerDirectionInputLabel

std::string	fShowerLengthOutputLabel

std::string	fShowerOpeningAngleOutputLabel

Additional Inherited Members
Protected Member Functions inherited from ShowerRecoTools::IShowerTool
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)

Detailed Description

Definition at line 17 of file ShowerLengthPercentile_tool.cc.

Constructor & Destructor Documentation

ShowerRecoTools::ShowerLengthPercentile::ShowerLengthPercentile ( const fhicl::ParameterSet & pset )

Definition at line 38 of file ShowerLengthPercentile_tool.cc.

     : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
     , fPercentile(pset.get<float>("Percentile"))
     , fPFParticleLabel(pset.get<art::InputTag>("PFParticleLabel"))
     , fVerbose(pset.get<int>("Verbose"))
     , fShowerStartPositionInputLabel(pset.get<std::string>("ShowerStartPositionInputLabel"))
     , fShowerDirectionInputLabel(pset.get<std::string>("ShowerDirectionInputLabel"))
     , fShowerLengthOutputLabel(pset.get<std::string>("ShowerLengthOutputLabel"))
     , fShowerOpeningAngleOutputLabel(pset.get<std::string>("ShowerOpeningAngleOutputLabel"))
   {}

Member Function Documentation

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

overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 50 of file ShowerLengthPercentile_tool.cc.

   {
 
     //Get the start position
     if (!ShowerEleHolder.CheckElement(fShowerStartPositionInputLabel)) {
       if (fVerbose)
         mf::LogError("ShowerLengthPercentile") << "Start position not set, returning " << std::endl;
       return 1;
     }
     //Only consider hits in the same tpcs as the vertex.
     TVector3 ShowerStartPosition = {-999, -999, -999};
     ShowerEleHolder.GetElement(fShowerStartPositionInputLabel, ShowerStartPosition);
 
     // Get the assocated pfParicle Handle
     auto const pfpHandle = Event.getValidHandle<std::vector<recob::PFParticle>>(fPFParticleLabel);
 
     // Get the spacepoint - PFParticle assn
     const art::FindManyP<recob::SpacePoint>& fmspp =
       ShowerEleHolder.GetFindManyP<recob::SpacePoint>(pfpHandle, Event, fPFParticleLabel);
 
     // Get the SpacePoints
     std::vector<art::Ptr<recob::SpacePoint>> spacePoints = fmspp.at(pfparticle.key());
     if (spacePoints.empty()) {
       if (fVerbose)
         mf::LogError("ShowerLengthPercentile") << "No Spacepoints, returning" << std::endl;
       return 1;
     }
 
     if (!ShowerEleHolder.CheckElement(fShowerDirectionInputLabel)) {
       if (fVerbose)
         mf::LogError("ShowerLengthPercentile") << "Direction not set, returning " << std::endl;
       return 1;
     }
 
     TVector3 ShowerDirection = {-999, -999, -999};
     ShowerEleHolder.GetElement(fShowerDirectionInputLabel, ShowerDirection);
 
     //Order the spacepoints
     IShowerTool::GetLArPandoraShowerAlg().OrderShowerSpacePoints(
       spacePoints, ShowerStartPosition, ShowerDirection);
 
     //Find the length as the value that contains % of the hits
     int lengthIter = fPercentile * spacePoints.size();
 
     //Find the length
     double ShowerLength = IShowerTool::GetLArPandoraShowerAlg().SpacePointProjection(
       spacePoints[lengthIter], ShowerStartPosition, ShowerDirection);
     double ShowerMaxProjection = IShowerTool::GetLArPandoraShowerAlg().SpacePointProjection(
       spacePoints[spacePoints.size() - 1], ShowerStartPosition, ShowerDirection);
 
     double ShowerLengthError = ShowerMaxProjection - ShowerLength;
 
     //Order the spacepoints in perpendicular
     IShowerTool::GetLArPandoraShowerAlg().OrderShowerSpacePointsPerpendicular(
       spacePoints, ShowerStartPosition, ShowerDirection);
 
     //Find the length as the value that contains % of the hits
     int perpIter = fPercentile * spacePoints.size();
 
     //Find the width of the shower
     double ShowerWidth = IShowerTool::GetLArPandoraShowerAlg().SpacePointPerpendicular(
       spacePoints[perpIter], ShowerStartPosition, ShowerDirection);
     // double ShowerMaxWidth = IShowerTool::GetLArPandoraShowerAlg().SpacePointPerpendicular(
     //     spacePoints[spacePoints.size() -1], ShowerStartPosition, ShowerDirection);
 
     double ShowerAngle = std::atan(ShowerWidth / ShowerLength);
     double ShowerAngleError = -999; //TODO: Do properly
 
     // Fill the shower element holder
     ShowerEleHolder.SetElement(ShowerLength, ShowerLengthError, fShowerLengthOutputLabel);
     ShowerEleHolder.SetElement(ShowerAngle, ShowerAngleError, fShowerOpeningAngleOutputLabel);
 
     return 0;
   }