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 &&
440 DECorr_[1][4]->
Fill(ne->regionId().ieta(),
j->regionId().ieta());
441 DECorr_[2][4]->
Fill(ne->regionId().iphi(),
j->regionId().iphi());
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)
473 std::vector<int> bit =
SortMinBiasBit(hfs->etSum(2), hfs->etSum(3));
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);