#include <DTNoiseTask.h>

Inheritance diagram for DTNoiseTask:

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

	DTNoiseTask (const edm::ParameterSet &ps)
	Constructor. More...

	~DTNoiseTask () override
	Destructor. More...

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

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

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

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

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

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

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

	~DQMEDAnalyzer () override=default

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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

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)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
	Analyze. More...

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

void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &setup) override

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)

Private Member Functions
void	bookHistos (DQMStore::IBooker &, DTChamberId chId)

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

Private Attributes
bool	doSegmentVeto

bool	doTimeBoxHistos

edm::EDGetTokenT< DTDigiCollection >	dtDigiToken_

edm::ESHandle< DTGeometry >	dtGeom

int	evtNumber

std::map< DTChamberId, int >	mapEvt

MonitorElement *	nEventMonitor

std::map< DTChamberId, MonitorElement * >	noiseHistos

edm::EDGetTokenT< DTRecSegment4DCollection >	recHits4DToken_

double	safeMargin

std::map< DTSuperLayerId, MonitorElement * >	tbHistos

edm::ESHandle< DTTtrig >	tTrigMap

std::map< DTChamberId, double >	tTrigStMap

Additional Inherited Members
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

No description available.

Authors: G. Mila , G. Cerminara - INFN Torino

Definition at line 41 of file DTNoiseTask.h.

Constructor & Destructor Documentation

DTNoiseTask::DTNoiseTask ( const edm::ParameterSet & ps )

Constructor.

Definition at line 41 of file DTNoiseTask.cc.

References doSegmentVeto, doTimeBoxHistos, dtDigiToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), recHits4DToken_, and safeMargin.

                                                : evtNumber(0) {
 
    LogVerbatim("DTNoiseTask") << "[DTNoiseTask]: Constructor"<<endl;
 
   //switch for timeBox booking
   doTimeBoxHistos = ps.getUntrackedParameter<bool>("doTbHistos", false);
 
   // The label to retrieve the digis
   dtDigiToken_ = consumes<DTDigiCollection>(
       ps.getParameter<InputTag>("dtDigiLabel"));
 
   // the name of the 4D rec hits collection
   recHits4DToken_ = consumes<DTRecSegment4DCollection>(
       edm::InputTag(ps.getParameter<string>("recHits4DLabel")));
 
   // switch for segment veto
   doSegmentVeto = ps.getUntrackedParameter<bool>("doSegmentVeto", false);
 
   // safe margin (ns) between ttrig and beginning of counting area
   safeMargin = ps.getUntrackedParameter<double>("safeMargin", 200.);
 
 }

DTNoiseTask::~DTNoiseTask ( )

override

Destructor.

Definition at line 67 of file DTNoiseTask.cc.

67 {}

Member Function Documentation

void DTNoiseTask::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotected

Analyze.

Definition at line 70 of file DTNoiseTask.cc.

References doSegmentVeto, doTimeBoxHistos, dtDigiToken_, evtNumber, MonitorElement::Fill(), edm::Event::getByToken(), mapEvt, nEventMonitor, noiseHistos, PostProcessor_cff::normalization, recHits4DToken_, safeMargin, edm::second(), tbHistos, tTrigStMap, and cuy::yBin.

                                                                    {
 
   evtNumber++;
   nEventMonitor->Fill(evtNumber);
 
   if(evtNumber%1000==0)
      LogVerbatim("DTNoiseTask") <<"[DTNoiseTask]: Analyzing evt number :"<<evtNumber<<endl;
 
   // map of the chambers with at least 1 segment
   std::map<DTChamberId, int> segmentsChId;
 
   // Get the 4D segment collection from the event
   edm::Handle<DTRecSegment4DCollection> all4DSegments;
   if(doSegmentVeto) {
     e.getByToken(recHits4DToken_, all4DSegments);
 
     // Loop over all chambers containing a segment and look for the number of segments
     DTRecSegment4DCollection::id_iterator chamberId;
     for (chamberId = all4DSegments->id_begin();
      chamberId != all4DSegments->id_end();
      ++chamberId){
       segmentsChId[*chamberId]=1;
     }
   }
 
   // Get the digis from the event
   edm::Handle<DTDigiCollection> dtdigis;
   e.getByToken(dtDigiToken_, dtdigis);
 
   // LOOP OVER ALL THE DIGIS OF THE EVENT
   DTDigiCollection::DigiRangeIterator dtLayerId_It;
   for (dtLayerId_It=dtdigis->begin(); dtLayerId_It!=dtdigis->end(); ++dtLayerId_It){
     for (DTDigiCollection::const_iterator digiIt = ((*dtLayerId_It).second).first;
      digiIt!=((*dtLayerId_It).second).second; ++digiIt){
 
       //Check the TDC trigger width
       int tdcTime = (*digiIt).countsTDC();
       double upperLimit = tTrigStMap[(*dtLayerId_It).first.superlayerId().chamberId()]-safeMargin;
       if(doTimeBoxHistos)
     tbHistos[(*dtLayerId_It).first.superlayerId()]->Fill(tdcTime);
       if(tdcTime>upperLimit)
     continue;
 
       //Check the chamber has no 4D segments (optional)
       if(doSegmentVeto &&
      segmentsChId.find((*dtLayerId_It).first.superlayerId().chamberId())!=segmentsChId.end())
     continue;
 
       // fill the occupancy histo
       // FIXME: needs to be optimized: no need to rescale the histo for each digi
       TH2F* noise_root = noiseHistos[(*dtLayerId_It).first.superlayerId().chamberId()]->getTH2F();
       double normalization=0;
       if(mapEvt.find((*dtLayerId_It).first.superlayerId().chamberId())!=mapEvt.end()) {
     LogVerbatim("DTNoiseTask")  << " Last fill: # of events: "
                     << mapEvt[(*dtLayerId_It).first.superlayerId().chamberId()]
                     << endl;
     normalization =  1e-9*upperLimit*mapEvt[(*dtLayerId_It).first.superlayerId().chamberId()];
     // revert back to # of entries
     noise_root->Scale(normalization);
       }
       int yBin=(*dtLayerId_It).first.layer()+(4*((*dtLayerId_It).first.superlayerId().superlayer()-1));
       noise_root->Fill((*digiIt).wire(),yBin);
       // normalize the occupancy histo
       mapEvt[(*dtLayerId_It).first.superlayerId().chamberId()] = evtNumber;
       LogVerbatim("DTNoiseTask")  << (*dtLayerId_It).first << " wire: " << (*digiIt).wire()
                   << " # counts: " << noise_root->GetBinContent((*digiIt).wire(),yBin)
                   << " Time interval: " << upperLimit
                   << " # of events: " << evtNumber << endl;;
       normalization = double( 1e-9*upperLimit*mapEvt[(*dtLayerId_It).first.superlayerId().chamberId()]);
       // update the rate
       noise_root->Scale(1./normalization);
       LogVerbatim("DTNoiseTask")  << "    noise rate: "
                   << noise_root->GetBinContent((*digiIt).wire(),yBin) << endl;
     }
   }
 }

