#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 42 of file L1TRPCTF.h.

Constructor & Destructor Documentation

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

override

Definition at line 36 of file L1TRPCTF.cc.

36 {}

Member Function Documentation

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_FULL_cff::distance, L1TRPCTF::BxDelays::eta_t, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), L1MuGMTReadoutCollection::getRecords(), mps_fire::i, edm::HandleBase::isValid(), dqmiolumiharvest::j, m_bxDiff, m_globBX, m_lastUsedBxInBxdiff, m_muonsEtaPhi, m_qualVsEta, min(), nev_, L1TRPCTF::BxDelays::phi_p, edm::Handle< T >::product(), rpctfbx, rpctfchargevalue, rpctfcratesynchro, rpctfetavalue, rpctfntrack_b, rpctfntrack_e, rpctfphivalue, rpctfptvalue, rpctfquality, rpctfSource_, 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";
 }

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(), makePileupJSON::bx, mps_fire::i, cmsLHEtoEOSManager::l, m_bxDiff, m_muonsEtaPhi, m_qualVsEta, nev_, nevRPC_, class-composition::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);
 }