|
|
|
#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 |
Definition at line 22 of file HiBremRecoveryClusterAlgo.h.
| 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;
}
| 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;
}
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().
Definition at line 83 of file HiBremRecoveryClusterAlgo.h.
Referenced by makeIslandSuperClusters(), and makeSuperClusters().
Definition at line 70 of file HiBremRecoveryClusterAlgo.h.
Referenced by HiBremRecoveryClusterAlgo(), makeIslandSuperClusters(), makeSuperClusters(), match(), and setVerbosity().
1.7.3