dt::Chunk Class Reference

Sequence of contiguous hits and dead cells all on the same plane. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/DiscreteTracker/Chunk.h"

Public Member Functions
	Chunk ()
	Chunk (int plane, int cellInLo, int cellInHi, int cellOutLo, int cellOutHi, bool up, art::PtrVector< rb::CellHit > hits)
int	Plane () const
int	EntryHitCell () const
int	ExitHitCell () const
int	EntryExtremeCell () const
int	ExitExtremeCell () const
int	TopHitCell () const
int	BottomHitCell () const
int	TopExtremeCell () const
int	BottomExtremeCell () const
bool	Up () const
void	SetUp (bool up)
bool	Matches (const Chunk &ch) const
bool	HasHits () const
bool	AllDead () const
art::PtrVector< rb::CellHit >	AllHits () const
art::PtrVector< rb::CellHit >	HitsOnLine (double grad, double icept) const
std::pair< Segment, Segment >	GetSegs () const
bool	IsWorthShowering () const
bool	IsShowerChunk () const
void	SetShowerChunk ()
bool	operator< (const Chunk &rhs) const

Protected Member Functions
void	CellPoints (const geo::PlaneGeo *plane, int cell, double &v0, double &v1, double &v2, double &v3, double &z0, double &z1) const
void	CellMeasurements (geo::View_t view, const geo::CellGeo *cell, double &v0, double &dv, double &z0, double &dz) const

Protected Attributes
int	fPlane
int	fEntryCellLo
int	fEntryCellHi
int	fExitCellLo
int	fExitCellHi
bool	fUp
bool	fShowerChunk
	"Shower chunk" can enter and exit anywhere More...
bool	fSegsDirty
std::pair< Segment, Segment >	fSegs
art::PtrVector< rb::CellHit >	fHits

Detailed Description

Sequence of contiguous hits and dead cells all on the same plane.

Definition at line 17 of file Chunk.h.

Constructor & Destructor Documentation

dt::Chunk::Chunk ( )

inline

Definition at line 20 of file Chunk.h.

References hits(), and NDAPDHVSetting::plane.

{}

dt::Chunk::Chunk	(	int	plane,
		int	cellInLo,
		int	cellInHi,
		int	cellOutLo,
		int	cellOutHi,
		bool	up,
		art::PtrVector< rb::CellHit >	hits
	)

Definition at line 21 of file Chunk.cxx.

References ana::assert().

    : fPlane(plane),
      fEntryCellLo(cellInLo), fEntryCellHi(cellInHi),
      fExitCellLo(cellOutLo), fExitCellHi(cellOutHi),
      fUp(up), fShowerChunk(false),
      fSegsDirty(true),
      fHits(hits)
  {
    assert(cellInLo <= cellInHi);
    assert(cellOutLo <= cellOutHi);
  }

Member Function Documentation

bool dt::Chunk::AllDead ( ) const

inline

Definition at line 36 of file Chunk.h.

References art::PtrVector< T >::empty(), and fHits.

{return  fHits.empty();}

art::PtrVector<rb::CellHit> dt::Chunk::AllHits ( ) const

inline

Definition at line 37 of file Chunk.h.

References fHits, GetSegs(), stan::math::grad(), HitsOnLine(), and IsWorthShowering().

Referenced by dt::DiscreteTracker::BestSeed().

{return fHits;}

int dt::Chunk::BottomExtremeCell ( ) const

inline

Definition at line 31 of file Chunk.h.

References fEntryCellLo, fExitCellLo, fUp, Matches(), SetUp(), and Up().

Referenced by GetSegs(), and Matches().

{return fUp ? fEntryCellLo : fExitCellLo; }

int dt::Chunk::BottomHitCell ( ) const

inline

Definition at line 29 of file Chunk.h.

References fEntryCellHi, fExitCellHi, and fUp.

Referenced by GetSegs(), Matches(), and dt::Chain::ToTrack().

{return fUp ? fEntryCellHi : fExitCellHi; }

