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 "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.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 override {
    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);
  ~FastjetJetProducer() override;
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
  static void fillDescriptionsFromFastJetProducer(edm::ParameterSetDescription& desc);
 
  void produce(edm::Event& iEvent, const edm::EventSetup& iSetup) override;
 
  // typedefs
  typedef std::shared_ptr<DynamicRfilt> DynamicRfiltPtr;
 
protected:
  //
  // member functions
  //
 
  virtual void produceTrackJets(edm::Event& iEvent, const edm::EventSetup& iSetup);
  void runAlgorithm(edm::Event& iEvent, const edm::EventSetup& iSetup) override;
 
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_;                
  bool correctShape_;               
  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 gridMaxRapidity_;  
  double gridSpacing_;      
 
  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