L1TRPCTF Class Reference

#include <L1TRPCTF.h>

Inheritance diagram for L1TRPCTF:

Classes
struct	BxDelays

Public Member Functions
	L1TRPCTF (const edm::ParameterSet &ps)
	~L1TRPCTF () 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
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	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 &ibooker, const edm::Run &, const edm::EventSetup &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
MonitorElement *	m_bxDiff
std::set< unsigned long long int >	m_globBX
long long int	m_lastUsedBxInBxdiff
MonitorElement *	m_muonsEtaPhi [3]
MonitorElement *	m_qualVsEta [3]
int	nev_
int	nevRPC_
std::string	output_dir_
MonitorElement *	rpctfbx
MonitorElement *	rpctfchargevalue [3]
MonitorElement *	rpctfcratesynchro [12]
MonitorElement *	rpctfetavalue [3]
MonitorElement *	rpctfntrack_b [3]
MonitorElement *	rpctfntrack_e [3]
MonitorElement *	rpctfphivalue [3]
MonitorElement *	rpctfptvalue [3]
MonitorElement *	rpctfquality [3]
edm::EDGetTokenT< L1MuGMTReadoutCollection >	rpctfSource_
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 41 of file L1TRPCTF.h.

Constructor & Destructor Documentation

L1TRPCTF()

L1TRPCTF::L1TRPCTF

(

const edm::ParameterSet &

ps

)

Definition at line 18 of file L1TRPCTF.cc.

References gather_cfg::cout, edm::ParameterSet::getUntrackedParameter(), and verbose_.

    : rpctfSource_(consumes<L1MuGMTReadoutCollection>(ps.getParameter<InputTag>("rpctfSource"))),
      //    digiSource_( ps.getParameter< InputTag >("rpctfRPCDigiSource") ),
      //   m_rpcDigiFine(false),
      //    m_useRpcDigi(true),
      m_lastUsedBxInBxdiff(0),
      output_dir_(ps.getUntrackedParameter<string>("output_dir"))
//    m_rpcDigiWithBX0(0),
//    m_rpcDigiWithBXnon0(0)

{
  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);

  if (verbose_)
    cout << "L1TRPCTF: constructor...." << endl;
}

~L1TRPCTF()

L1TRPCTF::~L1TRPCTF ( )

override

Definition at line 36 of file L1TRPCTF.cc.

{}

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 169 of file L1TRPCTF.cc.

