DTResolutionAnalysisTask Class Reference

#include <DTResolutionAnalysisTask.h>

Inheritance diagram for DTResolutionAnalysisTask:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
	To reset the MEs. More...
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	BookHistograms. More...
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &)
	BeginRun. More...
	DTResolutionAnalysisTask (const edm::ParameterSet &pset)
	Constructor. More...
virtual	~DTResolutionAnalysisTask ()
	Destructor. More...
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	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
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) 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 (DQMStore::IBooker &ibooker, DTSuperLayerId slId)
void	fillHistos (DTSuperLayerId slId, float distExtr, float residual)

Private Attributes
edm::ESHandle< DTGeometry >	dtGeom
std::map< DTSuperLayerId, std::vector< MonitorElement * > >	histosPerSL
int	prescaleFactor
edm::EDGetTokenT < DTRecSegment4DCollection >	recHits4DToken_
int	resetCycle
u_int32_t	thePhiHitsCut
u_int32_t	theZHitsCut
std::string	topHistoFolder

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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

Definition at line 34 of file DTResolutionAnalysisTask.h.

Constructor & Destructor Documentation

DTResolutionAnalysisTask::DTResolutionAnalysisTask

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 35 of file DTResolutionAnalysisTask.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_25ns14e33_v1_cff::InputTag, and dtDQMClient_cfg::prescaleFactor.

                                                                           {

  edm::LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask") << "[DTResolutionAnalysisTask] Constructor called!" << endl;

  // the name of the 4D rec hits collection
  recHits4DToken_ = consumes<DTRecSegment4DCollection>(
      edm::InputTag(pset.getParameter<string>("recHits4DLabel")));

  prescaleFactor  = pset.getUntrackedParameter<int>("diagnosticPrescale", 1);
  resetCycle = pset.getUntrackedParameter<int>("ResetCycle", -1);
  // top folder for the histograms in DQMStore
  topHistoFolder  = pset.getUntrackedParameter<string>("topHistoFolder","DT/02-Segments");

  thePhiHitsCut = pset.getUntrackedParameter<u_int32_t>("phiHitsCut",8);
  theZHitsCut = pset.getUntrackedParameter<u_int32_t>("zHitsCut",4);

}

DTResolutionAnalysisTask::~DTResolutionAnalysisTask ( )

virtual

Destructor.

Definition at line 54 of file DTResolutionAnalysisTask.cc.

                                                   {

  edm::LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask") << "[DTResolutionAnalysisTask] Destructor called!" << endl;

}

Member Function Documentation

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

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 105 of file DTResolutionAnalysisTask.cc.

