cvn::CVNProtoDUNEUtils Class Reference

#include <CVNProtoDUNEUtils.h>

Public Member Functions
	CVNProtoDUNEUtils ()
	~CVNProtoDUNEUtils ()
const std::vector< const recob::Hit * >	GetRecoSliceHits (const recob::Slice &slice, art::Event const &evt, const std::string sliceModule) const
const std::vector< const recob::Hit * >	GetRecoSliceHits (unsigned int sliceID, art::Event const &evt, const std::string sliceModule) const
const std::map< unsigned int, std::vector< const recob::Hit * > >	GetRecoSliceHitMap (art::Event const &evt, const std::string sliceModule) const
unsigned short	GetBeamSlice (art::Event const &evt, const std::string particleLabel) const
	Try to get the slice tagged as beam. Returns 9999 if no beam slice was found. More...
const std::map< unsigned int, std::vector< const recob::PFParticle * > >	GetPFParticleSliceMap (art::Event const &evt, const std::string particleLabel) const
	Get a map of slice index to the primary PFParticles within it. More...
const std::map< unsigned int, std::vector< const recob::PFParticle * > >	GetAllPFParticleSliceMap (art::Event const &evt, const std::string particleLabel) const
	Get a map of slice index to all of the PFParticles within it. More...
unsigned short	GetPFParticleSliceIndex (const recob::PFParticle &particle, art::Event const &evt, const std::string particleLabel) const
	Get the reconstructed slice number associated with a particle. More...
const recob::Slice *	GetPFParticleSlice (const recob::PFParticle &particle, art::Event const &evt, const std::string particleLabel) const
	Get the reconstructed slice associated with a particle. More...
const std::vector< const recob::SpacePoint * >	GetPFParticleSpacePoints (const recob::PFParticle &particle, art::Event const &evt, const std::string particleLabel) const
	Get the SpacePoints associated to the PFParticle. More...
bool	IsBeamParticle (const recob::PFParticle &particle, art::Event const &evt, const std::string particleLabel) const
	Use the pandora metadata to tell us if this is a beam particle or not. More...

Private Member Functions
const std::map< unsigned int, std::vector< const recob::PFParticle * > >	SliceMapHelper (art::Event const &evt, const std::string particleLabel, bool primaryOnly) const
	Helper to get the slice map and avoid code repetition. More...
const std::map< std::string, float >	GetPFParticleMetaData (const recob::PFParticle &particle, art::Event const &evt, const std::string particleLabel) const
	Get the metadata associated to a PFParticle from pandora. More...

Detailed Description

Definition at line 26 of file CVNProtoDUNEUtils.h.

Constructor & Destructor Documentation

cvn::CVNProtoDUNEUtils::CVNProtoDUNEUtils

(

)

Definition at line 8 of file CVNProtoDUNEUtils.cxx.

                                       {

}

cvn::CVNProtoDUNEUtils::~CVNProtoDUNEUtils

(

)

Definition at line 12 of file CVNProtoDUNEUtils.cxx.

                                        {

}

Member Function Documentation

const std::map< unsigned int, std::vector< const recob::PFParticle * > > cvn::CVNProtoDUNEUtils::GetAllPFParticleSliceMap	(	art::Event const &	evt,
		const std::string	particleLabel
	)		const

Get a map of slice index to all of the PFParticles within it.

Definition at line 86 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                                   {

  return SliceMapHelper(evt,particleLabel,false);

}

unsigned short cvn::CVNProtoDUNEUtils::GetBeamSlice	(	art::Event const &	evt,
		const std::string	particleLabel
	)		const

Try to get the slice tagged as beam. Returns 9999 if no beam slice was found.

Definition at line 62 of file CVNProtoDUNEUtils.cxx.

                                                                                                         {

  const std::map<unsigned int, std::vector<const recob::PFParticle*>> sliceMap = GetPFParticleSliceMap(evt,particleLabel);

  for(auto slice : sliceMap){
    for(auto particle : slice.second){
      if(IsBeamParticle(*particle,evt,particleLabel)){
        return slice.first;
      }
    }
  }

  return 9999;

}

const std::map< std::string, float > cvn::CVNProtoDUNEUtils::GetPFParticleMetaData ( const recob::PFParticle &  particle, art::Event const &  evt, const std::string  particleLabel  ) const

private

Get the metadata associated to a PFParticle from pandora.

Definition at line 169 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                                   {
  // Get the particles
  auto pfParticles = evt.getValidHandle<std::vector<recob::PFParticle>>(particleLabel);
  // And their meta data
  const art::FindManyP<larpandoraobj::PFParticleMetadata> findMetaData(pfParticles,evt,particleLabel);

  const larpandoraobj::PFParticleMetadata metaData = *((findMetaData.at(particle.Self())).at(0));

  return metaData.GetPropertiesMap();
}

const recob::Slice * cvn::CVNProtoDUNEUtils::GetPFParticleSlice	(	const recob::PFParticle &	particle,
		art::Event const &	evt,
		const std::string	particleLabel
	)		const

Get the reconstructed slice associated with a particle.

Definition at line 135 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                   {

  // Perhaps we should use the associations to do this? 
  auto pfParticles = evt.getValidHandle<std::vector<recob::PFParticle>>(particleLabel);
  const art::FindOneP<recob::Slice> findSlice(pfParticles,evt,particleLabel);

  const recob::Slice* slice = findSlice.at(particle.Self()).get();

  return slice;
}

unsigned short cvn::CVNProtoDUNEUtils::GetPFParticleSliceIndex	(	const recob::PFParticle &	particle,
		art::Event const &	evt,
		const std::string	particleLabel
	)		const

