PythiaDauVFilterMatchID.h

Go to the documentation of this file.

#ifndef PYTHIADAUVFILTERMATCHID_h
#define PYTHIADAUVFILTERMATCHID_h
// -*- C++ -*-
//
// Package:    PythiaDauVFilterMatchID
// Class:      PythiaDauVFilterMatchID
// 
//
// Original Author:  Daniele Pedrini
//         Created:  Apr 29 2008
//         Adapted by Stephan Wiederkehr
//
//


// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/global/EDFilter.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "Pythia8/Pythia.h"

//
// class decleration
//
namespace edm {
  class HepMCProduct;
}

struct decayTarget {
  int pdgID;
  double minPt,maxPt,minEta,maxEta;

};

class PythiaDauVFilterMatchID : public edm::global::EDFilter<> {
 public:
  explicit PythiaDauVFilterMatchID(const edm::ParameterSet&);
  ~PythiaDauVFilterMatchID() override;
  
  
  bool filter(edm::StreamID, edm::Event&, const edm::EventSetup&) const override;
 private:
  const int fVerbose;  
  const edm::EDGetTokenT<edm::HepMCProduct> token_;
  std::vector<int> dauIDs;
  const int particleID;
  const int motherID;
  const bool chargeconju; 
  const int ndaughters;
  std::vector<double> minptcut;
  const double maxptcut;
  std::vector<double> minetacut;
  std::vector<double> maxetacut;
  std::unique_ptr<Pythia8::Pythia> fLookupGen; // this instance is for accessing particleData information
};
#endif