#include <DTOccupancyCluster.h>

Public Member Functions
bool	addPoint (const DTOccupancyPoint &anotherPoint)

double	averageMean () const
	average cell occupancy of the layers in the cluster More...

double	averageRMS () const
	average RMS of the cell occpuancy distributions of the layers in the cluster More...

double	distance (const DTOccupancyPoint &point) const

	DTOccupancyCluster (const DTOccupancyPoint &firstPoint, const DTOccupancyPoint &secondPoint)
	Constructor. More...

	DTOccupancyCluster (const DTOccupancyPoint &singlePoint)
	Constructor. More...

TH2F *	getHisto (std::string histoName, int nBinsX, double minX, double maxX, int nBinsY, double minY, double maxY, int fillColor) const
	get a TH2F displaying the cluster More...

std::set< DTLayerId >	getLayerIDs () const

bool	isValid () const
	Check if the cluster candidate satisfies the quality requirements. More...

double	maxMean () const
	max average cell occupancy of the layers in the cluster More...

double	maxRMS () const
	max RMS of the cell occpuancy distributions of the layers in the cluster More...

int	nPoints () const

virtual	~DTOccupancyCluster ()
	Destructor. More...

Private Member Functions
void	computeRadius ()

bool	qualityCriterion (const DTOccupancyPoint &anotherPoint)

bool	qualityCriterion (const DTOccupancyPoint &firstPoint, const DTOccupancyPoint &secondPoint)

Private Attributes
double	radius

double	theMaxMean

double	theMaxRMS

double	theMeanSum

std::vector< DTOccupancyPoint >	thePoints

double	theRMSSum

bool	theValidity

Detailed Description

Cluster of DTOccupancyPoint used bt DTOccupancyTest to spot problematic layers. Layers are clusterized in the plane average cell occupancy - RMS of the cell occupancies.

Author: G. Cerminara - INFN Torino

Definition at line 20 of file DTOccupancyCluster.h.

Constructor & Destructor Documentation

◆ DTOccupancyCluster() [1/2]

DTOccupancyCluster::DTOccupancyCluster	(	const DTOccupancyPoint &	firstPoint,
		const DTOccupancyPoint &	secondPoint
	)

Constructor.

Definition at line 19 of file DTOccupancyCluster.cc.

     : radius(0), theMaxMean(-1.), theMaxRMS(-1.), theMeanSum(0.), theRMSSum(0.) {
   if (!qualityCriterion(firstPoint, secondPoint)) {
     theValidity = false;
   } else {
     // compute the cluster quantities
     thePoints.push_back(firstPoint);
     thePoints.push_back(secondPoint);
     if (firstPoint.mean() > secondPoint.mean())
       theMaxMean = firstPoint.mean();
     else
       theMaxMean = secondPoint.mean();
  
     if (firstPoint.rms() > secondPoint.rms())
       theMaxRMS = firstPoint.rms();
     else
       theMaxRMS = secondPoint.rms();
     theMeanSum += firstPoint.mean();
     theRMSSum += firstPoint.rms();
  
     theMeanSum += secondPoint.mean();
     theRMSSum += secondPoint.rms();
  
     computeRadius();
   }
 }

References computeRadius(), DTOccupancyPoint::mean(), qualityCriterion(), DTOccupancyPoint::rms(), theMaxMean, theMaxRMS, theMeanSum, thePoints, theRMSSum, and theValidity.

◆ DTOccupancyCluster() [2/2]

DTOccupancyCluster::DTOccupancyCluster ( const DTOccupancyPoint & singlePoint )

Constructor.

Definition at line 46 of file DTOccupancyCluster.cc.

     : radius(0),
       theMaxMean(singlePoint.mean()),
       theMaxRMS(singlePoint.rms()),
       theMeanSum(singlePoint.mean()),
       theRMSSum(singlePoint.rms()) {
   theValidity = true;
  
   // compute the cluster quantities
   thePoints.push_back(singlePoint);
 }

References thePoints, and theValidity.

◆ ~DTOccupancyCluster()

DTOccupancyCluster::~DTOccupancyCluster ( )

virtual

Destructor.

Definition at line 58 of file DTOccupancyCluster.cc.

58 {}

Member Function Documentation

◆ addPoint()

bool DTOccupancyCluster::addPoint ( const DTOccupancyPoint & anotherPoint )

Add a point to the cluster: returns false if the point does not satisfy the quality requirement

Definition at line 65 of file DTOccupancyCluster.cc.

                                                                       {
   LogTrace("DTDQM|DTMonitorClient|DTOccupancyTest|DTOccupancyCluster")
       << "   Add a point to the cluster: mean: " << anotherPoint.mean() << " rms: " << anotherPoint.rms() << endl;
   if (qualityCriterion(anotherPoint)) {
     LogTrace("DTDQM|DTMonitorClient|DTOccupancyTest|DTOccupancyCluster") << "   point is valid" << endl;
     thePoints.push_back(anotherPoint);
     // Compute the new cluster size
     computeRadius();
     // compute the max mean and RMS
     if (anotherPoint.mean() > theMaxMean) {
       theMaxMean = anotherPoint.mean();
     }
     if (anotherPoint.rms() > theMaxRMS) {
       theMaxRMS = anotherPoint.rms();
     }
     theMeanSum += anotherPoint.mean();
     theRMSSum += anotherPoint.rms();
     return true;
   }
   return false;
 }

