MMT_Road Class Reference

#include <MMT_Road.h>

Collaboration diagram for MMT_Road:

Public Member Functions
	MMT_Road (const char sector, const int roadSize, const int UpX, const int DownX, const int UpUV, const int DownUV, const int xthr, const int uvthr, const int iroadx, const int iroadu=-1, const int iroadv=-1)
	~MMT_Road ()=default
void	addHits (std::vector< std::shared_ptr< MMT_Hit > > &hits)
double	avgSofX () const
double	avgSofUV (const int uv1, const int uv2) const
double	avgZofUV (const int uv1, const int uv2) const
bool	checkCoincidences (const int bcwind) const
unsigned int	countHits () const
unsigned int	countRealHits () const
unsigned int	countUVHits (bool flag) const
unsigned int	countXHits (bool flag) const
bool	evaluateLowRes () const
bool	horizontalCheck () const
void	incrementAge (const int bcwind)
const std::vector< std::unique_ptr< MMT_Hit > > &	getHitVector () const
int	getRoadSize () const
int	getRoadSizeUpX () const
int	getRoadSizeDownX () const
int	getRoadSizeUpUV () const
int	getRoadSizeDownUV () const
char	getSector () const
int	getXthreshold () const
int	getUVthreshold () const
int	iRoadx () const
int	iRoadu () const
int	iRoadv () const
bool	matureCheck (const int bcwind) const
double	mxl () const
void	reset ()
bool	stereoCheck () const

Private Attributes
int	m_iroadx
int	m_iroadu
int	m_iroadv
char	m_sector
int	m_xthr
int	m_uvthr
int	m_roadSize
int	m_roadSizeUpX
int	m_roadSizeDownX
int	m_roadSizeUpUV
int	m_roadSizeDownUV
std::vector< std::unique_ptr< MMT_Hit > >	m_road_hits

Detailed Description

Definition at line 12 of file MMT_Road.h.

Constructor & Destructor Documentation

MMT_Road()

MMT_Road::MMT_Road	(	const char	sector,
		const int	roadSize,
		const int	UpX,
		const int	DownX,
		const int	UpUV,
		const int	DownUV,
		const int	xthr,
		const int	uvthr,
		const int	iroadx,
		const int	iroadu = -1,
		const int	iroadv = -1
	)

Definition at line 7 of file MMT_Road.cxx.

                                                                         {
  m_sector = sector;
  m_iroadx = iroadx;
  m_iroadu = (iroadu != -1) ? iroadu : iroadx;
  m_iroadv = (iroadv != -1) ? iroadv : iroadx;
  m_xthr = xthr;
  m_uvthr = uvthr;
 
  m_roadSize = roadSize;
  m_roadSizeUpX = UpX;
  m_roadSizeDownX = DownX;
  m_roadSizeUpUV = UpUV;
  m_roadSizeDownUV = DownUV;
}

~MMT_Road()

MMT_Road::~MMT_Road ( )

default

Member Function Documentation

addHits()

void MMT_Road::addHits

(

std::vector< std::shared_ptr< MMT_Hit > > &

hits

)

