SuperClusterToCandidate.h

Go to the documentation of this file.

#ifndef RecoAlgos_SuperClusterToCandidate_h
#define RecoAlgos_SuperClusterToCandidate_h
#include "CommonTools/RecoAlgos/interface/MassiveCandidateConverter.h"
#include "CommonTools/RecoAlgos/interface/CandidateProducer.h"
#include "DataFormats/EgammaReco/interface/SuperCluster.h"
#include "DataFormats/RecoCandidate/interface/RecoEcalCandidate.h"
#include "DataFormats/RecoCandidate/interface/RecoEcalCandidateFwd.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"

namespace converter {
  struct SuperClusterToCandidate : public MassiveCandidateConverter {
    typedef reco::SuperCluster value_type;
    typedef reco::SuperClusterCollection Components;
    typedef reco::RecoEcalCandidate Candidate;
    SuperClusterToCandidate(const edm::ParameterSet& cfg, edm::ConsumesCollector iC)
        : MassiveCandidateConverter(cfg, iC) {}
    void convert(reco::SuperClusterRef scRef, reco::RecoEcalCandidate& c) const {
      const reco::SuperCluster& sc = *scRef;
      math::XYZPoint v(0, 0, 0);  // this should be taken from something else...
      math::XYZVector p = sc.energy() * (sc.position() - v).unit();
      double t = sqrt(massSqr_ + p.mag2());
      c.setCharge(0);
      c.setVertex(v);
      c.setP4(reco::Candidate::LorentzVector(p.x(), p.y(), p.z(), t));
      c.setSuperCluster(scRef);
      c.setPdgId(particle_.pdgId());
    }
  };

  namespace helper {
    template <>
    struct CandConverter<reco::SuperCluster> {
      typedef SuperClusterToCandidate type;
    };
  }  // namespace helper
}  // namespace converter

#endif