References L1TRPCTF::BxDelays::bx, gather_cfg::cout, change_name::diff, HLT_2022v12_cff::distance, MillePedeFileConverter_cfg::e, L1TRPCTF::BxDelays::eta_t, dqm::impl::MonitorElement::Fill(), L1MuGMTReadoutCollection::getRecords(), mps_fire::i, edm::HandleBase::isValid(), dqmiolumiharvest::j, m_bxDiff, m_globBX, m_lastUsedBxInBxdiff, m_muonsEtaPhi, m_qualVsEta, SiStripPI::min, nev_, L1TRPCTF::BxDelays::phi_p, edm::Handle< T >::product(), rpctfbx, rpctfchargevalue, rpctfcratesynchro, rpctfetavalue, rpctfntrack_b, rpctfntrack_e, rpctfphivalue, rpctfptvalue, rpctfquality, rpctfSource_, hgcalTopologyTester_cfi::sector2, hgcalTowerProducer_cfi::tower, and verbose_.

                                                          {
  nev_++;
  if (verbose_)
    cout << "L1TRPCTF: analyze...." << endl;

  edm::Handle<L1MuGMTReadoutCollection> pCollection;
  e.getByToken(rpctfSource_, pCollection);

  if (!pCollection.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find L1MuGMTReadoutCollection";
    return;
  }

  L1MuGMTReadoutCollection const* gmtrc = pCollection.product();
  vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
  vector<L1MuGMTReadoutRecord>::const_iterator RRItr;

  std::vector<int> nrpctftrack_b(3, 0);
  std::vector<int> nrpctftrack_e(3, 0);

  vector<L1TRPCTF::BxDelays> all_bxdelays;

  for (RRItr = gmt_records.begin(); RRItr != gmt_records.end(); RRItr++) {
    if (verbose_)
      cout << "Readout Record " << RRItr->getBxInEvent() << endl;

    vector<vector<L1MuRegionalCand> > brlAndFwdCands;
    brlAndFwdCands.push_back(RRItr->getBrlRPCCands());
    brlAndFwdCands.push_back(RRItr->getFwdRPCCands());

    int beIndex = 0;
    vector<vector<L1MuRegionalCand> >::iterator RPCTFCands = brlAndFwdCands.begin();
    for (; RPCTFCands != brlAndFwdCands.end(); ++RPCTFCands) {
      for (vector<L1MuRegionalCand>::const_iterator ECItr = RPCTFCands->begin(); ECItr != RPCTFCands->end(); ++ECItr) {
        int bxindex = 1;  // bx == 0
        if (ECItr->bx() > 0)
          bxindex = 2;
        if (ECItr->bx() < 0)
          bxindex = 0;

        if (!ECItr->empty()) {
          if (beIndex == 0)
            ++nrpctftrack_b[bxindex];
          if (beIndex == 1)
            ++nrpctftrack_e[bxindex];

          if (verbose_)
            cout << "RPCTFCand bx " << ECItr->bx() << endl;

          int tower = ECItr->eta_packed();
          if (tower > 16) {
            tower = -((~tower & 63) + 1);
          }

          rpctfbx->Fill(ECItr->bx());

          rpctfetavalue[bxindex]->Fill(tower);
          if (verbose_)
            cout << "\tRPCTFCand eta value " << ECItr->etaValue() << endl;

          rpctfphivalue[bxindex]->Fill(ECItr->phi_packed());
          if (verbose_)
            cout << "\tRPCTFCand phi value " << ECItr->phiValue() << endl;

          rpctfptvalue[bxindex]->Fill(ECItr->ptValue());
          if (verbose_)
            cout << "\tRPCTFCand pt value " << ECItr->ptValue() << endl;

          rpctfchargevalue[bxindex]->Fill(ECItr->chargeValue());
          if (verbose_)
            cout << "\tRPCTFCand charge value " << ECItr->chargeValue() << endl;

          rpctfquality[bxindex]->Fill(ECItr->quality());
          if (verbose_)
            cout << "\tRPCTFCand quality " << ECItr->quality() << endl;

          m_qualVsEta[bxindex]->Fill(tower, ECItr->quality());
          m_muonsEtaPhi[bxindex]->Fill(tower, ECItr->phi_packed());

          BxDelays bx_del;
          bx_del.bx = ECItr->bx();
          bx_del.eta_t = tower;
          bx_del.phi_p = ECItr->phi_packed();
          all_bxdelays.push_back(bx_del);

        }  // if !empty
      }    // end candidates iteration
      ++beIndex;
    }  // end brl/endcap iteration
  }    // end GMT records iteration

  for (int bxI = 0; bxI < 3; ++bxI) {
    rpctfntrack_b[bxI]->Fill(nrpctftrack_b[bxI]);
    rpctfntrack_e[bxI]->Fill(nrpctftrack_e[bxI]);
  }

  for (unsigned int i = 0; i < all_bxdelays.size(); i++) {
    int sector = ((all_bxdelays[i].phi_p + 142) % 144) / 12;
    if (sector > 11 || sector < 0)
      continue;
    int eta_tower = all_bxdelays[i].eta_t;
    for (unsigned int j = 0; j < all_bxdelays.size(); j++) {
      if (i == j)
        continue;
      int sector2 = ((all_bxdelays[j].phi_p + 142) % 144) / 12;

      int distance_cut = 1;
      int distance = ((sector + 12) - sector2) % 12;
      distance = min(distance, 11 - distance);
      if (distance < distance_cut)
        continue;

      int bxDiff = all_bxdelays[i].bx - all_bxdelays[j].bx;
      rpctfcratesynchro[sector]->Fill(bxDiff, eta_tower);
    }
  }

  while (m_globBX.begin() != m_globBX.end()) {
    long long int diff = *m_globBX.begin() - m_lastUsedBxInBxdiff;  // first entry will go to overflow bin, ignore
    m_bxDiff->Fill(diff);
    m_lastUsedBxInBxdiff = *m_globBX.begin();
    m_globBX.erase(m_globBX.begin());
  }

  if (verbose_)
    cout << "L1TRPCTF: end job...." << endl;
  LogInfo("EndJob") << "analyzed " << nev_ << " events";
}

