File indexing completed on 2024-04-06 11:58:26
0001
0002
0003
0004
0005 #include "DTT0AbsoluteReferenceCorrection.h"
0006 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0007 #include "FWCore/Framework/interface/EventSetup.h"
0008 #include "FWCore/Framework/interface/ESHandle.h"
0009 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0010
0011 #include "Geometry/Records/interface/MuonGeometryRecord.h"
0012 #include "Geometry/DTGeometry/interface/DTGeometry.h"
0013 #include "Geometry/DTGeometry/interface/DTSuperLayer.h"
0014
0015 #include "DataFormats/MuonDetId/interface/DTWireId.h"
0016 #include "CondFormats/DTObjects/interface/DTT0.h"
0017 #include "CondFormats/DataRecord/interface/DTT0Rcd.h"
0018
0019 #include <string>
0020 #include <sstream>
0021
0022 using namespace std;
0023 using namespace edm;
0024
0025 namespace dtCalibration {
0026
0027 DTT0AbsoluteReferenceCorrection::DTT0AbsoluteReferenceCorrection(const ParameterSet& pset, edm::ConsumesCollector cc)
0028 : calibChamber_(pset.getParameter<string>("calibChamber")),
0029 reference_(pset.getParameter<double>("reference")),
0030 t0Token_(cc.esConsumes<edm::Transition::BeginRun>()) {
0031
0032 if (!calibChamber_.empty() && calibChamber_ != "None" && calibChamber_ != "All") {
0033 stringstream linestr;
0034 int selWheel, selStation, selSector;
0035 linestr << calibChamber_;
0036 linestr >> selWheel >> selStation >> selSector;
0037 chosenChamberId_ = DTChamberId(selWheel, selStation, selSector);
0038 LogVerbatim("Calibration") << "[DTT0AbsoluteReferenceCorrection] Chosen chamber: " << chosenChamberId_ << endl;
0039 }
0040
0041 }
0042
0043 DTT0AbsoluteReferenceCorrection::~DTT0AbsoluteReferenceCorrection() {}
0044
0045 void DTT0AbsoluteReferenceCorrection::setES(const EventSetup& setup) {
0046
0047 ESHandle<DTT0> t0H;
0048 t0H = setup.getHandle(t0Token_);
0049 t0Map_ = &*t0H;
0050 LogVerbatim("Calibration") << "[DTT0AbsoluteReferenceCorrection] T0 version: " << t0H->version();
0051 }
0052
0053 DTT0Data DTT0AbsoluteReferenceCorrection::correction(const DTWireId& wireId) {
0054
0055
0056 DTChamberId chamberId = wireId.layerId().superlayerId().chamberId();
0057
0058 if (calibChamber_.empty() || calibChamber_ == "None")
0059 return defaultT0(wireId);
0060 if (calibChamber_ != "All" && chamberId != chosenChamberId_)
0061 return defaultT0(wireId);
0062
0063
0064 float t0Mean, t0RMS;
0065 int status = t0Map_->get(wireId, t0Mean, t0RMS, DTTimeUnits::counts);
0066 if (status != 0)
0067 throw cms::Exception("[DTT0AbsoluteReferenceCorrection]")
0068 << "Could not find t0 entry in DB for" << wireId << endl;
0069
0070 float t0MeanNew = t0Mean - reference_;
0071 float t0RMSNew = t0RMS;
0072 return DTT0Data(t0MeanNew, t0RMSNew);
0073 }
0074
0075 DTT0Data DTT0AbsoluteReferenceCorrection::defaultT0(const DTWireId& wireId) {
0076
0077 float t0Mean, t0RMS;
0078 int status = t0Map_->get(wireId, t0Mean, t0RMS, DTTimeUnits::counts);
0079 if (!status) {
0080 return DTT0Data(t0Mean, t0RMS);
0081 } else {
0082
0083 throw cms::Exception("[DTT0AbsoluteReferenceCorrection]")
0084 << "Could not find t0 entry in DB for" << wireId << endl;
0085 }
0086 }
0087
0088 }