DTCalibrationMap Class Reference

#include <DTCalibrationMap.h>

Public Types
typedef std::vector< float >	CalibConsts
typedef std::map< Key, CalibConsts >::const_iterator	const_iterator
typedef DTWireId	Key

Public Member Functions
void	addCell (Key wireId, const CalibConsts &calibConst)
void	cleanTheConsts ()
	DTCalibrationMap (const edm::ParameterSet &pset)
	Constructor. More...
const CalibConsts *	getConsts (DTWireId wireId) const
float	getField (DTWireId wireId, int field) const
Key	getKey (DTWireId wireId) const
const_iterator	keyAndConsts_begin () const
const_iterator	keyAndConsts_end () const
float	kFactor (DTWireId wireId) const
	Return the kfactor for a particular wire. More...
float	meanVDrift (DTWireId wireId) const
	Return the mean drift velocity for a particular wire (cm/ns) More...
float	sigma_meanVDrift (DTWireId wireId) const
	Return the sigma of the mean drift velocity for a particular wire (cm/ns) More...
float	sigma_tTrig (DTWireId wireId) const
	Return the sigma of the t_trig (ns) for a particular wire. More...
float	tTrig (DTWireId wireId) const
	Return the t_trig (ns) for a particular wire. More...
void	writeConsts (const std::string &outputFileName) const
virtual	~DTCalibrationMap ()
	Destructor. More...

Private Types
enum	CalibGranularity { byChamber, bySL, byLayer, byWire }

Private Member Functions
bool	checkGranularity (Key aKey) const
void	readConsts (const std::string &inputFileName)

Private Attributes
std::string	calibConstFileName
std::string	calibConstGranularity
unsigned int	nFields
CalibGranularity	theGranularity
std::map< Key, CalibConsts >	theMap

Detailed Description

Allow saving and retrieving of calibration constants to/from txt file. This was originally provided for backward compatibility with the ORCA MuBarDigiParameters file. Can be used to save an arbitrary number of constants with the needed granularity and to retrieve them back using the wireId. Current field allocation: fields for each key are allocated to: –First block: legacy descriptors– [0] ttrig [1] sigma_ttrig (obsolete) [2] kfactor (obsolete) [3] vdrift
[4] sigma_vdrift (obsolete, was formerly hacked to include reoslution) [5] t0 [6] t0rms [7] noisy or dead flag [8-9] left for future usage –Second block (optional): free fields [10] Encoded information of free fields: (1000*version)+(100*type)+(number of fields); type is: ttrig=0, vdrift=1, uncertainties=3 [11-end] free fields

Author

G. Cerminara - INFN Torino

Definition at line 37 of file DTCalibrationMap.h.

Member Typedef Documentation

CalibConsts

typedef std::vector<float> DTCalibrationMap::CalibConsts

Definition at line 62 of file DTCalibrationMap.h.

const_iterator

typedef std::map<Key, CalibConsts>::const_iterator DTCalibrationMap::const_iterator

Definition at line 64 of file DTCalibrationMap.h.

Key

typedef DTWireId DTCalibrationMap::Key

Definition at line 63 of file DTCalibrationMap.h.

Member Enumeration Documentation

CalibGranularity

enum DTCalibrationMap::CalibGranularity

private

Enumerator
byChamber
bySL
byLayer
byWire

Definition at line 93 of file DTCalibrationMap.h.

{ byChamber, bySL, byLayer, byWire };

Constructor & Destructor Documentation

DTCalibrationMap()

DTCalibrationMap::DTCalibrationMap

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 23 of file DTCalibrationMap.cc.

References gather_cfg::cout, and muonDTDigis_cfi::pset.

                                                           {
  nFields = pset.getUntrackedParameter<int>("nFields", 5);
  calibConstFileName = pset.getUntrackedParameter<string>("calibConstFileName", "dummy.txt");
  calibConstGranularity = pset.getUntrackedParameter<string>("calibConstGranularity", "bySL");

  // Initialize correctly the enum which specify the granularity for the calibration
  if (calibConstGranularity == "byWire") {
    theGranularity = byWire;
  } else if (calibConstGranularity == "byLayer") {
    theGranularity = byLayer;
  } else if (calibConstGranularity == "bySL") {
    theGranularity = bySL;
  } else {
    theGranularity = byChamber;
    if (!(calibConstGranularity == "byChamber")) {
      cout << "[DTCalibrationMap]###Warning: Check parameter calibConstGranularity: " << calibConstGranularity
           << " options not available!" << endl;
    }
  }
  readConsts(calibConstFileName);
}

~DTCalibrationMap()

DTCalibrationMap::~DTCalibrationMap ( )

virtual

Destructor.

Definition at line 45 of file DTCalibrationMap.cc.

{}

Member Function Documentation

addCell()

void DTCalibrationMap::addCell	(	Key	wireId,
		const CalibConsts &	calibConst
	)

Definition at line 160 of file DTCalibrationMap.cc.

References Exception.

Referenced by DTVDriftCalibration::endJob().

                                                                        {
  if (!checkGranularity(theKey))
    throw cms::Exception("addCell") << "DTCalibrationMap:" << endl
                                    << "The added key is not compatible with the selected granularity" << endl;

  theMap[theKey] = calibConst;
}

checkGranularity()

bool DTCalibrationMap::checkGranularity ( Key  aKey ) const

private

Definition at line 180 of file DTCalibrationMap.cc.

References DTLayerId::layer(), runTheMatrix::ret, DTSuperLayerId::superlayer(), and DTWireId::wire().

                                                      {
  bool ret = true;

  // Check that the key is consistent with the given granularity
  if (theGranularity == byChamber) {
    if (aKey.superlayer() || aKey.layer() || aKey.wire()) {
      ret = false;
    }
  } else if (theGranularity == bySL) {
    if (aKey.layer() || aKey.wire()) {
      ret = false;
    }
  } else if (theGranularity == byLayer) {
    if (aKey.wire()) {
      ret = false;
    }
  } else if (theGranularity == byWire) {
    if (aKey.wire() == 0) {
      ret = false;
    }
  }
  return ret;
}

cleanTheConsts()

void DTCalibrationMap::cleanTheConsts ( )

inline

Definition at line 67 of file DTCalibrationMap.h.

References theMap.

{ theMap.clear(); }

getConsts()

const DTCalibrationMap::CalibConsts * DTCalibrationMap::getConsts

(

DTWireId

wireId

)

const

Definition at line 77 of file DTCalibrationMap.cc.

References utilities::cache().

Referenced by DTVDriftCalibration::endJob().

                                                                                    {
  // Create a cache
  static pair<Key, CalibConsts> cache;

  // Get the key
  Key theKey = getKey(wireId);

  // Check if the result is already cached
  if (theKey == cache.first) {
    return &(cache.second);
  }

  // Look for the given key into the map
  map<Key, CalibConsts>::const_iterator res = theMap.find(theKey);
  if (res != theMap.end()) {
    cache = (*res);
    return &((*res).second);
  } else {
    return nullptr;
  }
}

getField()

float DTCalibrationMap::getField	(	DTWireId	wireId,
		int	field
	)		const

Definition at line 101 of file DTCalibrationMap.cc.

References Exception.

                                                                 {
  const CalibConsts* cals = getConsts(wireId);
  if (cals == nullptr) {
    throw cms::Exception("NoCalibConsts") << "DTCalibrationMap:" << endl
                                          << "No parameters for wire: " << wireId << endl
                                          << "Check the " << calibConstFileName << " file!" << endl;
  } else {
    return (*(cals))[field];
  }
}

getKey()

DTCalibrationMap::Key DTCalibrationMap::getKey

(

DTWireId

wireId

)

const

Definition at line 64 of file DTCalibrationMap.cc.

References DTSuperLayerId::chamberId(), DTWireId::layerId(), and DTLayerId::superlayerId().

Referenced by DTVDriftCalibration::endJob(), MatrixToProcess.MatrixToProcess::getConfig(), MatrixToProcess.MatrixToProcess::getProcess(), MatrixToProcess.MatrixToProcess::identityTest(), and MatrixToProcess.MatrixToProcess::load().

                                                                  {
  if (theGranularity == byChamber) {
    return Key(wireId.chamberId(), 0, 0, 0);
  } else if (theGranularity == bySL) {
    return Key(wireId.superlayerId(), 0, 0);
  } else if (theGranularity == byLayer) {
    return Key(wireId.layerId(), 0);
  } else {
    return Key(wireId);
  }
}

keyAndConsts_begin()

const_iterator DTCalibrationMap::keyAndConsts_begin ( ) const

inline

Definition at line 86 of file DTCalibrationMap.h.

References theMap.

{ return theMap.begin(); }

keyAndConsts_end()

const_iterator DTCalibrationMap::keyAndConsts_end ( ) const

inline

Definition at line 88 of file DTCalibrationMap.h.

References theMap.

{ return theMap.end(); }

kFactor()

float DTCalibrationMap::kFactor

(

DTWireId

wireId

)

const

Return the kfactor for a particular wire.

Definition at line 54 of file DTCalibrationMap.cc.

References CaloRecHitAuxSetter::getField().

{ return getField(wireId, 2); }

meanVDrift()

float DTCalibrationMap::meanVDrift

(

DTWireId

wireId

)

