#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, CaloSubdetectorTopology &topology_p, double deltaEta, double minDeltaPhi, double maxDeltaPhi)
	PreshowerPhiClusterAlgo (float stripEnergyCut)
	PreshowerPhiClusterAlgo ()
	~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

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

Definition at line 24 of file PreshowerPhiClusterAlgo.h.

typedef math::XYZPoint PreshowerPhiClusterAlgo::Point

Definition at line 21 of file PreshowerPhiClusterAlgo.h.

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

Definition at line 23 of file PreshowerPhiClusterAlgo.h.

Constructor & Destructor Documentation

PreshowerPhiClusterAlgo::PreshowerPhiClusterAlgo ( ) [inline]

Definition at line 26 of file PreshowerPhiClusterAlgo.h.

                            : 
    esStripEnergyCut_(0.)
    {}

PreshowerPhiClusterAlgo::PreshowerPhiClusterAlgo ( float stripEnergyCut ) [inline]

Definition at line 30 of file PreshowerPhiClusterAlgo.h.

                                                :
    esStripEnergyCut_(stripEnergyCut)
    {}

PreshowerPhiClusterAlgo::~PreshowerPhiClusterAlgo ( ) [inline]

Definition at line 34 of file PreshowerPhiClusterAlgo.h.

{};

Member Function Documentation

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;
}

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

Definition at line 15 of file PreshowerPhiClusterAlgo.cc.

References edm::SortedCollection< T, SORT >::begin(), reco::deltaPhi(), edm::SortedCollection< T, SORT >::end(), PV3DBase< T, PVType, FrameType >::eta(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), goodStrip(), PV3DBase< T, PVType, FrameType >::phi(), ESDetId::plane(), pos, position, edm::SortedCollection< T, SORT >::push_back(), rechits_map, used_s, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{

  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)

  const CaloCellGeometry *refCell = geometry_p->getGeometry(strip);
  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;

    const CaloCellGeometry *thisCell = geometry_p->getGeometry(strip_it->first);
    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);
  used_strips = used_s;

  return cluster; 
}