|
|
|
#include <CosmicClusterAlgo.h>
Definition at line 28 of file CosmicClusterAlgo.h.
point in the space
Definition at line 63 of file CosmicClusterAlgo.h.
Definition at line 32 of file CosmicClusterAlgo.h.
| CosmicClusterAlgo::CosmicClusterAlgo | ( | ) | [inline] |
Definition at line 35 of file CosmicClusterAlgo.h.
{
}
| CosmicClusterAlgo::CosmicClusterAlgo | ( | double | ebst, |
| double | ebSt, | ||
| double | ebDt, | ||
| double | ebSp, | ||
| double | ecst, | ||
| double | ecSt, | ||
| double | ecDt, | ||
| double | ecSp, | ||
| const PositionCalc & | posCalc, | ||
| VerbosityLevel | the_verbosity = pERROR |
||
| ) | [inline] |
Definition at line 38 of file CosmicClusterAlgo.h.
References posCalculator_.
:
ecalBarrelSeedThreshold(ebst), ecalBarrelSingleThreshold(ebSt), ecalBarrelSecondThreshold(ebDt), ecalBarrelSupThreshold(ebSp), ecalEndcapSeedThreshold(ecst), ecalEndcapSingleThreshold(ecSt), ecalEndcapSecondThreshold(ecDt), ecalEndcapSupThreshold(ecSp), verbosity(the_verbosity) {
posCalculator_ = posCalc;
}
| virtual CosmicClusterAlgo::~CosmicClusterAlgo | ( | ) | [inline, virtual] |
Definition at line 43 of file CosmicClusterAlgo.h.
{
}
| void CosmicClusterAlgo::addCrystal | ( | const DetId & | det, |
| const bool | in9 | ||
| ) | [private] |
Definition at line 406 of file CosmicClusterAlgo.cc.
References current_v25, current_v9, edm::SortedCollection< T, SORT >::end(), edm::SortedCollection< T, SORT >::find(), recHits_, uncalibRecHits_, and used_s.
Referenced by prepareCluster().
{
EcalRecHitCollection::const_iterator thisIt = recHits_->find(det);
if ((thisIt != recHits_->end()) && (thisIt->id() != DetId(0)))
{
if ((used_s.find(thisIt->id()) == used_s.end()))
{
EcalUncalibratedRecHitCollection::const_iterator thisItu = uncalibRecHits_->find(det);
used_s.insert(det);
if ((thisIt->energy() >= -1.) && !(thisItu->chi2() < -1.))
{
if (in9) current_v9.push_back(det);
current_v25.push_back(det);
}
}
}
}
| bool CosmicClusterAlgo::checkMaxima | ( | CaloNavigator< DetId > & | navigator | ) | [private] |
Definition at line 307 of file CosmicClusterAlgo.cc.
References CaloNavigator< T >::east(), edm::SortedCollection< T, SORT >::end(), edm::SortedCollection< T, SORT >::find(), CaloNavigator< T >::home(), i, CaloNavigator< T >::north(), CaloNavigator< T >::pos(), recHits_, CaloNavigator< T >::south(), and CaloNavigator< T >::west().
Referenced by mainSearch().
{
bool maxima = true;
EcalRecHitCollection::const_iterator thisHit;
EcalRecHitCollection::const_iterator seedHit = recHits_->find(navigator.pos());
double thisEnergy = 0.;
double seedEnergy = seedHit->energy();
std::vector<DetId> swissCrossVec;
swissCrossVec.clear();
swissCrossVec.push_back(navigator.west());
navigator.home();
swissCrossVec.push_back(navigator.east());
navigator.home();
swissCrossVec.push_back(navigator.north());
navigator.home();
swissCrossVec.push_back(navigator.south());
navigator.home();
std::vector<DetId>::const_iterator detItr;
for (unsigned int i = 0; i < swissCrossVec.size(); ++i)
{
thisHit = recHits_->find(swissCrossVec[i]);
if ((swissCrossVec[i] == DetId(0)) || thisHit == recHits_->end()) thisEnergy = 0.0;
else thisEnergy = thisHit->energy();
if (thisEnergy > seedEnergy)
{
maxima = false;
break;
}
}
return maxima;
}
| void CosmicClusterAlgo::mainSearch | ( | const CaloSubdetectorGeometry * | geometry_p, |
| const CaloSubdetectorTopology * | topology_p, | ||
| const CaloSubdetectorGeometry * | geometryES_p, | ||
| EcalPart | ecalPart | ||
| ) | [private] |
Definition at line 154 of file CosmicClusterAlgo.cc.
References canSeed_s, checkMaxima(), gather_cfg::cout, current_v25, current_v25Sup, current_v9, makeCluster(), pINFO, prepareCluster(), seeds, used_s, and verbosity.
Referenced by makeClusters().
{
if (verbosity < pINFO)
{
std::cout << "Building clusters............" << std::endl;
}
// Loop over seeds:
std::vector<EcalRecHit>::iterator it;
for (it = seeds.begin(); it != seeds.end(); it++)
{
// check if this crystal is able to seed
// (event though it is already used)
bool usedButCanSeed = false;
if (canSeed_s.find(it->id()) != canSeed_s.end()) usedButCanSeed = true;
// make sure the current seed does not belong to a cluster already.
if ((used_s.find(it->id()) != used_s.end()) && (usedButCanSeed == false))
{
if (it == seeds.begin())
{
if (verbosity < pINFO)
{
std::cout << "##############################################################" << std::endl;
std::cout << "DEBUG ALERT: Highest energy seed already belongs to a cluster!" << std::endl;
std::cout << "##############################################################" << std::endl;
}
}
// seed crystal is used or is used and cannot seed a cluster
// so continue to the next seed crystal...
continue;
}
// clear the vector of hits in current clusterconst EcalRecHitCollection* hits,
current_v9.clear();
current_v25.clear();
current_v25Sup.clear();
// Create a navigator at the seed
CaloNavigator<DetId> navigator(it->id(), topology_p);
DetId seedId = navigator.pos();
navigator.setHome(seedId);
// Is the seed a local maximum?
bool localMaxima = checkMaxima(navigator);
if (localMaxima)
{
// build the 5x5 taking care over which crystals //JHaupt 4-27-08 3x3 is a good option...
// can seed new clusters and which can't
prepareCluster(navigator, geometry_p);
}
// If some crystals in the current vector then
// make them into a cluster
if (current_v25.size() > 0)
{
makeCluster(geometry_p, geometryES_p, seedId);
}
} // End loop on seed crystals
}
| void CosmicClusterAlgo::makeCluster | ( | const CaloSubdetectorGeometry * | geometry_p, |
| const CaloSubdetectorGeometry * | geometryES_p, | ||
| DetId | seedId | ||
| ) | [private] |
Definition at line 224 of file CosmicClusterAlgo.cc.
References EcalUncalibratedRecHit::amplitude(), PositionCalc::Calculate_Location(), clusters_v, gather_cfg::cout, current_v25, current_v25Sup, current_v9, ecalBarrelSecondThreshold, ecalBarrelSingleThreshold, ecalBarrelSupThreshold, ecalEndcapSecondThreshold, ecalEndcapSingleThreshold, ecalEndcapSupThreshold, edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), relval_parameters_module::energy, edm::SortedCollection< T, SORT >::find(), EcalRecHit::id(), EcalUncalibratedRecHit::id(), inEB, EcalRecHit::kGood, EcalRecHit::kOutOfTime, EcalRecHit::kPoorCalib, reco::CaloCluster::multi5x5, pINFO, posCalculator_, position, recHits_, EcalRecHit::recoFlag(), swap(), uncalibRecHits_, and verbosity.
Referenced by mainSearch().
{
double energy = 0;
double energySecond = 0.;//JHaupt 4-27-08 Added for the second crystal stream
double energyMax = 0.;//JHaupt 4-27-08 Added for the max crystal stream
double chi2 = 0;
DetId detFir;
DetId detSec;
//bool goodCluster = false; //JHaupt 4-27-08 Added so that some can be earased.. used another day Might not be needed as seeds are energy ordered...
std::vector<DetId>::iterator it;
for (it = current_v9.begin(); it != current_v9.end(); it++)
{
// Martina - check recoFlag for crystals sorrounding the good one
EcalRecHitCollection::const_iterator itt = recHits_->find(*it);
// double-check that iterator can be dereferenced
if(itt==recHits_->end()) continue;
EcalRecHit hit_p = *itt;
// if rechit affected by features other than these, do not allow 2 crystal seeding either
// features are hard coded, for now; add severity?
uint32_t rhFlag = (*itt).recoFlag();
if (!(
rhFlag == EcalRecHit::kGood ||
rhFlag == EcalRecHit::kOutOfTime ||
rhFlag == EcalRecHit::kPoorCalib
)
) continue;
EcalUncalibratedRecHitCollection::const_iterator ittu = uncalibRecHits_->find(*it);
EcalUncalibratedRecHit uhit_p = *ittu;
if (uhit_p.amplitude() > energySecond ) {energySecond = uhit_p.amplitude(); detSec = uhit_p.id();}
if (energySecond > energyMax ) {std::swap(energySecond,energyMax); std::swap(detFir,detSec);}
}
//if ((detFir.det() == DetId::Ecal) && (detFir.subdetId() == EcalEndcap)) inEB = false;
//We ignore the possiblity that the cluster may span into the EE
if ((energyMax < (inEB ? ecalBarrelSingleThreshold : ecalEndcapSingleThreshold)) && (energySecond < (inEB ? ecalBarrelSecondThreshold : ecalEndcapSecondThreshold) )) return;
for (it = current_v25.begin(); it != current_v25.end(); it++)
{
EcalRecHitCollection::const_iterator itt = recHits_->find(*it);
EcalRecHit hit_p = *itt;
EcalUncalibratedRecHitCollection::const_iterator ittu = uncalibRecHits_->find(*it);
EcalUncalibratedRecHit uhit_p = *ittu;
if ( uhit_p.amplitude() > ( inEB ? ecalBarrelSupThreshold : ecalEndcapSupThreshold))
{
// energy fraction = 1
current_v25Sup.push_back( std::pair<DetId, float>( hit_p.id(), 1.) );
energy += hit_p.energy(); //Keep the fully corrected energy
chi2 += 0;
}
}
Point position;
position = posCalculator_.Calculate_Location(current_v25Sup,recHits_, geometry, geometryES);
chi2 /= energy;
if (verbosity < pINFO)
{
std::cout << "JH******** NEW CLUSTER ********" << std::endl;
std::cout << "JHNo. of crystals = " << current_v25Sup.size() << std::endl;
std::cout << "JH Energy = " << energy << std::endl;
std::cout << "JH Phi = " << position.phi() << std::endl;
std::cout << "JH Eta = " << position.eta() << std::endl;
std::cout << "JH*****************************" << std::endl;
// specialize this
// std::cout << "JH****Emax**** "<<energyMax << " ieta " <<detFir.ieta() <<" iphi "<<detFir.ieta() << std::endl;
// std::cout << "JH****Esec**** "<<energySecond << " ieta " <<detSec.ieta() <<" iphi "<<detSec.ieta() << std::endl;
}
clusters_v.push_back(reco::BasicCluster(energy, position, reco::CaloID(), current_v25Sup, reco::CaloCluster::multi5x5, seedId));
}
| std::vector< reco::BasicCluster > CosmicClusterAlgo::makeClusters | ( | const EcalRecHitCollection * | hits, |
| const EcalUncalibratedRecHitCollection * | uncalibhits, | ||
| const CaloSubdetectorGeometry * | geometry, | ||
| const CaloSubdetectorTopology * | topology_p, | ||
| const CaloSubdetectorGeometry * | geometryES_p, | ||
| EcalPart | ecalPart, | ||
| bool | regional = false, |
||
| const std::vector< EcalEtaPhiRegion > & | regions = std::vector<EcalEtaPhiRegion>() |
||
| ) |
Definition at line 18 of file CosmicClusterAlgo.cc.
References EcalUncalibratedRecHit::amplitude(), barrel, edm::SortedCollection< T, SORT >::begin(), canSeed_s, clusters_v, gather_cfg::cout, ecalBarrelSeedThreshold, ecalEndcapSeedThreshold, edm::SortedCollection< T, SORT >::end(), endcap, edm::SortedCollection< T, SORT >::find(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), inEB, EcalRecHit::kGood, EcalRecHit::kOutOfTime, EcalRecHit::kPoorCalib, mainSearch(), pINFO, position, recHits_, seeds, python::multivaluedict::sort(), dtT0WireCalibration_cfg::threshold, uncalibRecHits_, used_s, and verbosity.
Referenced by CosmicClusterProducer::clusterizeECALPart().
{
seeds.clear();
used_s.clear();
canSeed_s.clear();
clusters_v.clear();
recHits_ = hits;
uncalibRecHits_ = uncalibhits;
inEB = true;
double threshold = 0;
std::string ecalPart_string;
if (ecalPart == endcap)
{
threshold = ecalEndcapSeedThreshold;
ecalPart_string = "EndCap";
inEB = false;
}
if (ecalPart == barrel)
{
threshold = ecalBarrelSeedThreshold;
ecalPart_string = "Barrel";
inEB = true;
}
if (verbosity < pINFO)
{
std::cout << "-------------------------------------------------------------" << std::endl;
std::cout << "Island algorithm invoked for ECAL" << ecalPart_string << std::endl;
std::cout << "Looking for seeds, threshold used = " << threshold << " ADC" <<std::endl;
}
int nregions=0;
if(regional) nregions=regions.size();
if(!regional || nregions) {
EcalRecHitCollection::const_iterator it;
for(it = recHits_->begin(); it != recHits_->end(); it++)
{
if (!uncalibRecHits_){
if (verbosity < pINFO)
{
std::cout << "-------------------------------------------------------------" << std::endl;
std::cout << "No Uncalibrated RecHits no Uncalibrated rec hit collection available" << std::endl;
}
continue;
}
// if rechit affected by features other than these, do not allow if seeding
// features are hard coded, for now; add severity?
uint32_t rhFlag = (*it).recoFlag();
if (!(
rhFlag == EcalRecHit::kGood ||
rhFlag == EcalRecHit::kOutOfTime ||
rhFlag == EcalRecHit::kPoorCalib
)
) continue;
EcalUncalibratedRecHitCollection::const_iterator itt = uncalibRecHits_->find(it->id());
if (itt == uncalibRecHits_->end()){
if (verbosity < pINFO)
{
std::cout << "-------------------------------------------------------------" << std::endl;
std::cout << "No Uncalibrated RecHit associated with the RecHit Probably no Uncalibrated rec hit collection available" << std::endl;
}
continue;
}
EcalUncalibratedRecHit uhit_p = *itt;
// looking for cluster seeds
if (uhit_p.amplitude() < (inEB ? ecalBarrelSeedThreshold : ecalEndcapSeedThreshold) ) continue; //
const CaloCellGeometry *thisCell = geometry_p->getGeometry(it->id());
GlobalPoint position = thisCell->getPosition();
// Require that RecHit is within clustering region in case
// of regional reconstruction
bool withinRegion = false;
if (regional) {
std::vector<EcalEtaPhiRegion>::const_iterator region;
for (region=regions.begin(); region!=regions.end(); region++) {
if (region->inRegion(position)) {
withinRegion = true;
break;
}
}
}
if (!regional || withinRegion) {
//float ET = it->energy() * sin(position.theta()); JHaupt Out 4-27-08 Et not needed for Cosmic Events...
// if (energy >= threshold)
seeds.push_back(*it); // JHaupt 4-27-2008 Et -> energy, most likely not needed as there is already a threshold requirement.
}
}
}
sort(seeds.begin(), seeds.end(), EcalRecHitLess());
if (verbosity < pINFO)
{
std::cout << "JH Total number of seeds found in event = " << seeds.size() << std::endl;
for (EcalRecHitCollection::const_iterator ji = seeds.begin(); ji != seeds.end(); ++ji)
{
//std::cout << "JH Seed Energy " << ji->energy() << " hashed " << ((EBDetId)ji->id()).hashedIndex() << std::endl;
}
}
mainSearch(geometry_p,topology_p,geometryES_p,ecalPart );
sort(clusters_v.rbegin(), clusters_v.rend(),ClusterEtLess());
if (verbosity < pINFO)
{
std::cout << "---------- end of main search. clusters have been sorted ----" << std::endl;
}
return clusters_v;
}
| void CosmicClusterAlgo::prepareCluster | ( | CaloNavigator< DetId > & | navigator, |
| const CaloSubdetectorGeometry * | geometry | ||
| ) | [private] |
Definition at line 345 of file CosmicClusterAlgo.cc.
References abs, addCrystal(), canSeed_s, CaloNavigator< T >::home(), and CaloNavigator< T >::offsetBy().
Referenced by mainSearch().
{
DetId thisDet;
std::set<DetId>::iterator setItr;
// now add the 5x5 taking care to mark the edges
// as able to seed and where overlapping in the central
// region with crystals that were previously able to seed
// change their status so they are not able to seed
//std::cout << std::endl;
for (int dx = -2; dx < 3; ++dx) //for (int dx = -2; dx < 3; ++dx)
{
for (int dy = -2; dy < 3; ++ dy) //for (int dy = -2; dy < 3; ++ dy)
{
// navigate in free steps forming
// a full 5x5
thisDet = navigator.offsetBy(dx, dy);
navigator.home();
// add the current crystal
//std::cout << "adding " << dx << ", " << dy << std::endl;
// now consider if we are in an edge (outer 16)
// or central (inner 9) region
if ((abs(dx) > 1) || (abs(dy) > 1))
{
// this is an "edge" so should be allowed to seed
// provided it is not already used
//std::cout << " setting can seed" << std::endl;
addCrystal(thisDet,false); //These are in the V25
canSeed_s.insert(thisDet);
} // end if "edge"
else
{
// or else we are in the central 3x3
// and must remove any of these crystals from the canSeed set
setItr = canSeed_s.find(thisDet);
addCrystal(thisDet,true); //These are in the V9
if (setItr != canSeed_s.end())
{
//std::cout << " unsetting can seed" << std::endl;
canSeed_s.erase(setItr);
}
} // end if "centre"
} // end loop on dy
} // end loop on dx
//std::cout << "*** " << std::endl;
//std::cout << " current_v contains " << current_v.size() << std::endl;
//std::cout << "*** " << std::endl;
}
| void CosmicClusterAlgo::setVerbosity | ( | VerbosityLevel | the_verbosity | ) | [inline] |
Definition at line 47 of file CosmicClusterAlgo.h.
References verbosity.
{
verbosity = the_verbosity;
}
std::set<DetId> CosmicClusterAlgo::canSeed_s [private] |
Definition at line 92 of file CosmicClusterAlgo.h.
Referenced by mainSearch(), makeClusters(), and prepareCluster().
std::vector<reco::BasicCluster> CosmicClusterAlgo::clusters_v [private] |
Definition at line 104 of file CosmicClusterAlgo.h.
Referenced by makeCluster(), and makeClusters().
std::vector<DetId> CosmicClusterAlgo::current_v25 [private] |
Definition at line 100 of file CosmicClusterAlgo.h.
Referenced by addCrystal(), mainSearch(), and makeCluster().
std::vector< std::pair<DetId, float> > CosmicClusterAlgo::current_v25Sup [private] |
Definition at line 101 of file CosmicClusterAlgo.h.
Referenced by mainSearch(), and makeCluster().
std::vector<DetId> CosmicClusterAlgo::current_v9 [private] |
Definition at line 99 of file CosmicClusterAlgo.h.
Referenced by addCrystal(), mainSearch(), and makeCluster().
double CosmicClusterAlgo::ecalBarrelSecondThreshold [private] |
Definition at line 74 of file CosmicClusterAlgo.h.
Referenced by makeCluster().
double CosmicClusterAlgo::ecalBarrelSeedThreshold [private] |
Definition at line 72 of file CosmicClusterAlgo.h.
Referenced by makeClusters().
double CosmicClusterAlgo::ecalBarrelSingleThreshold [private] |
Definition at line 73 of file CosmicClusterAlgo.h.
Referenced by makeCluster().
double CosmicClusterAlgo::ecalBarrelSupThreshold [private] |
Definition at line 75 of file CosmicClusterAlgo.h.
Referenced by makeCluster().
double CosmicClusterAlgo::ecalEndcapSecondThreshold [private] |
Definition at line 79 of file CosmicClusterAlgo.h.
Referenced by makeCluster().
double CosmicClusterAlgo::ecalEndcapSeedThreshold [private] |
Definition at line 77 of file CosmicClusterAlgo.h.
Referenced by makeClusters().
double CosmicClusterAlgo::ecalEndcapSingleThreshold [private] |
Definition at line 78 of file CosmicClusterAlgo.h.
Referenced by makeCluster().
double CosmicClusterAlgo::ecalEndcapSupThreshold [private] |
Definition at line 80 of file CosmicClusterAlgo.h.
Referenced by makeCluster().
bool CosmicClusterAlgo::inEB [private] |
Definition at line 96 of file CosmicClusterAlgo.h.
Referenced by makeCluster(), and makeClusters().
Definition at line 68 of file CosmicClusterAlgo.h.
Referenced by CosmicClusterAlgo(), and makeCluster().
const EcalRecHitCollection* CosmicClusterAlgo::recHits_ [private] |
Definition at line 83 of file CosmicClusterAlgo.h.
Referenced by addCrystal(), checkMaxima(), makeCluster(), and makeClusters().
std::vector<EcalRecHit> CosmicClusterAlgo::seeds [private] |
Definition at line 88 of file CosmicClusterAlgo.h.
Referenced by mainSearch(), and makeClusters().
const EcalUncalibratedRecHitCollection* CosmicClusterAlgo::uncalibRecHits_ [private] |
Definition at line 85 of file CosmicClusterAlgo.h.
Referenced by addCrystal(), makeCluster(), and makeClusters().
std::set<DetId> CosmicClusterAlgo::used_s [private] |
Definition at line 91 of file CosmicClusterAlgo.h.
Referenced by addCrystal(), mainSearch(), and makeClusters().
VerbosityLevel CosmicClusterAlgo::verbosity [private] |
Definition at line 107 of file CosmicClusterAlgo.h.
Referenced by mainSearch(), makeCluster(), makeClusters(), and setVerbosity().
1.7.3