CentralityDQM Class Reference

#include <CentralityDQM.h>

Inheritance diagram for CentralityDQM:

Public Member Functions
	CentralityDQM (const edm::ParameterSet &ps)
	~CentralityDQM () 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 (edm::Event const &e, edm::EventSetup const &eSetup) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	bookHistograms (DQMStore::IBooker &bei, edm::Run const &, edm::EventSetup const &) override

Private Attributes
edm::Handle< reco::Centrality >	centrality_
edm::Handle< int >	centralityBin_
edm::InputTag	centralityBinTag_
edm::EDGetTokenT< int >	centralityBinToken
edm::InputTag	centralityTag_
edm::EDGetTokenT< reco::Centrality >	centralityToken
edm::InputTag	eventplaneTag_
edm::EDGetTokenT< reco::EvtPlaneCollection >	eventplaneToken
MonitorElement *	h_cent_bin
MonitorElement *	h_ep_HFm2
MonitorElement *	h_ep_HFm3
MonitorElement *	h_ep_HFp2
MonitorElement *	h_ep_HFp3
MonitorElement *	h_ep_trackm2
MonitorElement *	h_ep_trackmid2
MonitorElement *	h_ep_trackmid3
MonitorElement *	h_ep_trackp2
MonitorElement *	h_hiEB
MonitorElement *	h_hiEE
MonitorElement *	h_hiEEminus
MonitorElement *	h_hiEEplus
MonitorElement *	h_hiET
MonitorElement *	h_hiHF
MonitorElement *	h_hiHFhit
MonitorElement *	h_hiHFhitMinus
MonitorElement *	h_hiHFhitPlus
MonitorElement *	h_hiHFminus
MonitorElement *	h_hiHFminusEta4
MonitorElement *	h_hiHFplus
MonitorElement *	h_hiHFplusEta4
MonitorElement *	h_hiNpix
MonitorElement *	h_hiNpixelTracks
MonitorElement *	h_hiNtracks
MonitorElement *	h_hiNtracksEtaCut
MonitorElement *	h_hiNtracksEtaPtCut
MonitorElement *	h_hiNtracksPtCut
MonitorElement *	h_hiZDC
MonitorElement *	h_hiZDCminus
MonitorElement *	h_hiZDCplus
MonitorElement *	h_vertex_x
MonitorElement *	h_vertex_y
MonitorElement *	h_vertex_z
edm::Handle< std::vector< reco::Vertex > >	vertex_
edm::InputTag	vertexTag_
edm::EDGetTokenT< std::vector< reco::Vertex > >	vertexToken

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 20 of file CentralityDQM.h.

Constructor & Destructor Documentation

CentralityDQM()

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

explicit

Definition at line 28 of file CentralityDQM.cc.

References edm::ParameterSet::getParameter().

                                                      {
  edm::LogInfo("CentralityDQM") << " Starting CentralityDQM "
                                << "\n";

  centralityTag_ = ps.getParameter<InputTag>("centralitycollection");
  centralityToken = consumes<reco::Centrality>(centralityTag_);

  centralityBinTag_ = (ps.getParameter<edm::InputTag>("centralitybincollection"));
  centralityBinToken = consumes<int>(centralityBinTag_);

  vertexTag_ = ps.getParameter<InputTag>("vertexcollection");
  vertexToken = consumes<std::vector<reco::Vertex> >(vertexTag_);

  eventplaneTag_ = ps.getParameter<InputTag>("eventplanecollection");
  eventplaneToken = consumes<reco::EvtPlaneCollection>(eventplaneTag_);

  // just to initialize
}

~CentralityDQM()

CentralityDQM::~CentralityDQM ( )

override

Definition at line 50 of file CentralityDQM.cc.

                              {
  edm::LogInfo("CentralityDQM") << " Deleting CentralityDQM "
                                << "\n";
}

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 110 of file CentralityDQM.cc.

