#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

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

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

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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 DQMEDAnalyzer
uint64_t	meId () const

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::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

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_FULL_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
	)

overrideprotectedvirtual

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

Reimplemented from DQMEDAnalyzer.

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(), L1MuGMTReadoutCollection::getRecords(), gmtSource_, gmtSourceToken_, edm::HandleBase::isValid(), edm::InputTag::label(), muonCollectionLabel_, muonCollectionToken_, patZpeak::muons, nev_, nev_dttf_, nev_dttf_track2_, numTracks, online_, phi, edm::InputTag::process(), edm::Handle< T >::product(), alignCSCRings::s, trackInputToken_, verbose_, and w.

                                                                             {
   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::implementation::IBooker::book1D(), dqm::implementation::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_, hlt_dqm_clientPB-live_cfg::nbins, nev_, nev_dttf_, nev_dttf_track2_, online_, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), setQualLabel(), setWheelLabel(), command_line::start, 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 &	r,
		edm::EventSetup const &	c
	)

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_FULL_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);
 }