PtAssignmentEngineAux2016 Class Reference

#include <PtAssignmentEngineAux2016.h>

Public Member Functions
int	getCLCT (int clct) const
int	getdEta (int dEta) const
int	getdPhiFromBin (int dPhiBin, int bits, int max=512) const
int	getdTheta (int dTheta) const
float	getEtaFromBin (int etaBin, int bits=5) const
float	getEtaFromEtaInt (int etaInt, int bits=5) const
float	getEtaFromThetaInt (int thetaInt, int bits=5) const
int	getEtaInt (float eta, int bits=5) const
int	getFRLUT (int sector, int station, int chamber) const
const int(*	getModeVariables () const)[6]
int	getNLBdPhi (int dPhi, int bits, int max=512) const
int	getNLBdPhiBin (int dPhi, int bits, int max=512) const

Detailed Description

Definition at line 10 of file PtAssignmentEngineAux2016.h.

Member Function Documentation

int PtAssignmentEngineAux2016::getCLCT

(

int

clct

)

const

Definition at line 270 of file PtAssignmentEngineAux2016.cc.

Referenced by PtAssignmentEngine2016::calculate_address().

                                                     {
  int clct_ = 0;
  int sign_ = 1;

  switch (clct) {
  case 10: clct_ = 0; sign_ = -1; break;
  case  9: clct_ = 1; sign_ =  1; break;
  case  8: clct_ = 1; sign_ = -1; break;
  case  7: clct_ = 2; sign_ =  1; break;
  case  6: clct_ = 2; sign_ = -1; break;
  case  5: clct_ = 3; sign_ =  1; break;
  case  4: clct_ = 3; sign_ = -1; break;
  case  3: clct_ = 3; sign_ =  1; break;
  case  2: clct_ = 3; sign_ = -1; break;
  case  1: clct_ = 3; sign_ = -1; break;  //case  1: clct_ = 3; sign_ =  1; break;
  case  0: clct_ = 3; sign_ = -1; break;
  default: clct_ = 3; sign_ = -1; break;
  }
  return (sign_ * clct_);
}

int PtAssignmentEngineAux2016::getdEta

(

int

dEta

)

const

Definition at line 313 of file PtAssignmentEngineAux2016.cc.

                                                     {
  int dEta_ = 0;

  if (dEta<=-5)
    dEta_ = 0;
  else if (dEta<=-2)
    dEta_ = 1;
  else if (dEta<=-1)
    dEta_ = 2;
  else if (dEta<=0)
    dEta_ = 3;
  else if (dEta<=1)
    dEta_ = 4;
  else if (dEta<=3)
    dEta_ = 5;
  else if (dEta<=6)
    dEta_ = 6;
  else
    dEta_ = 7;
  return (dEta_);
}

int PtAssignmentEngineAux2016::getdPhiFromBin	(	int	dPhiBin,
		int	bits,
		int	max = 512
	)		const

Definition at line 243 of file PtAssignmentEngineAux2016.cc.

References bits, photonValidator_cfi::dPhiBin, dPhiNLBMap_5bit_256Max, dPhiNLBMap_6bit_256Max, dPhiNLBMap_7bit_256Max, dPhiNLBMap_7bit_512Max, and dPhiNLBMap_8bit_512Max.

Referenced by PtAssignmentEngine2016::calculate_pt_xml().

                                                                                  {
  int dPhi_= (1<<bits)-1;

  if (dPhiBin>(1<<bits)-1)
    dPhiBin = (1<<bits)-1;

  if (max==256)
  {
    if (bits == 5)
      dPhi_ = dPhiNLBMap_5bit_256Max[dPhiBin];
    if (bits == 6)
      dPhi_ = dPhiNLBMap_6bit_256Max[dPhiBin];
    if (bits == 7)
      dPhi_ = dPhiNLBMap_7bit_256Max[dPhiBin];
  }

  else if (max==512)
  {
    if (bits == 7)
      dPhi_ = dPhiNLBMap_7bit_512Max[dPhiBin];
    if (bits == 8)
      dPhi_ = dPhiNLBMap_8bit_512Max[dPhiBin];
  }

  return (dPhi_);
}

int PtAssignmentEngineAux2016::getdTheta

(

int

dTheta

)

const

Definition at line 291 of file PtAssignmentEngineAux2016.cc.

