HiBremRecoveryClusterAlgo Class Reference

#include <HiBremRecoveryClusterAlgo.h>

Public Types
enum	VerbosityLevel { pDEBUG = 0, pWARNING = 1, pINFO = 2, pERROR = 3 }

Public Member Functions
	HiBremRecoveryClusterAlgo (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, double theSeedTransverseEnergyThreshold=0.40, double theBarrelBremEnergyThreshold=2.3, double theEndcapBremEnergyThreshold=5.7, VerbosityLevel the_verbosity=pERROR)
reco::SuperClusterCollection	makeSuperClusters (reco::CaloClusterPtrVector &clusters)
void	setVerbosity (VerbosityLevel the_verbosity)

Private Member Functions
void	makeIslandSuperClusters (reco::CaloClusterPtrVector &clusters_v, double etaRoad, double phiRoad)
bool	match (reco::CaloClusterPtr seed_p, reco::CaloClusterPtr cluster_p, double etaRoad, double phiRoad)

Private Attributes
double	BarrelBremEnergyThreshold
double	eb_rdeta_
double	eb_rdphi_
double	ec_rdeta_
double	ec_rdphi_
double	EndcapBremEnergyThreshold
double	seedTransverseEnergyThreshold
reco::SuperClusterCollection	superclusters_v
VerbosityLevel	verbosity

Detailed Description

Definition at line 22 of file HiBremRecoveryClusterAlgo.h.

Member Enumeration Documentation

enum HiBremRecoveryClusterAlgo::VerbosityLevel

Enumerator
pDEBUG
pWARNING
pINFO
pERROR

Definition at line 26 of file HiBremRecoveryClusterAlgo.h.

{ pDEBUG = 0, pWARNING = 1, pINFO = 2, pERROR = 3 }; 

Constructor & Destructor Documentation

HiBremRecoveryClusterAlgo::HiBremRecoveryClusterAlgo ( 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, double  theSeedTransverseEnergyThreshold = 0.40, double  theBarrelBremEnergyThreshold = 2.3, double  theEndcapBremEnergyThreshold = 5.7, VerbosityLevel  the_verbosity = pERROR  )

inline

Definition at line 28 of file HiBremRecoveryClusterAlgo.h.

References BarrelBremEnergyThreshold, eb_rdeta_, eb_rdphi_, ec_rdeta_, ec_rdphi_, EndcapBremEnergyThreshold, seedTransverseEnergyThreshold, and verbosity.

  {
      // 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;
      BarrelBremEnergyThreshold = theBarrelBremEnergyThreshold;
      EndcapBremEnergyThreshold = theEndcapBremEnergyThreshold;
      verbosity = the_verbosity;
  }

Member Function Documentation

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

private

Definition at line 47 of file HiBremRecoveryClusterAlgo.cc.

References edm::PtrVector< T >::begin(), edm::PtrVectorBase::clear(), gather_cfg::cout, edm::PtrVector< T >::end(), spr::find(), match(), pINFO, edm::PtrVector< T >::push_back(), seedTransverseEnergyThreshold, funct::sin(), superclusters_v, and verbosity.

Referenced by makeSuperClusters().

