#include <DTCalibValidation.h>

Inheritance diagram for DTCalibValidation:

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

void	beginJob ()
	BeginJob. More...

void	beginRun (const edm::Run &, const edm::EventSetup &)
	BeginRun. More...

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

void	endJob ()
	Endjob. More...

virtual	~DTCalibValidation ()
	Destructor. More...

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

template<typename type >
void	compute (const DTGeometry *dtGeom, const DTRecSegment4D &segment, const std::map< DTWireId, std::vector< type > > &recHitsPerWire, int step)

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

template<typename type >
const type *	findBestRecHit (const DTLayer *layer, DTWireId wireId, const std::vector< type > &recHits, const float simHitDist)

std::map< DTWireId, std::vector< DTRecHit1DPair > >	map1DRecHitsPerWire (const DTRecHitCollection *dt1DRecHitPairs)

std::map< DTWireId, std::vector< DTRecHit1D > >	map1DRecHitsPerWire (const DTRecSegment2DCollection *segment2Ds)

std::map< DTWireId, std::vector< DTRecHit1D > >	map1DRecHitsPerWire (const DTRecSegment4DCollection *segment4Ds)

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
bool	detailedAnalysis

edm::ESHandle< DTGeometry >	dtGeom

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

int	nevent

edm::ParameterSet	parameters

edm::EDGetTokenT < DTRecHitCollection >	recHits1DToken_

int	rightSegment

edm::EDGetTokenT < DTRecSegment2DCollection >	segment2DToken_

edm::EDGetTokenT < DTRecSegment4DCollection >	segment4DToken_

DQMStore *	theDbe

int	wrongSegment

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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

Detailed Description

Analysis on DT residuals to validate the kFactor

Author: G. Mila - INFN Torino

Definition at line 38 of file DTCalibValidation.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 34 of file DTCalibValidation.cc.

References cppFunctionSkipper::operator, and Parameters::parameters.

                                                              {
 
   //debug = pset.getUntrackedParameter<bool>("debug",false);
 
   // Get the DQM needed services
 
   // To remove into CMSSW versions before 20X
   theDbe = edm::Service<DQMStore>().operator->();
   // To add into CMSSW versions before 20X
   /*theDbe = edm::Service<DaqMonitorBEInterface>().operator->();
   theDbe->setVerbose(1);
   edm::Service<MonitorDaemon>().operator->();*/
 
   theDbe->setCurrentFolder("DT/DTCalibValidation");
 
   parameters = pset;
 }

DTCalibValidation::~DTCalibValidation ( )

virtual

Destructor.

Definition at line 53 of file DTCalibValidation.cc.

53 {

54 }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 124 of file DTCalibValidation.cc.

References bookConverter::compute(), LogTrace, nevent, and edm::Handle< T >::product().

                                                                                  {
 
   ++nevent;
   LogTrace("DTCalibValidation") << "[DTCalibValidation] Analyze #Run: " << event.id().run()
                                 << " #Event: " << nevent;
 
   // RecHit mapping at Step 1 -------------------------------
   map<DTWireId,vector<DTRecHit1DPair> > recHitsPerWire_1S;
 
   // RecHit mapping at Step 2 ------------------------------
   map<DTWireId,vector<DTRecHit1D> > recHitsPerWire_2S;
 
   if(detailedAnalysis){
      LogTrace("DTCalibValidation") << "  -- DTRecHit S1: begin analysis:";
      // Get the rechit collection from the event
      Handle<DTRecHitCollection> dtRecHits;
      event.getByToken(recHits1DToken_, dtRecHits);
      recHitsPerWire_1S = map1DRecHitsPerWire(dtRecHits.product());
 
      LogTrace("DTCalibValidation") << "  -- DTRecHit S2: begin analysis:";
      // Get the 2D rechits from the event
      Handle<DTRecSegment2DCollection> segment2Ds;
      event.getByToken(segment2DToken_, segment2Ds);
      recHitsPerWire_2S =  map1DRecHitsPerWire(segment2Ds.product());
   }
 
   // RecHit mapping at Step 3 ---------------------------------
   LogTrace("DTCalibValidation") << "  -- DTRecHit S3: begin analysis:";
   // Get the 4D rechits from the event
   Handle<DTRecSegment4DCollection> segment4Ds;
   event.getByToken(segment4DToken_, segment4Ds);
   map<DTWireId,vector<DTRecHit1D> > recHitsPerWire_3S =  map1DRecHitsPerWire(segment4Ds.product());
 
 
   // Loop over all 4D segments
   for(DTRecSegment4DCollection::const_iterator segment = segment4Ds->begin();
       segment != segment4Ds->end();
       ++segment) {
 
     if(detailedAnalysis){
        LogTrace("DTCalibValidation") << "Anlysis on recHit at step 1";
        compute(dtGeom.product(), (*segment), recHitsPerWire_1S, 1);
 
        LogTrace("DTCalibValidation") << "Anlysis on recHit at step 2";
        compute(dtGeom.product(), (*segment), recHitsPerWire_2S, 2);
     }
 
     LogTrace("DTCalibValidation") << "Anlysis on recHit at step 3";
     compute(dtGeom.product(), (*segment), recHitsPerWire_3S, 3);
   }
 
 }

