CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes
L1TDTTF Class Reference

#include <L1TDTTF.h>

Inheritance diagram for L1TDTTF:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 L1TDTTF (const edm::ParameterSet &ps)
 base services More...
 
virtual ~L1TDTTF ()
 
- Public Member Functions inherited from edm::EDAnalyzer
 EDAnalyzer ()
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndex indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Protected Member Functions

void analyze (const edm::Event &e, const edm::EventSetup &c)
 
void beginJob (void)
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
 
void endJob (void)
 
void endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
 
- Protected Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
CurrentProcessingContext const * currentContext () const
 
- 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 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 > *trackContainer, std::vector< L1MuRegionalCand > &gmtDttfCands)
 
void setQualLabel (MonitorElement *me, int axis)
 
void setWheelLabel (MonitorElement *me)
 

Private Attributes

DQMStoredbe_
 
MonitorElementdttf_bx [6][12]
 
MonitorElementdttf_bx_2ndTrack [6][12]
 
MonitorElementdttf_eta [6][12]
 
MonitorElementdttf_eta_fine_fraction [6][12]
 
MonitorElementdttf_eta_wheel_2ndTrack [6]
 
MonitorElementdttf_gmt_ghost
 
MonitorElementdttf_gmt_match
 
MonitorElementdttf_gmt_missed
 
MonitorElementdttf_nTracksPerEv [6][12]
 
MonitorElementdttf_nTracksPerEvent_integ
 
MonitorElementdttf_nTracksPerEvent_wheel [6]
 
MonitorElementdttf_phi [6][12]
 
MonitorElementdttf_phi_eta_coarse_wheel [6]
 
MonitorElementdttf_phi_eta_fine_wheel [6]
 
MonitorElementdttf_phi_eta_wheel_2ndTrack [6]
 
MonitorElementdttf_phi_wheel_2ndTrack [6]
 
MonitorElementdttf_pt [6][12]
 
MonitorElementdttf_pt_wheel_2ndTrack [6]
 
MonitorElementdttf_q [6][12]
 
MonitorElementdttf_q_wheel_2ndTrack [6]
 
MonitorElementdttf_qual [6][12]
 
MonitorElementdttf_quality_summary_wheel_2ndTrack [6]
 
MonitorElementdttf_quality_wheel_2ndTrack [6]
 
MonitorElementdttf_spare
 
edm::InputTag dttpgSource_
 
edm::InputTag gmtSource_
 
std::string l1tsubsystemfolder_
 
edm::InputTag muonCollectionLabel_
 
int nev_
 
int nev_dttf_
 
int nev_dttf_track2_
 
int numTracks [6][12]
 
bool online_
 
std::string outputFile_
 
edm::InputTag trackInputTag_
 
bool verbose_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
typedef WorkerT< EDAnalyzerWorkerType
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 

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 DQMStore?

Use ROOT Output?

Definition at line 211 of file L1TDTTF.cc.

References dbe_, dttpgSource_, edm::ParameterSet::getUntrackedParameter(), l1tsubsystemfolder_, edm::InputTag::label(), NULL, cppFunctionSkipper::operator, outputFile_, edm::InputTag::process(), DQMStore::setCurrentFolder(), DQMStore::setVerbose(), AlCaHLTBitMon_QueryRunRegistry::string, trackInputTag_, and verbose_.

212  : dttpgSource_( ps.getParameter< edm::InputTag >("dttpgSource") ),
213  gmtSource_( ps.getParameter< edm::InputTag >("gmtSource") ),
214  muonCollectionLabel_( ps.getParameter<edm::InputTag>("MuonCollection") ),
216  "L1T/L1TDTTF")),
217  online_( ps.getUntrackedParameter<bool>("online", true) ),
218  verbose_( ps.getUntrackedParameter<bool>("verbose", false) )
219 
220 {
221 
222  std::string trstring =
223  dttpgSource_.label() + ":DATA:" + dttpgSource_.process();
224  trackInputTag_ = edm::InputTag(trstring);
225 
227  if ( verbose_ ) edm::LogInfo("L1TDTTF: constructor") << "Verbose enabled";
228 
230  dbe_ = NULL;
231  if ( ps.getUntrackedParameter<bool>("DQMStore", false) ) {
233  dbe_->setVerbose(0);
235  }
236 
238  if ( ps.getUntrackedParameter<bool>("disableROOToutput", false) ) {
239 
240  outputFile_ = "";
241 
242  } else {
243 
244  outputFile_ = ps.getUntrackedParameter<std::string>("outputFile", "");
245  if ( ! outputFile_.empty() ) {
246  edm::LogInfo("L1TDTTF: constructor")
247  << "L1T Monitoring histograms will be saved to " << outputFile_;
248  }
249 
250  }
251 
252 }
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
std::string l1tsubsystemfolder_
Definition: L1TDTTF.h:70
edm::InputTag trackInputTag_
Definition: L1TDTTF.h:75
edm::InputTag muonCollectionLabel_
Definition: L1TDTTF.h:69
#define NULL
Definition: scimark2.h:8
DQMStore * dbe_
Definition: L1TDTTF.h:73
edm::InputTag gmtSource_
Definition: L1TDTTF.h:68
std::string outputFile_
Definition: L1TDTTF.h:74
edm::InputTag dttpgSource_
Definition: L1TDTTF.h:67
void setVerbose(unsigned level)
Definition: DQMStore.cc:398
bool verbose_
Definition: L1TDTTF.h:72
std::string const & label() const
Definition: InputTag.h:42
std::string const & process() const
Definition: InputTag.h:46
bool online_
Definition: L1TDTTF.h:71
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
L1TDTTF::~L1TDTTF ( )
virtual

Nothing to destroy

Definition at line 257 of file L1TDTTF.cc.

258 {
260 }

Member Function Documentation

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

counters

tracks handle

dttf counters

selection for offline

global muon selection plot

in case of problems accept all

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

Implements edm::EDAnalyzer.

Definition at line 621 of file L1TDTTF.cc.

References accept(), dttf_gmt_ghost, dttf_nTracksPerEv, dttf_nTracksPerEvent_integ, dttf_nTracksPerEvent_wheel, dttf_spare, dttpgSource_, eta(), MonitorElement::Fill(), fillMEs(), L1MuGMTReadoutCollection::getRecords(), gmtSource_, edm::HandleBase::isValid(), edm::InputTag::label(), muonCollectionLabel_, patZpeak::muons, nev_, nev_dttf_, nev_dttf_track2_, numTracks, online_, phi, edm::InputTag::process(), edm::Handle< T >::product(), alignCSCRings::s, trackInputTag_, verbose_, and w().

