FastjetJetProducer.h

Go to the documentation of this file.

#ifndef RecoJets_JetProducers_plugins_FastjetJetProducer_h
#define RecoJets_JetProducers_plugins_FastjetJetProducer_h

#include "RecoJets/JetProducers/interface/JetSpecific.h"

#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/RecoCandidate/interface/RecoChargedCandidate.h"

#include "RecoJets/JetProducers/plugins/VirtualJetProducer.h"

#include <fastjet/tools/Transformer.hh>

// a function returning
//   min(Rmax, deltaR_factor * deltaR(j1,j2))
// where j1 and j2 are the 2 subjets of j
// if the jet does not have exactly 2 pieces, Rmax is used.
class DynamicRfilt : public fastjet::FunctionOfPseudoJet<double>{
public:
  // default ctor
  DynamicRfilt(double Rmax, double deltaR_factor) : _Rmax(Rmax), _deltaR_factor(deltaR_factor){}

  // action of the function
  double result(const fastjet::PseudoJet &j) const{
    if (! j.has_pieces()) return _Rmax;

    std::vector<fastjet::PseudoJet> pieces = j.pieces();
    if (! pieces.size()==2) return _Rmax;

    double deltaR = pieces[0].delta_R(pieces[1]);
    return std::min(_Rmax, _deltaR_factor * deltaR);
  }

private:
  double _Rmax, _deltaR_factor;
};



class FastjetJetProducer : public VirtualJetProducer
{

public:
  // typedefs
  typedef fastjet::Transformer         transformer;
  typedef std::unique_ptr<transformer> transformer_ptr;
  typedef std::vector<transformer_ptr> transformer_coll;
  //
  // construction/destruction
  //
  explicit FastjetJetProducer(const edm::ParameterSet& iConfig);
  virtual ~FastjetJetProducer();

  virtual void produce( edm::Event & iEvent, const edm::EventSetup & iSetup );

  // typedefs
  typedef boost::shared_ptr<DynamicRfilt>  DynamicRfiltPtr;

protected:

  //
  // member functions
  //

  virtual void produceTrackJets( edm::Event & iEvent, const edm::EventSetup & iSetup );
  virtual void runAlgorithm( edm::Event& iEvent, const edm::EventSetup& iSetup );

 private:

  // trackjet clustering parameters
  bool useOnlyVertexTracks_;
  bool useOnlyOnePV_;
  float dzTrVtxMax_;
  float dxyTrVtxMax_;
  int minVtxNdof_;
  float maxVtxZ_;

  // jet trimming parameters
  bool useMassDropTagger_;    
  bool useFiltering_;         
  bool useDynamicFiltering_;  
  bool useTrimming_;          
  bool usePruning_;           
  bool useCMSBoostedTauSeedingAlgorithm_; 
  bool useKtPruning_;         
  bool useConstituentSubtraction_; 
  bool useSoftDrop_;          
  double muCut_;              
  double yCut_;               
  double rFilt_;              
  DynamicRfiltPtr rFiltDynamic_; 
  double rFiltFactor_;        
  int    nFilt_;              
  double trimPtFracMin_;      
  double zCut_;               
  double RcutFactor_;         
  double csRho_EtaMax_;       
  double csRParam_;           
  double beta_;               
  double R0_;                 


  double subjetPtMin_;        
  double muMin_;              
  double muMax_;              
  double yMin_;               
  double yMax_;               
  double dRMin_;              
  double dRMax_;              
  int    maxDepth_;           


  // tokens for the data access
  edm::EDGetTokenT<edm::View<reco::RecoChargedRefCandidate> > input_chrefcand_token_;
    
};


#endif