FastTimeDDDConstants Class Reference

#include <FastTimeDDDConstants.h>

Public Member Functions
	FastTimeDDDConstants (const FastTimeParameters *ftp)
int	getCells (int type) const
std::vector< GlobalPoint >	getCorners (int type, int izeta, int iphi, int zside) const
std::pair< int, int >	getEtaPhi (double r, double phi) const
GlobalPoint	getPosition (int type, int izeta, int iphi, int zside) const
double	getRin (int type) const
double	getRout (int type) const
double	getZHalf (int type) const
std::pair< int, int >	getZPhi (double z, double phi) const
double	getZPos (int type) const
bool	isValidXY (int type, int izeta, int iphi) const
int	numberEtaZ (int type) const
int	numberPhi (int type) const
	~FastTimeDDDConstants ()

Private Member Functions
void	initialize ()

Private Attributes
double	dEta_
double	dPhiBarrel_
double	dPhiEndcap_
double	dZBarrel_
double	etaMax_
double	etaMin_
const FastTimeParameters *	ftpar_
std::vector< double >	rLimits_

Detailed Description

this class reads the constant section of the numbering xml-file for fast timer device

Author

Sunanda Banerjee, SINP sunan.nosp@m.da.b.nosp@m.anerj.nosp@m.ee@c.nosp@m.ern.c.nosp@m.h

Definition at line 19 of file FastTimeDDDConstants.h.

Constructor & Destructor Documentation

FastTimeDDDConstants::FastTimeDDDConstants

(

const FastTimeParameters *

ftp

)

Definition at line 10 of file FastTimeDDDConstants.cc.

References initialize().

    : ftpar_(ft) {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") 
    << "FastTimeDDDConstants::FastTimeDDDConstants "
    << "( const FastTimeParameters* ft ) constructor";
#endif
  initialize();
}

FastTimeDDDConstants::~FastTimeDDDConstants

(

)

Definition at line 20 of file FastTimeDDDConstants.cc.

                                            {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "FastTimeDDDConstants:destructed!!!";
#endif
}

Member Function Documentation

int FastTimeDDDConstants::getCells

(

int

type

)

const

Definition at line 120 of file FastTimeDDDConstants.cc.

References ftpar_, FastTimeParameters::nEtaEndcap_, FastTimeParameters::nPhiBarrel_, FastTimeParameters::nPhiEndcap_, and FastTimeParameters::nZBarrel_.

Referenced by initialize().

                                                 {
  int numb(0);
  if (type == 1) {
    numb = (ftpar_->nZBarrel_) * (ftpar_->nPhiBarrel_);
  } else if (type == 2) {
    numb = (ftpar_->nEtaEndcap_) * (ftpar_->nPhiEndcap_);
  }
  return numb;
}

std::vector< GlobalPoint > FastTimeDDDConstants::getCorners	(	int	type,
		int	izeta,
		int	iphi,
		int	zside
	)		const

Definition at line 82 of file FastTimeDDDConstants.cc.

References funct::cos(), dPhiBarrel_, dPhiEndcap_, PVValHelper::dx, PVValHelper::dz, dZBarrel_, ftpar_, FastTimeParameters::geomParBarrel_, FastTimeParameters::geomParEndcap_, mps_fire::i, createfilelist::int, AlCaHLTBitMon_ParallelJobs::p, phi, RPCpg::pts, alignCSCRings::r, rLimits_, funct::sin(), x, y, and z.

Referenced by FastTimeGeometryLoader::buildGeom(), and FastTimeGeometry::getCorners().

                                                                           {
  double x(0), y(0), z(0), dx(0), dz(0), r(0), phi(0);
  if (type == 1) {
    phi = (iphi - 0.5) * dPhiBarrel_;
    r = ftpar_->geomParBarrel_[2];
    x = r * cos(phi);
    y = r * sin(phi);
    z = ftpar_->geomParBarrel_[0] + (izeta - 0.5) * dZBarrel_;
    dx = 0.5 * ftpar_->geomParBarrel_[3];
    dz = 0.5 * dZBarrel_;
  } else if (type == 2) {
    phi = (iphi - 0.5) * dPhiEndcap_;
    r = (izeta <= 0 || izeta >= (int)(rLimits_.size()))
            ? 0
            : 0.5 * (rLimits_[izeta - 1] + rLimits_[izeta]);
    x = (zside < 0) ? -r * cos(phi) : r * cos(phi);
    y = r * sin(phi);
    z = ftpar_->geomParEndcap_[2];
    dx = 0.5 * r * dPhiEndcap_;
    dz = 0.5 * ftpar_->geomParEndcap_[3];
  }
  if (zside < 0) {
    z = -z;
    dz = -dz;
  }
  static const int signx[8] = {-1, -1, 1, 1, -1, -1, 1, 1};
  static const int signy[8] = {-1, 1, 1, -1, -1, 1, 1, -1};
  static const int signz[8] = {-1, -1, -1, -1, 1, 1, 1, 1};
  std::vector<GlobalPoint> pts;
  for (unsigned int i = 0; i != 8; ++i) {
    GlobalPoint p(x + signx[i] * dx, y + signy[i] * dx, z + signz[i] * dz);
    pts.emplace_back(p);
  }
  return pts;
}

