#include <FastjetJetProducer.h>

Inheritance diagram for FastjetJetProducer:

Public Types
typedef boost::shared_ptr < DynamicRfilt >	DynamicRfiltPtr

typedef fastjet::Transformer	transformer

typedef std::vector < transformer_ptr >	transformer_coll

typedef std::unique_ptr < transformer >	transformer_ptr

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::RangeDefinition >	RangeDefPtr

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	FastjetJetProducer (const edm::ParameterSet &iConfig)

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

virtual	~FastjetJetProducer ()

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

	VirtualJetProducer (const edm::ParameterSet &iConfig)

virtual	~VirtualJetProducer ()

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
virtual void	produceTrackJets (edm::Event &iEvent, const edm::EventSetup &iSetup)

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

Protected Member Functions inherited from VirtualJetProducer
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)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Attributes
double	beta_
	for constituent subtraction : R parameter for KT alg in jet median background estimator More...

bool	correctShape_
	Soft drop. More...

double	csRho_EtaMax_
	for pruning: constituent dR * pt/2m < rcut_factor More...

double	csRParam_
	for constituent subtraction : maximum rapidity for ghosts More...

double	dRMax_
	for CMSBoostedTauSeedingAlgorithm : min dR More...

double	dRMin_
	for CMSBoostedTauSeedingAlgorithm : max asymmetry More...

float	dxyTrVtxMax_

float	dzTrVtxMax_

double	gridMaxRapidity_
	for soft drop : R0 (angular distance normalization - should be set to jet radius in most cases) More...

double	gridSpacing_
	for shape subtraction, get the fixed-grid rho More...

edm::EDGetTokenT< edm::View < reco::RecoChargedRefCandidate > >	input_chrefcand_token_
	for CMSBoostedTauSeedingAlgorithm : max depth for descending into clustering sequence More...

int	maxDepth_
	for CMSBoostedTauSeedingAlgorithm : max dR More...

float	maxVtxZ_

int	minVtxNdof_

double	muCut_
	Correct the shape of the jets. More...

double	muMax_
	for CMSBoostedTauSeedingAlgorithm : min mass-drop More...

double	muMin_
	for CMSBoostedTauSeedingAlgorithm : subjet pt min More...

int	nFilt_
	for dynamic filtering radius (as in arXiv:0802.2470) More...

double	R0_
	for soft drop : beta (angular exponent) More...

double	RcutFactor_
	for pruning OR soft drop: constituent minimum pt fraction of parent cluster More...

double	rFilt_
	for mass-drop tagging, symmetry cut: min(pt1^2,pt2^2) * dR(1,2) / mjet > ycut More...

DynamicRfiltPtr	rFiltDynamic_
	for filtering, trimming: dR scale of sub-clustering More...

double	rFiltFactor_
	for dynamic filtering radius (as in arXiv:0802.2470) More...

double	subjetPtMin_
	for shape subtraction, get the grid spacing More...

double	trimPtFracMin_
	for filtering, pruning: number of subjets expected More...

bool	useCMSBoostedTauSeedingAlgorithm_
	Jet pruning technique. More...

bool	useConstituentSubtraction_
	Use Kt clustering algorithm for pruning (default is Cambridge/Aachen) More...

bool	useDynamicFiltering_
	Jet filtering technique. More...

bool	useFiltering_
	Mass-drop tagging for boosted Higgs. More...

bool	useKtPruning_
	algorithm for seeding reconstruction of boosted Taus (similar to mass-drop tagging) More...

bool	useMassDropTagger_

bool	useOnlyOnePV_

bool	useOnlyVertexTracks_

bool	usePruning_
	Jet trimming technique. More...

bool	useSoftDrop_
	constituent subtraction technique More...

bool	useTrimming_
	Use dynamic filtering radius (as in arXiv:0802.2470) More...

double	yCut_
	for mass-drop tagging, m0/mjet (m0 = mass of highest mass subjet) More...

double	yMax_
	for CMSBoostedTauSeedingAlgorithm : min asymmetry More...

double	yMin_
	for CMSBoostedTauSeedingAlgorithm : max mass-drop More...

double	zCut_
	for trimming: constituent minimum pt fraction of full jet More...

