VirtualJetProducer Class Reference

#include <VirtualJetProducer.h>

Inheritance diagram for VirtualJetProducer:

Classes
struct	JetType

Public Types
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
std::string	jetType () const
virtual void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	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)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
virtual	~ProducerBase () noexcept(false)
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
static void	fillDescriptionsFromVirtualJetProducer (edm::ParameterSetDescription &desc)
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 Member Functions
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)
virtual void	runAlgorithm (edm::Event &iEvent, const edm::EventSetup &iSetup)=0
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 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)

Protected Attributes
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_
RangeDefPtr	fjRangeDef_
bool	fromHTTTopJetProducer_
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_

Private Attributes
std::auto_ptr< AnomalousTower >	anomalousTowerDef_
edm::EDGetTokenT< std::vector< edm::FwdPtr< reco::PFCandidate > > >	input_candidatefwdptr_token_
edm::EDGetTokenT< reco::CandidateView >	input_candidateview_token_
edm::EDGetTokenT< std::vector< edm::FwdPtr< pat::PackedCandidate > > >	input_packedcandidatefwdptr_token_

Detailed Description

Definition at line 35 of file VirtualJetProducer.h.

Member Typedef Documentation

typedef boost::shared_ptr<fastjet::GhostedAreaSpec> VirtualJetProducer::ActiveAreaSpecPtr

Definition at line 90 of file VirtualJetProducer.h.

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

Definition at line 91 of file VirtualJetProducer.h.

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

Definition at line 87 of file VirtualJetProducer.h.

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

Definition at line 89 of file VirtualJetProducer.h.

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

Definition at line 88 of file VirtualJetProducer.h.

typedef boost::shared_ptr<fastjet::RangeDefinition> VirtualJetProducer::RangeDefPtr

Definition at line 92 of file VirtualJetProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 127 of file VirtualJetProducer.cc.

References diMuonSkim_cfi::byName, Exception, reco::get(), and edm::ParameterSet::getParameter().

                                                                     {

    moduleLabel_        = iConfig.getParameter<string>  ("@module_label");
        src_                    = iConfig.getParameter<edm::InputTag>("src");
        srcPVs_                 = iConfig.getParameter<edm::InputTag>("srcPVs");
    jetType_            = iConfig.getParameter<string>  ("jetType");
    jetAlgorithm_       = iConfig.getParameter<string>  ("jetAlgorithm");
    rParam_             = iConfig.getParameter<double>  ("rParam");
    inputEtMin_         = iConfig.getParameter<double>  ("inputEtMin");
    inputEMin_          = iConfig.getParameter<double>  ("inputEMin");
    jetPtMin_           = iConfig.getParameter<double>  ("jetPtMin");
    doPVCorrection_     = iConfig.getParameter<bool>    ("doPVCorrection");
    doAreaFastjet_      = iConfig.getParameter<bool>    ("doAreaFastjet");
    doRhoFastjet_       = iConfig.getParameter<bool>    ("doRhoFastjet");
    jetCollInstanceName_    = iConfig.getParameter<string>  ("jetCollInstanceName");
    doPUOffsetCorr_     = iConfig.getParameter<bool>    ("doPUOffsetCorr");
    puSubtractorName_   = iConfig.getParameter<string>  ("subtractorName");
    useExplicitGhosts_  = iConfig.getParameter<bool>    ("useExplicitGhosts");  // use explicit ghosts in the fastjet clustering sequence?
    doAreaDiskApprox_   = iConfig.getParameter<bool>    ("doAreaDiskApprox");
    voronoiRfact_       = iConfig.getParameter<double>  ("voronoiRfact");   // Voronoi-based area calculation allows for an empirical scale factor
    rhoEtaMax_      = iConfig.getParameter<double>  ("Rho_EtaMax");         // do fasjet area / rho calcluation? => accept corresponding parameters
    ghostEtaMax_        = iConfig.getParameter<double>  ("Ghost_EtaMax");
    activeAreaRepeats_  = iConfig.getParameter<int>     ("Active_Area_Repeats");
    ghostArea_      = iConfig.getParameter<double>  ("GhostArea");
    restrictInputs_     = iConfig.getParameter<bool>    ("restrictInputs");     // restrict inputs to first "maxInputs" towers?
    maxInputs_          = iConfig.getParameter<unsigned int>("maxInputs");
    writeCompound_      = iConfig.getParameter<bool>    ("writeCompound");  // Check to see if we are writing compound jets for substructure and jet grooming
        writeJetsWithConst_     = iConfig.getParameter<bool>("writeJetsWithConst"); //write subtracted jet constituents
    doFastJetNonUniform_    = iConfig.getParameter<bool>    ("doFastJetNonUniform");
    puCenters_      = iConfig.getParameter<vector<double> >("puCenters");
    puWidth_        = iConfig.getParameter<double>  ("puWidth");
    nExclude_       = iConfig.getParameter<unsigned int>("nExclude");
    useDeterministicSeed_   = iConfig.getParameter<bool>    ("useDeterministicSeed");
    minSeed_        = iConfig.getParameter<unsigned int>("minSeed");
    verbosity_      = iConfig.getParameter<int> ("verbosity");

    anomalousTowerDef_ = auto_ptr<AnomalousTower>(new AnomalousTower(iConfig));

    input_vertex_token_ = consumes<reco::VertexCollection>(srcPVs_);
    input_candidateview_token_ = consumes<reco::CandidateView>(src_);
    input_candidatefwdptr_token_ = consumes<vector<edm::FwdPtr<reco::PFCandidate> > >(iConfig.getParameter<edm::InputTag>("src"));
    input_packedcandidatefwdptr_token_ = consumes<vector<edm::FwdPtr<pat::PackedCandidate> > >(iConfig.getParameter<edm::InputTag>("src"));
    //
    // additional parameters to think about:
    // - overlap threshold (set to 0.75 for the time being)
    // - p parameter for generalized kT (set to -2 for the time being)
    // - fastjet PU subtraction parameters (not yet considered)
    //
    if (jetAlgorithm_=="Kt") 
        fjJetDefinition_= JetDefPtr(new fastjet::JetDefinition(fastjet::kt_algorithm,rParam_));

    else if (jetAlgorithm_=="CambridgeAachen")
        fjJetDefinition_= JetDefPtr(new fastjet::JetDefinition(fastjet::cambridge_algorithm,rParam_) );

    else if (jetAlgorithm_=="AntiKt")
        fjJetDefinition_= JetDefPtr( new fastjet::JetDefinition(fastjet::antikt_algorithm,rParam_) );

    else if (jetAlgorithm_=="GeneralizedKt") 
        fjJetDefinition_= JetDefPtr( new fastjet::JetDefinition(fastjet::genkt_algorithm,rParam_,-2) );

    else if (jetAlgorithm_=="SISCone") {

        fjPlugin_ = PluginPtr( new fastjet::SISConePlugin(rParam_,0.75,0,0.0,false,fastjet::SISConePlugin::SM_pttilde) );
        fjJetDefinition_= JetDefPtr( new fastjet::JetDefinition(&*fjPlugin_) );

    } else if (jetAlgorithm_=="IterativeCone") {

        fjPlugin_ = PluginPtr(new fastjet::CMSIterativeConePlugin(rParam_,1.0));
        fjJetDefinition_= JetDefPtr(new fastjet::JetDefinition(&*fjPlugin_));

    } else if (jetAlgorithm_=="CDFMidPoint") {

        fjPlugin_ = PluginPtr(new fastjet::CDFMidPointPlugin(rParam_,0.75));
        fjJetDefinition_= JetDefPtr(new fastjet::JetDefinition(&*fjPlugin_));

    } else if (jetAlgorithm_=="ATLASCone") {

        fjPlugin_ = PluginPtr(new fastjet::ATLASConePlugin(rParam_));
        fjJetDefinition_= JetDefPtr(new fastjet::JetDefinition(&*fjPlugin_));

    } else {
        throw cms::Exception("Invalid jetAlgorithm") <<"Jet algorithm for VirtualJetProducer is invalid, Abort!\n";
    }

    jetTypeE=JetType::byName(jetType_);

    if ( doPUOffsetCorr_  ) {
        if(puSubtractorName_.empty()){
            LogWarning("VirtualJetProducer") << "Pile Up correction on; however, pile up type is not specified. Using default... \n";
            subtractor_ =  boost::shared_ptr<PileUpSubtractor>(new PileUpSubtractor(iConfig, consumesCollector()));
        } else subtractor_ =  boost::shared_ptr<PileUpSubtractor>(
                PileUpSubtractorFactory::get()->create( puSubtractorName_, iConfig, consumesCollector()));
    }

    // do approximate disk-based area calculation => warn if conflicting request
    if (doAreaDiskApprox_ && doAreaFastjet_)
        throw cms::Exception("Conflicting area calculations") << "Both the calculation of jet area via fastjet and via an analytical disk approximation have been requested. Please decide on one.\n";

    if ( doAreaFastjet_ || doRhoFastjet_ ) {

        if (voronoiRfact_ <= 0) {
            fjActiveArea_     = ActiveAreaSpecPtr(new fastjet::GhostedAreaSpec(ghostEtaMax_,activeAreaRepeats_,ghostArea_));
            fjActiveArea_->set_fj2_placement(true);

            if ( !useExplicitGhosts_ ) {
                fjAreaDefinition_ = AreaDefinitionPtr( new fastjet::AreaDefinition(fastjet::active_area, *fjActiveArea_ ) );
            } else {
                fjAreaDefinition_ = AreaDefinitionPtr( new fastjet::AreaDefinition(fastjet::active_area_explicit_ghosts, *fjActiveArea_ ) );
            }
        }
        fjRangeDef_ = RangeDefPtr( new fastjet::RangeDefinition(rhoEtaMax_) );
    } 

    if( ( doFastJetNonUniform_ ) && ( puCenters_.size() == 0 ) ) 
        throw cms::Exception("doFastJetNonUniform") << "Parameter puCenters for doFastJetNonUniform is not defined." << std::endl;
  
        // make the "produces" statements
        makeProduces( moduleLabel_, jetCollInstanceName_ );
    produces<vector<double> >("rhos");
    produces<vector<double> >("sigmas");
    produces<double>("rho");
    produces<double>("sigma");

  
}

