#include <L1TDTTF.h>

Inheritance diagram for L1TDTTF:

Public Member Functions
	L1TDTTF (const edm::ParameterSet &ps)
	base services More...

	~L1TDTTF () override

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

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

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

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

virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

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

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

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::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	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

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

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

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)

Private Member Functions
void	bookEta (int wh, int &nbins, float &start, float &stop)

void	fillMEs (std::vector< L1MuDTTrackCand > const *trackContainer, std::vector< L1MuRegionalCand > &gmtDttfCands)

void	setQualLabel (MonitorElement *me, int axis)

void	setWheelLabel (MonitorElement *me)

Private Attributes
MonitorElement *	dttf_bx [6][12]

MonitorElement *	dttf_bx_2ndTrack [6][12]

MonitorElement *	dttf_eta [6][12]

MonitorElement *	dttf_eta_fine_fraction [6][12]

MonitorElement *	dttf_eta_wheel_2ndTrack [6]

MonitorElement *	dttf_gmt_ghost

MonitorElement *	dttf_gmt_match

MonitorElement *	dttf_gmt_missed

MonitorElement *	dttf_nTracksPerEv [6][12]

MonitorElement *	dttf_nTracksPerEvent_integ

MonitorElement *	dttf_nTracksPerEvent_wheel [6]

MonitorElement *	dttf_phi [6][12]

MonitorElement *	dttf_phi_eta_coarse_wheel [6]

MonitorElement *	dttf_phi_eta_fine_wheel [6]

MonitorElement *	dttf_phi_eta_wheel_2ndTrack [6]

MonitorElement *	dttf_phi_wheel_2ndTrack [6]

MonitorElement *	dttf_pt [6][12]

MonitorElement *	dttf_pt_wheel_2ndTrack [6]

MonitorElement *	dttf_q [6][12]

MonitorElement *	dttf_q_wheel_2ndTrack [6]

MonitorElement *	dttf_qual [6][12]

MonitorElement *	dttf_quality_summary_wheel_2ndTrack [6]

MonitorElement *	dttf_quality_wheel_2ndTrack [6]

MonitorElement *	dttf_spare

edm::InputTag	dttpgSource_

edm::InputTag	gmtSource_

edm::EDGetTokenT< L1MuGMTReadoutCollection >	gmtSourceToken_

std::string	l1tsubsystemfolder_

edm::InputTag	muonCollectionLabel_

edm::EDGetTokenT< reco::MuonCollection >	muonCollectionToken_

int	nev_

int	nev_dttf_

int	nev_dttf_track2_

int	numTracks [6][12]

bool	online_

std::string	outputFile_

edm::InputTag	trackInputTag_

edm::EDGetTokenT< L1MuDTTrackContainer >	trackInputToken_

bool	verbose_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 30 of file L1TDTTF.h.

Constructor & Destructor Documentation

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

base services

DT input output tracks GMT GlobalMuon try

Verbose?

Use ROOT Output?

Definition at line 193 of file L1TDTTF.cc.

References dttpgSource_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), gmtSourceToken_, HLT_2018_cff::InputTag, edm::InputTag::label(), muonCollectionToken_, outputFile_, edm::InputTag::process(), AlCaHLTBitMon_QueryRunRegistry::string, trackInputTag_, trackInputToken_, and verbose_.

     : dttpgSource_(ps.getParameter<edm::InputTag>("dttpgSource")),
       gmtSource_(ps.getParameter<edm::InputTag>("gmtSource")),
       muonCollectionLabel_(ps.getParameter<edm::InputTag>("MuonCollection")),
       l1tsubsystemfolder_(ps.getUntrackedParameter<std::string>("l1tSystemFolder", "L1T/L1TDTTF")),
       online_(ps.getUntrackedParameter<bool>("online", true)),
       verbose_(ps.getUntrackedParameter<bool>("verbose", false))
 
 {
   std::string trstring = dttpgSource_.label() + ":DATA:" + dttpgSource_.process();
   trackInputTag_ = edm::InputTag(trstring);
 
   if (verbose_)
     edm::LogInfo("L1TDTTF: constructor") << "Verbose enabled";
 
   if (ps.getUntrackedParameter<bool>("disableROOToutput", false)) {
     outputFile_ = "";
 
   } else {
     outputFile_ = ps.getUntrackedParameter<std::string>("outputFile", "");
     if (!outputFile_.empty()) {
       edm::LogInfo("L1TDTTF: constructor") << "L1T Monitoring histograms will be saved to " << outputFile_;
     }
   }
 
   //set Token(-s)
   trackInputToken_ = consumes<L1MuDTTrackContainer>(trackInputTag_);
   muonCollectionToken_ = consumes<reco::MuonCollection>(ps.getParameter<edm::InputTag>("MuonCollection"));
   gmtSourceToken_ = consumes<L1MuGMTReadoutCollection>(ps.getParameter<edm::InputTag>("gmtSource"));
 }

L1TDTTF::~L1TDTTF ( )

override

Nothing to destroy

Definition at line 227 of file L1TDTTF.cc.

227 {

229 }

Member Function Documentation

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

overrideprotected

counters

tracks handle

dttf counters

selection for offline

global muon selection plot

in case of problems accept all

take only bx=0

in Gmt but not in DTTF

in phys values double phi= dttfCand->phiValue(); int sector = 1 + (phi + 15)/30; /// in phys values

Per event summaries

Definition at line 521 of file L1TDTTF.cc.

References accept(), dttf_gmt_ghost, dttf_nTracksPerEv, dttf_nTracksPerEvent_integ, dttf_nTracksPerEvent_wheel, dttf_spare, dttpgSource_, PVValHelper::eta, dqm::impl::MonitorElement::Fill(), fillMEs(), L1MuDTTrackContainer::getContainer(), L1MuGMTReadoutCollection::getRecords(), gmtSource_, gmtSourceToken_, edm::HandleBase::isValid(), edm::InputTag::label(), muonCollectionLabel_, muonCollectionToken_, PDWG_BPHSkim_cff::muons, nev_, nev_dttf_, nev_dttf_track2_, numTracks, online_, phi, edm::InputTag::process(), edm::Handle< T >::product(), alignCSCRings::s, trackInputToken_, verbose_, w, and makeMuonMisalignmentScenario::wheel.

                                                                             {
   if (verbose_)
     edm::LogInfo("L1TDTTF::Analyze::start") << "#################### START";
 
   ++nev_;
   memset(numTracks, 0, 72 * sizeof(int));
 
   edm::Handle<L1MuDTTrackContainer> myL1MuDTTrackContainer;
   try {
     event.getByToken(trackInputToken_, myL1MuDTTrackContainer);
   } catch (cms::Exception& iException) {
     edm::LogError("L1TDTTF::analyze::DataNotFound") << "can't getByToken L1MuDTTrackContainer with label "
                                                     << dttpgSource_.label() << ":DATA:" << dttpgSource_.process();
     return;
   }
 
   if (!myL1MuDTTrackContainer.isValid()) {
     edm::LogError("L1TDTTF::analyze::DataNotFound")
         << "can't find L1MuDTTrackContainer with label " << dttpgSource_.label() << ":DATA:" << dttpgSource_.process();
     return;
   }
 
   L1MuDTTrackContainer::TrackContainer const* trackContainer = myL1MuDTTrackContainer->getContainer();
 
   if (!trackContainer->empty()) {
     ++nev_dttf_;
     if (trackContainer->size() > 1)
       ++nev_dttf_track2_;
   }
 
   bool accept = true;
   if (!online_) {
     try {
       edm::Handle<reco::MuonCollection> muons;
       event.getByToken(muonCollectionToken_, muons);
       accept = false;
       if (muons.isValid()) {
         for (reco::MuonCollection::const_iterator recoMu = muons->begin(); recoMu != muons->end(); ++recoMu) {
           if (fabs(recoMu->eta()) < 1.4) {
             if (verbose_) {
               edm::LogInfo("L1TDTTFClient::Analyze:GM") << "Found a global muon!";
             }
             accept = true;
             break;
           }
         }
 
         if (!accept) {
           dttf_spare->Fill(!trackContainer->empty() ? 1 : 0);
 
           if (verbose_) {
             edm::LogInfo("L1TDTTFClient::Analyze:GM") << "No global muons in this event!";
           }
 
         } else {
           dttf_spare->Fill(!trackContainer->empty() ? 2 : 3);
         }
 
       } else {
         accept = true;
         edm::LogWarning("L1TDTTFClient::Analyze:GM")
             << "Invalid MuonCollection with label " << muonCollectionLabel_.label();
       }
 
     } catch (cms::Exception& iException) {
       accept = true;
       edm::LogError("DataNotFound") << "Unable to getByToken MuonCollection with label "
                                     << muonCollectionLabel_.label();
     }
   }
 
   std::vector<L1MuRegionalCand> gmtBx0DttfCandidates;
 
   try {
     edm::Handle<L1MuGMTReadoutCollection> pCollection;
     event.getByToken(gmtSourceToken_, pCollection);
 
     if (!pCollection.isValid()) {
       edm::LogError("DataNotFound") << "can't find L1MuGMTReadoutCollection with label " << gmtSource_.label();
     }
 
     // get GMT readout collection
     L1MuGMTReadoutCollection const* gmtrc = pCollection.product();
     std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
 
     std::vector<L1MuGMTReadoutRecord>::const_iterator RRItr;
 
     for (RRItr = gmt_records.begin(); RRItr != gmt_records.end(); ++RRItr) {
       std::vector<L1MuRegionalCand> dttfCands = RRItr->getDTBXCands();
       std::vector<L1MuRegionalCand>::iterator dttfCand;
 
       for (dttfCand = dttfCands.begin(); dttfCand != dttfCands.end(); ++dttfCand) {
         if (dttfCand->empty())
           continue;
         if (RRItr->getBxInEvent())
           continue;
 
         //       dttf_gmt_ghost_phys->Fill( dttfCand->eta_packed(),
         //               dttfCand->phi_packed() );
         gmtBx0DttfCandidates.push_back(*dttfCand);
       }
     }
 
   } catch (cms::Exception& iException) {
     edm::LogError("DataNotFound") << "Unable to getByToken L1MuGMTReadoutCollection with label " << gmtSource_.label();
   }
 
   // fill MEs if all selections are passed
   if (accept)
     fillMEs(trackContainer, gmtBx0DttfCandidates);
 
   std::vector<L1MuRegionalCand>::iterator dttfCand;
   for (dttfCand = gmtBx0DttfCandidates.begin(); dttfCand != gmtBx0DttfCandidates.end(); ++dttfCand) {
     if (dttfCand->empty())
       continue;
 
     int phi = dttfCand->phi_packed();
     int sector = 1 + (phi + 6) / 12;
     if (sector > 12)
       sector -= 12;
     double eta = dttfCand->etaValue();
 
     int wheel = -3;
     if (eta < -0.74) {
       wheel = -2;
     } else if (eta < -0.3) {
       wheel = -1;
 
     } else if (eta < 0.3) {
       wheel = 0;
 
     } else if (eta < 0.74) {
       wheel = 1;
     } else {
       wheel = 2;
     }
 
     dttf_gmt_ghost->Fill(wheel, sector);
     // dttf_gmt_ghost_phys->Fill( dttfCand->eta_packed(),
     //                            dttfCand->phi_packed() );
   }
 
   int numTracksInt = 0;
 
   for (int w = 0; w < 6; ++w) {
     int numTracks_wh = 0;
     for (int s = 0; s < 12; ++s) {
       dttf_nTracksPerEv[w][s]->Fill(numTracks[w][s]);
 
       numTracks_wh += numTracks[w][s];
     }
 
     numTracksInt += numTracks_wh;
     dttf_nTracksPerEvent_wheel[w]->Fill(numTracks_wh);
   }
 
   dttf_nTracksPerEvent_integ->Fill(numTracksInt);
 }

