EgammaReco/interface/HFEMClusterShape.h

0001
0002 #ifndef HFEMClusterShape_h
0003 #define HFEMClusterShape_h
0004
0005 #include <Rtypes.h>
0006 #include "DataFormats/EgammaReco/interface/HFEMClusterShapeFwd.h"
0007 #include "DataFormats/DetId/interface/DetId.h"
0008 /*class reco::HFEMClusterShape HFEMClusterShape.h DataFormats/EgammaReco/interface/HFEMClusterShape.h
0009  *
0010  * A Cluster Shape of a Possible EM cluster in the HF detector
0011  * Contains the DetId of its Seed
0012  *
0013  * \author Kevin Klapoetke, University of Minnesota
0014  *
0015  *
0016  */
0017
0018 namespace reco {
0019
0020   class HFEMClusterShape {
0021   public:
0022     HFEMClusterShape() {}
0023
0024     HFEMClusterShape(double eLong1x1,
0025                      double eShort1x1,
0026                      double eLong3x3,
0027                      double eShort3x3,
0028                      double eLong5x5,
0029                      double eShort5x5,
0030                      double eLongCore,
0031                      double CellEta,
0032                      double CellPhi,
0033                      DetId seed);
0034
0035     //energy in long or short fibers various cluster sizes
0036     double eLong1x1() const { return eLong1x1_; }
0037     double eShort1x1() const { return eShort1x1_; }
0038     double eLong3x3() const { return eLong3x3_; }
0039     double eShort3x3() const { return eShort3x3_; }
0040     double eLong5x5() const { return eLong5x5_; }
0041     double eShort5x5() const { return eShort5x5_; }
0042
0043     //total energy in various clusters
0044     double e1x1() const;
0045     double e3x3() const;
0046     double e5x5() const;
0047
0048     //Identification Variables
0049     //Longetudinal variable: E(3x3,short fibers)/E(3x3,long fibers)
0050     double eSeL() const;
0051     //Transverse Variable: E(Core of cluster)/E(3x3)
0052     double eCOREe9() const;
0053     //Shower Exclusion Variable: E(3x3)/E(5x5)
0054     double e9e25() const;
0055
0056     //energy in central highest energy cells (at least 50% energy of previous total energy startign with seed cell)
0057     double eCore() const { return eLongCore_; }
0058
0059     double CellEta() const { return CellEta_; }
0060     double CellPhi() const { return CellPhi_; }
0061
0062     //seed cell of cluster DetId
0063     DetId seed() const { return seed_; }
0064
0065   private:
0066     double eLong1x1_, eShort1x1_, eLong3x3_, eShort3x3_, eLong5x5_, eShort5x5_, eLongCore_, CellEta_, CellPhi_;
0067     DetId seed_;
0068   };
0069
0070 }  // namespace reco
0071
0072 #endif