VirtualJetProducer::~VirtualJetProducer ( )

virtual

Definition at line 254 of file VirtualJetProducer.cc.

{
} 

Member Function Documentation

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

inlineprotectedvirtual

Reimplemented in cms::HTTTopJetProducer.

Definition at line 129 of file VirtualJetProducer.h.

References convertSQLitetoXML_cfg::output.

                                                                       {};

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

protectedvirtual

Definition at line 468 of file VirtualJetProducer.cc.

References reco::CompositePtrCandidate::addDaughter(), mps_fire::i, and diffTreeTool::index.

{
  for (unsigned int i=0;i<fjConstituents.size();++i) { 
    int index = fjConstituents[i].user_index();
    if ( index >= 0 && static_cast<unsigned int>(index) < inputs_.size() )
      jet->addDaughter(inputs_[index]);
  }
}

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

static

Definition at line 960 of file VirtualJetProducer.cc.

References edm::ParameterSetDescription::add(), and edm::ConfigurationDescriptions::addDefault().

                                                                                    {

    edm::ParameterSetDescription desc;
        fillDescriptionsFromVirtualJetProducer(desc);
    desc.add<string>("jetCollInstanceName", ""  );

        // addDefault must be used here instead of add unless
        // all the classes that inherit from this class redefine
        // the fillDescriptions function. Otherwise, the autogenerated
        // cfi filenames are the same and conflict.
    descriptions.addDefault(desc);
}

void VirtualJetProducer::fillDescriptionsFromVirtualJetProducer ( edm::ParameterSetDescription &  desc )

static

Definition at line 973 of file VirtualJetProducer.cc.

References edm::ParameterSetDescription::add(), and funct::false.

Referenced by cms::CSJetProducer::fillDescriptions(), FastjetJetProducer::fillDescriptions(), cms::CATopJetProducer::fillDescriptions(), and cms::HTTTopJetProducer::fillDescriptions().

{
    desc.add<edm::InputTag>("src",      edm::InputTag("particleFlow") );
    desc.add<edm::InputTag>("srcPVs",   edm::InputTag("") );
    desc.add<string>("jetType",     "PFJet" );
    desc.add<string>("jetAlgorithm",    "AntiKt" );
    desc.add<double>("rParam",      0.4 );
    desc.add<double>("inputEtMin",      0.0 );
    desc.add<double>("inputEMin",       0.0 );
    desc.add<double>("jetPtMin",        5. );
    desc.add<bool>  ("doPVCorrection",  false );
    desc.add<bool>  ("doAreaFastjet",   false );
    desc.add<bool>  ("doRhoFastjet",    false );
    desc.add<bool>  ("doPUOffsetCorr",  false   );
    desc.add<double>("puPtMin",             10.);
        desc.add<double>("nSigmaPU",            1.0 );
        desc.add<double>("radiusPU",            0.5 );
    desc.add<string>("subtractorName",  ""  );
    desc.add<bool>  ("useExplicitGhosts",   false   );
    desc.add<bool>  ("doAreaDiskApprox",    false   );
    desc.add<double>("voronoiRfact",    -0.9    );
    desc.add<double>("Rho_EtaMax",      4.4     );
    desc.add<double>("Ghost_EtaMax",    5.  );
    desc.add<int>   ("Active_Area_Repeats", 1   );
    desc.add<double>("GhostArea",       0.01    );
    desc.add<bool>  ("restrictInputs",  false   );
    desc.add<unsigned int>  ("maxInputs",   1   );
    desc.add<bool>  ("writeCompound",   false   );
        desc.add<bool>  ("writeJetsWithConst",  false   );
    desc.add<bool>  ("doFastJetNonUniform", false   );
    desc.add<bool>  ("useDeterministicSeed",false   );
    desc.add<unsigned int>  ("minSeed",     14327   );
    desc.add<int>   ("verbosity",       0   );
    desc.add<double>("puWidth",     0.  );
    desc.add<unsigned int>("nExclude",  0   );
    desc.add<unsigned int>("maxBadEcalCells",   9999999 );
    desc.add<unsigned int>("maxBadHcalCells",   9999999 );
    desc.add<unsigned int>("maxProblematicEcalCells",   9999999 );
    desc.add<unsigned int>("maxProblematicHcalCells",   9999999 );
    desc.add<unsigned int>("maxRecoveredEcalCells", 9999999 );
    desc.add<unsigned int>("maxRecoveredHcalCells", 9999999 );
    vector<double>  puCentersDefault;
    desc.add<vector<double>>("puCenters",   puCentersDefault);
}

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