void dt::Chunk::CellMeasurements ( geo::View_t  view, const geo::CellGeo *  cell, double &  v0, double &  dv, double &  z0, double &  dz  ) const

protected

Parameters

v0	cell centre
v1	cell half-width

Definition at line 245 of file Chunk.cxx.

References geo::CellGeo::GetCenter(), geo::CellGeo::HalfD(), geo::CellGeo::HalfW(), and POTSpillRate::view.

Referenced by CellPoints(), and operator<().

  {
    double xyz[3];
    cell->GetCenter(xyz);
    v0 = xyz[view];
    dv = cell->HalfW();
    z0 = xyz[2];
    dz = cell->HalfD();
  }

void dt::Chunk::CellPoints ( const geo::PlaneGeo *  plane, int  cell, double &  v0, double &  v1, double &  v2, double &  v3, double &  z0, double &  z1  ) const

protected

Parameters

v0	lowest point not in a neighbouring cell
v1	lowest point in this cell
v2	highest point in this cell
v3	highest point not in a neighbouring cell

Definition at line 203 of file Chunk.cxx.

References geo::PlaneGeo::Cell(), CellMeasurements(), dz, POTSpillRate::view, and geo::PlaneGeo::View().

Referenced by GetSegs(), and operator<().

  {
    const geo::View_t view = plane->View();

    // Measurements of the cell itself
    double vc, dv, zc, dz;
    CellMeasurements(view, plane->Cell(cell), vc, dv, zc, dz);
    v1 = vc-dv;
    v2 = vc+dv;
    z0 = zc-dz;
    z1 = zc+dz;

    // Would handle bad channels by walking these until the channel they're on
    // is good.
    const geo::CellGeo* cellabove = plane->Cell(cell+1);
    const geo::CellGeo* cellbelow = plane->Cell(cell-1);

    double junk;
    if(cellbelow){
      // Lowest point not in another cell
      CellMeasurements(view, cellbelow, v0, dv, junk, junk);
      v0 += dv;
    }
    else{
      // No cell below, must be edge of the detector, could be anywhere low
      v0 = -1e3;
    }

    if(cellabove){
      // Highest point not in another cell
      CellMeasurements(view, cellabove, v3, dv, junk, junk);
      v3 -= dv;
    }
    else{
      // No cell above, must be edge of the detector, could be anywhere high
      v3 = +1e3;
    }
  }

int dt::Chunk::EntryExtremeCell ( ) const

inline

Definition at line 26 of file Chunk.h.

References fEntryCellHi, fEntryCellLo, and fUp.

{return fUp ? fEntryCellLo : fEntryCellHi;}

int dt::Chunk::EntryHitCell ( ) const

inline

Definition at line 24 of file Chunk.h.

References fEntryCellHi, fEntryCellLo, and fUp.

{return fUp ? fEntryCellHi : fEntryCellLo;}

int dt::Chunk::ExitExtremeCell ( ) const

inline

Definition at line 27 of file Chunk.h.

References fExitCellHi, fExitCellLo, and fUp.

{return fUp ? fExitCellHi  : fExitCellLo; }

int dt::Chunk::ExitHitCell ( ) const

inline

Definition at line 25 of file Chunk.h.

References fExitCellHi, fExitCellLo, and fUp.

{return fUp ? fExitCellLo  : fExitCellHi; }

std::pair< Segment, Segment > dt::Chunk::GetSegs

(

)

const

Definition at line 102 of file Chunk.cxx.

References BottomExtremeCell(), BottomHitCell(), CellPoints(), art::PtrVector< T >::empty(), fEntryCellHi, fEntryCellLo, fExitCellHi, fExitCellLo, fHits, fPlane, fSegs, fSegsDirty, fShowerChunk, fUp, geom(), geo::GeometryBase::Plane(), NDAPDHVSetting::plane, TopExtremeCell(), and TopHitCell().

