35 double deltacur=1000.;
37 bool MATCHEDSC =
false;
39 Direction candDir(position.eta(), position.phi());
44 for(reco::SuperClusterCollection::const_iterator scItr = superClusterCollectionH->begin(); scItr != superClusterCollectionH->end(); ++scItr){
48 if(supercluster->
seed()->algo() == 0){
49 deltacur = ROOT::Math::VectorUtil::DeltaR(supercluster->
position(),
position);
50 if (deltacur < delta1) {
52 matchedsupercluster = supercluster;
61 for(reco::BasicClusterCollection::const_iterator cItr = basicClusterCollectionH->begin(); cItr != basicClusterCollectionH->end(); ++cItr){
65 int ebc_bcalgo = cluster->algo();
66 double ebc_bce = cluster->energy();
67 double ebc_bceta = cluster->eta();
68 double ebc_bcet = ebc_bce*
sin(2*atan(
exp(ebc_bceta)));
71 if (ebc_bcet >
etMin_ && ebc_bcalgo == 0) {
75 bool inSuperCluster =
false;
80 for(;theEclust != matchedsupercluster->
clustersEnd(); ++theEclust) {
81 if ((**theEclust) == (*cluster) ) inSuperCluster =
true;
85 newDelta=ROOT::Math::VectorUtil::DeltaR(cluster->position(),
position);
const math::XYZPoint & position() const
cluster centroid position
void setVeto(const Veto &aVeto)
Set veto.
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Sin< T >::type sin(const T &t)
void addDeposit(double dr, double deposit)
Add deposit (ie. transverse energy or pT)
void addCandEnergy(double et)
Set energy or pT attached to cand trajectory.
XYZPointD XYZPoint
point in space with cartesian internal representation
static int position[264][3]
CaloCluster_iterator clustersBegin() const
fist iterator over BasicCluster constituents
const CaloClusterPtr & seed() const
seed BasicCluster
T const * get() const
Returns C++ pointer to the item.
CaloCluster_iterator clustersEnd() const
last iterator over BasicCluster constituents