GenFilters/plugins/PythiaAllDauVFilter.h

0001 #ifndef PYTHIAALLDAUVFILTER_h
0002 #define PYTHIAALLDAUVFILTER_h
0003 // -*- C++ -*-
0004 //
0005 // Package:    PythiaAllDauVFilter
0006 // Class:      PythiaAllDauVFilter
0007 //
0008 /**\class PythiaAllDauVFilter PythiaAllDauVFilter.cc
0009
0010  Description: Filter events using MotherId and ChildrenIds infos
0011         Accepts if event has a specified Mother with only specified daughters and all of the daughters  complies to respective pT and eta Cuts
0012
0013  Implementation:
0014      <Notes on implementation>
0015 */
0016 //
0017 // Original Author:  Daniele Pedrini, PythiaDauVFilter
0018 //                   Aravind T S    ,   Updated and Modified to PythiaAllDauVFilter
0019 //         Created:  Apr 29 2008
0020 //                   Apr 12 2021
0021 //
0022 //
0023
0024 // system include files
0025 #include <memory>
0026
0027 // user include files
0028 #include "FWCore/Framework/interface/Frameworkfwd.h"
0029 #include "FWCore/Framework/interface/global/EDFilter.h"
0030
0031 #include "FWCore/Framework/interface/Event.h"
0032 #include "FWCore/Framework/interface/MakerMacros.h"
0033
0034 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0035
0036 #include "Pythia8/Pythia.h"
0037
0038 //
0039 // class decleration
0040 //
0041 namespace edm {
0042   class HepMCProduct;
0043 }
0044
0045 class PythiaAllDauVFilter : public edm::global::EDFilter<> {
0046 public:
0047   explicit PythiaAllDauVFilter(const edm::ParameterSet&);
0048   ~PythiaAllDauVFilter() override;
0049
0050   bool filter(edm::StreamID, edm::Event&, const edm::EventSetup&) const override;
0051
0052 private:
0053   const int fVerbose;
0054   const edm::EDGetTokenT<edm::HepMCProduct> token_;
0055   std::vector<int> dauIDs, antiDauIDs;
0056   const int particleID;
0057   int antiParticleID;
0058   const int motherID;
0059   const bool chargeconju;
0060   const int ndaughters;
0061   std::vector<double> minptcut;
0062   const double maxptcut;
0063   std::vector<double> minetacut;
0064   std::vector<double> maxetacut;
0065   std::unique_ptr<Pythia8::Pythia> fLookupGen;  // this instance is for accessing particleData information
0066 };
0067 #endif