protectedvirtual

Definition at line 481 of file VirtualJetProducer.cc.

References mps_fire::i, diffTreeTool::index, and mps_fire::result.

Referenced by FastjetJetProducer::produceTrackJets(), and cms::SubEventGenJetProducer::runAlgorithm().

{
  vector<reco::CandidatePtr> result; result.reserve(fjConstituents.size()/2);
  for (unsigned int i=0;i<fjConstituents.size();i++) {
    auto index = fjConstituents[i].user_index();
    if ( index >= 0 && static_cast<unsigned int>(index) < inputs_.size() ) {
      result.emplace_back(inputs_[index]);
    }
  }
  return result;
}

void VirtualJetProducer::inputTowers ( )

protectedvirtual

Reimplemented in cms::CompoundJetProducer, SubjetFilterJetProducer, and cms::SubEventGenJetProducer.

Definition at line 409 of file VirtualJetProducer.cc.

References geometryDiff::epsilon, mps_fire::i, input, edm::isNotFinite(), and CaloTower::p4().

Referenced by SubjetFilterJetProducer::inputTowers(), cms::CompoundJetProducer::inputTowers(), and FastjetJetProducer::produceTrackJets().

{
  auto inBegin = inputs_.begin(),
    inEnd = inputs_.end(), i = inBegin;
  for (; i != inEnd; ++i ) {
    auto const & input = **i;
    // std::cout << "CaloTowerVI jets " << input->pt() << " " << input->et() << ' '<< input->energy() << ' ' << (isAnomalousTower(input) ? " bad" : " ok") << std::endl; 
    if (edm::isNotFinite(input.pt()))           continue;
    if (input.et()    <inputEtMin_)  continue;
    if (input.energy()<inputEMin_)   continue;
    if (isAnomalousTower(*i))      continue;
    // Change by SRR : this is no longer an error nor warning, this can happen with PU mitigation algos.
    // Also switch to something more numerically safe. (VI: 10^-42GeV????)
    if (input.pt() < 100 * std::numeric_limits<double>::epsilon() ) { 
      continue;
    }
    if (makeCaloJet(jetTypeE)&&doPVCorrection_) {
      const CaloTower & tower = dynamic_cast<const CaloTower &>(input);
      auto const &  ct = tower.p4(vertex_);  // very expensive as computed in eta/phi
      fjInputs_.emplace_back(ct.px(),ct.py(),ct.pz(),ct.energy());
      //std::cout << "tower:" << *tower << '\n';
    }
    else {
      /*
      if(makePFJet(jetTypeE)) {
    reco::PFCandidate& pfc = (reco::PFCandidate&)input;
    std::cout << "PF cand:" << pfc << '\n';
      }
      */
      fjInputs_.emplace_back(input.px(),input.py(),input.pz(),
                         input.energy());
    }
    fjInputs_.back().set_user_index(i - inBegin);
  }

  if ( restrictInputs_ && fjInputs_.size() > maxInputs_ ) {
    reco::helper::GreaterByPtPseudoJet   pTComparator;
    std::sort(fjInputs_.begin(), fjInputs_.end(), pTComparator);
    fjInputs_.resize(maxInputs_);
    edm::LogWarning("JetRecoTooManyEntries") << "Too many inputs in the event, limiting to first " << maxInputs_ << ". Output is suspect.";
  }
}

bool VirtualJetProducer::isAnomalousTower ( reco::CandidatePtr  input )

protectedvirtual

Definition at line 453 of file VirtualJetProducer.cc.

Referenced by cms::SubEventGenJetProducer::inputTowers().

{
  if (!makeCaloJet(jetTypeE)) 
      return false;
  else
      return (*anomalousTowerDef_)(*input);
}

std::string VirtualJetProducer::jetType ( ) const

inline

Definition at line 99 of file VirtualJetProducer.h.

References iEvent, input, metsig::jet, alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, and GlobalPosition_Frontier_DevDB_cff::tag.

{ return jetType_; }

bool VirtualJetProducer::makeBasicJet ( const JetType::Type &  fTag )

inlineprotected

Definition at line 72 of file VirtualJetProducer.h.

                                                    {
    return fTag == JetType::BasicJet;
  }

bool VirtualJetProducer::makeCaloJet ( const JetType::Type &  fTag )

inlineprotected

Definition at line 57 of file VirtualJetProducer.h.

References btagElecInJet_cfi::CaloJet.

                                                   {
    return fTag == JetType::CaloJet;
  }

bool VirtualJetProducer::makeGenJet ( const JetType::Type &  fTag )

inlineprotected

Definition at line 63 of file VirtualJetProducer.h.

                                                  {
    return fTag == JetType::GenJet;
  }

bool VirtualJetProducer::makePFClusterJet ( const JetType::Type &  fTag )

inlineprotected

Definition at line 69 of file VirtualJetProducer.h.

                                                        {
    return fTag == JetType::PFClusterJet;
  }

bool VirtualJetProducer::makePFJet ( const JetType::Type &  fTag )

inlineprotected

Definition at line 60 of file VirtualJetProducer.h.

                                                 {
    return fTag == JetType::PFJet;
  }

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

protectedvirtual

Definition at line 89 of file VirtualJetProducer.cc.

References GlobalPosition_Frontier_DevDB_cff::tag.

Referenced by SubjetFilterJetProducer::SubjetFilterJetProducer().

{


  if ( writeCompound_ ) {
    produces<reco::BasicJetCollection>();
  }

  if ( writeJetsWithConst_ ) {
    produces<reco::PFCandidateCollection>(tag).setBranchAlias(alias);
    produces<reco::PFJetCollection>();
  } else {
    if (makeCaloJet(jetTypeE)) {
      produces<reco::CaloJetCollection>(tag).setBranchAlias(alias);
    }
    else if (makePFJet(jetTypeE)) {
      produces<reco::PFJetCollection>(tag).setBranchAlias(alias);
    }
    else if (makeGenJet(jetTypeE)) {
      produces<reco::GenJetCollection>(tag).setBranchAlias(alias);
    }
    else if (makeTrackJet(jetTypeE)) {
      produces<reco::TrackJetCollection>(tag).setBranchAlias(alias);
    }
    else if (makePFClusterJet(jetTypeE)) {
      produces<reco::PFClusterJetCollection>(tag).setBranchAlias(alias);
    }
    else if (makeBasicJet(jetTypeE)) {
      produces<reco::BasicJetCollection>(tag).setBranchAlias(alias);
    }
  }
}

bool VirtualJetProducer::makeTrackJet ( const JetType::Type &  fTag )

inlineprotected

Definition at line 66 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::produce().

                                                    {
    return fTag == JetType::TrackJet;
  }

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

protected

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

protectedvirtual

Reimplemented in cms::CompoundJetProducer.

Definition at line 496 of file VirtualJetProducer.cc.

References btagElecInJet_cfi::CaloJet, Exception, iEvent, and metsig::jet.

