#include <DTLocalTriggerLutTask.h>

Inheritance diagram for DTLocalTriggerLutTask:

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

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

	~DTLocalTriggerLutTask () 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	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	BeginRun. More...

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

void	searchTMBestIn (std::vector< L1MuDTChambPhDigi > const *trigs)
	Find best (highest qual) TM trigger segments. More...

void	searchTMBestOut (std::vector< L1MuDTChambPhDigi > const *trigs)

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 Member Functions
void	bookHistos (DQMStore::IBooker &ibooker, DTChamberId chId)
	Book histos. More...

std::string &	topFolder ()
	Get the top folder. More...

Private Attributes
std::string	baseFolder

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

bool	detailedAnalysis

edm::ESHandle< DTGeometry >	muonGeom

int	nEvents

int	nLumis

int	nPhibBins

int	nPhiBins

bool	overUnderIn

edm::ParameterSet	parameters

double	rangePhi

double	rangePhiB

edm::EDGetTokenT< DTRecSegment4DCollection >	seg_Token_

std::string	theGeomLabel

edm::EDGetTokenT< L1MuDTChambPhContainer >	tm_TokenIn_

edm::EDGetTokenT< L1MuDTChambPhContainer >	tm_TokenOut_

DTArr3bool	track_ok

DTArr3Digi	trigBestIn

DTArr3Digi	trigBestOut

DTTrigGeomUtils *	trigGeomUtils

DTArr3int	trigQualBestIn

DTArr3int	trigQualBestOut

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

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 44 of file DTLocalTriggerLutTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 34 of file DTLocalTriggerLutTask.cc.

References baseFolder, detailedAnalysis, edm::ParameterSet::getUntrackedParameter(), LogTrace, nEvents, nLumis, nPhibBins, nPhiBins, overUnderIn, rangePhi, rangePhiB, seg_Token_, theGeomLabel, tm_TokenIn_, and tm_TokenOut_.

                                                                       : trigGeomUtils(nullptr) {
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerLutTask") << "[DTLocalTriggerLutTask]: Constructor" << endl;
 
   tm_TokenIn_ = consumes<L1MuDTChambPhContainer>(ps.getUntrackedParameter<InputTag>("inputTagTMin"));
   tm_TokenOut_ = consumes<L1MuDTChambPhContainer>(ps.getUntrackedParameter<InputTag>("inputTagTMout"));
   seg_Token_ = consumes<DTRecSegment4DCollection>(ps.getUntrackedParameter<InputTag>("inputTagSEG"));
 
   overUnderIn = ps.getUntrackedParameter<bool>("rebinOutFlowsInGraph");
   detailedAnalysis = ps.getUntrackedParameter<bool>("detailedAnalysis");
   theGeomLabel = ps.getUntrackedParameter<string>("geomLabel");
 
   if (detailedAnalysis) {
     nPhiBins = 401;
     rangePhi = 10.025;
     nPhibBins = 401;
     rangePhiB = 10.025;
   } else {
     nPhiBins = 51;
     rangePhi = 5.1;
     nPhibBins = 51;
     rangePhiB = 10.2;
   }
 
   baseFolder = "DT/03-LocalTrigger-TM/";
   parameters = ps;
 
   nEvents = 0;
   nLumis = 0;
 }

DTLocalTriggerLutTask::~DTLocalTriggerLutTask ( )

override

Destructor.

Definition at line 64 of file DTLocalTriggerLutTask.cc.

References LogTrace, nEvents, and trigGeomUtils.

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

Member Function Documentation

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

overrideprotected

Analyze.

Definition at line 175 of file DTLocalTriggerLutTask.cc.

References chHistos, DTTrigGeomUtils::computeSCCoordinates(), detailedAnalysis, edm::Event::getByToken(), L1MuDTChambPhContainer::getContainer(), SiStripPI::max, min(), nEvents, overUnderIn, rangePhi, rangePhiB, DetId::rawId(), searchTMBestIn(), searchTMBestOut(), seg_Token_, DTChamberId::station(), relativeConstraints::station, tm_TokenIn_, tm_TokenOut_, trigBestIn, trigBestOut, DTTrigGeomUtils::trigDir(), trigGeomUtils, DTTrigGeomUtils::trigPos(), trigQualBestIn, trigQualBestOut, DTTrigGeomUtils::trigToSeg(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                              {
   nEvents++;
 
   edm::Handle<L1MuDTChambPhContainer> trigHandleIn;
   e.getByToken(tm_TokenIn_, trigHandleIn);
   edm::Handle<L1MuDTChambPhContainer> trigHandleOut;
   e.getByToken(tm_TokenOut_, trigHandleOut);
 
   vector<L1MuDTChambPhDigi> const* trigsIn = trigHandleIn->getContainer();
   searchTMBestIn(trigsIn);
   vector<L1MuDTChambPhDigi> const* trigsOut = trigHandleOut->getContainer();
   searchTMBestOut(trigsOut);
 
   Handle<DTRecSegment4DCollection> segments4D;
   e.getByToken(seg_Token_, segments4D);
   DTRecSegment4DCollection::id_iterator chamberId;
 
   // Preliminary loop finds best 4D Segment and high quality ones
   vector<const DTRecSegment4D*> best4DSegments;
 
   for (chamberId = segments4D->id_begin(); chamberId != segments4D->id_end(); ++chamberId) {
     DTRecSegment4DCollection::range rangeInCh = segments4D->get(*chamberId);
     DTRecSegment4DCollection::const_iterator trackIt = rangeInCh.first;
     DTRecSegment4DCollection::const_iterator trackEnd = rangeInCh.second;
 
     const DTRecSegment4D* tmpBest = nullptr;
     int tmpdof = 0;
     int dof = 0;
 
     for (; trackIt != trackEnd; ++trackIt) {
       if (trackIt->hasPhi()) {
         dof = trackIt->phiSegment()->degreesOfFreedom();
         if (dof > tmpdof) {
           tmpBest = &(*trackIt);
           tmpdof = dof;
         }
       }
     }
 
     if (tmpBest)
       best4DSegments.push_back(tmpBest);
   }
 
   vector<const DTRecSegment4D*>::const_iterator bestTrackIt = best4DSegments.begin();
   vector<const DTRecSegment4D*>::const_iterator bestTrackEnd = best4DSegments.end();
 
   for (; bestTrackIt != bestTrackEnd; ++bestTrackIt) {
     if ((*bestTrackIt)->hasPhi()) {
       DTChamberId chId = (*bestTrackIt)->chamberId();
       int nHitsPhi = (*bestTrackIt)->phiSegment()->degreesOfFreedom() + 2;
 
       int wheel = chId.wheel();
       int station = chId.station();
       int scsector = 0;
       float trackPosPhi, trackPosEta, trackDirPhi, trackDirEta;
       trigGeomUtils->computeSCCoordinates((*bestTrackIt), scsector, trackPosPhi, trackDirPhi, trackPosEta, trackDirEta);
 
       map<string, MonitorElement*>& chMap = chHistos[chId.rawId()];
       // In part
       if (trigQualBestIn[wheel + 3][station][scsector] > 3 &&  // residuals only for correlate triggers
           trigQualBestIn[wheel + 3][station][scsector] < 7 && nHitsPhi >= 7) {
         float trigPos = trigGeomUtils->trigPos(trigBestIn[wheel + 3][station][scsector]);
         float trigDir = trigGeomUtils->trigDir(trigBestIn[wheel + 3][station][scsector]);
         trigGeomUtils->trigToSeg(station, trigPos, trackDirPhi);
 
         double deltaPos = trigPos - trackPosPhi;
         deltaPos = overUnderIn ? max(min(deltaPos, rangePhi - 0.01), -rangePhi + 0.01) : deltaPos;
         double deltaDir = trigDir - trackDirPhi;
         deltaDir = overUnderIn ? max(min(deltaDir, rangePhiB - 0.01), -rangePhiB + 0.01) : deltaDir;
         chMap.find("TM_PhiResidualIn")->second->Fill(deltaPos);
         chMap.find("TM_PhibResidualIn")->second->Fill(deltaDir);
 
         if (detailedAnalysis) {
           chMap.find("TM_PhitkvsPhitrig")->second->Fill(trigPos, trackPosPhi);
           chMap.find("TM_PhibtkvsPhibtrig")->second->Fill(trigDir, trackDirPhi);
           chMap.find("TM_PhibResidualvsTkPos")->second->Fill(trackPosPhi, trigDir - trackDirPhi);
           chMap.find("TM_PhiResidualvsTkPos")->second->Fill(trackPosPhi, trigPos - trackPosPhi);
         }
       }
 
       // Out part
       if (trigQualBestOut[wheel + 3][station][scsector] > 3 &&  // residuals only for correlate triggers
           trigQualBestOut[wheel + 3][station][scsector] < 7 && nHitsPhi >= 7) {
         float trigPos = trigGeomUtils->trigPos(trigBestOut[wheel + 3][station][scsector]);
         float trigDir = trigGeomUtils->trigDir(trigBestOut[wheel + 3][station][scsector]);
         trigGeomUtils->trigToSeg(station, trigPos, trackDirPhi);
 
         double deltaPos = trigPos - trackPosPhi;
         deltaPos = overUnderIn ? max(min(deltaPos, rangePhi - 0.01), -rangePhi + 0.01) : deltaPos;
         double deltaDir = trigDir - trackDirPhi;
         deltaDir = overUnderIn ? max(min(deltaDir, rangePhiB - 0.01), -rangePhiB + 0.01) : deltaDir;
         chMap.find("TM_PhiResidualOut")->second->Fill(deltaPos);
         chMap.find("TM_PhibResidualOut")->second->Fill(deltaDir);
       }
     }
   }
 }

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

overrideprotected

To reset the MEs.

Definition at line 156 of file DTLocalTriggerLutTask.cc.

References chHistos, LogTrace, and nLumis.

                                                                                                           {
   nLumis++;
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerLutTask") << "[DTLocalTriggerLutTask]: Begin of LS transition" << endl;
 
   if (nLumis % parameters.getUntrackedParameter<int>("ResetCycle") == 0) {
     LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerLutTask") << "[DTLocalTriggerLutTask]: Cleaning histos" << endl;
     map<uint32_t, map<string, MonitorElement*> >::const_iterator chambIt = chHistos.begin();
     map<uint32_t, map<string, MonitorElement*> >::const_iterator chambEnd = chHistos.end();
 
     for (; chambIt != chambEnd; ++chambIt) {
       map<string, MonitorElement*>::const_iterator histoIt = chambIt->second.begin();
       map<string, MonitorElement*>::const_iterator histoEnd = chambIt->second.end();
       for (; histoIt != histoEnd; ++histoIt) {
         histoIt->second->Reset();
       }
     }
   }
 }

void DTLocalTriggerLutTask::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	run,
		edm::EventSetup const &	context
	)

