#include <DTLocalTriggerTask.h>

Inheritance diagram for DTLocalTriggerTask:

Public Member Functions
	DTLocalTriggerTask (const edm::ParameterSet &ps)
	Constructor. More...

	~DTLocalTriggerTask () override
	Destructor. More...

Public Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

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

	DQMOneEDAnalyzer ()

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::Accumulator, Args... >
	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	analyze (const edm::Event &e, const edm::EventSetup &c) override
	Analyze. More...

void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context) override
	To reset the MEs. More...

void	bookBarrelHistos (DQMStore::IBooker &, std::string histoTag)
	Book the histograms. More...

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	Book the histograms. More...

void	bookHistos (DQMStore::IBooker &, const DTChamberId &dtCh, std::string folder, std::string histoTag)

void	bookWheelHistos (DQMStore::IBooker &, int wh, std::string histoTag)
	Book the histograms. More...

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

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

void	runSegmentAnalysis (edm::Handle< DTRecSegment4DCollection > &segments4D)
	Run analysis using DT 4D segments. More...

void	runTMAnalysis (std::vector< L1MuDTChambPhDigi > const phTrigs, std::vector< L1MuDTChambThDigi > const thTrigs)
	Run analysis on TM data. More...

void	setQLabels (MonitorElement *me, short int iaxis)
	Set Quality labels. More...

void	setQLabelsTheta (MonitorElement *me, short int iaxis)

std::string &	topFolder ()
	Get the Top folder (different between Physics and TP and TM) More...

void	triggerSource (const edm::Event &e)
	Get the L1A source. More...

Protected Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

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

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)

Protected Attributes
const int	wheelArrayShift = 3

Protected Attributes inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
edm::EDPutTokenT< DQMToken >	runToken_

Private Attributes
std::string	baseFolderTM

bool	detailedAnalysis

std::map< uint32_t, std::map< std::string, MonitorElement * > >	digiHistos

bool	doTMTheta

DTArr3PhDigi	iphbest

bool	isLocalRun

DTArr3ThDigi	ithbest

edm::EDGetTokenT< LTCDigiCollection >	ltcDigiCollectionToken_

DTArr3mapInt	mapDTTF

edm::ESHandle< DTGeometry >	muonGeom

int	nevents

edm::ParameterSet	parameters

DTArr3int	phcode_best

edm::EDGetTokenT< DTLocalTriggerCollection >	ros_Token_

edm::EDGetTokenT< DTRecSegment4DCollection >	seg_Token_

DTArr3int	thcode_best

MonitorElement *	tm_IDDataErrorPlot

edm::EDGetTokenT< L1MuDTChambPhContainer >	tm_Token_

edm::EDGetTokenT< L1MuDTChambThContainer >	tmTh_Token_

bool	tpMode

bool	track_ok [6][5][15]

DTTrigGeomUtils *	trigGeomUtils

std::string	trigsrc

bool	useSEG

bool	useTM

std::map< int, std::map< std::string, MonitorElement * > >	wheelHistos

Friends
class	DTMonitorModule

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
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)

Detailed Description

Definition at line 52 of file DTLocalTriggerTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 34 of file DTLocalTriggerTask.cc.

References baseFolderTM, detailedAnalysis, doTMTheta, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_2018_cff::InputTag, isLocalRun, LogTrace, ltcDigiCollectionToken_, ros_Token_, seg_Token_, tm_Token_, tmTh_Token_, and tpMode.

     : trigGeomUtils(nullptr), isLocalRun(ps.getUntrackedParameter<bool>("localrun", true)) {
   if (!isLocalRun) {
     ltcDigiCollectionToken_ = consumes<LTCDigiCollection>(ps.getParameter<edm::InputTag>("ltcDigiCollectionTag"));
   }
 
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask") << "[DTLocalTriggerTask]: Constructor" << endl;
 
   tpMode = ps.getUntrackedParameter<bool>("testPulseMode", false);
   detailedAnalysis = ps.getUntrackedParameter<bool>("detailedAnalysis", false);
   doTMTheta = ps.getUntrackedParameter<bool>("enableTMTheta", false);
   tm_Token_ = consumes<L1MuDTChambPhContainer>(
       ps.getUntrackedParameter<InputTag>("tm_label", InputTag("twinMuxStage2Digis:PhIn")));
   tmTh_Token_ = consumes<L1MuDTChambThContainer>(
       ps.getUntrackedParameter<edm::InputTag>("tmTh_label", InputTag("twinMuxStage2Digis:ThIn")));
 
   ros_Token_ =
       consumes<DTLocalTriggerCollection>(ps.getUntrackedParameter<InputTag>("ros_label", InputTag("dtunpacker")));
   seg_Token_ =
       consumes<DTRecSegment4DCollection>(ps.getUntrackedParameter<InputTag>("seg_label", InputTag("dt4DSegments")));
 
   if (tpMode) {
     baseFolderTM = "DT/11-LocalTriggerTP-TM/";
   } else {
     baseFolderTM = "DT/03-LocalTrigger-TM/";
   }
 
   parameters = ps;
 }

DTLocalTriggerTask::~DTLocalTriggerTask ( )

override

Destructor.