std::pair< int, int > FastTimeDDDConstants::getEtaPhi	(	double	r,
		double	phi
	)		const

Definition at line 40 of file FastTimeDDDConstants.cc.

References geant_units::operators::convertRadToDeg(), dPhiEndcap_, ftpar_, createfilelist::int, gen::k, FastTimeParameters::nEtaEndcap_, FastTimeParameters::nPhiEndcap_, geant_units::piRadians(), and rLimits_.

Referenced by FastTimeGeometry::getClosestCell(), and FastTimerSD::setDetUnitId().

                                                                      {
  if (phi < 0) phi += (2*geant_units::piRadians);
  int ir(ftpar_->nEtaEndcap_);
  for (unsigned int k = 1; k < rLimits_.size(); ++k) {
    if (r > rLimits_[k]) {
      ir = k;
      break;
    }
  }
  int iphi = (int)(phi / dPhiEndcap_) + 1;
  if (iphi > ftpar_->nPhiEndcap_) iphi = 1;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "FastTimeDDDConstants:Endcap r|phi " << r 
                << " " << convertRadToDeg(phi) << " ir|iphi " 
                << ir << " " << iphi;
#endif
  return std::pair<int, int>(ir, iphi);
}

GlobalPoint FastTimeDDDConstants::getPosition	(	int	type,
		int	izeta,
		int	iphi,
		int	zside
	)		const

Definition at line 60 of file FastTimeDDDConstants.cc.

References funct::cos(), dPhiBarrel_, dPhiEndcap_, dZBarrel_, ftpar_, FastTimeParameters::geomParBarrel_, FastTimeParameters::geomParEndcap_, createfilelist::int, AlCaHLTBitMon_ParallelJobs::p, phi, alignCSCRings::r, rLimits_, funct::sin(), x, y, and z.

Referenced by FastTimeGeometry::getPosition().

                                                               {
  double x(0), y(0), z(0);
  if (type == 1) {
    double phi = (iphi - 0.5) * dPhiBarrel_;
    x = ftpar_->geomParBarrel_[2] * cos(phi);
    y = ftpar_->geomParBarrel_[2] * sin(phi);
    z = ftpar_->geomParBarrel_[0] + (izeta - 0.5) * dZBarrel_;
  } else if (type == 2) {
    double phi = (iphi - 0.5) * dPhiEndcap_;
    double r = (izeta <= 0 || izeta >= (int)(rLimits_.size()))
                   ? 0
                   : 0.5 * (rLimits_[izeta - 1] + rLimits_[izeta]);
    x = (zside < 0) ? -r * cos(phi) : r * cos(phi);
    y = r * sin(phi);
    z = ftpar_->geomParEndcap_[2];
  }
  if (zside < 0) z = -z;
  GlobalPoint p(x, y, z);
  return p;
}

double FastTimeDDDConstants::getRin

(

int

type

)

const

Definition at line 130 of file FastTimeDDDConstants.cc.

References ftpar_, FastTimeParameters::geomParBarrel_, FastTimeParameters::geomParEndcap_, and relativeConstraints::value.

Referenced by FastTimeGeometryLoader::build(), and FastTimeGeometry::getSummary().

                                                  {
  double value(0);
  if (type == 1) {
    value = (ftpar_->geomParBarrel_[2]);
  } else if (type == 2) {
    value = (ftpar_->geomParEndcap_[0]);
  }
  return value;
}

double FastTimeDDDConstants::getRout

(

int

type

)

const

Definition at line 140 of file FastTimeDDDConstants.cc.

References ftpar_, FastTimeParameters::geomParBarrel_, FastTimeParameters::geomParEndcap_, and relativeConstraints::value.

Referenced by FastTimeGeometryLoader::build(), and FastTimeGeometry::getSummary().

                                                   {
  double value(0);
  if (type == 1) {
    value = (ftpar_->geomParBarrel_[2] + ftpar_->geomParBarrel_[3]);
  } else if (type == 2) {
    value = (ftpar_->geomParEndcap_[1]);
  }
  return value;
}

double FastTimeDDDConstants::getZHalf

(

int

type

)

const

Definition at line 150 of file FastTimeDDDConstants.cc.

References ftpar_, FastTimeParameters::geomParBarrel_, FastTimeParameters::geomParEndcap_, and relativeConstraints::value.

Referenced by FastTimeGeometryLoader::build(), and FastTimeGeometry::getSummary().

                                                    {
  double value(0);
  if (type == 1) {
    value = 0.5 * (ftpar_->geomParBarrel_[1] - ftpar_->geomParBarrel_[0]);
  } else if (type == 2) {
    value = (ftpar_->geomParEndcap_[3]);
  }
  return value;
}

std::pair< int, int > FastTimeDDDConstants::getZPhi	(	double	z,
		double	phi
	)		const

Definition at line 26 of file FastTimeDDDConstants.cc.

References geant_units::operators::convertRadToDeg(), dPhiBarrel_, dZBarrel_, ftpar_, createfilelist::int, FastTimeParameters::nPhiBarrel_, FastTimeParameters::nZBarrel_, and geant_units::piRadians().

Referenced by FastTimeGeometry::getClosestCell(), and FastTimerSD::setDetUnitId().

                                                                          {
  if (phi < 0) phi += (2*geant_units::piRadians);
  int iz = (int)(z / dZBarrel_) + 1;
  if (iz > ftpar_->nZBarrel_) iz = ftpar_->nZBarrel_;
  int iphi = (int)(phi / dPhiBarrel_) + 1;
  if (iphi > ftpar_->nPhiBarrel_) iphi = 1;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "FastTimeDDDConstants:Barrel z|phi " << z 
                << " " << convertRadToDeg(phi) << " iz|iphi "
                << iz << " " << iphi;
#endif
  return std::pair<int, int>(iz, iphi);
}

double FastTimeDDDConstants::getZPos

(

int

type

)

const

Definition at line 160 of file FastTimeDDDConstants.cc.

References ftpar_, FastTimeParameters::geomParBarrel_, FastTimeParameters::geomParEndcap_, and relativeConstraints::value.

Referenced by FastTimeGeometryLoader::build().

                                                   {
  double value(0);
  if (type == 1) {
    value = 0.5 * (ftpar_->geomParBarrel_[1] + ftpar_->geomParBarrel_[0]);
  } else if (type == 2) {
    value = (ftpar_->geomParEndcap_[2]);
  }
  return value;
}

void FastTimeDDDConstants::initialize ( void  )

private

Definition at line 202 of file FastTimeDDDConstants.cc.

References geant_units::operators::convertRadToDeg(), dEta_, dPhiBarrel_, dPhiEndcap_, dZBarrel_, EDM_ML_DEBUG, PVValHelper::eta, etaMax_, etaMin_, JetChargeProducer_cfi::exp, ftpar_, FastTimeParameters::geomParBarrel_, FastTimeParameters::geomParEndcap_, getCells(), gen::k, cmsBatch::log, FastTimeParameters::nEtaEndcap_, FastTimeParameters::nPhiBarrel_, FastTimeParameters::nPhiEndcap_, FastTimeParameters::nZBarrel_, geant_units::piRadians(), rLimits_, funct::tan(), and theta().

Referenced by FastTimeDDDConstants().

                                      {
  double thmin = atan(ftpar_->geomParEndcap_[0] / ftpar_->geomParEndcap_[2]);
  etaMax_ = -log(0.5 * thmin);
  double thmax = atan(ftpar_->geomParEndcap_[1] / ftpar_->geomParEndcap_[2]);
  etaMin_ = -log(0.5 * thmax);
  dEta_ = (etaMax_ - etaMin_) / ftpar_->nEtaEndcap_;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "Theta range " << convertRadToDeg(thmin) 
                << ":" << convertRadToDeg(thmax) << " Eta range "
                << etaMin_ << ":" << etaMax_ << ":" << dEta_;
#endif
  for (int k = 0; k <= ftpar_->nEtaEndcap_; ++k) {
    double eta = etaMin_ + k * dEta_;
    double theta = 2.0 * atan(exp(-eta));
    double rval = (ftpar_->geomParEndcap_[2]) * tan(theta);
    rLimits_.emplace_back(rval);
  }
  dZBarrel_ = ftpar_->geomParBarrel_[1] / ftpar_->nZBarrel_;
  dPhiBarrel_ = (2*geant_units::piRadians) / ftpar_->nPhiBarrel_;
  dPhiEndcap_ = (2*geant_units::piRadians) / ftpar_->nPhiEndcap_;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "FastTimeDDDConstants initialized with "
                << ftpar_->nZBarrel_ << ":" << ftpar_->nPhiBarrel_
                << ":" << getCells(1)
                << " cells for barrel; dz|dphi " << dZBarrel_ 
                << "|" << dPhiBarrel_ << " and " 
                << ftpar_->nEtaEndcap_ << ":" 
                << ftpar_->nPhiEndcap_ << ":" << getCells(2) 
                << " cells for endcap; dphi " << dPhiEndcap_
                << " The Limits in R are";
  for (unsigned int k = 0; k < rLimits_.size(); ++k)
    edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << rLimits_[k] << " ";
#endif
}