{
  // Write jets and constitutents. Will use fjJets_, inputs_
  // and fjClusterSeq_

  if ( writeCompound_ ) {
    // Write jets and subjets
    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:
      throw cms::Exception("InvalidInput") << "invalid jet type in CompoundJetProducer\n";
      break;
    };
  } else if ( writeJetsWithConst_ ) {
    // Write jets and new constituents.
    writeJetsWithConstituents<reco::PFJet>( iEvent, iSetup );
  } else {
    switch( jetTypeE ) {
    case JetType::CaloJet :
      writeJets<reco::CaloJet>( iEvent, iSetup);
      break;
    case JetType::PFJet :
      writeJets<reco::PFJet>( iEvent, iSetup);
      break;
    case JetType::GenJet :
      writeJets<reco::GenJet>( iEvent, iSetup);
      break;
    case JetType::TrackJet :
      writeJets<reco::TrackJet>( iEvent, iSetup);
      break;
    case JetType::PFClusterJet :
      writeJets<reco::PFClusterJet>( iEvent, iSetup);
      break;
    case JetType::BasicJet :
      writeJets<reco::BasicJet>( iEvent, iSetup);
      break;
    default:
           throw cms::Exception("InvalidInput") << "invalid jet type in VirtualJetProducer\n";
      break;
    };
  }
  
}

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

overridevirtual

Reimplemented in cms::HTTTopJetProducer, cms::CATopJetProducer, FastjetJetProducer, SubjetFilterJetProducer, cms::CSJetProducer, CMSInsideOutJetProducer, cms::SubEventGenJetProducer, and cms::SubJetProducer.

Definition at line 264 of file VirtualJetProducer.cc.

References edm::EventID::event(), edm::Event::getByToken(), mps_fire::i, edm::EventBase::id(), LogDebug, hpstanc_transforms::max, convertSQLitetoXML_cfg::output, edm::View< T >::ptrAt(), anotherprimaryvertexanalyzer_cfi::pvCollection, edm::EventID::run(), edm::View< T >::size(), and std::swap().

Referenced by cms::SubJetProducer::produce(), CMSInsideOutJetProducer::produce(), cms::CSJetProducer::produce(), SubjetFilterJetProducer::produce(), and FastjetJetProducer::produce().

{

  // If requested, set the fastjet random seed to a deterministic function
  // of the run/lumi/event. 
  // NOTE!!! The fastjet random number sequence is a global singleton.
  // Thus, we have to create an object and get access to the global singleton
  // in order to change it. 
  if ( useDeterministicSeed_ ) {
    fastjet::GhostedAreaSpec gas;
    std::vector<int> seeds(2);
    unsigned int runNum_uint = static_cast <unsigned int> (iEvent.id().run());
    unsigned int evNum_uint = static_cast <unsigned int> (iEvent.id().event()); 
    seeds[0] = std::max(runNum_uint,minSeed_ + 3) + 3 * evNum_uint;
    seeds[1] = std::max(runNum_uint,minSeed_ + 5) + 5 * evNum_uint;
    gas.set_random_status(seeds);
  }

  LogDebug("VirtualJetProducer") << "Entered produce\n";
  //determine signal vertex2
  vertex_=reco::Jet::Point(0,0,0);
  if ( (makeCaloJet(jetTypeE) || makePFJet(jetTypeE)) &&doPVCorrection_) {
    LogDebug("VirtualJetProducer") << "Adding PV info\n";
    edm::Handle<reco::VertexCollection> pvCollection;
    iEvent.getByToken(input_vertex_token_ , pvCollection);
    if (pvCollection->size()>0) vertex_=pvCollection->begin()->position();
  }

  // For Pileup subtraction using offset correction:
  // set up geometry map
  if ( doPUOffsetCorr_ ) {
     subtractor_->setupGeometryMap(iEvent, iSetup);
  }

  // clear data
  LogDebug("VirtualJetProducer") << "Clear data\n";
  fjInputs_.clear();
  fjJets_.clear();
  inputs_.clear();  
  
  // get inputs and convert them to the fastjet format (fastjet::PeudoJet)
  edm::Handle<reco::CandidateView> inputsHandle;
  
  edm::Handle< std::vector<edm::FwdPtr<reco::PFCandidate> > > pfinputsHandleAsFwdPtr; 
  edm::Handle< std::vector<edm::FwdPtr<pat::PackedCandidate> > > packedinputsHandleAsFwdPtr; 
  
  bool isView = iEvent.getByToken(input_candidateview_token_, inputsHandle);
  if ( isView ) {
    if ( inputsHandle->size() == 0) {
      output( iEvent, iSetup );
      return;
    }
    for (size_t i = 0; i < inputsHandle->size(); ++i) {
      inputs_.push_back(inputsHandle->ptrAt(i));
    }
  } else {
    bool isPF = iEvent.getByToken(input_candidatefwdptr_token_, pfinputsHandleAsFwdPtr);
    if ( isPF ) {
      if ( pfinputsHandleAsFwdPtr->size() == 0) {
    output( iEvent, iSetup );
    return;
      }
      for (size_t i = 0; i < pfinputsHandleAsFwdPtr->size(); ++i) {
    if ( (*pfinputsHandleAsFwdPtr)[i].ptr().isAvailable() ) {
      inputs_.push_back( (*pfinputsHandleAsFwdPtr)[i].ptr() );
    }
    else if ( (*pfinputsHandleAsFwdPtr)[i].backPtr().isAvailable() ) {
      inputs_.push_back( (*pfinputsHandleAsFwdPtr)[i].backPtr() );
    }
      }
    } else {
      iEvent.getByToken(input_packedcandidatefwdptr_token_, packedinputsHandleAsFwdPtr);
      if ( packedinputsHandleAsFwdPtr->size() == 0) {
    output( iEvent, iSetup );
    return;
      }
      for (size_t i = 0; i < packedinputsHandleAsFwdPtr->size(); ++i) {
    if ( (*packedinputsHandleAsFwdPtr)[i].ptr().isAvailable() ) {
      inputs_.push_back( (*packedinputsHandleAsFwdPtr)[i].ptr() );
    }
    else if ( (*packedinputsHandleAsFwdPtr)[i].backPtr().isAvailable() ) {
      inputs_.push_back( (*packedinputsHandleAsFwdPtr)[i].backPtr() );
    }
      }
    }
  }
  LogDebug("VirtualJetProducer") << "Got inputs\n";
  
  // Convert candidates to fastjet::PseudoJets.
  // Also correct to Primary Vertex. Will modify fjInputs_
  // and use inputs_
  fjInputs_.reserve(inputs_.size());
  inputTowers();
  LogDebug("VirtualJetProducer") << "Inputted towers\n";

  // For Pileup subtraction using offset correction:
  // Subtract pedestal. 
  if ( doPUOffsetCorr_ ) {
     subtractor_->setDefinition(fjJetDefinition_);
     subtractor_->reset(inputs_,fjInputs_,fjJets_);
     subtractor_->calculatePedestal(fjInputs_); 
     subtractor_->subtractPedestal(fjInputs_);    
     LogDebug("VirtualJetProducer") << "Subtracted pedestal\n";
  }
  // Run algorithm. Will modify fjJets_ and allocate fjClusterSeq_. 
  // This will use fjInputs_
  runAlgorithm( iEvent, iSetup );

  // if ( doPUOffsetCorr_ ) {
  //    subtractor_->setAlgorithm(fjClusterSeq_);
  // }

  LogDebug("VirtualJetProducer") << "Ran algorithm\n";
  // For Pileup subtraction using offset correction:
  // Now we find jets and need to recalculate their energy,
  // mark towers participated in jet,
  // remove occupied towers from the list and recalculate mean and sigma
  // put the initial towers collection to the jet,   
  // and subtract from initial towers in jet recalculated mean and sigma of towers 
  if ( doPUOffsetCorr_ ) {
    LogDebug("VirtualJetProducer") << "Do PUOffsetCorr\n";
    vector<fastjet::PseudoJet> orphanInput;
    subtractor_->calculateOrphanInput(orphanInput);
    subtractor_->calculatePedestal(orphanInput);
    subtractor_->offsetCorrectJets();
  }
  // Write the output jets.
  // This will (by default) call the member function template
  // "writeJets", but can be overridden. 
  // this will use inputs_
  output( iEvent, iSetup );
  LogDebug("VirtualJetProducer") << "Wrote jets\n";
  
  // 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_);  

  return;
}

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