void DTCalibValidation::beginJob ( void )

virtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 57 of file DTCalibValidation.cc.

References nevent, and Parameters::parameters.

                                 {
 
   // the name of the rechits collection at step 1
   recHits1DToken_ = consumes<DTRecHitCollection>(
       edm::InputTag(parameters.getUntrackedParameter<string>("recHits1DLabel")));
   // the name of the 2D segments
   segment2DToken_ = consumes<DTRecSegment2DCollection>(
       edm::InputTag(parameters.getUntrackedParameter<string>("segment2DLabel")));
   // the name of the 4D segments
   segment4DToken_ = consumes<DTRecSegment4DCollection>(
       edm::InputTag(parameters.getUntrackedParameter<string>("segment4DLabel")));
   // the counter of segments not used to compute residuals
   wrongSegment = 0;
   // the counter of segments used to compute residuals
   rightSegment = 0;
   // the analysis type
   detailedAnalysis = parameters.getUntrackedParameter<bool>("detailedAnalysis",false);
 
   nevent=0;
 
 }

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

virtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 79 of file DTCalibValidation.cc.

References bookHistos(), edm::EventSetup::get(), and relval_parameters_module::step.

                                                                         {
 
   // get the geometry
   setup.get<MuonGeometryRecord>().get(dtGeom);
 
   // 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) {
       DTSuperLayerId slId = (*sl_it)->id();
       if(detailedAnalysis){
          // Loop over the 3 steps
          for(int step = 1; step <= 3; ++step) bookHistos(slId,step);
       } else {
          // Only step 3
          bookHistos(slId,3);
       }
     }
   }
 
 }

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

private

Definition at line 481 of file DTCalibValidation.cc.

References mergeVDriftHistosByStation::histos, LogTrace, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, cmsPerfCommons::Step, relval_parameters_module::step, DTSuperLayerId::superlayer(), and DTChamberId::wheel().

                                                                 {
   LogTrace("DTCalibValidation") << "   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();
 
   theDbe->setCurrentFolder("DT/DTCalibValidation/Wheel" + wheel.str() +
                "/Station" + station.str() +
                "/Sector" + sector.str());
   // Create the monitor elements
   vector<MonitorElement *> histos;
   // Note hte order matters
   histos.push_back(theDbe->book1D("hResDist"+slHistoName,
                   "Residuals on the distance from wire (rec_hit - segm_extr) (cm)",
                   200, -0.4, 0.4));
   histos.push_back(theDbe->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));
   if(detailedAnalysis){
     histos.push_back(theDbe->book1D("hResPos"+slHistoName,
                     "Residuals on the position from wire (rec_hit - segm_extr) (cm)",
                     200, -0.4, 0.4));
     histos.push_back(theDbe->book2D("hResPosVsPos"+slHistoName,
                     "Residuals on the position (cm) from wire (rec_hit - segm_extr) vs distance  (cm)",
                     200, -2.5, 2.5, 200, -0.4, 0.4));
   }
 
   histosPerSL[make_pair(slId, step)] = histos;
 }

