Inheritance diagram for CTPPSCommonDQMSource:

Classes
struct	ArmPlots
	plots related to one arm More...

struct	GlobalPlots
	plots related to the whole system More...

Public Member Functions
	CTPPSCommonDQMSource (const edm::ParameterSet &ps)

	~CTPPSCommonDQMSource () override

Public Member Functions inherited from one::DQMEDAnalyzer< edm::LuminosityBlockCache< std::vector< int > > >
	DQMEDAnalyzer ()=default

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

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

	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override

void	analyzeCTPPSRecord (edm::Event const &event, edm::EventSetup const &eventSetup)

void	analyzeProtons (edm::Event const &event, edm::EventSetup const &eventSetup)

void	analyzeTracks (edm::Event const &event, edm::EventSetup const &eventSetup)

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

std::shared_ptr< std::vector< int > >	globalBeginLuminosityBlock (const edm::LuminosityBlock &iLumi, const edm::EventSetup &c) const override

void	globalEndLuminosityBlock (const edm::LuminosityBlock &iLumi, const edm::EventSetup &c) override

Private Attributes
std::map< unsigned int, ArmPlots >	armPlots

const edm::EDGetTokenT< CTPPSRecord >	ctppsRecordToken

int	currentLS

int	endLS

GlobalPlots	globalPlots

std::vector< int >	rpstate

const edm::EDGetTokenT< std::vector< CTPPSLocalTrackLite > >	tokenLocalTrackLite

const edm::EDGetTokenT< std::vector< reco::ForwardProton > >	tokenRecoProtons

const unsigned int	verbosity

Static Private Attributes
static int	MAX_LUMIS = 6000

static int	MAX_VBINS = 18

Detailed Description

Definition at line 30 of file CTPPSCommonDQMSource.cc.

Constructor & Destructor Documentation

CTPPSCommonDQMSource::CTPPSCommonDQMSource ( const edm::ParameterSet & ps )

Definition at line 249 of file CTPPSCommonDQMSource.cc.

References currentLS, endLS, and rpstate.

                                                                     :
   verbosity(ps.getUntrackedParameter<unsigned int>("verbosity", 0)),
   ctppsRecordToken(consumes<CTPPSRecord>(ps.getUntrackedParameter<edm::InputTag>("ctppsmetadata"))),
   tokenLocalTrackLite(consumes< vector<CTPPSLocalTrackLite> >(ps.getParameter<edm::InputTag>("tagLocalTrackLite"))),
   tokenRecoProtons(consumes<std::vector<reco::ForwardProton>>(ps.getParameter<InputTag>("tagRecoProtons")))
 {
   currentLS = 0;
   endLS = 0;
   rpstate.clear();
 }

CTPPSCommonDQMSource::~CTPPSCommonDQMSource ( )

override

Definition at line 262 of file CTPPSCommonDQMSource.cc.

263 {

264 }

Member Function Documentation

void CTPPSCommonDQMSource::analyze	(	edm::Event const &	e,
		edm::EventSetup const &	eSetup
	)

overrideprotected

Definition at line 282 of file CTPPSCommonDQMSource.cc.

References analyzeCTPPSRecord(), analyzeProtons(), and analyzeTracks().

 {
   analyzeCTPPSRecord(event, eventSetup);
   analyzeTracks(event, eventSetup);
   analyzeProtons(event, eventSetup);
 }

void CTPPSCommonDQMSource::analyzeCTPPSRecord	(	edm::Event const &	event,
		edm::EventSetup const &	eventSetup
	)

protected

Definition at line 291 of file CTPPSCommonDQMSource.cc.

References ctppsRecordToken, edm::Event::getLuminosityBlock(), mps_fire::i, edm::LuminosityBlock::index(), edm::HandleBase::isValid(), rpstate, CTPPSRecord::status(), and verbosity.

Referenced by analyze().

 {
   Handle<CTPPSRecord> hCTPPSRecord;
   event.getByToken(ctppsRecordToken, hCTPPSRecord);
 
   if (!hCTPPSRecord.isValid())
   {
     if (verbosity)
       LogProblem("CTPPSCommonDQMSource") << "ERROR in CTPPSCommonDQMSource::analyzeCTPPSRecord > input not available.";
 
     return;
   }
 
   auto& rpstate = *luminosityBlockCache(event.getLuminosityBlock().index());
   if (rpstate.empty())
   {
     rpstate.reserve(CTPPSRecord::RomanPot::Last);
     for (uint8_t i = 0; i < CTPPSRecord::RomanPot::Last; ++i)
       rpstate.push_back(hCTPPSRecord->status(i));
   }  
 }

