PreshowerPhiClusterAlgo Class Reference

#include <PreshowerPhiClusterAlgo.h>

Public Types
typedef std::set< DetId >	HitsID
typedef math::XYZPoint	Point
typedef std::map< DetId, EcalRecHit >	RecHitsMap

Public Member Functions
bool	goodStrip (RecHitsMap::iterator candidate_it)
reco::PreshowerCluster	makeOneCluster (ESDetId strip, HitsID *used_strips, RecHitsMap *rechits_map, const CaloSubdetectorGeometry *geometry_p, double deltaEta, double minDeltaPhi, double maxDeltaPhi)
	PreshowerPhiClusterAlgo ()
	PreshowerPhiClusterAlgo (float stripEnergyCut)
	~PreshowerPhiClusterAlgo ()

Private Attributes
float	esStripEnergyCut_
RecHitsMap *	rechits_map
std::vector< ESDetId >	road_2d
HitsID *	used_s

Detailed Description

Definition at line 17 of file PreshowerPhiClusterAlgo.h.

Member Typedef Documentation

HitsID

typedef std::set<DetId> PreshowerPhiClusterAlgo::HitsID

Definition at line 22 of file PreshowerPhiClusterAlgo.h.

Point

typedef math::XYZPoint PreshowerPhiClusterAlgo::Point

Definition at line 19 of file PreshowerPhiClusterAlgo.h.

RecHitsMap

typedef std::map<DetId, EcalRecHit> PreshowerPhiClusterAlgo::RecHitsMap

Definition at line 21 of file PreshowerPhiClusterAlgo.h.

Constructor & Destructor Documentation

PreshowerPhiClusterAlgo() [1/2]

PreshowerPhiClusterAlgo::PreshowerPhiClusterAlgo ( )

inline

Definition at line 24 of file PreshowerPhiClusterAlgo.h.

: esStripEnergyCut_(0.) {}

PreshowerPhiClusterAlgo() [2/2]

PreshowerPhiClusterAlgo::PreshowerPhiClusterAlgo ( float  stripEnergyCut )

inline

Definition at line 26 of file PreshowerPhiClusterAlgo.h.

: esStripEnergyCut_(stripEnergyCut) {}

~PreshowerPhiClusterAlgo()

PreshowerPhiClusterAlgo::~PreshowerPhiClusterAlgo ( )

inline

Definition at line 28 of file PreshowerPhiClusterAlgo.h.

{};

Member Function Documentation

goodStrip()

bool PreshowerPhiClusterAlgo::goodStrip

(

RecHitsMap::iterator

candidate_it

)

Definition at line 99 of file PreshowerPhiClusterAlgo.cc.

References esStripEnergyCut_, LogTrace, rechits_map, and used_s.

Referenced by makeOneCluster().

                                                                       {
  if (used_s->find(candidate_it->first) != used_s->end())
    LogTrace("PreShowerPhiClusterAlgo") << " This strip is in use";
  if (candidate_it == rechits_map->end())
    LogTrace("PreShowerPhiClusterAlgo") << " No such a strip in rechits_map";
  if (candidate_it->second.energy() <= esStripEnergyCut_)
    LogTrace("PreShowerPhiClusterAlgo") << "Strip energy" << candidate_it->second.energy() << "is below threshold";

  // crystal should not be included...
  if ((used_s->find(candidate_it->first) != used_s->end()) ||  // ...if it already belongs to a cluster
      (candidate_it == rechits_map->end()) ||                  // ...if it corresponds to a hit
      (candidate_it->second.energy() <= esStripEnergyCut_))    // ...if it has a negative or zero energy
  {
    return false;
  }
  return true;
}

makeOneCluster()

reco::PreshowerCluster PreshowerPhiClusterAlgo::makeOneCluster	(	ESDetId	strip,
		HitsID *	used_strips,
		RecHitsMap *	rechits_map,
		const CaloSubdetectorGeometry *	geometry_p,
		double	deltaEta,
		double	minDeltaPhi,
		double	maxDeltaPhi
	)

Definition at line 15 of file PreshowerPhiClusterAlgo.cc.

