Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 12:53:19

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