DTT0ChamberReferenceCorrection.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author A. Vilela Pereira
 */

#include "DTT0ChamberReferenceCorrection.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "CondFormats/DTObjects/interface/DTT0.h"
#include "CondFormats/DataRecord/interface/DTT0Rcd.h"

#include <string>
#include <sstream>

using namespace std;
using namespace edm;

namespace dtCalibration {

DTT0ChamberReferenceCorrection::DTT0ChamberReferenceCorrection(const ParameterSet& pset):
  calibChamber_( pset.getParameter<string>("calibChamber") ) {

  //DTChamberId chosenChamberId;
  if( calibChamber_ != "" && calibChamber_ != "None" && calibChamber_ != "All" ){
    stringstream linestr;
    int selWheel, selStation, selSector;
    linestr << calibChamber_;
    linestr >> selWheel >> selStation >> selSector;
    chosenChamberId_ = DTChamberId(selWheel, selStation, selSector);
    LogVerbatim("Calibration") << "[DTT0ChamberReferenceCorrection] Chosen chamber: " << chosenChamberId_ << endl;
  }
  //FIXME: Check if chosen chamber is valid.

}

DTT0ChamberReferenceCorrection::~DTT0ChamberReferenceCorrection() {
}

void DTT0ChamberReferenceCorrection::setES(const EventSetup& setup) {
  // Get t0 record from DB
  ESHandle<DTT0> t0H;
  setup.get<DTT0Rcd>().get(t0H);
  t0Map_ = &*t0H;
  LogVerbatim("Calibration") << "[DTT0ChamberReferenceCorrection] T0 version: " << t0H->version();
}

DTT0Data DTT0ChamberReferenceCorrection::correction(const DTWireId& wireId) {
  // Compute for selected chamber (or All) correction using as reference chamber mean
 
  DTChamberId chamberId = wireId.layerId().superlayerId().chamberId();

  if( calibChamber_ == "" || calibChamber_ == "None" )          return defaultT0(wireId);
  if( calibChamber_ != "All" && chamberId != chosenChamberId_ ) return defaultT0(wireId);

  // Access DB
  float t0Mean,t0RMS;
  int status = t0Map_->get(wireId,t0Mean,t0RMS,DTTimeUnits::counts);
  if(status != 0) 
     throw cms::Exception("[DTT0ChamberReferenceCorrection]") << "Could not find t0 entry in DB for" 
                                                              << wireId << endl;
  /*
    Leaving just the structure for future implementation
    ...
    ...
  */
  return DTT0Data(t0Mean,t0RMS);
}

DTT0Data DTT0ChamberReferenceCorrection::defaultT0(const DTWireId& wireId) {
  // Access default DB
  float t0Mean,t0RMS;
  int status = t0Map_->get(wireId,t0Mean,t0RMS,DTTimeUnits::counts);
  if(!status){
     return DTT0Data(t0Mean,t0RMS);
  } else{
     //... 
     throw cms::Exception("[DTT0ChamberReferenceCorrection]") << "Could not find t0 entry in DB for"
                                                      << wireId << endl;
  }
}

} // namespace