{
  // Vector of usedSeedEnergy, use the seed energy to check if this cluster is used.
  std::vector<double> usedSeedEnergy;
  usedSeedEnergy.clear();

  // Main brem recovery loop
  for ( reco::CaloCluster_iterator currentSeed = clusters_v.begin(); currentSeed != clusters_v.end(); ++currentSeed)
  {
     if (verbosity <= pINFO) {
        std::cout <<"Current Cluster: "<<(*currentSeed)->energy()<<" "<<(std::find(usedSeedEnergy.begin(), usedSeedEnergy.end(), (*currentSeed)->energy()) 
             != usedSeedEnergy.end())<<std::endl;
     }

     // check this seed was not already used
     if (std::find(usedSeedEnergy.begin(), usedSeedEnergy.end(), (*currentSeed)->energy()) 
             != usedSeedEnergy.end()) continue;

     // Does our highest energy cluster have high enough energy? If not, continue instead of break to be robust
     if ((*currentSeed)->energy() * sin((*currentSeed)->position().theta()) < seedTransverseEnergyThreshold) continue;

     // if yes, make it a seed for a new SuperCluster, the position of the SC is calculated by energy weighted sum:
     double energy_ = (*currentSeed)->energy();
     math::XYZVector position_((*currentSeed)->position().X(), 
                       (*currentSeed)->position().Y(), 
                       (*currentSeed)->position().Z());
     position_ *= energy_;
     usedSeedEnergy.push_back((*currentSeed)->energy());

     // Printout if verbose
     if (verbosity <= pINFO)
     {
       std::cout << "*****************************" << std::endl;
       std::cout << "******NEW SUPERCLUSTER*******" << std::endl;
       std::cout << "Seed R = " << (*currentSeed)->position().Rho() << std::endl;
       std::cout << "Seed Et = " << (*currentSeed)->energy()* sin((*currentSeed)->position().theta()) << std::endl;
     }

     // and add the matching clusters:
     reco::CaloClusterPtrVector constituentClusters;
     constituentClusters.push_back( *currentSeed );
     reco::CaloCluster_iterator currentCluster = currentSeed + 1;

     // loop over the clusters
     while (currentCluster != clusters_v.end())
     {
        // Print out the basic clusters
        if (verbosity <= pINFO) {
           std::cout <<"->Cluster: "<<(*currentCluster)->energy()
                     <<" Used = "<<(std::find(usedSeedEnergy.begin(), usedSeedEnergy.end(), (*currentCluster)->energy()) != usedSeedEnergy.end())
                     <<" Matched = "<<match(*currentSeed, *currentCluster, etaRoad, phiRoad)<<std::endl;
        }

        // if it's in the search window, and unused
        if (match(*currentSeed, *currentCluster, etaRoad, phiRoad)
            && (std::find(usedSeedEnergy.begin(), usedSeedEnergy.end(), (*currentCluster)->energy()) == usedSeedEnergy.end()))
        {
           // Add basic cluster
           constituentClusters.push_back(*currentCluster);
           energy_   += (*currentCluster)->energy();
           position_ += (*currentCluster)->energy() * math::XYZVector((*currentCluster)->position().X(), 
                                                                      (*currentCluster)->position().Y(), 
                                                                      (*currentCluster)->position().Z()); 
           // Add the cluster to the used list
           usedSeedEnergy.push_back((*currentCluster)->energy());

           if (verbosity <= pINFO) 
           {
              std::cout << "Cluster R = " << (*currentCluster)->position().Rho() << std::endl;
           }

        }
        ++currentCluster;
     }

     // Calculate the final position
     position_ /= energy_;

     if (verbosity <= pINFO)
     {
       std::cout << "Final SuperCluster R = " << position_.Rho() << std::endl;
     }

     // Add the supercluster to the new collection
     reco::SuperCluster newSuperCluster(energy_, 
                        math::XYZPoint(position_.X(), position_.Y(), position_.Z()), 
                        (*currentSeed), 
                        constituentClusters);

     superclusters_v.push_back(newSuperCluster);

     if (verbosity <= pINFO)
       {
         std::cout << "created a new supercluster of: " << std::endl;
         std::cout << "Energy = " << newSuperCluster.energy() << std::endl;
         std::cout << "Position in (R, phi, theta, eta) = (" 
               << newSuperCluster.position().Rho() << ", " 
               << newSuperCluster.position().phi() << ", "
               << newSuperCluster.position().theta() << ", "
               << newSuperCluster.position().eta() << ")" << std::endl;
       }
  }
  clusters_v.clear();
  usedSeedEnergy.clear();
}

reco::SuperClusterCollection HiBremRecoveryClusterAlgo::makeSuperClusters

(

reco::CaloClusterPtrVector &

clusters

)

Definition at line 5 of file HiBremRecoveryClusterAlgo.cc.

References BarrelBremEnergyThreshold, edm::PtrVector< T >::begin(), gather_cfg::cout, eb_rdeta_, eb_rdphi_, ec_rdeta_, ec_rdphi_, edm::PtrVector< T >::end(), EndcapBremEnergyThreshold, reco::CaloCluster::island, makeIslandSuperClusters(), pINFO, edm::PtrVector< T >::push_back(), superclusters_v, and verbosity.

Referenced by HiSuperClusterProducer::produceSuperclustersForECALPart().

