miscalibExample.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    miscalibExample
// Class:      miscalibExample
//
//
// Original Author:  Lorenzo AGOSTINO
//         Created:  Tue Jul 18 12:17:01 CEST 2006
//
//
 
// system include files
 
// user include files
#include "Calibration/EcalCalibAlgos/interface/miscalibExample.h"
 
//
 
#include "DataFormats/EgammaReco/interface/SuperCluster.h"
#include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include <iostream>
#include <stdexcept>
#include <vector>
 
// class declaration
//
/*
class miscalibExample : public edm::EDAnalyzer {
   public:
      explicit miscalibExample(const edm::ParameterSet&);
      ~miscalibExample();
 
      virtual void analyze(const edm::Event&, const edm::EventSetup&);
      virtual void beginJob();
      virtual void endJob();
   private:
 
 
      // ----------member data ---------------------------
      std::string rootfile_;
      std::string correctedHybridSuperClusterProducer_;
      std::string correctedHybridSuperClusterCollection_;
      std::string BarrelHitsCollection_;
      std::string ecalHitsProducer_ ;
      int read_events;
 
      TH1F* scEnergy;
};
 
*/
miscalibExample::miscalibExample(const edm::ParameterSet& iConfig) {
  rootfile_ = iConfig.getUntrackedParameter<std::string>("rootfile", "ecalSimpleTBanalysis.root");
  correctedHybridSuperClusterProducer_ = iConfig.getParameter<std::string>("correctedHybridSuperClusterProducer");
  correctedHybridSuperClusterCollection_ = iConfig.getParameter<std::string>("correctedHybridSuperClusterCollection");
}
 
miscalibExample::~miscalibExample() {}
 
//========================================================================
void miscalibExample::beginJob() {
  //========================================================================
 
  // Book histograms
  scEnergy = new TH1F("scEnergy", "SuperCluster energy", 100, 20., 80.);
  read_events = 0;
}
 
//========================================================================
void miscalibExample::endJob() {
  //========================================================================
 
  std::cout << "************* STATISTICS **************" << std::endl;
  std::cout << "Read Events: " << read_events << std::endl;
 
 
  TFile f(rootfile_.c_str(), "RECREATE");
 
  scEnergy->Write();
  f.Close();
}
 
//=================================================================================
void miscalibExample::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  //=================================================================================
  using namespace edm;
  using namespace std;
 
  read_events++;
 
  // Get hybrid super clusters after energy correction
 
  Handle<reco::SuperClusterCollection> pCorrectedHybridSuperClusters;
  iEvent.getByLabel(
      correctedHybridSuperClusterProducer_, correctedHybridSuperClusterCollection_, pCorrectedHybridSuperClusters);
  if (!pCorrectedHybridSuperClusters.isValid()) {
    LogError("EgammaSimpleAnalyzer") << "Error! can't get collection with label "
                                     << correctedHybridSuperClusterCollection_.c_str();
  }
  const reco::SuperClusterCollection* correctedHybridSuperClusters = pCorrectedHybridSuperClusters.product();
 
  reco::SuperClusterCollection::const_iterator superClusterIt;
  for (superClusterIt = correctedHybridSuperClusters->begin(); superClusterIt != correctedHybridSuperClusters->end();
       superClusterIt++) {
    scEnergy->Fill(superClusterIt->energy());
  }
}
 
//define this as a plug-in
//DEFINE_FWK_MODULE(miscalibExample);