#include <Multi5x5BremRecoveryClusterAlgo.h>

Public Member Functions
reco::SuperClusterCollection	makeSuperClusters (const reco::CaloClusterPtrVector &clusters)

	Multi5x5BremRecoveryClusterAlgo (const edm::ParameterSet &bremRecoveryPset, double eb_sc_road_etasize=0.06, double eb_sc_road_phisize=0.80, double ec_sc_road_etasize=0.14, double ec_sc_road_phisize=0.40, bool dynamicPhiRoad=true, double theSeedTransverseEnergyThreshold=0.40)

	~Multi5x5BremRecoveryClusterAlgo ()

Private Member Functions
void	makeIslandSuperClusters (reco::CaloClusterPtrVector &clusters_v, double etaRoad, double phiRoad)

Private Attributes
bool	dynamicPhiRoad_

double	eb_rdeta_

double	eb_rdphi_

double	ec_rdeta_

double	ec_rdphi_

BremRecoveryPhiRoadAlgo *	phiRoadAlgo_

double	seedTransverseEnergyThreshold

reco::SuperClusterCollection	superclusters_v

Detailed Description

Definition at line 24 of file Multi5x5BremRecoveryClusterAlgo.h.

Constructor & Destructor Documentation

Multi5x5BremRecoveryClusterAlgo::Multi5x5BremRecoveryClusterAlgo	(	const edm::ParameterSet &	bremRecoveryPset,
		double	eb_sc_road_etasize = `0.06`,
		double	eb_sc_road_phisize = `0.80`,
		double	ec_sc_road_etasize = `0.14`,
		double	ec_sc_road_phisize = `0.40`,
		bool	dynamicPhiRoad = `true`,
		double	theSeedTransverseEnergyThreshold = `0.40`
	)

inline

Definition at line 29 of file Multi5x5BremRecoveryClusterAlgo.h.

References cosmicSuperClusters_cfi::dynamicPhiRoad, dynamicPhiRoad_, eb_rdeta_, eb_rdphi_, ec_rdeta_, ec_rdphi_, phiRoadAlgo_, and seedTransverseEnergyThreshold.

     {
       // e*_rdeta_ and e*_rdphi_ are half the total window 
       // because they correspond to one direction (positive or negative)
       eb_rdeta_ = eb_sc_road_etasize / 2;
       eb_rdphi_ = eb_sc_road_phisize / 2;
       ec_rdeta_ = ec_sc_road_etasize / 2;
       ec_rdphi_ = ec_sc_road_phisize / 2;
 
       seedTransverseEnergyThreshold = theSeedTransverseEnergyThreshold;
       dynamicPhiRoad_ = dynamicPhiRoad;
       if (dynamicPhiRoad_) phiRoadAlgo_ = new BremRecoveryPhiRoadAlgo(bremRecoveryPset);
 
     }

Multi5x5BremRecoveryClusterAlgo::~Multi5x5BremRecoveryClusterAlgo ( )

inline

Definition at line 52 of file Multi5x5BremRecoveryClusterAlgo.h.

References fastPrimaryVertexProducer_cfi::clusters, dynamicPhiRoad_, makeIslandSuperClusters(), makeSuperClusters(), and phiRoadAlgo_.

   {
      if (dynamicPhiRoad_) delete phiRoadAlgo_;
   } 

Member Function Documentation

void Multi5x5BremRecoveryClusterAlgo::makeIslandSuperClusters	(	reco::CaloClusterPtrVector &	clusters_v,
		double	etaRoad,
		double	phiRoad
	)

private

Definition at line 41 of file Multi5x5BremRecoveryClusterAlgo.cc.

References funct::abs(), edm::PtrVectorBase::clear(), reco::dp, dynamicPhiRoad_, edm::PtrVectorBase::empty(), BremRecoveryPhiRoadAlgo::endcapPhiRoad(), reco::CaloCluster::energy(), stringResolutionProvider_cfi::et, PVValHelper::eta, objects.autophobj::float, mps_fire::i, LogTrace, M_PI, match(), p1, p2, phi, phiRoadAlgo_, reco::CaloCluster::position(), position, edm::PtrVector< T >::push_back(), seedTransverseEnergyThreshold, funct::sin(), edm::PtrVectorBase::size(), superclusters_v, theta(), mitigatedMETSequence_cff::U, X, DOFs::Y, and DOFs::Z.