References computeRadius(), LogTrace, DTOccupancyPoint::mean(), qualityCriterion(), DTOccupancyPoint::rms(), theMaxMean, theMaxRMS, theMeanSum, thePoints, and theRMSSum.

Referenced by DTOccupancyClusterBuilder::buildNewCluster().

◆ averageMean()

double DTOccupancyCluster::averageMean ( ) const

average cell occupancy of the layers in the cluster

Definition at line 101 of file DTOccupancyCluster.cc.

101 { return theMeanSum / (double)thePoints.size(); }

References theMeanSum, and thePoints.

Referenced by DTOccupancyClusterBuilder::buildNewCluster(), clusterIsLessThan(), and qualityCriterion().

◆ averageRMS()

double DTOccupancyCluster::averageRMS ( ) const

average RMS of the cell occpuancy distributions of the layers in the cluster

Definition at line 103 of file DTOccupancyCluster.cc.

103 { return theRMSSum / (double)thePoints.size(); }

References thePoints, and theRMSSum.

Referenced by DTOccupancyClusterBuilder::buildNewCluster(), clusterIsLessThan(), and qualityCriterion().

◆ computeRadius()

void DTOccupancyCluster::computeRadius ( )

private

Definition at line 153 of file DTOccupancyCluster.cc.

                                        {
   double radius_squared = 0;
   for (vector<DTOccupancyPoint>::const_iterator pt_i = thePoints.begin(); pt_i != thePoints.end(); ++pt_i) {
     for (vector<DTOccupancyPoint>::const_iterator pt_j = thePoints.begin(); pt_j != thePoints.end(); ++pt_j) {
       radius_squared += TMath::Power(pt_i->distance(*pt_j), 2);
     }
   }
   radius_squared = radius_squared / (2 * TMath::Power(thePoints.size() + 1, 2));
   radius = sqrt(radius_squared);
 }

References radius, mathSSE::sqrt(), and thePoints.

Referenced by addPoint(), and DTOccupancyCluster().

◆ distance()

double DTOccupancyCluster::distance ( const DTOccupancyPoint & point ) const

Compute the distance of a single point from the cluster (minimum distance with respect to the cluster points)

Definition at line 89 of file DTOccupancyCluster.cc.

                                                                        {
   double dist = 99999999;
   // compute the minimum distance from a point
   for (vector<DTOccupancyPoint>::const_iterator pt = thePoints.begin(); pt != thePoints.end(); ++pt) {
     double distance = point.distance(*pt);
     if (distance < dist) {
       dist = distance;
     }
   }
   return dist;
 }

References point, DiDispStaMuonMonitor_cfi::pt, and thePoints.

Referenced by DTOccupancyClusterBuilder::computeDistancesToCluster().

◆ getHisto()

TH2F * DTOccupancyCluster::getHisto	(	std::string	histoName,
		int	nBinsX,
		double	minX,
		double	maxX,
		int	nBinsY,
		double	minY,
		double	maxY,
		int	fillColor
	)		const

get a TH2F displaying the cluster

Definition at line 109 of file DTOccupancyCluster.cc.

                                                         {
   TH2F* histo = new TH2F(histoName.c_str(), histoName.c_str(), nBinsX, minX, maxX, nBinsY, minY, maxY);
   histo->SetFillColor(fillColor);
   for (vector<DTOccupancyPoint>::const_iterator pt = thePoints.begin(); pt != thePoints.end(); ++pt) {
     histo->Fill((*pt).mean(), (*pt).rms());
   }
   return histo;
 }

References RecoTauValidation_cfi::fillColor, timingPdfMaker::histo, HltBtagPostValidation_cff::histoName, HGVHistoProducerAlgoBlock_cfi::maxX, HGVHistoProducerAlgoBlock_cfi::maxY, HGVHistoProducerAlgoBlock_cfi::minX, HGVHistoProducerAlgoBlock_cfi::minY, DiDispStaMuonMonitor_cfi::pt, and thePoints.

◆ getLayerIDs()

set< DTLayerId > DTOccupancyCluster::getLayerIDs ( ) const

Definition at line 166 of file DTOccupancyCluster.cc.

                                                      {
   set<DTLayerId> ret;
   for (vector<DTOccupancyPoint>::const_iterator point = thePoints.begin(); point != thePoints.end(); ++point) {
     ret.insert((*point).layerId());
   }
   return ret;
 }

References point, runTheMatrix::ret, and thePoints.

◆ isValid()

bool DTOccupancyCluster::isValid ( void ) const

Check if the cluster candidate satisfies the quality requirements.

Definition at line 61 of file DTOccupancyCluster.cc.

61 { return theValidity; }

References theValidity.

Referenced by ntupleDataFormat._Object::_checkIsValid(), DTOccupancyClusterBuilder::buildNewCluster(), and core.AutoHandle.AutoHandle::ReallyLoad().

