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"))
117 l1GctHFRingETSum_ = ibooker.
book1D(
"centrality+centralityExt Et",
"centrality+centralityExt E_{T}", 8, -0.5, 7.5);
155 const std::string clabel[8] = {
"cenJet",
"forJet",
"single track",
"isoEm",
"nonIsoEm",
"EtSum",
"MET",
"HTT"};
156 const std::string olabel[3] = {
"rank",
"eta",
"phi"};
158 unsigned int Bin[3][8] = {{64,64,64,64,64,128,128,128},
161 float Min[3][8] = {{-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5},
162 {-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5},
163 {-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5}};
164 float Max[3][8] = {{63.5,63.5,63.5,63.5,63.5,1023.5,1023.5,1023.5},
165 {21.5,21.5,21.5,21.5,21.5,21.5,21.5,21.5},
166 {17.5,17.5,17.5,17.5,17.5,17.5,71.5,17.5}};
168 for(
int i=0;
i<3;
i++){
169 for(
int j=0;
j<8;
j++){
170 DECorr_[
i][
j] = ibooker.
book2D(clabel[
j]+olabel[
i]+
" data vs emul", clabel[
j]+olabel[
i]+
" data vs emul", Bin[
i][
j], Min[
i][j], Max[
i][j], Bin[
i][j], Min[
i][j], Max[
i][j]);
174 centralityCorr_=ibooker.
book2D(
"centrality data vs emul",
"centrality data vs emul",8,-0.5,7.5,8,-0.5,7.5);
175 centralityExtCorr_=ibooker.
book2D(
"centrality ext data vs emul",
"centrality ext data vs emul",8,-0.5,7.5,8,-0.5,7.5);
176 MinBiasCorr_ = ibooker.
book2D(
"Minimum Bias Trigger Data vs Emul",
"Minimum Bias Trigger Data vs Emul",6,-0.5,5.5,6,-0.5,5.5);
242 for (L1GctJetCandCollection::const_iterator cj = l1CenJets->begin();cj != l1CenJets->end(); cj++) {
253 for (L1GctJetCandCollection::const_iterator
j = l1CenJetsEmul->begin();
j != l1CenJetsEmul->end();
j++) {
257 DECorr_[1][0]->
Fill(cj->regionId().ieta(),
j->regionId().ieta());
258 DECorr_[2][0]->
Fill(cj->regionId().iphi(),
j->regionId().iphi());
262 for (L1GctJetCandCollection::const_iterator
j = l1CenJetsEmul->begin();
j != l1CenJetsEmul->end();
j++) {
266 for (L1GctJetCandCollection::const_iterator fj = l1ForJets->begin(); fj != l1ForJets->end(); fj++) {
277 for (L1GctJetCandCollection::const_iterator
j = l1ForJetsEmul->begin();
j != l1ForJetsEmul->end();
j++) {
280 DECorr_[1][1]->
Fill(fj->regionId().ieta(),
j->regionId().ieta());
281 DECorr_[2][1]->
Fill(fj->regionId().iphi(),
j->regionId().iphi());
286 for (L1GctJetCandCollection::const_iterator tj = l1TauJets->begin(); tj != l1TauJets->end(); tj++) {
297 for (L1GctJetCandCollection::const_iterator
j = l1TauJetsEmul->begin();
j!= l1TauJetsEmul->end();
j++) {
300 DECorr_[1][2]->
Fill(tj->regionId().ieta(),
j->regionId().ieta());
301 DECorr_[2][2]->
Fill(tj->regionId().iphi(),
j->regionId().iphi());
306 for (L1GctEtMissCollection::const_iterator
met = l1EtMiss->begin();
met != l1EtMiss->end();
met++) {
309 if (
met->overFlow() == 0 &&
met->et() > 0) {
317 for (L1GctEtMissCollection::const_iterator
j = l1EtMissEmul->begin();
j != l1EtMissEmul->end();
j++){
318 if (
met->bx()==0 &&
j->bx()==0){
325 for (L1GctEtHadCollection::const_iterator ht = l1EtHad->begin(); ht != l1EtHad->end(); ht++) {
332 for (L1GctEtHadCollection::const_iterator
j = l1EtHadEmul->begin();
j != l1EtHadEmul->end();
j++){
333 if (ht->bx()==0 &&
j->bx()==0){
341 for (L1GctEtTotalCollection::const_iterator et = l1EtTotal->begin(); et != l1EtTotal->end(); et++) {
348 for (L1GctEtTotalCollection::const_iterator
j = l1EtTotalEmul->begin();
j != l1EtTotalEmul->end();
j++){
349 if (et->bx()==0 &&
j->bx()==0){
357 for (L1GctEmCandCollection::const_iterator ie=l1IsoEm->begin(); ie!=l1IsoEm->end(); ie++) {
368 for (L1GctEmCandCollection::const_iterator
j = l1IsoEmEmul->begin();
j!= l1IsoEmEmul->end();
j++) {
371 DECorr_[1][3]->
Fill(ie->regionId().ieta(),
j->regionId().ieta());
372 DECorr_[2][3]->
Fill(ie->regionId().iphi(),
j->regionId().iphi());
377 for (L1GctEmCandCollection::const_iterator
ne=l1NonIsoEm->begin();
ne!=l1NonIsoEm->end();
ne++) {
388 for (L1GctEmCandCollection::const_iterator
j = l1NonIsoEmEmul->begin();
j!= l1NonIsoEmEmul->end();
j++) {
397 for (L1GctHFBitCountsCollection::const_iterator hfc=l1HFCounts->begin(); hfc!=l1HFCounts->end(); hfc++){
410 for (
unsigned i=0;
i<4;
i++){
415 for (L1GctHFRingEtSumsCollection::const_iterator hfs=l1HFSums->begin(); hfs!=l1HFSums->end(); hfs++){
423 std::vector<int> bit =
SortMinBiasBit(hfs->etSum(2),hfs->etSum(3));
424 for(std::vector<int>::const_iterator it=bit.begin(); it !=bit.end(); it++){
428 for (
unsigned i=0;
i<4;
i++){
431 for (L1GctHFRingEtSumsCollection::const_iterator
j=l1HFSumsEmul->begin();
j!=l1HFSumsEmul->end();
j++){
435 std::vector<int> dbit =
SortMinBiasBit(hfs->etSum(2),hfs->etSum(3));
441 for (L1CaloRegionCollection::const_iterator it=rgn->begin(); it!=rgn->end(); it++){
445 if (it->gctEta() < 4){
448 if (it->gctEta() > 17){
459 std::vector<int> Bit;
467 if (a==2 || a==3 || a==6 || a==7){
481 if (b==2 || b==3 || b==6 || b==7){
MonitorElement * l1GctEtTotal_
T getParameter(std::string const &) const
const unsigned int R6BINS
const unsigned int BXBINS
MonitorElement * l1GctHFRing1TowerCountPosEtaNegEta_
MonitorElement * HFPosEnergy_
const unsigned int OFBINS
MonitorElement * l1GctIsoEmOccEtaPhi_
edm::EDGetTokenT< L1GctHtMissCollection > l1HtMissEmulToken_
MonitorElement * l1GctEtHadOf_
const unsigned int EMETABINS
edm::EDGetTokenT< L1GctHtMissCollection > l1HtMissDataToken_
MonitorElement * l1GctIsoEmRankEtaPhi_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
MonitorElement * HFEnergy_
MonitorElement * l1GctHFRingTowerCountOccBx_
edm::EDGetTokenT< L1CaloRegionCollection > rctSource_L1CRCollection_
edm::EDGetTokenT< L1GctEtMissCollection > l1EtMissEmulToken_
MonitorElement * l1GctHFRing1ETSumPosEta_
void analyze(const edm::Event &e, const edm::EventSetup &c)
MonitorElement * l1GctTauJetsEtEtaPhi_
edm::EDGetTokenT< L1GctEmCandCollection > gctIsoEmSourceEmulToken_
edm::EDGetTokenT< L1GctJetCandCollection > gctCenJetsSourceEmulToken_
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_
virtual std::vector< int > SortMinBiasBit(uint16_t, uint16_t)
MonitorElement * centralityExtCorr_
MonitorElement * HFNegEnergy_
MonitorElement * l1GctCenJetsOccEtaPhi_
edm::EDGetTokenT< L1GctEmCandCollection > gctNonIsoEmSourceDataToken_
const unsigned int R3BINS
MonitorElement * l1GctForJetsEtEtaPhi_
virtual void dqmBeginRun(const edm::Run &, const edm::EventSetup &)
MonitorElement * l1GctAllEmOccRankBx_
const unsigned int JETETABINS
MonitorElement * l1GctEtMissPhi_
edm::EDGetTokenT< L1GctEtTotalCollection > l1EtTotalEmulToken_
MonitorElement * l1GctHFRingETDiff_
MonitorElement * book1D(Args &&...args)
Abs< T >::type abs(const T &t)
MonitorElement * l1GctEtMiss_
MonitorElement * l1GctForJetsOccEtaPhi_
edm::EDGetTokenT< L1GctEtTotalCollection > l1EtTotalDataToken_
virtual void beginLuminosityBlock(const edm::LuminosityBlock &, const edm::EventSetup &)
edm::EDGetTokenT< L1GctJetCandCollection > gctForJetsSourceDataToken_
MonitorElement * l1GctHFRing2TowerCountNegEta_
edm::EDGetTokenT< L1GctHFBitCountsCollection > l1HFCountsEmulToken_
edm::EDGetTokenT< L1GctHFRingEtSumsCollection > gctEnergySumsSourceEmulToken_
void setCurrentFolder(const std::string &fullpath)
MonitorElement * book2D(Args &&...args)
MonitorElement * l1GctEtTotalEtHadCorr_
MonitorElement * l1GctHFRingRatioPosEta_
const unsigned int R12BINS
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_
edm::EDGetTokenT< L1GctJetCandCollection > gctCenJetsSourceDataToken_
std::vector< L1CaloRegion > L1CaloRegionCollection
edm::EDGetTokenT< L1GctEtHadCollection > l1EtHadDataToken_
MonitorElement * l1GctNonIsoEmOccEtaPhi_
L1THIonImp(const edm::ParameterSet &ps)
virtual void bookHistograms(DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &) override
MonitorElement * l1GctEtTotalOf_
MonitorElement * l1GctNonIsoEmRank_
const unsigned int METPHIBINS
MonitorElement * l1GctNonIsoEmRankEtaPhi_
MonitorElement * l1GctCenJetsRank_