CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
ClusterShapeAlgo.h
Go to the documentation of this file.
1 #ifndef RecoEcal_EgammaCoreTools_ClusterShapeAlgo_h
2 #define RecoEcal_EgammaCoreTools_ClusterShapeAlgo_h
3 
13 #include <map>
14 
16 
26 
28 
31  double r;
32  double phi;
33 };
34 
36 public:
39  reco::ClusterShape Calculate(const reco::BasicCluster &passedCluster,
40  const EcalRecHitCollection *hits,
42  const CaloSubdetectorTopology *topology);
43 
44 private:
45  void Calculate_TopEnergy(const reco::BasicCluster &passedCluster, const EcalRecHitCollection *hits);
46  void Calculate_2ndEnergy(const reco::BasicCluster &passedCluster, const EcalRecHitCollection *hits);
47  void Create_Map(const EcalRecHitCollection *hits, const CaloSubdetectorTopology *topology);
48  void Calculate_e2x2();
49  void Calculate_e3x2();
50  void Calculate_e3x3();
51  void Calculate_e4x4();
52  void Calculate_e5x5();
53  void Calculate_e2x5Right();
54  void Calculate_e2x5Left();
55  void Calculate_e2x5Top();
56  void Calculate_e2x5Bottom();
57  void Calculate_Covariances(const reco::BasicCluster &passedCluster,
58  const EcalRecHitCollection *hits,
61  const EcalRecHitCollection *hits,
62  const int EtaPhi,
64  // defines a energy deposition topology in a reference system centered on the cluster
65  void Calculate_EnergyDepTopology(const reco::BasicCluster &passedCluster,
66  const EcalRecHitCollection *hits,
68  bool logW = true);
69  void Calculate_Polynomials(double rho);
70  double factorial(int n) const;
71  void Calculate_lat(const reco::BasicCluster &passedCluster);
72  void Calculate_ComplexZernikeMoments(const reco::BasicCluster &passedCluster);
73  // explicit implementation of polynomial part of
74  // Zernike-Functions for n<=5;
75  double f00(double r);
76  double f11(double r);
77  double f20(double r);
78  double f22(double r);
79  double f31(double r);
80  double f33(double r);
81  double f40(double r);
82  double f42(double r);
83  double f44(double r);
84  double f51(double r);
85  double f53(double r);
86  double f55(double r);
87  double absZernikeMoment(const reco::BasicCluster &passedCluster, int n, int m, double R0 = 6.6);
88  double fast_AbsZernikeMoment(const reco::BasicCluster &passedCluster, int n, int m, double R0);
89  // Calculation of Zernike-Moments for general values of (n,m)
90  double calc_AbsZernikeMoment(const reco::BasicCluster &passedCluster, int n, int m, double R0);
91 
93 
94  std::pair<DetId, double> energyMap_[5][5];
96 
100  double e3x2Ratio_;
101  double lat_;
102  double etaLat_;
103  double phiLat_;
104  double A20_, A42_;
105  std::vector<double> energyBasketFractionEta_;
106  std::vector<double> energyBasketFractionPhi_;
108  std::vector<EcalClusterEnergyDeposition> energyDistribution_;
109  std::vector<double> fcn_;
110 
111  enum { Eta, Phi };
112 };
113 
114 #endif
double f20(double r)
void Create_Map(const EcalRecHitCollection *hits, const CaloSubdetectorTopology *topology)
double fast_AbsZernikeMoment(const reco::BasicCluster &passedCluster, int n, int m, double R0)
std::vector< double > fcn_
std::vector< EcalClusterEnergyDeposition > energyDistribution_
edm::ParameterSet parameterSet_
double f00(double r)
double f44(double r)
void Calculate_BarrelBasketEnergyFraction(const reco::BasicCluster &passedCluster, const EcalRecHitCollection *hits, const int EtaPhi, const CaloSubdetectorGeometry *geometry)
double calc_AbsZernikeMoment(const reco::BasicCluster &passedCluster, int n, int m, double R0)
std::pair< DetId, double > energyMap_[5][5]
double factorial(int n) const
double f31(double r)
double absZernikeMoment(const reco::BasicCluster &passedCluster, int n, int m, double R0=6.6)
double f55(double r)
reco::ClusterShape Calculate(const reco::BasicCluster &passedCluster, const EcalRecHitCollection *hits, const CaloSubdetectorGeometry *geometry, const CaloSubdetectorTopology *topology)
void Calculate_Covariances(const reco::BasicCluster &passedCluster, const EcalRecHitCollection *hits, const CaloSubdetectorGeometry *geometry)
double f42(double r)
void Calculate_EnergyDepTopology(const reco::BasicCluster &passedCluster, const EcalRecHitCollection *hits, const CaloSubdetectorGeometry *geometry, bool logW=true)
double f40(double r)
double f33(double r)
void Calculate_TopEnergy(const reco::BasicCluster &passedCluster, const EcalRecHitCollection *hits)
void Calculate_Polynomials(double rho)
double f11(double r)
Definition: DetId.h:17
double f51(double r)
std::vector< double > energyBasketFractionPhi_
void Calculate_ComplexZernikeMoments(const reco::BasicCluster &passedCluster)
double f22(double r)
void Calculate_lat(const reco::BasicCluster &passedCluster)
std::vector< double > energyBasketFractionEta_
double f53(double r)
void Calculate_2ndEnergy(const reco::BasicCluster &passedCluster, const EcalRecHitCollection *hits)