#include <DTCalibValidationFromMuons.h>

Inheritance diagram for DTCalibValidationFromMuons:

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

void	dqmBeginRun (const edm::Run &r, const edm::EventSetup &c) override
	BeginRun. More...

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

	~DTCalibValidationFromMuons () override
	Destructor. More...

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

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

virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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

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::vector< ModuleDescription const * > &modules, 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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

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<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)

Private Member Functions
void	bookHistos (DTSuperLayerId slId, int step)

void	compute (const DTGeometry *dtGeom, const DTRecSegment4D &segment)

void	fillHistos (DTSuperLayerId slId, float distance, float residualOnDistance, float position, float residualOnPosition, int step)

float	recHitDistFromWire (const DTRecHit1DPair &hitPair, const DTLayer *layer)

float	recHitDistFromWire (const DTRecHit1D &recHit, const DTLayer *layer)

float	recHitPosition (const DTRecHit1DPair &hitPair, const DTLayer layer, const DTChamber chamber, float segmPos, int sl)

float	recHitPosition (const DTRecHit1D &recHit, const DTLayer layer, const DTChamber chamber, float segmPos, int sl)

Private Attributes
edm::ESHandle< DTGeometry >	dtGeom

std::map< std::pair< DTSuperLayerId, int >, std::vector< MonitorElement * > >	histosPerSL

edm::EDGetTokenT< reco::MuonCollection >	muonToken_

int	nevent

edm::ParameterSet	parameters

int	rightSegment

edm::EDGetTokenT< DTRecSegment4DCollection >	segment4DToken_

int	wrongSegment

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Analysis on DT residuals to validate the kFactor

Author: G. Mila - INFN Torino

Definition at line 44 of file DTCalibValidationFromMuons.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 32 of file DTCalibValidationFromMuons.cc.

References HLT_2018_cff::InputTag, nevent, and muonDTDigis_cfi::pset.

                                                                                {
   parameters = pset;
 
   // the name of the 4D segments
   segment4DToken_ =
       consumes<DTRecSegment4DCollection>(edm::InputTag(parameters.getUntrackedParameter<string>("segment4DLabel")));
   // muon collection for matching 4D segments to muons
   muonToken_ = consumes<reco::MuonCollection>(edm::InputTag(parameters.getUntrackedParameter<string>("muonLabel")));
   // the counter of segments not used to compute residuals
   wrongSegment = 0;
   // the counter of segments used to compute residuals
   rightSegment = 0;
 
   nevent = 0;
 }

DTCalibValidationFromMuons::~DTCalibValidationFromMuons ( )

override

Destructor.

Definition at line 48 of file DTCalibValidationFromMuons.cc.

                                                         {
   // FR the following was previously in the endJob
 
   LogVerbatim("DTCalibValidationFromMuons") << "Segments used to compute residuals: " << rightSegment;
   LogVerbatim("DTCalibValidationFromMuons") << "Segments not used to compute residuals: " << wrongSegment;
 }

Member Function Documentation

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

override

Definition at line 60 of file DTCalibValidationFromMuons.cc.