Definition at line 64 of file DTLocalTriggerTask.cc.

References LogTrace, nevents, and trigGeomUtils.

                                         {
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask")
       << "[DTLocalTriggerTask]: analyzed " << nevents << " events" << endl;
   if (trigGeomUtils) {
     delete trigGeomUtils;
   }
 }

Member Function Documentation

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

overrideprotected

Analyze.

Definition at line 194 of file DTLocalTriggerTask.cc.

References edm::Event::getByToken(), L1MuDTChambPhContainer::getContainer(), L1MuDTChambThContainer::getContainer(), edm::HandleBase::isValid(), nevents, runSegmentAnalysis(), runTMAnalysis(), seg_Token_, tm_Token_, tmTh_Token_, tpMode, triggerSource(), useSEG, and useTM.

                                                                           {
   if (!nevents) {
     edm::Handle<L1MuDTChambPhContainer> l1DTTPGPh;
     e.getByToken(tm_Token_, l1DTTPGPh);
     edm::Handle<L1MuDTChambThContainer> l1DTTPGTh;
     e.getByToken(tmTh_Token_, l1DTTPGTh);
 
     useTM = (l1DTTPGPh.isValid() || l1DTTPGTh.isValid()) && parameters.getUntrackedParameter<bool>("process_tm", true);
 
     Handle<DTRecSegment4DCollection> all4DSegments;
     e.getByToken(seg_Token_, all4DSegments);
     useSEG = all4DSegments.isValid() && parameters.getUntrackedParameter<bool>("process_seg", true);
   }
 
   nevents++;
 
   triggerSource(e);
 
   if (useTM) {
     edm::Handle<L1MuDTChambPhContainer> l1DTTPGPh;
     e.getByToken(tm_Token_, l1DTTPGPh);
     vector<L1MuDTChambPhDigi> const* l1PhTrig = l1DTTPGPh->getContainer();
 
     edm::Handle<L1MuDTChambThContainer> l1DTTPGTh;
     e.getByToken(tmTh_Token_, l1DTTPGTh);
 
     vector<L1MuDTChambThDigi> const* l1ThTrig = l1DTTPGTh->getContainer();
 
     runTMAnalysis(l1PhTrig, l1ThTrig);
   }
   if (!tpMode && useSEG) {
     Handle<DTRecSegment4DCollection> segments4D;
     e.getByToken(seg_Token_, segments4D);
 
     runSegmentAnalysis(segments4D);
   }
 }

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

overrideprotected

To reset the MEs.

Definition at line 179 of file DTLocalTriggerTask.cc.

References digiHistos, timingPdfMaker::histo, edm::LuminosityBlockBase::id(), LogTrace, edm::LuminosityBlockID::luminosityBlock(), and genParticles_cff::map.

                                                                                                        {
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask") << "[DTLocalTriggerTask]: Begin of LS transition" << endl;
 
   if (lumiSeg.id().luminosityBlock() % parameters.getUntrackedParameter<int>("ResetCycle", 3) == 0) {
     for (map<uint32_t, map<string, MonitorElement*> >::const_iterator histo = digiHistos.begin();
          histo != digiHistos.end();
          histo++) {
       for (map<string, MonitorElement*>::const_iterator ht = (*histo).second.begin(); ht != (*histo).second.end();
            ht++) {
         (*ht).second->Reset();
       }
     }
   }
 }

void DTLocalTriggerTask::bookBarrelHistos	(	DQMStore::IBooker &	,
		std::string	histoTag
	)

protected

Book the histograms.

Definition at line 232 of file DTLocalTriggerTask.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::impl::MonitorElement::setAxisTitle(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), tm_IDDataErrorPlot, and topFolder().

Referenced by bookHistograms().

                                                                                    {
   ibooker.setCurrentFolder(topFolder());
   if (histoTag == "TM_ErrorsChamberID") {
     tm_IDDataErrorPlot = ibooker.book1D(histoTag.c_str(), "TM Data ID Error", 5, -2, 3);
     tm_IDDataErrorPlot->setAxisTitle("wheel", 1);
   }
 
   return;
 }

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

overrideprotectedvirtual

Book the histograms.

Implements DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 78 of file DTLocalTriggerTask.cc.

