#include <CSJetProducer.h>

Inheritance diagram for cms::CSJetProducer:

Public Member Functions
	CSJetProducer (const edm::ParameterSet &ps)

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

	~CSJetProducer () override

Public Member Functions inherited from VirtualJetProducer
std::string	jetType () const

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

	VirtualJetProducer (const edm::ParameterSet &iConfig)

	~VirtualJetProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

static void	fillDescriptionsFromCSJetProducer (edm::ParameterSetDescription &desc)

Static Public Member Functions inherited from VirtualJetProducer
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

static void	fillDescriptionsFromVirtualJetProducer (edm::ParameterSetDescription &desc)

Protected Member Functions
void	runAlgorithm (edm::Event &iEvent, const edm::EventSetup &iSetup) override

Protected Member Functions inherited from VirtualJetProducer
virtual void	addHTTTopJetTagInfoCollection (edm::Event &iEvent, const edm::EventSetup &iSetup, edm::OrphanHandle< reco::BasicJetCollection > &oh)

virtual void	copyConstituents (const std::vector< fastjet::PseudoJet > &fjConstituents, reco::Jet *jet)

virtual std::vector< reco::CandidatePtr >	getConstituents (const std::vector< fastjet::PseudoJet > &fjConstituents)

virtual void	inputTowers ()

virtual bool	isAnomalousTower (reco::CandidatePtr input)

bool	makeBasicJet (const JetType::Type &fTag)

bool	makeCaloJet (const JetType::Type &fTag)

bool	makeGenJet (const JetType::Type &fTag)

bool	makePFClusterJet (const JetType::Type &fTag)

bool	makePFJet (const JetType::Type &fTag)

virtual void	makeProduces (std::string s, std::string tag="")

bool	makeTrackJet (const JetType::Type &fTag)

void	offsetCorrectJets (std::vector< fastjet::PseudoJet > &orphanInput)

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

template<typename T >
void	writeCompoundJets (edm::Event &iEvent, edm::EventSetup const &iSetup)
	function template to write out the outputs More...

template<typename T >
void	writeJets (edm::Event &iEvent, edm::EventSetup const &iSetup)

template<typename T >
void	writeJetsWithConstituents (edm::Event &iEvent, edm::EventSetup const &iSetup)
	function template to write out the outputs More...

Protected Attributes
double	csAlpha_
	for constituent subtraction : R parameter More...

double	csRParam_

edm::EDGetTokenT< std::vector< double > >	etaToken_
	for HI constituent subtraction : alpha (power of pt in metric) More...

edm::EDGetTokenT< std::vector< double > >	rhomToken_

edm::EDGetTokenT< std::vector< double > >	rhoToken_

Protected Attributes inherited from VirtualJetProducer
int	activeAreaRepeats_

bool	doAreaDiskApprox_

bool	doAreaFastjet_

bool	doFastJetNonUniform_

bool	doPUOffsetCorr_

bool	doPVCorrection_

bool	doRhoFastjet_

ActiveAreaSpecPtr	fjActiveArea_

AreaDefinitionPtr	fjAreaDefinition_

ClusterSequencePtr	fjClusterSeq_

std::vector< fastjet::PseudoJet >	fjInputs_

JetDefPtr	fjJetDefinition_

std::vector< fastjet::PseudoJet >	fjJets_

PluginPtr	fjPlugin_

SelectorPtr	fjSelector_

bool	fromHTTTopJetProducer_ = false

double	ghostArea_

double	ghostEtaMax_

edm::EDGetTokenT< reco::VertexCollection >	input_vertex_token_

double	inputEMin_

double	inputEtMin_

std::vector< edm::Ptr< reco::Candidate > >	inputs_

std::string	jetAlgorithm_

std::string	jetCollInstanceName_

double	jetPtMin_

std::string	jetType_

JetType::Type	jetTypeE

unsigned int	maxInputs_

unsigned int	minSeed_

std::string	moduleLabel_

unsigned int	nExclude_

std::vector< double >	puCenters_

std::string	puSubtractorName_

double	puWidth_

bool	restrictInputs_

double	rhoEtaMax_

double	rParam_

edm::InputTag	src_

edm::InputTag	srcPVs_

boost::shared_ptr< PileUpSubtractor >	subtractor_

bool	useDeterministicSeed_

bool	useExplicitGhosts_

int	verbosity_

reco::Particle::Point	vertex_

double	voronoiRfact_

bool	writeCompound_

bool	writeJetsWithConst_

Additional Inherited Members
Public Types inherited from VirtualJetProducer
typedef boost::shared_ptr< fastjet::GhostedAreaSpec >	ActiveAreaSpecPtr

typedef boost::shared_ptr< fastjet::AreaDefinition >	AreaDefinitionPtr

typedef boost::shared_ptr< fastjet::ClusterSequence >	ClusterSequencePtr

typedef boost::shared_ptr< fastjet::JetDefinition >	JetDefPtr

typedef boost::shared_ptr< fastjet::JetDefinition::Plugin >	PluginPtr

typedef boost::shared_ptr< fastjet::Selector >	SelectorPtr

Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 24 of file CSJetProducer.h.

Constructor & Destructor Documentation

CSJetProducer::CSJetProducer ( const edm::ParameterSet & ps )

Definition at line 13 of file CSJetProducer.cc.

References csAlpha_, csRParam_, etaToken_, edm::ParameterSet::getParameter(), rhomToken_, and rhoToken_.

                                                        :
   VirtualJetProducer( conf ),
   csRParam_(-1.0),
   csAlpha_(0.)
 {
   //get eta range, rho and rhom map
   etaToken_ = consumes<std::vector<double>>(conf.getParameter<edm::InputTag>( "etaMap" ));
   rhoToken_ = consumes<std::vector<double>>(conf.getParameter<edm::InputTag>( "rho" ));
   rhomToken_ = consumes<std::vector<double>>(conf.getParameter<edm::InputTag>( "rhom" ));
   csRParam_ = conf.getParameter<double>("csRParam");
   csAlpha_ = conf.getParameter<double>("csAlpha");
 }

cms::CSJetProducer::~CSJetProducer ( )

inlineoverride

Definition at line 30 of file CSJetProducer.h.

References fillDescriptions(), fillDescriptionsFromCSJetProducer(), iEvent, produce(), and runAlgorithm().

30 {}

Member Function Documentation

void CSJetProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 126 of file CSJetProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), fillDescriptionsFromCSJetProducer(), and VirtualJetProducer::fillDescriptionsFromVirtualJetProducer().

Referenced by ~CSJetProducer().

                                                                                {
 
   edm::ParameterSetDescription descCSJetProducer;
   fillDescriptionsFromCSJetProducer(descCSJetProducer);
   descCSJetProducer.add<string>("jetCollInstanceName", ""    );
   VirtualJetProducer::fillDescriptionsFromVirtualJetProducer(descCSJetProducer);
   descCSJetProducer.add<bool> ("sumRecHits", false);
   
   descriptions.add("CSJetProducer",descCSJetProducer);
   
 }

void CSJetProducer::fillDescriptionsFromCSJetProducer ( edm::ParameterSetDescription & desc )

static

Definition at line 141 of file CSJetProducer.cc.

References edm::ParameterSetDescription::add(), and DEFINE_FWK_MODULE.

Referenced by fillDescriptions(), and ~CSJetProducer().

                                                                                       {
 
   desc.add<double>("csRParam",-1.);
   desc.add<double>("csAlpha",2.);
 
   desc.add<edm::InputTag>("etaMap",edm::InputTag("hiFJRhoProducer","mapEtaEdges") );
   desc.add<edm::InputTag>("rho",edm::InputTag("hiFJRhoProducer","mapToRho") );
   desc.add<edm::InputTag>("rhom",edm::InputTag("hiFJRhoProducer","mapToRhoM") );
 
 }

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

override

Definition at line 26 of file CSJetProducer.cc.

References VirtualJetProducer::fjClusterSeq_, and VirtualJetProducer::produce().

Referenced by ~CSJetProducer().

 {
   // use the default production from one collection
   VirtualJetProducer::produce( iEvent, iSetup );
   //use runAlgorithm of this class
 
   //Delete allocated memory. It is allocated every time runAlgorithm is called
   fjClusterSeq_.reset();
 }

void CSJetProducer::runAlgorithm	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprotectedvirtual

Implements VirtualJetProducer.

Definition at line 37 of file CSJetProducer.cc.

