HcalCellType.cc

Go to the documentation of this file.

// File: HcalCellType.cc
// Description: Individual readout cell description for hadron calorimeter

#include "Geometry/HcalCommonData/interface/HcalCellType.h"

#include "CLHEP/Units/GlobalSystemOfUnits.h"
#include <iomanip>

HcalCellType::HcalCellType(HcalSubdetector detType, int etaBin, int phiBin, 
               int depthSegment,const HcalCellType::HcalCell& cell, 
               int readoutDirection, double samplingFactor,
               int numberZ, int nmodule, double halfSize, 
               int units) :
  theDetType(detType), theEtaBin(etaBin), theDepthSegment(depthSegment),
  theNumberOfZ(numberZ), theActualReadoutDirection(readoutDirection), 
  theUnitPhi(units), theSamplingFactor(samplingFactor){

  theEtaMin   = cell.eta - cell.deta;
  theEtaMax   = cell.eta + cell.deta;
  theRzFlag   = cell.flagrz;
  theDepthMin = (cell.rz  - cell.drz)/CLHEP::cm;
  theDepthMax = (cell.rz  + cell.drz)/CLHEP::cm;
  int nphi           = (int)(10*CLHEP::deg/cell.dphi);
  theNumberOfPhiBins = nphi*nmodule;
  double phimin      = cell.phi - cell.dphi;
  thePhiOffset       = (phimin - 2*(phiBin-1)*cell.dphi)/CLHEP::deg;
  thePhiBinWidth     = 2*(cell.dphi)/CLHEP::deg;
  theHalfSize        = halfSize/CLHEP::cm;
}

HcalCellType::HcalCellType(const HcalCellType &right) {

  theDetType                = right.theDetType;
  theEtaBin                 = right.theEtaBin;
  theDepthSegment           = right.theDepthSegment;
  theNumberOfPhiBins        = right.theNumberOfPhiBins;
  theNumberOfZ              = right.theNumberOfZ;
  theActualReadoutDirection = right.theActualReadoutDirection;
  theUnitPhi                = right.theUnitPhi;
  theRzFlag                 = right.theRzFlag;
  theEtaMin                 = right.theEtaMin;
  theEtaMax                 = right.theEtaMax;
  thePhiOffset              = right.thePhiOffset;
  thePhiBinWidth            = right.thePhiBinWidth;
  theDepthMin               = right.theDepthMin;
  theDepthMax               = right.theDepthMax;
  theHalfSize               = right.theHalfSize;
  theSamplingFactor         = right.theSamplingFactor;
  theMissingPhiPlus         = right.theMissingPhiPlus;
  theMissingPhiMinus        = right.theMissingPhiMinus;
}

const HcalCellType& HcalCellType::operator=(const HcalCellType &right) {

  theDetType                = right.theDetType;
  theEtaBin                 = right.theEtaBin;
  theDepthSegment           = right.theDepthSegment;
  theNumberOfPhiBins        = right.theNumberOfPhiBins;
  theNumberOfZ              = right.theNumberOfZ;
  theActualReadoutDirection = right.theActualReadoutDirection;
  theUnitPhi                = right.theUnitPhi;
  theRzFlag                 = right.theRzFlag;
  theEtaMin                 = right.theEtaMin;
  theEtaMax                 = right.theEtaMax;
  thePhiOffset              = right.thePhiOffset;
  thePhiBinWidth            = right.thePhiBinWidth;
  theDepthMin               = right.theDepthMin;
  theDepthMax               = right.theDepthMax;
  theHalfSize               = right.theHalfSize;
  theSamplingFactor         = right.theSamplingFactor;
  theMissingPhiPlus         = right.theMissingPhiPlus;
  theMissingPhiMinus        = right.theMissingPhiMinus;

  return *this;
}

HcalCellType::~HcalCellType() {}

void HcalCellType::setMissingPhi(std::vector<int>& v1, std::vector<int>& v2) {
  theMissingPhiPlus         = v1;
  theMissingPhiMinus        = v2;
}

int HcalCellType::nPhiMissingBins() const {
  int tmp = (int)(theMissingPhiPlus.size());
  if (theNumberOfZ > 1)  tmp += (int)(theMissingPhiMinus.size());
  return tmp;
}

void HcalCellType::setEta(int bin, double etamin, double etamax) {
  theEtaBin = bin;
  theEtaMin = etamin;
  theEtaMax = etamax;
}

void HcalCellType::setDepth(int bin, double dmin, double dmax) {
  theDepthSegment = bin;
  theDepthMin     = dmin;
  theDepthMax     = dmax;
}

void HcalCellType::setPhi(int bins, int unit, double dphi, double phioff) {
  theNumberOfPhiBins = bins;
  theUnitPhi         = unit;
  thePhiBinWidth     = dphi;
  thePhiOffset       = phioff;
}

std::ostream& operator<<(std::ostream& os, const HcalCellType& cell) {
  os << "Detector " << cell.detType() << " Eta " << cell.etaBin() << " (" 
     << cell.etaMin() << ":" << cell.etaMax() << ") Depth " 
     << cell.depthSegment() << " (" << cell.depthMin() << ":" 
     << cell.depthMax() << "; " << cell.depthType() << ") Phi " 
     << cell.nPhiBins() << " (" << cell.phiOffset() << ", "
     << cell.phiBinWidth() << ", " << cell.nPhiModule() << ", "
     << cell.unitPhi() << ") Halves " << cell.nHalves() << " Direction " 
     << cell.actualReadoutDirection() << " Half size " << cell.halfSize() 
     << " Sampling Factor " << cell.samplingFactor() << " # of missing cells "
     << cell.missingPhiPlus().size() << "/" << cell.missingPhiMinus().size();
  return os;
}