The Class for ECL Geometry Parameters. More...

#include <ECLGeometryPar.h>

Classes
struct	CrystalGeom_t
	crystal geometry More...

Public Member Functions
	ECLGeometryPar ()
	Constructor.

virtual	~ECLGeometryPar ()
	Destructor.

void	clear ()
	Clears.

void	Print () const
	Print some debug information.

void	read ()
	Gets geometry parameters from PhysicalVolumeStore.

void	Mapping (int cid)
	Mapping theta, phi Id.

int	GetCellID (int ThetaId, int PhiId)
	Get Cell Id.

int	GetCellID ()
	Get Cell Id.

int	GetThetaID ()
	Get Theta Id.

int	GetPhiID ()
	Get Phi Id.

int	ECLVolumeToCellID (const G4VTouchable *)
	Get Cell Id (LEP: new way)

int	TouchableToCellID (const G4VTouchable *)
	The same as above but without sanity checks.

int	TouchableDiodeToCellID (const G4VTouchable *)
	Mapping from G4VTouchable copyNumbers to Crystal CellID.

G4ThreeVector	getCrystalPos (int cid)
	The Position of crystal.

G4ThreeVector	getCrystalVec (int cid)
	The direction of crystal.

ROOT::Math::XYZVector	GetCrystalPos (int cid)
	The Position of crystal.

ROOT::Math::XYZVector	GetCrystalVec (int cid)
	The direction of crystal.

double	time2sensor (int cid, const G4ThreeVector &hit_pos)
	function to calculate flight time to diode sensor

Static Public Member Functions
static ECLGeometryPar *	Instance ()
	Static method to get a reference to the ECLGeometryPar instance.

Private Member Functions
void	InitCrystal (int cid)
	initialise the crystal

Private Attributes
G4Transform3D *	m_ECLForwardGlobalT = nullptr
	Global transformations for the forward part.

G4Transform3D *	m_ECLBackwardGlobalT = nullptr
	Global transformations for the backward part.

G4Transform3D *	m_ECLBarrelGlobalT = nullptr
	Global transformations for the barrel part.

std::vector< CrystalGeom_t >	m_crystals
	the crystals

CrystalGeom_t	m_current_crystal
	the current crystal

int	mPar_cellID
	The Cell ID information.

int	mPar_thetaID
	The Theta ID information.

int	mPar_phiID
	The Phi ID information.

int	m_ini_cid
	initial crystal ID

Static Private Attributes
static ECLGeometryPar *	m_B4ECLGeometryParDB = 0
	Pointer that saves the instance of this class.

Detailed Description

The Class for ECL Geometry Parameters.

This class provides ECL geometry parameters for simulation, reconstruction and so on. These parameters are got from geometry description

Definition at line 44 of file ECLGeometryPar.h.

Constructor & Destructor Documentation

◆ ECLGeometryPar()

ECLGeometryPar ( )

Constructor.

Definition at line 174 of file ECLGeometryPar.cc.

{
  clear();
  // delay crystal positions fetching to a moment when it actually needed and geometry is already built
  //  read();
}

◆ ~ECLGeometryPar()

~ECLGeometryPar ( )

virtual

Destructor.

Definition at line 181 of file ECLGeometryPar.cc.

{
  if (m_B4ECLGeometryParDB) {
    delete m_B4ECLGeometryParDB;
    B2DEBUG(150, "m_B4ECLGeometryParDB deleted ");
  }
  if (m_ECLForwardGlobalT) delete m_ECLForwardGlobalT;
  if (m_ECLBarrelGlobalT) delete m_ECLBarrelGlobalT;
  if (m_ECLBackwardGlobalT) delete m_ECLBackwardGlobalT;
}

Member Function Documentation

◆ clear()

void clear ( )

Clears.

Definition at line 192 of file ECLGeometryPar.cc.

{
  m_ini_cid = -1;
  mPar_cellID = 0;
  mPar_thetaID = 0;
  mPar_phiID = 0;
}

◆ ECLVolumeToCellID()

