EcalPreshowerNoiseDistrib.cc

Go to the documentation of this file.

/*
 * \file EcalPreshowerNoiseDistrib.cc
 *
*/
 
#include <Validation/EcalDigis/interface/EcalPreshowerNoiseDistrib.h>
 
EcalPreshowerNoiseDistrib::EcalPreshowerNoiseDistrib(const edm::ParameterSet& ps)
    : ESdigiCollectionToken_(consumes<ESDigiCollection>(ps.getParameter<edm::InputTag>("ESdigiCollection"))) {
  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
 
  // histos
  meESDigiMultiplicity_ = nullptr;
  for (int ii = 0; ii < 3; ii++) {
    meESDigiADC_[ii] = nullptr;
  }
}
 
void EcalPreshowerNoiseDistrib::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const&, edm::EventSetup const&) {
  Char_t histo[200];
 
  sprintf(histo, "multiplicity");
  meESDigiMultiplicity_ = ibooker.book1D(histo, histo, 1000, 0., 137728);
 
  for (int ii = 0; ii < 3; ii++) {
    sprintf(histo, "esRefHistos%02d", ii);
    meESDigiADC_[ii] = ibooker.book1D(histo, histo, 35, 983.5, 1018.5);
  }
 
  for (int ii = 0; ii < 3; ii++) {
    sprintf(histo, "esRefHistosCorr%02d", ii);
    meESDigiCorr_[ii] = ibooker.book2D(histo, histo, 35, 983.5, 1018.5, 35, 983.5, 1018.5);
  }
 
  meESDigi3D_ = ibooker.book3D("meESDigi3D_", "meESDigi3D_", 35, 983.5, 1018.5, 35, 983.5, 1018.5, 35, 983.5, 1018.5);
}
 
void EcalPreshowerNoiseDistrib::analyze(const edm::Event& e, const edm::EventSetup& c) {
  edm::Handle<ESDigiCollection> EcalDigiES;
 
  e.getByToken(ESdigiCollectionToken_, EcalDigiES);
 
  // retrun if no data
  if (!EcalDigiES.isValid())
    return;
 
  // loop over Digis
  const ESDigiCollection* preshowerDigi = EcalDigiES.product();
 
  std::vector<double> esADCCounts;
  esADCCounts.reserve(ESDataFrame::MAXSAMPLES);
 
  int nDigis = 0;
 
  for (unsigned int digis = 0; digis < EcalDigiES->size(); ++digis) {
    nDigis++;
    ESDataFrame esdf = (*preshowerDigi)[digis];
    int nrSamples = esdf.size();
    for (int sample = 0; sample < nrSamples; ++sample) {
      ESSample mySample = esdf[sample];
      if (meESDigiADC_[sample]) {
        meESDigiADC_[sample]->Fill(mySample.adc());
      }
    }
 
    // to study correlations
    if (meESDigiCorr_[0]) {
      meESDigiCorr_[0]->Fill(esdf[0].adc(), esdf[1].adc());
    }
    if (meESDigiCorr_[1]) {
      meESDigiCorr_[1]->Fill(esdf[0].adc(), esdf[2].adc());
    }
    if (meESDigiCorr_[2]) {
      meESDigiCorr_[2]->Fill(esdf[1].adc(), esdf[2].adc());
    }
 
    // reference histo: sample0, sample1, sample2
    if (meESDigi3D_)
      meESDigi3D_->Fill(esdf[0].adc(), esdf[1].adc(), esdf[2].adc());
  }
 
  if (meESDigiMultiplicity_)
    meESDigiMultiplicity_->Fill(nDigis);
}