◆ maxMean()

double DTOccupancyCluster::maxMean ( ) const

max average cell occupancy of the layers in the cluster

Definition at line 105 of file DTOccupancyCluster.cc.

105 { return theMaxMean; }

References theMaxMean.

Referenced by DTOccupancyClusterBuilder::buildClusters(), and DTOccupancyClusterBuilder::buildNewCluster().

◆ maxRMS()

double DTOccupancyCluster::maxRMS ( ) const

max RMS of the cell occpuancy distributions of the layers in the cluster

Definition at line 107 of file DTOccupancyCluster.cc.

107 { return theMaxRMS; }

References theMaxRMS.

Referenced by DTOccupancyClusterBuilder::buildClusters(), and DTOccupancyClusterBuilder::buildNewCluster().

◆ nPoints()

int DTOccupancyCluster::nPoints ( ) const

of layers belonging to the cluster

Definition at line 164 of file DTOccupancyCluster.cc.

164 { return thePoints.size(); }

References thePoints.

Referenced by DTOccupancyClusterBuilder::buildNewCluster(), and clusterIsLessThan().

◆ qualityCriterion() [1/2]

bool DTOccupancyCluster::qualityCriterion ( const DTOccupancyPoint & anotherPoint )

private

Definition at line 137 of file DTOccupancyCluster.cc.

                                                                               {
   double minrms = 0;
   if (anotherPoint.rms() < averageRMS())
     minrms = anotherPoint.rms();
   else
     minrms = averageRMS();
  
   if (fabs(averageMean() - anotherPoint.mean()) < averageRMS() &&
       fabs(averageRMS() - anotherPoint.rms()) < 2 * minrms / 3.) {
     theValidity = true;
     return true;
   }
   theValidity = false;
   return false;
 }

References averageMean(), averageRMS(), DTOccupancyPoint::mean(), DTOccupancyPoint::rms(), and theValidity.

◆ qualityCriterion() [2/2]

bool DTOccupancyCluster::qualityCriterion	(	const DTOccupancyPoint &	firstPoint,
		const DTOccupancyPoint &	secondPoint
	)

private

Definition at line 125 of file DTOccupancyCluster.cc.

                                                                                                                  {
   if (firstPoint.deltaMean(secondPoint) < computeAverageRMS(firstPoint, secondPoint) &&
       firstPoint.deltaRMS(secondPoint) < computeMinRMS(firstPoint, secondPoint)) {
     theValidity = true;
  
     return true;
   }
  
   theValidity = false;
   return false;
 }

References computeAverageRMS(), computeMinRMS(), DTOccupancyPoint::deltaMean(), DTOccupancyPoint::deltaRMS(), and theValidity.

Referenced by addPoint(), and DTOccupancyCluster().

Member Data Documentation

◆ radius

double DTOccupancyCluster::radius

private

Definition at line 82 of file DTOccupancyCluster.h.

Referenced by computeRadius().

◆ theMaxMean

double DTOccupancyCluster::theMaxMean

private

Definition at line 85 of file DTOccupancyCluster.h.

Referenced by addPoint(), DTOccupancyCluster(), and maxMean().

◆ theMaxRMS

double DTOccupancyCluster::theMaxRMS

private

Definition at line 86 of file DTOccupancyCluster.h.

Referenced by addPoint(), DTOccupancyCluster(), and maxRMS().

◆ theMeanSum

double DTOccupancyCluster::theMeanSum

private

Definition at line 87 of file DTOccupancyCluster.h.

Referenced by addPoint(), averageMean(), and DTOccupancyCluster().

◆ thePoints

std::vector<DTOccupancyPoint> DTOccupancyCluster::thePoints

private

Definition at line 83 of file DTOccupancyCluster.h.

Referenced by addPoint(), averageMean(), averageRMS(), computeRadius(), distance(), DTOccupancyCluster(), getHisto(), getLayerIDs(), and nPoints().

◆ theRMSSum

double DTOccupancyCluster::theRMSSum

private

Definition at line 88 of file DTOccupancyCluster.h.

Referenced by addPoint(), averageRMS(), and DTOccupancyCluster().

◆ theValidity

bool DTOccupancyCluster::theValidity

private

Definition at line 81 of file DTOccupancyCluster.h.

Referenced by DTOccupancyCluster(), isValid(), and qualityCriterion().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ DTOccupancyCluster() [1/2]

◆ DTOccupancyCluster() [2/2]

◆ ~DTOccupancyCluster()

Member Function Documentation

◆ addPoint()

◆ averageMean()

◆ averageRMS()

◆ computeRadius()

◆ distance()

◆ getHisto()

◆ getLayerIDs()

◆ isValid()

◆ maxMean()

◆ maxRMS()

◆ nPoints()

of layers belonging to the cluster

◆ qualityCriterion() [1/2]

◆ qualityCriterion() [2/2]

Member Data Documentation

◆ radius

◆ theMaxMean

◆ theMaxRMS

◆ theMeanSum

◆ thePoints

◆ theRMSSum

◆ theValidity