Referenced by AllHits(), dt::DiscreteTracker::BestSeed(), dt::Cand::MaybeFlip(), and dt::DiscreteTracker::PossibleNextSegs().

  {
    if(!fSegsDirty) return fSegs;
    fSegsDirty = false;

    art::ServiceHandle<geo::Geometry> geom;

    const geo::PlaneGeo* plane = geom->Plane(fPlane);

    double v0lo, v1lo, v2lo, v3lo, z1lo, z2lo;
    CellPoints(plane, BottomExtremeCell(),
               v0lo, v1lo, v2lo, v3lo,
               z1lo, z2lo);

    double v0hi, v1hi, v2hi, v3hi, z1hi, z2hi;
    CellPoints(plane, TopExtremeCell(),
               v0hi, v1hi, v2hi, v3hi,
               z1hi, z2hi);

    double junk;
    if(BottomHitCell() != BottomExtremeCell()){
      CellPoints(plane, BottomHitCell(),
                 junk, junk, v2lo, v3lo,
                 junk, junk);
    }
    if(TopHitCell() != TopExtremeCell()){
      CellPoints(plane, TopHitCell(),
                 v0hi, v1hi, junk, junk,
                 junk, junk);
    }

    // If this is a completely dead chunk, it's OK to enter in the crack part
    // on one side and exit in that same part without ever traversing any real
    // cell.
    if(fHits.empty()){
      v1lo = v0lo;
      v1hi = v0hi;
      v2lo = v3lo;
      v2hi = v3hi;
    }

    // If the track is upgoing it can enter the plane in plastic below the
    // lowest cell, and exit in the plastic above the highest cell. But it
    // mustn't enter higher than the scintillator of the first cell, or lower
    // than the scintillator of the last cell. Vice versa for downgoing.

    if(fShowerChunk){
      if(fUp){
        fSegs.first = Segment(z1lo, v0lo-kCornerR, v2hi,
                              z1lo+kCornerR, v0lo, v3hi,
                              true, fPlane, fEntryCellLo, fEntryCellHi);
        fSegs.second = Segment(z2hi-kCornerR, v0lo, v3hi,
                               z2hi, v1lo, v3hi+kCornerR,
                               false, fPlane, fExitCellLo, fExitCellHi);
      }
      else{
        fSegs.first = Segment(z1hi, v1lo, v3hi+kCornerR,
                              z1hi+kCornerR, v0lo, v3hi,
                              true, fPlane, fEntryCellLo, fEntryCellHi);
        fSegs.second = Segment(z2lo-kCornerR, v0lo, v3hi,
                               z2lo, v0lo-kCornerR, v2hi,
                               false, fPlane, fExitCellLo, fExitCellHi);
      }

      return fSegs;
    }

    if(fUp){
      // bottom left
      fSegs.first = Segment(z1lo, v0lo-kCornerR, v2lo,
                            z1lo+kCornerR, v0lo, v3hi,
                            true, fPlane, fEntryCellLo, fEntryCellHi);
      // top right
      fSegs.second = Segment(z2hi-kCornerR, v0lo, v3hi,
                             z2hi, v1hi, v3hi+kCornerR,
                             false, fPlane, fExitCellLo, fExitCellHi);
    }
    else{
      // top left
      fSegs.first = Segment(z1hi, v1hi, v3hi+kCornerR,
                            z1hi+kCornerR, v0lo, v3hi,
                            true, fPlane, fEntryCellLo, fEntryCellHi);
      // bottom right
      fSegs.second = Segment(z2lo-kCornerR, v0lo, v3hi,
                             z2lo, v0lo-kCornerR, v2lo,
                             false, fPlane, fExitCellLo, fExitCellHi);
    }

    return fSegs;
  }

bool dt::Chunk::HasHits ( ) const

inline

Definition at line 35 of file Chunk.h.

References art::PtrVector< T >::empty(), and fHits.

Referenced by dt::Cand::AddChunk().

{return !fHits.empty();}

art::PtrVector< rb::CellHit > dt::Chunk::HitsOnLine	(	double	grad,
		double	icept
	)		const

