HcalCellType.h

Go to the documentation of this file.

// File: HcalCellType.h
// Description: Hcal readout cell definition given by eta boundary and depth
#ifndef HcalCellType_h
#define HcalCellType_h

#include <vector>
#include <iostream>
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"

class HcalCellType {

public:

  struct HcalCell {
    bool   ok;
    double eta, deta, phi, dphi, rz, drz;
    bool   flagrz;
    HcalCell(bool fl=false, double et=0, double det=0, double fi=0,
             double dfi=0, double rzv=0, double drzv=0, bool frz=true) :
      ok(fl), eta(et), deta(det), phi(fi), dphi(dfi), rz(rzv), drz(drzv),
         flagrz(frz) {}
  };

  HcalCellType(HcalSubdetector detType, int etaBin, int phiBin, 
           int depthSegment, const HcalCell& cell, int readoutDirection=0,
           double samplingFactor=0, int numberZ=0, int nmodule=0,
           double halfSize=0, int units=0);
  HcalCellType(const HcalCellType &right);
  const HcalCellType& operator=(const HcalCellType &right);
  ~HcalCellType();

  void setMissingPhi(std::vector<int>&, std::vector<int>&);

  HcalSubdetector detType() const {return theDetType;}
                                                                               
  int etaBin() const {return theEtaBin;}
  void setEta(int bin, double etamin, double etamax);
                                                                               
  int depthSegment() const {return theDepthSegment;}
  void setDepth(int bin, double dmin, double dmax);

  int nPhiBins() const {return theNumberOfPhiBins;}
  int nPhiModule() const {return static_cast<int>(20./thePhiBinWidth);}
                                                                               
  double phiBinWidth() const {return thePhiBinWidth;}

  double phiOffset() const {return thePhiOffset;}
  int    unitPhi() const {return theUnitPhi;}
  void   setPhi(int bins, int unit, double dphi, double phioff);

  int nHalves() const {return theNumberOfZ;}
                                                                               
  int actualReadoutDirection() const {return theActualReadoutDirection;}
                                                                               
  double etaMin() const {return theEtaMin;}
                                                                               
  double etaMax() const {return theEtaMax;}

  double depth() const {return (theDepthMin+theDepthMax)/2;}
  double depthMin() const {return theDepthMin;}
  double depthMax() const {return theDepthMax;}
  bool   depthType() const {return theRzFlag;}
  double halfSize() const {return theHalfSize;}
                                                                               
  double samplingFactor() const {return theSamplingFactor;}

  std::vector<int> missingPhiPlus()  const {return theMissingPhiPlus;}
  std::vector<int> missingPhiMinus() const {return theMissingPhiMinus;}
  int nPhiMissingBins() const;

protected:
 
  HcalCellType();
                                                                               
private:

  HcalSubdetector  theDetType;
  int              theEtaBin;
  int              theDepthSegment;
  int              theNumberOfPhiBins;
  int              theNumberOfZ;
  int              theActualReadoutDirection;
  int              theUnitPhi;

  bool             theRzFlag;
                                                                               
  double           theEtaMin;
  double           theEtaMax;
  double           thePhiOffset;
  double           thePhiBinWidth;
  double           theDepthMin;
  double           theDepthMax;
  double           theHalfSize;
  double           theSamplingFactor;

  std::vector<int> theMissingPhiPlus;
  std::vector<int> theMissingPhiMinus;
};
                                                                               
std::ostream& operator<<(std::ostream&, const HcalCellType&);
#endif