References bookBarrelHistos(), bookHistos(), detailedAnalysis, doTMTheta, LogTrace, hgcalPlots::stat, and tpMode.

                                                                       {
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask") << "[DTLocalTriggerTask]: bookHistograms" << endl;
 
   if (parameters.getUntrackedParameter<bool>("staticBooking", true)) {  // Static histo booking
 
     vector<string> trigSources;
     if (parameters.getUntrackedParameter<bool>("localrun", true)) {
       trigSources.push_back("");
     } else {
       trigSources.push_back("_DTonly");
       trigSources.push_back("_NoDT");
       trigSources.push_back("_DTalso");
     }
     vector<string>::const_iterator trigSrcIt = trigSources.begin();
     vector<string>::const_iterator trigSrcEnd = trigSources.end();
 
     if (parameters.getUntrackedParameter<bool>("process_tm", true)) {
       bookBarrelHistos(ibooker, "TM_ErrorsChamberID");
     }
 
     if (tpMode) {
       for (int stat = 1; stat < 5; ++stat) {
         for (int wh = -2; wh < 3; ++wh) {
           for (int sect = 1; sect < 13; ++sect) {
             DTChamberId dtChId(wh, stat, sect);
 
             if (parameters.getUntrackedParameter<bool>("process_tm", true)) {  // TM data
               bookHistos(ibooker, dtChId, "LocalTriggerPhiIn", "TM_BXvsQual" + (*trigSrcIt));
               bookHistos(ibooker, dtChId, "LocalTriggerPhiIn", "TM_QualvsPhirad" + (*trigSrcIt));
             }
           }
         }
       }  // end of loop
     } else {
       for (; trigSrcIt != trigSrcEnd; ++trigSrcIt) {
         for (int wh = -2; wh < 3; ++wh) {
           for (int sect = 1; sect < 13; ++sect) {
             for (int stat = 1; stat < 5; ++stat) {
               DTChamberId dtChId(wh, stat, sect);
               if (parameters.getUntrackedParameter<bool>("process_tm", true)) {  // TM data
 
                 bookHistos(ibooker, dtChId, "LocalTriggerPhiIn", "TM_BXvsQual_In" + (*trigSrcIt));
                 if (detailedAnalysis) {
                   bookHistos(ibooker, dtChId, "LocalTriggerPhiIn", "TM_QualvsPhirad_In" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "LocalTriggerPhiIn", "TM_QualvsPhibend_In" + (*trigSrcIt));
                 }
                 bookHistos(ibooker, dtChId, "LocalTriggerPhiIn", "TM_Flag1stvsQual_In" + (*trigSrcIt));
                 bookHistos(ibooker, dtChId, "LocalTriggerPhiIn", "TM_FlagUpDownvsQual_In" + (*trigSrcIt));
 
                 bookHistos(ibooker, dtChId, "LocalTriggerPhiIn", "TM_BestQual_In" + (*trigSrcIt));
                 if (stat != 4 && doTMTheta) {
                   bookHistos(ibooker, dtChId, "LocalTriggerTheta", "TM_PositionvsBX" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "LocalTriggerTheta", "TM_QualityvsBX" + (*trigSrcIt));
 
                   bookHistos(ibooker, dtChId, "LocalTriggerTheta", "TM_PositionvsQual" + (*trigSrcIt));
 
                   bookHistos(ibooker, dtChId, "LocalTriggerTheta", "TM_ThetaBXvsQual" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "LocalTriggerTheta", "TM_ThetaBestQual" + (*trigSrcIt));
                 }
 
                 if (parameters.getUntrackedParameter<bool>("process_seg", true)) {  // TM + Segemnt
                   bookHistos(ibooker, dtChId, "Segment", "TM_PhitkvsPhitrig" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "Segment", "TM_PhibtkvsPhibtrig" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "Segment", "TM_PhiResidual" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "Segment", "TM_PhiResidualvsLUTPhi" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "Segment", "TM_PhibResidual" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "Segment", "TM_HitstkvsQualtrig" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "Segment", "TM_TrackPosvsAngle" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "Segment", "TM_TrackPosvsAngleandTrig" + (*trigSrcIt));
                   bookHistos(ibooker, dtChId, "Segment", "TM_TrackPosvsAngleandTrigHHHL" + (*trigSrcIt));
                   if (stat != 4) {
                     bookHistos(ibooker, dtChId, "Segment", "TM_TrackThetaPosvsAngle" + (*trigSrcIt));  // theta view
                     bookHistos(ibooker, dtChId, "Segment", "TM_TrackThetaPosvsAngleandTrig" + (*trigSrcIt));
                     bookHistos(ibooker,
                                dtChId,
                                "Segment",
                                "TM_TrackThetaPosvsAngleandTrigH" + (*trigSrcIt));  // TM theta quality not available!
                   }
                 }
               }
             }
           }
           for (int sect = 13; sect < 15; ++sect) {
             DTChamberId dtChId(wh, 4, sect);
             if (parameters.getUntrackedParameter<bool>("process_tm", true) &&
                 parameters.getUntrackedParameter<bool>("process_seg", true)) {  // TM+SEG LUTs data
               bookHistos(ibooker, dtChId, "Segment", "TM_PhitkvsPhitrig" + (*trigSrcIt));
               bookHistos(ibooker, dtChId, "Segment", "TM_PhibtkvsPhibtrig" + (*trigSrcIt));
               bookHistos(ibooker, dtChId, "Segment", "TM_PhiResidual" + (*trigSrcIt));
               bookHistos(ibooker, dtChId, "Segment", "TM_PhiResidualvsLUTPhi" + (*trigSrcIt));
               bookHistos(ibooker, dtChId, "Segment", "TM_PhibResidual" + (*trigSrcIt));
             }
           }
         }
       }  // end of loop
     }
   }
 }

void DTLocalTriggerTask::bookHistos	(	DQMStore::IBooker &	,
		const DTChamberId &	dtCh,
		std::string	folder,
		std::string	histoTag
	)

protected

