ThrParameters.cc

Go to the documentation of this file.

#include "RecoMuon/GlobalTrackingTools/interface/ThrParameters.h"
#include "FWCore/Framework/interface/ESHandle.h"
 
using namespace edm;
using namespace std;
 
ThrParameters::ThrParameters(const EventSetup* eSetup) {
  // Read the threshold DB
  ESHandle<DYTThrObject> dytThresholdsH;
 
  // This try catch is just temporary and
  // will be removed as soon as the DYTThrObject
  // record is included in a GT.
  // It is necessary here for testing
  try {
    eSetup->get<DYTThrObjectRcd>().get(dytThresholdsH);
    dytThresholds = dytThresholdsH.product();
    isValidThdDB_ = true;
  } catch (...) {
    isValidThdDB_ = false;
  }
 
  // Ape are always filled even they're null
  ESHandle<AlignmentErrorsExtended> dtAlignmentErrorsExtended;
  eSetup->get<DTAlignmentErrorExtendedRcd>().get(dtAlignmentErrorsExtended);
  for (std::vector<AlignTransformErrorExtended>::const_iterator it = dtAlignmentErrorsExtended->m_alignError.begin();
       it != dtAlignmentErrorsExtended->m_alignError.end();
       it++) {
    CLHEP::HepSymMatrix error = (*it).matrix();
    GlobalError glbErr(error[0][0], error[1][0], error[1][1], error[2][0], error[2][1], error[2][2]);
    DTChamberId DTid((*it).rawId());
    dtApeMap.insert(pair<DTChamberId, GlobalError>(DTid, glbErr));
  }
  ESHandle<AlignmentErrorsExtended> cscAlignmentErrorsExtended;
  eSetup->get<CSCAlignmentErrorExtendedRcd>().get(cscAlignmentErrorsExtended);
  for (std::vector<AlignTransformErrorExtended>::const_iterator it = cscAlignmentErrorsExtended->m_alignError.begin();
       it != cscAlignmentErrorsExtended->m_alignError.end();
       it++) {
    CLHEP::HepSymMatrix error = (*it).matrix();
    GlobalError glbErr(error[0][0], error[1][0], error[1][1], error[2][0], error[2][1], error[2][2]);
    CSCDetId CSCid((*it).rawId());
    cscApeMap.insert(pair<CSCDetId, GlobalError>(CSCid, glbErr));
  }
}
 
ThrParameters::~ThrParameters() {}