References angle(), SiStripBadComponentsDQMServiceTemplate_cfg::ep, reco::Centrality::EtEBSum(), reco::Centrality::EtEESum(), reco::Centrality::EtEESumMinus(), reco::Centrality::EtEESumPlus(), reco::Centrality::EtHFhitSum(), reco::Centrality::EtHFhitSumMinus(), reco::Centrality::EtHFhitSumPlus(), reco::Centrality::EtHFtowerSum(), reco::Centrality::EtHFtowerSumMinus(), reco::Centrality::EtHFtowerSumPlus(), reco::Centrality::EtHFtruncatedMinus(), reco::Centrality::EtHFtruncatedPlus(), reco::Centrality::EtMidRapiditySum(), hi::HFm2, hi::HFm3, hi::HFp2, hi::HFp3, iEvent, edm::HandleBase::isValid(), reco::Centrality::multiplicityPixel(), reco::Centrality::NpixelTracks(), reco::Centrality::Ntracks(), reco::Centrality::NtracksEtaCut(), reco::Centrality::NtracksEtaPtCut(), reco::Centrality::NtracksPtCut(), hi::trackm2, hi::trackmid2, hi::trackmid3, hi::trackp2, bphysicsOniaDQM_cfi::vertex, reco::Centrality::zdcSum(), reco::Centrality::zdcSumMinus(), and reco::Centrality::zdcSumPlus().

                                                                               {
  using namespace edm;
  edm::Handle<reco::Centrality> cent;
  iEvent.getByToken(centralityToken, cent);  //_centralitytag comes from the cfg

  edm::Handle<int> cbin;
  iEvent.getByToken(centralityBinToken, cbin);

  // as an inputTag and is
  //"hiCentrality"
  edm::Handle<reco::EvtPlaneCollection> ep;
  iEvent.getByToken(eventplaneToken, ep);

  if (cent.isValid()) {
    int hibin = -999;
    if (cbin.isValid())
      hibin = *cbin;

    //  std::cout<<  " ------------------------------------- "  << hibin << std::endl;

    h_cent_bin->Fill(hibin);
    h_hiNpix->Fill(cent->multiplicityPixel());
    h_hiNpixelTracks->Fill(cent->NpixelTracks());
    h_hiNtracks->Fill(cent->Ntracks());  //

    h_hiNtracksPtCut->Fill(cent->NtracksPtCut());
    h_hiNtracksEtaCut->Fill(cent->NtracksEtaCut());
    h_hiNtracksEtaPtCut->Fill(cent->NtracksEtaPtCut());

    h_hiHF->Fill(cent->EtHFtowerSum());
    h_hiHFplus->Fill(cent->EtHFtowerSumPlus());
    h_hiHFminus->Fill(cent->EtHFtowerSumMinus());
    h_hiHFplusEta4->Fill(cent->EtHFtruncatedPlus());
    h_hiHFminusEta4->Fill(cent->EtHFtruncatedMinus());

    h_hiHFhit->Fill(cent->EtHFhitSum());
    h_hiHFhitPlus->Fill(cent->EtHFhitSumPlus());
    h_hiHFhitMinus->Fill(cent->EtHFhitSumMinus());

    h_hiZDC->Fill(cent->zdcSum());
    h_hiZDCplus->Fill(cent->zdcSumPlus());
    h_hiZDCminus->Fill(cent->zdcSumMinus());

    h_hiEEplus->Fill(cent->EtEESumPlus());
    h_hiEEminus->Fill(cent->EtEESumMinus());
    h_hiEE->Fill(cent->EtEESum());
    h_hiEB->Fill(cent->EtEBSum());
    h_hiET->Fill(cent->EtMidRapiditySum());

    edm::Handle<std::vector<reco::Vertex> > vertex;
    iEvent.getByToken(vertexToken, vertex);
    h_vertex_x->Fill(vertex->begin()->x());
    h_vertex_y->Fill(vertex->begin()->y());
    h_vertex_z->Fill(vertex->begin()->z());
  }

  if (ep.isValid()) {
    EvtPlaneCollection::const_iterator rp = ep->begin();

    h_ep_HFm2->Fill(rp[HFm2].angle(0));
    h_ep_HFp2->Fill(rp[HFp2].angle(0));
    h_ep_trackmid2->Fill(rp[trackmid2].angle(0));
    h_ep_trackm2->Fill(rp[trackm2].angle(0));
    h_ep_trackp2->Fill(rp[trackp2].angle(0));

    h_ep_HFm3->Fill(rp[HFm3].angle(0));
    h_ep_HFp3->Fill(rp[HFp3].angle(0));
    h_ep_trackmid3->Fill(rp[trackmid3].angle(0));
  }
}

bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 58 of file CentralityDQM.cc.

References dqm::implementation::IBooker::book1D(), and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                            {
  // void CentralityDQM::bookHistograms(DQMStore* bei){

  bei.setCurrentFolder("Physics/Centrality/");

  h_hiNpix = bei.book1D("h_hiNpix", "h_hiNpix", 750, 0, 75000);
  h_hiNpixelTracks = bei.book1D("h_hiNpixelTracks", "hiNpixelTracks", 500, 0, 5000);
  h_hiNtracks = bei.book1D("h_hiNtracks", "h_hiNtracks", 500, 0, 5000);
  h_hiNtracksPtCut = bei.book1D("h_hiNtracksPtCut", "h_hiNtracksPtCut", 500, 0, 5000);
  h_hiNtracksEtaCut = bei.book1D("h_hiNtracksEtaCut", "h_hiNtracksEtaCut", 500, 0, 5000);
  h_hiNtracksEtaPtCut = bei.book1D("h_hiNtracksEtaPtCut", "h_hiNtracksEtaPtCut", 500, 0, 5000);

  h_hiHF = bei.book1D("h_hiHF", "h_hiHF", 900, 0, 9000);
  h_hiHFplus = bei.book1D("h_hiHFplus", "h_hiHFplus", 900, 0, 9000);
  h_hiHFminus = bei.book1D("h_hiHFminus", "h_hiHFminus", 900, 0, 9000);
  h_hiHFplusEta4 = bei.book1D("h_hiHFplusEta4", "h_hiHFplusEta4", 900, 0, 9000);
  h_hiHFminusEta4 = bei.book1D("h_hiHFminusEta4", "h_hiHFminusEta4", 900, 0, 9000);

  h_hiHFhit = bei.book1D("h_hiHFhit", "h_hiHFhit", 3000, 0, 300000);
  h_hiHFhitPlus = bei.book1D("h_hiHFhitPlus", "h_hiHFhitPlus", 2000, 0, 200000);
  h_hiHFhitMinus = bei.book1D("h_hiHFhitMinus", "h_hiHFhitMinus", 2000, 0, 200000);

  h_hiEB = bei.book1D("h_hiEB", "h_hiEB", 600, 0, 6000);
  h_hiET = bei.book1D("h_hiET", "h_hiET", 600, 0, 6000);
  h_hiEE = bei.book1D("h_hiEE", "h_hiEE", 600, 0, 6000);
  h_hiEEplus = bei.book1D("h_hiEEplus", "h_hiEEplus", 600, 0, 6000);
  h_hiEEminus = bei.book1D("h_hiEEminus", "h_hiEEminus", 600, 0, 6000);
  h_hiZDC = bei.book1D("h_hiZDC", "h_hiZDC", 600, 0, 6000);
  h_hiZDCplus = bei.book1D("h_hiZDCplus", "h_hiZDCplus", 600, 0, 6000);
  h_hiZDCminus = bei.book1D("h_hiZDCminus", "h_hiZDCminus", 600, 0, 6000);

  h_vertex_x = bei.book1D("h_vertex_x", "h_vertex_x", 400, -4, 4);
  h_vertex_y = bei.book1D("h_vertex_y", "h_vertex_y", 400, -4, 4);
  h_vertex_z = bei.book1D("h_vertex_z", "h_vertex_z", 400, -40, 40);

  h_cent_bin = bei.book1D("h_cent_bin", "h_cent_bin", 200, 0, 200);

  Double_t psirange = 4;
  bei.setCurrentFolder("Physics/Centrality/EventPlane/");
  h_ep_HFm2 = bei.book1D("h_ep_HFm2", "h_ep_HFm2", 800, -psirange, psirange);
  h_ep_HFp2 = bei.book1D("h_ep_HFp2", "h_ep_HFp2", 800, -psirange, psirange);
  h_ep_trackmid2 = bei.book1D("h_ep_trackmid2", "h_ep_trackmid2", 800, -psirange, psirange);
  h_ep_trackm2 = bei.book1D("h_ep_trackm2", "h_ep_trackm2", 800, -psirange, psirange);
  h_ep_trackp2 = bei.book1D("h_ep_trackp2", "h_ep_trackp2", 800, -psirange, psirange);
  h_ep_HFm3 = bei.book1D("h_ep_HFm3", "h_ep_HFm3", 800, -psirange, psirange);
  h_ep_HFp3 = bei.book1D("h_ep_HFp3", "h_ep_HFp3", 800, -psirange, psirange);
  h_ep_trackmid3 = bei.book1D("h_ep_trackmid3", "h_ep_trackmid3", 800, -psirange, psirange);
}

