#include <DTSegmentsTask.h>

Inheritance diagram for DTSegmentsTask:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	book the histos More...

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

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

virtual	~DTSegmentsTask ()
	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

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

	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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
bool	checkNoisyChannels

bool	debug

edm::ParameterSet	parameters

std::vector< MonitorElement * >	phiHistos

edm::EDGetTokenT< DTRecSegment4DCollection >	theRecHits4DLabel_

std::vector< MonitorElement * >	thetaHistos

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

DQM Analysis of 4D DT segments

Author: G. Mila - INFN Torino

Definition at line 26 of file DTSegmentsTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 28 of file DTSegmentsTask.cc.

References debug, edm::ParameterSet::getUntrackedParameter(), metProducer_cfi::parameters, muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                           {
 
   debug = pset.getUntrackedParameter<bool>("debug",false); 
   parameters = pset;
   
   // should be init from pset, but it was commented out...
   // better false than undefined
   checkNoisyChannels = false;
 
   // the name of the 4D rec hits collection
   theRecHits4DLabel_ = consumes<DTRecSegment4DCollection>(parameters.getParameter<std::string>("recHits4DLabel"));
 }

DTSegmentsTask::~DTSegmentsTask ( )

virtual

Destructor.

Definition at line 42 of file DTSegmentsTask.cc.

42 {

43 }

Member Function Documentation

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

override

book the histos

Definition at line 110 of file DTSegmentsTask.cc.

References DTStatusFlag::cellStatus(), gather_cfg::cout, debug, SoftLeptonByDistance_cfi::distance, edm::EventSetup::get(), findQualityFiles::size, and DTRecSegment2D::specificRecHits().

                                                                               {
 
   // Get the map of noisy channels
   //  bool checkNoisyChannels = parameters.getUntrackedParameter<bool>("checkNoisyChannels",false);
   ESHandle<DTStatusFlag> statusMap;
   if(checkNoisyChannels) {
     setup.get<DTStatusFlagRcd>().get(statusMap);
   } 
 
   // Get the 4D segment collection from the event
   edm::Handle<DTRecSegment4DCollection> all4DSegments;
   event.getByToken(theRecHits4DLabel_, all4DSegments);
 
   // 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);
     int nsegm = distance(range.first, range.second);
     if(debug)
       cout << "   Chamber: " << *chamberId << " has " << nsegm
        << " 4D segments" << endl;
     
         // Loop over the rechits of this ChamerId
     for (DTRecSegment4DCollection::const_iterator segment4D = range.first; segment4D!=range.second; ++segment4D){
 
       //FOR NOISY CHANNELS////////////////////////////////
       bool segmNoisy = false;
       if((*segment4D).hasPhi()){
     const DTChamberRecSegment2D* phiSeg = (*segment4D).phiSegment();
     vector<DTRecHit1D> phiHits = phiSeg->specificRecHits();
     map<DTSuperLayerId,vector<DTRecHit1D> > hitsBySLMap; 
     for(vector<DTRecHit1D>::const_iterator hit = phiHits.begin(); hit != phiHits.end(); ++hit) {
       DTWireId wireId = (*hit).wireId();
       
       // Check for noisy channels to skip them
       if(checkNoisyChannels) {
         bool isNoisy = false;
         bool isFEMasked = false;
         bool isTDCMasked = false;
         bool isTrigMask = false;
         bool isDead = false;
         bool isNohv = false;
         statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
         if(isNoisy) {
           if(debug)
         cout << "Wire: " << wireId << " is noisy, skipping!" << endl;
           segmNoisy = true;
         }      
       }
     }
       }
       
       if((*segment4D).hasZed()) {
     const DTSLRecSegment2D* zSeg = (*segment4D).zSegment();  // zSeg lives in the SL RF
     // Check for noisy channels to skip them
     vector<DTRecHit1D> zHits = zSeg->specificRecHits();
     for(vector<DTRecHit1D>::const_iterator hit = zHits.begin();
         hit != zHits.end(); ++hit) {
       DTWireId wireId = (*hit).wireId();
       if(checkNoisyChannels) {
         bool isNoisy = false;
         bool isFEMasked = false;
         bool isTDCMasked = false;
         bool isTrigMask = false;
         bool isDead = false;
         bool isNohv = false;
         //cout<<"wire id "<<wireId<<endl;
         statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
         if(isNoisy) {
           if(debug)
         cout << "Wire: " << wireId << " is noisy, skipping!" << endl;
           segmNoisy = true;
         }      
       }
     }
       } 
       
       if (segmNoisy) {
     if(debug)
      cout<<"skipping the segment: it contains noisy cells"<<endl;
     continue;
       }
       //END FOR NOISY CHANNELS////////////////////////////////
       
       // Fill the histos
       int nHits=0;
       if((*segment4D).hasPhi()){
     nHits = (((*segment4D).phiSegment())->specificRecHits()).size();
     if(debug)
       cout<<"Phi segment with number of hits: "<<nHits<<endl;
     phiHistos[0]->Fill((*chamberId).wheel(), nHits);
     phiHistos[1]->Fill((*chamberId).sector(), nHits);
     phiHistos[2]->Fill((*chamberId).station(), nHits);
       }
       if((*segment4D).hasZed()) {
     nHits = (((*segment4D).zSegment())->specificRecHits()).size();
     if(debug)
       cout<<"Zed segment with number of hits: "<<nHits<<endl;
     thetaHistos[0]->Fill((*chamberId).wheel(), nHits);
     thetaHistos[1]->Fill((*chamberId).sector(), nHits);
     thetaHistos[2]->Fill((*chamberId).station(), nHits);
       }
       
     } //loop over segments
   } // loop over chambers
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 45 of file DTSegmentsTask.cc.

