42 gctCenJetsDataSource_(ps.getParameter<
edm::
InputTag>(
"gctCentralJetsDataSource")),
43 gctForJetsDataSource_(ps.getParameter<
edm::
InputTag>(
"gctForwardJetsDataSource")),
44 gctTauJetsDataSource_(ps.getParameter<
edm::
InputTag>(
"gctTauJetsDataSource")),
45 gctEnergySumsDataSource_(ps.getParameter<
edm::
InputTag>(
"gctEnergySumsDataSource")),
46 gctIsoEmDataSource_(ps.getParameter<
edm::
InputTag>(
"gctIsoEmDataSource")),
47 gctNonIsoEmDataSource_(ps.getParameter<
edm::
InputTag>(
"gctNonIsoEmDataSource")),
51 gctCenJetsEmulSource_(ps.getParameter<
edm::
InputTag>(
"gctCentralJetsEmulSource")),
52 gctForJetsEmulSource_(ps.getParameter<
edm::
InputTag>(
"gctForwardJetsEmulSource")),
53 gctTauJetsEmulSource_(ps.getParameter<
edm::
InputTag>(
"gctTauJetsEmulSource")),
54 gctEnergySumsEmulSource_(ps.getParameter<
edm::
InputTag>(
"gctEnergySumsEmulSource")),
55 gctIsoEmEmulSource_(ps.getParameter<
edm::
InputTag>(
"gctIsoEmEmulSource")),
56 gctNonIsoEmEmulSource_(ps.getParameter<
edm::
InputTag>(
"gctNonIsoEmEmulSource")) {
138 "centrality and centrality ext E_{T} CORRELATION",
147 l1GctHFRingETSum_ = ibooker.
book1D(
"centrality+centralityExt Et",
"centrality+centralityExt E_{T}", 8, -0.5, 7.5);
153 ibooker.
book1D(
"centrality centralityExt ratio",
"centrality centralityExt ratio", 9, -0.5, 8.5);
191 const std::string clabel[8] = {
"cenJet",
"forJet",
"single track",
"isoEm",
"nonIsoEm",
"EtSum",
"MET",
"HTT"};
192 const std::string olabel[3] = {
"rank",
"eta",
"phi"};
194 unsigned int Bin[3][8] = {{64, 64, 64, 64, 64, 128, 128, 128},
197 float Min[3][8] = {{-0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5},
198 {-0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5},
199 {-0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5}};
200 float Max[3][8] = {{63.5, 63.5, 63.5, 63.5, 63.5, 1023.5, 1023.5, 1023.5},
201 {21.5, 21.5, 21.5, 21.5, 21.5, 21.5, 21.5, 21.5},
202 {17.5, 17.5, 17.5, 17.5, 17.5, 17.5, 71.5, 17.5}};
204 for (
int i = 0;
i < 3;
i++) {
205 for (
int j = 0;
j < 8;
j++) {
207 clabel[
j] + olabel[
i] +
" data vs emul",
217 centralityCorr_ = ibooker.
book2D(
"centrality data vs emul",
"centrality data vs emul", 8, -0.5, 7.5, 8, -0.5, 7.5);
219 ibooker.
book2D(
"centrality ext data vs emul",
"centrality ext data vs emul", 8, -0.5, 7.5, 8, -0.5, 7.5);
221 "Minimum Bias Trigger Data vs Emul",
"Minimum Bias Trigger Data vs Emul", 6, -0.5, 5.5, 6, -0.5, 5.5);
280 for (L1GctJetCandCollection::const_iterator cj = l1CenJets->begin(); cj != l1CenJets->end(); cj++) {
292 for (L1GctJetCandCollection::const_iterator
j = l1CenJetsEmul->begin();
j != l1CenJetsEmul->end();
j++) {
293 if (cj->bx() == 0 &&
j->bx() == 0 &&
297 DECorr_[1][0]->
Fill(cj->regionId().ieta(),
j->regionId().ieta());
298 DECorr_[2][0]->
Fill(cj->regionId().iphi(),
j->regionId().iphi());
302 for (L1GctJetCandCollection::const_iterator
j = l1CenJetsEmul->begin();
j != l1CenJetsEmul->end();
j++) {
306 for (L1GctJetCandCollection::const_iterator fj = l1ForJets->begin(); fj != l1ForJets->end(); fj++) {
318 for (L1GctJetCandCollection::const_iterator
j = l1ForJetsEmul->begin();
j != l1ForJetsEmul->end();
j++) {
319 if (fj->bx() == 0 &&
j->bx() == 0 &&
322 DECorr_[1][1]->
Fill(fj->regionId().ieta(),
j->regionId().ieta());
323 DECorr_[2][1]->
Fill(fj->regionId().iphi(),
j->regionId().iphi());
328 for (L1GctJetCandCollection::const_iterator tj = l1TauJets->begin(); tj != l1TauJets->end(); tj++) {
340 for (L1GctJetCandCollection::const_iterator
j = l1TauJetsEmul->begin();
j != l1TauJetsEmul->end();
j++) {
341 if (tj->bx() == 0 &&
j->bx() == 0 &&
344 DECorr_[1][2]->
Fill(tj->regionId().ieta(),
j->regionId().ieta());
345 DECorr_[2][2]->
Fill(tj->regionId().iphi(),
j->regionId().iphi());
350 for (L1GctEtMissCollection::const_iterator
met = l1EtMiss->begin();
met != l1EtMiss->end();
met++) {
352 if (
met->bx() == 0) {
353 if (
met->overFlow() == 0 &&
met->et() > 0) {
360 if (
met->overFlow() == 0 &&
met->et() > 0)
362 for (L1GctEtMissCollection::const_iterator
j = l1EtMissEmul->begin();
j != l1EtMissEmul->end();
j++) {
363 if (
met->bx() == 0 &&
j->bx() == 0) {
370 for (L1GctEtHadCollection::const_iterator ht = l1EtHad->begin(); ht != l1EtHad->end(); ht++) {
377 for (L1GctEtHadCollection::const_iterator
j = l1EtHadEmul->begin();
j != l1EtHadEmul->end();
j++) {
378 if (ht->bx() == 0 &&
j->bx() == 0) {
386 for (L1GctEtTotalCollection::const_iterator
et = l1EtTotal->begin();
et != l1EtTotal->end();
et++) {
393 for (L1GctEtTotalCollection::const_iterator
j = l1EtTotalEmul->begin();
j != l1EtTotalEmul->end();
j++) {
394 if (
et->bx() == 0 &&
j->bx() == 0) {
402 for (L1GctEmCandCollection::const_iterator ie = l1IsoEm->begin(); ie != l1IsoEm->end(); ie++) {
414 for (L1GctEmCandCollection::const_iterator
j = l1IsoEmEmul->begin();
j != l1IsoEmEmul->end();
j++) {
415 if (ie->bx() == 0 &&
j->bx() == 0 &&
418 DECorr_[1][3]->
Fill(ie->regionId().ieta(),
j->regionId().ieta());
419 DECorr_[2][3]->
Fill(ie->regionId().iphi(),
j->regionId().iphi());
424 for (L1GctEmCandCollection::const_iterator
ne = l1NonIsoEm->begin();
ne != l1NonIsoEm->end();
ne++) {
436 for (L1GctEmCandCollection::const_iterator
j = l1NonIsoEmEmul->begin();
j != l1NonIsoEmEmul->end();
j++) {
437 if (
ne->bx() == 0 &&
j->bx() == 0 &&
446 for (L1GctHFBitCountsCollection::const_iterator hfc = l1HFCounts->begin(); hfc != l1HFCounts->end(); hfc++) {
448 if (hfc->bx() == 0) {
459 for (
unsigned i = 0;
i < 4;
i++) {
464 for (L1GctHFRingEtSumsCollection::const_iterator hfs = l1HFSums->begin(); hfs != l1HFSums->end(); hfs++) {
465 if (hfs->bx() == 0) {
470 if (hfs->etSum(1) != 0)
474 for (std::vector<int>::const_iterator
it =
bit.begin();
it !=
bit.end();
it++) {
478 for (
unsigned i = 0;
i < 4;
i++) {
481 for (L1GctHFRingEtSumsCollection::const_iterator
j = l1HFSumsEmul->begin();
j != l1HFSumsEmul->end();
j++) {
482 if (hfs->bx() == 0 &&
j->bx() == 0 &&
486 std::vector<int> dbit =
SortMinBiasBit(hfs->etSum(2), hfs->etSum(3));
492 for (L1CaloRegionCollection::const_iterator
it = rgn->begin();
it != rgn->end();
it++) {
496 if (
it->gctEta() < 4) {
499 if (
it->gctEta() > 17) {
510 std::vector<int> Bit;
512 if ((
a + 1) / 4 > 0.5) {
518 if (
a == 2 ||
a == 3 ||
a == 6 ||
a == 7) {
524 Bit.push_back(
a % 2);
526 if ((
b + 1) / 4 > 0.5) {
532 if (
b == 2 ||
b == 3 ||
b == 6 ||
b == 7) {
538 Bit.push_back(
b % 2);
MonitorElement * l1GctEtTotal_
MonitorElement * l1GctHFRing1TowerCountPosEtaNegEta_
MonitorElement * HFPosEnergy_
T getParameter(std::string const &) const
MonitorElement * l1GctIsoEmOccEtaPhi_
edm::EDGetTokenT< L1GctHtMissCollection > l1HtMissEmulToken_
MonitorElement * l1GctEtHadOf_
virtual void setCurrentFolder(std::string const &fullpath)
edm::EDGetTokenT< L1GctHtMissCollection > l1HtMissDataToken_
MonitorElement * l1GctIsoEmRankEtaPhi_
MonitorElement * HFEnergy_
MonitorElement * l1GctHFRingTowerCountOccBx_
edm::EDGetTokenT< L1CaloRegionCollection > rctSource_L1CRCollection_
edm::EDGetTokenT< L1GctEtMissCollection > l1EtMissEmulToken_
MonitorElement * l1GctHFRing1ETSumPosEta_
MonitorElement * l1GctTauJetsEtEtaPhi_
edm::EDGetTokenT< L1GctEmCandCollection > gctIsoEmSourceEmulToken_
edm::EDGetTokenT< L1GctJetCandCollection > gctCenJetsSourceEmulToken_
const unsigned int EMETABINS
MonitorElement * l1GctForJetsRank_
MonitorElement * DECorr_[3][8]
MonitorElement * l1GctTauJetsRank_
MonitorElement * l1GctHFRing2TowerCountPosEta_
MonitorElement * l1GctEtHadOccBx_
MonitorElement * l1GctHFRing1PosEtaNegEta_
MonitorElement * l1GctEtMissOf_
edm::EDGetTokenT< L1GctJetCandCollection > gctTauJetsSourceEmulToken_
MonitorElement * l1GctHFRing1TowerCountNegEta_
MonitorElement * l1GctIsoEmRank_
edm::EDGetTokenT< L1GctEmCandCollection > gctIsoEmSourceDataToken_
edm::EDGetTokenT< L1GctHFBitCountsCollection > l1HFCountsDataToken_
edm::EDGetTokenT< L1GctEtHadCollection > l1EtHadEmulToken_
MonitorElement * l1GctAllJetsOccRankBx_
MonitorElement * l1GctHFRing2TowerCountPosEtaNegEta_
MonitorElement * l1GctHFRingETSumOccBx_
MonitorElement * l1GctEtMissOccBx_
edm::EDGetTokenT< L1GctEmCandCollection > gctNonIsoEmSourceEmulToken_
edm::EDGetTokenT< L1GctJetCandCollection > gctForJetsSourceEmulToken_
MonitorElement * centralityCorr_
edm::EDGetTokenT< L1GctEtMissCollection > l1EtMissDataToken_
const unsigned int R6BINS
virtual std::vector< int > SortMinBiasBit(uint16_t, uint16_t)
MonitorElement * centralityExtCorr_
MonitorElement * HFNegEnergy_
MonitorElement * l1GctCenJetsOccEtaPhi_
edm::EDGetTokenT< L1GctEmCandCollection > gctNonIsoEmSourceDataToken_
MonitorElement * l1GctForJetsEtEtaPhi_
MonitorElement * l1GctAllEmOccRankBx_
MonitorElement * l1GctEtMissPhi_
edm::EDGetTokenT< L1GctEtTotalCollection > l1EtTotalEmulToken_
MonitorElement * l1GctHFRingETDiff_
Abs< T >::type abs(const T &t)
MonitorElement * l1GctEtMiss_
MonitorElement * l1GctForJetsOccEtaPhi_
edm::EDGetTokenT< L1GctEtTotalCollection > l1EtTotalDataToken_
edm::EDGetTokenT< L1GctJetCandCollection > gctForJetsSourceDataToken_
MonitorElement * l1GctHFRing2TowerCountNegEta_
edm::EDGetTokenT< L1GctHFBitCountsCollection > l1HFCountsEmulToken_
edm::EDGetTokenT< L1GctHFRingEtSumsCollection > gctEnergySumsSourceEmulToken_
void analyze(const edm::Event &e, const edm::EventSetup &c) override
const unsigned int JETETABINS
const unsigned int METPHIBINS
const unsigned int R3BINS
MonitorElement * l1GctEtTotalEtHadCorr_
const unsigned int R12BINS
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
MonitorElement * l1GctHFRingRatioPosEta_
MonitorElement * l1GctCenJetsEtEtaPhi_
MonitorElement * l1GctHFRing1ETSumNegEta_
MonitorElement * l1GctEtTotalOccBx_
MonitorElement * l1GctHFRingETSum_
edm::EDGetTokenT< L1GctJetCandCollection > gctTauJetsSourceDataToken_
edm::EDGetTokenT< L1GctHFRingEtSumsCollection > gctEnergySumsSourceDataToken_
MonitorElement * l1GctEtHad_
MonitorElement * l1GctHFRing1TowerCountPosEta_
MonitorElement * l1GctTauJetsOccEtaPhi_
MonitorElement * MinBiasCorr_
MonitorElement * l1GctMinBiasBitHFEt_
const unsigned int BXBINS
edm::EDGetTokenT< L1GctJetCandCollection > gctCenJetsSourceDataToken_
std::vector< L1CaloRegion > L1CaloRegionCollection
edm::EDGetTokenT< L1GctEtHadCollection > l1EtHadDataToken_
MonitorElement * l1GctNonIsoEmOccEtaPhi_
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
const unsigned int PHIBINS
L1THIonImp(const edm::ParameterSet &ps)
void bookHistograms(DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &) override
MonitorElement * l1GctEtTotalOf_
MonitorElement * l1GctNonIsoEmRank_
const unsigned int OFBINS
MonitorElement * l1GctNonIsoEmRankEtaPhi_
MonitorElement * l1GctCenJetsRank_