HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA > Class Template Reference

Inheritance diagram for HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >:

Public Member Functions
	HGCalTriggerSimCluster (const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
virtual void	putInEvent (edm::Event &evt) override final
virtual void	reset () override final
virtual void	run (const l1t::HGCFETriggerDigiCollection &coll, const edm::EventSetup &es, edm::Event &evt)
virtual void	setProduces (edm::stream::EDProducer<> &prod) const override final
Public Member Functions inherited from HGCalTriggerBackend::Algorithm< FECODEC >
	Algorithm (const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
virtual void	setGeometry (const HGCalTriggerGeometryBase *const geom) override final
Public Member Functions inherited from HGCalTriggerBackendAlgorithmBase
	HGCalTriggerBackendAlgorithmBase (const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
const std::string &	name () const
virtual	~HGCalTriggerBackendAlgorithmBase ()

Private Member Functions
void	addToCluster (std::unordered_map< uint64_t, std::pair< int, l1t::HGCalCluster > > &cluster_container, uint64_t pid, int pdgid, float energy, float eta, float phi)
void	addToClusterShapes (std::unordered_map< uint64_t, std::pair< int, l1t::HGCalCluster > > &cluster_container, uint64_t pid, int pdgid, float energy, float eta, float phi, float r=0.0)

Private Attributes
HGCalTriggerCellCalibration	calibration_
std::unique_ptr< l1t::HGCalClusterBxCollection >	cluster_product_
std::string	HGCalEESensitive_
std::string	HGCalHESiliconSensitive_
edm::ESHandle< HGCalTopology >	hgceeTopoHandle_
edm::ESHandle< HGCalTopology >	hgchefTopoHandle_
edm::EDGetToken	inputbh_
edm::EDGetToken	inputee_
edm::EDGetToken	inputfh_
edm::Handle< std::vector< SimCluster > >	sim_handle_
edm::EDGetTokenT< std::vector< SimCluster > >	sim_token_

Additional Inherited Members
Protected Attributes inherited from HGCalTriggerBackend::Algorithm< FECODEC >
FECODEC	codec_
Protected Attributes inherited from HGCalTriggerBackendAlgorithmBase
const HGCalTriggerGeometryBase *	geometry_

Detailed Description

template<typename FECODEC, typename DATA>
class HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >

Definition at line 52 of file HGCalSimCluster.cc.

Constructor & Destructor Documentation

template<typename FECODEC , typename DATA >

HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::HGCalTriggerSimCluster ( const edm::ParameterSet &  conf, edm::ConsumesCollector &  cc  )

inline

Definition at line 117 of file HGCalSimCluster.cc.

References edm::ConsumesCollector::consumes(), and edm::ParameterSet::getParameter().

                                                                                        : 
                            Algorithm<FECODEC>(conf,cc),
                            HGCalEESensitive_(conf.getParameter<std::string>("HGCalEESensitive_tag")),
                            HGCalHESiliconSensitive_(conf.getParameter<std::string>("HGCalHESiliconSensitive_tag")),
                            calibration_(conf.getParameterSet("calib_parameters"))
            { 
                sim_token_ = cc.consumes< std::vector< SimCluster > >(conf.getParameter<edm::InputTag>("simcollection")); 
                inputee_ = cc.consumes<edm::PCaloHitContainer>( conf.getParameter<edm::InputTag>("simhitsee"));
                inputfh_ = cc.consumes<edm::PCaloHitContainer>( conf.getParameter<edm::InputTag>("simhitsfh"));
                // inputbh_ = cc.consumes<edm::PCaloHitContainer>(conf.getParameter<edm::InputTag>("g4SimHits:HGCHitsHEback")); 
                edm::LogWarning("HGCalTriggerSimCluster") <<"WARNING: BH simhits not loaded";
            }

Member Function Documentation

template<typename FECODEC , typename DATA >

void HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::addToCluster ( std::unordered_map< uint64_t, std::pair< int, l1t::HGCalCluster > > &  cluster_container, uint64_t  pid, int  pdgid, float  energy, float  eta, float  phi  )

inlineprivate

Definition at line 80 of file HGCalSimCluster.cc.

References stringResolutionProvider_cfi::eta, JetChargeProducer_cfi::exp, ResonanceBuilder::mass, p4, EnergyCorrector::pt, reco::return(), edm::second(), and lumiQTWidget::t.

Referenced by HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::run().

            {

                auto pair = cluster_container.emplace(pid, std::pair<int,l1t::HGCalCluster>(0,l1t::HGCalCluster() ) ) ;
                auto iterator = pair.first;
                math::PtEtaPhiMLorentzVectorD p4;
                p4.SetPt ( iterator -> second . second . pt()   ) ;
                p4.SetEta( iterator -> second . second . eta()  ) ;
                p4.SetPhi( iterator -> second . second . phi()  ) ;
                p4.SetM  ( iterator -> second . second . mass() ) ;
                math::PtEtaPhiMLorentzVectorD pp4; 
                float t = std::exp (- eta);
                pp4.SetPt ( energy * (2*t)/(1+t*t)  ) ;
                pp4.SetEta( eta ) ;
                pp4.SetPhi( phi ) ;
                pp4.SetM  (  0  ) ;
                p4 += pp4;
                iterator -> second . second . setP4(p4);
                //iterator -> second . second . shapes().Add( energy,eta,phi,r); // last is r, for 3d clusters
                return ;
            }

template<typename FECODEC , typename DATA >

void HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::addToClusterShapes ( std::unordered_map< uint64_t, std::pair< int, l1t::HGCalCluster > > &  cluster_container, uint64_t  pid, int  pdgid, float  energy, float  eta, float  phi, float  r = 0.0  )

inlineprivate

Definition at line 74 of file HGCalSimCluster.cc.

References alignCSCRings::r, and edm::second().

Referenced by HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::run().

                                                                                                                                                                                         {
                auto pair = cluster_container.emplace(pid, std::pair<int,l1t::HGCalCluster>(0,l1t::HGCalCluster() ) ) ;
                auto iterator = pair.first;
                iterator -> second . second . shapes().Add( energy,eta,phi,r); // last is r, for 3d clusters
            }

template<typename FECODEC , typename DATA >

virtual void HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::putInEvent ( edm::Event &  evt )

inlinefinaloverridevirtual

Implements HGCalTriggerBackendAlgorithmBase.

Definition at line 137 of file HGCalSimCluster.cc.

References eostools::move(), and HGCalTriggerBackendAlgorithmBase::name().

            {
                evt.put(std::move(cluster_product_),name());
            }

template<typename FECODEC , typename DATA >

virtual void HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::reset ( void  )

inlinefinaloverridevirtual

Implements HGCalTriggerBackendAlgorithmBase.

Definition at line 143 of file HGCalSimCluster.cc.

            {
                cluster_product_.reset( new l1t::HGCalClusterBxCollection );
            }

template<typename FECODEC , typename DATA >

virtual void HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::run ( const l1t::HGCFETriggerDigiCollection &  coll, const edm::EventSetup &  es, edm::Event &  evt  )

inlinevirtual

FIXME TODO

Implements HGCalTriggerBackendAlgorithmBase.

Definition at line 149 of file HGCalSimCluster.cc.

References HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::addToCluster(), HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::addToClusterShapes(), HGCalTriggerCellCalibration::calibrateInGeV(), HGCalTriggerBackend::Algorithm< FECODEC >::codec_, cscdqm::DATA, data, HGCalTopology::dddConstants(), HGCalTopology::detectorType(), stringResolutionProvider_cfi::eta, Exception, dedxEstimators_cff::fraction, HGCalTriggerBackendAlgorithmBase::geometry_, edm::EventSetup::get(), edm::Event::getByToken(), HGCalTriggerGeometryBase::getCellsFromTriggerCell(), photonIsolationHIProducer_cfi::hf, HGCEE, HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::hgceeTopoHandle_, HGCHEB, HGCHEF, HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::hgchefTopoHandle_, edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, reco::LeafCandidate::p4(), HadronAndPartonSelector_cfi::particles, sysUtil::pid, HGCalDDDConstants::simToReco(), DetId::subdetId(), and HGCalTestNumbering::unpackHexagonIndex().

            {
                //0.5. Get Digis, construct a map, detid -> energy
                
                edm::Handle<edm::PCaloHitContainer> ee_simhits_h;
                evt.getByToken(inputee_,ee_simhits_h);

                edm::Handle<edm::PCaloHitContainer> fh_simhits_h;
                evt.getByToken(inputfh_,fh_simhits_h);

                edm::Handle<edm::PCaloHitContainer> bh_simhits_h;
                //evt.getByToken(inputbh_,bh_simhits_h);

                
                if (not ee_simhits_h.isValid()){
                       throw cms::Exception("ContentError")<<"[HGCalTriggerSimCluster]::[run]::[ERROR] EE Digis from HGC not available"; 
                }
                if (not fh_simhits_h.isValid()){
                       throw cms::Exception("ContentError")<<"[HGCalTriggerSimCluster]::[run]::[ERROR] FH Digis from HGC not available"; 
                }

                std::unordered_map<uint64_t, double> hgc_simhit_energy;

                edm::ESHandle<HGCalTopology> topo_ee, topo_fh;
                es.get<IdealGeometryRecord>().get("HGCalEESensitive",topo_ee);
                es.get<IdealGeometryRecord>().get("HGCalHESiliconSensitive",topo_fh);

                if (ee_simhits_h.isValid()) 
                {
                    int layer=0,cell=0, sec=0, subsec=0, zp=0,subdet=0;
                    ForwardSubdetector mysubdet;
                    for (const auto& simhit : *ee_simhits_h)
                    {
                        HGCalTestNumbering::unpackHexagonIndex(simhit.id(), subdet, zp, layer, sec, subsec, cell); 
                        mysubdet = (ForwardSubdetector)(subdet);
                        std::pair<int,int> recoLayerCell = topo_ee->dddConstants().simToReco(cell,layer,sec,topo_ee->detectorType());
                        cell  = recoLayerCell.first;
                        layer = recoLayerCell.second;
                        if (layer<0 || cell<0) {
                          continue;
                        }
                        unsigned recoCell = HGCalDetId(mysubdet,zp,layer,subsec,sec,cell);

                        hgc_simhit_energy[recoCell] += simhit.energy();
                    }
                }
                if (fh_simhits_h.isValid())
                {
                    int layer=0,cell=0, sec=0, subsec=0, zp=0,subdet=0;
                    ForwardSubdetector mysubdet;
                    for (const auto& simhit : *fh_simhits_h)
                    {
                        HGCalTestNumbering::unpackHexagonIndex(simhit.id(), subdet, zp, layer, sec, subsec, cell); 
                        mysubdet = (ForwardSubdetector)(subdet);
                        std::pair<int,int> recoLayerCell = topo_fh->dddConstants().simToReco(cell,layer,sec,topo_fh->detectorType());
                        cell  = recoLayerCell.first;
                        layer = recoLayerCell.second;
                        if (layer<0 || cell<0) {
                          continue;
                        }
                        unsigned recoCell = HGCalDetId(mysubdet,zp,layer,subsec,sec,cell);
                        hgc_simhit_energy[recoCell] += simhit.energy();
                    }
                }

                if (bh_simhits_h.isValid() ) 
                {
                       throw cms::Exception("Not Implemented")<<"HGCalTriggerSimCluster: BH simhits not implemnted"; 

                }
                
                //1. construct a cluster container that hosts the cluster per truth-particle
                std::unordered_map<uint64_t,std::pair<int,l1t::HGCalCluster> > cluster_container;// PID-> bx,cluster
                evt.getByToken(sim_token_,sim_handle_);

                if (not sim_handle_.isValid()){
                       throw cms::Exception("ContentError")<<"[HGCalTriggerSimCluster]::[run]::[ERROR] Sim Cluster collection for HGC sim clustering not available"; 
                }
                // calibration
                es.get<IdealGeometryRecord>().get(HGCalEESensitive_, hgceeTopoHandle_);
                es.get<IdealGeometryRecord>().get(HGCalHESiliconSensitive_, hgchefTopoHandle_);

                // 1.5. pre-process the sim cluster to have easy accessible information
                // I want a map cell-> [ (pid, fraction),  ... 
                std::unordered_map<uint32_t, std::vector<std::pair< uint64_t, float > > > simclusters;
                for (auto& cluster : *sim_handle_)
                {
                    auto pid= cluster.particleId(); // not pdgId
                    const auto& hf = cluster.hits_and_fractions();
                    for (const auto & p : hf ) 
                    {
                        simclusters[p.first].push_back( std::pair<uint64_t, float>( pid,p.second) ) ;
                    }
                }
                
                //2. run on the digits,
                for( const auto& digi : coll ) 
                {
                    DATA data;
                    digi.decode(codec_,data);
                    //2.A get the trigger-cell information energy/id
                    //const HGCTriggerDetId& moduleId = digi.getDetId<HGCTriggerDetId>(); // this is a module Det Id

                    // there is a loss of generality here, due to the restriction imposed by the data formats
                    // it will work if inside a module there is a data.payload with an ordered list of all the energies
                    // one may think to add on top of it a wrapper if this stop to be the case for some of the data classes
                    for(const auto& triggercell : data.payload)
                    { 
                            if(triggercell.hwPt()<=0) continue;

                            const HGCalDetId tcellId(triggercell.detId());
                            // calbration
                            int subdet = tcellId.subdetId();
                            int cellThickness = 0;
                            
                            if( subdet == HGCEE ){ 
                                cellThickness = (hgceeTopoHandle_)->dddConstants().waferTypeL((unsigned int)tcellId.wafer() );
                            }else if( subdet == HGCHEF ){
                                cellThickness = (hgchefTopoHandle_)->dddConstants().waferTypeL((unsigned int)tcellId.wafer() );
                            }else if( subdet == HGCHEB ){
                                edm::LogWarning("DataNotFound") << "ATTENTION: the BH trgCells are not yet implemented !! ";
                            }

                            l1t::HGCalTriggerCell calibratedtriggercell(triggercell);
                            calibration_.calibrateInGeV(calibratedtriggercell, cellThickness); 
                            //uint32_t digiEnergy = data.payload; 
                            //auto digiEnergy=triggercell.p4().E();  
                            // using calibrated energy instead
                            auto calibratedDigiEnergy=calibratedtriggercell.p4().E();  
                            double eta=triggercell.p4().Eta();
                            double phi=triggercell.p4().Phi();
                            double z = triggercell.position().z(); // may be useful for cluster shapes
                            //2.B get the HGCAL-base-cell associated to it / geometry
                            
                            // normalization loop 
                            double norm=0.0;
                            map<unsigned, double> energy_for_cluster_shapes;
                            
                            for(const auto& cell : geometry_->getCellsFromTriggerCell( tcellId()) )  // HGCcell -- unsigned
                            {
                                HGCalDetId cellId(cell);

                                //2.C0 find energy of the hgc cell. default is very small value
                                double hgc_energy=1.0e-10; //average if not found / bh
                                const auto &it = hgc_simhit_energy.find(cell);
                                if (it != hgc_simhit_energy.end()) {  hgc_energy = it->second; }

                                //2.C get the particleId and energy fractions
                                const auto & iterator= simclusters.find(cellId);
                                if (iterator == simclusters.end() )  continue;
                                const auto & particles = iterator->second;
                                for ( const auto& p: particles ) 
                                {
                                    const auto & pid= p.first;
                                    const auto & fraction=p.second;
                                    norm += fraction * hgc_energy;
                                    energy_for_cluster_shapes[pid] += calibratedDigiEnergy *fraction *hgc_energy; // norm will be done later, with the above norm
                                }
                            }
                            
                            // second loop counting the energy
                            for(const auto& cell : geometry_->getCellsFromTriggerCell( tcellId()) )  // HGCcell -- unsigned
                            {
                                HGCalDetId cellId(cell);
                                double hgc_energy=1.0e-10; // 1 ->  average if not found / bh
                                const auto &it = hgc_simhit_energy.find(cell);
                                if (it != hgc_simhit_energy.end()) {  hgc_energy = it->second; }
                                //2.C get the particleId and energy fractions
                                const auto & iterator= simclusters.find(cellId);
                                if (iterator == simclusters.end() )  continue;
                                const auto & particles = iterator->second;
                                for ( const auto& p: particles ) 
                                {
                                    const auto & pid= p.first;
                                    const auto & fraction=p.second;
                                    //auto energy = fraction * calibratedDigiEnergy/norm;
                                    auto energy = fraction * hgc_energy* calibratedDigiEnergy/norm; // THIS IS WHAT I WANT
                                    //#warning FIXME_ENERGY
                                    //auto energy = fraction * hgc_energy* hgc_energy/norm;

                                    //2.D add to the corresponding cluster
                                    //eta/phi are the position of the trigger cell, to consider degradation 
                                    addToCluster(cluster_container, pid, 0,energy,eta,phi  ) ;
                                }
                            }

                            // third loop, cluster shapes to ensure the correct counts 
                            for(const auto& iterator :energy_for_cluster_shapes)
                            {
                                double energy = iterator.second / norm;
                                unsigned pid = iterator.first;
                                addToClusterShapes(cluster_container, pid, 0,energy,eta,phi,z  ) ;// only one for trigger cell 
                            }
                    } //end of for-loop
                }

                //3. Push the clusters in the cluster_product_
                for (auto&  p : cluster_container) 
                {
                    //std::unordered_map<uint64_t,std::pair<int,l1t::HGCalCluster> >
                    cluster_product_->push_back(p.second.first,p.second.second); // bx,cluster
                }

            } // end run

template<typename FECODEC , typename DATA >

virtual void HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::setProduces ( edm::stream::EDProducer<> &  prod ) const

inlinefinaloverridevirtual

Implements HGCalTriggerBackendAlgorithmBase.

Definition at line 131 of file HGCalSimCluster.cc.

References HGCalTriggerBackendAlgorithmBase::name(), and parseEventContent::prod.

            {
                prod.produces<l1t::HGCalClusterBxCollection>(name());
            }

Member Data Documentation

template<typename FECODEC , typename DATA >

HGCalTriggerCellCalibration HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::calibration_

private

Definition at line 66 of file HGCalSimCluster.cc.

template<typename FECODEC , typename DATA >

std::unique_ptr<l1t::HGCalClusterBxCollection> HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::cluster_product_

private

Definition at line 55 of file HGCalSimCluster.cc.

template<typename FECODEC , typename DATA >

std::string HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::HGCalEESensitive_

private

Definition at line 64 of file HGCalSimCluster.cc.

template<typename FECODEC , typename DATA >

std::string HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::HGCalHESiliconSensitive_

private

Definition at line 65 of file HGCalSimCluster.cc.

template<typename FECODEC , typename DATA >

edm::ESHandle<HGCalTopology> HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::hgceeTopoHandle_

private

Definition at line 59 of file HGCalSimCluster.cc.

Referenced by HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::run().

template<typename FECODEC , typename DATA >

edm::ESHandle<HGCalTopology> HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::hgchefTopoHandle_

private

Definition at line 60 of file HGCalSimCluster.cc.

Referenced by HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::run().

template<typename FECODEC , typename DATA >

edm::EDGetToken HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::inputbh_

private

Definition at line 71 of file HGCalSimCluster.cc.

template<typename FECODEC , typename DATA >

edm::EDGetToken HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::inputee_

private

Definition at line 71 of file HGCalSimCluster.cc.

template<typename FECODEC , typename DATA >

edm::EDGetToken HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::inputfh_

private

Definition at line 71 of file HGCalSimCluster.cc.

template<typename FECODEC , typename DATA >

edm::Handle< std::vector<SimCluster> > HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::sim_handle_

private

Definition at line 57 of file HGCalSimCluster.cc.

template<typename FECODEC , typename DATA >

edm::EDGetTokenT< std::vector<SimCluster> > HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::sim_token_

private

Definition at line 69 of file HGCalSimCluster.cc.