TPGCheck.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    TPGCheck
// Class:      TPGCheck
//
//
// Original Author:  Muriel Cerutti
//         Created:  Thu Oct 26 10:47:17 CEST 2006
//
 
// system include files
#include <memory>
 
// user include files
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include <string>
#include <vector>
 
#include "TFile.h"
#include "TH1I.h"
 
using namespace edm;
using namespace std;
 
//
// class declaration
//
 
class TPGCheck : public edm::EDAnalyzer {
public:
  explicit TPGCheck(const edm::ParameterSet &);
  ~TPGCheck() override;
 
private:
  void beginJob() override;
  void analyze(const edm::Event &, const edm::EventSetup &) override;
  void endJob() override;
 
  // ----------member data ---------------------------
  TH1I *ecal_et_[2];
  TH1I *ecal_tt_[2];
  TH1I *ecal_fgvb_[2];
 
  TFile *histFile_;
  std::string label_;
  std::string producer_;
  std::vector<std::string> ecal_parts_;
  // fix for consumes
  edm::EDGetTokenT<EcalTrigPrimDigiCollection> ecal_tp_token_;
};
 
//
// constructors and destructor
//
TPGCheck::TPGCheck(const edm::ParameterSet &iConfig) {
  // now do what ever initialization is needed
 
  ecal_parts_.push_back("Barrel");
  ecal_parts_.push_back("Endcap");
 
  histFile_ = new TFile("histos.root", "RECREATE");
  for (unsigned int i = 0; i < 2; ++i) {
    // Energy
    char t[30];
    sprintf(t, "%s_energy", ecal_parts_[i].c_str());
    ecal_et_[i] = new TH1I(t, "Et", 255, 0, 255);
 
    // Trigger Tower flag
    char titleTTF[30];
    sprintf(titleTTF, "%s_ttf", ecal_parts_[i].c_str());
    ecal_tt_[i] = new TH1I(titleTTF, "TTF", 10, 0, 10);
 
    // Fain Grain
    char titleFG[30];
    sprintf(titleFG, "%s_fgvb", ecal_parts_[i].c_str());
    ecal_fgvb_[i] = new TH1I(titleFG, "FGVB", 10, 0, 10);
  }
 
  label_ = iConfig.getParameter<std::string>("Label");
  producer_ = iConfig.getParameter<std::string>("Producer");
  ecal_tp_token_ = consumes<EcalTrigPrimDigiCollection>(edm::InputTag(label_, producer_));
}
 
TPGCheck::~TPGCheck() {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
 
  histFile_->Write();
  histFile_->Close();
}
 
//
// member functions
//
 
// ------------ method called to for each event  ------------
void TPGCheck::analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup) {
  // Get input
  edm::Handle<EcalTrigPrimDigiCollection> tp;
  iEvent.getByToken(ecal_tp_token_, tp);
  for (unsigned int i = 0; i < tp.product()->size(); i++) {
    EcalTriggerPrimitiveDigi d = (*(tp.product()))[i];
    int subdet = d.id().subDet() - 1;
    // for endcap, regroup double TP-s that are generated for the 2 interior
    // rings
    if (subdet == 0) {
      ecal_et_[subdet]->Fill(d.compressedEt());
    } else {
      if (d.id().ietaAbs() == 27 || d.id().ietaAbs() == 28) {
        if (i % 2)
          ecal_et_[subdet]->Fill(d.compressedEt() * 2.);
      } else
        ecal_et_[subdet]->Fill(d.compressedEt());
    }
    ecal_tt_[subdet]->Fill(d.ttFlag());
    ecal_fgvb_[subdet]->Fill(d.fineGrain());
  }
}
 
// ------------ method called once each job just before starting event loop
// ------------
void TPGCheck::beginJob() {}
 
// ------------ method called once each job just after ending the event loop
// ------------
void TPGCheck::endJob() {
  for (unsigned int i = 0; i < 2; ++i) {
    ecal_et_[i]->Write();
    ecal_tt_[i]->Write();
    ecal_fgvb_[i]->Write();
  }
}
 
// define this as a plug-in
DEFINE_FWK_MODULE(TPGCheck);