Get the reconstructed slice number associated with a particle.

Definition at line 147 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                     {

  // Try to use slices if we can
  try{
    const recob::Slice* slice = GetPFParticleSlice(particle,evt,particleLabel);
    return slice->ID();
  }
  // Otherwise fall back on metadata
  catch(...){
    std::map<std::string,float> mdMap = GetPFParticleMetaData(particle,evt,particleLabel);
    std::string search = "SliceIndex";
    if(mdMap.find(search) != mdMap.end()){
      return static_cast<unsigned short>(mdMap.at(search));
    }
    else{
//    std::cerr << "Object has no slice index... returning 9999" << std::endl;
      return 9999;
    }
  }

}

const std::map< unsigned int, std::vector< const recob::PFParticle * > > cvn::CVNProtoDUNEUtils::GetPFParticleSliceMap	(	art::Event const &	evt,
		const std::string	particleLabel
	)		const

Get a map of slice index to the primary PFParticles within it.

Definition at line 79 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                                {

  return SliceMapHelper(evt,particleLabel,true);

}

const std::vector< const recob::SpacePoint * > cvn::CVNProtoDUNEUtils::GetPFParticleSpacePoints	(	const recob::PFParticle &	particle,
		art::Event const &	evt,
		const std::string	particleLabel
	)		const

Get the SpacePoints associated to the PFParticle.

Definition at line 118 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                                               {

  // Get the particles and their associations
  auto particles = evt.getValidHandle<std::vector<recob::PFParticle>>(particleLabel);
  const art::FindManyP<recob::SpacePoint> findSpacePoints(particles,evt,particleLabel);
  const std::vector<art::Ptr<recob::SpacePoint>> pfpSpacePoints = findSpacePoints.at(particle.Self());

  // We don't want the art::Ptr so we need to get rid of it
  std::vector<const recob::SpacePoint*> sp;
  for(auto pointer : pfpSpacePoints){
    sp.push_back(pointer.get());
  }

  return sp;
}

const std::map< unsigned int, std::vector< const recob::Hit * > > cvn::CVNProtoDUNEUtils::GetRecoSliceHitMap	(	art::Event const &	evt,
		const std::string	sliceModule
	)		const

Definition at line 43 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                     {

  auto recoSlices = evt.getValidHandle<std::vector<recob::Slice> >(sliceModule);
  std::map<unsigned int, std::vector<const recob::Hit*>> hitMap;

  for(auto const slice : *recoSlices){

    const std::vector<const recob::Hit*> constvec = GetRecoSliceHits(slice.ID(),evt,sliceModule);
    for(auto const h : constvec){
      hitMap[slice.ID()].push_back(h);
    }

  }

  return hitMap;

}

const std::vector< const recob::Hit * > cvn::CVNProtoDUNEUtils::GetRecoSliceHits	(	const recob::Slice &	slice,
		art::Event const &	evt,
		const std::string	sliceModule
	)		const

Definition at line 17 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                      {

  return GetRecoSliceHits(slice.ID(),evt,sliceModule);  

}

const std::vector< const recob::Hit * > cvn::CVNProtoDUNEUtils::GetRecoSliceHits	(	unsigned int	sliceID,
		art::Event const &	evt,
		const std::string	sliceModule
	)		const

Definition at line 24 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                         {

  auto recoSlices = evt.getValidHandle<std::vector<recob::Slice> >(sliceModule);
  art::FindManyP<recob::Hit> findHits(recoSlices,evt,sliceModule);
  std::vector<art::Ptr<recob::Hit>> inputHits = findHits.at(sliceID);

  std::vector<const recob::Hit*> sliceHits;

  for(const art::Ptr<recob::Hit> hit : inputHits){

    sliceHits.push_back(hit.get());

  }

  return sliceHits;

}

bool cvn::CVNProtoDUNEUtils::IsBeamParticle	(	const recob::PFParticle &	particle,
		art::Event const &	evt,
		const std::string	particleLabel
	)		const

Use the pandora metadata to tell us if this is a beam particle or not.

Definition at line 181 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                  {
  std::map<std::string,float> mdMap = GetPFParticleMetaData(particle,evt,particleLabel);
  if(mdMap.find("IsTestBeam") != mdMap.end()){
    return true;
  }
  else{
    return false;
  }
}

const std::map< unsigned int, std::vector< const recob::PFParticle * > > cvn::CVNProtoDUNEUtils::SliceMapHelper ( art::Event const &  evt, const std::string  particleLabel, bool  primaryOnly  ) const

private

Helper to get the slice map and avoid code repetition.

Definition at line 93 of file CVNProtoDUNEUtils.cxx.

                                                                                                                                                                           {

  // Get the particles
  auto pfParticles = evt.getValidHandle<std::vector<recob::PFParticle>>(particleLabel);

  std::map<unsigned int, std::vector<const recob::PFParticle*>> sliceMap;

  for(unsigned int p = 0; p < pfParticles->size(); ++p){
    const recob::PFParticle* particle = &(pfParticles->at(p));

    //  Only the primary particles have the slice association
    if(primaryOnly && !particle->IsPrimary()) continue;

    unsigned int thisSlice = GetPFParticleSliceIndex(*particle,evt,particleLabel);

    if(thisSlice != 9999){
      sliceMap[thisSlice].push_back(particle);
    }
  }

  return sliceMap;

}

---

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

• dunereco/dunereco/CVN/func/CVNProtoDUNEUtils.h
• dunereco/dunereco/CVN/func/CVNProtoDUNEUtils.cxx