#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.

23 { 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.

30 { return m_road_hits; }

◆ getRoadSize()

int MMT_Road::getRoadSize ( ) const

inline

Definition at line 31 of file MMT_Road.h.

31 { return m_roadSize; }

◆ getRoadSizeDownUV()

int MMT_Road::getRoadSizeDownUV ( ) const

inline

Definition at line 35 of file MMT_Road.h.

35 { return m_roadSizeDownUV; }

◆ getRoadSizeDownX()

int MMT_Road::getRoadSizeDownX ( ) const

inline

Definition at line 33 of file MMT_Road.h.

33 { return m_roadSizeDownX; }

◆ getRoadSizeUpUV()

int MMT_Road::getRoadSizeUpUV ( ) const

inline

Definition at line 34 of file MMT_Road.h.

34 { return m_roadSizeUpUV; }

◆ getRoadSizeUpX()

int MMT_Road::getRoadSizeUpX ( ) const

inline

Definition at line 32 of file MMT_Road.h.

32 { return m_roadSizeUpX; }

◆ getSector()

char MMT_Road::getSector ( ) const

inline

Definition at line 36 of file MMT_Road.h.

36 { return m_sector; }

◆ getUVthreshold()

int MMT_Road::getUVthreshold ( ) const

inline

Definition at line 38 of file MMT_Road.h.

38 { return m_uvthr; }

◆ getXthreshold()

int MMT_Road::getXthreshold ( ) const

inline

Definition at line 37 of file MMT_Road.h.

37 { 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.

40 { return m_iroadu; }

◆ iRoadv()

int MMT_Road::iRoadv ( ) const

inline

Definition at line 41 of file MMT_Road.h.

41 { return m_iroadv; }

◆ iRoadx()

int MMT_Road::iRoadx ( ) const

inline

Definition at line 39 of file MMT_Road.h.

39 { 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.