Definition at line 67 of file Chunk.cxx.

References rb::CellHit::Cell(), geo::PlaneGeo::Cell(), dz, fHits, fPlane, geom(), geo::CellGeo::GetCenter(), geo::CellGeo::HalfD(), geo::CellGeo::HalfW(), getGoodRuns4SAM::n, geo::GeometryBase::Plane(), art::PtrVector< T >::push_back(), runNovaSAM::ret, art::PtrVector< T >::size(), std::swap(), and rb::CellHit::View().

Referenced by AllHits().

  {
    art::ServiceHandle<geo::Geometry> geom;

    const geo::PlaneGeo* geoplane = geom->Plane(fPlane);

    art::PtrVector<rb::CellHit> ret;

    for(unsigned int n = 0; n < fHits.size(); ++n){
      const art::Ptr<rb::CellHit>& hit = fHits[n];

      double xyz[3];
      const geo::CellGeo* geocell = geoplane->Cell(hit->Cell());
      geocell->GetCenter(xyz);
      const double z0 = xyz[2];
      const double dz = geocell->HalfD();
      const double z1 = z0-dz;
      const double z2 = z0+dz;
      const double v0 = xyz[hit->View()];
      const double dv = geocell->HalfW();
      const double vlo = v0-dv;
      const double vhi = v0+dv;

      double v1 = icept+grad*z1;
      double v2 = icept+grad*z2;
      if(v1 > v2) std::swap(v1, v2);

      if(v1 <= vhi && v2 >= vlo) ret.push_back(hit);
    } // end for n

    return ret;
  }

bool dt::Chunk::IsShowerChunk ( ) const

inline

Definition at line 43 of file Chunk.h.

References fShowerChunk.

Referenced by dt::DiscreteTracker::TryExtendCand().

{return fShowerChunk;}

bool dt::Chunk::IsWorthShowering

(

)

const

Definition at line 194 of file Chunk.cxx.

References fEntryCellHi, fEntryCellLo, fExitCellHi, and fExitCellLo.

Referenced by AllHits(), dt::DiscreteTracker::PossibleNextSegs(), and dt::DiscreteTracker::TryExtendCand().

  {
    // Single hit
    if(fEntryCellLo == fExitCellHi &&
       fEntryCellHi == fExitCellLo) return false;
    return true;
  }

bool dt::Chunk::Matches

(

const Chunk &

ch

)

const

Definition at line 35 of file Chunk.cxx.

References BottomExtremeCell(), BottomHitCell(), fPlane, TopExtremeCell(), and TopHitCell().

Referenced by BottomExtremeCell().

  {
    if(fPlane != ch.fPlane) return false;

    if(BottomExtremeCell() != ch.BottomExtremeCell()) return false;
    if(TopExtremeCell() != ch.TopExtremeCell()) return false;

    if(BottomHitCell() != ch.BottomHitCell()) return false;
    if(TopHitCell() != ch.TopHitCell()) return false;

    return true;
  }

bool dt::Chunk::operator< ( const Chunk &  rhs ) const

inline

Definition at line 46 of file Chunk.h.

References getBrightness::cell, CellMeasurements(), CellPoints(), dz, fPlane, and POTSpillRate::view.

{return fPlane < rhs.fPlane;}

int dt::Chunk::Plane ( ) const

inline

Definition at line 23 of file Chunk.h.

References fPlane.

Referenced by dt::Cand::AllHits(), dt::Cand::EstimateStraightLine(), dt::Chain::FirstPlane(), dt::Chain::LastPlane(), dt::Cand::operator<(), and dt::Chain::ToTrack().

{return fPlane;}

void dt::Chunk::SetShowerChunk ( )

inline

Definition at line 44 of file Chunk.h.

References fSegsDirty, and fShowerChunk.

Referenced by dt::DiscreteTracker::PossibleNextSegs().

{fSegsDirty = true; fShowerChunk = true;}

void dt::Chunk::SetUp

(

bool

up

)

