FindDistCone.h

Go to the documentation of this file.

#ifndef CalibrationIsolatedParticlesFindDistCone_h
#define CalibrationIsolatedParticlesFindDistCone_h

// system include files
#include <memory>
#include <cmath>
#include <string>
#include <map>
#include <vector>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"

#include "DataFormats/Common/interface/Handle.h"

#include "DataFormats/Candidate/interface/Candidate.h"

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"

#include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"

#include <cmath>

namespace spr {

  // Cone clustering core
  double getDistInPlaneTrackDir(const GlobalPoint& caloPoint, const GlobalVector& caloVector, const GlobalPoint& rechitPoint, bool debug=false);

  // Not used, but here for reference
  double getDistInCMatEcal(double eta1, double phi1, double eta2, double phi2, bool debug=false);
  double getDistInCMatHcal(double eta1, double phi1, double eta2, double phi2, bool debug=false);

  // get eta, phi, energy of rechits in collection
  void getEtaPhi(HBHERecHitCollection::const_iterator hit, std::vector<int>& RH_ieta, std::vector<int>& RH_iphi, std::vector<double>& RH_ene, bool debug=false);

  void getEtaPhi(edm::PCaloHitContainer::const_iterator hit, std::vector<int>& RH_ieta, std::vector<int>& RH_iphi, std::vector<double>& RH_ene, bool debug=false);

  void getEtaPhi(EcalRecHitCollection::const_iterator hit, std::vector<int>& RH_ieta, std::vector<int>& RH_iphi, std::vector<double>& RH_ene, bool debug=false);

  // get eta, phi of rechits in collection
  void getEtaPhi(HBHERecHitCollection::const_iterator hit,int& ieta,int& iphi, bool debug=false);
  void getEtaPhi(edm::PCaloHitContainer::const_iterator hit,int& ieta,int& iphi, bool debug=false);
  void getEtaPhi(EcalRecHitCollection::const_iterator hit,int& ieta,int& iphi, bool debug=false);

  double getEnergy(HBHERecHitCollection::const_iterator hit, int useRaw=0, bool debug=false);
  double getEnergy(edm::PCaloHitContainer::const_iterator hit, int useRaw=0, bool debug=false);
  double getEnergy(EcalRecHitCollection::const_iterator hit, int useRaw=0, bool debug=false);

  GlobalPoint getGpos(const CaloGeometry* geo, HBHERecHitCollection::const_iterator hit, bool debug=false);
  GlobalPoint getGpos(const CaloGeometry* geo, edm::PCaloHitContainer::const_iterator hit, bool debug=false);
  GlobalPoint getGpos(const CaloGeometry* geo, EcalRecHitCollection::const_iterator hit, bool debug=false);

  double getRawEnergy(HBHERecHitCollection::const_iterator hit, int useRaw=0);
  double getRawEnergy(edm::PCaloHitContainer::const_iterator hit, int useRaw=0);
  double getRawEnergy(EcalRecHitCollection::const_iterator hit, int useRaw=0);
  
}

#endif