#include <SubjetFilterJetProducer.h>

Inheritance diagram for SubjetFilterJetProducer:

Public Member Functions
void	endJob ()

void	inputTowers ()

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

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

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

	SubjetFilterJetProducer (const edm::ParameterSet &ps)

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

virtual	~SubjetFilterJetProducer ()

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

	VirtualJetProducer (const edm::ParameterSet &iConfig)

virtual	~VirtualJetProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Attributes
SubjetFilterAlgorithm	alg_

std::vector< CompoundPseudoJet >	fjCompoundJets_

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

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

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

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

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_

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_

reco::Particle::Point	vertex_

double	voronoiRfact_

bool	writeCompound_

Detailed Description

Definition at line 25 of file SubjetFilterJetProducer.h.

Constructor & Destructor Documentation

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

Definition at line 41 of file SubjetFilterJetProducer.cc.

References VirtualJetProducer::makeProduces(), and VirtualJetProducer::moduleLabel_.

   : VirtualJetProducer(iConfig)
   , alg_(iConfig.getParameter<string>  ("@module_label"),
      iConfig.getParameter<string>  ("jetAlgorithm"),
      iConfig.getParameter<unsigned>("nFatMax"),
      iConfig.getParameter<double>  ("rParam"),
      iConfig.getParameter<double>  ("rFilt"),
      iConfig.getParameter<double>  ("jetPtMin"),
      iConfig.getParameter<double>  ("massDropCut"),
      iConfig.getParameter<double>  ("asymmCut"),
      iConfig.getParameter<bool>    ("asymmCutLater"),
      iConfig.getParameter<bool>    ("doAreaFastjet"),
      iConfig.getParameter<double>  ("Ghost_EtaMax"),
      iConfig.getParameter<int>     ("Active_Area_Repeats"),
      iConfig.getParameter<double>  ("GhostArea"),
      iConfig.getUntrackedParameter<bool>("verbose",false))
 {
   produces<reco::BasicJetCollection>("fat");
   makeProduces(moduleLabel_,"sub");
   makeProduces(moduleLabel_,"filter");
 }

SubjetFilterJetProducer::~SubjetFilterJetProducer ( )

virtual

Definition at line 65 of file SubjetFilterJetProducer.cc.

 {
   
 }

Member Function Documentation

void SubjetFilterJetProducer::endJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 84 of file SubjetFilterJetProducer.cc.

References alg_, gather_cfg::cout, and SubjetFilterAlgorithm::summary().

 {
   cout<<alg_.summary()<<endl;
 }

void SubjetFilterJetProducer::inputTowers ( )

virtual

Reimplemented from VirtualJetProducer.

Definition at line 99 of file SubjetFilterJetProducer.cc.

References fjCompoundJets_, and VirtualJetProducer::inputTowers().

 {
   fjCompoundJets_.clear();
   VirtualJetProducer::inputTowers();
 }

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

virtual

Reimplemented from VirtualJetProducer.

Definition at line 107 of file SubjetFilterJetProducer.cc.

References VirtualJetProducer::JetType::BasicJet, VirtualJetProducer::JetType::CaloJet, VirtualJetProducer::JetType::GenJet, iEvent, VirtualJetProducer::jetTypeE, and VirtualJetProducer::JetType::PFJet.

 {
   // Write jets and constitutents. Will use fjCompoundJets_. 
   switch( jetTypeE ) {
   case JetType::CaloJet :
     writeCompoundJets<reco::CaloJet>( iEvent, iSetup );
     break;
   case JetType::PFJet :
     writeCompoundJets<reco::PFJet>( iEvent, iSetup );
     break;
   case JetType::GenJet :
     writeCompoundJets<reco::GenJet>( iEvent, iSetup );
     break;
   case JetType::BasicJet :
     writeCompoundJets<reco::BasicJet>( iEvent, iSetup );
     break;
   default:
     edm::LogError("InvalidInput")<<" invalid jet type in SubjetFilterJetProducer\n";
     break;
   };
   
 }

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

virtual

Reimplemented from VirtualJetProducer.

Definition at line 76 of file SubjetFilterJetProducer.cc.

References VirtualJetProducer::produce().

 {
   VirtualJetProducer::produce(iEvent,iSetup);
 }

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

virtual

Implements VirtualJetProducer.

Definition at line 91 of file SubjetFilterJetProducer.cc.

References alg_, fjCompoundJets_, VirtualJetProducer::fjInputs_, and SubjetFilterAlgorithm::run().

 {
   alg_.run(fjInputs_, fjCompoundJets_, iSetup);
 }

template<class T >

void SubjetFilterJetProducer::writeCompoundJets	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

Definition at line 134 of file SubjetFilterJetProducer.cc.