int ECLVolumeToCellID ( const G4VTouchable * touch )

Get Cell Id (LEP: new way)

Mapping from G4VTouchable copyNumbers to Crystal CellID

Definition at line 478 of file ECLGeometryPar.cc.

{
  int depth = touch->GetHistoryDepth();
  if ((depth != 3) && (depth != 5)) {
    B2WARNING("ECLGeometryPar::ECLVolumeToCellID: History depth = " << depth << " is out of range: should be 3 or 5.");
    return -1;
  }
  const G4String& vname = touch->GetVolume()->GetName();
  std::size_t pos0 = vname.find("lv_forward_crystal_");
  std::size_t pos1 = vname.find("lv_barrel_crystal_");
  std::size_t pos2 = vname.find("lv_backward_crystal_");
  if (pos0 == string::npos && pos1 == string::npos && pos2 == string::npos) {
    B2WARNING("ECLGeometryPar::ECLVolumeToCellID: Volume name does not match pattern. NAME=" << vname);
    return -1;
  }
  return TouchableToCellID(touch);
}

◆ GetCellID() [1/2]

int GetCellID ( )

inline

Get Cell Id.

Definition at line 75 of file ECLGeometryPar.h.

75{return mPar_cellID;};

◆ GetCellID() [2/2]

int GetCellID	(	int	ThetaId,
		int	PhiId
	)

Get Cell Id.

Definition at line 399 of file ECLGeometryPar.cc.

{
  return Mapping_t::CellID(ThetaId, PhiId);
}

◆ getCrystalPos()

G4ThreeVector getCrystalPos ( int cid )

inline

The Position of crystal.

Definition at line 87 of file ECLGeometryPar.h.

      {
        if (cid != m_ini_cid) InitCrystal(cid);
        return m_current_crystal.pos;
      }

◆ GetCrystalPos()

ROOT::Math::XYZVector GetCrystalPos ( int cid )

inline

The Position of crystal.

Definition at line 101 of file ECLGeometryPar.h.

      {
        if (cid != m_ini_cid) InitCrystal(cid);
        const G4ThreeVector& t = m_current_crystal.pos;
        return ROOT::Math::XYZVector(t.x(), t.y(), t.z());
      }

◆ getCrystalVec()

G4ThreeVector getCrystalVec ( int cid )

inline

The direction of crystal.

Definition at line 94 of file ECLGeometryPar.h.

      {
        if (cid != m_ini_cid) InitCrystal(cid);
        return m_current_crystal.dir;
      }

◆ GetCrystalVec()

ROOT::Math::XYZVector GetCrystalVec ( int cid )

inline

The direction of crystal.

Definition at line 109 of file ECLGeometryPar.h.

      {
        if (cid != m_ini_cid) InitCrystal(cid);
        const G4ThreeVector& t = m_current_crystal.dir;
        return ROOT::Math::XYZVector(t.x(), t.y(), t.z());
      }

◆ GetPhiID()

int GetPhiID ( )

inline

Get Phi Id.

Definition at line 79 of file ECLGeometryPar.h.

79{return mPar_phiID;};

◆ GetThetaID()

int GetThetaID ( )

inline

Get Theta Id.

Definition at line 77 of file ECLGeometryPar.h.

77{return mPar_thetaID;};

◆ InitCrystal()

void InitCrystal ( int cid )

private

initialise the crystal

Definition at line 362 of file ECLGeometryPar.cc.