References filterCSVwithJSON::copy, funct::cos(), edm::HandleBase::isValid(), DTWireId::layerId(), DTRecSegment2D::specificRecHits(), DTLayer::specificTopology(), DTChamber::superLayer(), DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), GeomDet::toGlobal(), GeomDet::toLocal(), DTWireId::wire(), DTTopology::wirePosition(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                        {

  edm::LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask") << "[DTResolutionAnalysisTask] Analyze #Run: " << event.id().run()
     << " #Event: " << event.id().event() << endl;


  // Get the 4D segment collection from the event
  edm::Handle<DTRecSegment4DCollection> all4DSegments;
  event.getByToken(recHits4DToken_, all4DSegments);

  // check the validity of the collection
  if(!all4DSegments.isValid()) return;

  // Loop over all chambers containing a segment
  DTRecSegment4DCollection::id_iterator chamberId;
  for (chamberId = all4DSegments->id_begin();
       chamberId != all4DSegments->id_end();
       ++chamberId) {
    // Get the range for the corresponding ChamerId
    DTRecSegment4DCollection::range  range = all4DSegments->get(*chamberId);

    // Get the chamber
    const DTChamber* chamber = dtGeom->chamber(*chamberId);

    // Loop over the rechits of this ChamerId
    for (DTRecSegment4DCollection::const_iterator segment4D = range.first;
     segment4D!=range.second;
     ++segment4D) {

      // If Statio != 4 skip RecHits with dimension != 4
      // For the Station 4 consider 2D RecHits
      if((*chamberId).station() != 4 && (*segment4D).dimension() != 4) {
    edm::LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask") << "[DTResolutionAnalysisTask]***Warning: RecSegment dimension is not 4 but "
           << (*segment4D).dimension() << "!" << endl;
    continue;
      } else if((*chamberId).station() == 4 && (*segment4D).dimension() != 2) {
    edm::LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask") << "[DTResolutionAnalysisTask]***Warning: RecSegment dimension is not 2 but "
           << (*segment4D).dimension() << "!" << endl;
    continue;
      }


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

      if((*segment4D).hasPhi()) { // has phi component
    const DTChamberRecSegment2D* phiSeg = (*segment4D).phiSegment();
    vector<DTRecHit1D> phiRecHits = phiSeg->specificRecHits();

    if(phiRecHits.size() < thePhiHitsCut) {
      continue;
    }
    copy(phiRecHits.begin(), phiRecHits.end(), back_inserter(recHits1D_S3));
      } else {

      }

      if((*segment4D).hasZed()) {
    const DTSLRecSegment2D* zSeg = (*segment4D).zSegment();
    vector<DTRecHit1D> zRecHits = zSeg->specificRecHits();
    if(zRecHits.size() < theZHitsCut) {

      continue;
    }
    copy(zRecHits.begin(), zRecHits.end(), back_inserter(recHits1D_S3));
      }

      // 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 = chamber->superLayer(wireId.superlayerId())->layer(wireId.layerId());
    float wireX = layer->specificTopology().wirePosition(wireId.wire());

    // Distance of the 1D rechit from the wire
    float distRecHitToWire = fabs(wireX - (*recHit1D).localPosition().x());

    // Extrapolate the segment to the z of the wire

    // Get wire position in chamber RF
    LocalPoint wirePosInLay(wireX,(*recHit1D).localPosition().y(),(*recHit1D).localPosition().z());
    GlobalPoint wirePosGlob = layer->toGlobal(wirePosInLay);
    LocalPoint wirePosInChamber = chamber->toLocal(wirePosGlob);

    // Segment position at Wire z in chamber local frame
    LocalPoint segPosAtZWire = (*segment4D).localPosition()
      + (*segment4D).localDirection()*wirePosInChamber.z()/cos((*segment4D).localDirection().theta());

    // Compute the distance of the segment from the wire
    int sl = wireId.superlayer();

    double distSegmToWire = -1;
    if(sl == 1 || sl == 3) {
      // RPhi SL
      distSegmToWire = fabs(wirePosInChamber.x() - segPosAtZWire.x());
    } else if(sl == 2) {
      // RZ SL
      distSegmToWire = fabs(wirePosInChamber.y() - segPosAtZWire.y());
    }

    if(distSegmToWire > 2.1)
      edm::LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask") << "  Warning: dist segment-wire: " << distSegmToWire << endl;

    double residual = distRecHitToWire - distSegmToWire;
    // FIXME: Fill the histos
    fillHistos(wireId.superlayerId(), distSegmToWire, residual);

      }// End of loop over 1D RecHit inside 4D segment
    }// End of loop over the rechits of this ChamerId
  }
  // -----------------------------------------------------------------------------
}

void DTResolutionAnalysisTask::beginLuminosityBlock ( edm::LuminosityBlock const &  lumiSeg, edm::EventSetup const &  context  )

virtual

To reset the MEs.

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 88 of file DTResolutionAnalysisTask.cc.

References timingPdfMaker::histo, i, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), python.multivaluedict::map(), and findQualityFiles::size.

                                                       {

  edm::LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask") << "[DTResolutionTask]: Begin of LS transition"<<endl;

  if(resetCycle != -1 && lumiSeg.id().luminosityBlock() % resetCycle == 0) {
    for(map<DTSuperLayerId, vector<MonitorElement*> > ::const_iterator histo = histosPerSL.begin();
    histo != histosPerSL.end();
    histo++) {
      int size = (*histo).second.size();
      for(int i=0; i<size; i++){
    (*histo).second[i]->Reset();
      }
    }
  }
}

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