Definition at line 242 of file DTLocalTriggerTask.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), digiHistos, printsummarytable::folder, fftjetimagerecorder_cfi::histoLabel, HltBtagPostValidation_cff::histoName, createfilelist::int, LogTrace, SiStripPI::max, simKBmtfStubs_cfi::maxBX, min(), simKBmtfStubs_cfi::minBX, LaserClient_cfi::nbins, DTTrigGeomUtils::phiRange(), DetId::rawId(), SimDataFormats::CaloAnalysis::sc, DTChamberId::sector(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), setQLabels(), setQLabelsTheta(), DTChamberId::station(), relativeConstraints::station, DTTrigGeomUtils::thetaRange(), topFolder(), trigGeomUtils, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

                                                      {
   int wh = dtCh.wheel();
   int sc = dtCh.sector();
   stringstream wheel;
   wheel << wh;
   stringstream station;
   station << dtCh.station();
   stringstream sector;
   sector << sc;
 
   double minBX = 0;
   double maxBX = 0;
   int rangeBX = 0;
 
   string histoType = histoTag.substr(3, histoTag.find("_", 3) - 3);
 
   ibooker.setCurrentFolder(topFolder() + "Wheel" + wheel.str() + "/Sector" + sector.str() + "/Station" + station.str() +
                            "/" + folder);
 
   string histoName = histoTag + "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
 
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask")
       << "[DTLocalTriggerTask]: booking " << topFolder() << "Wheel" << wheel.str() << "/Sector" << sector.str()
       << "/Station" << station.str() << "/" << folder << "/" << histoName << endl;
 
   if (histoType.find("BX") != string::npos) {
     if (histoTag.substr(0, 2) == "TM") {
       minBX = parameters.getUntrackedParameter<int>("minBXTM", -2) - 0.5;
       maxBX = parameters.getUntrackedParameter<int>("maxBXTM", 2) + 0.5;
     }
     rangeBX = (int)(maxBX - minBX);
   }
 
   if (folder == "LocalTriggerPhiIn") {
     if (histoType == "BXvsQual") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "BX vs trigger quality", 7, -0.5, 6.5, rangeBX, minBX, maxBX);
       setQLabels((digiHistos[dtCh.rawId()])[histoTag], 1);
       return;
     }
     if (histoType == "BestQual") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book1D(histoName, "Trigger quality of best primitives", 7, -0.5, 6.5);
       setQLabels((digiHistos[dtCh.rawId()])[histoTag], 1);
       return;
     }
     if (histoType == "QualvsPhirad") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "Trigger quality vs local position", 100, -500., 500., 7, -0.5, 6.5);
       setQLabels((digiHistos[dtCh.rawId()])[histoTag], 2);
       return;
     }
     if (histoType == "QualvsPhibend") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "Trigger quality vs local direction", 200, -40., 40., 7, -0.5, 6.5);
       setQLabels((digiHistos[dtCh.rawId()])[histoTag], 2);
       return;
     }
     if (histoType == "Flag1stvsQual") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "1st/2nd trig flag vs quality", 7, -0.5, 6.5, 2, -0.5, 1.5);
       setQLabels((digiHistos[dtCh.rawId()])[histoTag], 1);
       return;
     }
     if (histoType == "FlagUpDownvsQual") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "Up/Down trig flag vs quality", 7, -0.5, 6.5, 2, -0.5, 1.5);
       setQLabels((digiHistos[dtCh.rawId()])[histoTag], 1);
       return;
     }
 
   } else if (folder == "LocalTriggerTheta") {
     if (histoType == "PositionvsBX") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "Theta trigger position vs BX", rangeBX, minBX, maxBX, 7, -0.5, 6.5);
       return;
     }
     if (histoType == "PositionvsQual") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "Theta trigger position vs quality", 2, 0.5, 2.5, 7, -0.5, 6.5);
       setQLabelsTheta((digiHistos[dtCh.rawId()])[histoTag], 1);
       return;
     }
     if (histoType == "ThetaBXvsQual") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "BX vs trigger quality", 2, 0.5, 2.5, rangeBX, minBX, maxBX);
       setQLabelsTheta((digiHistos[dtCh.rawId()])[histoTag], 1);
     }
     if (histoType == "ThetaBestQual") {  //It does not make sense: 7 positions, how to define best?
       (digiHistos[dtCh.rawId()])[histoTag] = ibooker.book1D(
           histoName, "Trigger quality of best primitives (theta)", 2, 0.5, 2.5);  // 0 = not fired, 1 = L, 2 = H
       setQLabelsTheta((digiHistos[dtCh.rawId()])[histoTag], 1);
       return;
     }
 
   } else if (folder == "Segment") {
     if (histoType.find("TrackThetaPosvsAngle") == 0) {
       string histoLabel = "Position vs Angle (theta)";
       if (histoType.find("andTrigH") != string::npos)
         histoLabel += " for H triggers";
       else if (histoType.find("andTrig") != string::npos)
         histoLabel += " for triggers";
 
       float min, max;
       int nbins;
       trigGeomUtils->thetaRange(dtCh, min, max, nbins);
       (digiHistos[dtCh.rawId()])[histoTag] = ibooker.book2D(histoName, histoLabel, 16, -40., 40., nbins, min, max);
       return;
     }
     if (histoType.find("TrackPosvsAngle") == 0) {
       float min, max;
       int nbins;
       trigGeomUtils->phiRange(dtCh, min, max, nbins);
 
       string histoLabel = "Position vs Angle (phi)";
       if (histoType.find("andTrigHHHL") != string::npos)
         histoLabel += " for HH/HL triggers";
       else if (histoType.find("andTrig") != string::npos)
         histoLabel += " for triggers";
 
       (digiHistos[dtCh.rawId()])[histoTag] = ibooker.book2D(histoName, histoLabel, 16, -40., 40., nbins, min, max);
       return;
     }
     if (histoType == "PhitkvsPhitrig") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "Local position: segment vs trigger", 100, -500., 500., 100, -500., 500.);
       return;
     }
     if (histoType == "PhibtkvsPhibtrig") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "Local direction : segment vs trigger", 200, -40., 40., 200, -40., 40.);
       return;
     }
     if (histoType == "PhiResidual") {
       (digiHistos[dtCh.rawId()])[histoTag] = ibooker.book1D(
           histoName, "Trigger local position - Segment local position (correlated triggers)", 400, -10., 10.);
       return;
     }
     if (histoType == "PhibResidual") {
       (digiHistos[dtCh.rawId()])[histoTag] = ibooker.book1D(
           histoName, "Trigger local direction - Segment local direction (correlated triggers)", 500, -10., 10.);
       return;
     }
     if (histoType == "HitstkvsQualtrig") {
       (digiHistos[dtCh.rawId()])[histoTag] =
           ibooker.book2D(histoName, "Segment hits (phi) vs trigger quality", 7, -0.5, 6.5, 10, 0.5, 10.5);
       setQLabels((digiHistos[dtCh.rawId()])[histoTag], 1);
       return;
     }
   }
 }