protectedpure virtual

Implemented in cms::HTTTopJetProducer, cms::CATopJetProducer, FastjetJetProducer, SubjetFilterJetProducer, CMSInsideOutJetProducer, cms::CSJetProducer, cms::SubEventGenJetProducer, and cms::SubJetProducer.

template<class T >

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

protected

function template to write out the outputs

Definition at line 723 of file VirtualJetProducer.cc.

References gather_cfg::cout, metsig::jet, eostools::move(), point, edm::Event::put(), reco::Jet::setJetArea(), and reco::writeSpecific().

{
  if ( verbosity_ >= 1 ) { 
    std::cout << "<VirtualJetProducer::writeCompoundJets (moduleLabel = " << moduleLabel_ << ")>:" << std::endl;
  }

  // get a list of output jets
  auto jetCollection = std::make_unique<reco::BasicJetCollection>();
  // get a list of output subjets
  auto subjetCollection = std::make_unique<std::vector<T>>();

  // This will store the handle for the subjets after we write them
  edm::OrphanHandle< std::vector<T> > subjetHandleAfterPut;
  // this is the mapping of subjet to hard jet
  std::vector< std::vector<int> > indices;
  // this is the list of hardjet 4-momenta
  std::vector<math::XYZTLorentzVector> p4_hardJets;
  // this is the hardjet areas
  std::vector<double> area_hardJets;

  // Loop over the hard jets
  std::vector<fastjet::PseudoJet>::const_iterator it = fjJets_.begin(),
    iEnd = fjJets_.end(),
    iBegin = fjJets_.begin();
  indices.resize( fjJets_.size() );
  for ( ; it != iEnd; ++it ) {
    fastjet::PseudoJet const & localJet = *it;
    unsigned int jetIndex = it - iBegin;
    // Get the 4-vector for the hard jet
    p4_hardJets.push_back( math::XYZTLorentzVector(localJet.px(), localJet.py(), localJet.pz(), localJet.e() ));
    double localJetArea = 0.0;
    if ( doAreaFastjet_ && localJet.has_area() ) {
      localJetArea = localJet.area();
    }
    area_hardJets.push_back( localJetArea );

    // create the subjet list
    std::vector<fastjet::PseudoJet> constituents;
    if ( it->has_pieces() ) {
      constituents = it->pieces();
    } else if ( it->has_constituents() ) {
      constituents = it->constituents();
    }

    std::vector<fastjet::PseudoJet>::const_iterator itSubJetBegin = constituents.begin(),
      itSubJet = itSubJetBegin, itSubJetEnd = constituents.end();
    for (; itSubJet != itSubJetEnd; ++itSubJet ){

      fastjet::PseudoJet const & subjet = *itSubJet;      
      if ( verbosity_ >= 1 ) {
    std::cout << "subjet #" << (itSubJet - itSubJetBegin) << ": Pt = " << subjet.pt() << ", eta = " << subjet.eta() << ", phi = " << subjet.phi() << ", mass = " << subjet.m() 
          << " (#constituents = " << subjet.constituents().size() << ")" << std::endl;
    std::vector<fastjet::PseudoJet> subjet_constituents = subjet.constituents();
    int idx_constituent = 0;
    for ( std::vector<fastjet::PseudoJet>::const_iterator constituent = subjet_constituents.begin();
          constituent != subjet_constituents.end(); ++constituent ) {
      if ( constituent->pt() < 1.e-3 ) continue; // CV: skip ghosts
      std::cout << "  constituent #" << idx_constituent << ": Pt = " << constituent->pt() << ", eta = " << constituent->eta() << ", phi = " << constituent->phi() << "," 
            << " mass = " << constituent->m() << std::endl;
      ++idx_constituent;
    }
      }

      if ( verbosity_ >= 1 ) {
    std::cout << "subjet #" << (itSubJet - itSubJetBegin) << ": Pt = " << subjet.pt() << ", eta = " << subjet.eta() << ", phi = " << subjet.phi() << ", mass = " << subjet.m() 
          << " (#constituents = " << subjet.constituents().size() << ")" << std::endl;
    std::vector<fastjet::PseudoJet> subjet_constituents = subjet.constituents();
    int idx_constituent = 0;
    for ( std::vector<fastjet::PseudoJet>::const_iterator constituent = subjet_constituents.begin();
          constituent != subjet_constituents.end(); ++constituent ) {
      if ( constituent->pt() < 1.e-3 ) continue; // CV: skip ghosts
      std::cout << " constituent #" << idx_constituent << ": Pt = " << constituent->pt() << ", eta = " << constituent->eta() << ", phi = " << constituent->phi() << "," 
            << " mass = " << constituent->m() << std::endl;
      ++idx_constituent;
    }
      }

      math::XYZTLorentzVector p4Subjet(subjet.px(), subjet.py(), subjet.pz(), subjet.e() );
      reco::Particle::Point point(0,0,0);

      // This will hold ptr's to the subjets
      std::vector<reco::CandidatePtr> subjetConstituents;

      // Get the transient subjet constituents from fastjet
      std::vector<fastjet::PseudoJet> subjetFastjetConstituents = subjet.constituents();
      std::vector<reco::CandidatePtr> constituents =
    getConstituents(subjetFastjetConstituents );    

      indices[jetIndex].push_back( subjetCollection->size() );

      // Add the concrete subjet type to the subjet list to write to event record
      T jet;
      reco::writeSpecific( jet, p4Subjet, point, constituents, iSetup);
      double subjetArea = 0.0;
      if ( doAreaFastjet_ && itSubJet->has_area() ){
    subjetArea = itSubJet->area();
      }
      jet.setJetArea( subjetArea );
      subjetCollection->push_back( jet );
    }
  }

  // put subjets into event record
  subjetHandleAfterPut = iEvent.put(std::move(subjetCollection), jetCollInstanceName_);
  
  // Now create the hard jets with ptr's to the subjets as constituents
  std::vector<math::XYZTLorentzVector>::const_iterator ip4 = p4_hardJets.begin(),
    ip4Begin = p4_hardJets.begin(),
    ip4End = p4_hardJets.end();

  for ( ; ip4 != ip4End; ++ip4 ) {
    int p4_index = ip4 - ip4Begin;
    std::vector<int> & ind = indices[p4_index];
    std::vector<reco::CandidatePtr> i_hardJetConstituents;
    // Add the subjets to the hard jet
    for( std::vector<int>::const_iterator isub = ind.begin();
     isub != ind.end(); ++isub ) {
      reco::CandidatePtr candPtr( subjetHandleAfterPut, *isub, false );
      i_hardJetConstituents.push_back( candPtr );
    }   
    reco::Particle::Point point(0,0,0);
    reco::BasicJet toput( *ip4, point, i_hardJetConstituents);
    toput.setJetArea( area_hardJets[ip4 - ip4Begin] );
    jetCollection->push_back( toput );
  }

  // put hard jets into event record
  // Store the Orphan handle for adding HTT information
  edm::OrphanHandle<reco::BasicJetCollection>  oh = iEvent.put(std::move(jetCollection));

  if (fromHTTTopJetProducer_){
    addHTTTopJetTagInfoCollection( iEvent, iSetup, oh);
  }

}

