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);