#include <DTResidualCalibration.h>

Inheritance diagram for DTResidualCalibration:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

void	beginJob () override

void	beginRun (const edm::Run &, const edm::EventSetup &) override

	DTResidualCalibration (const edm::ParameterSet &pset)
	Constructor. More...

void	endJob () override

void	endRun (const edm::Run &, const edm::EventSetup &) override

	~DTResidualCalibration () override
	Destructor. More...

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	bookHistos (DTSuperLayerId slId)

void	bookHistos (DTLayerId slId)

void	fillHistos (DTSuperLayerId slId, float distance, float residualOnDistance)

void	fillHistos (DTLayerId slId, float distance, float residualOnDistance)

float	segmentToWireDistance (const DTRecHit1D &recHit1D, const DTRecSegment4D &segment)

Private Attributes
const bool	detailedAnalysis_

const DTGeometry *	dtGeom_

const edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeomToken_

std::map< DTLayerId, TH1F * >	histoMapPerLayerTH1F_

std::map< DTLayerId, TH2F * >	histoMapPerLayerTH2F_

std::map< DTSuperLayerId, TH1F * >	histoMapTH1F_

std::map< DTSuperLayerId, TH2F * >	histoMapTH2F_

const double	histRange_

unsigned int	nevent

const std::string	rootBaseDir_

TFile *	rootFile_

unsigned int	segmbad

const edm::EDGetTokenT< DTRecSegment4DCollection >	segment4DToken_

unsigned int	segmok

DTSegmentSelector *	select_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Extracts DT segment residuals

Definition at line 30 of file DTResidualCalibration.h.

Constructor & Destructor Documentation

◆ DTResidualCalibration()

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

Constructor.

Definition at line 34 of file DTResidualCalibration.cc.

References edm::EDConsumerBase::consumesCollector(), LogDebug, nevent, muonDTDigis_cfi::pset, rootFile_, CSCSkim_cfi::rootFileName, segmbad, segmok, select_, and AlCaHLTBitMon_QueryRunRegistry::string.

     : histRange_(pset.getParameter<double>("histogramRange")),
       segment4DToken_(consumes<DTRecSegment4DCollection>(pset.getParameter<edm::InputTag>("segment4DLabel"))),
       rootBaseDir_(pset.getUntrackedParameter<std::string>("rootBaseDir", "DT/Residuals")),
       detailedAnalysis_(pset.getUntrackedParameter<bool>("detailedAnalysis", false)),
       dtGeomToken_(esConsumes<edm::Transition::BeginRun>()) {
   edm::ConsumesCollector collector(consumesCollector());
   select_ = new DTSegmentSelector(pset, collector);
 
   LogDebug("Calibration") << "[DTResidualCalibration] Constructor called.";
   std::string rootFileName = pset.getUntrackedParameter<std::string>("rootFileName", "residuals.root");
   rootFile_ = new TFile(rootFileName.c_str(), "RECREATE");
   rootFile_->cd();
 
   segmok = 0;
   segmbad = 0;
   nevent = 0;
 }

◆ ~DTResidualCalibration()

DTResidualCalibration::~DTResidualCalibration ( )

override

Destructor.

Definition at line 53 of file DTResidualCalibration.cc.

References nevent, segmbad, segmok, and select_.

                                               {
   delete select_;
   edm::LogVerbatim("Calibration") << "[DTResidualCalibration] Destructor called.";
   edm::LogVerbatim("Calibration") << "[DTResidualCalibration] Analyzed events: " << nevent;
   edm::LogVerbatim("Calibration") << "[DTResidualCalibration] Good segments: " << segmok;
   edm::LogVerbatim("Calibration") << "[DTResidualCalibration] Bad segments: " << segmbad;
 }

Member Function Documentation

◆ analyze()

void DTResidualCalibration::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 87 of file DTResidualCalibration.cc.