References DQMStore::IBooker::book2D(), DQMStore::IBooker::cd(), and DQMStore::IBooker::setCurrentFolder().

                                               {
   
   ibooker.cd();
   ibooker.setCurrentFolder("Muons/DTSegmentsMonitor");
 
   // histos for phi segments
   phiHistos.push_back(ibooker.book2D("phiSegments_numHitsVsWheel", "phiSegments_numHitsVsWheel", 5, -2.5, 2.5, 20, 0, 20));
   phiHistos[0]->setBinLabel(1,"W-2",1);
   phiHistos[0]->setBinLabel(2,"W-1",1);
   phiHistos[0]->setBinLabel(3,"W0",1);
   phiHistos[0]->setBinLabel(4,"W1",1);
   phiHistos[0]->setBinLabel(5,"W2",1);
   phiHistos.push_back(ibooker.book2D("phiSegments_numHitsVsSector", "phiSegments_numHitsVsSector", 14, 0.5, 14.5, 20, 0, 20));
   phiHistos[1]->setBinLabel(1,"Sec1",1);
   phiHistos[1]->setBinLabel(2,"Sec2",1);
   phiHistos[1]->setBinLabel(3,"Sec3",1);
   phiHistos[1]->setBinLabel(4,"Sec4",1);
   phiHistos[1]->setBinLabel(5,"Sec5",1);
   phiHistos[1]->setBinLabel(6,"Sec6",1);
   phiHistos[1]->setBinLabel(7,"Sec7",1);
   phiHistos[1]->setBinLabel(8,"Sec8",1);
   phiHistos[1]->setBinLabel(9,"Sec9",1);
   phiHistos[1]->setBinLabel(10,"Sec10",1);
   phiHistos[1]->setBinLabel(11,"Sec11",1);
   phiHistos[1]->setBinLabel(12,"Sec12",1);
   phiHistos[1]->setBinLabel(13,"Sec13",1);
   phiHistos[1]->setBinLabel(14,"Sec14",1);
   phiHistos.push_back(ibooker.book2D("phiSegments_numHitsVsStation", "phiSegments_numHitsVsStation", 4, 0.5, 4.5, 20, 0, 20));
   phiHistos[2]->setBinLabel(1,"St1",1);
   phiHistos[2]->setBinLabel(2,"St2",1);
   phiHistos[2]->setBinLabel(3,"St3",1);
   phiHistos[2]->setBinLabel(4,"St4",1);
 
   // histos for theta segments
   thetaHistos.push_back(ibooker.book2D("thetaSegments_numHitsVsWheel", "thetaSegments_numHitsVsWheel", 5, -2.5, 2.5, 20, 0, 20));
   thetaHistos[0]->setBinLabel(1,"W-2",1);
   thetaHistos[0]->setBinLabel(2,"W-1",1);
   thetaHistos[0]->setBinLabel(3,"W0",1);
   thetaHistos[0]->setBinLabel(4,"W1",1);
   thetaHistos[0]->setBinLabel(5,"W2",1);
   thetaHistos.push_back(ibooker.book2D("thetaSegments_numHitsVsSector", "thetaSegments_numHitsVsSector", 14, 0.5, 14.5, 20, 0, 20));
   thetaHistos[1]->setBinLabel(1,"Sec1",1);
   thetaHistos[1]->setBinLabel(2,"Sec2",1);
   thetaHistos[1]->setBinLabel(3,"Sec3",1);
   thetaHistos[1]->setBinLabel(4,"Sec4",1);
   thetaHistos[1]->setBinLabel(5,"Sec5",1);
   thetaHistos[1]->setBinLabel(6,"Sec6",1);
   thetaHistos[1]->setBinLabel(7,"Sec7",1);
   thetaHistos[1]->setBinLabel(8,"Sec8",1);
   thetaHistos[1]->setBinLabel(9,"Sec9",1);
   thetaHistos[1]->setBinLabel(10,"Sec10",1);
   thetaHistos[1]->setBinLabel(11,"Sec11",1);
   thetaHistos[1]->setBinLabel(12,"Sec12",1);
   thetaHistos[1]->setBinLabel(13,"Sec13",1);
   thetaHistos[1]->setBinLabel(14,"Sec14",1);
   thetaHistos.push_back(ibooker.book2D("thetaSegments_numHitsVsStation", "thetaSegments_numHitsVsStation", 4, 0.5, 4.5, 20, 0, 20));
   thetaHistos[2]->setBinLabel(1,"St1",1);
   thetaHistos[2]->setBinLabel(2,"St2",1);
   thetaHistos[2]->setBinLabel(3,"St3",1);
   thetaHistos[2]->setBinLabel(4,"St4",1);
 
 }

Member Data Documentation

bool DTSegmentsTask::checkNoisyChannels

private

Definition at line 46 of file DTSegmentsTask.h.

bool DTSegmentsTask::debug

private

Definition at line 45 of file DTSegmentsTask.h.

Referenced by util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), util.rrapi.RRApi::get(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), and pkg.AbstractPkg::write().

edm::ParameterSet DTSegmentsTask::parameters

private

Definition at line 47 of file DTSegmentsTask.h.

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

std::vector<MonitorElement*> DTSegmentsTask::phiHistos

private

Definition at line 50 of file DTSegmentsTask.h.

edm::EDGetTokenT<DTRecSegment4DCollection> DTSegmentsTask::theRecHits4DLabel_

private

Definition at line 54 of file DTSegmentsTask.h.

std::vector<MonitorElement*> DTSegmentsTask::thetaHistos

private

Definition at line 51 of file DTSegmentsTask.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation