Main Page Namespace List Class Hierarchy Compound List File List Namespace Members Compound Members File Members

Atlfast::CellAssociatedClusterEnergy Class Reference

Collaboration diagram for Atlfast::CellAssociatedClusterEnergy:


Public Methods
	CellAssociatedClusterEnergy (ClusterCollection clusters, std::map< Cluster , double > &cluMap, double rConeBarrel)
void	operator() (IKinematic *cell)
Private Attributes
ClusterCollection *	m_clusters
std::map< Cluster *, double > &	m_cluMap
double	m_rConeBarrel

Constructor & Destructor Documentation

Atlfast::CellAssociatedClusterEnergy::CellAssociatedClusterEnergy	(	ClusterCollection *	clusters,
		std::map< Cluster *, double > &	cluMap,
		double	rConeBarrel
	)	`[inline]`

Definition at line 54 of file SharedConeStrategy.cxx.

References m_cluMap, m_clusters, and m_rConeBarrel.

00058                                  :
00059       m_clusters(clusters), m_cluMap(cluMap), m_rConeBarrel(rConeBarrel){
00060     }
    void operator()(IKinematic* cell){

Member Function Documentation

void Atlfast::CellAssociatedClusterEnergy::operator() ( IKinematic * cell ) [inline]

Definition at line 61 of file SharedConeStrategy.cxx.

References m_cluMap, m_clusters, and m_rConeBarrel.

00061                                      {
00062       
00063       
00064       // find all clusters within dR of the cell
00065       ClusterCollection::iterator endAssClu = 
00066         partition(m_clusters->begin(), 
00067                   m_clusters->end(), 
00068                   PartitionCondition::BelowThresholdDeltaR(*cell,m_rConeBarrel)
00069                   );
00070       
00071       //    assert(endAssClu != cluStart);
00072       
00073       //Find the sum of energies of all clusters within dR 
00074       // - will be used to calculate cluster weight
00075       //
00076       // in conversion from Fortran: use energy instead of 
00077       // eT*cosh(pseudoRapidity)
00078       
00079       double clusterEsum = 
00080         std::accumulate(m_clusters->begin(), endAssClu, 0., EsumIK() );
00081             
00082       //calculate the shared energy fractions
00083       
00084       ClusterCollection::iterator itr = m_clusters->begin();
00085       for(; itr != endAssClu; ++itr){
00086         
00087         double eCell         = cell->momentum().e();
00088         double eCluster      = ( *itr)->momentum().e();
00089         m_cluMap[*itr]      += eCell*eCluster/clusterEsum; 
00090         
00091       }
00092 
00093       return;
00094     }