References relativeConstraints::chamber, DTGeometry::chamber(), DTRecHitSegmentResidual::compute(), filterCSVwithJSON::copy, detailedAnalysis_, dtGeom_, edmPickEvents::event, fillHistos(), DTWireId::layerId(), LogTrace, nevent, FastTimerService_cff::range, rootFile_, segmbad, segment4DToken_, segmentToWireDistance(), segmok, select_, singleTopDQM_cfi::setup, DTRecSegment2D::specificRecHits(), and DTLayerId::superlayerId().

                                                                                      {
   rootFile_->cd();
   ++nevent;
 
   // Get the 4D rechits from the event
   const edm::Handle<DTRecSegment4DCollection>& segments4D = event.getHandle(segment4DToken_);
 
   // Loop over segments by chamber
   DTRecSegment4DCollection::id_iterator chamberIdIt;
   for (chamberIdIt = segments4D->id_begin(); chamberIdIt != segments4D->id_end(); ++chamberIdIt) {
     const DTChamber* chamber = dtGeom_->chamber(*chamberIdIt);
 
     // Get the range for the corresponding ChamberId
     DTRecSegment4DCollection::range range = segments4D->get((*chamberIdIt));
     // Loop over the rechits of this DetUnit
     for (DTRecSegment4DCollection::const_iterator segment = range.first; segment != range.second; ++segment) {
       LogTrace("Calibration") << "Segment local pos (in chamber RF): " << (*segment).localPosition()
                               << "\nSegment global pos: " << chamber->toGlobal((*segment).localPosition());
 
       if (!(*select_)(*segment, event, setup)) {
         segmbad++;
         continue;
       }
       segmok++;
 
       // Get all 1D RecHits at step 3 within the 4D segment
       std::vector<DTRecHit1D> recHits1D_S3;
 
       if ((*segment).hasPhi()) {
         const DTChamberRecSegment2D* phiSeg = (*segment).phiSegment();
         const std::vector<DTRecHit1D>& phiRecHits = phiSeg->specificRecHits();
         std::copy(phiRecHits.begin(), phiRecHits.end(), back_inserter(recHits1D_S3));
       }
 
       if ((*segment).hasZed()) {
         const DTSLRecSegment2D* zSeg = (*segment).zSegment();
         const std::vector<DTRecHit1D>& zRecHits = zSeg->specificRecHits();
         std::copy(zRecHits.begin(), zRecHits.end(), back_inserter(recHits1D_S3));
       }
 
       // Loop over 1D RecHit inside 4D segment
       for (std::vector<DTRecHit1D>::const_iterator recHit1D = recHits1D_S3.begin(); recHit1D != recHits1D_S3.end();
            ++recHit1D) {
         const DTWireId wireId = recHit1D->wireId();
 
         float segmDistance = segmentToWireDistance(*recHit1D, *segment);
         if (segmDistance > 2.1)
           LogTrace("Calibration") << "WARNING: segment-wire distance: " << segmDistance;
         else
           LogTrace("Calibration") << "segment-wire distance: " << segmDistance;
 
         float residualOnDistance = DTRecHitSegmentResidual().compute(dtGeom_, *recHit1D, *segment);
         LogTrace("Calibration") << "Wire Id " << wireId << " residual on distance: " << residualOnDistance;
 
         fillHistos(wireId.superlayerId(), segmDistance, residualOnDistance);
         if (detailedAnalysis_)
           fillHistos(wireId.layerId(), segmDistance, residualOnDistance);
       }
     }
   }
 }

◆ beginJob()

void DTResidualCalibration::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 61 of file DTResidualCalibration.cc.

61 { TH1::SetDefaultSumw2(true); }

◆ beginRun()

void DTResidualCalibration::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

override

Definition at line 63 of file DTResidualCalibration.cc.

References bookHistos(), DTGeometry::chambers(), detailedAnalysis_, dtGeom_, dtGeomToken_, histoMapTH1F_, hgcalTopologyTester_cfi::layers, edm::ESHandle< T >::product(), and singleTopDQM_cfi::setup.

                                                                                   {
   // get the geometry
   edm::ESHandle<DTGeometry> dtGeomH;
   dtGeomH = setup.getHandle(dtGeomToken_);
   dtGeom_ = dtGeomH.product();
 
   // Loop over all the chambers
   if (histoMapTH1F_.empty()) {
     for (auto ch_it : dtGeom_->chambers()) {
       // Loop over the SLs
       for (auto sl_it : ch_it->superLayers()) {
         DTSuperLayerId slId = (sl_it)->id();
         bookHistos(slId);
         if (detailedAnalysis_) {
           for (auto layer_it : (sl_it)->layers()) {
             DTLayerId layerId = (layer_it)->id();
             bookHistos(layerId);
           }
         }
       }
     }
   }
 }

