File indexing completed on 2022-11-01 03:29:36
0001 #include "DQM/Physics/src/CentralitypADQM.h"
0002
0003 #include <memory>
0004
0005
0006 #include "DQMServices/Core/interface/DQMStore.h"
0007 #include "DQMServices/Core/interface/DQMEDAnalyzer.h"
0008
0009
0010 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0011 #include "FWCore/Framework/interface/Event.h"
0012 #include "FWCore/Framework/interface/EventSetup.h"
0013 #include "FWCore/Framework/interface/LuminosityBlock.h"
0014 #include "FWCore/ServiceRegistry/interface/Service.h"
0015 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0016
0017
0018 #include "DataFormats/VertexReco/interface/Vertex.h"
0019
0020 using namespace edm;
0021 using namespace reco;
0022 using namespace hi;
0023
0024
0025
0026
0027 CentralitypADQM::CentralitypADQM(const edm::ParameterSet& ps) {
0028 edm::LogInfo("CentralitypADQM") << " Starting CentralitypADQM "
0029 << "\n";
0030
0031 centralityTag_ = ps.getParameter<InputTag>("centralitycollection");
0032 centralityToken = consumes<reco::Centrality>(centralityTag_);
0033
0034 vertexTag_ = ps.getParameter<InputTag>("vertexcollection");
0035 vertexToken = consumes<std::vector<reco::Vertex> >(vertexTag_);
0036
0037 }
0038
0039
0040
0041
0042 CentralitypADQM::~CentralitypADQM() {
0043 edm::LogInfo("CentralitypADQM") << " Deleting CentralitypADQM "
0044 << "\n";
0045 }
0046
0047
0048
0049
0050 void CentralitypADQM::bookHistograms(DQMStore::IBooker& bei, edm::Run const&, edm::EventSetup const&) {
0051
0052
0053 bei.setCurrentFolder("Physics/CentralitypA/");
0054
0055 h_hiNpix = bei.book1D("h_hiNpix", "h_hiNpix", 750, 0, 75000);
0056
0057 h_hiNtracks = bei.book1D("h_hiNtracks", "h_hiNtracks", 500, 0, 5000);
0058 h_hiNtracksPtCut = bei.book1D("h_hiNtracksPtCut", "h_hiNtracksPtCut", 500, 0, 5000);
0059 h_hiNtracksEtaCut = bei.book1D("h_hiNtracksEtaCut", "h_hiNtracksEtaCut", 500, 0, 5000);
0060 h_hiNtracksEtaPtCut = bei.book1D("h_hiNtracksEtaPtCut", "h_hiNtracksEtaPtCut", 500, 0, 5000);
0061
0062 h_hiHF = bei.book1D("h_hiHF", "h_hiHF", 900, 0, 9000);
0063 h_hiHFplus = bei.book1D("h_hiHFplus", "h_hiHFplus", 900, 0, 9000);
0064 h_hiHFminus = bei.book1D("h_hiHFminus", "h_hiHFminus", 900, 0, 9000);
0065 h_hiHFplusEta4 = bei.book1D("h_hiHFplusEta4", "h_hiHFplusEta4", 900, 0, 9000);
0066 h_hiHFminusEta4 = bei.book1D("h_hiHFminusEta4", "h_hiHFminusEta4", 900, 0, 9000);
0067
0068 h_hiHFhit = bei.book1D("h_hiHFhit", "h_hiHFhit", 3000, 0, 300000);
0069 h_hiHFhitPlus = bei.book1D("h_hiHFhitPlus", "h_hiHFhitPlus", 2000, 0, 200000);
0070 h_hiHFhitMinus = bei.book1D("h_hiHFhitMinus", "h_hiHFhitMinus", 2000, 0, 200000);
0071
0072 h_hiEB = bei.book1D("h_hiEB", "h_hiEB", 600, 0, 6000);
0073 h_hiET = bei.book1D("h_hiET", "h_hiET", 600, 0, 6000);
0074 h_hiEE = bei.book1D("h_hiEE", "h_hiEE", 600, 0, 6000);
0075 h_hiEEplus = bei.book1D("h_hiEEplus", "h_hiEEplus", 600, 0, 6000);
0076 h_hiEEminus = bei.book1D("h_hiEEminus", "h_hiEEminus", 600, 0, 6000);
0077 h_hiZDC = bei.book1D("h_hiZDC", "h_hiZDC", 600, 0, 6000);
0078 h_hiZDCplus = bei.book1D("h_hiZDCplus", "h_hiZDCplus", 600, 0, 6000);
0079 h_hiZDCminus = bei.book1D("h_hiZDCminus", "h_hiZDCminus", 600, 0, 6000);
0080
0081 h_vertex_x = bei.book1D("h_vertex_x", "h_vertex_x", 400, -4, 4);
0082 h_vertex_y = bei.book1D("h_vertex_y", "h_vertex_y", 400, -4, 4);
0083 h_vertex_z = bei.book1D("h_vertex_z", "h_vertex_z", 400, -40, 40);
0084 }
0085
0086
0087
0088 void CentralitypADQM::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
0089 using namespace edm;
0090 edm::Handle<reco::Centrality> cent;
0091 iEvent.getByToken(centralityToken, cent);
0092
0093 if (cent.isValid()) {
0094 h_hiNpix->Fill(cent->multiplicityPixel());
0095 h_hiNtracks->Fill(cent->Ntracks());
0096
0097 h_hiNtracksPtCut->Fill(cent->NtracksPtCut());
0098 h_hiNtracksEtaCut->Fill(cent->NtracksEtaCut());
0099 h_hiNtracksEtaPtCut->Fill(cent->NtracksEtaPtCut());
0100
0101 h_hiHF->Fill(cent->EtHFtowerSum());
0102 h_hiHFplus->Fill(cent->EtHFtowerSumPlus());
0103 h_hiHFminus->Fill(cent->EtHFtowerSumMinus());
0104 h_hiHFplusEta4->Fill(cent->EtHFtruncatedPlus());
0105 h_hiHFminusEta4->Fill(cent->EtHFtruncatedMinus());
0106
0107 h_hiHFhit->Fill(cent->EtHFhitSum());
0108 h_hiHFhitPlus->Fill(cent->EtHFhitSumPlus());
0109 h_hiHFhitMinus->Fill(cent->EtHFhitSumMinus());
0110
0111 h_hiZDC->Fill(cent->zdcSum());
0112 h_hiZDCplus->Fill(cent->zdcSumPlus());
0113 h_hiZDCminus->Fill(cent->zdcSumMinus());
0114
0115 h_hiEEplus->Fill(cent->EtEESumPlus());
0116 h_hiEEminus->Fill(cent->EtEESumMinus());
0117 h_hiEE->Fill(cent->EtEESum());
0118 h_hiEB->Fill(cent->EtEBSum());
0119 h_hiET->Fill(cent->EtMidRapiditySum());
0120
0121 edm::Handle<std::vector<reco::Vertex> > vertex;
0122 iEvent.getByToken(vertexToken, vertex);
0123 h_vertex_x->Fill(vertex->begin()->x());
0124 h_vertex_y->Fill(vertex->begin()->y());
0125 h_vertex_z->Fill(vertex->begin()->z());
0126 }
0127 }