CMSSW Reference Manual

00001 
00002 #ifndef HFEMClusterShape_h
00003 #define HFEMClusterShape_h
00004 
00005 #include <Rtypes.h>
00006 #include "DataFormats/EgammaReco/interface/HFEMClusterShapeFwd.h"
00007 #include "DataFormats/DetId/interface/DetId.h"
00008 /*class reco::HFEMClusterShape HFEMClusterShape.h DataFormats/EgammaReco/interface/HFEMClusterShape.h
00009  *  
00010  * A Cluster Shape of a Possible EM cluster in the HF detector
00011  * Contains the DetId of its Seed
00012  *
00013  * \author Kevin Klapoetke, University of Minnesota
00014  *
00015  * \version $Id: HFEMClusterShape.h,v 1.3 2007/10/08 18:52:13 futyand Exp $
00016  *
00017  */
00018 
00019 namespace reco {
00020 
00021   class HFEMClusterShape {
00022   public:
00023     HFEMClusterShape() { }
00024    
00025     HFEMClusterShape(double eLong1x1,double eShort1x1,
00026                      double eLong3x3,double eShort3x3,double eLong5x5,
00027                      double eShort5x5,double eLongCore,
00028                      double CellEta,double CellPhi,
00029                      DetId seed);
00030     
00031     
00032     //energy in long or short fibers various cluster sizes
00033     double eLong1x1() const {return eLong1x1_;}
00034     double eShort1x1() const {return eShort1x1_;}
00035     double eLong3x3() const {return eLong3x3_;}
00036     double eShort3x3() const {return eShort3x3_;}
00037     double eLong5x5() const {return eLong5x5_;}
00038     double eShort5x5() const {return eShort5x5_;}
00039     
00040     //total energy in various clusters
00041     double e1x1() const;
00042     double e3x3() const;
00043     double e5x5() const;
00044 
00045     //Identification Variables
00046     //Longetudinal variable: E(3x3,short fibers)/E(3x3,long fibers)
00047     double eSeL() const;
00048     //Transverse Variable: E(Core of cluster)/E(3x3)
00049     double eCOREe9() const;
00050     //Shower Exclusion Variable: E(3x3)/E(5x5)
00051     double e9e25() const;
00052 
00053     //energy in central highest energy cells (at least 50% energy of previous total energy startign with seed cell)
00054     double eCore() const {return eLongCore_;}
00055   
00056     double CellEta() const {return CellEta_;}
00057     double CellPhi() const {return CellPhi_;}
00058     
00059     //seed cell of cluster DetId
00060     DetId seed() const {return seed_;}
00061 
00062   private:
00063     double eLong1x1_, eShort1x1_, eLong3x3_,eShort3x3_, eLong5x5_, eShort5x5_,  eLongCore_,CellEta_,CellPhi_;
00064     DetId seed_;
00065  
00066   };
00067   
00068 }
00069 
00070 
00071 #endif