623 {
624 
625 
626  if ( verbose_ )
627  edm::LogInfo("L1TDTTF::Analyze::start") << "#################### START";
628 
630  ++nev_;
631  memset( numTracks, 0, 72 * sizeof(int) );
632 
634  edm::Handle<L1MuDTTrackContainer > myL1MuDTTrackContainer;
635  try {
636  event.getByLabel(trackInputTag_, myL1MuDTTrackContainer);
637  } catch (cms::Exception& iException) {
638  edm::LogError("L1TDTTF::analyze::DataNotFound")
639  << "can't getByLabel L1MuDTTrackContainer with label "
640  << dttpgSource_.label() << ":DATA:" << dttpgSource_.process();
641  return;
642  }
643 
644  if ( !myL1MuDTTrackContainer.isValid() ) {
645  edm::LogError("L1TDTTF::analyze::DataNotFound")
646  << "can't find L1MuDTTrackContainer with label "
647  << dttpgSource_.label() << ":DATA:" << dttpgSource_.process();
648  return;
649  }
650 
651  L1MuDTTrackContainer::TrackContainer * trackContainer =
652  myL1MuDTTrackContainer->getContainer();
653 
655  if ( trackContainer->size() > 0 ) {
656  ++nev_dttf_;
657  if( trackContainer->size() > 1 ) ++nev_dttf_track2_;
658  }
659 
663  bool accept = true;
664  if ( ! online_ ) {
665 
666  try {
667 
669  event.getByLabel(muonCollectionLabel_, muons);
670 
671  accept = false;
672  if ( muons.isValid() ) {
673  for (reco::MuonCollection::const_iterator recoMu = muons->begin();
674  recoMu!=muons->end(); ++recoMu ) {
675  if ( fabs( recoMu->eta() ) < 1.4 ) {
676  if ( verbose_ ) {
677  edm::LogInfo("L1TDTTFClient::Analyze:GM") << "Found a global muon!";
678  }
679  accept = true;
680  break;
681  }
682 
683  }
684 
686  if ( ! accept ) {
687  dttf_spare->Fill( trackContainer->size() ? 1 : 0 );
688 
689  if ( verbose_ ) {
690  edm::LogInfo("L1TDTTFClient::Analyze:GM")
691  << "No global muons in this event!";
692  }
693 
694  } else {
695  dttf_spare->Fill( trackContainer->size() ? 2 : 3 );
696  }
697 
698  } else {
700  accept = true;
701  edm::LogWarning("L1TDTTFClient::Analyze:GM")
702  << "Invalid MuonCollection with label "
704  }
705 
706 
707  } catch (cms::Exception& iException) {
709  accept = true;
710  edm::LogError("DataNotFound") << "Unable to getByLabel MuonCollection with label "
712  }
713 
714  }
715 
716 
718  std::vector<L1MuRegionalCand> gmtBx0DttfCandidates;
719 
720  try {
721 
723  event.getByLabel(gmtSource_, pCollection);
724 
725  if ( !pCollection.isValid() ) {
726  edm::LogError("DataNotFound") << "can't find L1MuGMTReadoutCollection with label "
727  << gmtSource_.label() ;
728  }
729 
730  // get GMT readout collection
731  L1MuGMTReadoutCollection const* gmtrc = pCollection.product();
732  std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
733 
734  std::vector<L1MuGMTReadoutRecord>::const_iterator RRItr;
735 
736  for ( RRItr = gmt_records.begin(); RRItr != gmt_records.end(); ++RRItr ) {
737 
738  std::vector<L1MuRegionalCand> dttfCands = RRItr->getDTBXCands();
739  std::vector<L1MuRegionalCand>::iterator dttfCand;
740 
741  for( dttfCand = dttfCands.begin(); dttfCand != dttfCands.end();
742  ++dttfCand ) {
743 
744  if(dttfCand->empty()) continue;
746  if ( RRItr->getBxInEvent() ) continue;
747 
748  // dttf_gmt_ghost_phys->Fill( dttfCand->eta_packed(),
749  // dttfCand->phi_packed() );
750  gmtBx0DttfCandidates.push_back( *dttfCand );
751 
752  }
753  }
754 
755  } catch (cms::Exception& iException) {
756  edm::LogError("DataNotFound") << "Unable to getByLabel L1MuGMTReadoutCollection with label "
757  << gmtSource_.label() ;
758  }
759 
760 
761  // fill MEs if all selections are passed
762  if ( accept ) fillMEs( trackContainer, gmtBx0DttfCandidates );
763 
765  std::vector<L1MuRegionalCand>::iterator dttfCand;
766  for( dttfCand = gmtBx0DttfCandidates.begin();
767  dttfCand != gmtBx0DttfCandidates.end(); ++dttfCand ) {
768  if( dttfCand->empty() ) continue;
769 
773  int phi= dttfCand->phi_packed();
774  int sector = 1 + (phi + 6)/12;
775  if (sector > 12 ) sector -= 12;
776  double eta = dttfCand->etaValue();
777 
778  int wheel = -3;
779  if ( eta < -0.74 ) {
780  wheel = -2;
781  } else if ( eta < -0.3 ) {
782  wheel = -1;
783 
784  } else if ( eta < 0.3 ) {
785  wheel = 0;
786 
787  } else if ( eta < 0.74 ) {
788  wheel = 1;
789  } else {
790  wheel = 2;
791  }
792 
793  dttf_gmt_ghost->Fill( wheel, sector );
794  // dttf_gmt_ghost_phys->Fill( dttfCand->eta_packed(),
795  // dttfCand->phi_packed() );
796  }
797 
798 
800  int numTracksInt = 0;
801 
802  for ( int w = 0; w < 6; ++w ) {
803 
804  int numTracks_wh = 0;
805  for ( int s = 0; s < 12; ++s ) {
806 
808 
809  numTracks_wh += numTracks[w][s];
810 
811  }
812 
813  numTracksInt += numTracks_wh;
814  dttf_nTracksPerEvent_wheel[w]->Fill( numTracks_wh );
815 
816  }
817 
818  dttf_nTracksPerEvent_integ->Fill( numTracksInt );
819 
820 
821 
822 }
edm::InputTag trackInputTag_
Definition: L1TDTTF.h:75
edm::InputTag muonCollectionLabel_
Definition: L1TDTTF.h:69
MonitorElement * dttf_nTracksPerEvent_wheel[6]
Definition: L1TDTTF.h:77
T eta() const
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:22
void Fill(long long x)
void fillMEs(std::vector< L1MuDTTrackCand > *trackContainer, std::vector< L1MuRegionalCand > &gmtDttfCands)
Definition: L1TDTTF.cc:828
MonitorElement * dttf_nTracksPerEvent_integ
Definition: L1TDTTF.h:98
edm::InputTag gmtSource_
Definition: L1TDTTF.h:68
edm::InputTag dttpgSource_
Definition: L1TDTTF.h:67
std::vector< L1MuDTTrackCand > TrackContainer
bool isValid() const
Definition: HandleBase.h:76
int nev_dttf_track2_
Definition: L1TDTTF.h:108
int numTracks[6][12]
Definition: L1TDTTF.h:109
bool verbose_
Definition: L1TDTTF.h:72
MonitorElement * dttf_gmt_ghost
Definition: L1TDTTF.h:103
int nev_dttf_
Definition: L1TDTTF.h:107
T const * product() const
Definition: Handle.h:74
std::string const & label() const
Definition: InputTag.h:42
std::string const & process() const
Definition: InputTag.h:46
tuple muons
Definition: patZpeak.py:38
bool online_
Definition: L1TDTTF.h:71
std::vector< L1MuGMTReadoutRecord > const & getRecords() const
MonitorElement * dttf_nTracksPerEv[6][12]
Definition: L1TDTTF.h:88
MonitorElement * dttf_spare
Definition: L1TDTTF.h:99
T w() const
int nev_
Definition: L1TDTTF.h:106
Definition: DDAxes.h:10
void L1TDTTF::beginJob ( void  )
protectedvirtual

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

ETA folder

integrated values: always packed

Only for online: occupancy summary - reset

Reimplemented from edm::EDAnalyzer.

Definition at line 265 of file L1TDTTF.cc.

References DQMStore::book1D(), DQMStore::book2D(), bookEta(), dbe_, 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::setCurrentFolder(), setQualLabel(), setWheelLabel(), dqm_diff::start, and AlCaHLTBitMon_QueryRunRegistry::string.

266 {
268  nev_ = 0;
269  nev_dttf_ = 0;
270  nev_dttf_track2_ = 0;
271 
272  // get hold of back-end interface
273 
274  if ( dbe_ ) {
275 
276  std::string dttf_trk_folder = l1tsubsystemfolder_;
277 
278  char hname[100];
279  char htitle[100];
280 
282  float start = 0;
283  float stop = 0;
284  int nbins = 0;
286 
288  dbe_->setCurrentFolder(dttf_trk_folder);
289 
290  std::string wheelpath[6] = { "/02-WHEEL_N2",
291  "/03-WHEEL_N1",
292  "/04-WHEEL_N0",
293  "/05-WHEEL_P0",
294  "/06-WHEEL_P1",
295  "/07-WHEEL_P2" };
296 
297 
298  char c_whn[6][3] = { "N2", "N1", "N0", "P0", "P1", "P2" };
299  // char bxn [3][3] = { "N1", "0", "P1" };
300  // char bxn[3][25] = {"/BX_NONZERO_ONLY/BX_N1", "", "/BX_NONZERO_ONLY/BX_P1"};
301 
302  for ( int iwh = 0; iwh < 6; ++iwh ) {
303 
304  bookEta( iwh, nbins, start, stop );
305 
309  std::string dttf_trk_folder_wheel = dttf_trk_folder + wheelpath[iwh];
310  dbe_->setCurrentFolder(dttf_trk_folder_wheel);
311 
313  sprintf(hname, "dttf_01_nTracksPerEvent_wh%s", c_whn[iwh]);
314  sprintf(htitle, "Wheel %s - Number Tracks Per Event", c_whn[iwh]);
315  dttf_nTracksPerEvent_wheel[iwh] = dbe_->book1D(hname, htitle,
316  10, 0.5, 10.5);
317  dttf_nTracksPerEvent_wheel[iwh]->setAxisTitle("# tracks/event", 1);
318 
320  sprintf(hname, "dttf_07_phi_vs_etaFine_wh%s", c_whn[iwh]);
321  sprintf(htitle, "Wheel %s - #eta-#phi DTTF Tracks occupancy (fine #eta only, unpacked values)", c_whn[iwh]);
322  dttf_phi_eta_fine_wheel[iwh] = dbe_->book2D(hname, htitle,
323  nbins, start-0.5, stop-0.5,
324  144, -6, 138);
325  // 144, -0.5, 143.5);
326 
327  dttf_phi_eta_fine_wheel[iwh]->setAxisTitle("#eta", 1);
328  dttf_phi_eta_fine_wheel[iwh]->setAxisTitle("#phi", 2);
329 
331  sprintf(hname, "dttf_08_phi_vs_etaCoarse_wh%s", c_whn[iwh]);
332  sprintf(htitle, "Wheel %s - #eta-#phi DTTF Tracks occupancy (coarse #eta only, unpacked values)", c_whn[iwh]);
333  dttf_phi_eta_coarse_wheel[iwh] = dbe_->book2D(hname, htitle,
334  nbins, start-0.5, stop-0.5,
335  144, -6, 138);
336  // 144, -0.5, 143.5);
337  dttf_phi_eta_coarse_wheel[iwh]->setAxisTitle("#eta", 1);
338  dttf_phi_eta_coarse_wheel[iwh]->setAxisTitle("#phi", 2);
339 
342  std::string dttf_trk_folder_wheel_2ndtrack =
343  dttf_trk_folder_wheel + "/2ND_TRACK_ONLY";
344  dbe_->setCurrentFolder(dttf_trk_folder_wheel_2ndtrack);
345 
346 
348  sprintf(hname, "dttf_04_quality_wh%s_2ndTrack", c_whn[iwh]);
349  sprintf(htitle, "Wheel %s - 2nd Tracks Quality distribution", c_whn[iwh]);
350  dttf_quality_wheel_2ndTrack[iwh] = dbe_->book1D(hname, htitle, 7, 1, 8);
352 
354  sprintf(hname, "dttf_05_quality_summary_wh%s_2ndTrack", c_whn[iwh]);
355  sprintf(htitle, "Wheel %s - 2nd Tracks - Quality", c_whn[iwh]);
356  dttf_quality_summary_wheel_2ndTrack[iwh] = dbe_->book2D(hname, htitle,
357  12, 1, 13, 7, 1, 8 );
360  // dttf_quality_summary_wheel_2ndTrack[iwh]->setAxisTitle("Quality", 2);
361 
363  sprintf(hname, "dttf_06_phi_vs_eta_wh%s_2ndTrack", c_whn[iwh]);
364  sprintf(htitle, "Wheel %s - #eta-#phi Distribution of DTTF 2nd Tracks",
365  c_whn[iwh]);
366 
367  dttf_phi_eta_wheel_2ndTrack[iwh] = dbe_->book2D(hname, htitle,
368  nbins, start-0.5,stop-0.5,
369  144, -6, 138);
370  // 144, -0.5, 143.5);
371  dttf_phi_eta_wheel_2ndTrack[iwh]->setAxisTitle("#eta", 1);
372  dttf_phi_eta_wheel_2ndTrack[iwh]->setAxisTitle("#phi", 2);
373 
374 
375 
377  sprintf(hname, "dttf_07_eta_wh%s_2ndTrack", c_whn[iwh]);
378  sprintf(htitle, "Wheel %s - DTTF 2nd Tracks #eta distribution (Packed values)",
379  c_whn[iwh]);
380  dttf_eta_wheel_2ndTrack[iwh] = dbe_->book1D(hname, htitle, 64, -0.5, 63.5);
381  dttf_eta_wheel_2ndTrack[iwh]->setAxisTitle("#eta", 1);
382 
384  sprintf(hname, "dttf_08_phi_wh%s_2ndTrack", c_whn[iwh]);
385  sprintf(htitle, "Wheel %s - DTTF 2nd Tracks Phi distribution (Packed values)",
386  c_whn[iwh]);
387  dttf_phi_wheel_2ndTrack[iwh] = dbe_->book1D(hname, htitle, 144, -6, 138. );
388  dttf_phi_wheel_2ndTrack[iwh]->setAxisTitle("#phi", 1);
389 
391  sprintf(hname, "dttf_09_pt_wh%s_2ndTrack", c_whn[iwh]);
392  sprintf(htitle, "Wheel %s - DTTF 2nd Tracks p_{T} distribution (Packed values)",
393  c_whn[iwh]);
394  dttf_pt_wheel_2ndTrack[iwh] = dbe_->book1D(hname, htitle, 32, -0.5, 31.5);
395  dttf_pt_wheel_2ndTrack[iwh]->setAxisTitle("p_{T}", 1);
396 
398  sprintf(hname, "dttf_10_charge_wh%s_2ndTrack", c_whn[iwh]);
399  sprintf(htitle, "Wheel %s - DTTF 2nd Tracks Charge distribution", c_whn[iwh]);
400  dttf_q_wheel_2ndTrack[iwh] = dbe_->book1D(hname, htitle, 2, -0.5, 1.5);
401  dttf_q_wheel_2ndTrack[iwh]->setAxisTitle("Charge", 1);
402 
403 
404 
405 
409 
411  std::string dttf_trk_folder_nTracksPerEvent = dttf_trk_folder_wheel + "/TracksPerEvent";
412  dbe_->setCurrentFolder(dttf_trk_folder_nTracksPerEvent);
413 
414  for(int ise = 0; ise < 12; ++ise) {
415  sprintf(hname, "dttf_nTracksPerEvent_wh%s_se%d", c_whn[iwh], ise+1);
416  sprintf(htitle, "Wheel %s Sector %d - Number of Tracks Per Event",
417  c_whn[iwh], ise+1);
418  dttf_nTracksPerEv[iwh][ise] = dbe_->book1D(hname, htitle, 2, 0.5, 2.5);
419  dttf_nTracksPerEv[iwh][ise]->setAxisTitle("# tracks/event", 1);
420  }
421 
422 
424  std::string dttf_trk_folder_wh_bxsec_all =
425  dttf_trk_folder_wheel + "/BX_BySector";
426  dbe_->setCurrentFolder(dttf_trk_folder_wh_bxsec_all);
427 
428  for(int ise = 0; ise < 12; ++ise ) {
429  sprintf(hname, "dttf_bx_wh%s_se%d", c_whn[iwh], ise+1);
430  sprintf(htitle, "Wheel %s Sector %d - BX Distribution",
431  c_whn[iwh], ise+1);
432  dttf_bx[iwh][ise] = dbe_->book1D(hname, htitle, 3, -1.5, 1.5);
433  dttf_bx[iwh][ise]->setAxisTitle("BX", 1);
434  }
435 
436  std::string dttf_trk_folder_wh_bxsec_trk2 =
437  dttf_trk_folder_wheel + "/BX_BySector/2ND_TRACK_ONLY";
438  dbe_->setCurrentFolder(dttf_trk_folder_wh_bxsec_trk2);
439 
440  for(int ise = 0; ise < 12; ++ise ) {
441  sprintf(hname, "dttf_bx_2ndTrack_wh%s_se%d", c_whn[iwh], ise+1);
442  sprintf(htitle, "Wheel %s Sector %d - BX 2nd Tracks only",
443  c_whn[iwh], ise+1);
444  dttf_bx_2ndTrack[iwh][ise] = dbe_->book1D(hname, htitle, 3, -1.5, 1.5);
445  dttf_bx_2ndTrack[iwh][ise]->setAxisTitle("BX", 1);
446  }
447 
449  std::string dttf_trk_folder_charge = dttf_trk_folder_wheel + "/Charge";
450  dbe_->setCurrentFolder(dttf_trk_folder_charge);
451 
452  for(int ise = 0; ise < 12; ++ise) {
453  sprintf(hname, "dttf_charge_wh%s_se%d", c_whn[iwh], ise+1);
454  sprintf(htitle, "Wheel %s Sector %d - Packed Charge", c_whn[iwh], ise+1);
455  dttf_q[iwh][ise] = dbe_->book1D(hname, htitle, 2, -0.5, 1.5);
456  dttf_q[iwh][ise]->setAxisTitle("Charge", 1);
457  }
458 
460  std::string dttf_trk_folder_pt = dttf_trk_folder_wheel + "/PT";
461  dbe_->setCurrentFolder(dttf_trk_folder_pt);
462 
463  for(int ise = 0; ise < 12; ++ise ) {
464  sprintf(hname, "dttf_pt_wh%s_se%d", c_whn[iwh], ise+1);
465  sprintf(htitle, "Wheel %s Sector %d - Packed p_{T}",
466  c_whn[iwh], ise + 1 );
467  dttf_pt[iwh][ise]= dbe_->book1D(hname, htitle, 32, -0.5, 31.5);
468  dttf_pt[iwh][ise]->setAxisTitle("p_{T}", 1);
469  }
470 
472  std::string dttf_trk_folder_phi = dttf_trk_folder_wheel + "/Phi";
473  dbe_->setCurrentFolder(dttf_trk_folder_phi);
474 
475  for(int ise = 0; ise < 12; ++ise ) {
476  sprintf(hname, "dttf_phi_wh%s_se%d", c_whn[iwh], ise+1);
477  sprintf(htitle, "Wheel %s Sector %d - Packed Phi", c_whn[iwh], ise+1);
478  dttf_phi[iwh][ise] = dbe_->book1D(hname, htitle, 144, -6, 138);
479  dttf_phi[iwh][ise]->setAxisTitle("#phi", 1);
480  //dttf_phi[iwh][ise] = dbe_->book1D(title,title, 32,-16.5, 15.5);
481  }
482 
484  std::string dttf_trk_folder_quality = dttf_trk_folder_wheel + "/Quality";
485  dbe_->setCurrentFolder(dttf_trk_folder_quality);
486 
487  for(int ise = 0; ise < 12; ++ise){
488  sprintf(hname, "dttf_qual_wh%s_se%d", c_whn[iwh], ise+1);
489  sprintf(htitle, "Wheel %s Sector %d - Packed Quality",
490  c_whn[iwh], ise+1);
491  dttf_qual[iwh][ise] = dbe_->book1D(hname, htitle, 7, 1, 8);
492  dttf_qual[iwh][ise]->setAxisTitle("Quality", 1);
493  setQualLabel( dttf_qual[iwh][ise], 1 );
494  }
495 
497  std::string dttf_trk_folder_eta = dttf_trk_folder_wheel + "/Eta";
498  dbe_->setCurrentFolder(dttf_trk_folder_eta);
499 
500  for (int ise = 0; ise < 12; ++ise ) {
501 
502  sprintf(hname, "dttf_eta_wh%s_se%d", c_whn[iwh], ise+1);
503  sprintf(htitle, "Wheel %s Sector %d - Packed #eta",
504  c_whn[iwh], ise+1);
505  dttf_eta[iwh][ise] = dbe_->book1D(hname, htitle, 64, -0.5, 63.5);
506  dttf_eta[iwh][ise]->setAxisTitle("#eta", 1);
507 
508  }
509 
511  dttf_trk_folder_eta = dttf_trk_folder_wheel + "/EtaFineFraction";
512  dbe_->setCurrentFolder(dttf_trk_folder_eta);
513 
514  for (int ise = 0; ise < 12; ++ise ) {
515 
516  sprintf(hname, "dttf_etaFine_fraction_wh%s_se%d", c_whn[iwh], ise+1);
517  sprintf(htitle, "Wheel %s Sector %d - Eta Fine Fraction",
518  c_whn[iwh], ise+1);
519  dttf_eta_fine_fraction[iwh][ise] = dbe_->book1D(hname, htitle, 2, 0, 2);
520  dttf_eta_fine_fraction[iwh][ise]->setAxisTitle("#eta", 1);
521  dttf_eta_fine_fraction[iwh][ise]->setBinLabel(1, "fine", 1);
522  dttf_eta_fine_fraction[iwh][ise]->setBinLabel(2, "coarse", 1);
523 
524  }
525 
526  }
527 
531  std::string dttf_trk_folder_inclusive = dttf_trk_folder + "/01-INCLUSIVE";
532  dbe_->setCurrentFolder(dttf_trk_folder_inclusive);
533 
534 
535  sprintf(hname, "dttf_01_nTracksPerEvent_integ");
536  sprintf(htitle, "Number of DTTF Tracks Per Event");
537  dttf_nTracksPerEvent_integ = dbe_->book1D(hname, htitle, 20, 0.5, 20.5);
538  dttf_nTracksPerEvent_integ->setAxisTitle("# tracks/event", 1);
539 
541  // sprintf(hname, "dttf_10_qual_eta_distr");
542  // sprintf(htitle, "DTTF Tracks Quality vs Eta Distribution");
543  // dttf_qual_eta_integ = dbe_->book2D(hname, htitle, 64, 0, 64, 7, 1, 8);
544  // setQualLabel( dttf_qual_eta_integ, 2);
545 
547  if ( online_ ) {
548  sprintf(hname, "dttf_04_tracks_occupancy_by_lumi");
549  sprintf(htitle, "DTTF Tracks in the last LumiSections");
550  dttf_spare = dbe_->book2D(hname, htitle, 6, 0, 6, 12, 1, 13);
552  dttf_spare->setAxisTitle("Sector", 2);
553  dttf_spare->getTH2F()->GetXaxis()->SetNdivisions(12);
554  } else {
555 
556  sprintf(hname, "dttf_04_global_muons_request");
557  sprintf(htitle, "Tracks compatible with a Global Muon in the Barrel");
558  dttf_spare = dbe_->book1D(hname, htitle, 4, -0.5, 3.5 );
559  dttf_spare->setBinLabel(1, "No tracks", 1);
560  dttf_spare->setBinLabel(2, "No tracks but GM", 1);
561  dttf_spare->setBinLabel(3, "Tracks wo GM", 1);
562  dttf_spare->setBinLabel(4, "Tracks w GM", 1);
563 
564  }
565 
566  std::string dttf_trk_folder_integrated_gmt =
567  dttf_trk_folder + "/08-GMT_MATCH";
568  dbe_->setCurrentFolder(dttf_trk_folder_integrated_gmt);
569 
570  sprintf(hname, "dttf_tracks_with_gmt_match");
571  sprintf(htitle, "DTTF Tracks With a Match in GMT");
572  dttf_gmt_match = dbe_->book2D(hname, htitle, 6, 0., 6., 12, 1., 13.);
574 
575  sprintf(hname, "dttf_tracks_without_gmt_match");
576  sprintf(htitle, "DTTF Tracks Without a Match in GMT");
577  dttf_gmt_missed = dbe_->book2D(hname, htitle, 6, 0., 6., 12, 1., 13.);
579 
580  sprintf(hname, "dttf_missing_tracks_in_gmt");
581  sprintf(htitle, "GMT Tracks Without a Corresponding Track in DTTF");
582  dttf_gmt_ghost = dbe_->book2D(hname, htitle, 5, -2, 3, 12, 1, 13.);
583 
584  dttf_gmt_ghost->setBinLabel(1, "N2", 1);
585  dttf_gmt_ghost->setBinLabel(2, "N1", 1);
586  dttf_gmt_ghost->setBinLabel(3, "N0/P0", 1);
587  dttf_gmt_ghost->setBinLabel(4, "P1", 1);
588  dttf_gmt_ghost->setBinLabel(5, "P2", 1);
589 
590 
591  // sprintf(hname, "dttf_eta_phi_missing_tracks_in_gmt");
592  // sprintf(htitle, "GMT Tracks Without a Corresponding Track in DTTF");
593  // dttf_gmt_ghost_phys = dbe_->book2D(hname, htitle, 64, 0., 64., 144, 0., 144. );
594 
595 
596  }
597 
598 }
std::string l1tsubsystemfolder_
Definition: L1TDTTF.h:70
tuple start
Check for commandline option errors.
Definition: dqm_diff.py:58
MonitorElement * dttf_phi_eta_coarse_wheel[6]
Definition: L1TDTTF.h:81
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:722
void bookEta(int wh, int &nbins, float &start, float &stop)
Definition: L1TDTTF.cc:1056
MonitorElement * dttf_nTracksPerEvent_wheel[6]
Definition: L1TDTTF.h:77
MonitorElement * dttf_phi_eta_fine_wheel[6]
Definition: L1TDTTF.h:80
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
void setWheelLabel(MonitorElement *me)
Definition: L1TDTTF.cc:1042
MonitorElement * dttf_gmt_match
Definition: L1TDTTF.h:101
MonitorElement * dttf_phi_wheel_2ndTrack[6]
Definition: L1TDTTF.h:84
MonitorElement * dttf_nTracksPerEvent_integ
Definition: L1TDTTF.h:98
DQMStore * dbe_
Definition: L1TDTTF.h:73
MonitorElement * dttf_pt_wheel_2ndTrack[6]
Definition: L1TDTTF.h:85
MonitorElement * dttf_eta_wheel_2ndTrack[6]
Definition: L1TDTTF.h:83
MonitorElement * dttf_qual[6][12]
Definition: L1TDTTF.h:91
MonitorElement * dttf_quality_wheel_2ndTrack[6]
Definition: L1TDTTF.h:78
MonitorElement * dttf_phi[6][12]
Definition: L1TDTTF.h:94
MonitorElement * dttf_gmt_missed
Definition: L1TDTTF.h:102
MonitorElement * dttf_quality_summary_wheel_2ndTrack[6]
Definition: L1TDTTF.h:79
MonitorElement * dttf_eta_fine_fraction[6][12]
Definition: L1TDTTF.h:92
MonitorElement * dttf_bx[6][12]
Definition: L1TDTTF.h:89
int nev_dttf_track2_
Definition: L1TDTTF.h:108
MonitorElement * dttf_gmt_ghost
Definition: L1TDTTF.h:103
MonitorElement * dttf_bx_2ndTrack[6][12]
Definition: L1TDTTF.h:90
int nev_dttf_
Definition: L1TDTTF.h:107
MonitorElement * dttf_q_wheel_2ndTrack[6]
Definition: L1TDTTF.h:86
MonitorElement * dttf_eta[6][12]
Definition: L1TDTTF.h:93
bool online_
Definition: L1TDTTF.h:71
MonitorElement * dttf_nTracksPerEv[6][12]
Definition: L1TDTTF.h:88
MonitorElement * dttf_pt[6][12]
Definition: L1TDTTF.h:95
MonitorElement * dttf_phi_eta_wheel_2ndTrack[6]
Definition: L1TDTTF.h:82
MonitorElement * dttf_spare
Definition: L1TDTTF.h:99
int nev_
Definition: L1TDTTF.h:106
TH2F * getTH2F(void) const
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:850
MonitorElement * dttf_q[6][12]
Definition: L1TDTTF.h:96
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
void setQualLabel(MonitorElement *me, int axis)
Definition: L1TDTTF.cc:1027
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
void L1TDTTF::beginLuminosityBlock ( edm::LuminosityBlock const &  lumiSeg,
edm::EventSetup const &  context 
)
inlineprotectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 51 of file L1TDTTF.h.