{
  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 (reco::CaloCluster_iterator it = clustersCollection.begin(); it != clustersCollection.end(); it++)
  {
    
     reco::CaloClusterPtr cluster_p = *it;
     if (cluster_p->algo() == reco::CaloCluster::island) 
     {
        if (verbosity <= pINFO)
        {
           std::cout <<"Basic Cluster: (eta,phi,energy) = "<<cluster_p->position().eta()<<" "<<cluster_p->position().phi()<<" "
                                      <<cluster_p->energy()<<std::endl;
        }

        // if the basic cluster pass the energy threshold -> fill it into the list
        if (fabs(cluster_p->position().eta()) < etaBorder)
        {
           if (cluster_p->energy() > BarrelBremEnergyThreshold) islandClustersBarrel_v.push_back(cluster_p);
        } else {
           if (cluster_p->energy() > EndcapBremEnergyThreshold) 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;
}

bool HiBremRecoveryClusterAlgo::match ( reco::CaloClusterPtr  seed_p, reco::CaloClusterPtr  cluster_p, double  etaRoad, double  phiRoad  )

private

Definition at line 155 of file HiBremRecoveryClusterAlgo.cc.

References funct::cos(), gather_cfg::cout, dPhi(), pINFO, and verbosity.

Referenced by makeIslandSuperClusters().

{
  math::XYZPoint clusterPosition = cluster_p->position();
  math::XYZPoint seedPosition = seed_p->position();

  double dPhi = acos(cos(seedPosition.phi() - clusterPosition.phi()));
 
  double dEta = fabs(seedPosition.eta() - clusterPosition.eta());
  if (verbosity <= pINFO) {
     std::cout <<"seed phi: "<<seedPosition.phi()<<" cluster phi: "<<clusterPosition.phi()<<" dphi = "<<dPhi<<" dphiMax = "<<dPhiMax<<std::endl;
     std::cout <<"seed eta: "<<seedPosition.eta()<<" cluster eta: "<<clusterPosition.eta()<<" deta = "<<dEta<<" detaMax = "<<dEtaMax<<std::endl;
  }
  if (dEta > dEtaMax) return false;
  if (dPhi > dPhiMax) return false;

  return true;
}

void HiBremRecoveryClusterAlgo::setVerbosity ( VerbosityLevel  the_verbosity )

inline

Definition at line 51 of file HiBremRecoveryClusterAlgo.h.

References verbosity.

  {
      verbosity = the_verbosity;
  }

Member Data Documentation

double HiBremRecoveryClusterAlgo::BarrelBremEnergyThreshold

private

Definition at line 80 of file HiBremRecoveryClusterAlgo.h.

Referenced by HiBremRecoveryClusterAlgo(), and makeSuperClusters().

double HiBremRecoveryClusterAlgo::eb_rdeta_

private

Definition at line 72 of file HiBremRecoveryClusterAlgo.h.

Referenced by HiBremRecoveryClusterAlgo(), and makeSuperClusters().

double HiBremRecoveryClusterAlgo::eb_rdphi_

private

Definition at line 73 of file HiBremRecoveryClusterAlgo.h.

Referenced by HiBremRecoveryClusterAlgo(), and makeSuperClusters().

double HiBremRecoveryClusterAlgo::ec_rdeta_

private

Definition at line 74 of file HiBremRecoveryClusterAlgo.h.

Referenced by HiBremRecoveryClusterAlgo(), and makeSuperClusters().

double HiBremRecoveryClusterAlgo::ec_rdphi_

private

Definition at line 75 of file HiBremRecoveryClusterAlgo.h.

Referenced by HiBremRecoveryClusterAlgo(), and makeSuperClusters().

double HiBremRecoveryClusterAlgo::EndcapBremEnergyThreshold

private

Definition at line 81 of file HiBremRecoveryClusterAlgo.h.

Referenced by HiBremRecoveryClusterAlgo(), and makeSuperClusters().

double HiBremRecoveryClusterAlgo::seedTransverseEnergyThreshold

private

Definition at line 77 of file HiBremRecoveryClusterAlgo.h.

Referenced by HiBremRecoveryClusterAlgo(), and makeIslandSuperClusters().

reco::SuperClusterCollection HiBremRecoveryClusterAlgo::superclusters_v

private

Definition at line 83 of file HiBremRecoveryClusterAlgo.h.

Referenced by makeIslandSuperClusters(), and makeSuperClusters().

VerbosityLevel HiBremRecoveryClusterAlgo::verbosity

private

Definition at line 70 of file HiBremRecoveryClusterAlgo.h.

Referenced by HiBremRecoveryClusterAlgo(), makeIslandSuperClusters(), makeSuperClusters(), match(), and setVerbosity().