bool FastTimeDDDConstants::isValidXY	(	int	type,
		int	izeta,
		int	iphi
	)		const

Definition at line 170 of file FastTimeDDDConstants.cc.

References convertSQLiteXML::ok.

Referenced by FastTimeTopology::valid().

                                                                        {
  bool ok(false);
  if (type == 1) {
    ok = ((izeta > 0) && (izeta <= ftpar_->nZBarrel_) && (iphi > 0) &&
          (iphi <= ftpar_->nPhiBarrel_));
  } else if (type == 2) {
    ok = ((izeta > 0) && (izeta <= ftpar_->nEtaEndcap_) && (iphi > 0) &&
          (iphi <= ftpar_->nPhiEndcap_));
  }
  return ok;
}

int FastTimeDDDConstants::numberEtaZ

(

int

type

)

const

Definition at line 182 of file FastTimeDDDConstants.cc.

References ftpar_, FastTimeParameters::nEtaEndcap_, and FastTimeParameters::nZBarrel_.

Referenced by FastTimeTopology::FastTimeTopology(), and FastTimeGeometry::newCell().

                                                   {
  int numb(0);
  if (type == 1) {
    numb = (ftpar_->nZBarrel_);
  } else if (type == 2) {
    numb = (ftpar_->nEtaEndcap_);
  }
  return numb;
}

int FastTimeDDDConstants::numberPhi

(

int

type

)

const

Definition at line 192 of file FastTimeDDDConstants.cc.

References ftpar_, FastTimeParameters::nPhiBarrel_, and FastTimeParameters::nPhiEndcap_.

Referenced by FastTimeTopology::FastTimeTopology(), and FastTimeGeometry::newCell().

                                                  {
  int numb(0);
  if (type == 1) {
    numb = (ftpar_->nPhiBarrel_);
  } else if (type == 2) {
    numb = (ftpar_->nPhiEndcap_);
  }
  return numb;
}

Member Data Documentation

double FastTimeDDDConstants::dEta_

private

Definition at line 42 of file FastTimeDDDConstants.h.

Referenced by initialize().

double FastTimeDDDConstants::dPhiBarrel_

private

Definition at line 43 of file FastTimeDDDConstants.h.

Referenced by getCorners(), getPosition(), getZPhi(), and initialize().

double FastTimeDDDConstants::dPhiEndcap_

private

Definition at line 43 of file FastTimeDDDConstants.h.

Referenced by getCorners(), getEtaPhi(), getPosition(), and initialize().

double FastTimeDDDConstants::dZBarrel_

private

Definition at line 43 of file FastTimeDDDConstants.h.

Referenced by getCorners(), getPosition(), getZPhi(), and initialize().

double FastTimeDDDConstants::etaMax_

private

Definition at line 42 of file FastTimeDDDConstants.h.

Referenced by initialize().

double FastTimeDDDConstants::etaMin_

private

Definition at line 42 of file FastTimeDDDConstants.h.

Referenced by initialize().

const FastTimeParameters* FastTimeDDDConstants::ftpar_

private

Definition at line 41 of file FastTimeDDDConstants.h.

Referenced by getCells(), getCorners(), getEtaPhi(), getPosition(), getRin(), getRout(), getZHalf(), getZPhi(), getZPos(), initialize(), numberEtaZ(), and numberPhi().

std::vector<double> FastTimeDDDConstants::rLimits_

private

Definition at line 44 of file FastTimeDDDConstants.h.

Referenced by getCorners(), getEtaPhi(), getPosition(), and initialize().