Referenced by ~CSJetProducer().

 {
   // run algorithm
   if ( !doAreaFastjet_ && !doRhoFastjet_) {
     fjClusterSeq_ = ClusterSequencePtr( new fastjet::ClusterSequence( fjInputs_, *fjJetDefinition_ ) );
   } else if (voronoiRfact_ <= 0) {
     fjClusterSeq_ = ClusterSequencePtr( new fastjet::ClusterSequenceArea( fjInputs_, *fjJetDefinition_ , *fjAreaDefinition_ ) );
   } else {
     fjClusterSeq_ = ClusterSequencePtr( new fastjet::ClusterSequenceVoronoiArea( fjInputs_, *fjJetDefinition_ , fastjet::VoronoiAreaSpec(voronoiRfact_) ) );
   }
 
   fjJets_.clear();
   std::vector<fastjet::PseudoJet> tempJets = fastjet::sorted_by_pt(fjClusterSeq_->inclusive_jets(jetPtMin_));
 
   //Get local rho and rhom map
   edm::Handle<std::vector<double>> etaRanges;
   edm::Handle<std::vector<double>> rhoRanges;
   edm::Handle<std::vector<double>> rhomRanges;
   
   iEvent.getByToken(etaToken_, etaRanges);
   iEvent.getByToken(rhoToken_, rhoRanges);
   iEvent.getByToken(rhomToken_, rhomRanges);
 
   //Check if size of eta and background density vectors is the same
   unsigned int bkgVecSize = etaRanges->size();
   if(bkgVecSize<1) { throw cms::Exception("WrongBkgInput") << "Producer needs vector with background estimates\n"; }
   if(bkgVecSize != (rhoRanges->size()+1) || bkgVecSize != (rhomRanges->size()+1) ) {
     throw cms::Exception("WrongBkgInput") << "Size of etaRanges (" << bkgVecSize << ") and rhoRanges (" << rhoRanges->size() << ") and/or rhomRanges (" << rhomRanges->size() << ") vectors inconsistent\n";
   }
 
   
   //Allow the background densities to change within the jet
 
   for(fastjet::PseudoJet& ijet : tempJets ) {
   
     //----------------------------------------------------------------------
     // sift ghosts and particles in the input jet
     std::vector<fastjet::PseudoJet> particles, ghosts;
     fastjet::SelectorIsPureGhost().sift(ijet.constituents(), ghosts, particles);
     unsigned long nParticles=particles.size();
     if(nParticles==0) continue; //don't subtract ghost jets
     
     //assign rho and rhom to ghosts according to local eta-dependent map
     double rho  = 1e-6;
     double rhom = 1e-6;
     for(fastjet::PseudoJet& ighost : ghosts) {
       
       if(ighost.eta()<=etaRanges->at(0) || bkgVecSize==1) {
         rho = rhoRanges->at(0);
         rhom = rhomRanges->at(0);
       } else if(ighost.eta()>=etaRanges->at(bkgVecSize-1)) {
         rho = rhoRanges->at(bkgVecSize-2);
         rhom = rhomRanges->at(bkgVecSize-2);
       } else {
         for(unsigned int ie = 0; ie<(bkgVecSize-1); ie++) {
           if(ighost.eta()>=etaRanges->at(ie) && ighost.eta()<etaRanges->at(ie+1)) {
             rho = rhoRanges->at(ie);
             rhom = rhomRanges->at(ie);
             break;
           }
         }
       }
       double pt = rho*ighost.area();
       double mass_squared=std::pow(rhom*ighost.area()+pt,2)-std::pow(pt,2);
       double mass=0;
       if (mass_squared>0) mass=sqrt(mass_squared);
       ighost.reset_momentum_PtYPhiM(pt,ighost.rap(),ighost.phi(),mass);
     }
 
     //----------------------------------------------------------------------
     //from here use official fastjet contrib package
     
     fastjet::contrib::ConstituentSubtractor subtractor;
     subtractor.set_distance_type(fastjet::contrib::ConstituentSubtractor::deltaR); // distance in eta-phi plane
     subtractor.set_max_distance(csRParam_); // free parameter for the maximal allowed distance between particle i and ghost k
     subtractor.set_alpha(csAlpha_);  // free parameter for the distance measure (the exponent of particle pt). Note that in older versions of the package alpha was multiplied by two but in newer versions this is not the case anymore
     subtractor.set_do_mass_subtraction(true);
     subtractor.set_remove_all_zero_pt_particles(true);
 
     std::vector<fastjet::PseudoJet> subtracted_particles = subtractor.do_subtraction(particles,ghosts);
 
     //Create subtracted jets
     fastjet::PseudoJet subtracted_jet=join(subtracted_particles);
     if(subtracted_jet.perp()>0.)
       fjJets_.push_back( subtracted_jet );
   }
   fjJets_ = fastjet::sorted_by_pt(fjJets_); 
 }