Inheritance diagram for VoronoiBackgroundProducer:

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

	~VoronoiBackgroundProducer ()

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

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_

edm::InputTag	src_

VoronoiAlgorithm *	voronoi_

std::vector < reco::VoronoiBackground >	vvm

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

typedef WorkerT< EDProducer >	WorkerType

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

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)

Detailed Description

Definition at line 26 of file VoronoiBackgroundProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 63 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")),
    isCalo_(iConfig.getParameter<bool>("isCalo")),
    etaBins_(iConfig.getParameter<int>("etaBins")),
    fourierOrder_(iConfig.getParameter<int>("fourierOrder"))
 {
 
    src_ = iConfig.getParameter<edm::InputTag>("src");
    //register your products
 
    produces<reco::VoronoiMap>();
    produces<std::vector<float> >();
 }

VoronoiBackgroundProducer::~VoronoiBackgroundProducer ( )

Definition at line 82 of file VoronoiBackgroundProducer.cc.

83 {

84 }

Member Function Documentation

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

static

Definition at line 143 of file VoronoiBackgroundProducer.cc.

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

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

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

privatevirtual

Implements edm::EDProducer.

Definition at line 93 of file VoronoiBackgroundProducer.cc.

References VoronoiAlgorithm::clear(), data, doEqualize_, equalizeR_, equalizeThreshold0_, equalizeThreshold1_, edm::helper::Filler< Map >::fill(), edm::Event::getByLabel(), i, edm::helper::Filler< Map >::insert(), isCalo_, edm::EventBase::isRealData(), LogDebug, VoronoiAlgorithm::particle_area(), VoronoiAlgorithm::perp_fourier(), VoronoiAlgorithm::push_back_particle(), edm::Event::put(), mathSSE::sqrt(), src_, VoronoiAlgorithm::subtracted_equalized_perp(), VoronoiAlgorithm::subtracted_unequalized_perp(), voronoi_, and vvm.

 {
    using namespace edm;
    if(voronoi_ == 0){
      bool data = iEvent.isRealData();
      voronoi_ = new VoronoiAlgorithm(equalizeR_,data,isCalo_,std::pair<double, double>(equalizeThreshold0_,equalizeThreshold1_),doEqualize_);
    }
 
    voronoi_->clear();
    vvm.clear();
 
    edm::Handle<reco::CandidateView> inputsHandle;
    iEvent.getByLabel(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();
    std::auto_ptr<std::vector<float> > vnout(new std::vector<float>(voronoi_vn.begin(), voronoi_vn.end()));
    std::auto_ptr<reco::VoronoiMap> mapout(new 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(vnout);
    iEvent.put(mapout);
  
 }