ShowerRecoTools::ShowerLinearEnergy Class Reference

Inheritance diagram for ShowerRecoTools::ShowerLinearEnergy:

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

Private Member Functions
double	CalculateEnergy (const detinfo::DetectorClocksData &clockData, const detinfo::DetectorPropertiesData &detProp, const std::vector< art::Ptr< recob::Hit >> &hits, const geo::PlaneID::PlaneID_t plane) const
Private 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)
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
unsigned int	fNumPlanes
std::vector< double >	fGradients
std::vector< double >	fIntercepts
art::InputTag	fPFParticleLabel
int	fVerbose
std::string	fShowerEnergyOutputLabel
std::string	fShowerBestPlaneOutputLabel
art::ServiceHandle< geo::Geometry >	fGeom

Detailed Description

Definition at line 23 of file ShowerLinearEnergy_tool.cc.

Constructor & Destructor Documentation

ShowerRecoTools::ShowerLinearEnergy::ShowerLinearEnergy

(

const fhicl::ParameterSet &

pset

)

Definition at line 54 of file ShowerLinearEnergy_tool.cc.

    : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
    , fGradients(pset.get<std::vector<double>>("Gradients"))
    , fIntercepts(pset.get<std::vector<double>>("Intercepts"))
    , fPFParticleLabel(pset.get<art::InputTag>("PFParticleLabel"))
    , fVerbose(pset.get<int>("Verbose"))
    , fShowerEnergyOutputLabel(pset.get<std::string>("ShowerEnergyOutputLabel"))
    , fShowerBestPlaneOutputLabel(pset.get<std::string>("ShowerBestPlaneOutputLabel"))
  {
    fNumPlanes = fGeom->Nplanes();
    if (fNumPlanes != fGradients.size() || fNumPlanes != fIntercepts.size()) {
      throw cet::exception("ShowerLinearEnergy")
        << "The number of planes does not match the size of the fcl parametes passed: Num Planes: "
        << fNumPlanes << ", Gradients size: " << fGradients.size()
        << ", Intercpts size: " << fIntercepts.size();
    }
  }

Member Function Documentation

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

overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 73 of file ShowerLinearEnergy_tool.cc.

  {

    // Get the assocated pfParicle vertex PFParticles
    auto const pfpHandle = Event.getValidHandle<std::vector<recob::PFParticle>>(fPFParticleLabel);

    //Get the clusters
    auto const clusHandle = Event.getValidHandle<std::vector<recob::Cluster>>(fPFParticleLabel);

    const art::FindManyP<recob::Cluster>& fmc =
      ShowerEleHolder.GetFindManyP<recob::Cluster>(pfpHandle, Event, fPFParticleLabel);
    // art::FindManyP<recob::Cluster> fmc(pfpHandle, Event, fPFParticleLabel);
    std::vector<art::Ptr<recob::Cluster>> clusters = fmc.at(pfparticle.key());

    //Get the hit association
    const art::FindManyP<recob::Hit>& fmhc =
      ShowerEleHolder.GetFindManyP<recob::Hit>(clusHandle, Event, fPFParticleLabel);
    // art::FindManyP<recob::Hit> fmhc(clusHandle, Event, fPFParticleLabel);

    std::map<geo::PlaneID::PlaneID_t, std::vector<art::Ptr<recob::Hit>>> planeHits;

    //Loop over the clusters in the plane and get the hits
    for (auto const& cluster : clusters) {

      //Get the hits
      std::vector<art::Ptr<recob::Hit>> hits = fmhc.at(cluster.key());

      //Get the plane.
      const geo::PlaneID::PlaneID_t plane(cluster->Plane().Plane);

      planeHits[plane].insert(planeHits[plane].end(), hits.begin(), hits.end());
    }

    // Calculate the energy for each plane && best plane
    geo::PlaneID::PlaneID_t bestPlane = std::numeric_limits<geo::PlaneID::PlaneID_t>::max();
    unsigned int bestPlaneNumHits = 0;

    //Holder for the final product
    std::vector<double> energyVec(fNumPlanes, -999.);
    std::vector<double> energyError(fNumPlanes, -999.);

    auto const clockData =
      art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(Event);
    auto const detProp =
      art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(Event, clockData);

    for (auto const& [plane, hits] : planeHits) {

      unsigned int planeNumHits = hits.size();

      //Calculate the Energy for
      double Energy = CalculateEnergy(clockData, detProp, hits, plane);
      // If the energy is negative, leave it at -999
      if (Energy > 0) energyVec.at(plane) = Energy;

      if (planeNumHits > bestPlaneNumHits) {
        bestPlane = plane;
        bestPlaneNumHits = planeNumHits;
      }
    }

    ShowerEleHolder.SetElement(energyVec, energyError, fShowerEnergyOutputLabel);
    // Only set the best plane if it has some hits in it
    if (bestPlane < fGeom->Nplanes()) {
      // Need to cast as an int for legacy default of -999
      // have to define a new variable as we pass-by-reference when filling
      int bestPlaneVal(bestPlane);
      ShowerEleHolder.SetElement(bestPlaneVal, fShowerBestPlaneOutputLabel);
    }

    return 0;
  }

double ShowerRecoTools::ShowerLinearEnergy::CalculateEnergy ( const detinfo::DetectorClocksData &  clockData, const detinfo::DetectorPropertiesData &  detProp, const std::vector< art::Ptr< recob::Hit >> &  hits, const geo::PlaneID::PlaneID_t  plane  ) const

private

Definition at line 151 of file ShowerLinearEnergy_tool.cc.

  {

    double totalCharge = 0, totalEnergy = 0;

    for (auto const& hit : hits) {
      totalCharge += (hit->Integral() * std::exp((sampling_rate(clockData) * hit->PeakTime()) /
                                                 (detProp.ElectronLifetime() * 1e3)));
    }

    totalEnergy = (totalCharge * fGradients.at(plane)) + fIntercepts.at(plane);

    return totalEnergy;
  }

Member Data Documentation

art::ServiceHandle<geo::Geometry> ShowerRecoTools::ShowerLinearEnergy::fGeom

private

Definition at line 51 of file ShowerLinearEnergy_tool.cc.

std::vector<double> ShowerRecoTools::ShowerLinearEnergy::fGradients

private

Definition at line 41 of file ShowerLinearEnergy_tool.cc.

std::vector<double> ShowerRecoTools::ShowerLinearEnergy::fIntercepts

private

Definition at line 42 of file ShowerLinearEnergy_tool.cc.

unsigned int ShowerRecoTools::ShowerLinearEnergy::fNumPlanes

private

Definition at line 40 of file ShowerLinearEnergy_tool.cc.

art::InputTag ShowerRecoTools::ShowerLinearEnergy::fPFParticleLabel

private

Definition at line 44 of file ShowerLinearEnergy_tool.cc.

std::string ShowerRecoTools::ShowerLinearEnergy::fShowerBestPlaneOutputLabel

private

Definition at line 48 of file ShowerLinearEnergy_tool.cc.

std::string ShowerRecoTools::ShowerLinearEnergy::fShowerEnergyOutputLabel

private

Definition at line 47 of file ShowerLinearEnergy_tool.cc.

int ShowerRecoTools::ShowerLinearEnergy::fVerbose

private

Definition at line 45 of file ShowerLinearEnergy_tool.cc.

---

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

• ShowerLinearEnergy_tool.cc