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 one::DQMEDAnalyzer< T >
	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

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_castor1
MonitorElement *	h_ep_castor2
MonitorElement *	h_ep_HFm1
MonitorElement *	h_ep_HFm2
MonitorElement *	h_ep_HFm3
MonitorElement *	h_ep_HFp1
MonitorElement *	h_ep_HFp2
MonitorElement *	h_ep_HFp3
MonitorElement *	h_ep_trackm1
MonitorElement *	h_ep_trackm2
MonitorElement *	h_ep_trackmid2
MonitorElement *	h_ep_trackmid3
MonitorElement *	h_ep_trackp1
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

Detailed Description

Definition at line 22 of file CentralityDQM.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 29 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 ( )

override

Definition at line 53 of file CentralityDQM.cc.

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

Member Function Documentation

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

overrideprotected

Definition at line 127 of file CentralityDQM.cc.

References angle(), hi::Castor1, hi::Castor2, 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(), edm::Event::getByToken(), hi::HFm1, hi::HFm2, hi::HFm3, hi::HFp1, hi::HFp2, hi::HFp3, edm::HandleBase::isValid(), reco::Centrality::multiplicityPixel(), reco::Centrality::NpixelTracks(), reco::Centrality::Ntracks(), reco::Centrality::NtracksEtaCut(), reco::Centrality::NtracksEtaPtCut(), reco::Centrality::NtracksPtCut(), hi::trackm1, hi::trackm2, hi::trackmid2, hi::trackmid3, hi::trackp1, hi::trackp2, 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_HFm1->Fill(rp[HFm1].angle(0));
  h_ep_HFp1->Fill(rp[HFp1].angle(0));
  h_ep_trackm1->Fill(rp[trackm1].angle(0));
  h_ep_trackp1->Fill(rp[trackp1].angle(0));
  h_ep_castor1->Fill(rp[Castor1].angle(0));


  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_castor2->Fill(rp[Castor2].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));

  }
}

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

overrideprivate

Definition at line 61 of file CentralityDQM.cc.

References DQMStore::IBooker::book1D(), and DQMStore::IBooker::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_HFm1 = bei.book1D("h_ep_HFm1", "h_ep_HFm1", 800,-psirange,psirange);
  h_ep_HFp1 = bei.book1D("h_ep_HFp1", "h_ep_HFp1", 800,-psirange,psirange);
  h_ep_trackm1 = bei.book1D("h_ep_trackm1", "h_ep_trackm1", 800,-psirange,psirange);
  h_ep_trackp1 = bei.book1D("h_ep_trackp1", "h_ep_trackp1", 800,-psirange,psirange);
  h_ep_castor1 = bei.book1D("h_ep_castor1", "h_ep_castor1", 800,-psirange,psirange);
  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_castor2 = bei.book1D("h_ep_castor2", "h_ep_castor2", 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

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

private

Definition at line 43 of file CentralityDQM.h.

edm::Handle<int> CentralityDQM::centralityBin_

private

Definition at line 54 of file CentralityDQM.h.

edm::InputTag CentralityDQM::centralityBinTag_

private

Definition at line 52 of file CentralityDQM.h.

edm::EDGetTokenT<int> CentralityDQM::centralityBinToken

private

Definition at line 53 of file CentralityDQM.h.

edm::InputTag CentralityDQM::centralityTag_

private

Definition at line 41 of file CentralityDQM.h.

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

private

Definition at line 42 of file CentralityDQM.h.

edm::InputTag CentralityDQM::eventplaneTag_

private

Definition at line 49 of file CentralityDQM.h.

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

private

Definition at line 50 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_cent_bin

private

Definition at line 89 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_castor1

private

Definition at line 95 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_castor2

private

Definition at line 102 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_HFm1

private

Definition at line 91 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_HFm2

private

Definition at line 97 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_HFm3

private

Definition at line 104 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_HFp1

private

Definition at line 92 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_HFp2

private

Definition at line 98 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_HFp3

private

Definition at line 105 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_trackm1

private

Definition at line 93 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_trackm2

private

Definition at line 100 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_trackmid2

private

Definition at line 99 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_trackmid3

private

Definition at line 106 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_trackp1

private

Definition at line 94 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_ep_trackp2

private

Definition at line 101 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiEB

private

Definition at line 76 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiEE

private

Definition at line 78 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiEEminus

private

Definition at line 80 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiEEplus

private

Definition at line 79 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiET

private

Definition at line 77 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiHF

private

Definition at line 68 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiHFhit

private

Definition at line 73 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiHFhitMinus

private

Definition at line 75 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiHFhitPlus

private

Definition at line 74 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiHFminus

private

Definition at line 70 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiHFminusEta4

private

Definition at line 72 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiHFplus

private

Definition at line 69 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiHFplusEta4

private

Definition at line 71 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiNpix

private

Definition at line 62 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiNpixelTracks

private

Definition at line 63 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiNtracks

private

Definition at line 64 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiNtracksEtaCut

private

Definition at line 66 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiNtracksEtaPtCut

private

Definition at line 67 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiNtracksPtCut

private

Definition at line 65 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiZDC

private

Definition at line 81 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiZDCminus

private

Definition at line 83 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_hiZDCplus

private

Definition at line 82 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_vertex_x

private

Definition at line 85 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_vertex_y

private

Definition at line 86 of file CentralityDQM.h.

MonitorElement* CentralityDQM::h_vertex_z

private

Definition at line 87 of file CentralityDQM.h.

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

private

Definition at line 47 of file CentralityDQM.h.

edm::InputTag CentralityDQM::vertexTag_

private

Definition at line 45 of file CentralityDQM.h.

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

private

Definition at line 46 of file CentralityDQM.h.