overridevirtual

bookHistograms

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

Definition at line 144 of file DTLocalTriggerLutTask.cc.

References bookHistos(), DTGeometry::chambers(), LogTrace, and muonGeom.

                                                                          {
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerLutTask") << "[DTLocalTriggerLutTask]: bookHistograms" << endl;
 
   std::vector<const DTChamber*>::const_iterator chambIt = muonGeom->chambers().begin();
   std::vector<const DTChamber*>::const_iterator chambEnd = muonGeom->chambers().end();
 
   for (; chambIt != chambEnd; ++chambIt)
     bookHistos(ibooker, (*chambIt)->id());
 }

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

private

Book histos.

Definition at line 78 of file DTLocalTriggerLutTask.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), chHistos, detailedAnalysis, nPhibBins, nPhiBins, rangePhi, rangePhiB, DetId::rawId(), DTChamberId::sector(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, topFolder(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms(), and topFolder().

                                                                                  {
   stringstream wheel;
   wheel << chId.wheel();
   stringstream sector;
   sector << chId.sector();
   stringstream station;
   station << chId.station();
 
   ibooker.setCurrentFolder(topFolder() + "Wheel" + wheel.str() + "/Sector" + sector.str() + "/Station" + station.str() +
                            "/Segment");
 
   string chTag = "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
   std::map<std::string, MonitorElement*>& chambMap = chHistos[chId.rawId()];
 
   string hName = "TM_PhiResidualIn";
 
   chambMap[hName] = ibooker.book1D(hName + chTag,
                                    "Trigger local position In - Segment local position (correlated triggers)",
                                    nPhiBins,
                                    -rangePhi,
                                    rangePhi);
 
   hName = "TM_PhiResidualOut";
 
   chambMap[hName] = ibooker.book1D(hName + chTag,
                                    "Trigger local position Out - Segment local position (correlated triggers)",
                                    nPhiBins,
                                    -rangePhi,
                                    rangePhi);
 
   hName = "TM_PhibResidualIn";
 
   chambMap[hName] = ibooker.book1D(hName + chTag,
                                    "Trigger local direction In - Segment local direction (correlated triggers)",
                                    nPhibBins,
                                    -rangePhiB,
                                    rangePhiB);
 
   hName = "TM_PhibResidualOut";
 
   chambMap[hName] = ibooker.book1D(hName + chTag,
                                    "Trigger local direction Out - Segment local direction (correlated triggers)",
                                    nPhibBins,
                                    -rangePhiB,
                                    rangePhiB);
 
   if (detailedAnalysis) {
     hName = "TM_PhitkvsPhitrig";
 
     chambMap[hName] =
         ibooker.book2D(hName + chTag, "Local position: segment vs trigger", 100, -500., 500., 100, -500., 500.);
     hName = "TM_PhibtkvsPhibtrig";
 
     chambMap[hName] =
         ibooker.book2D(hName + chTag, "Local direction : segment vs trigger", 200, -40., 40., 200, -40., 40.);
     hName = "TM_PhibResidualvsTkPos";
 
     chambMap[hName] =
         ibooker.book2D(hName + chTag, "Local direction residual vs Segment Position", 100, -500., 500., 200, -10., 10.);
     hName = "TM_PhiResidualvsTkPos";
 
     chambMap[hName] =
         ibooker.book2D(hName + chTag, "Local Position residual vs Segment Position", 100, -500., 500., 200, -10., 10.);
   }
 }

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

overrideprotected

BeginRun.

Definition at line 72 of file DTLocalTriggerLutTask.cc.

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

                                                                                        {
   // Get the geometry
   context.get<MuonGeometryRecord>().get(theGeomLabel, muonGeom);
   trigGeomUtils = new DTTrigGeomUtils(muonGeom);
 }

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

inlinefinalprotected

Definition at line 70 of file DTLocalTriggerLutTask.h.

70 {}

void DTLocalTriggerLutTask::searchTMBestIn ( std::vector< L1MuDTChambPhDigi > const * trigs )

protected

Find best (highest qual) TM trigger segments.

Definition at line 273 of file DTLocalTriggerLutTask.cc.

References dtResolutionTest_cfi::histoTag, fileinputsource_cfi::sec, trigBestIn, trigQualBestIn, and wheelArrayShift.

Referenced by analyze().

                                                                                     {
   string histoType;
   string histoTag;
 
   // define best quality trigger segment
   // start from 1 and zero is kept empty
   for (int st = 0; st <= 4; ++st)
     for (int wh = 0; wh <= 5; ++wh)
       for (int sec = 0; sec <= 12; ++sec)
         trigQualBestIn[wh][st][sec] = -1;
 
   vector<L1MuDTChambPhDigi>::const_iterator trigIt = trigs->begin();
   vector<L1MuDTChambPhDigi>::const_iterator trigEnd = trigs->end();
   for (; trigIt != trigEnd; ++trigIt) {
     int wh = trigIt->whNum();
     int sec = trigIt->scNum() + 1;  // DTTF -> DT sector range transform
     int st = trigIt->stNum();
     int qual = trigIt->code();
 
     if (qual > trigQualBestIn[wh + wheelArrayShift][st][sec] && qual < 7) {
       trigQualBestIn[wh + wheelArrayShift][st][sec] = qual;
       trigBestIn[wh + wheelArrayShift][st][sec] = &(*trigIt);
     }
   }
 }

void DTLocalTriggerLutTask::searchTMBestOut ( std::vector< L1MuDTChambPhDigi > const * trigs )

protected

Definition at line 299 of file DTLocalTriggerLutTask.cc.

References dtResolutionTest_cfi::histoTag, fileinputsource_cfi::sec, trigBestOut, trigQualBestOut, and wheelArrayShift.

Referenced by analyze().

                                                                                      {
   string histoType;
   string histoTag;
 
   // define best quality trigger segment
   //   // start from 1 and zero is kept empty
   for (int st = 0; st <= 4; ++st)
     for (int wh = 0; wh <= 5; ++wh)
       for (int sec = 0; sec <= 12; ++sec)
         trigQualBestOut[wh][st][sec] = -1;
 
   vector<L1MuDTChambPhDigi>::const_iterator trigIt = trigs->begin();
   vector<L1MuDTChambPhDigi>::const_iterator trigEnd = trigs->end();
   for (; trigIt != trigEnd; ++trigIt) {
     int wh = trigIt->whNum();
     int sec = trigIt->scNum() + 1;  // DTTF -> DT sector range transform
     int st = trigIt->stNum();
     int qual = trigIt->code();
 
     if (qual > trigQualBestOut[wh + wheelArrayShift][st][sec] && qual < 7) {
       trigQualBestOut[wh + wheelArrayShift][st][sec] = qual;
       trigBestOut[wh + wheelArrayShift][st][sec] = &(*trigIt);
     }
   }
 }