/data/refman/pasoursint/CMSSW_4_1_8_patch12/src/DataFormats/METReco/src/EcalHaloData.cc

Go to the documentation of this file.

#include "DataFormats/METReco/interface/EcalHaloData.h"
/*
  [class]:  EcalHaloData
  [authors]: R. Remington, The University of Florida
  [description]: See EcalHaloData.h
  [date]: October 15, 2009
*/

#include <cmath>
using namespace reco;
EcalHaloData::EcalHaloData()
{
}

int EcalHaloData::NumberOfHaloSuperClusters(float roundness, float angle ) const 
{

  int n = 0 ; 
  for( unsigned int i = 0 ; i < TheSuperClusterRefs.size() ; i++ )
    {
      edm::Ref<SuperClusterCollection> cluster(TheSuperClusterRefs, i) ;
      float r = ShowerShapes_Roundness[cluster];
      float a = ShowerShapes_Angle[cluster];
      
      if( (r > 0 && r < roundness) && ( a > 0 && a < angle ) ) 
        n++ ;
    }

  return  n ;
}