References fjCompoundJets_, VirtualJetProducer::inputs_, point, edm::Event::put(), and reco::writeSpecific().

 {
   auto_ptr<reco::BasicJetCollection> fatJets( new reco::BasicJetCollection() );
   auto_ptr<vector<T> >  subJets( new vector<T>() );
   auto_ptr<vector<T> >  filterJets( new vector<T>() );
 
   edm::OrphanHandle< vector<T> > subJetsAfterPut;
   edm::OrphanHandle< vector<T> > filterJetsAfterPut;
   
   vector< vector<int> > subIndices(fjCompoundJets_.size());
   vector< vector<int> > filterIndices(fjCompoundJets_.size());
   
   vector<math::XYZTLorentzVector> p4FatJets;
   vector<double>                  areaFatJets;
   
   vector<CompoundPseudoJet>::const_iterator itBegin(fjCompoundJets_.begin());
   vector<CompoundPseudoJet>::const_iterator itEnd(fjCompoundJets_.end());
   vector<CompoundPseudoJet>::const_iterator it(itBegin);
   
   for (;it!=itEnd;++it) {
 
     int jetIndex = it-itBegin;
     fastjet::PseudoJet fatJet = it->hardJet();
     p4FatJets.push_back(math::XYZTLorentzVector(fatJet.px(),fatJet.py(),
                         fatJet.pz(),fatJet.e()));
     areaFatJets.push_back(it->hardJetArea());
     
     vector<CompoundPseudoSubJet>::const_iterator itSubBegin(it->subjets().begin());
     vector<CompoundPseudoSubJet>::const_iterator itSubEnd(it->subjets().end());
     vector<CompoundPseudoSubJet>::const_iterator itSub(itSubBegin);
     
     for (; itSub!=itSubEnd;++itSub) {
       int subJetIndex = itSub-itSubBegin;
       fastjet::PseudoJet fjSubJet = itSub->subjet();
       math::XYZTLorentzVector p4SubJet(fjSubJet.px(),fjSubJet.py(),
                        fjSubJet.pz(),fjSubJet.e());
       reco::Particle::Point point(0,0,0);
       
       vector<reco::CandidatePtr> subJetConstituents;
       const vector<int>& subJetConstituentIndices = itSub->constituents();
       vector<int>::const_iterator itIndexBegin(subJetConstituentIndices.begin());
       vector<int>::const_iterator itIndexEnd(subJetConstituentIndices.end());
       vector<int>::const_iterator itIndex(itIndexBegin);
       for (;itIndex!=itIndexEnd;++itIndex)
     if ((*itIndex) < static_cast<int>(inputs_.size())) 
       subJetConstituents.push_back(inputs_[*itIndex]);
       
       T subJet;
       reco::writeSpecific(subJet,p4SubJet,point,subJetConstituents,iSetup);
       subJet.setJetArea(itSub->subjetArea());
       
       if (subJetIndex<2) {
     subIndices[jetIndex].push_back(subJets->size());
     subJets->push_back(subJet);
       }
       else {
     filterIndices[jetIndex].push_back(filterJets->size());
     filterJets->push_back(subJet);
       }
     }
   }
   
   subJetsAfterPut    = iEvent.put(subJets,   "sub");
   filterJetsAfterPut = iEvent.put(filterJets,"filter");
   
   vector<math::XYZTLorentzVector>::const_iterator itP4Begin(p4FatJets.begin());
   vector<math::XYZTLorentzVector>::const_iterator itP4End(p4FatJets.end());
   vector<math::XYZTLorentzVector>::const_iterator itP4(itP4Begin);
   for (;itP4!=itP4End;++itP4) {
     int fatIndex = itP4-itP4Begin;
     vector<int>& fatToSub    = subIndices[fatIndex];
     vector<int>& fatToFilter = filterIndices[fatIndex];
 
     vector<reco::CandidatePtr> i_fatJetConstituents;
 
     vector<int>::const_iterator itSubBegin(fatToSub.begin());
     vector<int>::const_iterator itSubEnd(fatToSub.end());
     vector<int>::const_iterator itSub(itSubBegin);
     for(;itSub!=itSubEnd;++itSub) {
       reco::CandidatePtr candPtr(subJetsAfterPut,(*itSub),false);
       i_fatJetConstituents.push_back(candPtr);
     }
 
     vector<int>::const_iterator itFilterBegin(fatToFilter.begin());
     vector<int>::const_iterator itFilterEnd(fatToFilter.end());
     vector<int>::const_iterator itFilter(itFilterBegin);
     for(;itFilter!=itFilterEnd;++itFilter) {
       reco::CandidatePtr candPtr(filterJetsAfterPut,(*itFilter),false);
       i_fatJetConstituents.push_back(candPtr);
     }
 
     reco::Particle::Point point(0,0,0);
     fatJets->push_back(reco::BasicJet((*itP4),point,i_fatJetConstituents));
     fatJets->back().setJetArea(areaFatJets[fatIndex]);
   }
   
   iEvent.put(fatJets,"fat");
 }

Member Data Documentation

SubjetFilterAlgorithm SubjetFilterJetProducer::alg_

private

Definition at line 52 of file SubjetFilterJetProducer.h.

Referenced by endJob(), and runAlgorithm().

std::vector<CompoundPseudoJet> SubjetFilterJetProducer::fjCompoundJets_

private

Definition at line 53 of file SubjetFilterJetProducer.h.

Referenced by inputTowers(), runAlgorithm(), and writeCompoundJets().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation