dtCalibration::DTVDriftSegment Class Reference

#include <DTVDriftSegment.h>

Inheritance diagram for dtCalibration::DTVDriftSegment:

Public Member Functions
DTVDriftData	compute (const DTSuperLayerId &) override
	DTVDriftSegment (edm::ParameterSet const &)
void	setES (const edm::EventSetup &setup) override
	~DTVDriftSegment () override
Public Member Functions inherited from dtCalibration::DTVDriftBaseAlgo
	DTVDriftBaseAlgo ()
virtual	~DTVDriftBaseAlgo ()

Private Member Functions
TH1F *	getHisto (const DTSuperLayerId &)
std::string	getHistoName (const DTSuperLayerId &)

Private Attributes
DTResidualFitter *	fitter_
const DTMtime *	mTimeMap_
unsigned int	nSigmas_
TFile *	rootFile_

Detailed Description

Definition at line 23 of file DTVDriftSegment.h.

Constructor & Destructor Documentation

DTVDriftSegment::DTVDriftSegment

(

edm::ParameterSet const &

)

Definition at line 34 of file DTVDriftSegment.cc.

References debug, fitter_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), rootFile_, and CSCSkim_cfi::rootFileName.

      : nSigmas_(pset.getUntrackedParameter<unsigned int>("nSigmasFitRange", 1)) {
    string rootFileName = pset.getParameter<string>("rootFileName");
    rootFile_ = new TFile(rootFileName.c_str(), "READ");

    bool debug = pset.getUntrackedParameter<bool>("debug", false);
    fitter_ = new DTResidualFitter(debug);
    //bool debug = pset.getUntrackedParameter<bool>("debug", false);
    //if(debug) fitter_->setVerbosity(1);
  }

DTVDriftSegment::~DTVDriftSegment ( )

override

Definition at line 45 of file DTVDriftSegment.cc.

References fitter_, and rootFile_.

                                    {
    rootFile_->Close();
    delete fitter_;
  }

Member Function Documentation

DTVDriftData DTVDriftSegment::compute ( const DTSuperLayerId &  )

overridevirtual

Implements dtCalibration::DTVDriftBaseAlgo.

Definition at line 57 of file DTVDriftSegment.cc.

References DTVelocityUnits::cm_per_ns, Exception, DTResidualFitResult::fitMean, DTResidualFitResult::fitMeanError, DTResidualFitter::fitResiduals(), DTResidualFitResult::fitSigma, DTResidualFitResult::fitSigmaError, fitter_, DTMtime::get(), getHisto(), LogTrace, mTimeMap_, nSigmas_, L1TObjectsTimingClient_cff::resolution, mps_update::status, and DTSuperLayerId::superLayer().

                                                                  {
    // Get original value from DB; vdrift is cm/ns , resolution is cm
    float vDrift = 0., resolution = 0.;
    int status = mTimeMap_->get(slId, vDrift, resolution, DTVelocityUnits::cm_per_ns);

    if (status != 0)
      throw cms::Exception("DTCalibration") << "Could not find vDrift entry in DB for" << slId << endl;

    // For RZ superlayers use original value
    if (slId.superLayer() == 2) {
      LogTrace("Calibration") << "[DTVDriftSegment]: RZ superlayer\n"
                              << "                   Will use original vDrift and resolution.";
      return DTVDriftData(vDrift, resolution);
    } else {
      TH1F* vDriftCorrHisto = getHisto(slId);
      // If empty histogram
      if (vDriftCorrHisto->GetEntries() == 0) {
        LogError("Calibration") << "[DTVDriftSegment]: Histogram " << vDriftCorrHisto->GetName() << " is empty.\n"
                                << "                   Will use original vDrift and resolution.";
        return DTVDriftData(vDrift, resolution);
      }

      LogTrace("Calibration") << "[DTVDriftSegment]: Fitting histogram " << vDriftCorrHisto->GetName();
      DTResidualFitResult fitResult = fitter_->fitResiduals(*vDriftCorrHisto, nSigmas_);
      LogTrace("Calibration") << "[DTVDriftSegment]: \n"
                              << "   Fit Mean  = " << fitResult.fitMean << " +/- " << fitResult.fitMeanError << "\n"
                              << "   Fit Sigma = " << fitResult.fitSigma << " +/- " << fitResult.fitSigmaError;

      float vDriftCorr = fitResult.fitMean;
      float vDriftNew = vDrift * (1. - vDriftCorr);
      float resolutionNew = resolution;
      return DTVDriftData(vDriftNew, resolutionNew);
    }
  }

TH1F * DTVDriftSegment::getHisto ( const DTSuperLayerId &  slId )

private

Definition at line 92 of file DTVDriftSegment.cc.

References Exception, getHistoName(), timingPdfMaker::histo, HltBtagPostValidation_cff::histoName, and rootFile_.

Referenced by compute().

                                                            {
    string histoName = getHistoName(slId);
    TH1F* histo = static_cast<TH1F*>(rootFile_->Get(histoName.c_str()));
    if (!histo)
      throw cms::Exception("DTCalibration") << "v-drift correction histogram not found:" << histoName << endl;
    return histo;
  }

string DTVDriftSegment::getHistoName ( const DTSuperLayerId &  slId )

private

Definition at line 100 of file DTVDriftSegment.cc.

References DTSuperLayerId::chamberId(), DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, DTSuperLayerId::superLayer(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by getHisto().

                                                                 {
    DTChamberId chId = slId.chamberId();

    // Compose the chamber name
    std::string wheel = std::to_string(chId.wheel());
    std::string station = std::to_string(chId.station());
    std::string sector = std::to_string(chId.sector());

    string chHistoName = "_W" + wheel + "_St" + station + "_Sec" + sector;

    return (slId.superLayer() != 2) ? ("hRPhiVDriftCorr" + chHistoName) : ("hRZVDriftCorr" + chHistoName);
  }

void DTVDriftSegment::setES ( const edm::EventSetup &  setup )

overridevirtual

Implements dtCalibration::DTVDriftBaseAlgo.

Definition at line 50 of file DTVDriftSegment.cc.

References edm::EventSetup::get(), and mTimeMap_.

                                                        {
    // Get the map of vdrift from the setup
    ESHandle<DTMtime> mTime;
    setup.get<DTMtimeRcd>().get(mTime);
    mTimeMap_ = &*mTime;
  }

Member Data Documentation

DTResidualFitter* dtCalibration::DTVDriftSegment::fitter_

private

Definition at line 39 of file DTVDriftSegment.h.

Referenced by compute(), DTVDriftSegment(), and ~DTVDriftSegment().

const DTMtime* dtCalibration::DTVDriftSegment::mTimeMap_

private

Definition at line 37 of file DTVDriftSegment.h.

Referenced by compute(), and setES().

unsigned int dtCalibration::DTVDriftSegment::nSigmas_

private

Definition at line 35 of file DTVDriftSegment.h.

Referenced by compute().

TFile* dtCalibration::DTVDriftSegment::rootFile_

private

Definition at line 38 of file DTVDriftSegment.h.

Referenced by DTVDriftSegment(), getHisto(), and ~DTVDriftSegment().