void CTPPSCommonDQMSource::analyzeProtons	(	edm::Event const &	event,
		edm::EventSetup const &	eventSetup
	)

protected

Definition at line 466 of file CTPPSCommonDQMSource.cc.

References armPlots, edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, nanoDQM_cff::plots, reco::ForwardProton::sector45, reco::ForwardProton::sector56, tokenRecoProtons, and verbosity.

Referenced by analyze().

 {
   // get event data
   Handle<vector<reco::ForwardProton>> hRecoProtons;
   event.getByToken(tokenRecoProtons, hRecoProtons);
 
   // check validity
   if (!hRecoProtons.isValid())
   {
     if (verbosity)
       LogProblem("CTPPSCommonDQMSource") << "ERROR in CTPPSCommonDQMSource::analyzeProtons > input not available.";
 
     return;
   }
 
   // loop over protons
   for (auto &p : *hRecoProtons)
   {
     if (!p.validFit())
       continue;
 
     signed int armIndex = -1;
     if (p.lhcSector() == reco::ForwardProton::LHCSector::sector45)
       armIndex = 0;
     if (p.lhcSector() == reco::ForwardProton::LHCSector::sector56)
       armIndex = 1;
     if (armIndex < 0)
       continue;
 
     auto &plots = armPlots[armIndex];
 
     plots.h_proton_xi->Fill(p.xi());
     plots.h_proton_t->Fill(fabs(p.t()));
     plots.h_proton_time->Fill(p.time());
   }
 }

void CTPPSCommonDQMSource::analyzeTracks	(	edm::Event const &	event,
		edm::EventSetup const &	eventSetup
	)

protected

Definition at line 315 of file CTPPSCommonDQMSource.cc.

References CTPPSDetId::arm(), armPlots, edm::EventBase::bunchCrossing(), CTPPSCommonDQMSource::GlobalPlots::events_per_bx, CTPPSCommonDQMSource::GlobalPlots::events_per_bx_short, MonitorElement::Fill(), globalPlots, CTPPSCommonDQMSource::GlobalPlots::h_trackCorr_hor, CTPPSCommonDQMSource::GlobalPlots::h_trackCorr_vert, training_settings::idx, edm::HandleBase::isValid(), TtSemiLepHitFitProducer_Electrons_cfi::mTop, AlCaHLTBitMon_ParallelJobs::p, nanoDQM_cff::plots, CTPPSDetId::rp(), tokenLocalTrackLite, and verbosity.