void DTNoiseTask::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	run,
		edm::EventSetup const &	setup
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 246 of file DTNoiseTask.cc.

References DQMStore::IBooker::bookFloat(), bookHistos(), DTGeometry::chambers(), doTimeBoxHistos, dtGeom, DTTtrig::get(), DTAnalyzerDetailed_cfi::kFactor, nEventMonitor, DTTimeUnits::ns, DQMStore::IBooker::setCurrentFolder(), tTrigMap, and tTrigStMap.

                                                                          {
 
   ibooker.setCurrentFolder("DT/EventInfo/Counters");
   nEventMonitor = ibooker.bookFloat("nProcessedEventsNoise");
 
 
 
   // 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) {
     DTChamberId chId = (*ch_it)->id();
     // histo booking
     bookHistos(ibooker,chId);
     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(doTimeBoxHistos)
     bookHistos(ibooker,slId);
       float tTrig, tTrigRMS, kFactor;
       tTrigMap->get(slId, tTrig, tTrigRMS,kFactor,DTTimeUnits::ns);
       // tTrig mapping per station
       // check that the ttrig is the lowest of the 3 SLs
       if(tTrigStMap.find(chId)==tTrigStMap.end() ||
      (tTrigStMap.find(chId)!=tTrigStMap.end() && tTrig < tTrigStMap[chId]))
     tTrigStMap[chId] = tTrig;
     }
   }
 
 
 }

void DTNoiseTask::bookHistos	(	DQMStore::IBooker &	ibooker,
		DTChamberId	chId
	)

private

Definition at line 147 of file DTNoiseTask.cc.