void L1TDTTF::bookEta	(	int	wh,
		int &	nbins,
		float &	start,
		float &	stop
	)

private

Definition at line 892 of file L1TDTTF.cc.

Referenced by bookHistograms().

                                                                    {
   switch (wh) {
     case 0:
       start = 0;
       stop = 18;
       nbins = 18;
       break;  // N2
     case 1:
       start = 8;
       stop = 28;
       nbins = 20;
       break;  // N1
     case 2:
       start = 22;
       stop = 32;
       nbins = 10;
       break;  // N0
     case 3:
       start = 22;
       stop = 42;
       nbins = 20;
       break;  // P0
     case 4:
       start = 36;
       stop = 56;
       nbins = 20;
       break;  // P1
     case 5:
       start = 46;
       stop = 64;
       nbins = 18;
       break;  // P2
     default:
       start = 0;
       stop = 0;
       nbins = 0;
       break;  // BOH
   }
 }

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

overrideprotectedvirtual

testing purposes

histo name

histo title

DTTF Output (6 wheels)

Per wheel summaries

number of tracks per event per wheel

phi vs etafine - for each wheel

phi vs etacoarse - for each wheel

Per wheel summaries : 2ND_TRACK_ONLY

DTTF Tracks Quality distribution

quality per wheel 2ND TRACK

phi vs eta - for each wheel 2ND TRACK

DTTF Tracks eta distribution (Packed values)

DTTF Tracks Phi distribution (Packed values)

DTTF Tracks p_{T} distribution (Packed values)

DTTF Tracks Charge distribution

Go in detailed subfolders

number of tracks per event folder

BX_SECTORS for each wheel

CHARGE folder

PT folder

PHI folder

QUALITY folder

ETA folder

integrated values: always packed

Only for online: occupancy summary - reset

Implements DQMEDAnalyzer.

