VBFForwardJetsFilter.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef GENERATORFILTERSVBFFORWARDJETSFILTER_H
#define GENERATORFILTERSVBFFORWARDJETSFILTER_H
 
#include "GeneratorModules/GenFilter.h"
#include "CLHEP/Vector/LorentzVector.h"
#include <vector>
#include "xAODJet/Jet.h"
 
 
class VBFForwardJetsFilter : public GenFilter {
public:
 
  VBFForwardJetsFilter(const std::string & name, ISvcLocator * pSvcLocator);
  virtual StatusCode filterInitialize();
  virtual StatusCode filterEvent();
 
private:
 
  double m_JetMinPt;
  double m_JetMaxEta;
  int m_NJets;
  double m_Jet1MinPt;
  double m_Jet1MaxEta;
  double m_Jet2MinPt;
  double m_Jet2MaxEta;
  bool m_UseOppositeSignEtaJet1Jet2;
  double m_DeltaEtaJJ;
  double m_MassJJ;
  bool m_UseLeadingJJ;
  double m_LGMinPt;
  double m_LGMaxEta;
  double m_DeltaRJLG;
  double m_RatioPtJLG;
 
  std::string m_TruthJetContainerName;
 
  CLHEP::HepLorentzVector sumDaughterNeutrinos(const HepMC::ConstGenParticlePtr& tau) const;
  void removePseudoJets(std::vector<const xAOD::Jet*>& jetList,
                        std::vector<HepMC::ConstGenParticlePtr>& MCTruthPhotonList,
                        std::vector<HepMC::ConstGenParticlePtr>& MCTruthElectronList,
                        std::vector<CLHEP::HepLorentzVector>  & MCTruthTauList);
  double getMinDeltaR(const xAOD::Jet* jet, std::vector<HepMC::ConstGenParticlePtr>& list) const;
  double getMinDeltaR(const xAOD::Jet* jet, std::vector<CLHEP::HepLorentzVector>& list) const;
};
 
#endif