Referenced by PtAssignmentEngine2016::calculate_address().

                                                         {
  int dTheta_ = 0;

  if (dTheta<=-3)
    dTheta_ = 0;
  else if (dTheta<=-2)
    dTheta_ = 1;
  else if (dTheta<=-1)
    dTheta_ = 2;
  else if (dTheta<=0)
    dTheta_ = 3;
  else if (dTheta<=1)
    dTheta_ = 4;
  else if (dTheta<=2)
    dTheta_ = 5;
  else if (dTheta<=3)
    dTheta_ = 6;
  else
    dTheta_ = 7;
  return (dTheta_);
}

float PtAssignmentEngineAux2016::getEtaFromBin	(	int	etaBin,
		int	bits = 5
	)		const

Definition at line 368 of file PtAssignmentEngineAux2016.cc.

References getEtaFromEtaInt().

                                                                         {
  // For backward compatibility
  return getEtaFromEtaInt(etaBin, bits);
}

float PtAssignmentEngineAux2016::getEtaFromEtaInt	(	int	etaInt,
		int	bits = 5
	)		const

Definition at line 357 of file PtAssignmentEngineAux2016.cc.

References bits, PVValHelper::eta, and edm::shift.

Referenced by PtAssignmentEngine2016::calculate_pt_xml(), and getEtaFromBin().

                                                                            {
  etaInt = (etaInt > 31) ? 31 : etaInt;  // 5-bit
  etaInt = (etaInt < 0) ? 0 : etaInt;
  bits = (bits > 5) ? 5 : bits;
  int shift = 5 - bits;
  etaInt <<= shift;
  // decode 5-bit etaInt to 0.9<abs(eta)<1.6
  float eta = ((0.5 + etaInt)*1.6)/32 + 0.9;
  return eta;
}

float PtAssignmentEngineAux2016::getEtaFromThetaInt	(	int	thetaInt,
		int	bits = 5
	)		const

Definition at line 347 of file PtAssignmentEngineAux2016.cc.

References PVValHelper::eta, cmsBatch::log, funct::tan(), and theta().

Referenced by PtAssignmentEngine2016::calculate_pt_xml().

                                                                                {
  thetaInt = (thetaInt > 127) ? 127 : thetaInt;  // 7-bit
  thetaInt = (thetaInt < 0) ? 0 : thetaInt;
  float theta = thetaInt;
  theta = (theta*0.2874016 + 8.5)*(3.14159265359/180);  // use this version to reproduce what happened when the pT LUT was made
  //theta = (theta*(45.0-8.5)/128. + 8.5) * M_PI/180.;
  float eta = -std::log(std::tan(theta/2));
  return eta;
}

int PtAssignmentEngineAux2016::getEtaInt	(	float	eta,
		int	bits = 5
	)		const

Definition at line 335 of file PtAssignmentEngineAux2016.cc.

References funct::abs(), bits, PVValHelper::eta, and edm::shift.

Referenced by PtAssignmentEngine2016::calculate_pt_xml().

                                                                  {
  eta = std::abs(eta);
  eta = (eta < 0.9) ? 0.9 : eta;
  bits = (bits > 5) ? 5 : bits;
  // encode 0.9<abs(eta)<1.6 in 5-bit (32 possible values)
  int etaInt = ((eta - 0.9)*32)/(1.6) - 0.5;
  int shift = 5 - bits;
  etaInt >>= shift;
  etaInt = (etaInt > 31) ? 31 : etaInt;
  return etaInt;
}

int PtAssignmentEngineAux2016::getFRLUT	(	int	sector,
		int	station,
		int	chamber
	)		const

Definition at line 381 of file PtAssignmentEngineAux2016.cc.

References bits, and relativeConstraints::chamber.

Referenced by PtAssignmentEngine2016::calculate_address().

                                                                                  {
  int bits = FRLUT[(sector-1)%2][station];
  bool isFront = bits & (1<<chamber);
  return isFront;
}

const int(* PtAssignmentEngineAux2016::getModeVariables

(

)

)[6]

Definition at line 87 of file PtAssignmentEngineAux2016.cc.

References ModeVariables_Scheme3.

Referenced by PtAssignmentEngine2016::calculate_pt_xml().

                                                                  {
  return ModeVariables_Scheme3;
}

int PtAssignmentEngineAux2016::getNLBdPhi	(	int	dPhi,
		int	bits,
		int	max = 512
	)		const

Definition at line 91 of file PtAssignmentEngineAux2016.cc.

