#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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

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_

Detailed Description

Definition at line 31 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 197 of file L1TDTTF.cc.

References dttpgSource_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), gmtSourceToken_, 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 239 of file L1TDTTF.cc.

240 {

242 }

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 581 of file L1TDTTF.cc.

References accept(), dttf_gmt_ghost, dttf_nTracksPerEv, dttf_nTracksPerEvent_integ, dttf_nTracksPerEvent_wheel, dttf_spare, dttpgSource_, PVValHelper::eta, MonitorElement::Fill(), fillMEs(), L1MuDTTrackContainer::getContainer(), L1MuGMTReadoutCollection::getRecords(), gmtSource_, gmtSourceToken_, edm::HandleBase::isValid(), edm::InputTag::label(), muonCollectionLabel_, muonCollectionToken_, extraflags_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 1015 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 &
	)

overrideprotected

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

Definition at line 251 of file L1TDTTF.cc.

References DQMStore::IBooker::book1D(), 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, MonitorElement::getTH2F(), l1tsubsystemfolder_, pileupCalc::nbins, nev_, nev_dttf_, nev_dttf_track2_, online_, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), 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 &
	)

overrideprotected

Definition at line 245 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 787 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, MonitorElement::Fill(), match(), numTracks, online_, HiIsolationCommonParameters_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 986 of file L1TDTTF.cc.

References MonitorElement::setAxisTitle(), and 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 1001 of file L1TDTTF.cc.

References MonitorElement::setAxisTitle(), and 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);
 }