template<typename type >

void DTCalibValidation::compute	(	const DTGeometry *	dtGeom,
		const DTRecSegment4D &	segment,
		const std::map< DTWireId, std::vector< type > > &	recHitsPerWire,
		int	step
	)

private

Definition at line 351 of file DTCalibValidation.cc.

References DTGeometry::chamber(), filterCSVwithJSON::copy, funct::cos(), DTRecSegment4D::dimension(), DTGeometry::layer(), DTRecSegment4D::localDirection(), DTRecSegment4D::localPosition(), LogTrace, DTRecSegment4D::phiSegment(), DTRecSegment2D::specificRecHits(), DTLayer::specificTopology(), DTSuperLayerId::superlayer(), 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
   // 8 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() != 8) {
       LogTrace("DTCalibValidation") << "[DTCalibValidation] Phi segments has: " << phiRecHits.size()
                                     << " hits, skipping"; // FIXME: info output
       computeResidual = false;
     }
     copy(phiRecHits.begin(), phiRecHits.end(), back_inserter(recHits1D_S3));
   }
   if(!phiSeg){
     LogTrace("DTCalibValidation") << " [DTCalibValidation] 4D segment has not the 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("DTCalibValidation") << "[DTCalibValidation] Theta segments has: " << zRecHits.size()
                                   << " hits, skipping"; // FIXME: info output
     computeResidual = false;
       }
       copy(zRecHits.begin(), zRecHits.end(), back_inserter(recHits1D_S3));
     }
     if(!zSeg){
       LogTrace("DTCalibValidation") << " [DTCalibValidation] 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("DTCalibValidation") << "SegmDistance: " << SegmDistance;
       } else if(sl == 2) {
     // RZ SL
     SegmDistance =  fabs(segPosAtZWire.y() - wirePosInChamber.y());
     LogTrace("DTCalibValidation") << "SegmDistance: " << SegmDistance;
       }
       if(SegmDistance > 2.1)
     LogTrace("DTCalibValidation") << "  Warning: dist segment-wire: " << SegmDistance;
 
       // Look for RecHits in the same cell
       if(recHitsPerWire.find(wireId) == recHitsPerWire.end()) {
         LogTrace("DTCalibValidation") << "   No RecHit found at Step: " << step << " in cell: " << wireId;
       } else {
     vector<type> recHits = recHitsPerWire.at(wireId);
     LogTrace("DTCalibValidation") << "   " << recHits.size() << " RecHits, Step " << step << " in channel: " << wireId;
 
     // Get the layer
     const DTLayer* layer = dtGeom->layer(wireId);
     // Find the best RecHits
     const type* theBestRecHit = findBestRecHit(layer, wireId, recHits, SegmDistance);
     // Compute the distance of the recHit from the wire
     float recHitWireDist =  recHitDistFromWire(*theBestRecHit, layer);
     LogTrace("DTCalibValidation") << "recHitWireDist: " << recHitWireDist;
 
     // Compute the residuals
     float residualOnDistance = recHitWireDist - SegmDistance;
     LogTrace("DTCalibValidation") << "WireId: " << wireId << "  ResidualOnDistance: " << residualOnDistance;
     float residualOnPosition = -1;
     float recHitPos = -1;
     if(sl == 1 || sl == 3) {
       recHitPos = recHitPosition(*theBestRecHit, layer, chamber,  segPosAtZWire.x(), sl);
       residualOnPosition = recHitPos - segPosAtZWire.x();
     }
     else{
       recHitPos = recHitPosition(*theBestRecHit, layer, chamber,  segPosAtZWire.y(), sl);
       residualOnPosition = recHitPos - segPosAtZWire.y();
     }
     LogTrace("DTCalibValidation") << "WireId: " << wireId << "  ResidualOnPosition: " << residualOnPosition;
 
     // Fill the histos
     if(sl == 1 || sl == 3)
       fillHistos(wireId.superlayerId(), SegmDistance, residualOnDistance, (wirePosInChamber.x() - segPosAtZWire.x()), residualOnPosition, step);
     else
       fillHistos(wireId.superlayerId(), SegmDistance, residualOnDistance, (wirePosInChamber.y() - segPosAtZWire.y()), residualOnPosition, step);
 
 
       }
     }
   }
 
 }