{
  if (m_crystals.size() == 0) read();
  int thetaid, phiid, nreplica, indx;
  Mapping_t::Mapping(cid, thetaid, phiid, nreplica, indx);
  //  cout<<cid<<" "<<thetaid<<" "<<phiid<<" "<<nreplica<<" "<<indx<<endl;
  const CrystalGeom_t& t = m_crystals[indx];
  double s, c;
  if (indx > 131)
    sincos<72>(ss72, nreplica, s, c);
  else
    sincos<16>(ss16, nreplica, s, c);
 
  G4Transform3D* T;
  if (ECLElementNumbers::isForward(cid + 1)) {
    T = m_ECLForwardGlobalT;
  } else if (ECLElementNumbers::isBarrel(cid + 1)) {
    T = m_ECLBarrelGlobalT;
  } else {
    T = m_ECLBackwardGlobalT;
    c = -c;
  }
  double xp = c * t.pos.x() - s * t.pos.y();
  double yp = s * t.pos.x() + c * t.pos.y();
  //  m_current_crystal.pos.set(xp, yp, t.pos.z());
 
  double xv = c * t.dir.x() - s * t.dir.y();
  double yv = s * t.dir.x() + c * t.dir.y();
  //  m_current_crystal.dir.set(xv, yv, t.dir.z());
 
  m_current_crystal.pos = (1 / CLHEP::cm) * (*T * G4Point3D(xp, yp, t.pos.z()));
  m_current_crystal.dir = *T * G4Vector3D(xv, yv, t.dir.z());
 
  //  cout<<t.pos<<" "<<t.dir<<" "<<m_current_crystal.pos<<" "<<m_current_crystal.dir<<endl;
  m_ini_cid = cid;
}

◆ Instance()

ECLGeometryPar * Instance ( )

static

Static method to get a reference to the ECLGeometryPar instance.

Returns: A reference to an instance of this class.

Definition at line 168 of file ECLGeometryPar.cc.

{
  if (!m_B4ECLGeometryParDB) m_B4ECLGeometryParDB = new ECLGeometryPar();
  return m_B4ECLGeometryParDB;
}

◆ Mapping()

void Mapping ( int cid )

Mapping theta, phi Id.

Definition at line 404 of file ECLGeometryPar.cc.

{
  mPar_cellID = cid;
  Mapping_t::Mapping(mPar_cellID, mPar_thetaID, mPar_phiID);
}

◆ read()

void read ( )

Gets geometry parameters from PhysicalVolumeStore.

Definition at line 253 of file ECLGeometryPar.cc.

{
  m_crystals.clear();
  m_crystals.reserve(46 * 2 + 132); // 10752 bytes
 
  // Endcap sectors are rotated since behaviour of G4Replica class. It
  // requires a physical volume during phi replication to be symmetric
  // around phi=0. So we need to rotate it by -((2*pi)/16)/2 angle at the
  // geometry description are return back here.
 
  G4Point3D p0(0, 0, 0); G4Vector3D v0(0, 0, 1);
 
  {
    m_ECLForwardGlobalT = new G4Transform3D(getT("ECLForwardPhysical"));
    G4Transform3D T =  getT("ECLForwardCrystalSectorPhysical_1");
    G4String tnamef("ECLForwardWrappedCrystal_Physical_");
    for (int i = 0; i < 72; i++) {
      G4String vname(tnamef); vname += to_string(i);
      G4Transform3D cT = T * getT(vname);
      CrystalGeom_t c = {cT * p0, cT * v0};
      m_crystals.push_back(c);
      //      G4cout << i << " " << c.pos << " " << c.dir << G4endl;
    }
  }
 
  {
    m_ECLBackwardGlobalT = new G4Transform3D(getT("ECLBackwardPhysical"));
    G4Transform3D T = getT("ECLBackwardCrystalSectorPhysical_0");
    G4String tnamef("ECLBackwardWrappedCrystal_Physical_");
    for (int i = 0; i < 60; i++) {
      G4String vname(tnamef); vname += to_string(i);
      G4Transform3D cT = T * getT(vname);
      CrystalGeom_t c = {cT * p0, cT * v0};
      m_crystals.push_back(c);
      //      G4cout << i << " " << c.pos << " " << c.dir << G4endl;
    }
  }
 
  {
    // get barrel sector (between two septums) transformation
 
    // extract global barrel movements assuming we know transformation of the sector in axial symmetric situation
    m_ECLBarrelGlobalT = new G4Transform3D(getT("ECLBarrelSectorPhysical_0") * G4RotateZ3D(M_PI / 2));
    G4Transform3D Ts = G4RotateZ3D(-M_PI / 2);
    // since barrel sector is symmetric around phi=0 we need to
    // translate crystal with negative phi back to positive rotating
    // crystal position by (2*M_PI/72) angle
    G4Transform3D Ts1 = G4RotateZ3D(M_PI / 36) * Ts;
 
    G4String tname("ECLBarrelWrappedCrystal_Physical_");
    for (int i = 0; i < 2 * 46; i++) {
      G4String vname(tname); vname += to_string(i);
      G4Transform3D cT = ((i % 2) ? Ts1 : Ts) * getT(vname);
      CrystalGeom_t c = {cT * p0, cT * v0};
      m_crystals.push_back(c);
      //      G4cout << i << " " << c.pos << " " <<c.dir<<G4endl;
    }
  }
 
  B2DEBUG(150, "ECLGeometryPar::read() initialized with " << m_crystals.size() << " crystals.");
}

