---

L2TauECALClustering Class Reference

#include <RecoTauTag/HLTProducers/interface/L2TauIsolationAlgs.h>

List of all members.

Public Member Functions
	L2TauECALClustering (double)
	L2TauECALClustering ()
void	run (const math::PtEtaPhiELorentzVectorCollection &, const CaloJet &, L2TauIsolationInfo &)
	~L2TauECALClustering ()
Private Member Functions
void	clusterize (const math::PtEtaPhiELorentzVectorCollection &)
std::vector< double >	clusterSeperation (const CaloJet &) const
Private Attributes
double	m_clusterRadius
L2TauECALClusterCollection	m_clusters

---

Detailed Description

Definition at line 59 of file L2TauIsolationAlgs.h.

---

Constructor & Destructor Documentation

L2TauECALClustering::L2TauECALClustering

(

)

Definition at line 49 of file L2TauIsolationAlgs.cc.

References m_clusterRadius.

{
  m_clusterRadius=0.08;
}

L2TauECALClustering::L2TauECALClustering

(

double

radius

)

Definition at line 54 of file L2TauIsolationAlgs.cc.

References m_clusterRadius.

{
  m_clusterRadius=radius;
}

L2TauECALClustering::~L2TauECALClustering

(

)

Definition at line 59 of file L2TauIsolationAlgs.cc.

{}

---

Member Function Documentation

void L2TauECALClustering::clusterize

(

const math::PtEtaPhiELorentzVectorCollection &

myRecHits

)

[private]

Definition at line 78 of file L2TauIsolationAlgs.cc.

References h, j, m_clusterRadius, m_clusters, and ptr.

Referenced by run().

{
   //If we have Hits do Clustering
   if(myRecHits.size()>0)
     {
       //Create the first Cluster by maximum Crystal
       m_clusters.push_back(L2TauECALCluster(*(myRecHits.begin())));

       //Loop on The Clusters if there are at least two hits
       if(myRecHits.size()>=2)
       for(math::PtEtaPhiELorentzVectorCollection::const_iterator h = myRecHits.begin()+1;h!=myRecHits.end();++h)
       {
          //These vars are used to find the nearest clusters to this hits
          double dR_min=100;  
          int ptr=0;
          int ptr_min=-1;
         
          for(L2TauECALClusterIt j=m_clusters.begin()+1;j!=m_clusters.end();j++)
          {
             if(ROOT::Math::VectorUtil::DeltaR(*h,j->p4())<m_clusterRadius)
                {
                  if(ROOT::Math::VectorUtil::DeltaR(*h,j->p4())<dR_min)
                    {
                      dR_min=ROOT::Math::VectorUtil::DeltaR(*h,j->p4());
                          ptr_min=ptr;
                    }
                }
                  ptr++;
          }
             
          //If it does not belong to cluster add a new one else add the Crystal to the Cluster
          if(ptr_min==-1)
            m_clusters.push_back(L2TauECALCluster(*h));
          else
            m_clusters[ptr_min].addCrystal(*h);

       }
              
     }

}

std::vector< double > L2TauECALClustering::clusterSeperation

(

const CaloJet &

jet

)

const [private]

Definition at line 122 of file L2TauIsolationAlgs.cc.

References c, DeltaPhi(), reco::Particle::eta(), m_clusters, reco::Particle::p4(), and funct::pow().

Referenced by run().

{
  double eta_rms=0.;
  double phi_rms=0.;
  double dr_rms=0.;

  double sumpt = 0.;

  std::vector<double> rmsVector; //declare the vector
  if(m_clusters.size()>0)
    {
      for(L2TauECALClusterIt c = m_clusters.begin();c!=m_clusters.end();++c) //loop on clusters
        {
          eta_rms+=c->p4().Pt()*pow(c->p4().Eta()-jet.eta(),2);
          phi_rms+=c->p4().Pt()*pow(ROOT::Math::VectorUtil::DeltaPhi(c->p4(),jet.p4().Vect()),2);
          dr_rms+=c->p4().Pt()*pow(ROOT::Math::VectorUtil::DeltaR(c->p4(),jet.p4().Vect()),2);
          sumpt+=c->p4().Pt();                     
        }
    }
  else
    {
      eta_rms=0.;
      phi_rms=0.;
      dr_rms =0.;
      sumpt=1.;
    }

  rmsVector.push_back(eta_rms/sumpt);
  rmsVector.push_back(phi_rms/sumpt);
  rmsVector.push_back(dr_rms/sumpt);

  return rmsVector;
}

void L2TauECALClustering::run	(	const math::PtEtaPhiELorentzVectorCollection &	hits,
		const CaloJet &	jet,
		L2TauIsolationInfo &	l2info
	)

Definition at line 63 of file L2TauIsolationAlgs.cc.

References clusterize(), clusterSeperation(), reco::L2TauIsolationInfo::ECALClusterDRRMS, reco::L2TauIsolationInfo::ECALClusterEtaRMS, reco::L2TauIsolationInfo::ECALClusterNClusters, reco::L2TauIsolationInfo::ECALClusterPhiRMS, and m_clusters.

Referenced by L2TauNarrowConeIsolationProducer::produce(), and L2TauIsolationProducer::produce().

{
  //Create Clusters
  clusterize(hits);


  //Fill info Class
  std::vector<double> rms  = clusterSeperation(jet);
  l2info.ECALClusterNClusters=m_clusters.size();
  l2info.ECALClusterEtaRMS=rms[0];
  l2info.ECALClusterPhiRMS=rms[1];
  l2info.ECALClusterDRRMS=rms[2];
}

---

Member Data Documentation

double L2TauECALClustering::m_clusterRadius [private]

Definition at line 75 of file L2TauIsolationAlgs.h.

Referenced by clusterize(), and L2TauECALClustering().

L2TauECALClusterCollection L2TauECALClustering::m_clusters [private]

Definition at line 78 of file L2TauIsolationAlgs.h.

Referenced by clusterize(), clusterSeperation(), and run().

---

The documentation for this class was generated from the following files:

• RecoTauTag/HLTProducers/interface/L2TauIsolationAlgs.h
• RecoTauTag/HLTProducers/src/L2TauIsolationAlgs.cc

---