|
|
|
Producer for particle flow rechits (PFRecHit) in ECAL. More...
#include <PFRecHitProducerECAL.h>
Public Member Functions | |
| PFRecHitProducerECAL (const edm::ParameterSet &) | |
| ~PFRecHitProducerECAL () | |
Private Member Functions | |
| reco::PFRecHit * | createEcalRecHit (const DetId &detid, double energy, PFLayer::Layer layer, const CaloSubdetectorGeometry *geom) |
| void | createRecHits (std::vector< reco::PFRecHit > &rechits, std::vector< reco::PFRecHit > &rechitsCleaned, edm::Event &, const edm::EventSetup &) |
| void | ecalNeighbArray (const EcalBarrelGeometry &barrelGeom, const CaloSubdetectorTopology &barrelTopo, const EcalEndcapGeometry &endcapGeom, const CaloSubdetectorTopology &endcapTopo) |
| bool | findEcalRecHitGeometry (const DetId &detid, const CaloSubdetectorGeometry *geom, math::XYZVector &position, math::XYZVector &axis) |
| find the position and the axis of the cell for a given rechit | |
| void | findRecHitNeighboursECAL (reco::PFRecHit &rh, const std::map< unsigned, unsigned > &sortedHits) |
| find rechit neighbours, using the hashed index | |
| DetId | move (DetId cell, const CaloDirection &dir) const |
| bool | stdmove (DetId &cell, const CaloDirection &dir, const CaloSubdetectorTopology &barrelTopo, const CaloSubdetectorTopology &endcapTopo, const EcalBarrelGeometry &barrelGeom, const EcalEndcapGeometry &endcapGeom) const |
| bool | stdsimplemove (DetId &cell, const CaloDirection &dir, const CaloSubdetectorTopology &barrelTopo, const CaloSubdetectorTopology &endcapTopo, const EcalBarrelGeometry &barrelGeom, const EcalEndcapGeometry &endcapGeom) const |
Private Attributes | |
| bool | crossBarrelEndcapBorder_ |
| if true, navigation will cross the barrel-endcap border | |
| edm::InputTag | inputTagEcalRecHitsEB_ |
| edm::InputTag | inputTagEcalRecHitsEE_ |
| bool | neighbourmapcalculated_ |
| set to true in ecalNeighbArray | |
| std::vector< std::vector< DetId > > | neighboursEB_ |
| for each ecal barrel rechit, keep track of the neighbours | |
| std::vector< std::vector< DetId > > | neighboursEE_ |
| for each ecal endcap rechit, keep track of the neighbours | |
| double | threshCleaningEB_ |
| double | threshCleaningEE_ |
| bool | timingCleaning_ |
| bool | topologicalCleaning_ |
Producer for particle flow rechits (PFRecHit) in ECAL.
Definition at line 41 of file PFRecHitProducerECAL.h.
| PFRecHitProducerECAL::PFRecHitProducerECAL | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 36 of file PFRecHitProducerECAL.cc.
References crossBarrelEndcapBorder_, edm::ParameterSet::getParameter(), inputTagEcalRecHitsEB_, inputTagEcalRecHitsEE_, neighbourmapcalculated_, threshCleaningEB_, threshCleaningEE_, timingCleaning_, and topologicalCleaning_.
: PFRecHitProducer(iConfig) { // access to the collections of rechits inputTagEcalRecHitsEB_ = iConfig.getParameter<InputTag>("ecalRecHitsEB"); inputTagEcalRecHitsEE_ = iConfig.getParameter<InputTag>("ecalRecHitsEE"); crossBarrelEndcapBorder_ = iConfig.getParameter<bool>("crossBarrelEndcapBorder"); timingCleaning_ = iConfig.getParameter<bool>("timing_Cleaning"); topologicalCleaning_ = iConfig.getParameter<bool>("topological_Cleaning"); threshCleaningEB_ = iConfig.getParameter<double>("thresh_Cleaning_EB"); threshCleaningEE_ = iConfig.getParameter<double>("thresh_Cleaning_EE"); neighbourmapcalculated_ = false; }
| PFRecHitProducerECAL::~PFRecHitProducerECAL | ( | ) |
Definition at line 69 of file PFRecHitProducerECAL.cc.
{}
| reco::PFRecHit * PFRecHitProducerECAL::createEcalRecHit | ( | const DetId & | detid, |
| double | energy, | ||
| PFLayer::Layer | layer, | ||
| const CaloSubdetectorGeometry * | geom | ||
| ) | [private] |
Definition at line 284 of file PFRecHitProducerECAL.cc.
References relval_parameters_module::energy, CaloCellGeometry::getCorners(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), position, DetId::rawId(), reco::PFRecHit::setNECorner(), reco::PFRecHit::setNWCorner(), reco::PFRecHit::setSECorner(), reco::PFRecHit::setSWCorner(), PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::z(), and z.
{
math::XYZVector position;
math::XYZVector axis;
const CaloCellGeometry *thisCell
= geom->getGeometry(detid);
// find rechit geometry
if(!thisCell) {
LogError("PFRecHitProducerECAL")
<<"warning detid "<<detid.rawId()
<<" not found in geometry"<<endl;
return 0;
}
position.SetCoordinates ( thisCell->getPosition().x(),
thisCell->getPosition().y(),
thisCell->getPosition().z() );
// the axis vector is the difference
const TruncatedPyramid* pyr
= dynamic_cast< const TruncatedPyramid* > (thisCell);
if( pyr ) {
axis.SetCoordinates( pyr->getPosition(1).x(),
pyr->getPosition(1).y(),
pyr->getPosition(1).z() );
math::XYZVector axis0( pyr->getPosition(0).x(),
pyr->getPosition(0).y(),
pyr->getPosition(0).z() );
axis -= axis0;
}
else return 0;
// if( !geomfound ) {
// LogError("PFRecHitProducerECAL")<<"cannor find geometry for detid "
// <<detid.rawId()<<" in layer "<<layer<<endl;
// return 0; // geometry not found, skip rechit
// }
reco::PFRecHit *rh
= new reco::PFRecHit( detid.rawId(), layer,
energy,
position.x(), position.y(), position.z(),
axis.x(), axis.y(), axis.z() );
const CaloCellGeometry::CornersVec& corners = thisCell->getCorners();
assert( corners.size() == 8 );
rh->setNECorner( corners[0].x(), corners[0].y(), corners[0].z() );
rh->setSECorner( corners[1].x(), corners[1].y(), corners[1].z() );
rh->setSWCorner( corners[2].x(), corners[2].y(), corners[2].z() );
rh->setNWCorner( corners[3].x(), corners[3].y(), corners[3].z() );
return rh;
}
| void PFRecHitProducerECAL::createRecHits | ( | std::vector< reco::PFRecHit > & | rechits, |
| std::vector< reco::PFRecHit > & | rechitsCleaned, | ||
| edm::Event & | , | ||
| const edm::EventSetup & | |||
| ) | [private, virtual] |
gets ecal barrel and endcap rechits, translate them to PFRecHits, which are stored in the rechits vector
Implements PFRecHitProducer.
| void PFRecHitProducerECAL::ecalNeighbArray | ( | const EcalBarrelGeometry & | barrelGeom, |
| const CaloSubdetectorTopology & | barrelTopo, | ||
| const EcalEndcapGeometry & | endcapGeom, | ||
| const CaloSubdetectorTopology & | endcapTopo | ||
| ) | [private] |
fill the vectors neighboursEB_ and neighboursEE_ which keep track of the neighbours of each rechit. to be called at the beginning of the run
Definition at line 454 of file PFRecHitProducerECAL.cc.
References EAST, DetId::Ecal, EcalBarrel, EcalEndcap, CaloSubdetectorTopology::getWindow(), EEDetId::hashedIndex(), EBDetId::hashedIndex(), recoMuon::in, LogDebug, neighbourmapcalculated_, neighboursEB_, neighboursEE_, NORTH, NORTHEAST, NORTHWEST, findQualityFiles::size, SOUTH, SOUTHEAST, SOUTHWEST, ntuplemaker::status, stdmove(), and WEST.
{
static const CaloDirection orderedDir[8]={SOUTHWEST,
SOUTH,
SOUTHEAST,
WEST,
EAST,
NORTHWEST,
NORTH,
NORTHEAST};
const unsigned nbarrel = 62000;
// Barrel first. The hashed index runs from 0 to 61199
neighboursEB_.resize(nbarrel);
//std::cout << " Building the array of neighbours (barrel) " ;
const std::vector<DetId>& vec(barrelGeom.getValidDetIds(DetId::Ecal,
EcalBarrel));
unsigned size=vec.size();
for(unsigned ic=0; ic<size; ++ic)
{
// We get the 9 cells in a square.
std::vector<DetId> neighbours(barrelTopo.getWindow(vec[ic],3,3));
// std::cout << " Cell " << EBDetId(vec[ic]) << std::endl;
unsigned nneighbours=neighbours.size();
unsigned hashedindex=EBDetId(vec[ic]).hashedIndex();
if(hashedindex>=nbarrel)
{
LogDebug("CaloGeometryTools") << " Array overflow " << std::endl;
}
// If there are 9 cells, it is easy, and this order is know:
// 6 7 8
// 3 4 5
// 0 1 2 (0 = SOUTHWEST)
if(nneighbours==9)
{
neighboursEB_[hashedindex].reserve(8);
for(unsigned in=0;in<nneighbours;++in)
{
// remove the centre
if(neighbours[in]!=vec[ic])
{
neighboursEB_[hashedindex].push_back(neighbours[in]);
// std::cout << " Neighbour " << in << " " << EBDetId(neighbours[in]) << std::endl;
}
}
}
else
{
DetId central(vec[ic]);
neighboursEB_[hashedindex].resize(8,DetId(0));
for(unsigned idir=0;idir<8;++idir)
{
DetId testid=central;
bool status=stdmove(testid,orderedDir[idir],
barrelTopo, endcapTopo,
barrelGeom, endcapGeom);
if(status) neighboursEB_[hashedindex][idir]=testid;
}
}
}
// Moved to the endcap
// std::cout << " done " << size << std::endl;
// std::cout << " Building the array of neighbours (endcap) " ;
// vec.clear();
const std::vector<DetId>& vecee=endcapGeom.getValidDetIds(DetId::Ecal,EcalEndcap);
size=vecee.size();
// There are some holes in the hashedIndex for the EE. Hence the array is bigger than the number
// of crystals
const unsigned nendcap=19960;
neighboursEE_.resize(nendcap);
for(unsigned ic=0; ic<size; ++ic)
{
// We get the 9 cells in a square.
std::vector<DetId> neighbours(endcapTopo.getWindow(vecee[ic],3,3));
unsigned nneighbours=neighbours.size();
// remove the centre
unsigned hashedindex=EEDetId(vecee[ic]).hashedIndex();
if(hashedindex>=nendcap)
{
LogDebug("CaloGeometryTools") << " Array overflow " << std::endl;
}
if(nneighbours==9)
{
neighboursEE_[hashedindex].reserve(8);
for(unsigned in=0;in<nneighbours;++in)
{
// remove the centre
if(neighbours[in]!=vecee[ic])
{
neighboursEE_[hashedindex].push_back(neighbours[in]);
}
}
}
else
{
DetId central(vecee[ic]);
neighboursEE_[hashedindex].resize(8,DetId(0));
for(unsigned idir=0;idir<8;++idir)
{
DetId testid=central;
bool status=stdmove(testid,orderedDir[idir],
barrelTopo, endcapTopo,
barrelGeom, endcapGeom);
if(status) neighboursEE_[hashedindex][idir]=testid;
}
}
}
// std::cout << " done " << size <<std::endl;
neighbourmapcalculated_ = true;
}
| bool PFRecHitProducerECAL::findEcalRecHitGeometry | ( | const DetId & | detid, |
| const CaloSubdetectorGeometry * | geom, | ||
| math::XYZVector & | position, | ||
| math::XYZVector & | axis | ||
| ) | [private] |
find the position and the axis of the cell for a given rechit
Definition at line 354 of file PFRecHitProducerECAL.cc.
References CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), DetId::rawId(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
const CaloCellGeometry *thisCell
= geom->getGeometry(detid);
// find rechit geometry
if(!thisCell) {
LogError("PFRecHitProducerECAL")
<<"warning detid "<<detid.rawId()
<<" not found in geometry"<<endl;
return false;
}
position.SetCoordinates ( thisCell->getPosition().x(),
thisCell->getPosition().y(),
thisCell->getPosition().z() );
// the axis vector is the difference
const TruncatedPyramid* pyr
= dynamic_cast< const TruncatedPyramid* > (thisCell);
if( pyr ) {
axis.SetCoordinates( pyr->getPosition(1).x(),
pyr->getPosition(1).y(),
pyr->getPosition(1).z() );
math::XYZVector axis0( pyr->getPosition(0).x(),
pyr->getPosition(0).y(),
pyr->getPosition(0).z() );
axis -= axis0;
return true;
}
else return false;
}
| void PFRecHitProducerECAL::findRecHitNeighboursECAL | ( | reco::PFRecHit & | rh, |
| const std::map< unsigned, unsigned > & | sortedHits | ||
| ) | [private] |
find rechit neighbours, using the hashed index
Definition at line 401 of file PFRecHitProducerECAL.cc.
References reco::PFRecHit::add4Neighbour(), reco::PFRecHit::add8Neighbour(), reco::PFRecHit::detId(), east, EAST, i, north, NORTH, NORTHEAST, NORTHWEST, DetId::rawId(), SOUTH, south, SOUTHEAST, SOUTHWEST, west, and WEST.
{
DetId center( rh.detId() );
DetId north = move( center, NORTH );
DetId northeast = move( center, NORTHEAST );
DetId northwest = move( center, NORTHWEST );
DetId south = move( center, SOUTH );
DetId southeast = move( center, SOUTHEAST );
DetId southwest = move( center, SOUTHWEST );
DetId east = move( center, EAST );
DetId west = move( center, WEST );
IDH i = sortedHits.find( north.rawId() );
if(i != sortedHits.end() )
rh.add4Neighbour( i->second );
i = sortedHits.find( northeast.rawId() );
if(i != sortedHits.end() )
rh.add8Neighbour( i->second );
i = sortedHits.find( south.rawId() );
if(i != sortedHits.end() )
rh.add4Neighbour( i->second );
i = sortedHits.find( southwest.rawId() );
if(i != sortedHits.end() )
rh.add8Neighbour( i->second );
i = sortedHits.find( east.rawId() );
if(i != sortedHits.end() )
rh.add4Neighbour( i->second );
i = sortedHits.find( southeast.rawId() );
if(i != sortedHits.end() )
rh.add8Neighbour( i->second );
i = sortedHits.find( west.rawId() );
if(i != sortedHits.end() )
rh.add4Neighbour( i->second );
i = sortedHits.find( northwest.rawId() );
if(i != sortedHits.end() )
rh.add8Neighbour( i->second );
}
| DetId PFRecHitProducerECAL::move | ( | DetId | cell, |
| const CaloDirection & | dir | ||
| ) | const [private] |
Definition at line 759 of file PFRecHitProducerECAL.cc.
References dir, EcalBarrel, EEDetId::hashedIndex(), EBDetId::hashedIndex(), neighbourmapcalculated_, neighboursEB_, neighboursEE_, NONE, query::result, and DetId::subdetId().
{
DetId originalcell = cell;
if(dir==NONE || cell==DetId(0)) return false;
// Conversion CaloDirection and index in the table
// CaloDirection :NONE,SOUTH,SOUTHEAST,SOUTHWEST,EAST,WEST, NORTHEAST,NORTHWEST,NORTH
// Table : SOUTHWEST,SOUTH,SOUTHEAST,WEST,EAST,NORTHWEST,NORTH, NORTHEAST
static const int calodirections[9]={-1,1,2,0,4,3,7,5,6};
assert(neighbourmapcalculated_);
DetId result = (originalcell.subdetId()==EcalBarrel) ?
neighboursEB_[EBDetId(originalcell).hashedIndex()][calodirections[dir]]:
neighboursEE_[EEDetId(originalcell).hashedIndex()][calodirections[dir]];
return result;
}
| bool PFRecHitProducerECAL::stdmove | ( | DetId & | cell, |
| const CaloDirection & | dir, | ||
| const CaloSubdetectorTopology & | barrelTopo, | ||
| const CaloSubdetectorTopology & | endcapTopo, | ||
| const EcalBarrelGeometry & | barrelGeom, | ||
| const EcalEndcapGeometry & | endcapGeom | ||
| ) | const [private] |
Definition at line 666 of file PFRecHitProducerECAL.cc.
References EAST, NORTH, NORTHEAST, NORTHWEST, query::result, SOUTH, SOUTHEAST, SOUTHWEST, stdsimplemove(), and WEST.
Referenced by ecalNeighbArray().
{
bool result;
if(dir==NORTH) {
result = stdsimplemove(cell,NORTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
return result;
}
else if(dir==SOUTH) {
result = stdsimplemove(cell,SOUTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
return result;
}
else if(dir==EAST) {
result = stdsimplemove(cell,EAST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
return result;
}
else if(dir==WEST) {
result = stdsimplemove(cell,WEST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
return result;
}
// One has to try both paths
else if(dir==NORTHEAST)
{
result = stdsimplemove(cell,NORTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
if(result)
return stdsimplemove(cell,EAST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
else
{
result = stdsimplemove(cell,EAST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
if(result)
return stdsimplemove(cell,NORTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
else
return false;
}
}
else if(dir==NORTHWEST)
{
result = stdsimplemove(cell,NORTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
if(result)
return stdsimplemove(cell,WEST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
else
{
result = stdsimplemove(cell,WEST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
if(result)
return stdsimplemove(cell,NORTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
else
return false;
}
}
else if(dir == SOUTHEAST)
{
result = stdsimplemove(cell,SOUTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
if(result)
return stdsimplemove(cell,EAST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
else
{
result = stdsimplemove(cell,EAST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
if(result)
return stdsimplemove(cell,SOUTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
else
return false;
}
}
else if(dir == SOUTHWEST)
{
result = stdsimplemove(cell,SOUTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
if(result)
return stdsimplemove(cell,WEST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
else
{
result = stdsimplemove(cell,SOUTH, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
if(result)
return stdsimplemove(cell,WEST, barrelTopo, endcapTopo, barrelGeom, endcapGeom );
else
return false;
}
}
cell = DetId(0);
return false;
}
| bool PFRecHitProducerECAL::stdsimplemove | ( | DetId & | cell, |
| const CaloDirection & | dir, | ||
| const CaloSubdetectorTopology & | barrelTopo, | ||
| const CaloSubdetectorTopology & | endcapTopo, | ||
| const EcalBarrelGeometry & | barrelGeom, | ||
| const EcalEndcapGeometry & | endcapGeom | ||
| ) | const [private] |
Definition at line 588 of file PFRecHitProducerECAL.cc.
References crossBarrelEndcapBorder_, EcalBarrel, EcalEndcap, CaloSubdetectorTopology::getNeighbours(), EBDetId::ietaAbs(), EEDetId::iPhiOuterRing(), EBDetId::MAX_IETA, and DetId::subdetId().
Referenced by stdmove().
{
std::vector<DetId> neighbours;
// BARREL CASE
if(cell.subdetId()==EcalBarrel) {
EBDetId ebDetId = cell;
neighbours = barrelTopo.getNeighbours(ebDetId,dir);
// first try to move according to the standard navigation
if(neighbours.size()>0 && !neighbours[0].null()) {
cell = neighbours[0];
return true;
}
// failed.
if(crossBarrelEndcapBorder_) {
// are we on the outer ring ?
const int ietaAbs ( ebDetId.ietaAbs() ) ; // abs value of ieta
if( EBDetId::MAX_IETA == ietaAbs ) {
// get ee nbrs for for end of barrel crystals
// yes we are
const EcalBarrelGeometry::OrderedListOfEEDetId&
ol( * barrelGeom.getClosestEndcapCells( ebDetId ) ) ;
// take closest neighbour on the other side, that is in the barrel.
cell = *(ol.begin() );
return true;
}
}
}
// ENDCAP CASE
else if(cell.subdetId()==EcalEndcap) {
EEDetId eeDetId = cell;
neighbours= endcapTopo.getNeighbours(eeDetId,dir);
if(neighbours.size()>0 && !neighbours[0].null()) {
cell = neighbours[0];
return true;
}
// failed.
if(crossBarrelEndcapBorder_) {
// are we on the outer ring ?
const int iphi ( eeDetId.iPhiOuterRing() ) ;
if( iphi!= 0) {
// yes we are
const EcalEndcapGeometry::OrderedListOfEBDetId&
ol( * endcapGeom.getClosestBarrelCells( eeDetId ) ) ;
// take closest neighbour on the other side, that is in the barrel.
cell = *(ol.begin() );
return true;
}
}
}
// everything failed
cell = DetId(0);
return false;
}
bool PFRecHitProducerECAL::crossBarrelEndcapBorder_ [private] |
if true, navigation will cross the barrel-endcap border
Definition at line 117 of file PFRecHitProducerECAL.h.
Referenced by PFRecHitProducerECAL(), and stdsimplemove().
Definition at line 120 of file PFRecHitProducerECAL.h.
Referenced by PFRecHitProducerECAL().
Definition at line 121 of file PFRecHitProducerECAL.h.
Referenced by PFRecHitProducerECAL().
bool PFRecHitProducerECAL::neighbourmapcalculated_ [private] |
set to true in ecalNeighbArray
Definition at line 114 of file PFRecHitProducerECAL.h.
Referenced by ecalNeighbArray(), move(), and PFRecHitProducerECAL().
std::vector<std::vector<DetId> > PFRecHitProducerECAL::neighboursEB_ [private] |
for each ecal barrel rechit, keep track of the neighbours
Definition at line 108 of file PFRecHitProducerECAL.h.
Referenced by ecalNeighbArray(), and move().
std::vector<std::vector<DetId> > PFRecHitProducerECAL::neighboursEE_ [private] |
for each ecal endcap rechit, keep track of the neighbours
Definition at line 111 of file PFRecHitProducerECAL.h.
Referenced by ecalNeighbArray(), and move().
double PFRecHitProducerECAL::threshCleaningEB_ [private] |
Definition at line 125 of file PFRecHitProducerECAL.h.
Referenced by PFRecHitProducerECAL().
double PFRecHitProducerECAL::threshCleaningEE_ [private] |
Definition at line 126 of file PFRecHitProducerECAL.h.
Referenced by PFRecHitProducerECAL().
bool PFRecHitProducerECAL::timingCleaning_ [private] |
Definition at line 124 of file PFRecHitProducerECAL.h.
Referenced by PFRecHitProducerECAL().
bool PFRecHitProducerECAL::topologicalCleaning_ [private] |
Definition at line 129 of file PFRecHitProducerECAL.h.
Referenced by PFRecHitProducerECAL().
1.7.3