References bits, particleFlow_cfi::dPhi, dPhiNLBMap_5bit_256Max, dPhiNLBMap_6bit_256Max, dPhiNLBMap_7bit_256Max, dPhiNLBMap_7bit_512Max, dPhiNLBMap_8bit_512Max, and SiStripPI::max.

                                                                           {
  int dPhi_= max;
  int sign_ = 1;
  if (dPhi<0)
    sign_ = -1;
  dPhi = sign_ * dPhi;

  if (max==256)
  {
    if (bits == 5)
    {
      dPhi_ = dPhiNLBMap_5bit_256Max[(1<<bits)-1];
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_5bit_256Max[edge]<=dPhi && dPhiNLBMap_5bit_256Max[edge+1]>dPhi)
        {
          dPhi_ = dPhiNLBMap_5bit_256Max[edge];
          break;
        }
      }
    }
    if (bits == 6)
    {
      dPhi_ = dPhiNLBMap_6bit_256Max[(1<<bits)-1];
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_6bit_256Max[edge]<=dPhi && dPhiNLBMap_6bit_256Max[edge+1]>dPhi)
        {
          dPhi_ = dPhiNLBMap_6bit_256Max[edge];
          break;
        }
      }
    }
    if (bits == 7)
    {
      dPhi_ = dPhiNLBMap_7bit_256Max[(1<<bits)-1];
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_7bit_256Max[edge]<=dPhi && dPhiNLBMap_7bit_256Max[edge+1]>dPhi)
        {
          dPhi_ = dPhiNLBMap_7bit_256Max[edge];
          break;
        }
      }
    }
  }

  else if (max==512)
  {
    if (bits == 7)
    {
      dPhi_ = dPhiNLBMap_7bit_512Max[(1<<bits)-1];
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_7bit_512Max[edge]<=dPhi && dPhiNLBMap_7bit_512Max[edge+1]>dPhi)
        {
          dPhi_ = dPhiNLBMap_7bit_512Max[edge];
          break;
        }
      }
    }
    if (bits == 8)
    {
      dPhi_ = dPhiNLBMap_8bit_512Max[(1<<bits)-1];
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_8bit_512Max[edge]<=dPhi && dPhiNLBMap_8bit_512Max[edge+1]>dPhi)
        {
          dPhi_ = dPhiNLBMap_8bit_512Max[edge];
          break;
        }
      }
    }
  }

  if (dPhi>=max) dPhi_ = max;
  return (sign_ * dPhi_);
}

int PtAssignmentEngineAux2016::getNLBdPhiBin	(	int	dPhi,
		int	bits,
		int	max = 512
	)		const

Definition at line 170 of file PtAssignmentEngineAux2016.cc.

References bits, particleFlow_cfi::dPhi, dPhiNLBMap_5bit_256Max, dPhiNLBMap_6bit_256Max, dPhiNLBMap_7bit_256Max, dPhiNLBMap_7bit_512Max, and dPhiNLBMap_8bit_512Max.

Referenced by PtAssignmentEngine2016::calculate_address().

                                                                              {
  int dPhiBin_= (1<<bits)-1;
  int sign_ = 1;
  if (dPhi<0)
    sign_ = -1;
  dPhi = sign_ * dPhi;

  if (max==256)
  {
    if (bits == 5)
    {
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_5bit_256Max[edge]<=dPhi && dPhiNLBMap_5bit_256Max[edge+1]>dPhi)
        {
          dPhiBin_ = edge;
          break;
        }
      }
    }
    if (bits == 6)
    {
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_6bit_256Max[edge]<=dPhi && dPhiNLBMap_6bit_256Max[edge+1]>dPhi)
        {
          dPhiBin_ = edge;
          break;
        }
      }
    }
    if (bits == 7)
    {
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_7bit_256Max[edge]<=dPhi && dPhiNLBMap_7bit_256Max[edge+1]>dPhi)
        {
          dPhiBin_ = edge;
          break;
        }
      }
    }
  }

  else if (max==512)
  {
    if (bits == 7)
    {
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_7bit_512Max[edge]<=dPhi && dPhiNLBMap_7bit_512Max[edge+1]>dPhi)
        {
          dPhiBin_ = edge;
          break;
        }
      }
    }
    if (bits == 8)
    {
      for (int edge=0; edge<(1<<bits)-1; edge++)
      {
        if (dPhiNLBMap_8bit_512Max[edge]<=dPhi && dPhiNLBMap_8bit_512Max[edge+1]>dPhi)
        {
          dPhiBin_ = edge;
          break;
        }
      }
    }
  }

  return (dPhiBin_);
}