template<typename T >

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

protected

Definition at line 563 of file VirtualJetProducer.cc.

References reco::deltaR2(), SoftLeptonByDistance_cfi::distance, PVValHelper::eta, Exception, plotBeamSpotDB::first, reco::helper::VirtualJetProducerHelper::intersection(), edm::isNotFinite(), metsig::jet, fwrapper::jets, gen::k, M_PI, eostools::move(), DetIdAssociatorESProducer_cff::nEta, edm::Event::put(), fastjet::BackgroundEstimator::rho(), edm::second(), fastjet::BackgroundEstimator::set_excluded_jets(), fastjet::BackgroundEstimator::sigma(), mathSSE::sqrt(), and reco::writeSpecific().

{
  // std::cout << "writeJets " << typeid(T).name() 
  //          << (doRhoFastjet_ ? " doRhoFastjet " : "")
  //          << (doAreaFastjet_ ? " doAreaFastjet " : "")
  //          << (doAreaDiskApprox_ ? " doAreaDiskApprox " : "")
  //          << std::endl;

  if (doRhoFastjet_) {
    // declare jet collection without the two jets, 
    // for unbiased background estimation.
    std::vector<fastjet::PseudoJet> fjexcluded_jets;
    fjexcluded_jets=fjJets_;
    
    if(fjexcluded_jets.size()>2) fjexcluded_jets.resize(nExclude_);
    
    if(doFastJetNonUniform_){
      auto rhos = std::make_unique<std::vector<double>>();
      auto sigmas = std::make_unique<std::vector<double>>();
      int nEta = puCenters_.size();
      rhos->reserve(nEta);
      sigmas->reserve(nEta);
      fastjet::ClusterSequenceAreaBase const* clusterSequenceWithArea =
        dynamic_cast<fastjet::ClusterSequenceAreaBase const *> ( &*fjClusterSeq_ );

      if (clusterSequenceWithArea ==nullptr ){
    if (fjJets_.size() > 0) {
      throw cms::Exception("LogicError")<<"fjClusterSeq is not initialized while inputs are present\n ";
    }
      } else {
    for(int ie = 0; ie < nEta; ++ie){
      double eta = puCenters_[ie];
      double etamin=eta-puWidth_;
      double etamax=eta+puWidth_;
      fastjet::RangeDefinition range_rho(etamin,etamax);
      fastjet::BackgroundEstimator bkgestim(*clusterSequenceWithArea,range_rho);
      bkgestim.set_excluded_jets(fjexcluded_jets);
      rhos->push_back(bkgestim.rho());
      sigmas->push_back(bkgestim.sigma());
    }
      }
      iEvent.put(std::move(rhos),"rhos");
      iEvent.put(std::move(sigmas),"sigmas");
    }else{
      auto rho = std::make_unique<double>(0.0);
      auto sigma = std::make_unique<double>(0.0);
      double mean_area = 0;
      
      fastjet::ClusterSequenceAreaBase const* clusterSequenceWithArea =
        dynamic_cast<fastjet::ClusterSequenceAreaBase const *> ( &*fjClusterSeq_ );
      /*
    const double nemptyjets = clusterSequenceWithArea->n_empty_jets(*fjRangeDef_);
    if(( nemptyjets  < -15 ) || ( nemptyjets > fjRangeDef_->area()+ 15)) {
    edm::LogWarning("StrangeNEmtpyJets") << "n_empty_jets is : " << clusterSequenceWithArea->n_empty_jets(*fjRangeDef_) << " with range " << fjRangeDef_->description() << ".";
    }
      */
      if (clusterSequenceWithArea ==nullptr ){
    if (fjJets_.size() > 0) {
      throw cms::Exception("LogicError")<<"fjClusterSeq is not initialized while inputs are present\n ";
    }
      } else {
    clusterSequenceWithArea->get_median_rho_and_sigma(*fjRangeDef_,false,*rho,*sigma,mean_area);
    if((*rho < 0)|| (edm::isNotFinite(*rho))) {
      edm::LogError("BadRho") << "rho value is " << *rho << " area:" << mean_area << " and n_empty_jets: " << clusterSequenceWithArea->n_empty_jets(*fjRangeDef_) << " with range " << fjRangeDef_->description()
                  <<". Setting rho to rezo.";
      *rho = 0;
    }
      }
      iEvent.put(std::move(rho),"rho");
      iEvent.put(std::move(sigma),"sigma");
    }
  } // doRhoFastjet_
  
  // produce output jet collection
  
  using namespace reco;

  // allocate fjJets_.size() Ts in vector  
  auto jets = std::make_unique<std::vector<T>>(fjJets_.size());
  
  // Distance between jet centers and overlap area -- for disk-based area calculation
  using RIJ = std::pair<double,double>; 
  std::vector<RIJ>   rijStorage(fjJets_.size()*(fjJets_.size()/2));
  RIJ * rij[fjJets_.size()];
  unsigned int k=0;
  for (unsigned int ijet=0;ijet<fjJets_.size();++ijet) {
     rij[ijet] = &rijStorage[k]; k+=ijet;
  }

  float etaJ[fjJets_.size()],  phiJ[fjJets_.size()];

  auto orParam_ = 1./rParam_;
  // fill jets 
  for (unsigned int ijet=0;ijet<fjJets_.size();++ijet) {
    auto & jet = (*jets)[ijet];
    // get the fastjet jet
    const fastjet::PseudoJet& fjJet = fjJets_[ijet];
    // get the constituents from fastjet
    std::vector<fastjet::PseudoJet> const & fjConstituents = fastjet::sorted_by_pt(fjJet.constituents());
    // convert them to CandidatePtr vector
    std::vector<CandidatePtr> const & constituents = getConstituents(fjConstituents);

    // write the specifics to the jet (simultaneously sets 4-vector, vertex).
    // These are overridden functions that will call the appropriate
    // specific allocator.
    writeSpecific(jet,
                  Particle::LorentzVector(fjJet.px(),
                                          fjJet.py(),
                                          fjJet.pz(),
                                          fjJet.E()),
                  vertex_,
                  constituents, iSetup);
    phiJ[ijet] = jet.phi();
    etaJ[ijet] = jet.eta();
  }

   // calcuate the jet area
  for (unsigned int ijet=0;ijet<fjJets_.size();++ijet) {
    // calcuate the jet area
    double jetArea=0.0;
    // get the fastjet jet
    const auto & fjJet = fjJets_[ijet];
    if ( doAreaFastjet_ && fjJet.has_area() ) {
      jetArea = fjJet.area();
    }
    else if ( doAreaDiskApprox_ ) {
      // Here it is assumed that fjJets_ is in decreasing order of pT, 
      // which should happen in FastjetJetProducer::runAlgorithm() 
      jetArea   = M_PI;
        RIJ *  distance  = rij[ijet];
        for (unsigned jJet = 0; jJet < ijet; ++jJet) {
          distance[jJet].first      = std::sqrt(reco::deltaR2(etaJ[ijet],phiJ[ijet], etaJ[jJet],phiJ[jJet]))*orParam_;
          distance[jJet].second = reco::helper::VirtualJetProducerHelper::intersection(distance[jJet].first);
          jetArea            -=distance[jJet].second;
          for (unsigned kJet = 0; kJet < jJet; ++kJet) {
            jetArea          += reco::helper::VirtualJetProducerHelper::intersection(distance[jJet].first, distance[kJet].first, rij[jJet][kJet].first, 
                                                                                     distance[jJet].second, distance[kJet].second, rij[jJet][kJet].second);
          } // end loop over harder jets
        } // end loop over harder jets
      jetArea  *= (rParam_*rParam_);
    } 
    auto & jet = (*jets)[ijet]; 
    jet.setJetArea (jetArea);
    
    if(doPUOffsetCorr_){
      jet.setPileup(subtractor_->getPileUpEnergy(ijet));
    }else{
      jet.setPileup (0.0);
    }
        
   // std::cout << "area " << ijet << " " << jetArea << " " << Area<T>::get(jet) << std::endl;
   // std::cout << "JetVI " << ijet << ' '<< jet.pt() << " " << jet.et() << ' '<< jet.energy() << ' '<< jet.mass() << std::endl;

  }
  // put the jets in the collection
  iEvent.put(std::move(jets),jetCollInstanceName_);
}