References bookConverter::compute(), flavorHistoryFilter_cfi::dr, MuonSubdetId::DT, PVValHelper::dx, PVValHelper::dy, LogTrace, DetId::Muon, PDWG_BPHSkim_cff::muons, nevent, edm::Handle< T >::product(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                                                           {
   ++nevent;
   LogTrace("DTCalibValidationFromMuons") << "[DTCalibValidationFromMuons] Analyze #Run: " << event.id().run()
                                          << " #Event: " << nevent;
 
   // RecHit mapping at Step 3 ---------------------------------
   LogTrace("DTCalibValidationFromMuons") << "  -- DTRecHit S3: begin analysis:";
   // Get the 4D rechits from the event
   Handle<reco::MuonCollection> muonH;
   event.getByToken(muonToken_, muonH);
   const vector<reco::Muon> *muons = muonH.product();
 
   // Get the 4D rechits from the event
   Handle<DTRecSegment4DCollection> segment4Ds;
   event.getByToken(segment4DToken_, segment4Ds);
 
   vector<const DTRecSegment4D *> selectedSegment4Ds;
 
   for (auto &imuon : *muons) {
     for (const auto &ch : imuon.matches()) {
       DetId chId(ch.id.rawId());
       if (chId.det() != DetId::Muon)
         continue;
       if (chId.subdetId() != MuonSubdetId::DT)
         continue;
       if (imuon.pt() < 15)
         continue;
       if (!imuon.isGlobalMuon())
         continue;
 
       int nsegs = ch.segmentMatches.size();
       if (!nsegs)
         continue;
 
       // get the DT segments that were used to construct the muon
       DTChamberId matchId = ch.id();
       DTRecSegment4DCollection::range segs = segment4Ds->get(matchId);
       for (DTRecSegment4DCollection::const_iterator segment = segs.first; segment != segs.second; ++segment) {
         LocalPoint posHit = segment->localPosition();
         float dx = (posHit.x() ? posHit.x() - ch.x : 0);
         float dy = (posHit.y() ? posHit.y() - ch.y : 0);
         float dr = sqrt(dx * dx + dy * dy);
         if (dr < 5)
           selectedSegment4Ds.push_back(&(*segment));
       }
     }
   }
 
   // Loop over all 4D segments
   for (auto segment : selectedSegment4Ds) {
     LogTrace("DTCalibValidationFromMuons") << "Anlysis on recHit at step 3";
     compute(dtGeom.product(), *segment);
   }
 }

void DTCalibValidationFromMuons::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 288 of file DTCalibValidationFromMuons.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), combine::histos, LogTrace, DTChamberId::sector(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                              {
   // FR substitute the DQMStore instance by ibooker
   ibooker.setCurrentFolder("DT/DTCalibValidationFromMuons");
 
   DTSuperLayerId slId;
 
   // Loop over all the chambers
   vector<const DTChamber *>::const_iterator ch_it = dtGeom->chambers().begin();
   vector<const DTChamber *>::const_iterator ch_end = dtGeom->chambers().end();
   for (; ch_it != ch_end; ++ch_it) {
     vector<const DTSuperLayer *>::const_iterator sl_it = (*ch_it)->superLayers().begin();
     vector<const DTSuperLayer *>::const_iterator sl_end = (*ch_it)->superLayers().end();
     // Loop over the SLs
     for (; sl_it != sl_end; ++sl_it) {
       slId = (*sl_it)->id();
 
       // TODO! fix this is a leftover
       int firstStep = 3;
       // Loop over the 3 steps
       for (int step = firstStep; step <= 3; ++step) {
         LogTrace("DTCalibValidationFromMuons") << "   Booking histos for SL: " << slId;
 
         // Compose the chamber name
         stringstream wheel;
         wheel << slId.wheel();
         stringstream station;
         station << slId.station();
         stringstream sector;
         sector << slId.sector();
         stringstream superLayer;
         superLayer << slId.superlayer();
         // Define the step
         stringstream Step;
         Step << step;
 
         string slHistoName = "_STEP" + Step.str() + "_W" + wheel.str() + "_St" + station.str() + "_Sec" + sector.str() +
                              "_SL" + superLayer.str();
 
         ibooker.setCurrentFolder("DT/DTCalibValidationFromMuons/Wheel" + wheel.str() + "/Station" + station.str() +
                                  "/Sector" + sector.str());
         // Create the monitor elements
         vector<MonitorElement *> histos;
         // Note the order matters
         histos.push_back(ibooker.book1D(
             "hResDist" + slHistoName, "Residuals on the distance from wire (rec_hit - segm_extr) (cm)", 200, -0.4, 0.4));
         histos.push_back(ibooker.book2D("hResDistVsDist" + slHistoName,
                                         "Residuals on the distance (cm) from wire (rec_hit "
                                         "- segm_extr) vs distance  (cm)",
                                         100,
                                         0,
                                         2.5,
                                         200,
                                         -0.4,
                                         0.4));
 
         histosPerSL[make_pair(slId, step)] = histos;
       }
     }
   }
 }

void DTCalibValidationFromMuons::bookHistos	(	DTSuperLayerId	slId,
		int	step
	)

private

void DTCalibValidationFromMuons::compute	(	const DTGeometry *	dtGeom,
		const DTRecSegment4D &	segment
	)

private

Definition at line 167 of file DTCalibValidationFromMuons.cc.

References relativeConstraints::chamber, DTGeometry::chamber(), filterCSVwithJSON::copy, funct::cos(), DTRecSegment4D::dimension(), DTGeometry::layer(), DTRecSegment4D::localDirection(), DTRecSegment4D::localPosition(), LogTrace, DTRecSegment4D::phiSegment(), DTRecSegment2D::specificRecHits(), DTLayer::specificTopology(), DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), PV3DBase< T, PVType, FrameType >::theta(), GeomDet::toGlobal(), GeomDet::toLocal(), DTWireId::wire(), DTTopology::wirePosition(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and DTRecSegment4D::zSegment().

                                                                                                 {
   bool computeResidual = true;
 
   // Get all 1D RecHits at step 3 within the 4D segment
   vector<DTRecHit1D> recHits1D_S3;
 
   // Get 1D RecHits at Step 3 and select only events with
   // >=7 hits in phi and 4 hits in theta (if any)
   const DTChamberRecSegment2D *phiSeg = segment.phiSegment();
   if (phiSeg) {
     vector<DTRecHit1D> phiRecHits = phiSeg->specificRecHits();
     if (phiRecHits.size() < 7) {
       LogTrace("DTCalibValidationFromMuons") << "[DTCalibValidationFromMuons] Phi segments has: " << phiRecHits.size()
                                              << " hits, skipping";  // FIXME: info output
       computeResidual = false;
     }
     copy(phiRecHits.begin(), phiRecHits.end(), back_inserter(recHits1D_S3));
   }
   if (!phiSeg) {
     LogTrace("DTCalibValidationFromMuons") << " [DTCalibValidationFromMuons] 4D segment has no phi segment! ";
     computeResidual = false;
   }
 
   if (segment.dimension() == 4) {
     const DTSLRecSegment2D *zSeg = segment.zSegment();
     if (zSeg) {
       vector<DTRecHit1D> zRecHits = zSeg->specificRecHits();
       if (zRecHits.size() != 4) {
         LogTrace("DTCalibValidationFromMuons") << "[DTCalibValidationFromMuons] Theta segments has: " << zRecHits.size()
                                                << " hits, skipping";  // FIXME: info output
         computeResidual = false;
       }
       copy(zRecHits.begin(), zRecHits.end(), back_inserter(recHits1D_S3));
     }
     if (!zSeg) {
       LogTrace("DTCalibValidationFromMuons") << " [DTCalibValidationFromMuons] 4D segment has not the z segment! ";
       computeResidual = false;
     }
   }
 
   if (!computeResidual)
     ++wrongSegment;
 
   if (computeResidual) {
     ++rightSegment;
 
     // Loop over 1D RecHit inside 4D segment
     for (vector<DTRecHit1D>::const_iterator recHit1D = recHits1D_S3.begin(); recHit1D != recHits1D_S3.end();
          ++recHit1D) {
       const DTWireId wireId = (*recHit1D).wireId();
 
       // Get the layer and the wire position
       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((*recHit1D).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) {
         // RPhi SL
         SegmDistance = fabs(wirePosInChamber.x() - segPosAtZWire.x());
         LogTrace("DTCalibValidationFromMuons") << "SegmDistance: " << SegmDistance;
       } else if (sl == 2) {
         // RZ SL
         SegmDistance = fabs(segPosAtZWire.y() - wirePosInChamber.y());
         LogTrace("DTCalibValidationFromMuons") << "SegmDistance: " << SegmDistance;
       }
 
       if (SegmDistance > 2.1)
         LogTrace("DTCalibValidationFromMuons") << "  Warning: dist segment-wire: " << SegmDistance;
 
       // Compute the distance of the recHit from the wire
       float recHitWireDist = recHitDistFromWire(*recHit1D, layer);
       LogTrace("DTCalibValidationFromMuons") << "recHitWireDist: " << recHitWireDist;
 
       // Compute the residuals
       float residualOnDistance = recHitWireDist - SegmDistance;
       LogTrace("DTCalibValidationFromMuons") << "WireId: " << wireId << "  ResidualOnDistance: " << residualOnDistance;
       float residualOnPosition = -1;
       float recHitPos = -1;
       if (sl == 1 || sl == 3) {
         recHitPos = recHitPosition(*recHit1D, layer, chamber, segPosAtZWire.x(), sl);
         residualOnPosition = recHitPos - segPosAtZWire.x();
       } else {
         recHitPos = recHitPosition(*recHit1D, layer, chamber, segPosAtZWire.y(), sl);
         residualOnPosition = recHitPos - segPosAtZWire.y();
       }
       LogTrace("DTCalibValidationFromMuons") << "WireId: " << wireId << "  ResidualOnPosition: " << residualOnPosition;
 
       // Fill the histos
       if (sl == 1 || sl == 3)
         fillHistos(wireId.superlayerId(),
                    SegmDistance,
                    residualOnDistance,
                    (wirePosInChamber.x() - segPosAtZWire.x()),
                    residualOnPosition,
                    3);
       else
         fillHistos(wireId.superlayerId(),
                    SegmDistance,
                    residualOnDistance,
                    (wirePosInChamber.y() - segPosAtZWire.y()),
                    residualOnPosition,
                    3);
     }
   }
 }

void DTCalibValidationFromMuons::dqmBeginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

override

BeginRun.

Definition at line 55 of file DTCalibValidationFromMuons.cc.

References edm::EventSetup::get().

                                                                                           {
   // get the geometry
   setup.get<MuonGeometryRecord>().get(dtGeom);
 }

void DTCalibValidationFromMuons::fillHistos	(	DTSuperLayerId	slId,
		float	distance,
		float	residualOnDistance,
		float	position,
		float	residualOnPosition,
		int	step
	)

private

Definition at line 352 of file DTCalibValidationFromMuons.cc.

References combine::histos.

                                                                                                                        {
   // FIXME: optimization of the number of searches
   vector<MonitorElement *> histos = histosPerSL[make_pair(slId, step)];
   histos[0]->Fill(residualOnDistance);
   histos[1]->Fill(distance, residualOnDistance);
 }

float DTCalibValidationFromMuons::recHitDistFromWire	(	const DTRecHit1DPair &	hitPair,
		const DTLayer *	layer
	)

private

Definition at line 116 of file DTCalibValidationFromMuons.cc.

References DTEnums::Left, DTRecHit1DPair::localPosition(), DTEnums::Right, and PV3DBase< T, PVType, FrameType >::x().

                                                                                                         {
   return fabs(hitPair.localPosition(DTEnums::Left).x() - hitPair.localPosition(DTEnums::Right).x()) / 2.;
 }

float DTCalibValidationFromMuons::recHitDistFromWire	(	const DTRecHit1D &	recHit,
		const DTLayer *	layer
	)

private

Definition at line 121 of file DTCalibValidationFromMuons.cc.

References DTRecHit1D::localPosition(), DTLayer::specificTopology(), DTWireId::wire(), DTRecHit1D::wireId(), DTTopology::wirePosition(), and PV3DBase< T, PVType, FrameType >::x().

                                                                                                    {
   return fabs(recHit.localPosition().x() - layer->specificTopology().wirePosition(recHit.wireId().wire()));
 }

float DTCalibValidationFromMuons::recHitPosition	(	const DTRecHit1DPair &	hitPair,
		const DTLayer *	layer,
		const DTChamber *	chamber,
		float	segmPos,
		int	sl
	)

private

Definition at line 126 of file DTCalibValidationFromMuons.cc.

References DTEnums::Left, DTRecHit1DPair::localPosition(), DTEnums::Right, GeomDet::toGlobal(), GeomDet::toLocal(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                                                                              {
   // Get the layer and the wire position
   GlobalPoint hitPosGlob_right = layer->toGlobal(hitPair.localPosition(DTEnums::Right));
   LocalPoint hitPosInChamber_right = chamber->toLocal(hitPosGlob_right);
   GlobalPoint hitPosGlob_left = layer->toGlobal(hitPair.localPosition(DTEnums::Left));
   LocalPoint hitPosInChamber_left = chamber->toLocal(hitPosGlob_left);
 
   float recHitPos = -1;
   if (sl != 2) {
     if (fabs(hitPosInChamber_left.x() - segmentPos) < fabs(hitPosInChamber_right.x() - segmentPos))
       recHitPos = hitPosInChamber_left.x();
     else
       recHitPos = hitPosInChamber_right.x();
   } else {
     if (fabs(hitPosInChamber_left.y() - segmentPos) < fabs(hitPosInChamber_right.y() - segmentPos))
       recHitPos = hitPosInChamber_left.y();
     else
       recHitPos = hitPosInChamber_right.y();
   }
 
   return recHitPos;
 }

float DTCalibValidationFromMuons::recHitPosition	(	const DTRecHit1D &	recHit,
		const DTLayer *	layer,
		const DTChamber *	chamber,
		float	segmPos,
		int	sl
	)

private

Definition at line 151 of file DTCalibValidationFromMuons.cc.

References DTRecHit1D::localPosition(), GeomDet::toGlobal(), GeomDet::toLocal(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                                                                         {
   // Get the layer and the wire position
   GlobalPoint recHitPosGlob = layer->toGlobal(recHit.localPosition());
   LocalPoint recHitPosInChamber = chamber->toLocal(recHitPosGlob);
 
   float recHitPos = -1;
   if (sl != 2)
     recHitPos = recHitPosInChamber.x();
   else
     recHitPos = recHitPosInChamber.y();
 
   return recHitPos;
 }