CSJetProducer.h

Go to the documentation of this file.

#ifndef RecoJets_JetProducers_CSJetProducer_h
#define RecoJets_JetProducers_CSJetProducer_h
 
/* *********************************************************
  \class CSJetProducer
 
  \brief Jet producer to produce CMS-style constituent subtracted jets
 
  \author   Marta Verweij
  \modified for granular eta map Chris McGinn
  \version
 
         Notes on implementation:
 
         Constituent subtraction using fastjet contrib package
         The background densities change within the jet as function of eta.
 
 ************************************************************/
 
#include <vector>
#include "RecoJets/JetProducers/plugins/VirtualJetProducer.h"
 
namespace cms {
  class CSJetProducer : public VirtualJetProducer {
  public:
    CSJetProducer(const edm::ParameterSet& ps);
 
    ~CSJetProducer() override {}
    static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
    static void fillDescriptionsFromCSJetProducer(edm::ParameterSetDescription& desc);
 
    void produce(edm::Event& iEvent, const edm::EventSetup& iSetup) override;
 
  protected:
    void runAlgorithm(edm::Event& iEvent, const edm::EventSetup& iSetup) override;
 
    double getModulatedRhoFactor(double phi, double eventPlane2, double eventPlane3, double par1, double par2);
 
    double csRParam_;  
    double csAlpha_;   
 
    bool useModulatedRho_;    
    double minFlowChi2Prob_;  
    double maxFlowChi2Prob_;  
    //input rho and rho_m + eta map
    edm::EDGetTokenT<std::vector<double>> etaToken_;
    edm::EDGetTokenT<std::vector<double>> rhoToken_;
    edm::EDGetTokenT<std::vector<double>> rhomToken_;
    edm::EDGetTokenT<std::vector<double>> rhoFlowFitParamsToken_;
  };
}  // namespace cms
#endif