◆ time2sensor()

double time2sensor	(	int	cid,
		const G4ThreeVector &	hit_pos
	)

function to calculate flight time to diode sensor

Definition at line 496 of file ECLGeometryPar.cc.

{
  if (cid != m_ini_cid) InitCrystal(cid);
  double z = 15. - (hit_pos - m_current_crystal.pos) * m_current_crystal.dir; // position along the vector of crystal axis
  double dt = 6.05 + z * (0.0749 - z * 0.00112); // flight time to diode sensor in nanoseconds
  return dt;
}

◆ TouchableDiodeToCellID()

int TouchableDiodeToCellID ( const G4VTouchable * touch )

Mapping from G4VTouchable copyNumbers to Crystal CellID.

Definition at line 410 of file ECLGeometryPar.cc.

{
  return TouchableToCellID(touch);
}

◆ TouchableToCellID()

int TouchableToCellID ( const G4VTouchable * touch )

The same as above but without sanity checks.

Definition at line 415 of file ECLGeometryPar.cc.

{
  //  touch->GetCopyNumber() is a virtual call, avoid it by using
  //  directly G4NavigationHistory so we will have only one virtual
  //  call instead of three here
  const G4NavigationHistory* h = touch->GetHistory();
  int hd = h->GetDepth();
  int i1 = h->GetReplicaNo(hd - 1); // index of each volume is set at physical volume creation
  int i2 = h->GetReplicaNo(hd - 2); // go up in volume hierarchy
 
  int ThetaId = Mapping_t::Indx2ThetaId(i1);
  int NCryst  = Mapping_t::ThetaId2NCry(ThetaId); // the number of crystals in a sector at given ThetaId
  int Offset  = Mapping_t::Offset(ThetaId);
 
  int ik = i1 - Offset;
  int PhiId;
  if (NCryst == 2) {
    PhiId = (i2 - (ik % 2)) * NCryst + ik;
    if (PhiId < 0) PhiId += 144;
  } else {
    int i3 = h->GetReplicaNo(hd - 3); // go up in volume hierarchy
    if (ThetaId < 13)
      PhiId = (i2 + 2 * i3) * NCryst + ik;
    else {
      // int tt[] = {3,2,1,0,7,6,5,4};
      // if(i3r!=tt[i3]){
      //  cout<<i3<<" "<<i3<<" "<<tt[i3]<<endl;
      //  exit(0);
      // }
      int i3r = (3 - (i3 % 4)) + 4 * (i3 / 4);
      PhiId = (2 * i3r + (1 - i2)) * NCryst + ik;
    }
  }
 
  int cellID = Offset * 16 + PhiId;
  //    cout<<"ECLGeometryPar::TouchableToCellID "<<h->GetVolume(hd-1)->GetName()<<" "<<i1<<" "<<i2<<" "<<h->GetReplicaNo(hd - 3)<<" "<<ThetaId<<" "<<NCryst<<" "<<Offset<<" "<<PhiId<<endl;
 
  // test of the position and direction of crystal
  if (0) {
    G4AffineTransform t = h->GetTopTransform();
    G4ThreeVector o(0, 0, 0), n(0, 0, 1);
    G4ThreeVector ro = t.Inverse().TransformPoint(o);
    G4ThreeVector rn = t.Inverse().TransformAxis(n);
 
    InitCrystal(cellID);
    ro *= 1 / CLHEP::cm;
 
    G4ThreeVector dr = m_current_crystal.pos - ro, dn = m_current_crystal.dir - rn;
    if (dr.mag() > 1e-10 || dn.mag() > 1e-10) {
      cout << "Mismatch " << h->GetVolume(hd - 1)->GetName() << " " << cellID << " " << ThetaId << " " << PhiId << " " << NCryst << " "
           << hd << " " << i2 << " " << i1 << " " << m_current_crystal.pos << " " << ro << " " << rn << " " << dr << " " << dn << endl;
 
      for (int i = 0; i < 144; i++) {
        int ci = Mapping_t::CellID(ThetaId, i);
        InitCrystal(ci);
        dr = m_current_crystal.pos - ro;
        if (dr.mag() < 1e-10) cout << "best PhiId = " << i << endl;
      }
    }
  }
  return cellID;
}

