pizero::Cone Class Reference

#include <Cone.h>

Public Member Functions
	Cone (const art::Event &evt, detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, const TVector3 &newStart, const TVector3 &newDir, const double newLen)
TVector3	start () const
TVector3	end () const
TVector3	direction () const
double	length () const
double	width () const
std::vector< const recob::Hit * >	hits () const
std::vector< const recob::SpacePoint * >	spacepoints () const
double	energy (calo::CalorimetryAlg caloAlg) const
double	completeness (const simb::MCParticle &mcpart, const std::string hitLabel="hitpdune") const
double	purity (const simb::MCParticle &mcpart) const

Private Attributes
const art::Event *	m_evt
detinfo::DetectorClocksData const *	m_clockdata
detinfo::DetectorPropertiesData const *	m_detprop
TVector3	m_start
TVector3	m_direction
double	m_length
double	m_ang = 15.0/360 * 2*3.1415926
double	m_width = 1
double	m_startoffset = 10
std::vector< const recob::Hit * >	m_hits
std::vector< const recob::SpacePoint * >	m_spacepoints

Detailed Description

Definition at line 22 of file Cone.h.

Constructor & Destructor Documentation

pizero::Cone::Cone ( const art::Event &  evt, detinfo::DetectorClocksData const &  clockData, detinfo::DetectorPropertiesData const &  detProp, const TVector3 &  newStart, const TVector3 &  newDir, const double  newLen  )

inline

Definition at line 38 of file Cone.h.

                           :
    m_evt(&evt),
    m_clockdata(&clockData),
    m_detprop(&detProp),
    m_start(newStart), m_direction(newDir.Unit()), m_length(newLen)
  {
    // Move the start back by a constant to catch all hits.
    m_start = m_start - m_direction.Unit()*m_startoffset;
    m_length += m_startoffset;
    // Determine the width of the cone at the end.
    m_width = m_length * tan(m_ang);

    // Get all spacepoints in the event.
    std::string spsLabel = "pandora";
    auto spHandle = m_evt->getHandle<std::vector<recob::SpacePoint> >(spsLabel);
    if (!spHandle) { return; }
    // Get associations between the spacepoints and hits.
    const art::FindManyP<recob::Hit> sphits(spHandle, *m_evt, spsLabel);

    // Loop through spacepoints to test whether they're in the cone.
    for(const recob::SpacePoint& sp : *spHandle) {
      const TVector3 pos(sp.XYZ());
      // 0 <= projection <= cone length
      const double proj = (pos-m_start).Dot(m_direction);
      if(proj < 0 || proj > m_length) continue;
      // Cone radius at projection.
      const double cradius = (proj/m_length)*m_width/2;
      // rejection < cone radius
      const double rej = ((pos-m_start) - proj*m_direction).Mag();
      if(rej > cradius) continue;

      // This spacepoint is inside the cone. Put it among the results.
      for(const art::Ptr<recob::Hit>& hit : sphits.at(sp.ID())) {
        m_hits.push_back(&*hit);
      }
      m_spacepoints.push_back(&sp);
    }
  }

Member Function Documentation

double pizero::Cone::completeness ( const simb::MCParticle &  mcpart, const std::string  hitLabel = "hitpdune"  ) const

inline

Definition at line 100 of file Cone.h.

                                                          {
  // Get all hits in the event.
  auto hitHandle = m_evt->getHandle<std::vector<recob::Hit> >(hitLabel);
  if(!hitHandle) {
    std::cout << "pizero::Cone::completeness: could not find hits in event.\n";
    return 0;
  }

  // Check for each hit whether it came from the MCParticle and if so,
  // whether it fell inside the cone.
  art::ServiceHandle<cheat::BackTrackerService> bt_serv;
  double sharedHits = 0;
  double MChits = 0;
  for(const recob::Hit& hit : *hitHandle) {
    for(const int trackId : bt_serv->HitToTrackIds(*m_clockdata, hit)) {
      if(std::abs(trackId) == std::abs(mcpart.TrackId())) {
        // Hit is in MCParticle.
        ++MChits;
        // Check if hit in cone.
        for(const recob::Hit* chit : m_hits) {
          // Compare the time and place of hit for lack of ID
          if(chit->PeakTime() - hit.PeakTime() < 1e-5 &&
             chit->Channel() == hit.Channel()) {
            ++sharedHits;
            // std::cout << "Found shared hit!\n";
            break;
          }
        }
        break;
      } // if hit is from MCP
    } // for trackID in hit
  } // for hit in event

  return MChits==0? 0: sharedHits/MChits;
} // Cone::completeness

TVector3 pizero::Cone::direction ( ) const

inline

Definition at line 85 of file Cone.h.

{ return m_direction; }

TVector3 pizero::Cone::end ( void  ) const

inline

Definition at line 84 of file Cone.h.

{ return m_start + m_direction*m_length; }

double pizero::Cone::energy ( calo::CalorimetryAlg  caloAlg ) const

inline

Definition at line 91 of file Cone.h.

                                                  {
    protoana::ProtoDUNEShowerUtils shUtil;
    return shUtil.EstimateEnergyFromHitCharge(*m_clockdata, *m_detprop, m_hits, caloAlg)[2];
  }

std::vector<const recob::Hit*> pizero::Cone::hits ( void  ) const

inline

Definition at line 88 of file Cone.h.

{ return m_hits; }

double pizero::Cone::length ( ) const

inline

Definition at line 86 of file Cone.h.

{ return m_length; }

double pizero::Cone::purity ( const simb::MCParticle &  mcpart ) const

inline

Definition at line 137 of file Cone.h.

                                                             {
  unsigned sharedHits = 0;
  art::ServiceHandle<cheat::BackTrackerService> bt_serv;
  for(const recob::Hit* hit : m_hits) {
    for(const int trackId : bt_serv->HitToTrackIds(*m_clockdata, *hit)) {
      if(std::abs(trackId) == std::abs(mcpart.TrackId())) {
        ++sharedHits;
        break;
      }
    }
  }

  return (double)sharedHits / m_hits.size();
}

std::vector<const recob::SpacePoint*> pizero::Cone::spacepoints ( ) const

inline

Definition at line 89 of file Cone.h.

{ return m_spacepoints; }

TVector3 pizero::Cone::start ( void  ) const

inline

Definition at line 83 of file Cone.h.

{ return m_start; }

double pizero::Cone::width ( ) const

inline

Definition at line 87 of file Cone.h.

{ return m_width; }

Member Data Documentation

double pizero::Cone::m_ang = 15.0/360 * 2*3.1415926

private

Definition at line 30 of file Cone.h.

detinfo::DetectorClocksData const* pizero::Cone::m_clockdata

private

Definition at line 25 of file Cone.h.

detinfo::DetectorPropertiesData const* pizero::Cone::m_detprop

private

Definition at line 26 of file Cone.h.

TVector3 pizero::Cone::m_direction

private

Definition at line 28 of file Cone.h.

const art::Event* pizero::Cone::m_evt

private

Definition at line 24 of file Cone.h.

std::vector<const recob::Hit*> pizero::Cone::m_hits

private

Definition at line 34 of file Cone.h.

double pizero::Cone::m_length

private

Definition at line 29 of file Cone.h.

std::vector<const recob::SpacePoint*> pizero::Cone::m_spacepoints

private

Definition at line 35 of file Cone.h.

TVector3 pizero::Cone::m_start

private

Definition at line 27 of file Cone.h.

double pizero::Cone::m_startoffset = 10

private

Definition at line 32 of file Cone.h.

double pizero::Cone::m_width = 1

private

Definition at line 31 of file Cone.h.

---

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

• protoduneana/protoduneana/singlephase/EMTaskForce/Cone.h