void DTLocalTriggerTask::bookWheelHistos	(	DQMStore::IBooker &	,
		int	wh,
		std::string	histoTag
	)

protected

Book the histograms.

Definition at line 397 of file DTLocalTriggerTask.cc.

References HltBtagPostValidation_cff::histoName, LogTrace, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), topFolder(), and makeMuonMisalignmentScenario::wheel.

                                                                                           {
   stringstream wheel;
   wheel << wh;
 
   string histoType = histoTag.substr(3, histoTag.find("_", 3) - 3);
 
   ibooker.setCurrentFolder(topFolder() + "Wheel" + wheel.str() + "/");
 
   string histoName = histoTag + "_W" + wheel.str();
 
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask")
       << "[DTLocalTriggerTask]: booking " << topFolder() << "Wheel" << wheel.str() << "/" << histoName << endl;
 }

void DTLocalTriggerTask::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	context
	)

overrideprotected

Beginrun.

Definition at line 72 of file DTLocalTriggerTask.cc.

References edm::EventSetup::get(), muonGeom, nevents, and trigGeomUtils.

                                                                                     {
   nevents = 0;
   context.get<MuonGeometryRecord>().get(muonGeom);
   trigGeomUtils = new DTTrigGeomUtils(muonGeom);
 }

void DTLocalTriggerTask::endLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	context
	)

inlinefinalprotected

Definition at line 92 of file DTLocalTriggerTask.h.

References MillePedeFileConverter_cfg::e, and triggerSource().

92 {}

void DTLocalTriggerTask::runSegmentAnalysis ( edm::Handle< DTRecSegment4DCollection > & segments4D )

protected

Run analysis using DT 4D segments.

Definition at line 539 of file DTLocalTriggerTask.cc.

References DTTrigGeomUtils::computeSCCoordinates(), digiHistos, iphbest, phcode_best, FastTimerService_cff::range, DetId::rawId(), relativeConstraints::station, thcode_best, HLT_2018_cff::track, track_ok, DTTrigGeomUtils::trigDir(), trigGeomUtils, DTTrigGeomUtils::trigPos(), trigsrc, DTTrigGeomUtils::trigToSeg(), useTM, and makeMuonMisalignmentScenario::wheel.

