#include <ClusterTools.h>

Public Member Functions
	ClusterTools (const edm::ParameterSet &, edm::ConsumesCollector &)

float	getClusterHadronFraction (const reco::CaloCluster &) const

void	getEvent (const edm::Event &)

void	getEventSetup (const edm::EventSetup &)

double	getMultiClusterEnergy (const reco::HGCalMultiCluster &) const

math::XYZPoint	getMultiClusterPosition (const reco::HGCalMultiCluster &, double vz=0.) const

	~ClusterTools ()

Private Attributes
const HGCRecHitCollection *	bhrh_

const edm::EDGetTokenT < HGCRecHitCollection >	bhtok

const HGCRecHitCollection *	eerh_

const edm::EDGetTokenT < HGCRecHitCollection >	eetok

const HGCRecHitCollection *	fhrh_

const edm::EDGetTokenT < HGCRecHitCollection >	fhtok

RecHitTools	rhtools_

Detailed Description

Definition at line 26 of file ClusterTools.h.

Constructor & Destructor Documentation

ClusterTools::ClusterTools	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	sumes
	)

Definition at line 18 of file ClusterTools.cc.

                                                        :
   eetok( sumes.consumes<HGCRecHitCollection>(conf.getParameter<edm::InputTag>("HGCEEInput")) ),
   fhtok( sumes.consumes<HGCRecHitCollection>(conf.getParameter<edm::InputTag>("HGCFHInput")) ),
   bhtok( sumes.consumes<HGCRecHitCollection>(conf.getParameter<edm::InputTag>("HGCBHInput")) ) {
 }

hgcal::ClusterTools::~ClusterTools ( )

inline

Definition at line 29 of file ClusterTools.h.

29 {}

Member Function Documentation

float ClusterTools::getClusterHadronFraction ( const reco::CaloCluster & clus ) const

Definition at line 40 of file ClusterTools.cc.

References bhrh_, eerh_, relval_parameters_module::energy, Exception, f, fhrh_, edm::SortedCollection< T, SORT >::find(), DetId::Forward, HLT_FULL_cff::fraction, DetId::Hcal, HcalEndcap, HGCEE, HGCHEF, reco::CaloCluster::hitsAndFractions(), and groupFilesInBlocks::temp.

                                                                               {
   float energy=0.f, energyHad=0.f;
   const auto& hits = clus.hitsAndFractions();
   for( const auto& hit : hits ) {
     const auto& id = hit.first;
     const float fraction = hit.second;
     if( id.det() == DetId::Forward ) {
       switch( id.subdetId() ) {
       case HGCEE:
         energy += eerh_->find(id)->energy()*fraction;
         break;
       case HGCHEF:
         {
           const float temp = fhrh_->find(id)->energy();
           energy += temp*fraction;
           energyHad += temp*fraction;
         }
         break;
       default:
         throw cms::Exception("HGCalClusterTools")
           << " Cluster contains hits that are not from HGCal! " << std::endl;
       }
     } else if ( id.det() == DetId::Hcal && id.subdetId() == HcalEndcap ) {
       const float temp = bhrh_->find(id)->energy();
       energy += temp*fraction;
       energyHad += temp*fraction;
     } else {
       throw cms::Exception("HGCalClusterTools")
         << " Cluster contains hits that are not from HGCal! " << std::endl;
     }    
   }
   float fraction = -1.f;
   if( energy > 0.f ) {
     fraction = energyHad/energy;
   }
   return fraction;
 }

void ClusterTools::getEvent ( const edm::Event & ev )

Definition at line 25 of file ClusterTools.cc.

References bhrh_, bhtok, eerh_, eetok, fhrh_, fhtok, edm::Event::getByToken(), hgcal::RecHitTools::getEvent(), edm::Handle< T >::product(), rhtools_, and groupFilesInBlocks::temp.

                                               {
   rhtools_.getEvent(ev);
   edm::Handle<HGCRecHitCollection> temp;
   ev.getByToken(eetok, temp);
   eerh_ = temp.product();
   ev.getByToken(fhtok, temp);
   fhrh_ = temp.product();
   ev.getByToken(bhtok, temp);
   bhrh_ = temp.product();
 }

void ClusterTools::getEventSetup ( const edm::EventSetup & es )

Definition at line 36 of file ClusterTools.cc.

References hgcal::RecHitTools::getEventSetup(), and rhtools_.

                                                         {
   rhtools_.getEventSetup(es);
 }

double ClusterTools::getMultiClusterEnergy ( const reco::HGCalMultiCluster & clu ) const

Definition at line 105 of file ClusterTools.cc.

References reco::HGCalMultiCluster::clusters().

                                                                                  {
   double acc = 0.0;
   for(const auto& ptr : clu.clusters() ) {
     acc += ptr->energy();
   }
   return acc;
 }

math::XYZPoint ClusterTools::getMultiClusterPosition	(	const reco::HGCalMultiCluster &	clu,
		double	vz = `0.`
	)		const

Definition at line 79 of file ClusterTools.cc.

References assert(), reco::HGCalMultiCluster::clusters(), f, mathSSE::sqrt(), groupFilesInBlocks::temp, histoStyle::weight, and x().

                                                                                                     {
   if( clu.clusters().size() == 0 ) return math::XYZPoint();
   double acc_rho = 0.0;
   double acc_eta = 0.0;
   double acc_phi = 0.0;
   double totweight = 0.;
   for( const auto& ptr : clu.clusters() ) {   
     const double x = ptr->x();    
     const double y = ptr->y();
     const double point_r2 = (x*x + y*y);    
     const double point_z = ptr->z()-vz;
     const double point_h = std::sqrt(point_r2 + point_z*point_z);
     const double weight = ptr->energy() * ptr->size();
     assert((y != 0. || x != 0.) && "Cluster position somehow in beampipe.");
     assert(point_z != 0.f && "Layer-cluster position given as reference point.");
     const double point_r = std::sqrt(point_r2);
     acc_rho += point_r * weight;
     acc_phi += vdt::fast_atan2(y,x) * weight;
     acc_eta += -1. * vdt::fast_log(point_r/(point_z + point_h)) * weight;
     totweight += weight;
   }
   const double invweight = 1.0/totweight;
   reco::PFCluster::REPPoint temp(acc_rho*invweight,acc_eta*invweight,acc_phi*invweight);
   return math::XYZPoint(temp.x(),temp.y(),temp.z());
 }