◆ bookHistos() [1/2]

void DTResidualCalibration::bookHistos ( DTSuperLayerId slId )

private

Definition at line 197 of file DTResidualCalibration.cc.

References ALCARECODTCalibSynchDQM_cff::baseDir, histoMapTH1F_, histoMapTH2F_, histRange_, LogDebug, rootBaseDir_, rootFile_, DTChamberId::sector(), DTChamberId::station(), AlCaHLTBitMon_QueryRunRegistry::string, DTSuperLayerId::superlayer(), cond::impl::to_string(), and DTChamberId::wheel().

Referenced by beginRun().

                                                           {
   TH1AddDirectorySentry addDir;
   rootFile_->cd();
 
   LogDebug("Calibration") << "[DTResidualCalibration] Booking histos for SL: " << slId;
 
   // Compose the chamber name
   // Define the step
   int step = 3;
 
   std::string wheelStr = std::to_string(slId.wheel());
   std::string stationStr = std::to_string(slId.station());
   std::string sectorStr = std::to_string(slId.sector());
 
   std::string slHistoName = "_STEP" + std::to_string(step) + "_W" + wheelStr + "_St" + stationStr + "_Sec" + sectorStr +
                             "_SL" + std::to_string(slId.superlayer());
 
   LogDebug("Calibration") << "Accessing " << rootBaseDir_;
   TDirectory* baseDir = rootFile_->GetDirectory(rootBaseDir_.c_str());
   if (!baseDir)
     baseDir = rootFile_->mkdir(rootBaseDir_.c_str());
   LogDebug("Calibration") << "Accessing " << ("Wheel" + wheelStr);
   TDirectory* wheelDir = baseDir->GetDirectory(("Wheel" + wheelStr).c_str());
   if (!wheelDir)
     wheelDir = baseDir->mkdir(("Wheel" + wheelStr).c_str());
   LogDebug("Calibration") << "Accessing " << ("Station" + stationStr);
   TDirectory* stationDir = wheelDir->GetDirectory(("Station" + stationStr).c_str());
   if (!stationDir)
     stationDir = wheelDir->mkdir(("Station" + stationStr).c_str());
   LogDebug("Calibration") << "Accessing " << ("Sector" + sectorStr);
   TDirectory* sectorDir = stationDir->GetDirectory(("Sector" + sectorStr).c_str());
   if (!sectorDir)
     sectorDir = stationDir->mkdir(("Sector" + sectorStr).c_str());
 
   sectorDir->cd();
 
   // Create the monitor elements
   TH1F* histosTH1F = new TH1F(("hResDist" + slHistoName).c_str(),
                               "Residuals on the distance from wire (rec_hit - segm_extr) (cm)",
                               200,
                               -histRange_,
                               histRange_);
   TH2F* histosTH2F = new TH2F(("hResDistVsDist" + slHistoName).c_str(),
                               "Residuals on the dist. (cm) from wire (rec_hit - segm_extr) vs dist. (cm)",
                               100,
                               0,
                               2.5,
                               200,
                               -histRange_,
                               histRange_);
   histoMapTH1F_[slId] = histosTH1F;
   histoMapTH2F_[slId] = histosTH2F;
 }

◆ bookHistos() [2/2]

void DTResidualCalibration::bookHistos ( DTLayerId slId )

private

Definition at line 251 of file DTResidualCalibration.cc.