52  {};
void L1TDTTF::bookEta ( int  wh,
int &  nbins,
float &  start,
float &  stop 
)
private

Definition at line 1056 of file L1TDTTF.cc.

Referenced by beginJob().

1057 {
1058 
1059  switch ( wh ) {
1060  case 0 : start = 0; stop = 18; nbins = 18; break; // N2
1061  case 1 : start = 8; stop = 28; nbins = 20; break; // N1
1062  case 2 : start = 22; stop = 32; nbins = 10; break; // N0
1063  case 3 : start = 22; stop = 42; nbins = 20; break; // P0
1064  case 4 : start = 36; stop = 56; nbins = 20; break; // P1
1065  case 5 : start = 46; stop = 64; nbins = 18; break; // P2
1066  default : start = 0; stop = 0; nbins = 0; break; // BOH
1067  }
1068 
1069 }
tuple start
Check for commandline option errors.
Definition: dqm_diff.py:58
void L1TDTTF::endJob ( void  )
protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 603 of file L1TDTTF.cc.

References dbe_, nev_, nev_dttf_, nev_dttf_track2_, outputFile_, DQMStore::save(), and verbose_.

604 {
605  if (verbose_) {
606  edm::LogInfo("EndJob") << "L1TDTTF: end job....";
607  edm::LogInfo("EndJob") << "analyzed " << nev_ << " events";
608  edm::LogInfo("EndJob") << "containing at least one dttf track : "
609  << nev_dttf_;
610  edm::LogInfo("EndJob") << "containing two dttf tracks : "
611  << nev_dttf_track2_;
612  }
613 
614  if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
615 
616 }
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:2118
DQMStore * dbe_
Definition: L1TDTTF.h:73
std::string outputFile_
Definition: L1TDTTF.h:74
int nev_dttf_track2_
Definition: L1TDTTF.h:108
bool verbose_
Definition: L1TDTTF.h:72
int nev_dttf_
Definition: L1TDTTF.h:107
int nev_
Definition: L1TDTTF.h:106
void L1TDTTF::endLuminosityBlock ( edm::LuminosityBlock const &  lumiSeg,
edm::EventSetup const &  context 
)
inlineprotectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 54 of file L1TDTTF.h.

