DTVDriftMeanTimer.cc

Go to the documentation of this file.

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

#include "DTVDriftMeanTimer.h"

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

#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "CondFormats/DTObjects/interface/DTMtime.h"

#include "CalibMuon/DTCalibration/interface/DTMeanTimerFitter.h"
#include "CalibMuon/DTCalibration/interface/DTCalibDBUtils.h"

#include <string>
#include <vector>

#include "TFile.h"
#include "TString.h"

using namespace std;
using namespace edm;

namespace dtCalibration {

DTVDriftMeanTimer::DTVDriftMeanTimer(const ParameterSet& pset) {
  string rootFileName = pset.getParameter<string>("rootFileName");
  rootFile_ = new TFile(rootFileName.c_str(), "READ");
  fitter_ = new DTMeanTimerFitter(rootFile_);
  bool debug = pset.getUntrackedParameter<bool>("debug", false);
  if(debug) fitter_->setVerbosity(1);
}

DTVDriftMeanTimer::~DTVDriftMeanTimer() {
  rootFile_->Close();
  delete fitter_;
}

void DTVDriftMeanTimer::setES(const edm::EventSetup& setup) {}

DTVDriftData DTVDriftMeanTimer::compute(DTSuperLayerId const& slId) {

  // Evaluate v_drift and sigma from the TMax histograms
  DTWireId wireId(slId, 0, 0);
  TString N = ( ( ( ( (TString)"TMax" + (long)wireId.wheel() )
                                      + (long)wireId.station() )
                              + (long)wireId.sector() ) + (long)wireId.superLayer() );
  vector<float> vDriftAndReso = fitter_->evaluateVDriftAndReso(N);

  // Don't write the constants for the SL if the vdrift was not computed
  if(vDriftAndReso.front() == -1)
     throw cms::Exception("DTCalibration") << "Could not compute valid vDrift value for SL " << slId << endl;

  return DTVDriftData(vDriftAndReso[0],vDriftAndReso[1]);
}

} // namespace