Referenced by analyze().

                                                                                         {
   DTRecSegment4DCollection::const_iterator track;
 
   // Find best tracks & good tracks
   memset(track_ok, false, 450 * sizeof(bool));
 
   DTRecSegment4DCollection::id_iterator chamberId;
   vector<const DTRecSegment4D*> best4DSegments;
 
   // Preliminary loop finds best 4D Segment and high quality ones
   for (chamberId = segments4D->id_begin(); chamberId != segments4D->id_end(); ++chamberId) {
     DTRecSegment4DCollection::range range = segments4D->get(*chamberId);
     const DTRecSegment4D* tmpBest = nullptr;
     int tmpdof = 0;
     int dof = 0;
 
     for (track = range.first; track != range.second; ++track) {
       if ((*track).hasPhi()) {
         dof = (*track).phiSegment()->degreesOfFreedom();
         if (dof > tmpdof) {
           tmpBest = &(*track);
           tmpdof = dof;
 
           int wheel = (*track).chamberId().wheel();
           int sector = (*track).chamberId().sector();
           int station = (*track).chamberId().station();
           if (sector == 13) {
             sector = 4;
           } else if (sector == 14) {
             sector = 10;
           }
           track_ok[wheel + 3][station][sector] = (!track_ok[wheel + 3][station][sector] && dof >= 2);
         }
       }
     }
     if (tmpBest)
       best4DSegments.push_back(tmpBest);
   }
 
   vector<const DTRecSegment4D*>::const_iterator btrack;
 
   for (btrack = best4DSegments.begin(); btrack != best4DSegments.end(); ++btrack) {
     if ((*btrack)->hasPhi()) {  // Phi component
 
       int wheel = (*btrack)->chamberId().wheel();
       int station = (*btrack)->chamberId().station();
       int sector = (*btrack)->chamberId().sector();
       int scsector = 0;
       float x_track, y_track, x_angle, y_angle;
       trigGeomUtils->computeSCCoordinates((*btrack), scsector, x_track, x_angle, y_track, y_angle);
       int nHitsPhi = (*btrack)->phiSegment()->degreesOfFreedom() + 2;
 
       DTChamberId dtChId(wheel, station, sector);  // get chamber for LUTs histograms (Sectors 1 to 14)
       uint32_t indexCh = dtChId.rawId();
       map<string, MonitorElement*>& innerMECh = digiHistos[indexCh];
 
       DTChamberId dtChIdSC =
           DTChamberId(wheel, station, scsector);  // get chamber for histograms SC granularity (sectors 1 to 12)
       indexCh = dtChIdSC.rawId();
       map<string, MonitorElement*>& innerME = digiHistos[indexCh];
 
       if (useTM && phcode_best[wheel + 3][station][scsector] > -1 && phcode_best[wheel + 3][station][scsector] < 7) {
         innerME.find("TM_HitstkvsQualtrig" + trigsrc)->second->Fill(phcode_best[wheel + 3][station][scsector], nHitsPhi);
 
         if (phcode_best[wheel + 3][station][scsector] > 3 && nHitsPhi >= 7) {
           float x_trigger = trigGeomUtils->trigPos(iphbest[wheel + 3][station][scsector]);
           float angle_trigger = trigGeomUtils->trigDir(iphbest[wheel + 3][station][scsector]);
           trigGeomUtils->trigToSeg(station, x_trigger, x_angle);
 
           innerMECh.find("TM_PhitkvsPhitrig" + trigsrc)->second->Fill(x_trigger, x_track);
           innerMECh.find("TM_PhibtkvsPhibtrig" + trigsrc)->second->Fill(angle_trigger, x_angle);
           innerMECh.find("TM_PhiResidual" + trigsrc)->second->Fill(x_trigger - x_track);
           innerMECh.find("TM_PhibResidual" + trigsrc)->second->Fill(angle_trigger - x_angle);
         }
       }
 
       if (useTM) {
         // check for triggers elsewhere in the sector
         bool trigFlagTM = false;
         for (int ist = 1; ist < 5; ist++) {
           if (ist != station && phcode_best[wheel + 3][ist][scsector] >= 2 &&
               phcode_best[wheel + 3][ist][scsector] < 7 && track_ok[wheel + 3][ist][scsector] == true) {
             trigFlagTM = true;
             break;
           }
         }
 
         if (trigFlagTM && fabs(x_angle) < 40. && nHitsPhi >= 7) {
           // position vs angle of track for reconstruced tracks (denom. for trigger efficiency)
           innerME.find("TM_TrackPosvsAngle" + trigsrc)->second->Fill(x_angle, x_track);
           if (phcode_best[wheel + 3][station][scsector] >= 2 && phcode_best[wheel + 3][station][scsector] < 7) {
             innerME.find("TM_TrackPosvsAngleandTrig" + trigsrc)->second->Fill(x_angle, x_track);
             if (phcode_best[wheel + 3][station][scsector] > 4) {  //HH & HL Triggers
               innerME.find("TM_TrackPosvsAngleandTrigHHHL" + trigsrc)->second->Fill(x_angle, x_track);
             }
           }
         }
 
         if ((*btrack)->hasZed() && trigFlagTM && fabs(y_angle) < 40. &&
             (*btrack)->zSegment()->degreesOfFreedom() >= 1) {
           // position va angle of track for reconstruced tracks (denom. for trigger efficiency) along theta direction
           innerME.find("TM_TrackThetaPosvsAngle" + trigsrc)->second->Fill(y_angle, y_track);
           if (thcode_best[wheel + 3][station][scsector] > 0) {
             innerME.find("TM_TrackThetaPosvsAngleandTrig" + trigsrc)->second->Fill(y_angle, y_track);
             if (thcode_best[wheel + 3][station][scsector] == 2) {
               innerME.find("TM_TrackThetaPosvsAngleandTrigH" + trigsrc)->second->Fill(y_angle, y_track);
             }
           }
         }
       }
     }
   }
 }

void DTLocalTriggerTask::runTMAnalysis	(	std::vector< L1MuDTChambPhDigi > const *	phTrigs,
		std::vector< L1MuDTChambThDigi > const *	thTrigs
	)

protected

Run analysis on TM data.

Definition at line 411 of file DTLocalTriggerTask.cc.