Member Data Documentation

centrality_

edm::Handle<reco::Centrality> CentralityDQM::centrality_

private

Definition at line 39 of file CentralityDQM.h.

centralityBin_

edm::Handle<int> CentralityDQM::centralityBin_

private

Definition at line 50 of file CentralityDQM.h.

centralityBinTag_

edm::InputTag CentralityDQM::centralityBinTag_

private

Definition at line 48 of file CentralityDQM.h.

centralityBinToken

edm::EDGetTokenT<int> CentralityDQM::centralityBinToken

private

Definition at line 49 of file CentralityDQM.h.

centralityTag_

edm::InputTag CentralityDQM::centralityTag_

private

Definition at line 37 of file CentralityDQM.h.

centralityToken

edm::EDGetTokenT<reco::Centrality> CentralityDQM::centralityToken

private

Definition at line 38 of file CentralityDQM.h.

eventplaneTag_

edm::InputTag CentralityDQM::eventplaneTag_

private

Definition at line 45 of file CentralityDQM.h.

eventplaneToken

edm::EDGetTokenT<reco::EvtPlaneCollection> CentralityDQM::eventplaneToken

private

Definition at line 46 of file CentralityDQM.h.

h_cent_bin

MonitorElement* CentralityDQM::h_cent_bin

private

Definition at line 84 of file CentralityDQM.h.

h_ep_HFm2

MonitorElement* CentralityDQM::h_ep_HFm2

private

Definition at line 86 of file CentralityDQM.h.

h_ep_HFm3

MonitorElement* CentralityDQM::h_ep_HFm3

private

Definition at line 92 of file CentralityDQM.h.

h_ep_HFp2

MonitorElement* CentralityDQM::h_ep_HFp2

private

Definition at line 87 of file CentralityDQM.h.

h_ep_HFp3

MonitorElement* CentralityDQM::h_ep_HFp3

private

Definition at line 93 of file CentralityDQM.h.

h_ep_trackm2

MonitorElement* CentralityDQM::h_ep_trackm2

private

Definition at line 89 of file CentralityDQM.h.

h_ep_trackmid2

MonitorElement* CentralityDQM::h_ep_trackmid2

private

Definition at line 88 of file CentralityDQM.h.

h_ep_trackmid3

MonitorElement* CentralityDQM::h_ep_trackmid3

private

Definition at line 94 of file CentralityDQM.h.

h_ep_trackp2

MonitorElement* CentralityDQM::h_ep_trackp2

private

Definition at line 90 of file CentralityDQM.h.

h_hiEB

MonitorElement* CentralityDQM::h_hiEB

private

Definition at line 71 of file CentralityDQM.h.

h_hiEE

MonitorElement* CentralityDQM::h_hiEE

private

Definition at line 73 of file CentralityDQM.h.

h_hiEEminus

MonitorElement* CentralityDQM::h_hiEEminus

private

Definition at line 75 of file CentralityDQM.h.

h_hiEEplus

MonitorElement* CentralityDQM::h_hiEEplus

private

Definition at line 74 of file CentralityDQM.h.

h_hiET

MonitorElement* CentralityDQM::h_hiET

private

Definition at line 72 of file CentralityDQM.h.

h_hiHF

MonitorElement* CentralityDQM::h_hiHF

private

Definition at line 63 of file CentralityDQM.h.

h_hiHFhit

MonitorElement* CentralityDQM::h_hiHFhit

private

Definition at line 68 of file CentralityDQM.h.

h_hiHFhitMinus

MonitorElement* CentralityDQM::h_hiHFhitMinus

private

Definition at line 70 of file CentralityDQM.h.

h_hiHFhitPlus

MonitorElement* CentralityDQM::h_hiHFhitPlus

private

Definition at line 69 of file CentralityDQM.h.

h_hiHFminus

MonitorElement* CentralityDQM::h_hiHFminus

private

Definition at line 65 of file CentralityDQM.h.

h_hiHFminusEta4

MonitorElement* CentralityDQM::h_hiHFminusEta4

private

Definition at line 67 of file CentralityDQM.h.

h_hiHFplus

MonitorElement* CentralityDQM::h_hiHFplus

private

Definition at line 64 of file CentralityDQM.h.

h_hiHFplusEta4

MonitorElement* CentralityDQM::h_hiHFplusEta4

private

Definition at line 66 of file CentralityDQM.h.

h_hiNpix

MonitorElement* CentralityDQM::h_hiNpix

private

Definition at line 57 of file CentralityDQM.h.

h_hiNpixelTracks

MonitorElement* CentralityDQM::h_hiNpixelTracks

private

Definition at line 58 of file CentralityDQM.h.

h_hiNtracks

MonitorElement* CentralityDQM::h_hiNtracks

private

Definition at line 59 of file CentralityDQM.h.

h_hiNtracksEtaCut

MonitorElement* CentralityDQM::h_hiNtracksEtaCut

private

Definition at line 61 of file CentralityDQM.h.

h_hiNtracksEtaPtCut

MonitorElement* CentralityDQM::h_hiNtracksEtaPtCut

private

Definition at line 62 of file CentralityDQM.h.

h_hiNtracksPtCut

MonitorElement* CentralityDQM::h_hiNtracksPtCut

private

Definition at line 60 of file CentralityDQM.h.

h_hiZDC

MonitorElement* CentralityDQM::h_hiZDC

private

Definition at line 76 of file CentralityDQM.h.

h_hiZDCminus

MonitorElement* CentralityDQM::h_hiZDCminus

private

Definition at line 78 of file CentralityDQM.h.

h_hiZDCplus

MonitorElement* CentralityDQM::h_hiZDCplus

private

Definition at line 77 of file CentralityDQM.h.

h_vertex_x

MonitorElement* CentralityDQM::h_vertex_x

private

Definition at line 80 of file CentralityDQM.h.

h_vertex_y

MonitorElement* CentralityDQM::h_vertex_y

private

Definition at line 81 of file CentralityDQM.h.

h_vertex_z

MonitorElement* CentralityDQM::h_vertex_z

private

Definition at line 82 of file CentralityDQM.h.

vertex_

edm::Handle<std::vector<reco::Vertex> > CentralityDQM::vertex_

private

Definition at line 43 of file CentralityDQM.h.

vertexTag_

edm::InputTag CentralityDQM::vertexTag_

private

Definition at line 41 of file CentralityDQM.h.

vertexToken

edm::EDGetTokenT<std::vector<reco::Vertex> > CentralityDQM::vertexToken

private

Definition at line 42 of file CentralityDQM.h.