template<class T >

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

protected

function template to write out the outputs

Definition at line 861 of file VirtualJetProducer.cc.

References gather_cfg::cout, hcalTTPDigis_cfi::id, metsig::jet, eostools::move(), point, edm::Event::put(), reco::Jet::setJetArea(), reco::LeafCandidate::setP4(), reco::PFCandidate::setSourceCandidatePtr(), reco::PFCandidate::translatePdgIdToType(), and reco::writeSpecific().

{
  if ( verbosity_ >= 1 ) {
    std::cout << "<VirtualJetProducer::writeJetsWithConstituents (moduleLabel = " << moduleLabel_ << ")>:" << std::endl;
  }

  // get a list of output jets  MV: make this compatible with template
  auto jetCollection = std::make_unique<reco::PFJetCollection>();
  
  // this is the mapping of jet to constituents
  std::vector< std::vector<int> > indices;
  // this is the list of jet 4-momenta
  std::vector<math::XYZTLorentzVector> p4_Jets;
  // this is the jet areas
  std::vector<double> area_Jets;

  // get a list of output constituents
  auto constituentCollection = std::make_unique<reco::PFCandidateCollection>();
  
  // This will store the handle for the constituents after we write them
  edm::OrphanHandle<reco::PFCandidateCollection> constituentHandleAfterPut;
    
  // Loop over the jets and extract constituents
  std::vector<fastjet::PseudoJet> constituentsSub;
  std::vector<fastjet::PseudoJet>::const_iterator it = fjJets_.begin(),
    iEnd = fjJets_.end(),
    iBegin = fjJets_.begin();
  indices.resize( fjJets_.size() );

  for ( ; it != iEnd; ++it ) {
    fastjet::PseudoJet const & localJet = *it;
    unsigned int jetIndex = it - iBegin;
    // Get the 4-vector for the hard jet
    p4_Jets.push_back( math::XYZTLorentzVector(localJet.px(), localJet.py(), localJet.pz(), localJet.e() ));
    double localJetArea = 0.0;
    if ( doAreaFastjet_ && localJet.has_area() ) {
      localJetArea = localJet.area();
    }
    area_Jets.push_back( localJetArea );

    // create the constituent list
    std::vector<fastjet::PseudoJet> constituents,ghosts;
    if ( it->has_pieces() )
      constituents = it->pieces();
    else if ( it->has_constituents() )
      fastjet::SelectorIsPureGhost().sift(it->constituents(), ghosts, constituents); //filter out ghosts

    //loop over constituents of jet (can be subjets or normal constituents)
    indices[jetIndex].reserve(constituents.size());
    constituentsSub.reserve(constituentsSub.size()+constituents.size());
    for (fastjet::PseudoJet const& constit : constituents) {
      indices[jetIndex].push_back( constituentsSub.size() );
      constituentsSub.push_back(constit);
    }
  }

  //Loop over constituents and store in the event
  static const reco::PFCandidate dummySinceTranslateIsNotStatic;
  for (fastjet::PseudoJet const& constit : constituentsSub) {
    auto orig = inputs_[constit.user_index()];
    auto id = dummySinceTranslateIsNotStatic.translatePdgIdToType(orig->pdgId());
    reco::PFCandidate pCand( reco::PFCandidate(orig->charge(), orig->p4(), id) );
    math::XYZTLorentzVector pVec;
    pVec.SetPxPyPzE(constit.px(),constit.py(),constit.pz(),constit.e());
    pCand.setP4(pVec);
    pCand.setSourceCandidatePtr( orig->sourceCandidatePtr(0) );
    constituentCollection->push_back(pCand);
  }
  // put constituents into event record
  constituentHandleAfterPut = iEvent.put(std::move(constituentCollection), jetCollInstanceName_ );

  // Now create the jets with ptr's to the constituents
  std::vector<math::XYZTLorentzVector>::const_iterator ip4 = p4_Jets.begin(),
    ip4Begin = p4_Jets.begin(),
    ip4End = p4_Jets.end();

  for ( ; ip4 != ip4End; ++ip4 ) {
    int p4_index = ip4 - ip4Begin;
    std::vector<int> & ind = indices[p4_index];
    std::vector<reco::CandidatePtr> i_jetConstituents;
    // Add the constituents to the jet
    for( std::vector<int>::const_iterator iconst = ind.begin(); iconst != ind.end(); ++iconst ) {
      reco::CandidatePtr candPtr( constituentHandleAfterPut, *iconst, false );
      i_jetConstituents.push_back( candPtr );
    }
    if(i_jetConstituents.size()>0) { //only keep jets which have constituents after subtraction
      reco::Particle::Point point(0,0,0);
      reco::PFJet jet;
      reco::writeSpecific(jet,*ip4,point,i_jetConstituents,iSetup);
      jet.setJetArea( area_Jets[ip4 - ip4Begin] );
      jetCollection->emplace_back( jet );
    }
  }

  // put jets into event record
  iEvent.put(std::move(jetCollection));
}

Member Data Documentation

int VirtualJetProducer::activeAreaRepeats_

protected

Definition at line 190 of file VirtualJetProducer.h.

std::auto_ptr<AnomalousTower> VirtualJetProducer::anomalousTowerDef_

private

Definition at line 225 of file VirtualJetProducer.h.

bool VirtualJetProducer::doAreaDiskApprox_

protected

Definition at line 182 of file VirtualJetProducer.h.

bool VirtualJetProducer::doAreaFastjet_

protected

Definition at line 180 of file VirtualJetProducer.h.

Referenced by cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

bool VirtualJetProducer::doFastJetNonUniform_

