SiStripDetSummary Class Reference

#include <SiStripDetSummary.h>

Classes
struct	Values

Public Member Functions
void	add (DetId detid, float value)
	Used to compute the mean value of the value variable divided by subdetector, layer and mono/stereo. More...
void	add (DetId detid)
	Used to compute the number of entries 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 ( const TrackerTopology *  tTopo )

inlineexplicit

Definition at line 31 of file SiStripDetSummary.h.

References add(), LayerTriplets::layers(), and valueMap_.

                                                           : 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;
      }
    }
  }

Member Function Documentation

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 5 of file SiStripDetSummary.cc.

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(), counter.Counter::register(), SequenceTypes.Task::remove(), SequenceTypes.Task::replace(), and SiStripDetSummary().

{
  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;
}

void SiStripDetSummary::add ( DetId  detid )

inline

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

Definition at line 54 of file SiStripDetSummary.h.

References add(), computeMean_, SiStripPI::mean, and print().

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

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

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

inline

Definition at line 73 of file SiStripDetSummary.h.

References valueMap_.

  {
    return valueMap_;
  }

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.

References computeMean_, KineDebug3::count(), gamEcalExtractorBlocks_cff::detector, SiStripPI::rms, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, and valueMap_.

Referenced by add().

{
  // 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;
  }
}

Member Data Documentation

bool SiStripDetSummary::computeMean_

protected

Definition at line 80 of file SiStripDetSummary.h.

Referenced by add(), and print().

const TrackerTopology* SiStripDetSummary::trackerTopo_

private

Definition at line 82 of file SiStripDetSummary.h.

Referenced by add().

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

protected

Definition at line 79 of file SiStripDetSummary.h.

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