Additional Inherited Members
Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from VirtualJetProducer
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_

RangeDefPtr	fjRangeDef_

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	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_

Detailed Description

Definition at line 39 of file FastjetJetProducer.h.

Member Typedef Documentation

typedef boost::shared_ptr<DynamicRfilt> FastjetJetProducer::DynamicRfiltPtr

Definition at line 56 of file FastjetJetProducer.h.

typedef fastjet::Transformer FastjetJetProducer::transformer

Definition at line 44 of file FastjetJetProducer.h.

typedef std::vector<transformer_ptr> FastjetJetProducer::transformer_coll

Definition at line 46 of file FastjetJetProducer.h.

typedef std::unique_ptr<transformer> FastjetJetProducer::transformer_ptr

Definition at line 45 of file FastjetJetProducer.h.

Constructor & Destructor Documentation

FastjetJetProducer::FastjetJetProducer ( const edm::ParameterSet & iConfig )

explicit

Can't use these together

Definition at line 59 of file FastjetJetProducer.cc.

References beta_, correctShape_, csRho_EtaMax_, csRParam_, dRMax_, dRMin_, dxyTrVtxMax_, dzTrVtxMax_, Exception, edm::ParameterSet::exists(), VirtualJetProducer::fjAreaDefinition_, edm::ParameterSet::getParameter(), gridMaxRapidity_, gridSpacing_, input_chrefcand_token_, maxDepth_, maxVtxZ_, minVtxNdof_, muCut_, muMax_, muMin_, nFilt_, R0_, RcutFactor_, rFilt_, rFiltDynamic_, rFiltFactor_, VirtualJetProducer::src_, subjetPtMin_, trimPtFracMin_, useCMSBoostedTauSeedingAlgorithm_, useConstituentSubtraction_, useDynamicFiltering_, VirtualJetProducer::useExplicitGhosts_, useFiltering_, useKtPruning_, useMassDropTagger_, useOnlyOnePV_, useOnlyVertexTracks_, usePruning_, useSoftDrop_, useTrimming_, yCut_, yMax_, yMin_, and zCut_.

   : VirtualJetProducer( iConfig ),
     useMassDropTagger_(false),
     useFiltering_(false),
     useDynamicFiltering_(false),
     useTrimming_(false),
     usePruning_(false),
     useCMSBoostedTauSeedingAlgorithm_(false),
     useKtPruning_(false),
     useConstituentSubtraction_(false),
     useSoftDrop_(false),
     correctShape_(false),
     muCut_(-1.0),
     yCut_(-1.0),
     rFilt_(-1.0),
     rFiltFactor_(-1.0),
     nFilt_(-1),
     trimPtFracMin_(-1.0),
     zCut_(-1.0),
     RcutFactor_(-1.0),
     csRho_EtaMax_(-1.0),
     csRParam_(-1.0),
     beta_(-1.0),
     R0_(-1.0),
     gridMaxRapidity_(-1.0), // For fixed-grid rho
     gridSpacing_(-1.0)  // For fixed-grid rho
 {
 
   if ( iConfig.exists("UseOnlyVertexTracks") )
     useOnlyVertexTracks_ = iConfig.getParameter<bool>("UseOnlyVertexTracks");
   else 
     useOnlyVertexTracks_ = false;
   
   if ( iConfig.exists("UseOnlyOnePV") )
     useOnlyOnePV_        = iConfig.getParameter<bool>("UseOnlyOnePV");
   else
     useOnlyOnePV_ = false;
 
   if ( iConfig.exists("DzTrVtxMax") )
     dzTrVtxMax_          = iConfig.getParameter<double>("DzTrVtxMax");
   else
     dzTrVtxMax_ = 999999.;
   if ( iConfig.exists("DxyTrVtxMax") )
     dxyTrVtxMax_          = iConfig.getParameter<double>("DxyTrVtxMax");
   else
     dxyTrVtxMax_ = 999999.;
   if ( iConfig.exists("MinVtxNdof") )
     minVtxNdof_ = iConfig.getParameter<int>("MinVtxNdof");
   else
     minVtxNdof_ = 5;
   if ( iConfig.exists("MaxVtxZ") )
     maxVtxZ_ = iConfig.getParameter<double>("MaxVtxZ");
   else
     maxVtxZ_ = 15;
 
 
   if ( iConfig.exists("useFiltering") ||
        iConfig.exists("useTrimming") ||
        iConfig.exists("usePruning") ||
        iConfig.exists("useMassDropTagger") ||
        iConfig.exists("useCMSBoostedTauSeedingAlgorithm") ||
        iConfig.exists("useConstituentSubtraction") ||
        iConfig.exists("useSoftDrop")
        ) {
     useMassDropTagger_=false;
     useFiltering_=false;
     useDynamicFiltering_=false;
     useTrimming_=false;
     usePruning_=false;
     useCMSBoostedTauSeedingAlgorithm_=false;
     useKtPruning_=false;
     useConstituentSubtraction_=false;
     useSoftDrop_ = false;
     rFilt_=-1.0;
     rFiltFactor_=-1.0;
     nFilt_=-1;
     trimPtFracMin_=-1.0;
     zCut_=-1.0;
     RcutFactor_=-1.0;
     muCut_=-1.0;
     yCut_=-1.0;
     subjetPtMin_ = -1.0;
     muMin_ = -1.0;
     muMax_ = -1.0;
     yMin_ = -1.0;
     yMax_ = -1.0;
     dRMin_ = -1.0;
     dRMax_ = -1.0;
     maxDepth_ = -1;
     csRho_EtaMax_ = -1.0;
     csRParam_ = -1.0;
     beta_ = -1.0;
     R0_ = -1.0;
     useExplicitGhosts_ = true;
 
     if ( iConfig.exists("useMassDropTagger") ) {
       useMassDropTagger_ = iConfig.getParameter<bool>("useMassDropTagger");
       muCut_ = iConfig.getParameter<double>("muCut");
       yCut_ = iConfig.getParameter<double>("yCut");
     }
 
     if ( iConfig.exists("useFiltering") ) {
       useFiltering_ = iConfig.getParameter<bool>("useFiltering");
       rFilt_ = iConfig.getParameter<double>("rFilt");
       nFilt_ = iConfig.getParameter<int>("nFilt");
       if ( iConfig.exists("useDynamicFiltering") ) {
         useDynamicFiltering_ = iConfig.getParameter<bool>("useDynamicFiltering");
         rFiltFactor_ = iConfig.getParameter<double>("rFiltFactor");
         if ( useDynamicFiltering_ )
           rFiltDynamic_ = DynamicRfiltPtr(new DynamicRfilt(rFilt_, rFiltFactor_));
       }
     }
   
     if ( iConfig.exists("useTrimming") ) {
       useTrimming_ = iConfig.getParameter<bool>("useTrimming");
       rFilt_ = iConfig.getParameter<double>("rFilt");
       trimPtFracMin_ = iConfig.getParameter<double>("trimPtFracMin");
     }
 
     if ( iConfig.exists("usePruning") ) {
       usePruning_ = iConfig.getParameter<bool>("usePruning");
       zCut_ = iConfig.getParameter<double>("zcut");
       RcutFactor_ = iConfig.getParameter<double>("rcut_factor");
       nFilt_ = iConfig.getParameter<int>("nFilt");
       if ( iConfig.exists("useKtPruning") )
         useKtPruning_ = iConfig.getParameter<bool>("useKtPruning");
     }
 
     if ( iConfig.exists("useCMSBoostedTauSeedingAlgorithm") ) {
       useCMSBoostedTauSeedingAlgorithm_ = iConfig.getParameter<bool>("useCMSBoostedTauSeedingAlgorithm");
       subjetPtMin_ = iConfig.getParameter<double>("subjetPtMin");
       muMin_ = iConfig.getParameter<double>("muMin");
       muMax_ = iConfig.getParameter<double>("muMax");
       yMin_ = iConfig.getParameter<double>("yMin");
       yMax_ = iConfig.getParameter<double>("yMax");
       dRMin_ = iConfig.getParameter<double>("dRMin");
       dRMax_ = iConfig.getParameter<double>("dRMax");
       maxDepth_ = iConfig.getParameter<int>("maxDepth");
     }
 
     if ( iConfig.exists("useConstituentSubtraction") ) {
 
       if ( fjAreaDefinition_.get() == 0 ) {
     throw cms::Exception("AreaMustBeSet") << "Logic error. The area definition must be set if you use constituent subtraction." << std::endl;
       }
 
       useConstituentSubtraction_ = iConfig.getParameter<bool>("useConstituentSubtraction");
       csRho_EtaMax_ = iConfig.getParameter<double>("csRho_EtaMax");
       csRParam_ = iConfig.getParameter<double>("csRParam");
     }
 
     if ( iConfig.exists("useSoftDrop") ) {
       if ( usePruning_ ) {   
     throw cms::Exception("PruningAndSoftDrop") << "Logic error. Soft drop is a generalized pruning, do not run them together." << std::endl;
       }
       useSoftDrop_ = iConfig.getParameter<bool>("useSoftDrop");
       zCut_ = iConfig.getParameter<double>("zcut");
       beta_ = iConfig.getParameter<double>("beta");
       R0_ = iConfig.getParameter<double>("R0");
     }
 
   }
 
   if ( iConfig.exists("correctShape") ) {
     correctShape_ = iConfig.getParameter<bool>("correctShape");
     gridMaxRapidity_ = iConfig.getParameter<double>("gridMaxRapidity");
     gridSpacing_ = iConfig.getParameter<double>("gridSpacing");
     useExplicitGhosts_ = true;
   }
 
   input_chrefcand_token_ = consumes<edm::View<reco::RecoChargedRefCandidate> >(src_);
 
 }