55  {};
void L1TDTTF::fillMEs ( std::vector< L1MuDTTrackCand > *  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

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 828 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_, and verbose_.

Referenced by analyze().

830 {
831 
832  L1MuDTTrackContainer::TrackContainer::const_iterator track
833  = trackContainer->begin();
834  L1MuDTTrackContainer::TrackContainer::const_iterator trackEnd
835  = trackContainer->end();
836 
837  for ( ; track != trackEnd; ++track ) {
838 
839  if ( verbose_ ) {
840  edm::LogInfo("L1TDTTF::Analyze") << "bx = " << track->bx();
841  edm::LogInfo("L1TDTTF::Analyze") << "quality (packed) = "
842  << track->quality_packed();
843  edm::LogInfo("L1TDTTF::Analyze") << "pt (packed) = "
844  << track->pt_packed()
845  << " , pt (GeV) = " << track->ptValue();
846  edm::LogInfo("L1TDTTF::Analyze") << "phi (packed) = "
847  << track->phi_packed()
848  << " , phi (rad) = " << track->phiValue();
849  edm::LogInfo("L1TDTTF::Analyze") << "charge (packed) = "
850  << track->charge_packed();
851  }
852 
853 
855  if ( ( track->whNum() == -1 ) && ! track->eta_packed() ) {
856  edm::LogInfo("L1TDTTF::Analyze") << "Skipping N0 with zero eta value";
857 
858  continue;
859  }
860 
861 
862  int bxindex = track->bx() + 1;
863  int se = track->scNum();
864  int sector = se + 1;
865  int whindex = track->whNum();
866 
867  whindex = ( whindex < 0 ) ? whindex + 3 : whindex + 2;
868 
869  if ( whindex < 0 || whindex > 5 ) {
870  edm::LogError("L1TDTTF::Analyze::WHEEL_ERROR") << track->whNum()
871  << "(" << whindex << ")";
872  continue;
873  }
874 
875  if ( se < 0 || se > 11 ) {
876  edm::LogError("L1TDTTF::Analyze::SECTOR_ERROR") << se;
877  continue;
878  }
879 
881 
883  // int phi_local = track->phi_packed();//range: 0 < phi_local < 31
884  // if ( phi_local > 15 ) phi_local -= 32; //range: -16 < phi_local < 15
885 
886  // int phi_global = phi_local + se * 12; //range: -16 < phi_global < 147
887  // if(phi_global < 0) phi_global += 144; //range: 0 < phi_global < 147
888  // if(phi_global > 143) phi_global -= 144; //range: 0 < phi_global < 143
889  // // float phi_phys = phi_global * 2.5 + 1.25;
890 
892  int phi_global = track->phi_packed();
893  phi_global = (phi_global > 15 ? phi_global - 32 : phi_global ) + se * 12;
894  if ( phi_global < -6 ) phi_global += 144; //range: 0 < phi_global < 147
895  if ( phi_global > 137 ) phi_global -= 144; //range: 0 < phi_global < 143
896 
897  // int eta_global = track->eta_packed();
898  // int eta_global = track->eta_packed() - 32;
899  // dttf_eta[bxindex][whindex][se]->Fill(eta_global);
900  // float eta_phys = eta_global / 2.4 ;
901 
905 
907  dttf_bx[whindex][se]->Fill(track->bx());
908 
910  if( track->TrkTag() == 1 ) {
911 
913  dttf_bx_2ndTrack[whindex][se]->Fill(track->bx());
914 
915  }
916 
920 
921 
922  if ( bxindex == 1 ) {
923 
925  ++numTracks[whindex][se];
926 
928  dttf_phi[whindex][se]->Fill(phi_global);
929 
931  dttf_qual[whindex][se]->Fill(track->quality_packed());
932 
934  dttf_pt[whindex][se]->Fill(track->pt_packed());
935 
937  dttf_q[whindex][se]->Fill(track->charge_packed());
938 
939 
941  dttf_eta[whindex][se]->Fill( track->eta_packed() );
942 
943  if( track->isFineHalo() ) {
944 
945  dttf_eta_fine_fraction[whindex][se]->Fill( 0 );
946 
948  dttf_phi_eta_fine_wheel[whindex]->Fill( track->eta_packed(), phi_global );
949 
950  } else {
951 
952  dttf_eta_fine_fraction[whindex][se]->Fill( 1 );
953 
955  dttf_phi_eta_coarse_wheel[whindex]->Fill( track->eta_packed(), phi_global );
956  }
957 
959  if ( online_ ) {
960  dttf_spare->Fill( whindex, sector );
961  }
962 
964  // dttf_qual_eta_integ->Fill(track->eta_packed(), track->quality_packed());
965 
967  if ( track->TrkTag() == 1 ) {
968 
970  dttf_phi_wheel_2ndTrack[whindex]->Fill(phi_global);
971 
973  dttf_pt_wheel_2ndTrack[whindex]->Fill(track->pt_packed());
974 
976  dttf_eta_wheel_2ndTrack[whindex]->Fill(track->eta_packed());
977 
979  dttf_quality_wheel_2ndTrack[whindex]->Fill(track->quality_packed());
980 
982  dttf_q_wheel_2ndTrack[whindex]->Fill(track->charge_packed());
983 
985  dttf_quality_summary_wheel_2ndTrack[whindex]->Fill( sector, track->quality_packed() );
986 
988  dttf_phi_eta_wheel_2ndTrack[whindex]->Fill( track->eta_packed(), phi_global );
989 
990  }
991 
993  bool match = false;
994  std::vector<L1MuRegionalCand>::iterator dttfCand;
996  unsigned int gmt_phi = ( phi_global < 0 ? phi_global + 144 : phi_global );
997 
998  for ( dttfCand = gmtDttfCands.begin(); dttfCand != gmtDttfCands.end();
999  ++dttfCand ) {
1000 
1002  if ( dttfCand->empty() ) continue;
1003  if ( ( dttfCand->phi_packed() == gmt_phi ) &&
1004  dttfCand->quality_packed() == track->quality_packed() ) {
1005  match = true;
1006  dttfCand->reset();
1007  break;
1008  }
1009 
1010 
1011  }
1012 
1013  if ( match ) {
1014  dttf_gmt_match->Fill( whindex, sector );
1015  } else {
1016  dttf_gmt_missed->Fill( whindex, sector );
1017  }
1018 
1019  }
1020 
1021  }
1022 
1023 }
MonitorElement * dttf_phi_eta_coarse_wheel[6]
Definition: L1TDTTF.h:81
MonitorElement * dttf_phi_eta_fine_wheel[6]
Definition: L1TDTTF.h:80
MonitorElement * dttf_gmt_match
Definition: L1TDTTF.h:101
MonitorElement * dttf_phi_wheel_2ndTrack[6]
Definition: L1TDTTF.h:84
void Fill(long long x)
MonitorElement * dttf_pt_wheel_2ndTrack[6]
Definition: L1TDTTF.h:85
MonitorElement * dttf_eta_wheel_2ndTrack[6]
Definition: L1TDTTF.h:83
MonitorElement * dttf_qual[6][12]
Definition: L1TDTTF.h:91
MonitorElement * dttf_quality_wheel_2ndTrack[6]
Definition: L1TDTTF.h:78
MonitorElement * dttf_phi[6][12]
Definition: L1TDTTF.h:94
MonitorElement * dttf_gmt_missed
Definition: L1TDTTF.h:102
MonitorElement * dttf_quality_summary_wheel_2ndTrack[6]
Definition: L1TDTTF.h:79
MonitorElement * dttf_eta_fine_fraction[6][12]
Definition: L1TDTTF.h:92
MonitorElement * dttf_bx[6][12]
Definition: L1TDTTF.h:89
int numTracks[6][12]
Definition: L1TDTTF.h:109
bool verbose_
Definition: L1TDTTF.h:72
MonitorElement * dttf_bx_2ndTrack[6][12]
Definition: L1TDTTF.h:90
MonitorElement * dttf_q_wheel_2ndTrack[6]
Definition: L1TDTTF.h:86
MonitorElement * dttf_eta[6][12]
Definition: L1TDTTF.h:93
bool online_
Definition: L1TDTTF.h:71
MonitorElement * dttf_pt[6][12]
Definition: L1TDTTF.h:95
MonitorElement * dttf_phi_eta_wheel_2ndTrack[6]
Definition: L1TDTTF.h:82
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:6
MonitorElement * dttf_spare
Definition: L1TDTTF.h:99
MonitorElement * dttf_q[6][12]
Definition: L1TDTTF.h:96
void L1TDTTF::setQualLabel ( MonitorElement me,
int  axis 
)
private

Definition at line 1027 of file L1TDTTF.cc.

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

Referenced by beginJob().

1028 {
1029 
1030  if( axis == 1 )
1031  me->setAxisTitle("Quality", axis);
1032  me->setBinLabel(1, "T34", axis);
1033  me->setBinLabel(2, "T23/24", axis);
1034  me->setBinLabel(3, "T12/13/14", axis);
1035  me->setBinLabel(4, "T234", axis);
1036  me->setBinLabel(5, "T134", axis);
1037  me->setBinLabel(6, "T123/124", axis);
1038  me->setBinLabel(7, "T1234", axis);
1039 }
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
void L1TDTTF::setWheelLabel ( MonitorElement me)
private

Definition at line 1042 of file L1TDTTF.cc.

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

Referenced by beginJob().

1043 {
1044  me->setAxisTitle("Wheel", 1);
1045  me->setBinLabel(1, "N2", 1);
1046  me->setBinLabel(2, "N1", 1);
1047  me->setBinLabel(3, "N0", 1);
1048  me->setBinLabel(4, "P0", 1);
1049  me->setBinLabel(5, "P1", 1);
1050  me->setBinLabel(6, "P2", 1);
1051 }
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)

