L1AnalysisRecoCluster.cc

Go to the documentation of this file.

#include "DataFormats/Common/interface/Handle.h"
#include "L1Trigger/L1TNtuples/interface/L1AnalysisRecoCluster.h"
#include "DataFormats/EgammaReco/interface/SuperCluster.h"

L1Analysis::L1AnalysisRecoCluster::L1AnalysisRecoCluster()
{
}


L1Analysis::L1AnalysisRecoCluster::~L1AnalysisRecoCluster()
{
}

void L1Analysis::L1AnalysisRecoCluster::Set(const reco::CaloClusterCollection & caloClusters, unsigned maxCl)
{
  recoCluster_.nClusters=recoCluster_.eta.size();
  for(reco::CaloClusterCollection::const_iterator it=caloClusters.begin(); it!=caloClusters.end() && recoCluster_.nClusters<maxCl; it++)
  {
     recoCluster_.eta.push_back( it->eta() );
     recoCluster_.phi.push_back( it->phi() );
     recoCluster_.et.push_back ( it->energy() * sin( it->position().theta() ) );
     recoCluster_.e.push_back  ( it->energy() );
     recoCluster_.nClusters++;
  }

}

void L1Analysis::L1AnalysisRecoCluster::Set(const reco::SuperClusterCollection & superClusters, unsigned maxCl)
{
   recoCluster_.nClusters=recoCluster_.eta.size();

  for(reco::SuperClusterCollection::const_iterator it=superClusters.begin(); it!=superClusters.end() && recoCluster_.nClusters<maxCl; it++)
  {
     recoCluster_.eta.push_back(it->eta());
     recoCluster_.phi.push_back(it->phi());
     recoCluster_.et.push_back (it->energy() * sin( it->position().theta() ) );
     recoCluster_.e.push_back  (it->energy());
     recoCluster_.nClusters++;
   }

}