FastjetJetProducer::~FastjetJetProducer ( )

virtual

Definition at line 235 of file FastjetJetProducer.cc.

236 {

237 }

Member Function Documentation

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

virtual

Reimplemented from VirtualJetProducer.

Reimplemented in cms::CATopJetProducer.

Definition at line 244 of file FastjetJetProducer.cc.

References VirtualJetProducer::fjClusterSeq_, VirtualJetProducer::jetTypeE, VirtualJetProducer::makeTrackJet(), VirtualJetProducer::produce(), and produceTrackJets().

Referenced by cms::CATopJetProducer::produce().

 {
 
   // for everything but track jets
   if (!makeTrackJet(jetTypeE)) {
  
      // use the default production from one collection
      VirtualJetProducer::produce( iEvent, iSetup );
 
   } else { // produce trackjets from tracks grouped per primary vertex
 
     produceTrackJets(iEvent, iSetup);
   
   }
 
   // fjClusterSeq_ retains quite a lot of memory - about 1 to 7Mb at 200 pileup
   // depending on the exact configuration; and there are 24 FastjetJetProducers in the
   // sequence so this adds up to about 60 Mb. It's allocated every time runAlgorithm
   // is called, so safe to delete here.
   fjClusterSeq_.reset();
 }

void FastjetJetProducer::produceTrackJets	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