Member Data Documentation

DQMStore* L1TDTTF::dbe_
private

Definition at line 73 of file L1TDTTF.h.

Referenced by beginJob(), endJob(), and L1TDTTF().

MonitorElement* L1TDTTF::dttf_bx[6][12]
private

Definition at line 89 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_bx_2ndTrack[6][12]
private

Definition at line 90 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_eta[6][12]
private

Definition at line 93 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_eta_fine_fraction[6][12]
private

Definition at line 92 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_eta_wheel_2ndTrack[6]
private

Definition at line 83 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_gmt_ghost
private

Definition at line 103 of file L1TDTTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTF::dttf_gmt_match
private

Definition at line 101 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_gmt_missed
private

Definition at line 102 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_nTracksPerEv[6][12]
private

Definition at line 88 of file L1TDTTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTF::dttf_nTracksPerEvent_integ
private

Definition at line 98 of file L1TDTTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTF::dttf_nTracksPerEvent_wheel[6]
private

Definition at line 77 of file L1TDTTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTF::dttf_phi[6][12]
private

Definition at line 94 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_phi_eta_coarse_wheel[6]
private

Definition at line 81 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_phi_eta_fine_wheel[6]
private

Definition at line 80 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_phi_eta_wheel_2ndTrack[6]
private

Definition at line 82 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_phi_wheel_2ndTrack[6]
private