References ALCARECODTCalibSynchDQM_cff::baseDir, histoMapPerLayerTH1F_, histoMapPerLayerTH2F_, histRange_, DTLayerId::layer(), LogDebug, rootBaseDir_, rootFile_, DTChamberId::sector(), DTChamberId::station(), AlCaHLTBitMon_QueryRunRegistry::string, DTSuperLayerId::superlayer(), cond::impl::to_string(), and DTChamberId::wheel().

                                                         {
   TH1AddDirectorySentry addDir;
   rootFile_->cd();
 
   LogDebug("Calibration") << "[DTResidualCalibration] Booking histos for layer: " << layerId;
 
   // Compose the chamber name
   std::string wheelStr = std::to_string(layerId.wheel());
   std::string stationStr = std::to_string(layerId.station());
   std::string sectorStr = std::to_string(layerId.sector());
   std::string superLayerStr = std::to_string(layerId.superlayer());
   std::string layerStr = std::to_string(layerId.layer());
   // Define the step
   int step = 3;
 
   std::string layerHistoName = "_STEP" + std::to_string(step) + "_W" + wheelStr + "_St" + stationStr + "_Sec" +
                                sectorStr + "_SL" + superLayerStr + "_Layer" + layerStr;
 
   LogDebug("Calibration") << "Accessing " << rootBaseDir_;
   TDirectory* baseDir = rootFile_->GetDirectory(rootBaseDir_.c_str());
   if (!baseDir)
     baseDir = rootFile_->mkdir(rootBaseDir_.c_str());
   LogDebug("Calibration") << "Accessing " << ("Wheel" + wheelStr);
   TDirectory* wheelDir = baseDir->GetDirectory(("Wheel" + wheelStr).c_str());
   if (!wheelDir)
     wheelDir = baseDir->mkdir(("Wheel" + wheelStr).c_str());
   LogDebug("Calibration") << "Accessing " << ("Station" + stationStr);
   TDirectory* stationDir = wheelDir->GetDirectory(("Station" + stationStr).c_str());
   if (!stationDir)
     stationDir = wheelDir->mkdir(("Station" + stationStr).c_str());
   LogDebug("Calibration") << "Accessing " << ("Sector" + sectorStr);
   TDirectory* sectorDir = stationDir->GetDirectory(("Sector" + sectorStr).c_str());
   if (!sectorDir)
     sectorDir = stationDir->mkdir(("Sector" + sectorStr).c_str());
   LogDebug("Calibration") << "Accessing " << ("SL" + superLayerStr);
   TDirectory* superLayerDir = sectorDir->GetDirectory(("SL" + superLayerStr).c_str());
   if (!superLayerDir)
     superLayerDir = sectorDir->mkdir(("SL" + superLayerStr).c_str());
 
   superLayerDir->cd();
   // Create histograms
   TH1F* histosTH1F = new TH1F(("hResDist" + layerHistoName).c_str(),
                               "Residuals on the distance from wire (rec_hit - segm_extr) (cm)",
                               200,
                               -histRange_,
                               histRange_);
   TH2F* histosTH2F = new TH2F(("hResDistVsDist" + layerHistoName).c_str(),
                               "Residuals on the dist. (cm) from wire (rec_hit - segm_extr) vs dist. (cm)",
                               100,
                               0,
                               2.5,
                               200,
                               -histRange_,
                               histRange_);
   histoMapPerLayerTH1F_[layerId] = histosTH1F;
   histoMapPerLayerTH2F_[layerId] = histosTH2F;
 }

◆ endJob()

void DTResidualCalibration::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 178 of file DTResidualCalibration.cc.

References LogDebug, and rootFile_.

Referenced by o2olib.O2ORunMgr::executeJob().

                                    {
   LogDebug("Calibration") << "[DTResidualCalibration] Writing histos to file.";
   rootFile_->cd();
   rootFile_->Write();
   rootFile_->Close();
 
   /*std::map<DTSuperLayerId, TH1F* >::const_iterator itSlHistos = histoMapTH1F_.begin();
   std::map<DTSuperLayerId, TH1F* >::const_iterator itSlHistos_end = histoMapTH1F_.end(); 
   for(; itSlHistos != itSlHistos_end; ++itSlHistos){
      std::vector<TH1F*>::const_iterator itHistTH1F = (*itSlHistos).second.begin();
      std::vector<TH1F*>::const_iterator itHistTH1F_end = (*itSlHistos).second.end();
      for(; itHistTH1F != itHistTH1F_end; ++itHistTH1F) (*itHistTH1F)->Write();
 
      std::vector<TH2F*>::const_iterator itHistTH2F = histoMapTH2F_[(*itSlHistos).first].begin();
      std::vector<TH2F*>::const_iterator itHistTH2F_end = histoMapTH2F_[(*itSlHistos).first].end();
      for(; itHistTH2F != itHistTH2F_end; ++itHistTH2F) (*itHistTH2F)->Write();
   }*/
 }

