ESPedestalTask.cc

Go to the documentation of this file.

#include <memory>
#include <fstream>
#include <iostream>

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/ESDetId.h"
#include "DataFormats/EcalDigi/interface/ESDataFrame.h"
#include "DataFormats/EcalRawData/interface/EcalRawDataCollections.h"

#include "DQMServices/Core/interface/DQMStore.h"

#include "DQM/EcalPreshowerMonitorModule/interface/ESPedestalTask.h"

using namespace cms;
using namespace edm;
using namespace std;

ESPedestalTask::ESPedestalTask(const edm::ParameterSet& ps) {
  digitoken_ = consumes<ESDigiCollection>(ps.getParameter<InputTag>("DigiLabel"));
  lookup_ = ps.getUntrackedParameter<FileInPath>("LookupTable");
  outputFile_ = ps.getUntrackedParameter<string>("OutputFile", "");
  prefixME_ = ps.getUntrackedParameter<string>("prefixME", "EcalPreshower");

  for (int i = 0; i < 2; ++i)
    for (int j = 0; j < 2; ++j)
      for (int k = 0; k < 40; ++k)
        for (int l = 0; l < 40; ++l)
          senCount_[i][j][k][l] = -1;

  ievt_ = 0;
}

void ESPedestalTask::bookHistograms(DQMStore::IBooker& iBooker, Run const&, EventSetup const&) {
  int iz, ip, ix, iy, fed, kchip, pace, bundle, fiber, optorx;
  int senZ_[4288], senP_[4288], senX_[4288], senY_[4288];

  // read in look-up table
  ifstream file(lookup_.fullPath().c_str());
  if (!file.is_open())
    throw cms::Exception("FileNotFound") << lookup_.fullPath();

  file >> nLines_;

  for (int i = 0; i < nLines_; ++i) {
    file >> iz >> ip >> ix >> iy >> fed >> kchip >> pace >> bundle >> fiber >> optorx;

    senZ_[i] = iz;
    senP_[i] = ip;
    senX_[i] = ix;
    senY_[i] = iy;

    iz = (senZ_[i] == 1) ? 0 : 1;
    senCount_[iz][senP_[i] - 1][senX_[i] - 1][senY_[i] - 1] = i;
  }

  char hname[300];

  iBooker.setCurrentFolder(prefixME_ + "/ESPedestalTask");

  for (int i = 0; i < nLines_; ++i) {
    for (int is = 0; is < 32; ++is) {
      sprintf(hname, "ADC Z %d P %d X %d Y %d Str %d", senZ_[i], senP_[i], senX_[i], senY_[i], is + 1);
      meADC_[i][is] = iBooker.book1D(hname, hname, 1000, 899.5, 1899.5);
    }
  }
}

void ESPedestalTask::endJob(void) { LogInfo("ESPedestalTask") << "analyzed " << ievt_ << " events"; }

void ESPedestalTask::analyze(const edm::Event& e, const edm::EventSetup& iSetup) {
  ievt_++;
  runNum_ = e.id().run();

  Handle<ESDigiCollection> digis;
  e.getByToken(digitoken_, digis);

  runtype_ = 1;  // Let runtype_ = 1

  // Digis
  int zside, plane, ix, iy, strip, iz;
  if (digis.isValid()) {
    for (ESDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr) {
      ESDataFrame dataframe = (*digiItr);
      ESDetId id = dataframe.id();

      zside = id.zside();
      plane = id.plane();
      ix = id.six();
      iy = id.siy();
      strip = id.strip();
      iz = (zside == 1) ? 0 : 1;

      if (meADC_[senCount_[iz][plane - 1][ix - 1][iy - 1]][strip - 1]) {
        if (runtype_ == 1) {
          meADC_[senCount_[iz][plane - 1][ix - 1][iy - 1]][strip - 1]->Fill(dataframe.sample(0).adc());
          meADC_[senCount_[iz][plane - 1][ix - 1][iy - 1]][strip - 1]->Fill(dataframe.sample(1).adc());
          meADC_[senCount_[iz][plane - 1][ix - 1][iy - 1]][strip - 1]->Fill(dataframe.sample(2).adc());
        } else if (runtype_ == 3) {
          meADC_[senCount_[iz][plane - 1][ix - 1][iy - 1]][strip - 1]->Fill(dataframe.sample(1).adc());
        }
      }
    }
  }
}

DEFINE_FWK_MODULE(ESPedestalTask);