protectedvirtual

Definition at line 267 of file FastjetJetProducer.cc.

Referenced by produce().

 {
 
     // read in the track candidates
   edm::Handle<edm::View<reco::RecoChargedRefCandidate> > inputsHandle;
     iEvent.getByToken(input_chrefcand_token_, inputsHandle);
 
     // make collection with pointers so we can play around with it
     std::vector<edm::Ptr<reco::RecoChargedRefCandidate> > allInputs;
     std::vector<edm::Ptr<reco::Candidate> > origInputs;
     for (size_t i = 0; i < inputsHandle->size(); ++i) {
       allInputs.push_back(inputsHandle->ptrAt(i));
       origInputs.push_back(inputsHandle->ptrAt(i));
     }
 
     // read in the PV collection
     edm::Handle<reco::VertexCollection> pvCollection;
     iEvent.getByToken(input_vertex_token_, pvCollection);
     // define the overall output jet container
     std::auto_ptr<std::vector<reco::TrackJet> > jets(new std::vector<reco::TrackJet>() );
 
     // loop over the good vertices, clustering for each vertex separately
     for (reco::VertexCollection::const_iterator itVtx = pvCollection->begin(); itVtx != pvCollection->end(); ++itVtx) {
       if (itVtx->isFake() || itVtx->ndof() < minVtxNdof_ || fabs(itVtx->z()) > maxVtxZ_) continue;
 
       // clear the intermediate containers
       inputs_.clear();
       fjInputs_.clear();
       fjJets_.clear();
 
       // if only vertex-associated tracks should be used
       if (useOnlyVertexTracks_) {
         // loop over the tracks associated to the vertex
         for (reco::Vertex::trackRef_iterator itTr = itVtx->tracks_begin(); itTr != itVtx->tracks_end(); ++itTr) {
           // whether a match was found in the track candidate input
           bool found = false;
           // loop over input track candidates
           for (std::vector<edm::Ptr<reco::RecoChargedRefCandidate> >::iterator itIn = allInputs.begin(); itIn != allInputs.end(); ++itIn) {
             // match the input track candidate to the track from the vertex
             reco::TrackRef trref(itTr->castTo<reco::TrackRef>());
             // check if the tracks match
             if ((*itIn)->track() == trref) {
               found = true;
               // add this track candidate to the input for clustering
               inputs_.push_back(*itIn);
               // erase the track candidate from the total list of input, so we don't reuse it later
               allInputs.erase(itIn);
               // found the candidate track corresponding to the vertex track, so stop the loop
               break;
             } // end if match found
           } // end loop over input tracks
           // give an info message in case no match is found (can happen if candidates are subset of tracks used for clustering)
           if (!found) edm::LogInfo("FastjetTrackJetProducer") << "Ignoring a track at vertex which is not in input track collection!";
         } // end loop over tracks associated to vertex
       // if all inpt track candidates should be used
       } else {
         // loop over input track candidates
         for (std::vector<edm::Ptr<reco::RecoChargedRefCandidate> >::iterator itIn = allInputs.begin(); itIn != allInputs.end(); ++itIn) {
           // check if the track is close enough to the vertex
           float dz = (*itIn)->track()->dz(itVtx->position());
           float dxy = (*itIn)->track()->dxy(itVtx->position());
           if (fabs(dz) > dzTrVtxMax_) continue;
           if (fabs(dxy) > dxyTrVtxMax_) continue;
           bool closervtx = false;
           // now loop over the good vertices a second time
           for (reco::VertexCollection::const_iterator itVtx2 = pvCollection->begin(); itVtx2 != pvCollection->end(); ++itVtx2) {
             if (itVtx->isFake() || itVtx->ndof() < minVtxNdof_ || fabs(itVtx->z()) > maxVtxZ_) continue;
             // and check this track is closer to any other vertex (if more than 1 vertex considered)
             if (!useOnlyOnePV_ &&
                 itVtx != itVtx2 &&
                 fabs((*itIn)->track()->dz(itVtx2->position())) < fabs(dz)) {
               closervtx = true;
               break; // 1 closer vertex makes the track already not matched, so break
             }
           }
           // don't add this track if another vertex is found closer
           if (closervtx) continue;
           // add this track candidate to the input for clustering
           inputs_.push_back(*itIn);
           // erase the track candidate from the total list of input, so we don't reuse it later
           allInputs.erase(itIn);
           // take a step back in the loop since we just erased
           --itIn;
         }
       }
 
       // convert candidates in inputs_ to fastjet::PseudoJets in fjInputs_
       fjInputs_.reserve(inputs_.size());
       inputTowers();
       LogDebug("FastjetTrackJetProducer") << "Inputted towers\n";
 
       // run algorithm, using fjInputs_, modifying fjJets_ and allocating fjClusterSeq_
       runAlgorithm(iEvent, iSetup);
       LogDebug("FastjetTrackJetProducer") << "Ran algorithm\n";
 
       // convert our jets and add to the overall jet vector
       for (unsigned int ijet=0;ijet<fjJets_.size();++ijet) {
         // get the constituents from fastjet
         std::vector<fastjet::PseudoJet> fjConstituents = sorted_by_pt(fjClusterSeq_->constituents(fjJets_[ijet]));
         // convert them to CandidatePtr vector
         std::vector<reco::CandidatePtr> constituents = getConstituents(fjConstituents);
         // fill the trackjet
         reco::TrackJet jet;
         // write the specifics to the jet (simultaneously sets 4-vector, vertex).
         writeSpecific( jet,
                        reco::Particle::LorentzVector(fjJets_[ijet].px(), fjJets_[ijet].py(), fjJets_[ijet].pz(), fjJets_[ijet].E()),
                        vertex_, constituents, iSetup);
         jet.setJetArea(0);
         jet.setPileup(0);
         jet.setPrimaryVertex(edm::Ref<reco::VertexCollection>(pvCollection, (int) (itVtx-pvCollection->begin())));
         jet.setVertex(itVtx->position());
         jets->push_back(jet);
       }
 
       if (useOnlyOnePV_) break; // stop vertex loop if only one vertex asked for
     } // end loop over vertices
 
     // put the jets in the collection
     LogDebug("FastjetTrackJetProducer") << "Put " << jets->size() << " jets in the event.\n";
     iEvent.put(jets);
 
     // Clear the work vectors so that memory is free for other modules.
     // Use the trick of swapping with an empty vector so that the memory
     // is actually given back rather than silently kept.
     decltype(fjInputs_)().swap(fjInputs_);
     decltype(fjJets_)().swap(fjJets_);
     decltype(inputs_)().swap(inputs_);  
 }

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

