SingleCellClusterAlgo.cc

Go to the documentation of this file.

#include "L1Trigger/L1THGCal/interface/HGCalTriggerBackendAlgorithmBase.h"
#include "L1Trigger/L1THGCal/interface/fe_codecs/HGCalBestChoiceCodec.h"
#include "DataFormats/ForwardDetId/interface/HGCalDetId.h"

#include "DataFormats/L1THGCal/interface/HGCalCluster.h"

using namespace HGCalTriggerBackend;

class SingleCellClusterAlgo : public Algorithm<HGCalBestChoiceCodec> 
{
    public:

        SingleCellClusterAlgo(const edm::ParameterSet& conf):
            Algorithm<HGCalBestChoiceCodec>(conf),
            cluster_product_( new l1t::HGCalClusterBxCollection ){}

        virtual void setProduces(edm::EDProducer& prod) const override final 
        {
            prod.produces<l1t::HGCalClusterBxCollection>(name());
        }

        virtual void run(const l1t::HGCFETriggerDigiCollection& coll,
                const std::unique_ptr<HGCalTriggerGeometryBase>& geom) override final;

        virtual void putInEvent(edm::Event& evt) override final 
        {
            evt.put(std::move(cluster_product_),name());
        }

        virtual void reset() override final 
        {
            cluster_product_.reset( new l1t::HGCalClusterBxCollection );
        }

    private:
        std::unique_ptr<l1t::HGCalClusterBxCollection> cluster_product_;

};

/*****************************************************************/
void SingleCellClusterAlgo::run(const l1t::HGCFETriggerDigiCollection& coll,
        const std::unique_ptr<HGCalTriggerGeometryBase>& geom) 
/*****************************************************************/
{
    for( const auto& digi : coll ) 
    {
        HGCalBestChoiceCodec::data_type data;
        data.reset();
        const HGCalDetId& moduleId = digi.getDetId<HGCalDetId>();
        digi.decode(codec_, data);
        int i = 0;
        for(const auto& value : data.payload)
        {
            if(value>0)
            {
                GlobalPoint point = geom->getModulePosition(moduleId);
                math::PtEtaPhiMLorentzVector p4((double)value/cosh(point.eta()), point.eta(), point.phi(), 0.);
                // index in module stored as hwEta
                l1t::HGCalCluster cluster( 
                        reco::LeafCandidate::LorentzVector(),
                        value, i, 0);
                cluster.setP4(p4);
                cluster.setModule(moduleId.wafer());
                cluster.setLayer(moduleId.layer());
                cluster.setSubDet(moduleId.subdetId());
                cluster_product_->push_back(0,cluster);
            }
            i++;
        }

    }
}

DEFINE_EDM_PLUGIN(HGCalTriggerBackendAlgorithmFactory, 
        SingleCellClusterAlgo,
        "SingleCellClusterAlgo");