void DTCalibValidation::endJob ( void )

virtual

Endjob.

Reimplemented from edm::EDAnalyzer.

Definition at line 106 of file DTCalibValidation.cc.

References dumpDBToFile_GT_ttrig_cfg::outputFileName, Parameters::parameters, and AlCaHLTBitMon_QueryRunRegistry::string.

                               {
 
  LogVerbatim("DTCalibValidation") << "Segments used to compute residuals: " << rightSegment;
  LogVerbatim("DTCalibValidation") << "Segments not used to compute residuals: " << wrongSegment;
 
  //theDbe->showDirStructure();
  bool outputMEsInRootFile = parameters.getParameter<bool>("OutputMEsInRootFile");
  std::string outputFileName = parameters.getParameter<std::string>("OutputFileName");
  if(outputMEsInRootFile){
    theDbe->showDirStructure();
    theDbe->save(outputFileName);
  }
 
  theDbe->rmdir("DT/DTCalibValidation");
 }

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

private

Definition at line 527 of file DTCalibValidation.cc.

References mergeVDriftHistosByStation::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);
   if(detailedAnalysis){
     histos[2]->Fill(residualOnPosition);
     histos[3]->Fill(position, residualOnPosition);
   }
 }

template<typename type >

const type * DTCalibValidation::findBestRecHit	(	const DTLayer *	layer,
		DTWireId	wireId,
		const std::vector< type > &	recHits,
		const float	simHitDist
	)

private

Definition at line 244 of file DTCalibValidation.cc.

                                                       {
   float res = 99999;
   const type* theBestRecHit = 0;
   // Loop over RecHits within the cell
   for(typename vector<type>::const_iterator recHit = recHits.begin();
       recHit != recHits.end();
       ++recHit) {
     float distTmp = recHitDistFromWire(*recHit, layer);
     if(fabs(distTmp-segmDist) < res) {
       res = fabs(distTmp-segmDist);
       theBestRecHit = &(*recHit);
     }
   } // End of loop over RecHits within the cell
 
   return theBestRecHit;
 }

map< DTWireId, vector< DTRecHit1DPair > > DTCalibValidation::map1DRecHitsPerWire ( const DTRecHitCollection * dt1DRecHitPairs )

private

Definition at line 180 of file DTCalibValidation.cc.

References run_regression::ret.

                                                                                 {
   map<DTWireId, vector<DTRecHit1DPair> > ret;
 
   for(DTRecHitCollection::const_iterator rechit = dt1DRecHitPairs->begin();
       rechit != dt1DRecHitPairs->end(); ++rechit) {
     ret[(*rechit).wireId()].push_back(*rechit);
   }
 
   return ret;
 }

map< DTWireId, vector< DTRecHit1D > > DTCalibValidation::map1DRecHitsPerWire ( const DTRecSegment2DCollection * segment2Ds )

private

Definition at line 194 of file DTCalibValidation.cc.

References run_regression::ret.

                                                                                  {
   map<DTWireId, vector<DTRecHit1D> > ret;
 
   // Loop over all 2D segments
   for(DTRecSegment2DCollection::const_iterator segment = segment2Ds->begin();
       segment != segment2Ds->end();
       ++segment) {
     vector<DTRecHit1D> component1DHits= (*segment).specificRecHits();
     // Loop over all component 1D hits
     for(vector<DTRecHit1D>::const_iterator hit = component1DHits.begin();
         hit != component1DHits.end(); ++hit) {
       ret[(*hit).wireId()].push_back(*hit);
     }
   }
   return ret;
 }

map< DTWireId, std::vector< DTRecHit1D > > DTCalibValidation::map1DRecHitsPerWire ( const DTRecSegment4DCollection * segment4Ds )

private

