EcalTBHodoscopeRawInfoDumper.cc

Go to the documentation of this file.

#include "RecoTBCalo/EcalTBHodoscopeReconstructor/interface/EcalTBHodoscopeRawInfoDumper.h"
#include "TBDataFormats/EcalTBObjects/interface/EcalTBHodoscopeRawInfo.h"
#include "DataFormats/Common/interface/EDCollection.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"

#include <TFile.h>
#include <FWCore/Framework/interface/MakerMacros.h>
#include "FWCore/MessageLogger/interface/MessageLogger.h"

EcalTBHodoscopeRawInfoDumper::EcalTBHodoscopeRawInfoDumper(edm::ParameterSet const& ps) {
  rawInfoCollection_ = ps.getParameter<std::string>("rawInfoCollection");
  rawInfoProducer_ = ps.getParameter<std::string>("rawInfoProducer");
  rootfile_ = ps.getUntrackedParameter<std::string>("rootfile", "ecalHodoscopeRawInfoPlots.root");
}

EcalTBHodoscopeRawInfoDumper::~EcalTBHodoscopeRawInfoDumper() {}

//========================================================================
void EcalTBHodoscopeRawInfoDumper::beginJob() {
  //========================================================================

  char histoName[100];
  char histoTitle[100];

  for (int i = 0; i < 4; i++) {
    sprintf(histoName, "h_numberOfFiredHits_%d", i);
    sprintf(histoTitle, "NumberOfFiredHits Plane %d", i);
    h_numberOfFiredHits_[i] = new TH1F(histoName, histoTitle, 10, 0., 10.);
  }

  for (int i = 0; i < 4; i++) {
    sprintf(histoName, "h_firedHits_%d", i);
    sprintf(histoTitle, "firedHits Plane %d", i);
    h_firedHits_[i] = new TH1F(histoName, histoTitle, 64, -0.5, 63.5);
  }
}

//========================================================================
void EcalTBHodoscopeRawInfoDumper::endJob() {
  //========================================================================

  TFile f(rootfile_.c_str(), "RECREATE");

  for (int i = 0; i < 4; i++)
    h_numberOfFiredHits_[i]->Write();

  for (int i = 0; i < 4; i++)
    h_firedHits_[i]->Write();

  f.Close();
}

void EcalTBHodoscopeRawInfoDumper::analyze(const edm::Event& e, const edm::EventSetup& es) {
  // Get input
  edm::Handle<EcalTBHodoscopeRawInfo> ecalRawHodoscope;
  const EcalTBHodoscopeRawInfo* hodoscopeRawInfo = nullptr;
  //evt.getByLabel( digiProducer_, digiCollection_, pDigis);
  e.getByLabel(rawInfoProducer_, ecalRawHodoscope);
  if (!ecalRawHodoscope.isValid()) {
    edm::LogError("EcalTBHodoscopeRecInfoError") << "Error! can't get the product " << rawInfoCollection_.c_str();
  } else {
    hodoscopeRawInfo = ecalRawHodoscope.product();
  }

  if (hodoscopeRawInfo)
    for (int i = 0; i < 4; i++) {
      std::ostringstream str;
      str << " Hits ";
      std::vector<int> firedHits;
      h_numberOfFiredHits_[i]->Fill((*hodoscopeRawInfo)[i].numberOfFiredHits());
      for (int j = 0; j < 64; j++)
        if ((*hodoscopeRawInfo)[i][j]) {
          h_firedHits_[i]->Fill(j);
          firedHits.push_back(j);
          str << j << " ";
        }
      LogDebug("EcalTBHodoscope") << "Looking plane " << i << " number of hits "
                                  << (*hodoscopeRawInfo)[i].numberOfFiredHits() << str.str();
    }
  // Create empty output
}