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

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

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