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_
bool	readLegacyVDriftDB
TFile *	rootFile_
const DTRecoConditions *	vDriftMap_

Detailed Description

Definition at line 24 of file DTVDriftSegment.h.

Constructor & Destructor Documentation

DTVDriftSegment()

DTVDriftSegment::DTVDriftSegment

(

edm::ParameterSet const &

pset

)

Definition at line 36 of file DTVDriftSegment.cc.

      : 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");
  }

References debug, fitter_, muonDTDigis_cfi::pset, readLegacyVDriftDB, rootFile_, and CSCSkim_cfi::rootFileName.

~DTVDriftSegment()

DTVDriftSegment::~DTVDriftSegment ( )

override

Definition at line 50 of file DTVDriftSegment.cc.

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

References fitter_, and rootFile_.

Member Function Documentation

compute()

DTVDriftData DTVDriftSegment::compute ( const DTSuperLayerId &  slId )

overridevirtual

Implements dtCalibration::DTVDriftBaseAlgo.

Definition at line 73 of file DTVDriftSegment.cc.

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

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_.

getHisto()

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

private

Definition at line 113 of file DTVDriftSegment.cc.

                                                            {
    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;
  }

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

Referenced by compute().

getHistoName()

string DTVDriftSegment::getHistoName ( const DTSuperLayerId &  slId )

private

Definition at line 121 of file DTVDriftSegment.cc.

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

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

Referenced by getHisto().

setES()

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

overridevirtual

Implements dtCalibration::DTVDriftBaseAlgo.

Definition at line 55 of file DTVDriftSegment.cc.

                                                        {
    // Get the map of vdrift from the setup
    if (readLegacyVDriftDB) {
      ESHandle<DTMtime> mTime;
      setup.get<DTMtimeRcd>().get(mTime);
      mTimeMap_ = &*mTime;
    } else {
      ESHandle<DTRecoConditions> hVdrift;
      setup.get<DTRecoConditionsVdriftRcd>().get(hVdrift);
      vDriftMap_ = &*hVdrift;
      // 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";
      }
    }
  }

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

Member Data Documentation

fitter_

DTResidualFitter* dtCalibration::DTVDriftSegment::fitter_

private

Definition at line 42 of file DTVDriftSegment.h.

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

mTimeMap_

const DTMtime* dtCalibration::DTVDriftSegment::mTimeMap_

private

Definition at line 38 of file DTVDriftSegment.h.

Referenced by compute(), and setES().

nSigmas_

unsigned int dtCalibration::DTVDriftSegment::nSigmas_

private

Definition at line 36 of file DTVDriftSegment.h.

Referenced by compute().

readLegacyVDriftDB

bool dtCalibration::DTVDriftSegment::readLegacyVDriftDB

private

Definition at line 40 of file DTVDriftSegment.h.

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

rootFile_

TFile* dtCalibration::DTVDriftSegment::rootFile_

private

Definition at line 41 of file DTVDriftSegment.h.

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

vDriftMap_

const DTRecoConditions* dtCalibration::DTVDriftSegment::vDriftMap_

private

Definition at line 39 of file DTVDriftSegment.h.

Referenced by compute(), and setES().