protectedvirtual

Implements VirtualJetProducer.

Reimplemented in cms::CATopJetProducer.

Definition at line 398 of file FastjetJetProducer.cc.

Referenced by produceTrackJets().

 {
   // run algorithm
   /*
   fjInputs_.clear();
   double px, py , pz, E;
   string line;
   std::ifstream fin("dump3.txt");
   while (getline(fin, line)){
     if (line == "#END") break;
     if (line.substr(0,1) == "#") {continue;}
     istringstream istr(line);
     istr >> px >> py >> pz >> E;
     // create a fastjet::PseudoJet with these components and put it onto
     // back of the input_particles vector
     fastjet::PseudoJet j(px,py,pz,E);
     //if ( fabs(j.rap()) < inputEtaMax )
     fjInputs_.push_back(fastjet::PseudoJet(px,py,pz,E)); 
   }
   fin.close();
   */
 
   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_) ) );
   }
 
   if ( !(useMassDropTagger_ || useCMSBoostedTauSeedingAlgorithm_ || useTrimming_ || useFiltering_ || usePruning_ || useSoftDrop_ || useConstituentSubtraction_ ) ) {
     fjJets_ = fastjet::sorted_by_pt(fjClusterSeq_->inclusive_jets(jetPtMin_));
   } else {
     fjJets_.clear();
 
 
     transformer_coll transformers;
 
 
     std::vector<fastjet::PseudoJet> tempJets = fastjet::sorted_by_pt(fjClusterSeq_->inclusive_jets(jetPtMin_));
 
     unique_ptr<fastjet::JetMedianBackgroundEstimator> bge_rho;
     if ( useConstituentSubtraction_ ) {
       fastjet::Selector rho_range =  fastjet::SelectorAbsRapMax(csRho_EtaMax_);
       bge_rho = unique_ptr<fastjet::JetMedianBackgroundEstimator> (new  fastjet::JetMedianBackgroundEstimator(rho_range, fastjet::JetDefinition(fastjet::kt_algorithm, csRParam_), *fjAreaDefinition_) );
       bge_rho->set_particles(fjInputs_);
       fastjet::contrib::ConstituentSubtractor * constituentSubtractor = new fastjet::contrib::ConstituentSubtractor(bge_rho.get());
       // this sets the same background estimator to be used for deltaMass density, rho_m, as for pt density, rho:
       constituentSubtractor->use_common_bge_for_rho_and_rhom(true); // for massless input particles it does not make any difference (rho_m is always zero)
 
       transformers.push_back( transformer_ptr(constituentSubtractor) );
     };
     if ( useMassDropTagger_ ) {
       fastjet::MassDropTagger * md_tagger = new fastjet::MassDropTagger ( muCut_, yCut_ );
       transformers.push_back( transformer_ptr(md_tagger) );
     }
     if ( useCMSBoostedTauSeedingAlgorithm_ ) {
       fastjet::contrib::CMSBoostedTauSeedingAlgorithm * tau_tagger = 
     new fastjet::contrib::CMSBoostedTauSeedingAlgorithm ( subjetPtMin_, muMin_, muMax_, yMin_, yMax_, dRMin_, dRMax_, maxDepth_, verbosity_ );
       transformers.push_back( transformer_ptr(tau_tagger ));
     }
     if ( useTrimming_ ) {
       fastjet::Filter * trimmer = new fastjet::Filter(fastjet::JetDefinition(fastjet::kt_algorithm, rFilt_), fastjet::SelectorPtFractionMin(trimPtFracMin_));
       transformers.push_back( transformer_ptr(trimmer) );
     } 
     if ( (useFiltering_) && (!useDynamicFiltering_) ) {
       fastjet::Filter * filter = new fastjet::Filter(fastjet::JetDefinition(fastjet::cambridge_algorithm, rFilt_), fastjet::SelectorNHardest(nFilt_));
       transformers.push_back( transformer_ptr(filter));
     } 
 
     if ( (usePruning_)  && (!useKtPruning_) ) {
       fastjet::Pruner * pruner = new fastjet::Pruner(fastjet::cambridge_algorithm, zCut_, RcutFactor_);
       transformers.push_back( transformer_ptr(pruner ));
     }
 
     if ( useDynamicFiltering_ ){
       fastjet::Filter * filter = new fastjet::Filter( fastjet::Filter(&*rFiltDynamic_, fastjet::SelectorNHardest(nFilt_)));
       transformers.push_back( transformer_ptr(filter));
     }
 
     if ( useKtPruning_ ) {
       fastjet::Pruner * pruner = new fastjet::Pruner(fastjet::kt_algorithm, zCut_, RcutFactor_);
       transformers.push_back( transformer_ptr(pruner ));
     }
 
     if ( useSoftDrop_ ) {
       fastjet::contrib::SoftDrop * sd = new fastjet::contrib::SoftDrop(beta_, zCut_, R0_ );
       transformers.push_back( transformer_ptr(sd) );
     }
 
     unique_ptr<fastjet::Subtractor> subtractor;
     unique_ptr<fastjet::GridMedianBackgroundEstimator> bge_rho_grid;
     if ( correctShape_ ) {
       bge_rho_grid = unique_ptr<fastjet::GridMedianBackgroundEstimator> (new  fastjet::GridMedianBackgroundEstimator(gridMaxRapidity_, gridSpacing_) );
       bge_rho_grid->set_particles(fjInputs_);
       subtractor = unique_ptr<fastjet::Subtractor>( new fastjet::Subtractor(  bge_rho_grid.get()) );
       subtractor->set_use_rho_m();
       //subtractor->use_common_bge_for_rho_and_rhom(true);
     }
 
 
     for ( std::vector<fastjet::PseudoJet>::const_iterator ijet = tempJets.begin(),
         ijetEnd = tempJets.end(); ijet != ijetEnd; ++ijet ) {
 
       fastjet::PseudoJet transformedJet = *ijet;
       bool passed = true;
       for ( transformer_coll::const_iterator itransf = transformers.begin(),
           itransfEnd = transformers.end(); itransf != itransfEnd; ++itransf ) {
     if ( transformedJet != 0 ) {
       transformedJet = (**itransf)(transformedJet);
     } else {
       passed=false;
     }
       }
 
       if ( correctShape_ ) {
     transformedJet = (*subtractor)(transformedJet);
       }
 
       if ( passed ) {
     fjJets_.push_back( transformedJet );
       }
     }
   }
 
 }