Referenced by makeSuperClusters(), and ~Multi5x5BremRecoveryClusterAlgo().

 {
   if(clusters_v.empty()) return;
 
   const auto clustersSize = clusters_v.size();
   assert(clustersSize>0);
 
   bool usedSeed[clustersSize];
   for (auto ic=0U; ic<clustersSize; ++ic) usedSeed[ic]=false;
  
   float eta[clustersSize], phi[clustersSize], et[clustersSize];
   for (auto ic=0U; ic<clustersSize; ++ic) {
     eta[ic]=clusters_v[ic]->eta();
     phi[ic]=clusters_v[ic]->phi();
     et[ic]=clusters_v[ic]->energy() * sin(clusters_v[ic]->position().theta());
   }
 
   for (auto is=0U; is<clustersSize; ++is) {
     // check this seed was not already used
     if (usedSeed[is]) continue;
     auto const & currentSeed = clusters_v[is];
     
     // Does our highest energy cluster have high enough energy?
     // changed this to continue from break (to be robust against the order of sorting of the seed clusters)
     if (et[is] < seedTransverseEnergyThreshold) continue;
 
     // if yes, make it a seed for a new SuperCluster:
     double energy = (currentSeed)->energy();
     math::XYZVector position_((currentSeed)->position().X(), 
                   (currentSeed)->position().Y(), 
                   (currentSeed)->position().Z());
     position_ *= energy;
     usedSeed[is]=true;
 
     LogTrace("EcalClusters") << "*****************************";
     LogTrace("EcalClusters") << "******NEW SUPERCLUSTER*******";
     LogTrace("EcalClusters") << "Seed R = " << (currentSeed)->position().Rho();
     
     reco::CaloClusterPtrVector constituentClusters;
     constituentClusters.push_back(currentSeed);
     auto ic = is + 1;
         
     while (ic < clustersSize) {
       auto const & currentCluster = clusters_v[ic];
   
       // if dynamic phi road is enabled then compute the phi road for a cluster
       // of energy of existing clusters + the candidate cluster.
       if (dynamicPhiRoad_)
     phiRoad = phiRoadAlgo_->endcapPhiRoad(energy + (currentCluster)->energy());
 
       auto dphi = [](float p1,float p2) {
     auto dp=std::abs(p1-p2); if (dp>float(M_PI)) dp-=float(2*M_PI);
     return std::abs(dp);
       };
 
       auto match = [&](int i, int j) {
     return (dphi(phi[i],phi[j])< phiRoad) & (std::abs(eta[i]-eta[j])< etaRoad);
       };
       
       // does the cluster match the phi road for this candidate supercluster
       if (!usedSeed[ic] && match(is,ic)) {
     
     // add basic cluster to supercluster constituents
     constituentClusters.push_back(currentCluster);
     energy   += (currentCluster)->energy();
     position_ += (currentCluster)->energy() * math::XYZVector((currentCluster)->position().X(),
                                   (currentCluster)->position().Y(), 
                                   (currentCluster)->position().Z()
                                   );
     // remove cluster from vector of available clusters
     usedSeed[ic]=true;
     LogTrace("EcalClusters") << "Cluster R = " << (currentCluster)->position().Rho();
       }
       ++ic;
     }
     
     position_ /= energy;
 
     LogTrace("EcalClusters") <<"Final SuperCluster R = " << position_.Rho();
     
     reco::SuperCluster newSuperCluster(energy, 
                        math::XYZPoint(position_.X(), position_.Y(), position_.Z()),
                        currentSeed, 
                        constituentClusters);
     
     superclusters_v.push_back(newSuperCluster);
     LogTrace("EcalClusters") << "created a new supercluster of: ";
     LogTrace("EcalClusters") << "Energy = " << newSuperCluster.energy();
     LogTrace("EcalClusters") << "Position in (R, phi, theta) = ("
                  << newSuperCluster.position().Rho() << ", "
                  << newSuperCluster.position().phi() << ", "
                  << newSuperCluster.position().theta() << ")" ;
     
     
   }
   
   clusters_v.clear();
   
 }

reco::SuperClusterCollection Multi5x5BremRecoveryClusterAlgo::makeSuperClusters ( const reco::CaloClusterPtrVector & clusters )

Definition at line 5 of file Multi5x5BremRecoveryClusterAlgo.cc.

References eb_rdeta_, eb_rdphi_, ec_rdeta_, ec_rdphi_, makeIslandSuperClusters(), reco::CaloCluster::multi5x5, edm::PtrVector< T >::push_back(), and superclusters_v.

Referenced by Multi5x5SuperClusterProducer::produceSuperclustersForECALPart(), and ~Multi5x5BremRecoveryClusterAlgo().

 {
 
     const float etaBorder = 1.479;
     superclusters_v.clear();
 
     // create vectors of references to clusters of a specific origin...
     reco::CaloClusterPtrVector islandClustersBarrel_v;
     reco::CaloClusterPtrVector islandClustersEndCap_v;
 
     // ...and populate them:
     for (auto const& cluster_p : clustersCollection)
     {
         if (cluster_p->algo() == reco::CaloCluster::multi5x5) 
         {
             if (fabs(cluster_p->position().eta()) < etaBorder)
             {
                 islandClustersBarrel_v.push_back(cluster_p);
             }
             else
             {
                 islandClustersEndCap_v.push_back(cluster_p);
             }
         }
     }
 
     // make the superclusters from the Barrel clusters - Island
     makeIslandSuperClusters(islandClustersBarrel_v, eb_rdeta_, eb_rdphi_);
     // make the superclusters from the EndCap clusters - Island
     makeIslandSuperClusters(islandClustersEndCap_v, ec_rdeta_, ec_rdphi_);
 
     return superclusters_v;
 }