Definition at line 237 of file L1TDTTF.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), bookEta(), dttf_bx, dttf_bx_2ndTrack, dttf_eta, dttf_eta_fine_fraction, dttf_eta_wheel_2ndTrack, dttf_gmt_ghost, dttf_gmt_match, dttf_gmt_missed, dttf_nTracksPerEv, dttf_nTracksPerEvent_integ, dttf_nTracksPerEvent_wheel, dttf_phi, dttf_phi_eta_coarse_wheel, dttf_phi_eta_fine_wheel, dttf_phi_eta_wheel_2ndTrack, dttf_phi_wheel_2ndTrack, dttf_pt, dttf_pt_wheel_2ndTrack, dttf_q, dttf_q_wheel_2ndTrack, dttf_qual, dttf_quality_summary_wheel_2ndTrack, dttf_quality_wheel_2ndTrack, dttf_spare, dqm::impl::MonitorElement::getTH2F(), l1tsubsystemfolder_, LaserClient_cfi::nbins, nev_, nev_dttf_, nev_dttf_track2_, online_, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), setQualLabel(), setWheelLabel(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                           {
   nev_ = 0;
   nev_dttf_ = 0;
   nev_dttf_track2_ = 0;
 
   // get hold of back-end interface
 
   std::string dttf_trk_folder = l1tsubsystemfolder_;
 
   char hname[100];   
   char htitle[100];  
 
   float start = 0;
   float stop = 0;
   int nbins = 0;
 
   ibooker.setCurrentFolder(dttf_trk_folder);
 
   std::string wheelpath[6] = {
       "/02-WHEEL_N2", "/03-WHEEL_N1", "/04-WHEEL_N0", "/05-WHEEL_P0", "/06-WHEEL_P1", "/07-WHEEL_P2"};
 
   char c_whn[6][3] = {"N2", "N1", "N0", "P0", "P1", "P2"};
   // char bxn [3][3] = { "N1", "0", "P1" };
   // char bxn[3][25] = {"/BX_NONZERO_ONLY/BX_N1", "", "/BX_NONZERO_ONLY/BX_P1"};
 
   for (int iwh = 0; iwh < 6; ++iwh) {
     bookEta(iwh, nbins, start, stop);  
 
     std::string dttf_trk_folder_wheel = dttf_trk_folder + wheelpath[iwh];
     ibooker.setCurrentFolder(dttf_trk_folder_wheel);
 
     sprintf(hname, "dttf_01_nTracksPerEvent_wh%s", c_whn[iwh]);
     sprintf(htitle, "Wheel %s - Number Tracks Per Event", c_whn[iwh]);
     dttf_nTracksPerEvent_wheel[iwh] = ibooker.book1D(hname, htitle, 10, 0.5, 10.5);
     dttf_nTracksPerEvent_wheel[iwh]->setAxisTitle("# tracks/event", 1);
 
     sprintf(hname, "dttf_07_phi_vs_etaFine_wh%s", c_whn[iwh]);
     sprintf(htitle, "Wheel %s -   #eta-#phi DTTF Tracks occupancy (fine #eta only, unpacked values)", c_whn[iwh]);
     dttf_phi_eta_fine_wheel[iwh] = ibooker.book2D(hname, htitle, nbins, start - 0.5, stop - 0.5, 144, -6, 138);
     // 144, -0.5, 143.5);
 
     dttf_phi_eta_fine_wheel[iwh]->setAxisTitle("#eta", 1);
     dttf_phi_eta_fine_wheel[iwh]->setAxisTitle("#phi", 2);
 
     sprintf(hname, "dttf_08_phi_vs_etaCoarse_wh%s", c_whn[iwh]);
     sprintf(htitle, "Wheel %s -   #eta-#phi DTTF Tracks occupancy (coarse #eta only, unpacked values)", c_whn[iwh]);
     dttf_phi_eta_coarse_wheel[iwh] = ibooker.book2D(hname, htitle, nbins, start - 0.5, stop - 0.5, 144, -6, 138);
     // 144, -0.5, 143.5);
     dttf_phi_eta_coarse_wheel[iwh]->setAxisTitle("#eta", 1);
     dttf_phi_eta_coarse_wheel[iwh]->setAxisTitle("#phi", 2);
 
     std::string dttf_trk_folder_wheel_2ndtrack = dttf_trk_folder_wheel + "/2ND_TRACK_ONLY";
     ibooker.setCurrentFolder(dttf_trk_folder_wheel_2ndtrack);
 
     sprintf(hname, "dttf_04_quality_wh%s_2ndTrack", c_whn[iwh]);
     sprintf(htitle, "Wheel %s - 2nd Tracks Quality distribution", c_whn[iwh]);
     dttf_quality_wheel_2ndTrack[iwh] = ibooker.book1D(hname, htitle, 7, 1, 8);
     setQualLabel(dttf_quality_wheel_2ndTrack[iwh], 1);
 
     sprintf(hname, "dttf_05_quality_summary_wh%s_2ndTrack", c_whn[iwh]);
     sprintf(htitle, "Wheel %s - 2nd Tracks - Quality", c_whn[iwh]);
     dttf_quality_summary_wheel_2ndTrack[iwh] = ibooker.book2D(hname, htitle, 12, 1, 13, 7, 1, 8);
     dttf_quality_summary_wheel_2ndTrack[iwh]->setAxisTitle("Sector", 1);
     setQualLabel(dttf_quality_summary_wheel_2ndTrack[iwh], 2);
     // dttf_quality_summary_wheel_2ndTrack[iwh]->setAxisTitle("Quality", 2);
 
     sprintf(hname, "dttf_06_phi_vs_eta_wh%s_2ndTrack", c_whn[iwh]);
     sprintf(htitle, "Wheel %s -   #eta-#phi Distribution of DTTF 2nd Tracks", c_whn[iwh]);
 
     dttf_phi_eta_wheel_2ndTrack[iwh] = ibooker.book2D(hname, htitle, nbins, start - 0.5, stop - 0.5, 144, -6, 138);
     // 144, -0.5, 143.5);
     dttf_phi_eta_wheel_2ndTrack[iwh]->setAxisTitle("#eta", 1);
     dttf_phi_eta_wheel_2ndTrack[iwh]->setAxisTitle("#phi", 2);
 
     sprintf(hname, "dttf_07_eta_wh%s_2ndTrack", c_whn[iwh]);
     sprintf(htitle, "Wheel %s - DTTF 2nd Tracks #eta distribution (Packed values)", c_whn[iwh]);
     dttf_eta_wheel_2ndTrack[iwh] = ibooker.book1D(hname, htitle, 64, -0.5, 63.5);
     dttf_eta_wheel_2ndTrack[iwh]->setAxisTitle("#eta", 1);
 
     sprintf(hname, "dttf_08_phi_wh%s_2ndTrack", c_whn[iwh]);
     sprintf(htitle, "Wheel %s - DTTF 2nd Tracks Phi distribution (Packed values)", c_whn[iwh]);
     dttf_phi_wheel_2ndTrack[iwh] = ibooker.book1D(hname, htitle, 144, -6, 138.);
     dttf_phi_wheel_2ndTrack[iwh]->setAxisTitle("#phi", 1);
 
     sprintf(hname, "dttf_09_pt_wh%s_2ndTrack", c_whn[iwh]);
     sprintf(htitle, "Wheel %s - DTTF 2nd Tracks p_{T} distribution (Packed values)", c_whn[iwh]);
     dttf_pt_wheel_2ndTrack[iwh] = ibooker.book1D(hname, htitle, 32, -0.5, 31.5);
     dttf_pt_wheel_2ndTrack[iwh]->setAxisTitle("p_{T}", 1);
 
     sprintf(hname, "dttf_10_charge_wh%s_2ndTrack", c_whn[iwh]);
     sprintf(htitle, "Wheel %s - DTTF 2nd Tracks Charge distribution", c_whn[iwh]);
     dttf_q_wheel_2ndTrack[iwh] = ibooker.book1D(hname, htitle, 2, -0.5, 1.5);
     dttf_q_wheel_2ndTrack[iwh]->setAxisTitle("Charge", 1);
 
 
     std::string dttf_trk_folder_nTracksPerEvent = dttf_trk_folder_wheel + "/TracksPerEvent";
     ibooker.setCurrentFolder(dttf_trk_folder_nTracksPerEvent);
 
     for (int ise = 0; ise < 12; ++ise) {
       sprintf(hname, "dttf_nTracksPerEvent_wh%s_se%d", c_whn[iwh], ise + 1);
       sprintf(htitle, "Wheel %s Sector %d - Number of Tracks Per Event", c_whn[iwh], ise + 1);
       dttf_nTracksPerEv[iwh][ise] = ibooker.book1D(hname, htitle, 2, 0.5, 2.5);
       dttf_nTracksPerEv[iwh][ise]->setAxisTitle("# tracks/event", 1);
     }
 
     std::string dttf_trk_folder_wh_bxsec_all = dttf_trk_folder_wheel + "/BX_BySector";
     ibooker.setCurrentFolder(dttf_trk_folder_wh_bxsec_all);
 
     for (int ise = 0; ise < 12; ++ise) {
       sprintf(hname, "dttf_bx_wh%s_se%d", c_whn[iwh], ise + 1);
       sprintf(htitle, "Wheel %s Sector %d - BX Distribution", c_whn[iwh], ise + 1);
       dttf_bx[iwh][ise] = ibooker.book1D(hname, htitle, 3, -1.5, 1.5);
       dttf_bx[iwh][ise]->setAxisTitle("BX", 1);
     }
 
     std::string dttf_trk_folder_wh_bxsec_trk2 = dttf_trk_folder_wheel + "/BX_BySector/2ND_TRACK_ONLY";
     ibooker.setCurrentFolder(dttf_trk_folder_wh_bxsec_trk2);
 
     for (int ise = 0; ise < 12; ++ise) {
       sprintf(hname, "dttf_bx_2ndTrack_wh%s_se%d", c_whn[iwh], ise + 1);
       sprintf(htitle, "Wheel %s Sector %d - BX 2nd Tracks only", c_whn[iwh], ise + 1);
       dttf_bx_2ndTrack[iwh][ise] = ibooker.book1D(hname, htitle, 3, -1.5, 1.5);
       dttf_bx_2ndTrack[iwh][ise]->setAxisTitle("BX", 1);
     }
 
     std::string dttf_trk_folder_charge = dttf_trk_folder_wheel + "/Charge";
     ibooker.setCurrentFolder(dttf_trk_folder_charge);
 
     for (int ise = 0; ise < 12; ++ise) {
       sprintf(hname, "dttf_charge_wh%s_se%d", c_whn[iwh], ise + 1);
       sprintf(htitle, "Wheel %s Sector %d - Packed Charge", c_whn[iwh], ise + 1);
       dttf_q[iwh][ise] = ibooker.book1D(hname, htitle, 2, -0.5, 1.5);
       dttf_q[iwh][ise]->setAxisTitle("Charge", 1);
     }
 
     std::string dttf_trk_folder_pt = dttf_trk_folder_wheel + "/PT";
     ibooker.setCurrentFolder(dttf_trk_folder_pt);
 
     for (int ise = 0; ise < 12; ++ise) {
       sprintf(hname, "dttf_pt_wh%s_se%d", c_whn[iwh], ise + 1);
       sprintf(htitle, "Wheel %s Sector %d - Packed p_{T}", c_whn[iwh], ise + 1);
       dttf_pt[iwh][ise] = ibooker.book1D(hname, htitle, 32, -0.5, 31.5);
       dttf_pt[iwh][ise]->setAxisTitle("p_{T}", 1);
     }
 
     std::string dttf_trk_folder_phi = dttf_trk_folder_wheel + "/Phi";
     ibooker.setCurrentFolder(dttf_trk_folder_phi);
 
     for (int ise = 0; ise < 12; ++ise) {
       sprintf(hname, "dttf_phi_wh%s_se%d", c_whn[iwh], ise + 1);
       sprintf(htitle, "Wheel %s Sector %d - Packed Phi", c_whn[iwh], ise + 1);
       dttf_phi[iwh][ise] = ibooker.book1D(hname, htitle, 144, -6, 138);
       dttf_phi[iwh][ise]->setAxisTitle("#phi", 1);
       //dttf_phi[iwh][ise] = dbe_->book1D(title,title, 32,-16.5, 15.5);
     }
 
     std::string dttf_trk_folder_quality = dttf_trk_folder_wheel + "/Quality";
     ibooker.setCurrentFolder(dttf_trk_folder_quality);
 
     for (int ise = 0; ise < 12; ++ise) {
       sprintf(hname, "dttf_qual_wh%s_se%d", c_whn[iwh], ise + 1);
       sprintf(htitle, "Wheel %s Sector %d - Packed Quality", c_whn[iwh], ise + 1);
       dttf_qual[iwh][ise] = ibooker.book1D(hname, htitle, 7, 1, 8);
       dttf_qual[iwh][ise]->setAxisTitle("Quality", 1);
       setQualLabel(dttf_qual[iwh][ise], 1);
     }
 
     std::string dttf_trk_folder_eta = dttf_trk_folder_wheel + "/Eta";
     ibooker.setCurrentFolder(dttf_trk_folder_eta);
 
     for (int ise = 0; ise < 12; ++ise) {
       sprintf(hname, "dttf_eta_wh%s_se%d", c_whn[iwh], ise + 1);
       sprintf(htitle, "Wheel %s Sector %d - Packed #eta", c_whn[iwh], ise + 1);
       dttf_eta[iwh][ise] = ibooker.book1D(hname, htitle, 64, -0.5, 63.5);
       dttf_eta[iwh][ise]->setAxisTitle("#eta", 1);
     }
 
     dttf_trk_folder_eta = dttf_trk_folder_wheel + "/EtaFineFraction";
     ibooker.setCurrentFolder(dttf_trk_folder_eta);
 
     for (int ise = 0; ise < 12; ++ise) {
       sprintf(hname, "dttf_etaFine_fraction_wh%s_se%d", c_whn[iwh], ise + 1);
       sprintf(htitle, "Wheel %s Sector %d - Eta Fine Fraction", c_whn[iwh], ise + 1);
       dttf_eta_fine_fraction[iwh][ise] = ibooker.book1D(hname, htitle, 2, 0, 2);
       dttf_eta_fine_fraction[iwh][ise]->setAxisTitle("#eta", 1);
       dttf_eta_fine_fraction[iwh][ise]->setBinLabel(1, "fine", 1);
       dttf_eta_fine_fraction[iwh][ise]->setBinLabel(2, "coarse", 1);
     }
   }
 
   std::string dttf_trk_folder_inclusive = dttf_trk_folder + "/01-INCLUSIVE";
   ibooker.setCurrentFolder(dttf_trk_folder_inclusive);
 
   sprintf(hname, "dttf_01_nTracksPerEvent_integ");
   sprintf(htitle, "Number of DTTF Tracks Per Event");
   dttf_nTracksPerEvent_integ = ibooker.book1D(hname, htitle, 20, 0.5, 20.5);
   dttf_nTracksPerEvent_integ->setAxisTitle("# tracks/event", 1);
 
   // sprintf(hname, "dttf_10_qual_eta_distr");
   // sprintf(htitle, "DTTF Tracks Quality vs Eta Distribution");
   // dttf_qual_eta_integ = dbe_->book2D(hname, htitle, 64, 0, 64, 7, 1, 8);
   // setQualLabel( dttf_qual_eta_integ, 2);
 
   if (online_) {
     sprintf(hname, "dttf_04_tracks_occupancy_by_lumi");
     sprintf(htitle, "DTTF Tracks in the last LumiSections");
     dttf_spare = ibooker.book2D(hname, htitle, 6, 0, 6, 12, 1, 13);
     setWheelLabel(dttf_spare);
     dttf_spare->setAxisTitle("Sector", 2);
     dttf_spare->getTH2F()->GetXaxis()->SetNdivisions(12);
   } else {
     sprintf(hname, "dttf_04_global_muons_request");
     sprintf(htitle, "Tracks compatible with a Global Muon in the Barrel");
     dttf_spare = ibooker.book1D(hname, htitle, 4, -0.5, 3.5);
     dttf_spare->setBinLabel(1, "No tracks", 1);
     dttf_spare->setBinLabel(2, "No tracks but GM", 1);
     dttf_spare->setBinLabel(3, "Tracks wo GM", 1);
     dttf_spare->setBinLabel(4, "Tracks w GM", 1);
   }
 
   std::string dttf_trk_folder_integrated_gmt = dttf_trk_folder + "/08-GMT_MATCH";
   ibooker.setCurrentFolder(dttf_trk_folder_integrated_gmt);
 
   sprintf(hname, "dttf_tracks_with_gmt_match");
   sprintf(htitle, "DTTF Tracks With a Match in GMT");
   dttf_gmt_match = ibooker.book2D(hname, htitle, 6, 0., 6., 12, 1., 13.);
   setWheelLabel(dttf_gmt_match);
 
   sprintf(hname, "dttf_tracks_without_gmt_match");
   sprintf(htitle, "DTTF Tracks Without a Match in GMT");
   dttf_gmt_missed = ibooker.book2D(hname, htitle, 6, 0., 6., 12, 1., 13.);
   setWheelLabel(dttf_gmt_missed);
 
   sprintf(hname, "dttf_missing_tracks_in_gmt");
   sprintf(htitle, "GMT Tracks Without a Corresponding Track in DTTF");
   dttf_gmt_ghost = ibooker.book2D(hname, htitle, 5, -2, 3, 12, 1, 13.);
 
   dttf_gmt_ghost->setBinLabel(1, "N2", 1);
   dttf_gmt_ghost->setBinLabel(2, "N1", 1);
   dttf_gmt_ghost->setBinLabel(3, "N0/P0", 1);
   dttf_gmt_ghost->setBinLabel(4, "P1", 1);
   dttf_gmt_ghost->setBinLabel(5, "P2", 1);
 
   // sprintf(hname, "dttf_eta_phi_missing_tracks_in_gmt");
   // sprintf(htitle, "GMT Tracks Without a Corresponding Track in DTTF");
   // dttf_gmt_ghost_phys = dbe_->book2D(hname, htitle, 64, 0., 64., 144, 0., 144. );
 }

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 232 of file L1TDTTF.cc.

                                                                  {
   //empty
 }

void L1TDTTF::fillMEs	(	std::vector< L1MuDTTrackCand > const *	trackContainer,
		std::vector< L1MuRegionalCand > &	gmtDttfCands
	)

private

Forget N0 with zero eta value for physical values

from 0 to 11

from 1 to 12

wh has possible values {-3,-2,-1,1,2,3}

make wh2 go from 0 to 5

useful conversions

calculate phi in physical coordinates: keep it int, set labels later

new attempt

Fill per sector bx WHEEL_s/dttf_bx_whs

Fill per sector 2nd bx

WHEEL_s/BX_SECTORS/TRACK_2_ONLY/dttf_bx_2ndTrack_whs_sed

COUNTERS global

Fill per sector phi: WHEEL_s/BX_d/dttf_phi_whs_sed

Fill per sector quality WHEEL_s/BX_d/dttf_qual_whs_sed

Fill per sector pt WHEEL_s/BX_d/dttf_pt_whs_sed

Fill per sector charge WHEEL_s/BX_d/dttf_q_whs_sed

Fill per sector eta WHEEL_s/BX_d/dttf_eta_whs_sed

WHEEL_s/dttf_phi_eta_whs

Only for online: INCLUSIVE/dttf_occupancy_summary_r

second track summary

WHEEL_s/dttf_phi_integ

WHEEL_s/dttf_pt_integ

WHEEL_s/dttf_eta_integ

WHEEL_s/dttf_qual_integ

WHEEL_s/dttf_q_integ

WHEEL_s/dttf_quality_whs

WHEEL_s/dttf_phi_eta_whs

gmt phi_packed() goes from 0 to 143

calculate phi in physical coordinates: keep it int, set labels later

Definition at line 697 of file L1TDTTF.cc.

References dttf_bx, dttf_bx_2ndTrack, dttf_eta, dttf_eta_fine_fraction, dttf_eta_wheel_2ndTrack, dttf_gmt_match, dttf_gmt_missed, dttf_phi, dttf_phi_eta_coarse_wheel, dttf_phi_eta_fine_wheel, dttf_phi_eta_wheel_2ndTrack, dttf_phi_wheel_2ndTrack, dttf_pt, dttf_pt_wheel_2ndTrack, dttf_q, dttf_q_wheel_2ndTrack, dttf_qual, dttf_quality_summary_wheel_2ndTrack, dttf_quality_wheel_2ndTrack, dttf_spare, dqm::impl::MonitorElement::Fill(), match(), numTracks, online_, HLT_2018_cff::track, and verbose_.

Referenced by analyze().

                                                                                                                  {
   L1MuDTTrackContainer::TrackContainer::const_iterator track = trackContainer->begin();
   L1MuDTTrackContainer::TrackContainer::const_iterator trackEnd = trackContainer->end();
 
   for (; track != trackEnd; ++track) {
     if (verbose_) {
       edm::LogInfo("L1TDTTF::Analyze") << "bx = " << track->bx();
       edm::LogInfo("L1TDTTF::Analyze") << "quality (packed) = " << track->quality_packed();
       edm::LogInfo("L1TDTTF::Analyze") << "pt      (packed) = " << track->pt_packed()
                                        << "  , pt  (GeV) = " << track->ptValue();
       edm::LogInfo("L1TDTTF::Analyze") << "phi     (packed) = " << track->phi_packed()
                                        << " , phi (rad) = " << track->phiValue();
       edm::LogInfo("L1TDTTF::Analyze") << "charge  (packed) = " << track->charge_packed();
     }
 
     if ((track->whNum() == -1) && !track->eta_packed()) {
       edm::LogInfo("L1TDTTF::Analyze") << "Skipping N0 with zero eta value";
 
       continue;
     }
 
     int bxindex = track->bx() + 1;
     int se = track->scNum();       
     int sector = se + 1;           
     int whindex = track->whNum();  
 
     whindex = (whindex < 0) ? whindex + 3 : whindex + 2;  
 
     if (whindex < 0 || whindex > 5) {
       edm::LogError("L1TDTTF::Analyze::WHEEL_ERROR") << track->whNum() << "(" << whindex << ")";
       continue;
     }
 
     if (se < 0 || se > 11) {
       edm::LogError("L1TDTTF::Analyze::SECTOR_ERROR") << se;
       continue;
     }
 
 
     // int phi_local = track->phi_packed();//range: 0 < phi_local < 31
     // if ( phi_local > 15 ) phi_local -= 32; //range: -16 < phi_local < 15
 
     // int phi_global = phi_local + se * 12; //range: -16 < phi_global < 147
     // if(phi_global < 0) phi_global += 144; //range: 0 < phi_global < 147
     // if(phi_global > 143) phi_global -= 144; //range: 0 < phi_global < 143
     // // float phi_phys = phi_global * 2.5 + 1.25;
 
     int phi_global = track->phi_packed();
     phi_global = (phi_global > 15 ? phi_global - 32 : phi_global) + se * 12;
     if (phi_global < -6)
       phi_global += 144;  //range: 0 < phi_global < 147
     if (phi_global > 137)
       phi_global -= 144;  //range: 0 < phi_global < 143
 
     // int eta_global = track->eta_packed();
     // int eta_global = track->eta_packed() - 32;
     // dttf_eta[bxindex][whindex][se]->Fill(eta_global);
     // float eta_phys = eta_global / 2.4 ;
 
 
     dttf_bx[whindex][se]->Fill(track->bx());
 
     if (track->TrkTag() == 1) {
       dttf_bx_2ndTrack[whindex][se]->Fill(track->bx());
     }
 
 
     if (bxindex == 1) {
       ++numTracks[whindex][se];
 
       dttf_phi[whindex][se]->Fill(phi_global);
 
       dttf_qual[whindex][se]->Fill(track->quality_packed());
 
       dttf_pt[whindex][se]->Fill(track->pt_packed());
 
       dttf_q[whindex][se]->Fill(track->charge_packed());
 
       dttf_eta[whindex][se]->Fill(track->eta_packed());
 
       if (track->isFineHalo()) {
         dttf_eta_fine_fraction[whindex][se]->Fill(0);
 
         dttf_phi_eta_fine_wheel[whindex]->Fill(track->eta_packed(), phi_global);
 
       } else {
         dttf_eta_fine_fraction[whindex][se]->Fill(1);
 
         dttf_phi_eta_coarse_wheel[whindex]->Fill(track->eta_packed(), phi_global);
       }
 
       if (online_) {
         dttf_spare->Fill(whindex, sector);
       }
 
       // dttf_qual_eta_integ->Fill(track->eta_packed(), track->quality_packed());
 
       if (track->TrkTag() == 1) {
         dttf_phi_wheel_2ndTrack[whindex]->Fill(phi_global);
 
         dttf_pt_wheel_2ndTrack[whindex]->Fill(track->pt_packed());
 
         dttf_eta_wheel_2ndTrack[whindex]->Fill(track->eta_packed());
 
         dttf_quality_wheel_2ndTrack[whindex]->Fill(track->quality_packed());
 
         dttf_q_wheel_2ndTrack[whindex]->Fill(track->charge_packed());
 
         dttf_quality_summary_wheel_2ndTrack[whindex]->Fill(sector, track->quality_packed());
 
         dttf_phi_eta_wheel_2ndTrack[whindex]->Fill(track->eta_packed(), phi_global);
       }
 
       bool match = false;
       std::vector<L1MuRegionalCand>::iterator dttfCand;
       unsigned int gmt_phi = (phi_global < 0 ? phi_global + 144 : phi_global);
 
       for (dttfCand = gmtDttfCands.begin(); dttfCand != gmtDttfCands.end(); ++dttfCand) {
         if (dttfCand->empty())
           continue;
         if ((dttfCand->phi_packed() == gmt_phi) && dttfCand->quality_packed() == track->quality_packed()) {
           match = true;
           dttfCand->reset();
           break;
         }
       }
 
       if (match) {
         dttf_gmt_match->Fill(whindex, sector);
       } else {
         dttf_gmt_missed->Fill(whindex, sector);
       }
     }
   }
 }

void L1TDTTF::setQualLabel	(	MonitorElement *	me,
		int	axis
	)

private

Definition at line 868 of file L1TDTTF.cc.

References dqm::impl::MonitorElement::setAxisTitle(), and dqm::impl::MonitorElement::setBinLabel().

Referenced by bookHistograms().

                                                        {
   if (axis == 1)
     me->setAxisTitle("Quality", axis);
   me->setBinLabel(1, "T34", axis);
   me->setBinLabel(2, "T23/24", axis);
   me->setBinLabel(3, "T12/13/14", axis);
   me->setBinLabel(4, "T234", axis);
   me->setBinLabel(5, "T134", axis);
   me->setBinLabel(6, "T123/124", axis);
   me->setBinLabel(7, "T1234", axis);
 }

void L1TDTTF::setWheelLabel ( MonitorElement * me )

private

Definition at line 881 of file L1TDTTF.cc.

References dqm::impl::MonitorElement::setAxisTitle(), and dqm::impl::MonitorElement::setBinLabel().

Referenced by bookHistograms().

                                               {
   me->setAxisTitle("Wheel", 1);
   me->setBinLabel(1, "N2", 1);
   me->setBinLabel(2, "N1", 1);
   me->setBinLabel(3, "N0", 1);
   me->setBinLabel(4, "P0", 1);
   me->setBinLabel(5, "P1", 1);
   me->setBinLabel(6, "P2", 1);
 }