References DQMStore::IBooker::book2D(), DTGeometry::chamber(), dtGeom, LayerTriplets::layers(), noiseHistos, DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, DTChamber::superLayers(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

                                                                        {
 
   // set the folder
   stringstream wheel; wheel << chId.wheel();
   stringstream station; station << chId.station();
   stringstream sector; sector << chId.sector();
 
   ibooker.setCurrentFolder("DT/05-Noise/Wheel" + wheel.str() +
 //          "/Station" + station.str() +
             "/Sector" + sector.str());
 
   // Build the histo name
   string histoName = string("NoiseRate")
     + "_W" + wheel.str()
     + "_St" + station.str()
     + "_Sec" + sector.str() ;
 
    LogVerbatim("DTNoiseTask") << "[DTNoiseTask]: booking chamber histo:"<<endl;
    LogVerbatim("DTNoiseTask") << "              folder "<< "DT/05-Noise/Wheel" + wheel.str() +
 //     "/Station" + station.str() +
     "/Sector" + sector.str() + "/"<<endl;
    LogVerbatim("DTNoiseTask") << "              histoName "<<histoName<<endl;
 
   // Get the chamber from the geometry
   int nWires_max = 0;
   const DTChamber* dtchamber = dtGeom->chamber(chId);
   const vector<const DTSuperLayer*>& superlayers = dtchamber->superLayers();
 
   // Loop over layers and find the max # of wires
   for(vector<const DTSuperLayer*>::const_iterator sl = superlayers.begin();
       sl != superlayers.end(); ++sl) { // loop over SLs
     vector<const DTLayer*> layers = (*sl)->layers();
     for(vector<const DTLayer*>::const_iterator lay = layers.begin();
     lay != layers.end(); ++lay) { // loop over layers
       int nWires = (*lay)->specificTopology().channels();
       if(nWires > nWires_max) nWires_max = nWires;
     }
   }
 
   noiseHistos[chId] = ibooker.book2D(histoName,"Noise rate (Hz) per channel", nWires_max,1, nWires_max+1,12,1,13);
   noiseHistos[chId]->setAxisTitle("wire number",1);
   noiseHistos[chId]->setBinLabel(1,"SL1-L1",2);
   noiseHistos[chId]->setBinLabel(2,"SL1-L2",2);
   noiseHistos[chId]->setBinLabel(3,"SL1-L3",2);
   noiseHistos[chId]->setBinLabel(4,"SL1-L4",2);
   noiseHistos[chId]->setBinLabel(5,"SL2-L1",2);
   noiseHistos[chId]->setBinLabel(6,"SL2-L2",2);
   noiseHistos[chId]->setBinLabel(7,"SL2-L3",2);
   noiseHistos[chId]->setBinLabel(8,"SL2-L4",2);
   noiseHistos[chId]->setBinLabel(9,"SL3-L1",2);
   noiseHistos[chId]->setBinLabel(10,"SL3-L2",2);
   noiseHistos[chId]->setBinLabel(11,"SL3-L3",2);
   noiseHistos[chId]->setBinLabel(12,"SL3-L4",2);
 
 }

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

private

Definition at line 204 of file DTNoiseTask.cc.

References DQMStore::IBooker::book1D(), DTSuperLayerId::chamberId(), DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, DTSuperLayerId::superlayer(), tbHistos, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                           {
 
   // set the folder
   stringstream wheel; wheel << slId.chamberId().wheel();
   stringstream station; station << slId.chamberId().station();
   stringstream sector; sector << slId.chamberId().sector();
   stringstream superlayer; superlayer << slId.superlayer();
 
   ibooker.setCurrentFolder("DT/05-Noise/Wheel" + wheel.str() +
             "/Station" + station.str() +
             "/Sector" + sector.str());
 
   // Build the histo name
   string histoName = string("TimeBox")
     + "_W" + wheel.str()
     + "_St" + station.str()
     + "_Sec" + sector.str()
     + "_SL" + superlayer.str();
 
    LogVerbatim("DTNoiseTask") <<"[DTNoiseTask]: booking SL histo:"<<endl;
    LogVerbatim("DTNoiseTask") <<"              folder "<< "DT/05-Noise/Wheel" + wheel.str() +
     "/Station" + station.str() +
     "/Sector" + sector.str() + "/" << endl;
    LogVerbatim("DTNoiseTask") <<"              histoName "<<histoName<<endl;
 
   tbHistos[slId] = ibooker.book1D(histoName,"Time Box (TDC counts)", 1000, 0, 6000);
 
 }

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

overrideprotected

Definition at line 233 of file DTNoiseTask.cc.

References dtGeom, edm::EventSetup::get(), and tTrigMap.

                                                                      {
 
   LogVerbatim("DTNoiseTask") <<"[DTNoiseTask]: Begin of run"<<endl;
 
   // tTrig Map
   setup.get<DTTtrigRcd>().get(tTrigMap);
 
   // get the geometry
   setup.get<MuonGeometryRecord>().get(dtGeom);
 
 }

void DTNoiseTask::endLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	setup
	)

overrideprotected

Definition at line 282 of file DTNoiseTask.cc.

References MillePedeFileConverter_cfg::e, evtNumber, mapEvt, noiseHistos, PostProcessor_cff::normalization, safeMargin, and tTrigStMap.

                                                                                             {
   LogVerbatim("DTNoiseTask") << "[DTNoiseTask]: End LS, update rates in all histos" << endl;
 
   // update the rate of all histos (usefull for histos with few entries: they are not updated very often
   for(map<DTChamberId, MonitorElement*>::const_iterator meAndChamber = noiseHistos.begin();
       meAndChamber != noiseHistos.end(); ++meAndChamber) {
     DTChamberId chId = (*meAndChamber).first;
     TH2F* noise_root = (*meAndChamber).second->getTH2F();
     double upperLimit = tTrigStMap[chId]-safeMargin;
 
     double normalization=0;
     if(mapEvt.find(chId) != mapEvt.end()) {
       LogVerbatim("DTNoiseTask")  << " Ch: " << chId << " Last fill: # of events: " << mapEvt[chId] << endl;
     normalization =  1e-9*upperLimit*mapEvt[chId];
     // revert back to # of entries
     noise_root->Scale(normalization);
     }
     //check that event analyzed != 0 might happen oline
     if (evtNumber) {
       // set the # of events analyzed until this update
       LogVerbatim("DTNoiseTask")  << "          Update for events: " << evtNumber << endl;
       mapEvt[chId] = evtNumber;
       // update the rate
       normalization = double( 1e-9*upperLimit*evtNumber);
       noise_root->Scale(1./normalization);
     }
   }
 }