Definition at line 23 of file MMT_Road.cxx.

                                                               {
  for (const auto &hit_i : hits) {
    if (m_sector != hit_i->getSector()) continue;
 
    int iroad = 0;
    unsigned short int olow = 0, ohigh = 0;
    if (hit_i->isX()) {
      iroad = m_iroadx;
      olow  = m_roadSizeDownX;
      ohigh = m_roadSizeUpX;
    }
    else if (hit_i->isU()) {
      iroad = m_iroadu;
      olow  = m_roadSizeDownUV;
      ohigh = m_roadSizeUpUV;
    }
    else if (hit_i->isV()) {
      iroad = m_iroadv;
      olow  = m_roadSizeDownUV;
      ohigh = m_roadSizeUpUV;
    }
    else continue;
 
    double val = hit_i->getShift();
    double slow  = val + (m_roadSize*iroad     + 0.5 - olow )*hit_i->getPitchOverZ();
    double shigh = val + (m_roadSize*(iroad+1) + 0.5 + ohigh)*hit_i->getPitchOverZ();
 
    val = hit_i->getRZSlope();
    bool has_hit = (val > 0.) ? (val > slow && val < shigh) : (val > shigh && val < slow);
    if (!has_hit) continue;
 
    has_hit = false;
    unsigned short int bo = hit_i->getPlane();
    auto it = std::find_if(m_road_hits.begin(), m_road_hits.end(), [&bo](const auto &hit) { return (hit->getPlane() == bo); });
    if (it != m_road_hits.end()) {
      has_hit = true;
      if (!hit_i->isNoise() && (*it)->isNoise()) {
        m_road_hits.erase(it);
        has_hit = false;
      }
    }
 
    if (has_hit) continue;
    m_road_hits.emplace_back(std::make_unique<MMT_Hit>(hit_i.get()));
    m_road_hits.back()->setAge(0);
  }
}

avgSofUV()

double MMT_Road::avgSofUV	(	const int	uv1,
		const int	uv2
	)		const

Definition at line 81 of file MMT_Road.cxx.

                                                            {
  std::vector<double> sl;
  for (const auto &hit : m_road_hits) {
    int bo = hit->getPlane();
    if (bo == uv1 || bo == uv2) sl.push_back(hit->getRZSlope());
  }
  double avg_uv = std::accumulate(sl.begin(), sl.end(), 0.0)/(double)sl.size();
  return avg_uv;
}

avgSofX()

double MMT_Road::avgSofX

(

)

const

Definition at line 71 of file MMT_Road.cxx.

                               {
  std::vector<double> sl;
  for (const auto &hit : m_road_hits) {
    int bo = hit->getPlane();
    if (bo < 2 || bo > 5) sl.push_back(hit->getRZSlope());
  }
  double avg_x = std::accumulate(sl.begin(), sl.end(), 0.0)/(double)sl.size();
  return avg_x;
}

avgZofUV()

double MMT_Road::avgZofUV	(	const int	uv1,
		const int	uv2
	)		const

Definition at line 91 of file MMT_Road.cxx.

                                                            {
  std::vector<double> zs;
  for (const auto &hit : m_road_hits) {
    int bo = hit->getPlane();
    if (bo == uv1 || bo == uv2) zs.push_back(hit->getZ());
  }
  double avg_z = std::accumulate(zs.begin(), zs.end(), 0.0)/(double)zs.size();
  return avg_z;
}

checkCoincidences()

bool MMT_Road::checkCoincidences

(

const int

bcwind

)

const

Definition at line 101 of file MMT_Road.cxx.

                                                       {
  bool passHorizontalCheck = this->horizontalCheck();
  if (!passHorizontalCheck) return false;
  bool passStereoCheck = this->stereoCheck();
  if (!passStereoCheck) return false;
  bool passMatureCheck = this->matureCheck(bcwind);
  return (passHorizontalCheck && passStereoCheck && passMatureCheck);
}

countHits()

unsigned int MMT_Road::countHits ( ) const

inline

Definition at line 23 of file MMT_Road.h.

{ return m_road_hits.size(); }

countRealHits()

unsigned int MMT_Road::countRealHits

(

)

const

Definition at line 110 of file MMT_Road.cxx.

                                           {
  int nreal = 0;
  for (const auto &hit : m_road_hits) {
    if (hit->isNoise() == false) nreal++;
  }
  return nreal;
}

countUVHits()

unsigned int MMT_Road::countUVHits

(

bool

flag

)

const

Definition at line 118 of file MMT_Road.cxx.

                                                  {
  unsigned int nuv = 0;
  for (const auto &hit : m_road_hits) {
    if (hit->getPlane() == 2 || hit->getPlane() == 4) {
      if (hit->isNoise() == flag) nuv++;
    }
    if (hit->getPlane() == 3 || hit->getPlane() == 5) {
      if (hit->isNoise() == flag) nuv++;
    }
  }
  return nuv;
}

countXHits()

unsigned int MMT_Road::countXHits

(

bool

flag

)

const

Definition at line 131 of file MMT_Road.cxx.

                                                 {
  unsigned int nx = 0;
  for (const auto &hit : m_road_hits) {
    if (hit->getPlane() < 2 || hit->getPlane() > 5) {
      if (hit->isNoise() == flag) nx++;
    }
  }
  return nx;
}

evaluateLowRes()

bool MMT_Road::evaluateLowRes

(

)

const

Definition at line 141 of file MMT_Road.cxx.

                                    {
  unsigned int nhits1 = 0, nhits2 = 0;
  for (const auto &hit : m_road_hits) {
    if (hit->getPlane() < 4 && !hit->isNoise()) nhits1++;
    else if (hit->getPlane() > 3 && !hit->isNoise()) nhits2++;
  }
  return (nhits1 < 4 || nhits2 < 4);
}

getHitVector()

const std::vector<std::unique_ptr<MMT_Hit> >& MMT_Road::getHitVector ( ) const

inline

Definition at line 30 of file MMT_Road.h.

{ return m_road_hits; }

getRoadSize()

int MMT_Road::getRoadSize ( ) const

inline

Definition at line 31 of file MMT_Road.h.

{ return m_roadSize; }

getRoadSizeDownUV()

int MMT_Road::getRoadSizeDownUV ( ) const

inline

Definition at line 35 of file MMT_Road.h.

{ return m_roadSizeDownUV; }

getRoadSizeDownX()

int MMT_Road::getRoadSizeDownX ( ) const

inline

Definition at line 33 of file MMT_Road.h.

{ return m_roadSizeDownX; }

getRoadSizeUpUV()

int MMT_Road::getRoadSizeUpUV ( ) const

inline

Definition at line 34 of file MMT_Road.h.

{ return m_roadSizeUpUV; }

getRoadSizeUpX()

int MMT_Road::getRoadSizeUpX ( ) const

inline

Definition at line 32 of file MMT_Road.h.

{ return m_roadSizeUpX; }

getSector()

char MMT_Road::getSector ( ) const

inline

Definition at line 36 of file MMT_Road.h.

{ return m_sector; }

getUVthreshold()

int MMT_Road::getUVthreshold ( ) const

inline

Definition at line 38 of file MMT_Road.h.

{ return m_uvthr; }

getXthreshold()

int MMT_Road::getXthreshold ( ) const

inline

Definition at line 37 of file MMT_Road.h.

{ return m_xthr; }

horizontalCheck()

bool MMT_Road::horizontalCheck

(

)

const

Definition at line 150 of file MMT_Road.cxx.

                                     {
  int nx1 = 0, nx2 = 0;
  for (const auto &hit : m_road_hits) {
    if (hit->getPlane() >-1 && hit->getPlane() < 2) nx1++;
    if (hit->getPlane() > 5 && hit->getPlane() < 8) nx2++;
  }
  return (nx1 > 0 && nx2 > 0 && (nx1+nx2) >= m_xthr);
}

incrementAge()

void MMT_Road::incrementAge

(

const int

bcwind

)

Definition at line 159 of file MMT_Road.cxx.

                                            {
  std::vector<unsigned int> old_ihits;
  for (unsigned int j = 0; j < m_road_hits.size(); j++) {
    m_road_hits[j]->setAge(m_road_hits[j]->getAge() +1);
    if (m_road_hits[j]->getAge() > (bcwind-1)) old_ihits.push_back(j);
  }
  for (int j = old_ihits.size()-1; j > -1; j--) m_road_hits.erase(m_road_hits.begin()+j);
}

iRoadu()