Member Data Documentation

◆ m_B4ECLGeometryParDB

ECLGeometryPar * m_B4ECLGeometryParDB = 0

staticprivate

Pointer that saves the instance of this class.

Definition at line 145 of file ECLGeometryPar.h.

◆ m_crystals

std::vector<CrystalGeom_t> m_crystals

private

the crystals

Definition at line 133 of file ECLGeometryPar.h.

◆ m_current_crystal

CrystalGeom_t m_current_crystal

private

the current crystal

Definition at line 135 of file ECLGeometryPar.h.

◆ m_ECLBackwardGlobalT

G4Transform3D* m_ECLBackwardGlobalT = nullptr

private

Global transformations for the backward part.

Definition at line 122 of file ECLGeometryPar.h.

◆ m_ECLBarrelGlobalT

G4Transform3D* m_ECLBarrelGlobalT = nullptr

private

Global transformations for the barrel part.

Definition at line 124 of file ECLGeometryPar.h.

◆ m_ECLForwardGlobalT

G4Transform3D* m_ECLForwardGlobalT = nullptr

private

Global transformations for the forward part.

Definition at line 120 of file ECLGeometryPar.h.

◆ m_ini_cid

int m_ini_cid

private

initial crystal ID

Definition at line 143 of file ECLGeometryPar.h.

◆ mPar_cellID

int mPar_cellID

private

The Cell ID information.

Definition at line 137 of file ECLGeometryPar.h.

◆ mPar_phiID

int mPar_phiID

private

The Phi ID information.

Definition at line 141 of file ECLGeometryPar.h.

◆ mPar_thetaID

int mPar_thetaID

private

The Theta ID information.

Definition at line 139 of file ECLGeometryPar.h.

The documentation for this class was generated from the following files:

ecl/geometry/include/ECLGeometryPar.h
ecl/geometry/src/ECLGeometryPar.cc

Classes

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ECLGeometryPar()

◆ ~ECLGeometryPar()

Member Function Documentation

◆ clear()

◆ ECLVolumeToCellID()

◆ GetCellID() [1/2]

◆ GetCellID() [2/2]

◆ getCrystalPos()

◆ GetCrystalPos()

◆ getCrystalVec()

◆ GetCrystalVec()

◆ GetPhiID()

◆ GetThetaID()

◆ InitCrystal()

◆ Instance()

◆ Mapping()

◆ read()

◆ time2sensor()

◆ TouchableDiodeToCellID()

◆ TouchableToCellID()

Member Data Documentation

◆ m_B4ECLGeometryParDB

◆ m_crystals

◆ m_current_crystal

◆ m_ECLBackwardGlobalT

◆ m_ECLBarrelGlobalT

◆ m_ECLForwardGlobalT

◆ m_ini_cid

◆ mPar_cellID

◆ mPar_phiID

◆ mPar_thetaID