HcalTB02HcalNumberingScheme Class Reference

#include <SimG4CMS/HcalTestBeam/interface/HcalTB02HcalNumberingScheme.h>

Inheritance diagram for HcalTB02HcalNumberingScheme:

Public Member Functions
int	getetaID (int sID) const
int	getetaScaleF () const
int	getlayerID (int sID) const
int	getphiID (int sID) const
int	getphiScaleF () const
int	getUnitID (const G4Step *aStep) const override
	HcalTB02HcalNumberingScheme ()
	~HcalTB02HcalNumberingScheme () override
Public Member Functions inherited from HcalTB02NumberingScheme
	HcalTB02NumberingScheme ()
virtual	~HcalTB02NumberingScheme ()

Private Attributes
int	etaScale
int	phiScale

Detailed Description

Description: Numbering scheme for hadron calorimeter in 2002 test beam

Usage: Sets up unique identifier for HB towers in 2002 test beam

Definition at line 24 of file HcalTB02HcalNumberingScheme.h.

Constructor & Destructor Documentation

HcalTB02HcalNumberingScheme()

HcalTB02HcalNumberingScheme::HcalTB02HcalNumberingScheme

(

)

Definition at line 25 of file HcalTB02HcalNumberingScheme.cc.

    : HcalTB02NumberingScheme(), phiScale(1000000), etaScale(10000) {
  edm::LogVerbatim("HcalTBSim") << "Creating HcalTB02HcalNumberingScheme";
}

~HcalTB02HcalNumberingScheme()

HcalTB02HcalNumberingScheme::~HcalTB02HcalNumberingScheme ( )

override

Definition at line 30 of file HcalTB02HcalNumberingScheme.cc.

                                                          {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalTBSim") << "Deleting HcalTB02HcalNumberingScheme";
#endif
}

Member Function Documentation

getetaID()

int HcalTB02HcalNumberingScheme::getetaID

(

int

sID

)

const

Definition at line 110 of file HcalTB02HcalNumberingScheme.cc.

                                                       {
  sID = abs(sID);
  int aux = sID - int(float(sID) / float(phiScale)) * phiScale;
  int etaID = int(float(aux) / float(etaScale));
 
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalTBSim") << "HcalTB02HcalNumberingScheme:: scintID " << sID << " eta = " << etaID;
#endif
  return etaID;
}

References funct::abs(), printConversionInfo::aux, etaScale, createfilelist::int, and phiScale.

getetaScaleF()

int HcalTB02HcalNumberingScheme::getetaScaleF ( ) const

inline

Definition at line 31 of file HcalTB02HcalNumberingScheme.h.

{ return etaScale; }

References etaScale.

getlayerID()

int HcalTB02HcalNumberingScheme::getlayerID

(

int

sID

)

const

Definition at line 81 of file HcalTB02HcalNumberingScheme.cc.

                                                         {
  sID = abs(sID);
  int layerID = sID;
  if ((layerID != 17) && (layerID != 18))
    layerID = sID - int(float(sID) / float(etaScale)) * etaScale;
 
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalTBSim") << "HcalTB02HcalNumberingScheme:: scintID " << sID << " layer = " << layerID;
#endif
  return layerID;
}

References funct::abs(), etaScale, and createfilelist::int.

getphiID()

int HcalTB02HcalNumberingScheme::getphiID

(

int

sID

)

const

Definition at line 93 of file HcalTB02HcalNumberingScheme.cc.

                                                       {
  float IDsign = 1.;
  if (sID < 0)
    IDsign = -1;
  sID = abs(sID);
  int phiID = int(float(sID) / float(phiScale));
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalTBSim") << "HcalTB02HcalNumberingScheme:: scintID " << sID << " phi = " << phiID;
#endif
  if (IDsign > 0) {
    phiID += 4;
  } else {
    phiID = abs(phiID - 3);
  }
  return phiID;
}

References funct::abs(), createfilelist::int, and phiScale.

getphiScaleF()

int HcalTB02HcalNumberingScheme::getphiScaleF ( ) const

inline

Definition at line 30 of file HcalTB02HcalNumberingScheme.h.

{ return phiScale; }

References phiScale.

getUnitID()

int HcalTB02HcalNumberingScheme::getUnitID ( const G4Step *  aStep ) const

overridevirtual

Implements HcalTB02NumberingScheme.

Definition at line 40 of file HcalTB02HcalNumberingScheme.cc.

                                                                    {
  int scintID = 0;
 
  G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
  const G4ThreeVector& hitPoint = preStepPoint->GetPosition();
  float hx = hitPoint.x();
  float hy = hitPoint.y();
  float hz = hitPoint.z();
  float hr = std::sqrt(pow(hx, 2) + pow(hy, 2));
 
  // Check if hit happened in first HO layer or second.
 
  if ((hr > 3. * m) && (hr < 3.830 * m))
    return scintID = 17;
  if (hr > 3.830 * m)
    return scintID = 18;
 
  // Compute the scintID in the HB.
 
  float hR = hitPoint.mag();  //sqrt( pow(hx,2)+pow(hy,2)+pow(hz,2) );
  float htheta = (hR == 0. ? 0. : acos(std::max(std::min(hz / hR, float(1.)), float(-1.))));
  float hsintheta = sin(htheta);
  float hphi = (hR * hsintheta == 0. ? 0. : acos(std::max(std::min(hx / (hR * hsintheta), float(1.)), float(-1.))));
  float heta = (fabs(hsintheta) == 1. ? 0. : -log(fabs(tan(htheta / 2.))));
  int eta = int(heta / 0.087);
  int phi = int(hphi / (5. * degree));
 
  G4VPhysicalVolume* thePV = preStepPoint->GetPhysicalVolume();
  int ilayer = ((thePV->GetCopyNo()) / 10) % 100;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalTBSim") << "HcalTB02HcalNumberingScheme:: Layer " << thePV->GetName()
                                << " found at phi = " << phi << " eta = " << eta << " lay = " << thePV->GetCopyNo()
                                << " " << ilayer;
#endif
  scintID = phiScale * phi + etaScale * eta + ilayer;
  if (hy < 0.)
    scintID = -scintID;
 
  return scintID;
}

References PVValHelper::eta, etaScale, createfilelist::int, dqm-mbProfile::log, visualization-live-secondInstance_cfg::m, SiStripPI::max, min(), phi, phiScale, funct::pow(), funct::sin(), mathSSE::sqrt(), and funct::tan().

Member Data Documentation

etaScale

int HcalTB02HcalNumberingScheme::etaScale

private

Definition at line 39 of file HcalTB02HcalNumberingScheme.h.

Referenced by getetaID(), getetaScaleF(), getlayerID(), and getUnitID().

phiScale

int HcalTB02HcalNumberingScheme::phiScale

private

Definition at line 38 of file HcalTB02HcalNumberingScheme.h.

Referenced by getetaID(), getphiID(), getphiScaleF(), and getUnitID().