PythiaHepMCFilterGammaGamma.h

Go to the documentation of this file.

#ifndef PYTHIAHEPMCFILTERGAMMAGAMMA_h
#define PYTHIAHEPMCFILTERGAMMAGAMMA_h

//
// Package:    GeneratorInterface/GenFilters
// Class:      PythiaHepMCFilterGammaGamma
// 
// Original Author:  Matteo Sani
//
//

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "GeneratorInterface/Core/interface/BaseHepMCFilter.h"
#include "FWCore/Framework/interface/EDFilter.h"

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

namespace edm {
      class HepMCProduct;
}

class PythiaHepMCFilterGammaGamma : public BaseHepMCFilter {
 public:
  explicit PythiaHepMCFilterGammaGamma(const edm::ParameterSet&);
  ~PythiaHepMCFilterGammaGamma() override;
  
  bool filter(const HepMC::GenEvent* myGenEvent) override;
 private:

  const edm::EDGetTokenT<edm::HepMCProduct> token_;

  //----------
  // filter parameters
  //----------

  const double ptSeedThr, etaSeedThr;

  const double ptGammaThr, etaGammaThr;

  const double ptTkThr, etaTkThr;

  const double ptElThr, etaElThr;

  const double dRTkMax;

  const double dRSeedMax;

  const double dPhiSeedMax, dEtaSeedMax;

  const double dRNarrowCone;

  const double pTMinCandidate1, pTMinCandidate2;

  const double etaMaxCandidate;

  const double invMassMin, invMassMax;

  const double energyCut;

  const int nTkConeMax;

  const int nTkConeSum;

  const bool acceptPrompts;

  const double promptPtThreshold;

};
#endif