Definition at line 55 of file Chunk.cxx.

References fEntryCellHi, fEntryCellLo, fExitCellHi, fExitCellLo, fSegsDirty, fUp, and std::swap().

Referenced by BottomExtremeCell(), dt::Cand::MaybeFlip(), and dt::DiscreteTracker::PossibleNextSegs().

  {
    fSegsDirty = true;

    if(up != fUp){
      std::swap(fEntryCellLo, fExitCellLo);
      std::swap(fEntryCellHi, fExitCellHi);
    }
    fUp = up;
  }

int dt::Chunk::TopExtremeCell ( ) const

inline

Definition at line 30 of file Chunk.h.

References fEntryCellHi, fExitCellHi, and fUp.

Referenced by GetSegs(), and Matches().

{return fUp ? fExitCellHi  : fEntryCellHi;}

int dt::Chunk::TopHitCell ( ) const

inline

Definition at line 28 of file Chunk.h.

References fEntryCellLo, fExitCellLo, and fUp.

Referenced by GetSegs(), Matches(), and dt::Chain::ToTrack().

{return fUp ? fExitCellLo  : fEntryCellLo;}

bool dt::Chunk::Up

(

)

const

Definition at line 49 of file Chunk.cxx.

References fUp.

Referenced by dt::Cand::AddChunk(), dt::DiscreteTracker::BestSeed(), BottomExtremeCell(), and dt::DiscreteTracker::PossibleNextSegs().

  {
    return fUp;
  }

Member Data Documentation

int dt::Chunk::fEntryCellHi

protected

Definition at line 61 of file Chunk.h.

Referenced by BottomHitCell(), EntryExtremeCell(), EntryHitCell(), GetSegs(), IsWorthShowering(), SetUp(), and TopExtremeCell().

int dt::Chunk::fEntryCellLo

protected

Definition at line 61 of file Chunk.h.

Referenced by BottomExtremeCell(), EntryExtremeCell(), EntryHitCell(), GetSegs(), IsWorthShowering(), SetUp(), and TopHitCell().

int dt::Chunk::fExitCellHi

protected

Definition at line 61 of file Chunk.h.

Referenced by BottomHitCell(), ExitExtremeCell(), ExitHitCell(), GetSegs(), IsWorthShowering(), SetUp(), and TopExtremeCell().

int dt::Chunk::fExitCellLo

protected

Definition at line 61 of file Chunk.h.

Referenced by BottomExtremeCell(), ExitExtremeCell(), ExitHitCell(), GetSegs(), IsWorthShowering(), SetUp(), and TopHitCell().

art::PtrVector<rb::CellHit> dt::Chunk::fHits

protected

Definition at line 68 of file Chunk.h.

Referenced by AllDead(), AllHits(), GetSegs(), HasHits(), and HitsOnLine().

int dt::Chunk::fPlane

protected

Definition at line 60 of file Chunk.h.

Referenced by GetSegs(), HitsOnLine(), Matches(), operator<(), and Plane().

std::pair<Segment, Segment> dt::Chunk::fSegs

mutableprotected

Definition at line 66 of file Chunk.h.

Referenced by GetSegs().

bool dt::Chunk::fSegsDirty

mutableprotected

Definition at line 65 of file Chunk.h.

Referenced by GetSegs(), SetShowerChunk(), and SetUp().

bool dt::Chunk::fShowerChunk

protected

"Shower chunk" can enter and exit anywhere

Definition at line 63 of file Chunk.h.

Referenced by GetSegs(), IsShowerChunk(), and SetShowerChunk().

bool dt::Chunk::fUp

protected

Definition at line 62 of file Chunk.h.

Referenced by BottomExtremeCell(), BottomHitCell(), EntryExtremeCell(), EntryHitCell(), ExitExtremeCell(), ExitHitCell(), GetSegs(), SetUp(), TopExtremeCell(), TopHitCell(), and Up().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/DiscreteTracker/Chunk.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/DiscreteTracker/Chunk.cxx