Definition at line 84 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_pt[6][12]
private

Definition at line 95 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_pt_wheel_2ndTrack[6]
private

Definition at line 85 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_q[6][12]
private

Definition at line 96 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_q_wheel_2ndTrack[6]
private

Definition at line 86 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_qual[6][12]
private

Definition at line 91 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_quality_summary_wheel_2ndTrack[6]
private

Definition at line 79 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_quality_wheel_2ndTrack[6]
private

Definition at line 78 of file L1TDTTF.h.

Referenced by beginJob(), and fillMEs().

MonitorElement* L1TDTTF::dttf_spare
private

Definition at line 99 of file L1TDTTF.h.

Referenced by analyze(), beginJob(), and fillMEs().

edm::InputTag L1TDTTF::dttpgSource_
private

Definition at line 67 of file L1TDTTF.h.

Referenced by analyze(), and L1TDTTF().

edm::InputTag L1TDTTF::gmtSource_
private

Definition at line 68 of file L1TDTTF.h.

Referenced by analyze().

std::string L1TDTTF::l1tsubsystemfolder_
private

Definition at line 70 of file L1TDTTF.h.

Referenced by beginJob(), and L1TDTTF().

edm::InputTag L1TDTTF::muonCollectionLabel_
private

Definition at line 69 of file L1TDTTF.h.

Referenced by analyze().

int L1TDTTF::nev_
private

Definition at line 106 of file L1TDTTF.h.

Referenced by analyze(), beginJob(), and endJob().

int L1TDTTF::nev_dttf_
private

Definition at line 107 of file L1TDTTF.h.

Referenced by analyze(), beginJob(), and endJob().

int L1TDTTF::nev_dttf_track2_
private

Definition at line 108 of file L1TDTTF.h.

Referenced by analyze(), beginJob(), and endJob().

int L1TDTTF::numTracks[6][12]
private

Definition at line 109 of file L1TDTTF.h.

Referenced by analyze(), and fillMEs().

bool L1TDTTF::online_
private

Definition at line 71 of file L1TDTTF.h.

Referenced by analyze(), beginJob(), and fillMEs().

std::string L1TDTTF::outputFile_
private

Definition at line 74 of file L1TDTTF.h.

Referenced by endJob(), and L1TDTTF().

edm::InputTag L1TDTTF::trackInputTag_
private

Definition at line 75 of file L1TDTTF.h.

Referenced by analyze(), and L1TDTTF().

bool L1TDTTF::verbose_
private

Definition at line 72 of file L1TDTTF.h.

Referenced by analyze(), endJob(), fillMEs(), and L1TDTTF().