SiStripDetSummary Class Reference

#include <SiStripDetSummary.h>

Classes
struct	Values

Public Member Functions
void	add (DetId detid)
	Used to compute the number of entries divided by subdetector, layer and mono/stereo. More...
void	add (DetId detid, float value)
	Used to compute the mean value of the value variable divided by subdetector, layer and mono/stereo. More...
std::map< unsigned int, Values >	getCounts ()
void	print (std::stringstream &ss, const bool mean=true) const
	SiStripDetSummary (const TrackerTopology *tTopo)

Protected Attributes
bool	computeMean_
std::map< unsigned int, Values >	valueMap_

Private Attributes
const TrackerTopology *	trackerTopo_

Detailed Description

Author

M. De Mattia

Date

26/3/2009 Class to compute and print summary information.

If values are passed together with DetIds (method add( detId, value)), it computes the mean value and rms of a given quantity and is able to print a summary divided by layer/disk for each subdetector.
If instead only DetIds are passed (method add( detId )), it prints the count divided by layer/disk for each subdetector.

Note: consider the possibility to move this class inside SiStripBaseObject as a protected member class.

Definition at line 28 of file SiStripDetSummary.h.

Constructor & Destructor Documentation

SiStripDetSummary()

SiStripDetSummary::SiStripDetSummary ( const TrackerTopology *  tTopo )

inlineexplicit

Definition at line 30 of file SiStripDetSummary.h.

                                                           : computeMean_(true), trackerTopo_(tTopo) {
    // Initialize valueMap_ with zeros
    // WARNING: this initialization is strongly connected with how the map is filled in the add method
    // TIB: layers = 4, stereo = the first 2
    // TOB: layers = 6, stereo = the first 2
    // TEC: wheels = 9, stereo = 9
    // TID: wheels = 3, stereo = 3
    unsigned int layers[] = {4, 6, 9, 3};
    unsigned int stereo[] = {2, 2, 9, 3};
    Values initValues;
    for (unsigned int subDet = 0; subDet < 4; ++subDet) {
      // Layers start from 1
      for (unsigned int layer = 1; layer <= layers[subDet]; ++layer) {
        valueMap_[1000 * (subDet + 1) + layer * 10] = initValues;
        if (layer <= stereo[subDet])
          valueMap_[1000 * (subDet + 1) + layer * 10 + 1] = initValues;
      }
    }
  }

References hgcalTopologyTester_cfi::layers, and valueMap_.

Member Function Documentation

add() [1/2]

void SiStripDetSummary::add ( DetId  detid )

inline

Used to compute the number of entries divided by subdetector, layer and mono/stereo.

Definition at line 53 of file SiStripDetSummary.h.

                               {
    computeMean_ = false;
    add(detid, 0);
  }

References add(), and computeMean_.

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

add() [2/2]

void SiStripDetSummary::add	(	DetId	detid,
		float	value
	)

Used to compute the mean value of the value variable divided by subdetector, layer and mono/stereo.

Definition at line 6 of file SiStripDetSummary.cc.

                                                    {
  int layer = 0;
  int stereo = 0;
  int detNum = 0;
 
  // Using the operator[] if the element does not exist it is created with the default value. That is 0 for integral types.
  switch (detid.subdetId()) {
    case StripSubdetector::TIB:
      layer = trackerTopo_->tibLayer(detid);
      stereo = trackerTopo_->tibStereo(detid);
      detNum = 1000;
      break;
    case StripSubdetector::TOB:
      layer = trackerTopo_->tobLayer(detid);
      stereo = trackerTopo_->tobStereo(detid);
      detNum = 2000;
      break;
    case StripSubdetector::TEC:
      // is this module in TEC+ or TEC-?
      layer = trackerTopo_->tecWheel(detid);
      stereo = trackerTopo_->tecStereo(detid);
      detNum = 3000;
      break;
    case StripSubdetector::TID:
      // is this module in TID+ or TID-?
      layer = trackerTopo_->tidWheel(detid);
      stereo = trackerTopo_->tidStereo(detid);
      detNum = 4000;
      break;
  }
  detNum += layer * 10 + stereo;
  // string name( detector + boost::lexical_cast<string>(layer) + boost::lexical_cast<string>(stereo) );
  valueMap_[detNum].mean += value;
  valueMap_[detNum].rms += value * value;
  valueMap_[detNum].count += 1;
}

References DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecStereo(), TrackerTopology::tecWheel(), StripSubdetector::TIB, TrackerTopology::tibLayer(), TrackerTopology::tibStereo(), StripSubdetector::TID, TrackerTopology::tidStereo(), TrackerTopology::tidWheel(), StripSubdetector::TOB, TrackerTopology::tobLayer(), TrackerTopology::tobStereo(), trackerTopo_, relativeConstraints::value, and valueMap_.

Referenced by add(), SiStripPI::fillNoiseDetSummary(), counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

getCounts()

std::map<unsigned int, Values> SiStripDetSummary::getCounts ( )

inline

Definition at line 70 of file SiStripDetSummary.h.

{ return valueMap_; }

References valueMap_.

print()

void SiStripDetSummary::print	(	std::stringstream &	ss,
		const bool	mean = true
	)		const

Method used to write the output. By default mean == true and it writes the mean value. If mean == false it will write the count.

Definition at line 43 of file SiStripDetSummary.cc.

                                                                        {
  // Compute the mean for each detector and for each layer.
  std::map<unsigned int, Values>::const_iterator valueMapIt = valueMap_.begin();
 
  ss << "subDet" << std::setw(15) << "layer" << std::setw(16) << "mono/stereo" << std::setw(20);
  if (mean)
    ss << "mean +- rms" << std::endl;
  else
    ss << "count" << std::endl;
 
  std::string detector;
  std::string oldDetector;
 
  for (; valueMapIt != valueMap_.end(); ++valueMapIt) {
    int count = valueMapIt->second.count;
    double mean = 0.;
    double rms = 0.;
    if (computeMean_ && count != 0) {
      mean = (valueMapIt->second.mean) / count;
      rms = (valueMapIt->second.rms) / count - mean * mean;
      if (rms <= 0)
        rms = 0;
      else
        rms = sqrt(rms);
    }
    // Detector type
    switch ((valueMapIt->first) / 1000) {
      case 1:
        detector = "TIB ";
        break;
      case 2:
        detector = "TOB ";
        break;
      case 3:
        detector = "TEC ";
        break;
      case 4:
        detector = "TID ";
        break;
    }
    if (detector != oldDetector) {
      ss << std::endl << detector;
      oldDetector = detector;
    } else
      ss << "    ";
    // Layer number
    int layer = (valueMapIt->first) / 10 - (valueMapIt->first) / 1000 * 100;
    int stereo = valueMapIt->first - layer * 10 - (valueMapIt->first) / 1000 * 1000;
 
    ss << std::setw(15) << layer << std::setw(13) << stereo << std::setw(18);
    if (computeMean_)
      ss << mean << " +- " << rms << std::endl;
    else
      ss << count << std::endl;
  }
}

References computeMean_, KineDebug3::count(), hgcalTestNeighbor_cfi::detector, SiStripPI::mean, SiStripPI::rms, mathSSE::sqrt(), contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, and valueMap_.

Member Data Documentation

computeMean_

bool SiStripDetSummary::computeMean_

protected

Definition at line 75 of file SiStripDetSummary.h.

Referenced by add(), and print().

trackerTopo_

const TrackerTopology* SiStripDetSummary::trackerTopo_

private

Definition at line 78 of file SiStripDetSummary.h.

Referenced by add().

valueMap_

std::map<unsigned int, Values> SiStripDetSummary::valueMap_

protected

Definition at line 74 of file SiStripDetSummary.h.

Referenced by add(), getCounts(), print(), and SiStripDetSummary().