protected

Definition at line 185 of file VirtualJetProducer.h.

bool VirtualJetProducer::doPUOffsetCorr_

protected

Definition at line 194 of file VirtualJetProducer.h.

bool VirtualJetProducer::doPVCorrection_

protected

Definition at line 173 of file VirtualJetProducer.h.

bool VirtualJetProducer::doRhoFastjet_

protected

Definition at line 184 of file VirtualJetProducer.h.

Referenced by cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

ActiveAreaSpecPtr VirtualJetProducer::fjActiveArea_

protected

Definition at line 202 of file VirtualJetProducer.h.

AreaDefinitionPtr VirtualJetProducer::fjAreaDefinition_

protected

Definition at line 203 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::FastjetJetProducer(), cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

ClusterSequencePtr VirtualJetProducer::fjClusterSeq_

protected

Definition at line 199 of file VirtualJetProducer.h.

Referenced by cms::CSJetProducer::produce(), FastjetJetProducer::produce(), FastjetJetProducer::produceTrackJets(), cms::SubEventGenJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

std::vector<fastjet::PseudoJet> VirtualJetProducer::fjInputs_

protected

Definition at line 205 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::produce(), FastjetJetProducer::produceTrackJets(), cms::SubJetProducer::runAlgorithm(), cms::SubEventGenJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), CMSInsideOutJetProducer::runAlgorithm(), SubjetFilterJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

JetDefPtr VirtualJetProducer::fjJetDefinition_

protected

Definition at line 200 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

std::vector<fastjet::PseudoJet> VirtualJetProducer::fjJets_

protected

Definition at line 206 of file VirtualJetProducer.h.

Referenced by cms::HTTTopJetProducer::addHTTTopJetTagInfoCollection(), cms::SubEventGenJetProducer::produce(), FastjetJetProducer::produceTrackJets(), cms::SubEventGenJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), CMSInsideOutJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

PluginPtr VirtualJetProducer::fjPlugin_

protected

Definition at line 201 of file VirtualJetProducer.h.

RangeDefPtr VirtualJetProducer::fjRangeDef_

protected

Definition at line 204 of file VirtualJetProducer.h.

bool VirtualJetProducer::fromHTTTopJetProducer_

protected

Definition at line 222 of file VirtualJetProducer.h.

Referenced by cms::HTTTopJetProducer::HTTTopJetProducer().

double VirtualJetProducer::ghostArea_

protected

Definition at line 191 of file VirtualJetProducer.h.

double VirtualJetProducer::ghostEtaMax_

protected

Definition at line 189 of file VirtualJetProducer.h.

edm::EDGetTokenT<std::vector<edm::FwdPtr<reco::PFCandidate> > > VirtualJetProducer::input_candidatefwdptr_token_

private

Definition at line 229 of file VirtualJetProducer.h.

edm::EDGetTokenT<reco::CandidateView> VirtualJetProducer::input_candidateview_token_

private

Definition at line 228 of file VirtualJetProducer.h.

edm::EDGetTokenT<std::vector<edm::FwdPtr<pat::PackedCandidate> > > VirtualJetProducer::input_packedcandidatefwdptr_token_

private

Definition at line 230 of file VirtualJetProducer.h.

edm::EDGetTokenT<reco::VertexCollection> VirtualJetProducer::input_vertex_token_

protected

Definition at line 233 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::produceTrackJets().

double VirtualJetProducer::inputEMin_

protected

Definition at line 171 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::inputTowers().

double VirtualJetProducer::inputEtMin_

protected

Definition at line 170 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::inputTowers().

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

protected

Definition at line 197 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::inputTowers(), cms::SubEventGenJetProducer::produce(), FastjetJetProducer::produceTrackJets(), SubjetFilterJetProducer::writeCompoundJets(), and cms::CompoundJetProducer::writeCompoundJets().

std::string VirtualJetProducer::jetAlgorithm_

protected

Definition at line 168 of file VirtualJetProducer.h.

std::string VirtualJetProducer::jetCollInstanceName_

protected

Definition at line 213 of file VirtualJetProducer.h.

Referenced by cms::CompoundJetProducer::writeCompoundJets().

double VirtualJetProducer::jetPtMin_

protected

Definition at line 172 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), and FastjetJetProducer::runAlgorithm().

std::string VirtualJetProducer::jetType_

protected

Definition at line 167 of file VirtualJetProducer.h.

JetType::Type VirtualJetProducer::jetTypeE

protected

Definition at line 55 of file VirtualJetProducer.h.

Referenced by SubjetFilterJetProducer::output(), cms::CompoundJetProducer::output(), and FastjetJetProducer::produce().

unsigned int VirtualJetProducer::maxInputs_

protected

Definition at line 177 of file VirtualJetProducer.h.

unsigned int VirtualJetProducer::minSeed_

protected

Definition at line 219 of file VirtualJetProducer.h.

std::string VirtualJetProducer::moduleLabel_

protected

Definition at line 164 of file VirtualJetProducer.h.

Referenced by Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), Modules.ESProducer::nameInProcessDesc_(), and SubjetFilterJetProducer::SubjetFilterJetProducer().

unsigned int VirtualJetProducer::nExclude_

protected

Definition at line 211 of file VirtualJetProducer.h.

std::vector<double> VirtualJetProducer::puCenters_

protected

Definition at line 209 of file VirtualJetProducer.h.

std::string VirtualJetProducer::puSubtractorName_

protected

Definition at line 195 of file VirtualJetProducer.h.

double VirtualJetProducer::puWidth_

protected

Definition at line 210 of file VirtualJetProducer.h.

bool VirtualJetProducer::restrictInputs_

protected

Definition at line 176 of file VirtualJetProducer.h.

double VirtualJetProducer::rhoEtaMax_

protected

Definition at line 188 of file VirtualJetProducer.h.

double VirtualJetProducer::rParam_

protected

Definition at line 169 of file VirtualJetProducer.h.

edm::InputTag VirtualJetProducer::src_

protected

Definition at line 165 of file VirtualJetProducer.h.

Referenced by cms::CATopJetProducer::CATopJetProducer(), and cms::SubEventGenJetProducer::SubEventGenJetProducer().

edm::InputTag VirtualJetProducer::srcPVs_

protected

Definition at line 166 of file VirtualJetProducer.h.

boost::shared_ptr<PileUpSubtractor> VirtualJetProducer::subtractor_

protected

Definition at line 216 of file VirtualJetProducer.h.

bool VirtualJetProducer::useDeterministicSeed_

protected

Definition at line 218 of file VirtualJetProducer.h.

bool VirtualJetProducer::useExplicitGhosts_

protected

Definition at line 181 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::FastjetJetProducer().

int VirtualJetProducer::verbosity_

protected

Definition at line 221 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::runAlgorithm().

reco::Particle::Point VirtualJetProducer::vertex_

protected

Definition at line 198 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::produceTrackJets(), and cms::SubEventGenJetProducer::runAlgorithm().

double VirtualJetProducer::voronoiRfact_

protected

Definition at line 186 of file VirtualJetProducer.h.

Referenced by cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

bool VirtualJetProducer::writeCompound_

protected

Definition at line 214 of file VirtualJetProducer.h.

bool VirtualJetProducer::writeJetsWithConst_

protected

Definition at line 215 of file VirtualJetProducer.h.