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 20 of file FastTimeDDDConstants.h.

Constructor & Destructor Documentation

FastTimeDDDConstants::FastTimeDDDConstants

(

const FastTimeParameters *

ftp

)

Definition at line 11 of file FastTimeDDDConstants.cc.

References gather_cfg::cout, and initialize().

                                                                       : ftpar_(ft) {

#ifdef EDM_ML_DEBUG
  std::cout << "FastTimeDDDConstants::FastTimeDDDConstants ( const FastTimeParameters* ft ) constructor\n";
#endif
  initialize();

}

FastTimeDDDConstants::~FastTimeDDDConstants

(

)

Definition at line 20 of file FastTimeDDDConstants.cc.

References gather_cfg::cout.

                                            { 
#ifdef EDM_ML_DEBUG
  std::cout << "FastTimeDDDConstants:destructed!!!" << std::endl;
#endif
}

Member Function Documentation

int FastTimeDDDConstants::getCells

(

int

type

)

const

Definition at line 115 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 80 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 gather_cfg::cout, dPhiEndcap_, ftpar_, createfilelist::int, gen::k, FastTimeParameters::nEtaEndcap_, FastTimeParameters::nPhiEndcap_, and rLimits_.

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

                                                                           {

  if (phi < 0) phi += CLHEP::twopi;
  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
  std::cout << "FastTimeDDDConstants:Endcap r|phi " << r << " " 
        << phi/CLHEP::deg << " ir|iphi " << ir << " " << iphi << std::endl;
#endif
  return std::pair<int,int>(ir,iphi);
}

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

Definition at line 58 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 125 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 135 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 145 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 gather_cfg::cout, dPhiBarrel_, dZBarrel_, ftpar_, createfilelist::int, FastTimeParameters::nPhiBarrel_, and FastTimeParameters::nZBarrel_.

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

                                                                         {

  if (phi < 0) phi += CLHEP::twopi;
  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
  std::cout << "FastTimeDDDConstants:Barrel z|phi " << z << " " 
        << phi/CLHEP::deg << " iz|iphi " << iz << " " << iphi << std::endl;
#endif
  return std::pair<int,int>(iz,iphi);
}

double FastTimeDDDConstants::getZPos

(

int

type

)

const

Definition at line 155 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 197 of file FastTimeDDDConstants.cc.

References gather_cfg::cout, 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_, 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
  std::cout << "Theta range " << thmin/CLHEP::deg << ":" << thmax/CLHEP::deg
        << " Eta range " << etaMin_ << ":" << etaMax_ << ":" << dEta_ 
        << std::endl;
#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_ = CLHEP::twopi/ftpar_->nPhiBarrel_;
  dPhiEndcap_ = CLHEP::twopi/ftpar_->nPhiEndcap_;
#ifdef EDM_ML_DEBUG
  std::cout << "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" << std::endl;
  for (unsigned int k=0; k<rLimits_.size(); ++k) {
    std::cout << "[" << k << "] " << rLimits_[k] << " ";
    if (k%8 == 7) std::cout << std::endl;
  }
  if ((rLimits_.size()-1)%8 != 7) std::cout << std::endl;
#endif
}

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

Definition at line 165 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 177 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 187 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 46 of file FastTimeDDDConstants.h.

Referenced by initialize().

double FastTimeDDDConstants::dPhiBarrel_

private

Definition at line 47 of file FastTimeDDDConstants.h.

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

double FastTimeDDDConstants::dPhiEndcap_

private

Definition at line 47 of file FastTimeDDDConstants.h.

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

double FastTimeDDDConstants::dZBarrel_

private

Definition at line 47 of file FastTimeDDDConstants.h.

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

double FastTimeDDDConstants::etaMax_

private

Definition at line 46 of file FastTimeDDDConstants.h.

Referenced by initialize().

double FastTimeDDDConstants::etaMin_

private

Definition at line 46 of file FastTimeDDDConstants.h.

Referenced by initialize().

const FastTimeParameters* FastTimeDDDConstants::ftpar_

private

Definition at line 45 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 48 of file FastTimeDDDConstants.h.

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