const

Return the mean drift velocity for a particular wire (cm/ns)

Definition at line 57 of file DTCalibrationMap.cc.

References CaloRecHitAuxSetter::getField().

{ return getField(wireId, 3); }

readConsts()

void DTCalibrationMap::readConsts ( const std::string &  inputFileName )

private

Definition at line 113 of file DTCalibrationMap.cc.

References filterCSVwithJSON::copy, gather_cfg::cout, geometryDiff::file, InefficientDoubleROC::inputFileName, and mps_splice::line.

                                                             {
  ifstream file(inputFileName.c_str());
  // Check if the file exists
  if (!file) {
    cout << "[DTCalibrationMap]***Warning: File: " << inputFileName << " not found in current directory!!!" << endl;
  }

  string line;

  // The numbers to be read to build the key
  int wheel_id = 0;
  int station_id = 0;
  int sector_id = 0;
  int superlayer_id = 0;
  int layer_id = 0;
  int wire_id = 0;

  // Read all the lines
  while (getline(file, line)) {
    if (line.empty() || line[0] == '#')
      continue;  // Skip comments and empty lines
    stringstream linestr;
    linestr << line;

    pair<Key, CalibConsts> wireCalib;

    linestr >> wheel_id >> station_id >> sector_id >> superlayer_id >> layer_id >> wire_id;

    // Build the key
    wireCalib.first = Key(wheel_id, station_id, sector_id, superlayer_id, layer_id, wire_id);

    if (!checkGranularity(wireCalib.first))
      cout << "[DTCalibrationMap]***Warning: the CalibConstFile is not consistent with the selected granularity!"
           << endl;

    // Read the calibration constants
    copy(istream_iterator<float>(linestr), istream_iterator<float>(), back_inserter(wireCalib.second));

    if (wireCalib.second.size() < nFields) {
      cout << "[DTCalibrationMap]***Warning: the CalibConstFile is not consistent with the number of fields!" << endl;
    }

    theMap.insert(wireCalib);
  }
}

sigma_meanVDrift()

float DTCalibrationMap::sigma_meanVDrift

(

DTWireId

wireId

)

const

Return the sigma of the mean drift velocity for a particular wire (cm/ns)

Definition at line 60 of file DTCalibrationMap.cc.

References CaloRecHitAuxSetter::getField().

{ return getField(wireId, 4); }

sigma_tTrig()

float DTCalibrationMap::sigma_tTrig

(

DTWireId

wireId

)

const

Return the sigma of the t_trig (ns) for a particular wire.

Definition at line 51 of file DTCalibrationMap.cc.

References CaloRecHitAuxSetter::getField().

{ return getField(wireId, 1); }

tTrig()

float DTCalibrationMap::tTrig

(

DTWireId

wireId

)

const

Return the t_trig (ns) for a particular wire.

Definition at line 48 of file DTCalibrationMap.cc.

References CaloRecHitAuxSetter::getField().

{ return getField(wireId, 0); }

writeConsts()

void DTCalibrationMap::writeConsts

(

const std::string &

outputFileName

)

const

Definition at line 169 of file DTCalibrationMap.cc.

References filterCSVwithJSON::copy, MillePedeFileConverter_cfg::out, and reco_skim_cfg_mod::outputFileName.

Referenced by DTVDriftCalibration::endJob().

                                                                     {
  ofstream out(outputFileName.c_str());
  for (map<Key, CalibConsts>::const_iterator iter = theMap.begin(); iter != theMap.end(); ++iter) {
    out << (*iter).first.wheel() << ' ' << (*iter).first.station() << ' ' << (*iter).first.sector() << ' '
        << (*iter).first.superlayer() << ' ' << (*iter).first.layer() << ' ' << (*iter).first.wire() << ' ';
    copy((*iter).second.begin(), (*iter).second.end(), ostream_iterator<float>(out, " "));
    out << endl;
  }
}

Member Data Documentation

calibConstFileName

std::string DTCalibrationMap::calibConstFileName

private

Definition at line 106 of file DTCalibrationMap.h.

calibConstGranularity

std::string DTCalibrationMap::calibConstGranularity

private

Definition at line 109 of file DTCalibrationMap.h.

nFields

unsigned int DTCalibrationMap::nFields

private

Definition at line 103 of file DTCalibrationMap.h.

theGranularity

CalibGranularity DTCalibrationMap::theGranularity

private

Definition at line 94 of file DTCalibrationMap.h.

theMap

std::map<Key, CalibConsts> DTCalibrationMap::theMap

private

Definition at line 112 of file DTCalibrationMap.h.

Referenced by cleanTheConsts(), keyAndConsts_begin(), and keyAndConsts_end().