Member Data Documentation

double FastjetJetProducer::beta_

private

for constituent subtraction : R parameter for KT alg in jet median background estimator

Definition at line 99 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::correctShape_

private

Soft drop.

Definition at line 87 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::csRho_EtaMax_

private

for pruning: constituent dR * pt/2m < rcut_factor

Definition at line 97 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::csRParam_

private

for constituent subtraction : maximum rapidity for ghosts

Definition at line 98 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::dRMax_

private

for CMSBoostedTauSeedingAlgorithm : min dR

Definition at line 111 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::dRMin_

private

for CMSBoostedTauSeedingAlgorithm : max asymmetry

Definition at line 110 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

float FastjetJetProducer::dxyTrVtxMax_

private

Definition at line 73 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and produceTrackJets().

float FastjetJetProducer::dzTrVtxMax_

private

Definition at line 72 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and produceTrackJets().

double FastjetJetProducer::gridMaxRapidity_

private

for soft drop : R0 (angular distance normalization - should be set to jet radius in most cases)

Definition at line 101 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::gridSpacing_

private

for shape subtraction, get the fixed-grid rho

Definition at line 102 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

edm::EDGetTokenT<edm::View<reco::RecoChargedRefCandidate> > FastjetJetProducer::input_chrefcand_token_

private

