dtCalibration::DTVDriftSegment Class Reference

#include <DTVDriftSegment.h>

Inheritance diagram for dtCalibration::DTVDriftSegment:

Public Member Functions
DTVDriftData	compute (const DTSuperLayerId &) override
	DTVDriftSegment (edm::ParameterSet const &, edm::ConsumesCollector)
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_
edm::ESGetToken< DTMtime, DTMtimeRcd >	mTimeMapToken_
unsigned int	nSigmas_
bool	readLegacyVDriftDB
TFile *	rootFile_
const DTRecoConditions *	vDriftMap_
edm::ESGetToken< DTRecoConditions, DTRecoConditionsVdriftRcd >	vDriftMapToken_

Detailed Description

Definition at line 30 of file DTVDriftSegment.h.

Constructor & Destructor Documentation

DTVDriftSegment()

DTVDriftSegment::DTVDriftSegment	(	edm::ParameterSet const &	pset,
		edm::ConsumesCollector	cc
	)

Definition at line 36 of file DTVDriftSegment.cc.

References edm::BeginRun, gpuPixelDoublets::cc, debug, fitter_, mTimeMapToken_, muonDTDigis_cfi::pset, readLegacyVDriftDB, rootFile_, CSCSkim_cfi::rootFileName, and vDriftMapToken_.

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

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

    readLegacyVDriftDB = pset.getParameter<bool>("readLegacyVDriftDB");
    if (readLegacyVDriftDB) {
      mTimeMapToken_ = cc.esConsumes<edm::Transition::BeginRun>();
    } else {
      vDriftMapToken_ = cc.esConsumes<edm::Transition::BeginRun>();
    }
  }

~DTVDriftSegment()

DTVDriftSegment::~DTVDriftSegment ( )

override

Definition at line 55 of file DTVDriftSegment.cc.

References fitter_, and rootFile_.

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

Member Function Documentation

compute()

DTVDriftData DTVDriftSegment::compute ( const DTSuperLayerId &  slId )

overridevirtual

Implements dtCalibration::DTVDriftBaseAlgo.

Definition at line 74 of file DTVDriftSegment.cc.

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

                                                                  {
    // Get original value from DB; vdrift is cm/ns , resolution is cm
    // Note that resolution is irrelevant as it is no longer used anywhere in reconstruction.

    float vDrift = 0., resolution = 0.;
    if (readLegacyVDriftDB) {  // Legacy format
      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;
    } else {  // New DB format
      vDrift = vDriftMap_->get(DTWireId(slId.rawId()));
    }

    // 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);
    }
  }

getHisto()

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

private

Definition at line 114 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;
  }

getHistoName()

string DTVDriftSegment::getHistoName ( const DTSuperLayerId &  slId )

private

Definition at line 122 of file DTVDriftSegment.cc.

References DTSuperLayerId::chamberId(), hgcalTBTopologyTester_cfi::sector, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, DTSuperLayerId::superLayer(), cond::impl::to_string(), 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);
  }

setES()

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

overridevirtual

Implements dtCalibration::DTVDriftBaseAlgo.

Definition at line 60 of file DTVDriftSegment.cc.

References Exception, mTimeMap_, mTimeMapToken_, readLegacyVDriftDB, singleTopDQM_cfi::setup, vDriftMap_, vDriftMapToken_, DTRecoConditions::version(), and BeamSplash_cfg::version.

                                                        {
    // Get the map of vdrift from the setup
    if (readLegacyVDriftDB) {
      mTimeMap_ = &setup.getData(mTimeMapToken_);
    } else {
      vDriftMap_ = &setup.getData(vDriftMapToken_);
      // Consistency check: no parametrization is implemented for the time being
      int version = vDriftMap_->version();
      if (version != 1) {
        throw cms::Exception("Configuration") << "only version 1 is presently supported for VDriftDB";
      }
    }
  }

Member Data Documentation

fitter_

DTResidualFitter* dtCalibration::DTVDriftSegment::fitter_

private

Definition at line 48 of file DTVDriftSegment.h.

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

mTimeMap_

const DTMtime* dtCalibration::DTVDriftSegment::mTimeMap_

private

Definition at line 44 of file DTVDriftSegment.h.

Referenced by compute(), and setES().

mTimeMapToken_

edm::ESGetToken<DTMtime, DTMtimeRcd> dtCalibration::DTVDriftSegment::mTimeMapToken_

private

Definition at line 50 of file DTVDriftSegment.h.

Referenced by DTVDriftSegment(), and setES().

nSigmas_

unsigned int dtCalibration::DTVDriftSegment::nSigmas_

private

Definition at line 42 of file DTVDriftSegment.h.

Referenced by compute().

readLegacyVDriftDB

bool dtCalibration::DTVDriftSegment::readLegacyVDriftDB

private

Definition at line 46 of file DTVDriftSegment.h.

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

rootFile_

TFile* dtCalibration::DTVDriftSegment::rootFile_

private

Definition at line 47 of file DTVDriftSegment.h.

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

vDriftMap_

const DTRecoConditions* dtCalibration::DTVDriftSegment::vDriftMap_

private

Definition at line 45 of file DTVDriftSegment.h.

Referenced by compute(), and setES().

vDriftMapToken_

edm::ESGetToken<DTRecoConditions, DTRecoConditionsVdriftRcd> dtCalibration::DTVDriftSegment::vDriftMapToken_

private

Definition at line 51 of file DTVDriftSegment.h.

Referenced by DTVDriftSegment(), and setES().