int MMT_Road::iRoadu ( ) const

inline

Definition at line 40 of file MMT_Road.h.

{ return m_iroadu; }

iRoadv()

int MMT_Road::iRoadv ( ) const

inline

Definition at line 41 of file MMT_Road.h.

{ return m_iroadv; }

iRoadx()

int MMT_Road::iRoadx ( ) const

inline

Definition at line 39 of file MMT_Road.h.

{ return m_iroadx; }

matureCheck()

bool MMT_Road::matureCheck

(

const int

bcwind

)

const

Definition at line 168 of file MMT_Road.cxx.

                                                 {
  for (const auto &hit : m_road_hits) {
    if (hit->getAge() == (bcwind - 1)) return true;
  }
  return false;
}

mxl()

double MMT_Road::mxl

(

)

const

Definition at line 175 of file MMT_Road.cxx.

                           {
  std::vector<double> ys, zs;
  for (const auto &hit : m_road_hits) {
    int bo = hit->getPlane();
    if (bo < 2 || bo > 5) {
      ys.push_back(hit->getR());
      zs.push_back(hit->getZ());
    }
  }
  double mxl = 0;
  double avg_z = std::accumulate(zs.begin(), zs.end(), 0.0)/(double)zs.size();
  double sum_sq_z = std::inner_product(zs.begin(), zs.end(), zs.begin(), 0.0);
  for (unsigned int i = 0; i < ys.size(); i++) mxl += ys[i]*( (zs[i]-avg_z) / (sum_sq_z - zs.size()*std::pow(avg_z,2)) );
 
  return mxl;
}

reset()

void MMT_Road::reset

(

)

Definition at line 192 of file MMT_Road.cxx.

                     {
  if (!m_road_hits.empty()) m_road_hits.clear();
}

stereoCheck()

bool MMT_Road::stereoCheck

(

)

const

Definition at line 196 of file MMT_Road.cxx.

                                 {
 
  if (this->getUVthreshold() == 0) return true;
 
  int nu = 0, nv = 0;
  for (const auto &hit : m_road_hits) {
    if (hit->getPlane() == 2 || hit->getPlane() == 4) nu++;
    if (hit->getPlane() == 3 || hit->getPlane() == 5) nv++;
  }
 
  return (nu > 0 && nv > 0 && (nu+nv) >= m_uvthr);
}

Member Data Documentation

m_iroadu

int MMT_Road::m_iroadu

private

Definition at line 49 of file MMT_Road.h.

m_iroadv

int MMT_Road::m_iroadv

private

Definition at line 50 of file MMT_Road.h.

m_iroadx

int MMT_Road::m_iroadx

private

Definition at line 48 of file MMT_Road.h.

m_road_hits

std::vector<std::unique_ptr<MMT_Hit> > MMT_Road::m_road_hits

private

Definition at line 54 of file MMT_Road.h.

m_roadSize

int MMT_Road::m_roadSize

private

Definition at line 53 of file MMT_Road.h.

m_roadSizeDownUV

int MMT_Road::m_roadSizeDownUV

private

Definition at line 53 of file MMT_Road.h.

m_roadSizeDownX

int MMT_Road::m_roadSizeDownX

private

Definition at line 53 of file MMT_Road.h.

m_roadSizeUpUV

int MMT_Road::m_roadSizeUpUV

private

Definition at line 53 of file MMT_Road.h.

m_roadSizeUpX

int MMT_Road::m_roadSizeUpX

private

Definition at line 53 of file MMT_Road.h.

m_sector

char MMT_Road::m_sector

private

Definition at line 51 of file MMT_Road.h.

m_uvthr

int MMT_Road::m_uvthr

private

Definition at line 52 of file MMT_Road.h.

m_xthr

int MMT_Road::m_xthr

private

Definition at line 52 of file MMT_Road.h.

---

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

• MMT_Road.h
• MMT_Road.cxx