Inheritance diagram for VoronoiBackgroundProducer:

Public Member Functions
	VoronoiBackgroundProducer (const edm::ParameterSet &)

	~VoronoiBackgroundProducer ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

	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 Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
bool	doEqualize_

double	equalizeR_

double	equalizeThreshold0_

double	equalizeThreshold1_

int	etaBins_

int	fourierOrder_

bool	isCalo_

bool	jetCorrectorFormat_

edm::EDGetTokenT< reco::CandidateView >	src_

std::string	tableLabel_

bool	useTextTable_

VoronoiAlgorithm *	voronoi_

std::vector< reco::VoronoiBackground >	vvm

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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)

Detailed Description

Definition at line 37 of file VoronoiBackgroundProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 77 of file VoronoiBackgroundProducer.cc.

References edm::ParameterSet::getParameter(), and src_.

                                                                                   :
    voronoi_(0),
    doEqualize_(iConfig.getParameter<bool>("doEqualize")),
    equalizeThreshold0_(iConfig.getParameter<double>("equalizeThreshold0")),
    equalizeThreshold1_(iConfig.getParameter<double>("equalizeThreshold1")),
    equalizeR_(iConfig.getParameter<double>("equalizeR")),
    useTextTable_(iConfig.getParameter<bool>("useTextTable")),
    jetCorrectorFormat_(iConfig.getParameter<bool>("jetCorrectorFormat")),
    isCalo_(iConfig.getParameter<bool>("isCalo")),
    tableLabel_(iConfig.getParameter<std::string>("tableLabel")),
    etaBins_(iConfig.getParameter<int>("etaBins")),
    fourierOrder_(iConfig.getParameter<int>("fourierOrder"))
 {
 
    src_ = consumes<reco::CandidateView>(iConfig.getParameter<edm::InputTag>("src"));
    //register your products
 
    produces<reco::VoronoiMap>();
    produces<std::vector<float> >();
 }

VoronoiBackgroundProducer::~VoronoiBackgroundProducer ( )

Definition at line 99 of file VoronoiBackgroundProducer.cc.

100 {

101 }

Member Function Documentation

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

static

Definition at line 191 of file VoronoiBackgroundProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, and edm::ParameterSetDescription::setUnknown().

                                                                                       {
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

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

privatevirtual

Definition at line 110 of file VoronoiBackgroundProducer.cc.

References HiRecoJets_cff::bkg, VoronoiAlgorithm::clear(), doEqualize_, equalizeR_, equalizeThreshold0_, equalizeThreshold1_, objects.autophobj::filler, edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, isCalo_, edm::EventBase::isRealData(), jetCorrectorFormat_, JetCorrectorParametersCollection::L1Offset, LogDebug, eostools::move(), NULL, VoronoiAlgorithm::particle_area(), jetCorrFactors_cfi::payload, VoronoiAlgorithm::perp_fourier(), edm::ESHandle< T >::product(), VoronoiAlgorithm::push_back_particle(), edm::Event::put(), edm::View< T >::size(), mathSSE::sqrt(), src_, VoronoiAlgorithm::subtracted_equalized_perp(), VoronoiAlgorithm::subtracted_unequalized_perp(), tableLabel_, useTextTable_, UETable::values, voronoi_, and vvm.

 {
    using namespace edm;
    if(voronoi_ == 0){
     UECalibration *ue = NULL;
     if (useTextTable_) {
      const bool isData = iEvent.isRealData();
      const char *calibrationFile = NULL;
       if(isCalo_){
         if(isData) calibrationFile = "RecoHI/HiJetAlgos/data/ue_calibrations_hermite_calo_data.txt";
         if(!isData) calibrationFile = "RecoHI/HiJetAlgos/data/ue_calibrations_hermite_calo_mc.txt";
       }else{
         if(isData) calibrationFile = "RecoHI/HiJetAlgos/data/ue_calibrations_hermite_pf_data.txt";
         if(!isData) calibrationFile = "RecoHI/HiJetAlgos/data/ue_calibrations_hermite_pf_mc.txt";
       }
     ue = new UECalibration(calibrationFile);
     }
     else if (jetCorrectorFormat_) {
         edm::ESHandle<JetCorrectorParametersCollection> ueHandle;
 
         iSetup.get<JetCorrectionsRecord>().get(tableLabel_,ueHandle);
 
         const JetCorrectorParametersCollection *payload = ueHandle.product();
         std::vector<float> ue_vec = (*payload)[JetCorrectorParametersCollection::L1Offset].record(0).parameters();
 
         ue = new UECalibration(ue_vec);
     }
     else {
      edm::ESHandle<UETable> ueHandle;
 
      iSetup.get<HeavyIonUERcd>().get(tableLabel_,ueHandle);
 
      const UETable *ueTable = ueHandle.product();
      ue = new UECalibration(ueTable->values);
     }
 
      voronoi_ = new VoronoiAlgorithm(ue,equalizeR_,std::pair<double, double>(equalizeThreshold0_,equalizeThreshold1_),doEqualize_);
    }
 
    voronoi_->clear();
    vvm.clear();
 
    edm::Handle<reco::CandidateView> inputsHandle;
    iEvent.getByToken(src_,inputsHandle);
 
    for(unsigned int i = 0; i < inputsHandle->size(); ++i){
       reco::CandidateViewRef ref(inputsHandle,i);
       voronoi_->push_back_particle(ref->pt(),ref->eta(),ref->phi(),0);
    }
 
    std::vector<double> subtracted_momenta = voronoi_->subtracted_unequalized_perp();
    std::vector<double> equalized_momenta = voronoi_->subtracted_equalized_perp();
    std::vector<double> particle_area = voronoi_->particle_area();
    std::vector<double> voronoi_vn = voronoi_->perp_fourier();
    auto vnout = std::make_unique<std::vector<float>>(voronoi_vn.begin(), voronoi_vn.end());
    auto mapout = std::make_unique<reco::VoronoiMap>();
    reco::VoronoiMap::Filler filler(*mapout);
 
    for(unsigned int i = 0; i < inputsHandle->size(); ++i){
       reco::CandidateViewRef ref(inputsHandle,i);
       const double pre_eq_pt = subtracted_momenta[i];
       const double post_eq_pt = equalized_momenta[i];
       const double area = particle_area[i];
       const double mass_square = ref->massSqr();
       const double pre_eq_mt = sqrt(mass_square + pre_eq_pt * pre_eq_pt);
       const double post_eq_mt = sqrt(mass_square + post_eq_pt * post_eq_pt);
 
       reco::VoronoiBackground bkg(pre_eq_pt,post_eq_pt,pre_eq_mt,post_eq_mt,area);
       LogDebug("VoronoiBackgroundProducer")<<"Subtraction --- oldpt : "<<ref->pt()<<" --- newpt : "<<post_eq_pt<<endl;
       vvm.push_back(bkg);
    }
 
    filler.insert(inputsHandle,vvm.begin(),vvm.end());
    filler.fill();
    iEvent.put(std::move(vnout));
    iEvent.put(std::move(mapout));
  
 }