◆ endRun()

void DTResidualCalibration::endRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

inlineoverride

Definition at line 40 of file DTResidualCalibration.h.

40 {};

◆ fillHistos() [1/2]

void DTResidualCalibration::fillHistos	(	DTSuperLayerId	slId,
		float	distance,
		float	residualOnDistance
	)

private

Definition at line 310 of file DTResidualCalibration.cc.

References HLT_2022v15_cff::distance, histoMapTH1F_, and histoMapTH2F_.

Referenced by analyze().

                                                                                                     {
   histoMapTH1F_[slId]->Fill(residualOnDistance);
   histoMapTH2F_[slId]->Fill(distance, residualOnDistance);
 }

◆ fillHistos() [2/2]

void DTResidualCalibration::fillHistos	(	DTLayerId	slId,
		float	distance,
		float	residualOnDistance
	)

private

Definition at line 316 of file DTResidualCalibration.cc.

References HLT_2022v15_cff::distance, histoMapPerLayerTH1F_, and histoMapPerLayerTH2F_.

                                                                                                   {
   histoMapPerLayerTH1F_[layerId]->Fill(residualOnDistance);
   histoMapPerLayerTH2F_[layerId]->Fill(distance, residualOnDistance);
 }

◆ segmentToWireDistance()

float DTResidualCalibration::segmentToWireDistance	(	const DTRecHit1D &	recHit1D,
		const DTRecSegment4D &	segment
	)

private

Definition at line 149 of file DTResidualCalibration.cc.

References relativeConstraints::chamber, DTGeometry::chamber(), DTSuperLayerId::chamberId(), funct::cos(), dtGeom_, DTGeometry::layer(), phase1PixelTopology::layer, DTWireId::layerId(), DTRecSegment4D::localDirection(), DTRecHit1D::localPosition(), DTRecSegment4D::localPosition(), DTSuperLayerId::superlayer(), PV3DBase< T, PVType, FrameType >::theta(), DTWireId::wire(), DTRecHit1D::wireId(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

                                                                                                             {
   // Get the layer and the wire position
   const DTWireId wireId = recHit1D.wireId();
   const DTLayer* layer = dtGeom_->layer(wireId);
   float wireX = layer->specificTopology().wirePosition(wireId.wire());
 
   // Extrapolate the segment to the z of the wire
   // Get wire position in chamber RF
   // (y and z must be those of the hit to be coherent in the transf. of RF in case of rotations of the layer alignment)
   LocalPoint wirePosInLay(wireX, recHit1D.localPosition().y(), recHit1D.localPosition().z());
   GlobalPoint wirePosGlob = layer->toGlobal(wirePosInLay);
   const DTChamber* chamber = dtGeom_->chamber(wireId.layerId().chamberId());
   LocalPoint wirePosInChamber = chamber->toLocal(wirePosGlob);
 
   // Segment position at Wire z in chamber local frame
   LocalPoint segPosAtZWire =
       segment.localPosition() + segment.localDirection() * wirePosInChamber.z() / cos(segment.localDirection().theta());
 
   // Compute the distance of the segment from the wire
   int sl = wireId.superlayer();
   float segmDistance = -1;
   if (sl == 1 || sl == 3)
     segmDistance = fabs(wirePosInChamber.x() - segPosAtZWire.x());
   else if (sl == 2)
     segmDistance = fabs(segPosAtZWire.y() - wirePosInChamber.y());
 
   return segmDistance;
 }

Member Data Documentation

◆ detailedAnalysis_

const bool DTResidualCalibration::detailedAnalysis_

private

Definition at line 61 of file DTResidualCalibration.h.

Referenced by analyze(), and beginRun().

◆ dtGeom_

const DTGeometry* DTResidualCalibration::dtGeom_

private

Definition at line 64 of file DTResidualCalibration.h.

Referenced by analyze(), beginRun(), and segmentToWireDistance().

◆ dtGeomToken_

const edm::ESGetToken<DTGeometry, MuonGeometryRecord> DTResidualCalibration::dtGeomToken_

private

Definition at line 65 of file DTResidualCalibration.h.

Referenced by beginRun().

◆ histoMapPerLayerTH1F_

std::map<DTLayerId, TH1F*> DTResidualCalibration::histoMapPerLayerTH1F_

private

Definition at line 71 of file DTResidualCalibration.h.

Referenced by bookHistos(), and fillHistos().

◆ histoMapPerLayerTH2F_

std::map<DTLayerId, TH2F*> DTResidualCalibration::histoMapPerLayerTH2F_

private

Definition at line 72 of file DTResidualCalibration.h.

Referenced by bookHistos(), and fillHistos().

◆ histoMapTH1F_

std::map<DTSuperLayerId, TH1F*> DTResidualCalibration::histoMapTH1F_

private

Definition at line 68 of file DTResidualCalibration.h.

Referenced by beginRun(), bookHistos(), and fillHistos().

◆ histoMapTH2F_

std::map<DTSuperLayerId, TH2F*> DTResidualCalibration::histoMapTH2F_

private

Definition at line 69 of file DTResidualCalibration.h.

Referenced by bookHistos(), and fillHistos().

◆ histRange_

const double DTResidualCalibration::histRange_

private

Definition at line 57 of file DTResidualCalibration.h.

Referenced by bookHistos().

◆ nevent

unsigned int DTResidualCalibration::nevent

private

Definition at line 45 of file DTResidualCalibration.h.

Referenced by analyze(), DTResidualCalibration(), and ~DTResidualCalibration().

◆ rootBaseDir_

const std::string DTResidualCalibration::rootBaseDir_

private

Definition at line 59 of file DTResidualCalibration.h.

Referenced by bookHistos().

◆ rootFile_

TFile* DTResidualCalibration::rootFile_

private

Definition at line 62 of file DTResidualCalibration.h.

Referenced by analyze(), bookHistos(), DTResidualCalibration(), and endJob().

◆ segmbad

unsigned int DTResidualCalibration::segmbad

private

Definition at line 46 of file DTResidualCalibration.h.

Referenced by analyze(), DTResidualCalibration(), and ~DTResidualCalibration().

◆ segment4DToken_

const edm::EDGetTokenT<DTRecSegment4DCollection> DTResidualCalibration::segment4DToken_

private

Definition at line 58 of file DTResidualCalibration.h.

Referenced by analyze().

◆ segmok

unsigned int DTResidualCalibration::segmok

private

Definition at line 46 of file DTResidualCalibration.h.

Referenced by analyze(), DTResidualCalibration(), and ~DTResidualCalibration().

◆ select_

DTSegmentSelector* DTResidualCalibration::select_

private

Definition at line 56 of file DTResidualCalibration.h.

Referenced by analyze(), DTResidualCalibration(), and ~DTResidualCalibration().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTResidualCalibration()

◆ ~DTResidualCalibration()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginRun()

◆ bookHistos() [1/2]

◆ bookHistos() [2/2]

◆ endJob()

◆ endRun()

◆ fillHistos() [1/2]

◆ fillHistos() [2/2]

◆ segmentToWireDistance()

Member Data Documentation

◆ detailedAnalysis_

◆ dtGeom_

◆ dtGeomToken_

◆ histoMapPerLayerTH1F_

◆ histoMapPerLayerTH2F_

◆ histoMapTH1F_

◆ histoMapTH2F_

◆ histRange_

◆ nevent

◆ rootBaseDir_

◆ rootFile_

◆ segmbad

◆ segment4DToken_

◆ segmok

◆ select_