EcalShowerContainmentCorrections.h

Go to the documentation of this file.

#ifndef CondFormats_EcalShowerContainmentCorrections_h
#define CondFormats_EcalShowerContainmentCorrections_h
// -*- C++ -*-
//
// Package:     CondFormats
// Class  :     EcalShowerContainmentCorrections
// 
#include "CondFormats/Serialization/interface/Serializable.h"

#include <vector>
#include <algorithm>
#include <map>
#include <boost/cstdint.hpp>

#include <DataFormats/Math/interface/Point3D.h>

class EBDetId;

class EcalShowerContainmentCorrections {

 public:


  struct Coefficients{

    Coefficients(){for(unsigned int i=0; i<Coefficients::kSize; ++i) data[i]=0;}
    Coefficients(const Coefficients& coeff){
      std::copy(coeff.data,
        coeff.data+Coefficients::kSize,
        data);
    }
 
    Coefficients& operator=(const Coefficients& coeff){
      if (this == &coeff) return *this; 
      std::copy(coeff.data,coeff.data+Coefficients::kSize,data);
      return *this;
    }

    static const int kPolynomialDegree = 4;

    static const int kNTypes           = 8; 
    static const unsigned int kSize    = kPolynomialDegree* kNTypes  ;

    double data[kSize];

  
 COND_SERIALIZABLE;
};



  const Coefficients correctionCoefficients(const EBDetId& centerxtal) const;


  void fillCorrectionCoefficients(const EBDetId& xtal, 
                                  int   group, 
                  const Coefficients&  coefficients);


  void fillCorrectionCoefficients(const int supermodule, const int module, 
                  const Coefficients&  coefficients);



  const double correction3x3(const EBDetId& xtal, 
                 const  math::XYZPoint& pos) const;



  const double correction5x5(const EBDetId& xtal, 
                 const  math::XYZPoint& pos) const;



 private:
  enum  Direction{eX,eY};
  enum  Type{e3x3,e5x5};

  const double correctionXY(const EBDetId& xtal, 
                double position, 
                Direction dir,
                Type type) const ;


  typedef std::map<EBDetId,int> GroupMap;

  GroupMap     groupmap_;


  std::vector<Coefficients> coefficients_;  


 COND_SERIALIZABLE;
};

#endif