overridevirtual

BookHistograms.

Implements DQMEDAnalyzer.

Definition at line 67 of file DTResolutionAnalysisTask.cc.

References bookHistos(), chambers, and DTChamberId::station().

                                                                      {

  // Book the histograms
  vector<const DTChamber*> chambers = dtGeom->chambers();
  for(vector<const DTChamber*>::const_iterator chamber = chambers.begin();
      chamber != chambers.end(); ++chamber) {  // Loop over all chambers
    DTChamberId dtChId = (*chamber)->id();
    for(int sl = 1; sl <= 3; ++sl) { // Loop over SLs
      if(dtChId.station() == 4 && sl == 2) continue;
      const  DTSuperLayerId dtSLId(dtChId,sl);
      bookHistos(ibooker,dtSLId);
    }
  }

}

void DTResolutionAnalysisTask::bookHistos ( DQMStore::IBooker &  ibooker, DTSuperLayerId  slId  )

private

Definition at line 227 of file DTResolutionAnalysisTask.cc.

References DQMStore::IBooker::book1D(), mergeVDriftHistosByStation::histos, DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), and DTChamberId::wheel().

                                                                                        {

  edm::LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask") << "   Booking histos for SL: " << slId << endl;

  // Compose the chamber name
  stringstream wheel; wheel << slId.wheel();
  stringstream station; station << slId.station();
  stringstream sector; sector << slId.sector();
  stringstream superLayer; superLayer << slId.superlayer();


  string slHistoName =
    "_W" + wheel.str() +
    "_St" + station.str() +
    "_Sec" + sector.str() +
    "_SL" + superLayer.str();

  ibooker.setCurrentFolder(topHistoFolder + "/Wheel" + wheel.str() +
               "/Sector" + sector.str() +
               "/Station" + station.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));
  histosPerSL[slId] = histos;
}

void DTResolutionAnalysisTask::dqmBeginRun ( const edm::Run &  run, const edm::EventSetup &  setup  )

virtual

BeginRun.

Reimplemented from DQMEDAnalyzer.

Definition at line 60 of file DTResolutionAnalysisTask.cc.

References edm::EventSetup::get().

                                                                                  {

  // Get the DT Geometry
  setup.get<MuonGeometryRecord>().get(dtGeom);

}

void DTResolutionAnalysisTask::fillHistos ( DTSuperLayerId  slId, float  distExtr, float  residual  )

private

Definition at line 258 of file DTResolutionAnalysisTask.cc.

References mergeVDriftHistosByStation::histos.

                                      {
  vector<MonitorElement *> histos =  histosPerSL[slId];
  histos[0]->Fill(residual);
}

Member Data Documentation

edm::ESHandle<DTGeometry> DTResolutionAnalysisTask::dtGeom

private

Definition at line 59 of file DTResolutionAnalysisTask.h.

std::map<DTSuperLayerId, std::vector<MonitorElement*> > DTResolutionAnalysisTask::histosPerSL

private

Definition at line 77 of file DTResolutionAnalysisTask.h.

int DTResolutionAnalysisTask::prescaleFactor

private

Definition at line 61 of file DTResolutionAnalysisTask.h.

edm::EDGetTokenT<DTRecSegment4DCollection> DTResolutionAnalysisTask::recHits4DToken_

private

Definition at line 68 of file DTResolutionAnalysisTask.h.

int DTResolutionAnalysisTask::resetCycle

private

Definition at line 62 of file DTResolutionAnalysisTask.h.

u_int32_t DTResolutionAnalysisTask::thePhiHitsCut

private

Definition at line 64 of file DTResolutionAnalysisTask.h.

u_int32_t DTResolutionAnalysisTask::theZHitsCut

private

Definition at line 65 of file DTResolutionAnalysisTask.h.

std::string DTResolutionAnalysisTask::topHistoFolder

private

Definition at line 80 of file DTResolutionAnalysisTask.h.