for CMSBoostedTauSeedingAlgorithm : max depth for descending into clustering sequence

Definition at line 116 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and produceTrackJets().

int FastjetJetProducer::maxDepth_

private

for CMSBoostedTauSeedingAlgorithm : max dR

Definition at line 112 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

float FastjetJetProducer::maxVtxZ_

private

Definition at line 75 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and produceTrackJets().

int FastjetJetProducer::minVtxNdof_

private

Definition at line 74 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and produceTrackJets().

double FastjetJetProducer::muCut_

private

Correct the shape of the jets.

Definition at line 88 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::muMax_

private

for CMSBoostedTauSeedingAlgorithm : min mass-drop

Definition at line 107 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::muMin_

private

for CMSBoostedTauSeedingAlgorithm : subjet pt min

Definition at line 106 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

int FastjetJetProducer::nFilt_

private

for dynamic filtering radius (as in arXiv:0802.2470)

Definition at line 93 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::R0_

private

for soft drop : beta (angular exponent)

Definition at line 100 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::RcutFactor_

private

for pruning OR soft drop: constituent minimum pt fraction of parent cluster

Definition at line 96 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::rFilt_

private

for mass-drop tagging, symmetry cut: min(pt1^2,pt2^2) * dR(1,2) / mjet > ycut

