#include <DTTTrigResidualCorrection.h>

Inheritance diagram for dtCalibration::DTTTrigResidualCorrection:

Public Member Functions
DTTTrigData	correction (const DTSuperLayerId &) override

	DTTTrigResidualCorrection (const edm::ParameterSet &)

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

	~DTTTrigResidualCorrection () override

Public Member Functions inherited from dtCalibration::DTTTrigBaseCorrection
	DTTTrigBaseCorrection ()

virtual	~DTTTrigBaseCorrection ()

Private Member Functions
const TH1F *	getHisto (const DTSuperLayerId &)

std::string	getHistoName (const DTSuperLayerId &slID)

Private Attributes
std::string	dbLabel_

DTResidualFitter *	fitter_

const DTMtime *	mTimeMap_

std::string	rootBaseDir_

TFile *	rootFile_

const DTTtrig *	tTrigMap_

bool	useFit_

bool	useSlopesCalib_

double	vDriftEff_ [5][14][4][3]

Detailed Description

Definition at line 28 of file DTTTrigResidualCorrection.h.

Constructor & Destructor Documentation

DTTTrigResidualCorrection::DTTTrigResidualCorrection ( const edm::ParameterSet & pset )

Definition at line 43 of file DTTTrigResidualCorrection.cc.

References debug, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

                                                                              {
   string residualsRootFile = pset.getParameter<string>("residualsRootFile");
   rootFile_ = new TFile(residualsRootFile.c_str(),"READ");
   rootBaseDir_ = pset.getUntrackedParameter<string>("rootBaseDir","/DQMData/DT/DTCalibValidation"); 
   useFit_ = pset.getParameter<bool>("useFitToResiduals");
   //useConstantvDrift_ = pset.getParameter<bool>("useConstantDriftVelocity");
   dbLabel_  = pset.getUntrackedParameter<string>("dbLabel", "");
   useSlopesCalib_ = pset.getUntrackedParameter<bool>("useSlopesCalib",false);
 
   // Load external slopes
   if(useSlopesCalib_){ 
      ifstream fileSlopes; 
      fileSlopes.open( pset.getParameter<FileInPath>("slopesFileName").fullPath().c_str() );
 
      int tmp_wheel = 0, tmp_sector = 0, tmp_station = 0, tmp_SL = 0;
      double tmp_ttrig = 0., tmp_t0 = 0., tmp_kfact = 0.;
      int tmp_a = 0, tmp_b = 0, tmp_c = 0, tmp_d = 0;
      double tmp_v_eff = 0.;
      while(!fileSlopes.eof()){
         fileSlopes >> tmp_wheel >> tmp_sector >> tmp_station >> tmp_SL  >> tmp_a >> tmp_b >>
                       tmp_ttrig >> tmp_t0 >> tmp_kfact >> tmp_c >> tmp_d >> tmp_v_eff;
         vDriftEff_[tmp_wheel+2][tmp_sector-1][tmp_station-1][tmp_SL-1] = -tmp_v_eff;
      }
      fileSlopes.close();
   } 
 
   bool debug = pset.getUntrackedParameter<bool>("debug",false);
   fitter_ = new DTResidualFitter(debug);
 }

DTTTrigResidualCorrection::~DTTTrigResidualCorrection ( )

override

Definition at line 73 of file DTTTrigResidualCorrection.cc.

                                                       {
   delete rootFile_;
   delete fitter_;
 }

Member Function Documentation

DTTTrigData DTTTrigResidualCorrection::correction ( const DTSuperLayerId & slId )

overridevirtual

Implements dtCalibration::DTTTrigBaseCorrection.

Definition at line 91 of file DTTTrigResidualCorrection.cc.

References DTVelocityUnits::cm_per_ns, Exception, DTResidualFitResult::fitMean, DTResidualFitResult::fitSigma, DTLinearDriftAlgo_cfi::hitResolution, DTAnalyzerDetailed_cfi::kFactor, LogTrace, DTTimeUnits::ns, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, mps_update::status, DTSuperLayerId::superLayer(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                             {
 
   float tTrigMean,tTrigSigma,kFactor;
   int status = tTrigMap_->get(slId,tTrigMean,tTrigSigma,kFactor,DTTimeUnits::ns);
   if(status != 0) throw cms::Exception("[DTTTrigResidualCorrection]") << "Could not find tTrig entry in DB for"
                                                                       << slId << endl;
 
   float vDrift,hitResolution;
   status = mTimeMap_->get(slId,vDrift,hitResolution,DTVelocityUnits::cm_per_ns);
   if(status != 0) throw cms::Exception("[DTTTrigResidualCorrection]") << "Could not find vDrift entry in DB for"
                                                                       << slId << endl;
   TH1F residualHisto = *(getHisto(slId));
   LogTrace("Calibration") << "[DTTTrigResidualCorrection]: \n"
                           << "   Mean, RMS     = " << residualHisto.GetMean() << ", " << residualHisto.GetRMS();
 
   double fitMean = -1.;
   double fitSigma = -1.;
   if(useFit_){
     LogTrace("Calibration") << "[DTTTrigResidualCorrection]: Fitting histogram " << residualHisto.GetName(); 
     int nSigmas = 2;
 
     DTResidualFitResult fitResult = fitter_->fitResiduals(residualHisto,nSigmas);
     fitMean = fitResult.fitMean;
     fitSigma = fitResult.fitSigma;
 
     LogTrace("Calibration") << "[DTTTrigResidualCorrection]: \n"
                             << "   Fit Mean      = " << fitMean << "\n"
                             << "   Fit Sigma     = " << fitSigma;
   }
   double resMean = (useFit_) ? fitMean : residualHisto.GetMean();
 
   int wheel = slId.wheel();
   int sector = slId.sector();
   int station = slId.station();
   int superLayer = slId.superLayer();
   double resTime = 0.;
   if(useSlopesCalib_){
      double vdrift_eff = vDriftEff_[wheel+2][sector-1][station-1][superLayer-1];
      if(vdrift_eff == 0) vdrift_eff = vDrift;
 
      if(vdrift_eff) resTime = resMean/vdrift_eff;
 
      LogTrace("Calibration") << "[DTTTrigResidualCorrection]: Effective vDrift, correction to tTrig = "
                              << vdrift_eff << ", " << resTime;
   } else{
      if(vDrift) resTime = resMean/vDrift;
 
      LogTrace("Calibration") << "[DTTTrigResidualCorrection]: vDrift from DB, correction to tTrig = "
                              << vDrift << ", " << resTime;
   }
  
   double corrMean = tTrigMean;
   double corrSigma = (tTrigSigma != 0.) ? tTrigSigma : 1.;
   double corrKFact = (tTrigSigma != 0.) ? (kFactor + resTime/tTrigSigma) : resTime;
 
   return DTTTrigData(corrMean,corrSigma,corrKFact);  
 }

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

private

Definition at line 149 of file DTTTrigResidualCorrection.cc.

References Exception, mergeVDriftHistosByStation::getHistoName(), trackerHits::histo, and LogTrace.

                                                                           {
   string histoName = getHistoName(slId);
   LogTrace("Calibration") << "[DTTTrigResidualCorrection]: Accessing histogram " << histoName.c_str();
   TH1F* histo = static_cast<TH1F*>(rootFile_->Get(histoName.c_str()));
   if(!histo) throw cms::Exception("[DTTTrigResidualCorrection]") << "residual histogram not found:"
                                                                  << histoName << endl; 
   return histo;
 }

string DTTTrigResidualCorrection::getHistoName ( const DTSuperLayerId & slID )

private

Definition at line 158 of file DTTTrigResidualCorrection.cc.

References DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, cmsPerfCommons::Step, AlCaHLTBitMon_QueryRunRegistry::string, DTSuperLayerId::superlayer(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                          {
 
   int step = 3;
 
   std::string wheel = std::to_string(slId.wheel());
   std::string station = std::to_string(slId.station());
   std::string sector = std::to_string(slId.sector());
   std::string superLayer = std::to_string(slId.superlayer());
   std::string Step = std::to_string(step);
 
   string histoName =
     rootBaseDir_ + "/Wheel" + wheel + "/Station" + station + "/Sector" + sector + 
     "/hResDist_STEP" + Step + "_W" + wheel + "_St" + station + "_Sec" + sector + "_SL" + superLayer;
 
   return histoName;
 }

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

overridevirtual

Implements dtCalibration::DTTTrigBaseCorrection.

Definition at line 78 of file DTTTrigResidualCorrection.cc.

References edm::EventSetup::get().

                                                              {
   // Get tTrig record from DB
   ESHandle<DTTtrig> tTrig;
   //setup.get<DTTtrigRcd>().get(tTrig);
   setup.get<DTTtrigRcd>().get(dbLabel_,tTrig);
   tTrigMap_ = &*tTrig;
 
   // Get vDrift record
   ESHandle<DTMtime> mTimeHandle;
   setup.get<DTMtimeRcd>().get(mTimeHandle);
   mTimeMap_ = &*mTimeHandle;
 }