References edm::SortedCollection< T, SORT >::begin(), spr::deltaEta, reco::deltaPhi(), edm::SortedCollection< T, SORT >::end(), PV3DBase< T, PVType, FrameType >::eta(), CaloSubdetectorGeometry::getGeometry(), goodStrip(), CommPDSkim_cfg::maxDeltaPhi, VBFGenJetFilter_cfi::minDeltaPhi, PV3DBase< T, PVType, FrameType >::phi(), ESDetId::plane(), position, edm::SortedCollection< T, SORT >::push_back(), rechits_map, digitizers_cfi::strip, and used_s.

                                                                                   {
  rechits_map = the_rechitsMap_p;
  used_s = used_strips;

  int plane = strip.plane();

  // create null-cluster
  std::vector<std::pair<DetId, float> > dummy;
  Point posi(0, 0, 0);
  reco::PreshowerCluster nullcluster = reco::PreshowerCluster(0., posi, dummy, plane);
  if (strip == ESDetId(0))
    return nullcluster;  // works in case of no intersected strip found (e.g. in the Barrel)

  auto refCell = geometry_p->getGeometry(strip);
  const GlobalPoint& refpos = refCell->getPosition();
  double refEta = refpos.eta();
  double refPhi = refpos.phi();

  // Collection of cluster strips
  EcalRecHitCollection clusterRecHits;

  double x_pos = 0;
  double y_pos = 0;
  double z_pos = 0;

  RecHitsMap::iterator strip_it;
  for (strip_it = rechits_map->begin(); strip_it != rechits_map->end(); ++strip_it) {
    if (!goodStrip(strip_it))
      continue;

    ESDetId mystrip = (strip_it->first == DetId(0)) ? ESDetId(0) : ESDetId(strip_it->first);
    if (mystrip.plane() != plane)
      continue;

    auto thisCell = geometry_p->getGeometry(strip_it->first);
    const GlobalPoint& position = thisCell->getPosition();

    if (fabs(position.eta() - refEta) < deltaEta) {
      //std::cout<<"all strips : "<<mystrip.plane()<<" "<<position.phi()<<" "<<reco::deltaPhi(position.phi(), refPhi)<<std::endl;

      if (reco::deltaPhi(position.phi(), refPhi) > 0 && reco::deltaPhi(position.phi(), refPhi) > maxDeltaPhi)
        continue;
      if (reco::deltaPhi(position.phi(), refPhi) < 0 && reco::deltaPhi(position.phi(), refPhi) < minDeltaPhi)
        continue;

      //std::cout<<mystrip.zside()<<" "<<mystrip.plane()<<" "<<mystrip.six()<<" "<<mystrip.siy()<<" "<<mystrip.strip()<<" "<<position.eta()<<" "<<position.phi()<<" "<<strip_it->second.energy()<<" "<<strip_it->second.recoFlag()<<" "<<refEta<<" "<<refPhi<<" "<<reco::deltaPhi(position.phi(), refPhi)<<std::endl;

      clusterRecHits.push_back(strip_it->second);
      used_s->insert(strip_it->first);

      x_pos += strip_it->second.energy() * position.x();
      y_pos += strip_it->second.energy() * position.y();
      z_pos += strip_it->second.energy() * position.z();
    }
  }

  EcalRecHitCollection::iterator it;
  double Eclust = 0;

  std::vector<std::pair<DetId, float> > usedHits;
  for (it = clusterRecHits.begin(); it != clusterRecHits.end(); it++) {
    Eclust += it->energy();
    usedHits.push_back(std::pair<DetId, float>(it->id(), 1.));
  }

  if (Eclust > 0.) {
    x_pos /= Eclust;
    y_pos /= Eclust;
    z_pos /= Eclust;
  }
  Point pos(x_pos, y_pos, z_pos);

  reco::PreshowerCluster cluster = reco::PreshowerCluster(Eclust, pos, usedHits, plane);

  return cluster;
}

Member Data Documentation

esStripEnergyCut_

float PreshowerPhiClusterAlgo::esStripEnergyCut_

private

Definition at line 41 of file PreshowerPhiClusterAlgo.h.

Referenced by goodStrip().

rechits_map

RecHitsMap* PreshowerPhiClusterAlgo::rechits_map

private

Definition at line 46 of file PreshowerPhiClusterAlgo.h.

Referenced by goodStrip(), and makeOneCluster().

road_2d

std::vector<ESDetId> PreshowerPhiClusterAlgo::road_2d

private

Definition at line 43 of file PreshowerPhiClusterAlgo.h.

used_s

HitsID* PreshowerPhiClusterAlgo::used_s

private

Definition at line 49 of file PreshowerPhiClusterAlgo.h.

Referenced by goodStrip(), and makeOneCluster().