Definition at line 90 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

DynamicRfiltPtr FastjetJetProducer::rFiltDynamic_

private

for filtering, trimming: dR scale of sub-clustering

Definition at line 91 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::rFiltFactor_

private

for dynamic filtering radius (as in arXiv:0802.2470)

Definition at line 92 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer().

double FastjetJetProducer::subjetPtMin_

private

for shape subtraction, get the grid spacing

Definition at line 105 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::trimPtFracMin_

private

for filtering, pruning: number of subjets expected

Definition at line 94 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::useCMSBoostedTauSeedingAlgorithm_

private

Jet pruning technique.

Definition at line 83 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::useConstituentSubtraction_

private

Use Kt clustering algorithm for pruning (default is Cambridge/Aachen)

Definition at line 85 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::useDynamicFiltering_

private

Jet filtering technique.

Definition at line 80 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::useFiltering_

private

Mass-drop tagging for boosted Higgs.

Definition at line 79 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::useKtPruning_

private

algorithm for seeding reconstruction of boosted Taus (similar to mass-drop tagging)

Definition at line 84 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::useMassDropTagger_

private

Definition at line 78 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::useOnlyOnePV_

private

Definition at line 71 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and produceTrackJets().

bool FastjetJetProducer::useOnlyVertexTracks_

private

Definition at line 70 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and produceTrackJets().

bool FastjetJetProducer::usePruning_

private

Jet trimming technique.

Definition at line 82 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::useSoftDrop_

private

constituent subtraction technique

Definition at line 86 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

bool FastjetJetProducer::useTrimming_

private

Use dynamic filtering radius (as in arXiv:0802.2470)

Definition at line 81 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::yCut_

private

for mass-drop tagging, m0/mjet (m0 = mass of highest mass subjet)

Definition at line 89 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::yMax_

private

for CMSBoostedTauSeedingAlgorithm : min asymmetry

Definition at line 109 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::yMin_

private

for CMSBoostedTauSeedingAlgorithm : max mass-drop

Definition at line 108 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

double FastjetJetProducer::zCut_

private

for trimming: constituent minimum pt fraction of full jet

Definition at line 95 of file FastjetJetProducer.h.

Referenced by FastjetJetProducer(), and runAlgorithm().

Public Types

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation