L1RecoTreeProducer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    L1Trigger/L1TNtuples
// Class:      L1RecoTreeProducer
//
// system include files
#include <memory>

// framework
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/ESHandle.h"

// ROOT output stuff
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"
#include "TH1.h"
#include "TTree.h"
#include "TF1.h"

// EDM formats
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"

//local  data formats
#include "L1Trigger/L1TNtuples/interface/L1AnalysisRecoVertexDataFormat.h"


//
// class declaration
//

class L1RecoTreeProducer : public edm::EDAnalyzer {
public:
  explicit L1RecoTreeProducer(const edm::ParameterSet&);
  ~L1RecoTreeProducer();


private:
  virtual void beginJob(void) ;
  virtual void analyze(const edm::Event&, const edm::EventSetup&);
  virtual void endJob();

public:
  L1Analysis::L1AnalysisRecoVertexDataFormat*              vtxData_;

private:

  // output file
  edm::Service<TFileService> fs_;

  // tree
  TTree * tree_;

  // EDM input tags
  edm::EDGetTokenT<reco::VertexCollection>     vtxToken_;

  unsigned int maxVtx_;

};



L1RecoTreeProducer::L1RecoTreeProducer(const edm::ParameterSet& iConfig)
{
  
  vtxToken_ = consumes<reco::VertexCollection>(iConfig.getUntrackedParameter("vtxToken",edm::InputTag("offlinePrimaryVertices")));

  maxVtx_         = iConfig.getParameter<unsigned int>("maxVtx");

  vtxData_ = new L1Analysis::L1AnalysisRecoVertexDataFormat();

  // set up output
  tree_=fs_->make<TTree>("RecoTree", "RecoTree");
  tree_->Branch("Vertex", "L1Analysis::L1AnalysisRecoVertexDataFormat", &vtxData_, 32000, 3);

}


L1RecoTreeProducer::~L1RecoTreeProducer()
{

  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)

}


//
// member functions
//

// ------------ method called to for each event  ------------
void L1RecoTreeProducer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{

  vtxData_->Reset();

  // get vertices
  edm::Handle<reco::VertexCollection> vertices;
  iEvent.getByToken(vtxToken_, vertices);

  if (vertices.isValid()) {

    for(reco::VertexCollection::const_iterator it=vertices->begin();
    it!=vertices->end() && vtxData_->nVtx < maxVtx_;
    ++it) {
      
      if (!it->isFake()) {
    vtxData_->NDoF.push_back(it->ndof());
    vtxData_->Z.push_back(it->z());
    vtxData_->Rho.push_back(it->position().rho());  
    vtxData_->nVtx++;
      }
      
    }
    tree_->Fill();
  }

}


// ------------ method called once each job just before starting event loop  ------------
void
L1RecoTreeProducer::beginJob(void)
{
}

// ------------ method called once each job just after ending the event loop  ------------
void
L1RecoTreeProducer::endJob() {
}

//define this as a plug-in
DEFINE_FWK_MODULE(L1RecoTreeProducer);