References angle(), detailedAnalysis, digiHistos, doTMTheta, dqm::impl::MonitorElement::Fill(), dtResolutionTest_cfi::histoTag, mps_fire::i, triggerObjects_cff::id, iphbest, ithbest, dqmiolumiharvest::j, dqmdumpme::k, phcode_best, DetId::rawId(), SimDataFormats::CaloAnalysis::sc, thcode_best, tm_IDDataErrorPlot, tpMode, DTTrigGeomUtils::trigDir(), trigGeomUtils, DTTrigGeomUtils::trigPos(), trigsrc, wheelArrayShift, and x.

Referenced by analyze().

                                                                                     {
   //exit(0);
 
   string histoType;
   string histoTag;
 
   // define best quality trigger segment (phi and theta)
   // in any station start from 1 and zero is kept empty
   for (int i = 0; i < 5; ++i) {
     for (int j = 0; j < 6; ++j) {
       for (int k = 0; k < 13; ++k) {
         phcode_best[j][i][k] = -1;
         thcode_best[j][i][k] = -1;
       }
     }
   }
 
   vector<L1MuDTChambPhDigi>::const_iterator iph = phTrigs->begin();
   vector<L1MuDTChambPhDigi>::const_iterator iphe = phTrigs->end();
   for (; iph != iphe; ++iph) {
     int phwheel = iph->whNum();
     int phsec = iph->scNum() + 1;  // SM The track finder goes from 0 to 11. I need them from 1 to 12 !!!!!
     int phst = iph->stNum();
     int phbx = iph->bxNum();
     int phcode = iph->code();
     int phi1st = iph->Ts2Tag();
     int updown = iph->UpDownTag();
 
     // FIXME: workaround for TM data with station ID
     if (phst == 0) {
       tm_IDDataErrorPlot->Fill(phwheel);
       continue;
     }
 
     if (phcode > phcode_best[phwheel + 3][phst][phsec] && phcode < 7) {
       phcode_best[phwheel + 3][phst][phsec] = phcode;
       iphbest[phwheel + 3][phst][phsec] = &(*iph);
     }
 
     DTChamberId dtChId(phwheel, phst, phsec);
 
     float x = trigGeomUtils->trigPos(&(*iph));
     float angle = trigGeomUtils->trigDir(&(*iph));
     uint32_t indexCh = dtChId.rawId();
 
     map<string, MonitorElement*>& innerME = digiHistos[indexCh];
 
     if (tpMode) {
       innerME.find("TM_BXvsQual" + trigsrc)->second->Fill(phcode, phbx - phi1st);  // SM BX vs Qual Phi view (1st tracks)
       innerME.find("TM_QualvsPhirad" + trigsrc)->second->Fill(x, phcode);          // SM Qual vs radial angle Phi view
     } else {
       innerME.find("TM_BXvsQual_In" + trigsrc)
           ->second->Fill(phcode, phbx - phi1st);  // SM BX vs Qual Phi view (1st tracks)
       innerME.find("TM_Flag1stvsQual_In" + trigsrc)
           ->second->Fill(phcode, phi1st);  // SM Qual 1st/2nd track flag Phi view
       innerME.find("TM_FlagUpDownvsQual_In" + trigsrc)
           ->second->Fill(phcode, updown);  // SM Qual Up/Down track flag Phi view
       if (detailedAnalysis) {
         innerME.find("TM_QualvsPhirad_In" + trigsrc)->second->Fill(x, phcode);       // SM Qual vs radial angle Phi view
         innerME.find("TM_QualvsPhibend_In" + trigsrc)->second->Fill(angle, phcode);  // SM Qual vs bending Phi view
       }
     }
   }
 
   if (doTMTheta) {
     int thcode[7];
     vector<L1MuDTChambThDigi>::const_iterator ith = thTrigs->begin();
     vector<L1MuDTChambThDigi>::const_iterator ithe = thTrigs->end();
     for (; ith != ithe; ++ith) {
       int thwheel = ith->whNum();
       int thsec = ith->scNum() + 1;  // SM The track finder goes from 0 to 11. I need them from 1 to 12 !!!!!
       int thst = ith->stNum();
       int thbx = ith->bxNum();
 
       for (int pos = 0; pos < 7; pos++) {
         thcode[pos] = ith->code(pos);
 
         if (thcode[pos] > thcode_best[thwheel + 3][thst][thsec]) {
           thcode_best[thwheel + 3][thst][thsec] = thcode[pos];
           ithbest[thwheel + 3][thst][thsec] = &(*ith);
         }
       }
 
       DTChamberId dtChId(thwheel, thst, thsec);
       uint32_t indexCh = dtChId.rawId();
 
       map<string, MonitorElement*>& innerME = digiHistos[indexCh];
 
       for (int pos = 0; pos < 7; pos++) {  //SM fill position for non zero position bit in theta view
         if (thcode[pos] > 0) {             //Fired
           innerME.find("TM_PositionvsBX" + trigsrc)->second->Fill(thbx, pos);  // SM BX vs Position Theta view
           innerME.find("TM_PositionvsQual" + trigsrc)
               ->second->Fill(thcode[pos], pos);  //code = pos + qual; so 0, 1, 2 for 0, L, H resp.
           innerME.find("TM_ThetaBXvsQual" + trigsrc)
               ->second->Fill(thcode[pos], thbx);  //code = pos + qual; so 0, 1, 2 for 0, L, H resp.
         }
       }
     }
   }
 
   // Fill Quality plots with best TM triggers in phi & theta
   if (!tpMode) {
     for (int st = 1; st < 5; ++st) {
       for (int wh = -2; wh < 3; ++wh) {
         for (int sc = 1; sc < 13; ++sc) {
           if (phcode_best[wh + wheelArrayShift][st][sc] > -1 && phcode_best[wh + wheelArrayShift][st][sc] < 7) {
             DTChamberId id(wh, st, sc);
             uint32_t indexCh = id.rawId();
             map<string, MonitorElement*>& innerME = digiHistos[indexCh];
 
             innerME.find("TM_BestQual_In" + trigsrc)
                 ->second->Fill(phcode_best[wh + wheelArrayShift][st][sc]);  // Best Qual Trigger Phi view
           }
           if (thcode_best[wh + wheelArrayShift][st][sc] > 0 && thcode_best[wh + wheelArrayShift][st][sc] < 3) {
             DTChamberId id(wh, st, sc);
             uint32_t indexCh = id.rawId();
             map<string, MonitorElement*>& innerME = digiHistos[indexCh];
 
             innerME.find("TM_ThetaBestQual" + trigsrc)
                 ->second->Fill(thcode_best[wh + wheelArrayShift][st][sc]);  // Best Qual Trigger Theta view
           }
         }
       }
     }
   }
 }

void DTLocalTriggerTask::setQLabels	(	MonitorElement *	me,
		short int	iaxis
	)

protected

Set Quality labels.

Definition at line 653 of file DTLocalTriggerTask.cc.

References dqm::impl::MonitorElement::getTH1(), timingPdfMaker::histo, mps_fire::i, and SummaryClient_cfi::labels.

Referenced by bookHistos().

                                                                        {
   TH1* histo = me->getTH1();
   if (!histo)
     return;
 
   TAxis* axis = nullptr;
   if (iaxis == 1) {
     axis = histo->GetXaxis();
   } else if (iaxis == 2) {
     axis = histo->GetYaxis();
   }
   if (!axis)
     return;
 
   string labels[7] = {"LI", "LO", "HI", "HO", "LL", "HL", "HH"};
   int istart = axis->GetXmin() < -1 ? 2 : 1;
   for (int i = 0; i < 7; i++) {
     axis->SetBinLabel(i + istart, labels[i].c_str());
   }
 }

void DTLocalTriggerTask::setQLabelsTheta	(	MonitorElement *	me,
		short int	iaxis
	)

protected

Definition at line 674 of file DTLocalTriggerTask.cc.

References dqm::impl::MonitorElement::getTH1(), timingPdfMaker::histo, mps_fire::i, and SummaryClient_cfi::labels.

Referenced by bookHistos().

                                                                             {
   TH1* histo = me->getTH1();
   if (!histo)
     return;
 
   TAxis* axis = nullptr;
   if (iaxis == 1) {
     axis = histo->GetXaxis();
   } else if (iaxis == 2) {
     axis = histo->GetYaxis();
   }
   if (!axis)
     return;
 
   string labels[2] = {"L", "H"};
   int istart = axis->GetXmin() < -1 ? 2 : 1;
   for (int i = 0; i < 2; i++) {
     axis->SetBinLabel(i + istart, labels[i].c_str());
   }
 }

std::string& DTLocalTriggerTask::topFolder ( )

inlineprotected

Get the Top folder (different between Physics and TP and TM)

Definition at line 98 of file DTLocalTriggerTask.h.

References baseFolderTM.

Referenced by bookBarrelHistos(), bookHistos(), and bookWheelHistos().

98 { return baseFolderTM; }

DTLocalTriggerTask::baseFolderTM

std::string baseFolderTM

Definition: DTLocalTriggerTask.h:113

void DTLocalTriggerTask::triggerSource ( const edm::Event & e )

protected

Get the L1A source.

Definition at line 695 of file DTLocalTriggerTask.cc.

References edm::Event::getByToken(), mps_fire::i, isLocalRun, ltcDigiCollectionToken_, and trigsrc.

Referenced by analyze(), and endLuminosityBlock().

                                                         {
   if (!isLocalRun) {
     Handle<LTCDigiCollection> ltcdigis;
     e.getByToken(ltcDigiCollectionToken_, ltcdigis);
 
     for (std::vector<LTCDigi>::const_iterator ltc_it = ltcdigis->begin(); ltc_it != ltcdigis->end(); ltc_it++) {
       size_t otherTriggerSum = 0;
       for (size_t i = 1; i < 6; i++) {
         otherTriggerSum += size_t((*ltc_it).HasTriggered(i));
       }
       if ((*ltc_it).HasTriggered(0) && otherTriggerSum == 0)
         trigsrc = "_DTonly";
       else if (!(*ltc_it).HasTriggered(0))
         trigsrc = "_NoDT";
       else if ((*ltc_it).HasTriggered(0) && otherTriggerSum > 0)
         trigsrc = "_DTalso";
     }
     return;
   }
 
   trigsrc = "";
   return;
 }