ggPFESClusters Class Reference

#include <ggPFESClusters.h>

Public Member Functions
std::map< unsigned int, unsigned int >	getClosestEECls (vector< reco::PreshowerCluster > PFPS, vector< reco::CaloCluster > PF)
double	getLinkDist (reco::PreshowerCluster clusterPS, reco::CaloCluster clusterECAL)
virtual vector < reco::PreshowerCluster >	getPFESClusters (reco::SuperCluster sc)
	ggPFESClusters (edm::Handle< EcalRecHitCollection > &ESRecHits, const CaloSubdetectorGeometry *geomEnd)
virtual	~ggPFESClusters ()

Private Attributes
Handle< EcalRecHitCollection >	ESRecHits_
const CaloSubdetectorGeometry *	geomEnd_

Detailed Description

Definition at line 15 of file ggPFESClusters.h.

Constructor & Destructor Documentation

ggPFESClusters::ggPFESClusters ( edm::Handle< EcalRecHitCollection > &  ESRecHits, const CaloSubdetectorGeometry *  geomEnd  )

explicit

Definition at line 11 of file ggPFESClusters.cc.

                    :
  ESRecHits_(ESRecHits),
  geomEnd_(geomEnd)
{
  
}

ggPFESClusters::~ggPFESClusters ( )

virtual

Definition at line 21 of file ggPFESClusters.cc.

                               {

}

Member Function Documentation

std::map< unsigned int, unsigned int > ggPFESClusters::getClosestEECls	(	vector< reco::PreshowerCluster >	PFPS,
		vector< reco::CaloCluster >	PF
	)

Definition at line 79 of file ggPFESClusters.cc.

References getLinkDist(), i, and j.

                                                                                                                              {
  std::map<unsigned int,unsigned int> linkedIndices;
  std::map<double, unsigned int> links;

  for(unsigned int j=0; j<PFPS.size(); ++j){
    for(unsigned int i=0; i<PF.size(); ++i){
      double link_ = getLinkDist(PFPS[j], PF[i]);
      if(link_!=-1.)
    links.insert(std::make_pair(link_,i));
    }
    std::map<double,unsigned int>::iterator itlinks = links.begin();
    if((*itlinks).first>=0.)
      linkedIndices.insert(std::make_pair(j, (*itlinks).second));
    links.clear();
  }

  return linkedIndices;
}

double ggPFESClusters::getLinkDist	(	reco::PreshowerCluster	clusterPS,
		reco::CaloCluster	clusterECAL
	)

Definition at line 100 of file ggPFESClusters.cc.

References first, geomEnd_, CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), reco::CaloCluster::hitsAndFractions(), AlCaHLTBitMon_ParallelJobs::p, reco::PreshowerCluster::plane(), reco::CaloCluster::position(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), vdt::x, detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by getClosestEECls().

                                                                                             {

  static double resPSpitch = 0.19/sqrt(12.);
  static double resPSlength = 6.1/sqrt(12.);

  // ECAL cluster position
  double zECAL  = clusterECAL.position().Z(); // cluster centroid position 
  double xECAL  = clusterECAL.position().X();
  double yECAL  = clusterECAL.position().Y();

  // PS cluster position 
  double zPS = clusterPS.position().Z();
  double xPS = clusterPS.position().X();
  double yPS = clusterPS.position().Y();

  // Check that zEcal and zPs have the same sign
  if (zECAL*zPS < 0.) return -1.;
  double deltaX = 0.;
  double deltaY = 0.;
  double sqr12 = std::sqrt(12.);
  switch (clusterPS.plane()) {
  case 1:
    // vertical strips, measure x with pitch precision
    deltaX = resPSpitch * sqr12;
    deltaY = resPSlength * sqr12;
    break;
  case 2:
    // horizontal strips, measure y with pitch precision
    deltaY = resPSpitch * sqr12;
    deltaX = resPSlength * sqr12;
    break;
  default:
    break;
  }

  // Get the rechits
  const std::vector< std::pair<DetId, float> > fracs = clusterECAL.hitsAndFractions();
  bool linkedbyrechit = false;
  //loop rechits
  for(unsigned int rhit = 0; rhit < fracs.size(); rhit++){
    double fraction = fracs[rhit].second;
    if(fraction < 1E-4) continue;
    if(fracs[rhit].first.null()) continue;

    //getting rechit center position
    DetId eehitid=DetId(fracs[rhit].first.rawId());
    const CaloCellGeometry *xtal = geomEnd_->getGeometry(eehitid);
    if(!xtal) continue;
    
    //getting rechit corners
    const CaloCellGeometry::CornersVec& corners_vec = xtal->getCorners();
    math::XYZPoint corners[4];
    // They have this order just for compatibility with the original function. 
    corners[0] = math::XYZPoint(corners_vec[3].x(), corners_vec[3].y(), corners_vec[3].z() );
    corners[1] = math::XYZPoint(corners_vec[2].x(), corners_vec[2].y(), corners_vec[2].z() );
    corners[2] = math::XYZPoint(corners_vec[1].x(), corners_vec[1].y(), corners_vec[1].z() );
    corners[3] = math::XYZPoint(corners_vec[0].x(), corners_vec[0].y(), corners_vec[0].z() );
    const GlobalPoint p = geomEnd_->getGeometry(eehitid)->getPosition();

    // Extrapolating from the xtal's surface (because of working with DetId's) 
    const math::XYZPoint& posxyz = math::XYZPoint(p.x()*zPS/p.z(), p.y()*zPS/p.z(), p.z()*zPS/p.z());

    double x[5];
    double y[5];
    for ( unsigned jc=0; jc<4; ++jc ) {
      // corner position projected onto the preshower
      math::XYZPoint cornerpos = corners[jc] * zPS/p.z();
      // Inflate the size by the size of the PS strips, and by 5% to include ECAL cracks.
      x[jc] = cornerpos.X() + (cornerpos.X()-posxyz.X()) * (0.05 +1.0/fabs((cornerpos.X()-posxyz.X()))*deltaX/2.);
      y[jc] = cornerpos.Y() + (cornerpos.Y()-posxyz.Y()) * (0.05 +1.0/fabs((cornerpos.Y()-posxyz.Y()))*deltaY/2.);
    }//loop corners

    //need to close the polygon in order to use the TMath::IsInside fonction from root lib
    x[4] = x[0];
    y[4] = y[0];

    //Check if the extrapolation point of the track falls within the rechit boundaries
    bool isinside = TMath::IsInside(xPS,yPS,5,x,y);
      
    if( isinside ){
      linkedbyrechit = true;
      break;
    }

  }//loop rechits
  
  if( linkedbyrechit ) {
    double dist = std::sqrt( (xECAL/1000. - xPS/1000.)*(xECAL/1000. - xPS/1000.) 
              + (yECAL/1000. - yPS/1000.)*(yECAL/1000. - yPS/1000.) );
    return dist;
  } else { 
    return -1.;
  }
  
}

vector< reco::PreshowerCluster > ggPFESClusters::getPFESClusters ( reco::SuperCluster  sc )

virtual

Definition at line 26 of file ggPFESClusters.cc.

References gather_cfg::cout, ESRecHits_, first, findQualityFiles::jj, GetRecoTauVFromDQM_MC_cff::kk, ESDetId::plane(), reco::SuperCluster::preshowerClustersBegin(), reco::SuperCluster::preshowerClustersEnd(), and alignCSCRings::s.

Referenced by ggPFPhotons::fillPFClusters().

                                   {
  std::vector<PreshowerCluster>PFPreShowerClust;
  int jj = 0;
  float energies[100];
  for(reco::CaloCluster_iterator ps=sc.preshowerClustersBegin(); ps!=sc.preshowerClustersEnd(); ++ps){
    if(jj>99)
      cout << endl << "Attention! More than 100 ES clusters!" << endl;
    else {
      energies[jj] = (*ps)->energy();
      bool doubl = false;
      for(int kk=0; kk<jj; kk++){
        if(fabs(energies[kk]-(*ps)->energy())<0.0000001){
          //cout << endl << "Duplicate cluster";
          doubl = true;
      break;
        }
      }
      jj++;
      if(doubl) continue;
    }

    std::vector< std::pair<DetId, float> > psCells=(*ps)->hitsAndFractions();
    float PS1E=0;
    float PS2E=0;   
    int plane=0;
    for(unsigned int s=0; s<psCells.size(); ++s){
      for(EcalRecHitCollection::const_iterator es=ESRecHits_->begin(); es!= ESRecHits_->end(); ++es){
    if(es->id().rawId()==psCells[s].first.rawId()){ 
      ESDetId id=ESDetId(psCells[s].first.rawId());
      plane=id.plane();
      if(id.plane()==1)PS1E=PS1E + es->energy() * psCells[s].second;
      if(id.plane()==2)PS2E=PS2E + es->energy() * psCells[s].second;
      break;
    }         
      }
    }
    //make PreShower object storing plane PSEnergy and Position
    if(plane==1){
      PreshowerCluster PS1=PreshowerCluster(PS1E,(*ps)->position(),(*ps)->hitsAndFractions(),plane);
      PFPreShowerClust.push_back(PS1);  
    }
    if(plane==2){
      PreshowerCluster PS2=PreshowerCluster(PS2E,(*ps)->position(),(*ps)->hitsAndFractions(),plane);
      PFPreShowerClust.push_back(PS2);      
    }
  }  
  return PFPreShowerClust;  
}

Member Data Documentation

Handle<EcalRecHitCollection> ggPFESClusters::ESRecHits_

private

Definition at line 29 of file ggPFESClusters.h.

Referenced by getPFESClusters().

const CaloSubdetectorGeometry* ggPFESClusters::geomEnd_

private

Definition at line 30 of file ggPFESClusters.h.

Referenced by getLinkDist().