Referenced by analyze().

 {
   // get event data
   Handle< vector<CTPPSLocalTrackLite> > hTracks;
   event.getByToken(tokenLocalTrackLite, hTracks);
 
   // check validity
   if (!hTracks.isValid())
   {
     if (verbosity)
       LogProblem("CTPPSCommonDQMSource") << "ERROR in CTPPSCommonDQMSource::analyzeTracks > input not available.";
 
     return;
   }
 
   //------------------------------
   // collect indeces of RP with tracks, for each correlation plot
   set<signed int> s_rp_idx_global_hor, s_rp_idx_global_vert;
   map<unsigned int, set<signed int>> ms_rp_idx_arm;
 
   for (auto &tr : *hTracks)
   {
     const CTPPSDetId rpId(tr.getRPId());
     const unsigned int arm = rpId.arm();
     const unsigned int stNum = rpId.station();
     const unsigned int rpNum = rpId.rp();
     const unsigned int stRPNum = stNum * 10 + rpNum;
 
     {
       signed int idx = -1;
       if (stRPNum ==  3) idx = 0;
       if (stRPNum == 23) idx = 1;
       if (stRPNum == 16) idx = 2;
 
       if (idx >= 0)
         s_rp_idx_global_hor.insert(3*arm + idx);
     }
 
     {
       signed int idx = -1;
       if (stRPNum ==  4) idx = 0;
       if (stRPNum ==  5) idx = 1;
       if (stRPNum == 24) idx = 2;
       if (stRPNum == 25) idx = 3;
 
       if (idx >= 0)
         s_rp_idx_global_vert.insert(4*arm + idx);
     }
 
     {
       signed int idx = -1;
       if (stRPNum ==  3) idx = 0;
       if (stRPNum ==  4) idx = 1;
       if (stRPNum ==  5) idx = 2;
       if (stRPNum == 16) idx = 3;
       if (stRPNum == 23) idx = 4;
       if (stRPNum == 24) idx = 5;
       if (stRPNum == 25) idx = 6;
 
       const signed int hor = ((rpNum == 2) || (rpNum == 3) || (rpNum == 6)) ? 1 : 0;
 
       if (idx >= 0)
         ms_rp_idx_arm[arm].insert(idx * 10 + hor);
     }
   }
 
   //------------------------------
   // Global Plots
 
   globalPlots.events_per_bx->Fill(event.bunchCrossing());
   globalPlots.events_per_bx_short->Fill(event.bunchCrossing());
 
   for (const auto &idx1 : s_rp_idx_global_hor)
     for (const auto &idx2 : s_rp_idx_global_hor)
       globalPlots.h_trackCorr_hor->Fill(idx1, idx2);
 
   for (const auto &idx1 : s_rp_idx_global_vert)
     for (const auto &idx2 : s_rp_idx_global_vert)
       globalPlots.h_trackCorr_vert->Fill(idx1, idx2);
 
   //------------------------------
   // Arm Plots
 
   map<unsigned int, set<unsigned int>> mTop, mHor, mBot;
 
   for (auto &tr : *hTracks)
   {
     CTPPSDetId rpId(tr.getRPId());
     const unsigned int rpNum = rpId.rp();
     const unsigned int armIdx = rpId.arm();
 
     if (rpNum == 0 || rpNum == 4)
       mTop[armIdx].insert(rpId);
     if (rpNum == 2 || rpNum == 3 || rpNum == 6)
       mHor[armIdx].insert(rpId);
     if (rpNum == 1 || rpNum == 5)
       mBot[armIdx].insert(rpId);
 
     auto &ap = armPlots[rpId.arm()];
     unsigned int rpDecId = rpId.arm()*100 + rpId.station()*10 + rpId.rp();
 
     // fill in reference tracking-RP plots
     {
       auto it = ap.trackingRPPlots.find(rpDecId);
       if (it != ap.trackingRPPlots.end())
       {
         it->second.h_x->Fill(tr.getX());
         it->second.h_y->Fill(tr.getY());
       }
     }
 
     // fill in reference timing-RP plots
     {
       auto it = ap.timingRPPlots.find(rpDecId);
       if (it != ap.timingRPPlots.end())
       {
         it->second.h_x->Fill(tr.getX());
         it->second.h_time->Fill(tr.getTime());
       }
     }
   }
 
   for (auto &p : armPlots)
   {
     p.second.h_numRPWithTrack_top->Fill(mTop[p.first].size());
     p.second.h_numRPWithTrack_hor->Fill(mHor[p.first].size());
     p.second.h_numRPWithTrack_bot->Fill(mBot[p.first].size());
   }
 
   //------------------------------
   // Correlation plots
 
   for (const auto &ap : ms_rp_idx_arm)
   {
     auto &plots = armPlots[ap.first];
 
     for (const auto &idx1 : ap.second)
     {
       for (const auto &idx2 : ap.second)
       {
         plots.h_trackCorr->Fill(idx1/10, idx2/10);
 
         if ((idx1 % 10) != (idx2 % 10))
           plots.h_trackCorr_overlap->Fill(idx1/10, idx2/10);
       }
     }
   }
 }

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

overrideprotected

Definition at line 268 of file CTPPSCommonDQMSource.cc.

References armPlots, globalPlots, and CTPPSCommonDQMSource::GlobalPlots::Init().

 {
   // global plots
   globalPlots.Init(ibooker);
 
   // loop over arms
   for (unsigned int arm = 0; arm < 2; arm++)
   {
     armPlots[arm] = ArmPlots(ibooker, arm);
   }
 }

std::shared_ptr< std::vector< int > > CTPPSCommonDQMSource::globalBeginLuminosityBlock	(	const edm::LuminosityBlock &	iLumi,
		const edm::EventSetup &	c
	)		const

overrideprotected

Definition at line 505 of file CTPPSCommonDQMSource.cc.

 {
   return std::make_shared<std::vector<int>>();
 }

void CTPPSCommonDQMSource::globalEndLuminosityBlock	(	const edm::LuminosityBlock &	iLumi,
		const edm::EventSetup &	c
	)

overrideprotected

Definition at line 512 of file CTPPSCommonDQMSource.cc.

References currentLS, DEFINE_FWK_MODULE, globalPlots, mps_fire::i, edm::LuminosityBlockBase::id(), edm::LuminosityBlock::index(), edm::LuminosityBlockID::luminosityBlock(), rpstate, CTPPSCommonDQMSource::GlobalPlots::RPState, and MonitorElement::setBinContent().

 {
   auto const& rpstate = *luminosityBlockCache(iLumi.index());
   auto currentLS = iLumi.id().luminosityBlock();
   for(std::vector<int>::size_type i=0; i<rpstate.size();i++)
     globalPlots.RPState->setBinContent(currentLS, i+1, rpstate[i]);
 }