bookHistograms()

void L1TRPCTF::bookHistograms ( DQMStore::IBooker &  ibooker, const edm::Run &  , const edm::EventSetup &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 38 of file L1TRPCTF.cc.

References b, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), simKBmtfDigis_cfi::bx, mps_fire::i, createfilelist::int, cmsLHEtoEOSManager::l, m_bxDiff, m_muonsEtaPhi, m_qualVsEta, nev_, nevRPC_, EcalTangentSkim_cfg::o, output_dir_, rpctfbx, rpctfchargevalue, rpctfcratesynchro, rpctfetavalue, rpctfntrack_b, rpctfntrack_e, rpctfphivalue, rpctfptvalue, rpctfquality, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), and contentValuesCheck::ss.

                                                                                           {
  nev_ = 0;
  nevRPC_ = 0;

  ostringstream oDir;
  oDir << output_dir_ << "/CrateSynchroHistograms/";
  ibooker.setCurrentFolder(oDir.str());
  for (unsigned int i = 0; i < 12; i++) {
    ostringstream o;
    o << "RPCTF_crate_" << i << "_synchro";
    rpctfcratesynchro[i] = ibooker.book2D(o.str(), o.str(), 5, -2.5, 2.5, 33, -16.5, 16.5);
    for (int bx = -2; bx < 3; ++bx) {
      ostringstream b;
      b << "BX=" << bx;
      rpctfcratesynchro[i]->setBinLabel(bx + 3, b.str(), 1);
    }
    rpctfcratesynchro[i]->setAxisTitle("Tower", 2);
  }
  ibooker.setCurrentFolder(output_dir_);

  rpctfetavalue[1] = ibooker.book1D("RPCTF_eta_value_bx0", "RPCTF eta value bx=0", 33, -16.5, 16.5);
  rpctfetavalue[2] = ibooker.book1D("RPCTF_eta_value_bx+", "RPCTF eta value bx>0", 33, -16.5, 16.5);
  rpctfetavalue[0] = ibooker.book1D("RPCTF_eta_value_bx-", "RPCTF eta value bx<0", 33, -16.5, 16.5);

  rpctfphivalue[1] = ibooker.book1D("RPCTF_phi_value_bx0", "RPCTF phi value bx=0", 144, -0.5, 143.5);
  rpctfphivalue[2] = ibooker.book1D("RPCTF_phi_value_bx+", "RPCTF phi value bx>0", 144, -0.5, 143.5);
  rpctfphivalue[0] = ibooker.book1D("RPCTF_phi_value_bx-", "RPCTF phi value bx<0", 144, -0.5, 143.5);

  rpctfptvalue[1] = ibooker.book1D("RPCTF_pt_value_bx0", "RPCTF pt value bx=0", 160, -0.5, 159.5);
  rpctfptvalue[2] = ibooker.book1D("RPCTF_pt_value_bx+", "RPCTF pt value bx>0", 160, -0.5, 159.5);
  rpctfptvalue[0] = ibooker.book1D("RPCTF_pt_value_bx-", "RPCTF pt value bx<0", 160, -0.5, 159.5);

  rpctfchargevalue[1] = ibooker.book1D("RPCTF_charge_value_bx0", "RPCTF charge value bx=0", 3, -1.5, 1.5);
  rpctfchargevalue[2] = ibooker.book1D("RPCTF_charge_value_bx+", "RPCTF charge value bx>0", 3, -1.5, 1.5);
  rpctfchargevalue[0] = ibooker.book1D("RPCTF_charge_value_bx-", "RPCTF charge value bx<01", 3, -1.5, 1.5);

  rpctfquality[1] = ibooker.book1D("RPCTF_quality", "RPCTF quality bx=0", 6, -0.5, 5.5);
  rpctfquality[2] = ibooker.book1D("RPCTF_quality_bx+", "RPCTF quality bx>0", 6, -0.5, 5.5);
  rpctfquality[0] = ibooker.book1D("RPCTF_quality_bx-", "RPCTF quality bx<0", 6, -0.5, 5.5);

  rpctfntrack_b[1] = ibooker.book1D("RPCTF_ntrack_brl_bx0", "RPCTF number of tracks - barrel, bx=0", 5, -0.5, 4.5);
  rpctfntrack_b[2] = ibooker.book1D("RPCTF_ntrack_brl_bx+", "RPCTF number of tracks - barrel, bx>0", 5, -0.5, 4.5);
  rpctfntrack_b[0] = ibooker.book1D("RPCTF_ntrack_brl_bx-", "RPCTF number of tracks - barrel, bx<0", 5, -0.5, 4.5);

  rpctfntrack_e[1] = ibooker.book1D("RPCTF_ntrack_fwd_bx0", "RPCTF number of tracks - endcap, bx=0", 5, -0.5, 4.5);
  rpctfntrack_e[2] = ibooker.book1D("RPCTF_ntrack_fwd_bx+", "RPCTF number of tracks - endcap, bx>0", 5, -0.5, 4.5);
  rpctfntrack_e[0] = ibooker.book1D("RPCTF_ntrack_fwd_bx-", "RPCTF number of tracks - endcap, bx<0", 5, -0.5, 4.5);

  m_qualVsEta[1] = ibooker.book2D("RPCTF_quality_vs_eta_bx0",
                                  "RPCTF quality vs eta, bx=0",
                                  33,
                                  -16.5,
                                  16.5,
                                  6,
                                  -0.5,
                                  5.5);  // Currently only 0...3 quals are possible
  m_qualVsEta[2] = ibooker.book2D("RPCTF_quality_vs_eta_bx+",
                                  "RPCTF quality vs eta, bx>0",
                                  33,
                                  -16.5,
                                  16.5,
                                  6,
                                  -0.5,
                                  5.5);  // Currently only 0...3 quals are possible
  m_qualVsEta[0] = ibooker.book2D("RPCTF_quality_vs_eta_bx-",
                                  "RPCTF quality vs eta, bx<0",
                                  33,
                                  -16.5,
                                  16.5,
                                  6,
                                  -0.5,
                                  5.5);  // Currently only 0...3 quals are possible

  m_muonsEtaPhi[1] =
      ibooker.book2D("RPCTF_muons_eta_phi_bx0", "RPCTF occupancy(eta,phi), bx=0", 33, -16.5, 16.5, 144, -0.5, 143.5);
  m_muonsEtaPhi[2] =
      ibooker.book2D("RPCTF_muons_eta_phi_bx+", "RPCTF occupancy(eta,phi), bx>0", 33, -16.5, 16.5, 144, -0.5, 143.5);
  m_muonsEtaPhi[0] =
      ibooker.book2D("RPCTF_muons_eta_phi_bx-", "RPCTF occupancy(eta,phi), bx<0", 33, -16.5, 16.5, 144, -0.5, 143.5);

  rpctfbx = ibooker.book1D("RPCTF_bx", "RPCTF bx distribiution", 7, -3.5, 3.5);

  //axis labels
  for (int l = 0; l < 3; ++l) {
    m_muonsEtaPhi[l]->setAxisTitle("tower", 1);
    m_qualVsEta[l]->setAxisTitle("tower");
    rpctfetavalue[l]->setAxisTitle("tower");

    m_muonsEtaPhi[l]->setAxisTitle("phi", 2);
    rpctfphivalue[l]->setAxisTitle("phi");
  }

  // set phi bin labels
  for (int i = 0; i < 12; ++i) {
    //float lPhi  = (30./360)*i*2*3.14;
    int lPhi = 30 * i;
    int lBin = int((30. / 360) * i * 144) + 1;
    std::stringstream ss;
    ss << "phi=" << lPhi;
    for (int l = 0; l < 3; ++l) {
      rpctfphivalue[l]->setBinLabel(lBin, ss.str());
      m_muonsEtaPhi[l]->setBinLabel(lBin, ss.str(), 2);
    }
  }

  /*
    // set TC numbers on phi axis
    for (int tc = 0; tc < 12 ; ++tc ){
      int lBin  = (tc*12+3+1)%144;
      std::stringstream ss;
      ss << "TC" <<tc;
      for (int l = 0; l<3; ++l){
        rpctfphivalue[l]->setBinLabel(lBin,ss.str());
        m_muonsEtaPhi[l]->setBinLabel(lBin,ss.str(), 2);
      }
  }*/

  // set eta bin labels
  for (int i = -16; i < 17; ++i) {
    std::stringstream ss;
    ss << i;
    for (int l = 0; l < 3; ++l) {
      rpctfetavalue[l]->setBinLabel(i + 17, ss.str());
      m_muonsEtaPhi[l]->setBinLabel(i + 17, ss.str(), 1);
      m_qualVsEta[l]->setBinLabel(i + 17, ss.str());
    }
  }

  m_bxDiff = ibooker.book1D("RPCTF_bx_diff", "RPCTrigger - bx difference", 12000, -.5, 11999.5);
}

Member Data Documentation

m_bxDiff

MonitorElement* L1TRPCTF::m_bxDiff

private

Definition at line 71 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

m_globBX

std::set<unsigned long long int> L1TRPCTF::m_globBX

private

Definition at line 74 of file L1TRPCTF.h.

Referenced by analyze().

m_lastUsedBxInBxdiff

long long int L1TRPCTF::m_lastUsedBxInBxdiff

private

Definition at line 82 of file L1TRPCTF.h.

Referenced by analyze().

m_muonsEtaPhi

MonitorElement* L1TRPCTF::m_muonsEtaPhi[3]

private

Definition at line 68 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

m_qualVsEta

MonitorElement* L1TRPCTF::m_qualVsEta[3]

private

Definition at line 67 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

nev_

int L1TRPCTF::nev_

private

Definition at line 78 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

nevRPC_

int L1TRPCTF::nevRPC_

private

Definition at line 79 of file L1TRPCTF.h.

Referenced by bookHistograms().

output_dir_

std::string L1TRPCTF::output_dir_

private

Definition at line 83 of file L1TRPCTF.h.

Referenced by bookHistograms().

rpctfbx

MonitorElement* L1TRPCTF::rpctfbx

private

Definition at line 66 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

rpctfchargevalue

MonitorElement* L1TRPCTF::rpctfchargevalue[3]

private

Definition at line 62 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

rpctfcratesynchro

MonitorElement* L1TRPCTF::rpctfcratesynchro[12]

private

Definition at line 72 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

rpctfetavalue

MonitorElement* L1TRPCTF::rpctfetavalue[3]

private

Definition at line 59 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

rpctfntrack_b

MonitorElement* L1TRPCTF::rpctfntrack_b[3]

private

Definition at line 64 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

rpctfntrack_e

MonitorElement* L1TRPCTF::rpctfntrack_e[3]

private

Definition at line 65 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

rpctfphivalue

MonitorElement* L1TRPCTF::rpctfphivalue[3]

private

Definition at line 60 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

rpctfptvalue

MonitorElement* L1TRPCTF::rpctfptvalue[3]

private

Definition at line 61 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

rpctfquality

MonitorElement* L1TRPCTF::rpctfquality[3]

private

Definition at line 63 of file L1TRPCTF.h.

Referenced by analyze(), and bookHistograms().

rpctfSource_

edm::EDGetTokenT<L1MuGMTReadoutCollection> L1TRPCTF::rpctfSource_

private

Definition at line 76 of file L1TRPCTF.h.

Referenced by analyze().

verbose_

bool L1TRPCTF::verbose_

private

Definition at line 80 of file L1TRPCTF.h.

Referenced by analyze(), and L1TRPCTF().