Definition at line 213 of file DTCalibValidation.cc.

References run_regression::ret, and DTRecHit1D::wireId().

                                                                                  {
   map<DTWireId, vector<DTRecHit1D> > ret;
   // Loop over all 4D segments
   for(DTRecSegment4DCollection::const_iterator segment = segment4Ds->begin();
       segment != segment4Ds->end();
       ++segment) {
     // Get component 2D segments
     vector<const TrackingRecHit*> segment2Ds = (*segment).recHits();
     // Loop over 2D segments:
     for(vector<const TrackingRecHit*>::const_iterator segment2D = segment2Ds.begin();
         segment2D != segment2Ds.end();
         ++segment2D) {
       // Get 1D component rechits
       vector<const TrackingRecHit*> hits = (*segment2D)->recHits();
       // Loop over them
       for(vector<const TrackingRecHit*>::const_iterator hit = hits.begin();
           hit != hits.end(); ++hit) {
         const DTRecHit1D* hit1D = dynamic_cast<const DTRecHit1D*>(*hit);
         ret[hit1D->wireId()].push_back(*hit1D);
       }
     }
   }
 
   return ret;
 }

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

private

Definition at line 267 of file DTCalibValidation.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 DTCalibValidation::recHitDistFromWire	(	const DTRecHit1D &	recHit,
		const DTLayer *	layer
	)

private

Definition at line 275 of file DTCalibValidation.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()));
 
   //to check the compatibility position / distance
   /* // Get the recHit and the wire position
   const DTChamber* chamber = (*layer).chamber();
   GlobalPoint recHitPosGlob = layer->toGlobal(recHit.localPosition());
   LocalPoint recHitPosInChamber = chamber->toLocal(recHitPosGlob);
   float wireX = layer->specificTopology().wirePosition(recHit.wireId().wire());
   LocalPoint wirePosInLay(wireX,recHit.localPosition().y(),0);
   GlobalPoint wirePosGlob = layer->toGlobal(wirePosInLay);
   LocalPoint wirePosInChamber = chamber->toLocal(wirePosGlob);
 
   float recHitDist = -1;
   if(recHit.wireId().layerId().superlayerId().superlayer() != 2)
     recHitDist = fabs(recHitPosInChamber.x()-wirePosInChamber.x());
   else
     recHitDist = fabs(recHitPosInChamber.y()-wirePosInChamber.y());
 
     return recHitDist; */
 
 }

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

private

Definition at line 303 of file DTCalibValidation.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 DTCalibValidation::recHitPosition	(	const DTRecHit1D &	recHit,
		const DTLayer *	layer,
		const DTChamber *	chamber,
		float	segmPos,
		int	sl
	)

private

Definition at line 331 of file DTCalibValidation.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;
 
 }

Member Data Documentation

bool DTCalibValidation::detailedAnalysis

private

Definition at line 75 of file DTCalibValidation.h.

edm::ESHandle<DTGeometry> DTCalibValidation::dtGeom

private

Definition at line 77 of file DTCalibValidation.h.

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

private

Definition at line 131 of file DTCalibValidation.h.

int DTCalibValidation::nevent

private

Definition at line 73 of file DTCalibValidation.h.

edm::ParameterSet DTCalibValidation::parameters

private

Definition at line 70 of file DTCalibValidation.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

edm::EDGetTokenT<DTRecHitCollection> DTCalibValidation::recHits1DToken_

private

Definition at line 80 of file DTCalibValidation.h.

int DTCalibValidation::rightSegment

private

Definition at line 72 of file DTCalibValidation.h.

edm::EDGetTokenT<DTRecSegment2DCollection> DTCalibValidation::segment2DToken_

private

Definition at line 82 of file DTCalibValidation.h.

edm::EDGetTokenT<DTRecSegment4DCollection> DTCalibValidation::segment4DToken_

private

Definition at line 84 of file DTCalibValidation.h.

DQMStore* DTCalibValidation::theDbe

private

Definition at line 64 of file DTCalibValidation.h.

int DTCalibValidation::wrongSegment

private

Definition at line 71 of file DTCalibValidation.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation