Modules
	arich data objects

	arich modules

Classes
class	ARICHCalibrationChecker
	ARICH calibration checker. More...

class	ARICHChannelMaskMaker
	ARICH channel mask calibration algorithm. More...

class	ARICHDatabaseImporter
	ARICH database importer. More...

class	ARICHAerogelInfo
	The Class for ARICH Aerogel Parameters. More...

class	ARICHAerogelMap
	Map of the Aerogel placement. More...

class	ARICHAerogelRayleighScatteringFit
	The Class for ARICH aerogel Rayleigh scattering fit parameters. More...

class	ARICHAeroTilesAlignment
	Alignment parameters for aerogel tiles. More...

class	ARICHAeroTilesInfo
	The Class for ARICH aerogel tiles properties for even reconstruction. More...

class	ARICHAsicChannelMask
	ARICH Channel Mask. More...

class	ARICHAsicInfo
	Tested ASIC chips. More...

class	ARICHB2LinkInfo
	List of Belle2Link boards. More...

class	ARICHBadChannels
	ARICH bad channels. More...

class	ARICHBiasCablesMapping
	The Class for ARICH mapping of bias cables to modules. More...

class	ARICHBiasChannelsMapping
	This class provides ARICH mapping of HAPD modules to bias power supply channels. More...

class	ARICHBiasCrateCableMapping
	This class provides ARICH bias voltages at gain 40. More...

class	ARICHBiasVoltages
	The Class for ARICH bias voltages at gain 40. More...

class	ARICHCableInfo
	List of cables. More...

class	ARICHChannelMapping
	The Class for ARICH HAPD channel mapping. More...

class	ARICHChannelMask
	The Class for ARICH HAPD channel mask. More...

class	ARICHComponentTest
	Test of the ARICH component. More...

class	ARICHCopperMapping
	The Class for ARICH mapping of merger to copper boards. More...

class	ARICHFEBoardInfo
	Front End Electronics Boards. More...

class	ARICHFebTest
	ARICH FEB test data. More...

class	ARICHFirmware
	ARICH firmware. More...

class	ARICHGeoAerogelPlane
	Geometry parameters of HAPD. More...

class	ARICHGeoBase
	Base class for geometry parameters. More...

class	ARICHGeoCablesEnvelope
	Geometry parameters of cable envelope. More...

class	ARICHGeoCooling
	Geometry parameters of cooling system. More...

class	ARICHGeoDetectorPlane
	Geometry parameters of ARICH photon detector plane. More...

class	ARICHGeoFEBCooling
	Geometry parameters of Cooling System - version2 (v2). More...

class	ARICHGeoGlobalDisplacement
	Global displacement parameters for ARICH. More...

class	ARICHGeoHAPD
	Geometry parameters of HAPD. More...

class	ARICHGeoMasterVolume
	Geometry parameters of ARICH Master volume (envelope) More...

class	ARICHGeoMerger
	Geometry parameters of Merger PCB. More...

class	ARICHGeoMergerCooling
	Geometry parameters of Merger Cooling System - version2 (v2). More...

class	ARICHGeometryConfig
	The Class for ARICH Geometry Parameters. More...

class	ARICHGeoMirrorDisplacement
	Mirror displacement parameters for ARICH. More...

class	ARICHGeoMirrors
	Geometry parameters of HAPD. More...

class	ARICHGeoSupport
	Geometry parameters of ARICH support structures and neutron shield. More...

class	ARICHGlobalAlignment
	Geometry parameters of ARICH Master volume (envelope). More...

class	ARICHHapdChipInfo
	Contains manufacturer data of one of the 4 photo sensors chips. More...

class	ARICHHapdInfo
	Contains manufacturer data of the photo sensor - HAPD. More...

class	ARICHHapdMask
	ARICH HAPD Sensor Mask: map of operational/nonoperational HAPDs : one bit per sensor. More...

class	ARICHHapdQA
	Contains manufacturer data of the photo sensor - HAPD. More...

class	ARICHHapdQE
	Contains manufacturer data of the photo sensor - HAPD. More...

class	ARICHHighVoltageConfig
	Configuration parameters of the High Voltage Channel. More...

class	ARICHHvCablesMapping
	The Class for ARICH mapping of bias cables to modules. More...

class	ARICHHvChannelsMapping
	The Class for ARICH mapping of bias power supply channels to modules. More...

class	ARICHHvCrateCableMapping
	The Class for ARICH bias voltages at gain 40. More...

class	ARICHMagnetTest
	Contains fractions of dead times measured in a magnet. More...

class	ARICHMergerInfo
	List of MergerBoards. More...

class	ARICHMergerMap
	Mapping of the merger board to the detector. More...

class	ARICHMergerMapping
	The Class for ARICH mapping of modules to merger boards. More...

class	ARICHMirrorAlignment
	Mirror alignment parameters for ARICH. More...

class	ARICHMirrorInfo
	List of MergerBoards. More...

class	ARICHMirrorMap
	Mapping of the mirrors on the detector. More...

class	ARICHModulesInfo
	The Class for information on HAPD modules installed in ARICH. More...

class	ARICHModuleTest
	ARICH module test data. More...

class	ARICHPositionElement
	Position element for ARICH. More...

class	ARICHReconstructionPar
	The Class for ARICH reconstruction parameters. More...

class	ARICHSensorModuleInfo
	Sensor-module information. More...

class	ARICHSensorModuleMap
	Mapping of the sensor-board connections to the detector. More...

class	ARICHSimulationPar
	The Class for ARICH simulation parameters. More...

struct	tessellatedSolidStr
	Structure which holds apexes of the tessellation volumes. More...

class	ARICHTracking
	Beamtest ARICH Geometry Tracking Class. More...

class	ARICHBtestGeometryPar
	The Class for ARICH Beamtest Geometry Parameters. More...

class	ARICHGeometryPar
	The Class for ARICH Geometry Parameters. More...

class	ARICHAerogelHist
	Base class for geometry parameters. More...

class	ARICHChannelHist
	ARICH histogram with HAPD plane 3 options for bin segmentation are available type 0 - one bin per HAPD, type 1 - one bin per channel, type 2 - one bin per APD chip HAPD numbering is 1-420 (software numbering), channel is ASIC channel, APD is ASIC channel / 36. More...

Functions
double	getDetectorZPosition () const
	get z position of detector plane (starting z of HAPDs)

int	doBeamBackgroundStudy () const
	returns 1 if beam background study (to add additional sensitive modules, detect neutrons, ...), 0 else

double	getColEffi () const
	get HAPD collection efficiency.

double	getSensitiveSurfaceSize () const
	get size of detector sensitive surface (size of two chips + gap between)

int	getNMCopies () const
	get the total number of HAPD modules

ROOT::Math::XYVector	getChipLocPos (int chipID)
	get center position of chipID-th chip of detector module (in detector module local coordinates)

bool	isInit ()
	returns true if parametrisation is initialized, false otherwise

bool	isSimple ()
	returns true if parametrisation is initialized with simple geometry (beamtest)

double	getAerogelTransmissionLength (int ilayer)
	get transmission length of "ilayer" aerogel layer

double	getAerogelRefIndex (int ilayer)
	get refractive index of "ilayer" aerogel layer

double	getAerogelThickness (int ilayer)
	get thickness of "ilayer" aerogel layer

double	getAerogelZPosition (int ilayer)
	get z position of "ilayer" aerogel layer

int	getNMirrors ()
	get number of mirror plates

double	getDetectorWindowThickness ()
	get thickness of detector module window

double	getDetectorWindowRefIndex ()
	get refractive index of detector module window

int	getNumberOfAerogelRadiators ()
	get number of aerogel layers

double	getDetectorPadSize ()
	get detectors module pad size

int	getDetectorXPadNumber ()
	get number of pads of detector module (in one direction)

double	getMirrorsStartAngle ()
	get the phi angle of first corner of mirror polygon

double	getMirrorsZPosition ()
	get the z position of mirror plates (starting z)

double	getDetectorZPosition ()
	get z position of detector plane (starting z)

double	getQEScaling ()
	get QE scaling factor for photons internally reflected in HAPD window

double	getWindowAbsorbtion ()
	get absorbtion probability for photons internally reflected in HAPD window

double	getChipNegativeCrosstalk ()
	get factor to suppress p.e. detection efficiency due to negative polarity crosstalk among chip channels

int	getNrTiles ()
	get number of aerogel tiles in radial direction (number of rings)

int	getNphiTiles (int nRing)
	get number of aerogel tiles in each radial ring

double	getTileGap ()
	get size of the gap between aerogel tiles

double	getAeroTubeRin ()
	get the inner size of aerogel tube

double	getAeroTubeRout ()
	get the outer size of aerogel tube

static ARICHBtestGeometryPar *	Instance ()
	Static method to get a reference to the ARICHBtestGeometryPar instance.

	ARICHBtestGeometryPar ()
	Constructor.

	ARICHBtestGeometryPar (const ARICHBtestGeometryPar &arichPar)
	copy constructor

virtual	~ARICHBtestGeometryPar ()
	Destructor.

void	setMwpc (ARICHTracking *m_mwpc)
	Set the pointer of the tracking MWPCs.

ARICHTracking *	getMwpc ()
	Get the pointer of the tracking MWPCs.

std::pair< double, double >	GetHapdChannelPosition (int)
	Get the position of the HAPD channel.

std::pair< int, int >	GetHapdElectronicMap (int)
	Get the mapping of the electronic channel to the HAPD module nr and the channel number.

int	AddHapdChannelPositionPair (double, double)
	Set the position of the HAPD channel.

int	AddHapdElectronicMapPair (int, int)
	Set the mapping of the electronic channel to the HAPD module nr and the channel number.

void	setTrackingShift (const ROOT::Math::XYZVector &)
	Set the tracking shift.

ROOT::Math::XYZVector	getTrackingShift ()
	Get the tracking shift.

void	setRotationCenter (const ROOT::Math::XYZVector &)
	Set the rotation center of the Aerogel RICH frame.

ROOT::Math::XYZVector	getRotationCenter ()
	Get the rotation center of the Aerogel RICH frame.

void	setFrameRotation (double)
	Set the rotation angle of the Aerogel RICH frame.

ROOT::Math::Rotation3D	getFrameRotation ()
	Get the rotation matrix of the Aerogel RICH frame.

void	setAverageAgel (bool)
	Set the flag for the reconstruction by using the average aerogel refractive index.

bool	getAverageAgel ()
	Get the flag for the reconstruction by using the average aerogel refractive index.

void	clear (void)
	Clears.

void	setOffset (const ROOT::Math::XYZVector &)
	Set of the setup global offset.

ROOT::Math::XYZVector	getOffset ()
	Get the setup global offset.

void	Print (void) const
	Print some debug information.

static ARICHGeometryPar *	Instance ()
	Static method to get a reference to the ARICHGeometryPar instance.

	ARICHGeometryPar ()
	Constructor.

virtual	~ARICHGeometryPar ()
	Destructor.

void	Initialize (const GearDir &content, const GearDir &mirrorinfo)
	calculates detector parameters needed for geometry build and reconstruction - the mirrorinfo can be in separate directory

void	Initialize (const GearDir &content)
	calculates detector parameters needed for geometry build and reconstruction.

void	clear (void)
	Clears.

void	read (const GearDir &content)
	gets geometry parameters from gearbox.

void	frontEndMapping (const GearDir &content)
	reads front-end board to merger to cooper mapping from an xml file

int	getMergerFromCooper (int cooperID, int finesse)
	returns merger board ID from cooperID and finesse

int	getBoardFromMerger (int mergerID, int slot)
	returns front-end board ID from merger ID and slot

int	getNBoardsOnMerger (int mergerID)
	returns number of front-end boards connected to the merger

void	readModuleInfo (const GearDir &content)
	read parameters of each module from gearbox.

double	QE (double e) const
	get photocathode quantum efficiency at energy e.

void	Print (void) const
	Print some debug information.

double	getChannelQE (int moduleID, int channelID)
	get channel quantum efficiency

int	getChannelID (ROOT::Math::XYVector hit)
	get ID number of channel containing point "hit" (hit is in detector module local coordinates)

void	modulesPosition (const GearDir &content)
	calculates the positions of HAPD modules, with the parameters from xml.

void	modulesPositionSimple (const GearDir &content)
	gets modules positions directly from xml file (for simple "beamtest" geometry).

int	getCopyNo (const ROOT::Math::XYZVector &hit)
	get the copy number of HAPD module containing point "hit"

ROOT::Math::XYZVector	getOrigin (int copyno)
	get the position of copyno-th HAPD module origin

G4ThreeVector	getOriginG4 (int copyNo)
	get the position of copyNo-th HAPD module origin (returns G4ThreeVector)

double	getModAngle (int copyno)
	get the angle of copyno-th HAPD rotation

void	chipLocPosition ()
	calculates the centers of chips in detector module local coordinates

int	getChipID (ROOT::Math::XYVector locpos)
	get ID number of chip containing point "locpos"

ROOT::Math::XYZVector	getChannelCenterGlob (int modID, int chanID)
	get center of chanID channel of modID detector module (in global coordinates)

ROOT::Math::XYVector	getChannelCenterLoc (int chID)
	get center position of chID channel (in detector module local coordinates)

void	padPositions ()
	calculates the centers of channels in local (detector module) and global coordinates

void	mirrorPositions (const GearDir &content)
	calculates parameters of all mirror planes (normal vector and point on plane)

void	mirrorPositionSimple (const GearDir &content)
	Gets mirrors positions directly from xml file (in case of simple "beamtest" geometry).

ROOT::Math::XYZVector	getMirrorNormal (int mirID)
	get normal vector of mirID-th mirror plate

ROOT::Math::XYZVector	getMirrorPoint (int mirID)
	get one point lying on mirID-th mirror plate

void	setAeroTransLength (int ilayer, double trlen)
	set transmission length of "ilayer" aerogel layer

void	setAeroRefIndex (int ilayer, double n)
	set refractive index of "ilayer" aerogel layer

void	setAerogelThickness (int ilayer, double thick)
	set thickness of "ilayer" aerogel layer

void	setAerogelZPosition (int ilayer, double zPos)
	set z position of "ilayer" aerogel layer

void	setWindowRefIndex (double refInd)
	set detector module window refractive index

void	initDetectorMask (int nmodules)
	initialize detector mask

void	setActive (int module, int channel, bool val)
	set the channel on/off

bool	isActive (int module, int channel)
	check the activity of the channel

void	readMirrorAlignment (const GearDir &content)
	Reads mirror plates alignment parameters.

int	getAerogelTileID (ROOT::Math::XYVector locpos)
	returns ID number of aerogel tile containing locpos (x-y) point.

Variables
static ARICHBtestGeometryPar *	p_B4ARICHBtestGeometryParDB = 0
	Pointer that saves the instance of this class.

static ARICHGeometryPar *	p_B4ARICHGeometryParDB = 0
	Pointer that saves the instance of this class.

Detailed Description

Function Documentation

◆ AddHapdChannelPositionPair()

int AddHapdChannelPositionPair	(	double	x,
		double	y
	)

Set the position of the HAPD channel.

Definition at line 93 of file ARICHBtestGeometryPar.cc.

  {
    m_hapdmap.push_back(std::make_pair(x, y));
    return m_hapdmap.size();
  }

◆ AddHapdElectronicMapPair()

int AddHapdElectronicMapPair	(	int	x,
		int	y
	)

Set the mapping of the electronic channel to the HAPD module nr and the channel number.

Definition at line 99 of file ARICHBtestGeometryPar.cc.

  {
    m_hapdeid.push_back(std::make_pair(x, y));
    return m_hapdeid.size();
  }

◆ ARICHBtestGeometryPar() [1/2]

ARICHBtestGeometryPar ( )

Constructor.

Definition at line 37 of file ARICHBtestGeometryPar.cc.

                                              :
    m_mwpc(new ARICHTracking),
    m_trackingShift(ROOT::Math::XYZVector()),
    m_rotationCenter(ROOT::Math::XYZVector()),
    m_frameRotation(ROOT::Math::Rotation3D()),
    m_averageAgel(0),
    m_nPads(0),
    m_init(0),
    m_simple(0)
  {
    clear();
  }

◆ ARICHBtestGeometryPar() [2/2]

ARICHBtestGeometryPar ( const ARICHBtestGeometryPar & arichPar )

copy constructor

Definition at line 51 of file ARICHBtestGeometryPar.cc.

                                                                                    :
    m_trackingShift(arichPar.m_trackingShift),
    m_rotationCenter(arichPar.m_rotationCenter),
    m_frameRotation(arichPar.m_frameRotation),
    m_hapdmap(arichPar.m_hapdmap),
    m_hapdeid(arichPar.m_hapdeid)
  {
    m_averageAgel = arichPar.m_averageAgel;
    m_nPads = arichPar.m_nPads;
    m_init = arichPar.m_init;
    m_simple = arichPar.m_simple;
    m_mwpc = new ARICHTracking();
    *m_mwpc = *arichPar.m_mwpc;
    p_B4ARICHBtestGeometryParDB = this;
  }

◆ ARICHGeometryPar()

ARICHGeometryPar ( )

Constructor.

Definition at line 41 of file ARICHGeometryPar.cc.

  {
    clear();
  }

◆ chipLocPosition()

void chipLocPosition ( )

private

calculates the centers of chips in detector module local coordinates

Definition at line 410 of file ARICHGeometryPar.cc.

  {
    double xycenter =  m_padSize * m_nPadX / 4. + m_chipGap / 2.;
    m_chipLocPos.push_back(ROOT::Math::XYVector(xycenter - m_padSize * m_nPadX / 4., xycenter - m_padSize * m_nPadX / 4.));
    m_chipLocPos.push_back(ROOT::Math::XYVector(xycenter - m_padSize * m_nPadX / 4., -xycenter - m_padSize * m_nPadX / 4.));
    m_chipLocPos.push_back(ROOT::Math::XYVector(-xycenter - m_padSize * m_nPadX / 4., xycenter - m_padSize * m_nPadX / 4.));
    m_chipLocPos.push_back(ROOT::Math::XYVector(-xycenter - m_padSize * m_nPadX / 4., -xycenter - m_padSize * m_nPadX / 4.));
  }

◆ clear() [1/2]

void clear ( void )

Clears.

Definition at line 156 of file ARICHBtestGeometryPar.cc.

  {
    m_init = false;
    m_simple = false;
    m_trackingShift = ROOT::Math::XYZVector();
    m_frameRotation   = ROOT::Math::Rotation3D();
    m_rotationCenter  = ROOT::Math::XYZVector();
  }

◆ clear() [2/2]

void clear ( void )

Clears.

Definition at line 76 of file ARICHGeometryPar.cc.

  {
    m_chipGap = 0.0;
    m_detZpos = 0.0;
    m_modXSize = 0.0;
    m_modZSize = 0.0;
    m_winThick = 0.0;
    m_winRefInd = 0.0;
    m_mirrorOuterRad = 0.0;
    m_mirrorThickness = 0.0;
    m_mirrorStartAng = 0.0;
    m_nMirrors = 0;
    m_mirrorZPos = 0.0;
    m_nPadX = 0;
    m_padSize = 0.0;
    m_detInnerRadius = 0.0;
    m_detOuterRadius = 0.0;
    m_nrow = 0;
    m_nRad = 0;
    m_init = false;
    m_simple = false;
    m_ncol.clear(); m_fDFi.clear(); m_fDR.clear(); m_fR.clear();
    m_fFi.clear(); m_fFiMod.clear(); m_chipLocPos.clear(); m_padWorldPositions.clear(); m_mirrornorm.clear(); m_mirrorpoint.clear();
    m_ChannelQE.clear();
    for (int i = 0; i < MAX_N_ALAYERS; i++) {
      m_aeroTrLength[i] = 0; m_aeroRefIndex[i] = 0;
      m_aeroZPosition[i] = 0; m_aeroThickness[i] = 0;
    }
    m_nPads = 0;
    m_ColEffi = 0.;
    m_LambdaFirst = 0.;
    m_LambdaStep = 0.;
    m_NpointsQE = 0;
    m_qeScale = 0.;
    m_windowAbsorbtion = 0.;
    m_chipNegativeCrosstalk = 0;
    for (int i = 0; i < MAXPTS_QE; i++) {m_QE[i] = 0.;}
 
    m_tileNr = 0;
    m_tileGap = 0.;
    m_aeroRin = 0.;
    m_aeroRout = 0.;
    for (int i = 0; i < 5; i++) m_tileNphi[i] = 0;
 
  }

◆ doBeamBackgroundStudy()

int doBeamBackgroundStudy ( ) const

inline

returns 1 if beam background study (to add additional sensitive modules, detect neutrons, ...), 0 else

Definition at line 289 of file ARICHGeometryConfig.h.

  {
    return m_bbstudy;
  }

◆ frontEndMapping()

void frontEndMapping ( const GearDir & content )

private

reads front-end board to merger to cooper mapping from an xml file

Definition at line 184 of file ARICHGeometryPar.cc.

  {
 
    GearDir mapping(content, "FrontEndMapping");
 
    for (const GearDir& merger : mapping.getNodes("Merger")) {
      unsigned mergerID = (unsigned) merger.getInt("@id");
 
      auto testMer = m_mergerIDs.find(mergerID);
      if (testMer != m_mergerIDs.end()) {
        B2ERROR(mapping.getPath() << "/MergerID " << mergerID <<
                " ***input already used");
      }
 
      m_mergerIDs.insert(mergerID);
 
      std::vector<unsigned> boardIDs;
      for (const GearDir& board : merger.getNodes("FEboards/FEboard")) {
        unsigned boardID = board.getInt();
        auto testBor = m_boardIDs.find(boardID);
        if (testBor != m_boardIDs.end()) {
          B2ERROR(mapping.getPath() << "/FEboardID " << boardID <<
                  " ***input already used");
        }
        boardIDs.push_back(boardID);
        m_boardIDs.insert(boardID);
      }
      m_merger2feb.insert(std::pair<int, std::vector<unsigned>>(mergerID, boardIDs));
      unsigned copperID = (unsigned) merger.getInt("COPPERid");
      string finesseSlot = merger.getString("FinesseSlot");
      int finesse = 0;
      if (finesseSlot == "A") {finesse = 0;}
      else if (finesseSlot == "B") {finesse = 1;}
      else if (finesseSlot == "C") {finesse = 2;}
      else if (finesseSlot == "D") {finesse = 3;}
      else {
        B2ERROR(merger.getPath() << "/FinesseSlot " << finesseSlot <<
                " ***invalid slot (valid are A, B, C, D)");
        continue;
      }
 
      m_copperIDs.insert(copperID);
 
      std::pair<unsigned, int> copfin(copperID, finesse);
      m_copper2merger.insert(std::pair<std::pair<unsigned, int>, unsigned>(copfin, mergerID));
 
    }
 
  }

◆ getAerogelRefIndex()

double getAerogelRefIndex ( int ilayer )

inline

get refractive index of "ilayer" aerogel layer

Definition at line 313 of file ARICHGeometryPar.h.

  {
    return m_aeroRefIndex[layer];
  }

◆ getAerogelThickness()

double getAerogelThickness ( int ilayer )

inline

get thickness of "ilayer" aerogel layer

Definition at line 318 of file ARICHGeometryPar.h.

  {
    return m_aeroThickness[layer];
  }

◆ getAerogelTileID()

int getAerogelTileID ( ROOT::Math::XYVector locpos )

returns ID number of aerogel tile containing locpos (x-y) point.

Definition at line 585 of file ARICHGeometryPar.cc.

  {
 
    double size = (m_aeroRout - m_aeroRin) / double(m_tileNr);
    double r = locpos.R();
    double phi = ROOT::Math::VectorUtil::Phi_0_2pi(locpos.Phi());
    int nr = int((r - m_aeroRin) / size);
 
    if (r < m_aeroRin + nr * size + m_tileGap / 2. || r >  m_aeroRin + (nr + 1)*size - m_tileGap / 2.) return 0;
 
    double dphi = 2.*M_PI / double(m_tileNphi[nr]);
    double gapPhi = m_tileGap / (m_aeroRin + (nr + 1) * size - m_tileGap / 2.);
 
    int nphi = int(phi / dphi);
    if (phi < nphi * dphi + gapPhi / 2. || phi > (nphi + 1)*dphi - gapPhi / 2.) return 0;
 
    int tileID = nphi;
    for (int i = 0; i < nr; i++) tileID += m_tileNphi[i];
 
    return tileID + 1;
 
  }

◆ getAerogelTransmissionLength()

double getAerogelTransmissionLength ( int ilayer )

inline

get transmission length of "ilayer" aerogel layer

Definition at line 308 of file ARICHGeometryPar.h.

  {
    return m_aeroTrLength[layer];
  }

◆ getAerogelZPosition()

double getAerogelZPosition ( int ilayer )

inline

get z position of "ilayer" aerogel layer

Definition at line 323 of file ARICHGeometryPar.h.

  {
    return m_aeroZPosition[layer];
  }

◆ getAeroTubeRin()

double getAeroTubeRin ( )

inline

get the inner size of aerogel tube

Definition at line 404 of file ARICHGeometryPar.h.

  {
    return m_aeroRin;
  }

◆ getAeroTubeRout()

double getAeroTubeRout ( )

inline

get the outer size of aerogel tube

Definition at line 409 of file ARICHGeometryPar.h.

  {
    return m_aeroRout;
  }

◆ getAverageAgel()

bool getAverageAgel ( )

Get the flag for the reconstruction by using the average aerogel refractive index.

Definition at line 149 of file ARICHBtestGeometryPar.cc.

  {
    return m_averageAgel;
  }

◆ getBoardFromMerger()

int getBoardFromMerger	(	int	mergerID,
		int	slot
	)

returns front-end board ID from merger ID and slot

Definition at line 246 of file ARICHGeometryPar.cc.

  {
    auto boards = m_merger2feb.find(mergerID);
    if (boards == m_merger2feb.end()) {
      B2ERROR("getBoardFromMerger: " << " merger " << mergerID << " is not mapped");
      return -1;
    }
    if ((boards->second).size() <= unsigned(slot)) {
      B2ERROR("getBoardFromMerger: " << " merger " << mergerID << " slot " << slot << " is not assigned to FE board.");
      return -1;
    }
    return (boards->second).at(slot);
  }

◆ getChannelCenterGlob()

ROOT::Math::XYZVector getChannelCenterGlob	(	int	modID,
		int	chanID
	)

get center of chanID channel of modID detector module (in global coordinates)

Definition at line 431 of file ARICHGeometryPar.cc.

  {
    int id = (modID - 1) * m_nPads + chanID;
    return ROOT::Math::XYZVector(m_padWorldPositions.at(id).X(), m_padWorldPositions.at(id).Y(), m_detZpos + m_winThick);
  }

◆ getChannelCenterLoc()

ROOT::Math::XYVector getChannelCenterLoc ( int chID )

get center position of chID channel (in detector module local coordinates)

Definition at line 437 of file ARICHGeometryPar.cc.

  {
    return m_padLocPositions[chID];
  }

◆ getChannelID()

int getChannelID ( ROOT::Math::XYVector hit )

get ID number of channel containing point "hit" (hit is in detector module local coordinates)

Definition at line 313 of file ARICHGeometryPar.cc.

  {
    int ChipID = getChipID(position);
    int Npad = int(m_nPadX / 2);
    ROOT::Math::XYVector chipPos = getChipLocPos(ChipID);
    ROOT::Math::XYVector locloc = position - chipPos;
    int ix = int(locloc.X() / m_padSize);
    int iy = int(locloc.Y() / m_padSize);
    if (locloc.X() < 0 || locloc.Y() < 0) return -1;
    if (ix > Npad - 1 || iy > Npad - 1) return -1;
    int chID = ChipID * Npad * Npad + iy + ix * Npad;
    return chID;
  }

◆ getChannelQE()

double getChannelQE	(	int	moduleID,
		int	channelID
	)

get channel quantum efficiency

Definition at line 307 of file ARICHGeometryPar.cc.

  {
    int id = (moduleID - 1) * m_nPads + channelID;
    return m_ChannelQE.at(id) / 100.0;
  }

◆ getChipID()

int getChipID ( ROOT::Math::XYVector locpos )

get ID number of chip containing point "locpos"

Definition at line 420 of file ARICHGeometryPar.cc.

  {
    if (locpos.X() > 0) {
      if (locpos.Y() > 0) return 0;
      return 1;
    }
    if (locpos.Y() > 0) return 2;
    return 3;
  }

◆ getChipLocPos()

ROOT::Math::XYVector getChipLocPos ( int chipID )

inline

get center position of chipID-th chip of detector module (in detector module local coordinates)

Definition at line 292 of file ARICHGeometryPar.h.

  {
    return m_chipLocPos.at(chipID);
  }

◆ getChipNegativeCrosstalk()

double getChipNegativeCrosstalk ( )

inline

get factor to suppress p.e. detection efficiency due to negative polarity crosstalk among chip channels

Definition at line 383 of file ARICHGeometryPar.h.

  {
    return m_chipNegativeCrosstalk;
  }

◆ getColEffi()

double getColEffi ( ) const

inline

get HAPD collection efficiency.

Definition at line 277 of file ARICHGeometryPar.h.

  {
    return m_ColEffi;
  }

◆ getCopyNo()

int getCopyNo ( const ROOT::Math::XYZVector & hit )

get the copy number of HAPD module containing point "hit"

Definition at line 375 of file ARICHGeometryPar.cc.

  {
    double x = hit.X();
    double y = hit.Y();
    double r = sqrt(x * x + y * y);
    double fi = atan2(y, x);
    if (fi < 0) fi += 2 * M_PI;
    int ntot = 0;
    for (int i = 0; i < m_nrow; i++) {
      int nfi = int(fi / m_fDFi[i]);
      int copyno = ntot + nfi;
      if (fabs(r - m_fR[copyno]) <  m_modXSize / 2.) return copyno + 1;
      ntot += m_ncol[i];
    }
    return -1;
  }

◆ getDetectorPadSize()

double getDetectorPadSize ( )

inline

get detectors module pad size

Definition at line 348 of file ARICHGeometryPar.h.

  {
    return m_padSize;
  }

◆ getDetectorWindowRefIndex()

double getDetectorWindowRefIndex ( )

inline

get refractive index of detector module window

Definition at line 338 of file ARICHGeometryPar.h.

  {
    return m_winRefInd;
  }

◆ getDetectorWindowThickness()

double getDetectorWindowThickness ( )

inline

get thickness of detector module window

Definition at line 333 of file ARICHGeometryPar.h.

  {
    return m_winThick;
  }

◆ getDetectorXPadNumber()

int getDetectorXPadNumber ( )

inline

get number of pads of detector module (in one direction)

Definition at line 353 of file ARICHGeometryPar.h.

  {
    return m_nPadX;
  }

◆ getDetectorZPosition() [1/2]

double getDetectorZPosition ( )

inline

get z position of detector plane (starting z)

Definition at line 368 of file ARICHGeometryPar.h.

  {
    return m_detZpos;
  }

◆ getDetectorZPosition() [2/2]

double getDetectorZPosition ( ) const

inline

get z position of detector plane (starting z of HAPDs)

Definition at line 284 of file ARICHGeometryConfig.h.

  {
    return m_detectorPlane.getPosition().Z() - m_hapd.getModuleSizeZ() / 2.;
  }

◆ getFrameRotation()

ROOT::Math::Rotation3D getFrameRotation ( )

Get the rotation matrix of the Aerogel RICH frame.

Definition at line 139 of file ARICHBtestGeometryPar.cc.

  {
    return m_frameRotation;
  }

◆ GetHapdChannelPosition()

std::pair< double, double > GetHapdChannelPosition ( int id )

Get the position of the HAPD channel.

Definition at line 83 of file ARICHBtestGeometryPar.cc.

  {
    return m_hapdmap[id];
  }

◆ GetHapdElectronicMap()

std::pair< int, int > GetHapdElectronicMap ( int id )

Get the mapping of the electronic channel to the HAPD module nr and the channel number.

Definition at line 88 of file ARICHBtestGeometryPar.cc.

  {
    return m_hapdeid[id];
  }

◆ getMergerFromCooper()

int getMergerFromCooper	(	int	cooperID,
		int	finesse
	)

returns merger board ID from cooperID and finesse

Definition at line 234 of file ARICHGeometryPar.cc.

  {
    auto merger = m_copper2merger.find(std::pair<unsigned, int>(copperID, finesse));
    if (merger == m_copper2merger.end()) {
      // B2INFO("getMergerFromCooper: " << " copper " << copperID << ", finesse "
      //       << finesse << " is not assigned to any merger board");
      return 0;
    }
 
    return merger->second;
  }

◆ getMirrorNormal()

ROOT::Math::XYZVector getMirrorNormal ( int mirID )

get normal vector of mirID-th mirror plate

Definition at line 503 of file ARICHGeometryPar.cc.

  {
    return m_mirrornorm[mirID];
  }

◆ getMirrorPoint()

ROOT::Math::XYZVector getMirrorPoint ( int mirID )

get one point lying on mirID-th mirror plate

Definition at line 508 of file ARICHGeometryPar.cc.

  {
    return m_mirrorpoint[mirID];
  }

◆ getMirrorsStartAngle()

double getMirrorsStartAngle ( )

inline

get the phi angle of first corner of mirror polygon

Definition at line 358 of file ARICHGeometryPar.h.

  {
    return m_mirrorStartAng;
  }

◆ getMirrorsZPosition()

double getMirrorsZPosition ( )

inline

get the z position of mirror plates (starting z)

Definition at line 363 of file ARICHGeometryPar.h.

  {
    return m_mirrorZPos;
  }

◆ getModAngle()

double getModAngle ( int copyno )

get the angle of copyno-th HAPD rotation

Definition at line 405 of file ARICHGeometryPar.cc.

  {
    return m_fFiMod[copyno - 1];
  }

◆ getMwpc()

ARICHTracking * getMwpc ( )

Get the pointer of the tracking MWPCs.

Definition at line 78 of file ARICHBtestGeometryPar.cc.

  {
    return m_mwpc;
  }

◆ getNBoardsOnMerger()

int getNBoardsOnMerger ( int mergerID )

returns number of front-end boards connected to the merger

Definition at line 260 of file ARICHGeometryPar.cc.

  {
    auto boards = m_merger2feb.find(mergerID);
    if (boards == m_merger2feb.end()) {
      B2ERROR("getNBoardsOnMerger: " << " merger " << mergerID << " is not mapped");
    }
    return (boards->second).size();
  }

◆ getNMCopies()

int getNMCopies ( ) const

inline

get the total number of HAPD modules

Definition at line 287 of file ARICHGeometryPar.h.

  {
    return m_fR.size();
  }

◆ getNMirrors()

int getNMirrors ( )

inline

get number of mirror plates

Definition at line 328 of file ARICHGeometryPar.h.

  {
    return m_nMirrors;
  }

◆ getNphiTiles()

int getNphiTiles ( int nRing )

inline

get number of aerogel tiles in each radial ring

Definition at line 393 of file ARICHGeometryPar.h.

  {
    if (nRing > 4) return 0;
    return m_tileNphi[nRing];
  }

◆ getNrTiles()

int getNrTiles ( )

inline

get number of aerogel tiles in radial direction (number of rings)

Definition at line 388 of file ARICHGeometryPar.h.

  {
    return m_tileNr;
  }

◆ getNumberOfAerogelRadiators()

int getNumberOfAerogelRadiators ( )

inline

get number of aerogel layers

Definition at line 343 of file ARICHGeometryPar.h.

  {
    return m_nRad;
  }

◆ getOffset()

ROOT::Math::XYZVector getOffset ( )

Get the setup global offset.

Definition at line 172 of file ARICHBtestGeometryPar.cc.

  {
    return m_offset;
  }

◆ getOrigin()

ROOT::Math::XYZVector getOrigin ( int copyno )

get the position of copyno-th HAPD module origin

Definition at line 392 of file ARICHGeometryPar.cc.

  {
    const double amag = std::abs(m_fR[copyNo - 1]);
    const double phi = m_fFi[copyNo - 1];
    return ROOT::Math::XYZVector(amag * std::cos(phi), amag * std::sin(phi), m_detZpos + m_modZSize / 2.);
  }

◆ getOriginG4()

G4ThreeVector getOriginG4 ( int copyNo )

get the position of copyNo-th HAPD module origin (returns G4ThreeVector)

Definition at line 399 of file ARICHGeometryPar.cc.

  {
    ROOT::Math::XYZVector origin = getOrigin(copyNo);
    return G4ThreeVector(origin.X() / Unit::mm, origin.Y() / Unit::mm, origin.Z() / Unit::mm);
  }

◆ getQEScaling()

double getQEScaling ( )

inline

get QE scaling factor for photons internally reflected in HAPD window

Definition at line 373 of file ARICHGeometryPar.h.

  {
    return m_qeScale;
  }

◆ getRotationCenter()

ROOT::Math::XYZVector getRotationCenter ( )

Get the rotation center of the Aerogel RICH frame.

Definition at line 127 of file ARICHBtestGeometryPar.cc.

  {
    return m_rotationCenter;
  }

◆ getSensitiveSurfaceSize()

double getSensitiveSurfaceSize ( ) const

inline

get size of detector sensitive surface (size of two chips + gap between)

Definition at line 282 of file ARICHGeometryPar.h.

  {
    return m_nPadX * m_padSize + m_chipGap;
  }

◆ getTileGap()

double getTileGap ( )

inline

get size of the gap between aerogel tiles

Definition at line 399 of file ARICHGeometryPar.h.

  {
    return m_tileGap;
  }

◆ getTrackingShift()

ROOT::Math::XYZVector getTrackingShift ( )

Get the tracking shift.

Definition at line 114 of file ARICHBtestGeometryPar.cc.

  {
    return m_trackingShift;
  }

◆ getWindowAbsorbtion()

double getWindowAbsorbtion ( )

inline

get absorbtion probability for photons internally reflected in HAPD window

Definition at line 378 of file ARICHGeometryPar.h.

  {
    return m_windowAbsorbtion;
  }

◆ initDetectorMask()

void initDetectorMask ( int nmodules )

initialize detector mask

Definition at line 539 of file ARICHGeometryPar.cc.

  {
    int m_DetectorMaskSize = m_nPads * nmodules / 32 + 1;
    for (int i = 0; i < m_DetectorMaskSize; i++) {
      m_DetectorMask.push_back(0xFFFFFFFF);
    }
    B2INFO("DetectorMask initialized size=" << m_nPads << " * " << nmodules << " +1 =" << m_DetectorMaskSize);
 
  }

◆ Initialize() [1/2]

void Initialize ( const GearDir & content )

calculates detector parameters needed for geometry build and reconstruction.

Definition at line 71 of file ARICHGeometryPar.cc.

  {
    Initialize(content, content);
  }

◆ Initialize() [2/2]

void Initialize	(	const GearDir &	content,
		const GearDir &	mirrorinfo
	)

calculates detector parameters needed for geometry build and reconstruction - the mirrorinfo can be in separate directory

Definition at line 50 of file ARICHGeometryPar.cc.

  {
 
    read(content);
    string Type = content.getString("@type", "");
    if (Type == "beamtest") {
      m_simple = true;
      modulesPositionSimple(content);
      mirrorPositionSimple(mirrorcontent);
    } else {
      modulesPosition(content);
      mirrorPositions(content);
      frontEndMapping(content);
    }
    chipLocPosition();
    padPositions();
    initDetectorMask(m_fR.size());
    readModuleInfo(content);
    m_init = true;
  }

◆ Instance() [1/2]

ARICHBtestGeometryPar * Instance ( )

static

Static method to get a reference to the ARICHBtestGeometryPar instance.

Returns: A reference to an instance of this class.

Definition at line 29 of file ARICHBtestGeometryPar.cc.

  {
    if (!p_B4ARICHBtestGeometryParDB) {
      p_B4ARICHBtestGeometryParDB = new ARICHBtestGeometryPar();
    }
    return p_B4ARICHBtestGeometryParDB;
  }

◆ Instance() [2/2]

ARICHGeometryPar * Instance ( )

static

Static method to get a reference to the ARICHGeometryPar instance.

Returns: A reference to an instance of this class.

Definition at line 33 of file ARICHGeometryPar.cc.

  {
    if (!p_B4ARICHGeometryParDB) {
      p_B4ARICHGeometryParDB = new ARICHGeometryPar();
    }
    return p_B4ARICHGeometryParDB;
  }

◆ isActive()

bool isActive	(	int	module,
		int	channel
	)

check the activity of the channel

Definition at line 561 of file ARICHGeometryPar.cc.

  {
    int ch  = (module - 1) * m_nPads + channel;
    int bit = ch % 32;
    int idx = ch / 32;
    return    m_DetectorMask[idx] & (1 << bit);
  }

◆ isInit()

bool isInit ( )

inline

returns true if parametrisation is initialized, false otherwise

Definition at line 297 of file ARICHGeometryPar.h.

  {
    return m_init;
  }

◆ isSimple()

bool isSimple ( )

inline

returns true if parametrisation is initialized with simple geometry (beamtest)

Definition at line 302 of file ARICHGeometryPar.h.

  {
    return m_simple;
  }

◆ mirrorPositions()

void mirrorPositions ( const GearDir & content )

private

calculates parameters of all mirror planes (normal vector and point on plane)

Calculates mirror positions (normal vectors and on point of every mirror plate) and stores them.

Definition at line 477 of file ARICHGeometryPar.cc.

  {
    double rmir = m_mirrorOuterRad * cos(M_PI / m_nMirrors) - m_mirrorThickness;
    for (int i = 0; i < m_nMirrors; i++) {
      ROOT::Math::XYZVector norm(cos(2.*M_PI / double(m_nMirrors) * (i + 0.5) + m_mirrorStartAng),
                                 sin(2.*M_PI / double(m_nMirrors) * (i + 0.5) + m_mirrorStartAng), 0);
      m_mirrornorm.push_back(norm);
      m_mirrorpoint.push_back(rmir * norm);
    }
    readMirrorAlignment(content);
  }

◆ mirrorPositionSimple()

void mirrorPositionSimple ( const GearDir & content )

private

Gets mirrors positions directly from xml file (in case of simple "beamtest" geometry).

Definition at line 489 of file ARICHGeometryPar.cc.

  {
    double thick = content.getLength("Mirrors/thickness");
    BOOST_FOREACH(const GearDir & mirror, content.getNodes("Mirrors/Mirror")) {
      double angle = mirror.getAngle("angle");
      ROOT::Math::XYZVector point(mirror.getLength("xPos") - cos(angle)*thick / 2., mirror.getLength("yPos") - sin(angle)*thick / 2., 0);
      ROOT::Math::XYZVector norm(cos(angle), sin(angle), 0);
      m_mirrorpoint.push_back(point);
      m_mirrornorm.push_back(norm);
      m_mirrorZPos = mirror.getLength("zPos");
      m_nMirrors++;
    }
  }

◆ modulesPosition()

void modulesPosition ( const GearDir & content )

private

calculates the positions of HAPD modules, with the parameters from xml.

Definition at line 327 of file ARICHGeometryPar.cc.

  {
 
    GearDir detParams(content, "Detector/Plane/Rings");
 
    double r = m_detInnerRadius;
 
    BOOST_FOREACH(const GearDir & ring, detParams.getNodes("Ring")) {
      double dR = ring.getLength("dR");
      r += dR;
      double rcenter = r + m_modXSize / 2.;
      if (rcenter + m_modXSize * sqrt(2) / 2. > m_detOuterRadius) {
        B2WARNING(m_ncol.size() + 1  << "th ring of ARICH photon detectors will not be placed (out of detector tube).");
        break;
      }
      m_nrow++;
      int nSeg = ring.getInt("nSegments") ;
      double dFi = ring.getLength("dFi");
      m_fDR.push_back(dR);
      double f = 2.*atan2((m_modXSize + dFi) / 2., r);
      int blaa = int(2.*M_PI / f / nSeg) * nSeg;
      m_ncol.push_back(blaa);
      f = 2.*M_PI / double(blaa);
      m_fDFi.push_back(f);
      B2INFO(blaa << " modules of " << m_ncol.size() << "th ring of ARICH photon detectors will be placed at r = " << rcenter << "cm. ");
      for (int nv = 0; nv < blaa; ++nv) {
        m_fR.push_back(rcenter);
        double fi = f * (nv + 0.5);
        m_fFi.push_back(fi);
        m_fFiMod.push_back(fi);
      }
      r += (m_modXSize  + r * (1 - cos(f / 2.)));
    }
    B2INFO("Altogether " << m_fR.size() << " ARICH photon detector modules will be placed.");
  }

◆ modulesPositionSimple()

void modulesPositionSimple ( const GearDir & content )

private

gets modules positions directly from xml file (for simple "beamtest" geometry).

Definition at line 363 of file ARICHGeometryPar.cc.

  {
    BOOST_FOREACH(const GearDir & module, content.getNodes("Detector/Plane/Modules/Module")) {
      ROOT::Math::XYVector position(module.getLength("xPos"), module.getLength("yPos"));
      double angle = module.getAngle("angle") / Unit::rad;
      m_fFi.push_back(position.Phi());
      m_fR.push_back(position.R());
      m_fFiMod.push_back(angle);
    }
    B2INFO("Altogether " << m_fR.size() << " ARICH photon detector modules will be placed.");
  }

◆ padPositions()

void padPositions ( )

private

calculates the centers of channels in local (detector module) and global coordinates

Definition at line 443 of file ARICHGeometryPar.cc.

  {
    int Npad = int(m_nPadX / 2.);
    ROOT::Math::XYVector xstart(m_padSize / 2., m_padSize / 2.);
    for (int chipID = 0; chipID < 4; chipID++) {
      ROOT::Math::XYVector chipPos = getChipLocPos(chipID);
      for (int ix = 0; ix < Npad; ix++) {
        for (int iy = 0; iy < Npad; iy++) {
          int chanID = chipID * Npad * Npad + ix * Npad + iy;
          ROOT::Math::XYVector center(m_padSize / 2. + ix * m_padSize, m_padSize / 2. + iy * m_padSize);
          center = center + chipPos;
          m_padLocPositions[chanID] = center;
        }
      }
    }
    for (int iMod = 0; iMod < getNMCopies(); iMod++) {
      const double magnitude = std::abs(m_fR[iMod]);
      const double phi = m_fFi[iMod];
      const double phiMod = m_fFiMod[iMod];
      const double iModCenterX = magnitude * std::cos(phi);
      const double iModCenterY = magnitude * std::sin(phi);
 
      for (unsigned int iChan = 0; iChan < m_padLocPositions.size(); iChan++) {
        const ROOT::Math::XYVector& iChanCenter = m_padLocPositions[iChan];
        const double cosPhiMod = std::cos(phiMod);
        const double sinPhiMod = std::sin(phiMod);
        const double iChanCenterRotatedX = iChanCenter.X() * cosPhiMod - iChanCenter.Y() * sinPhiMod;
        const double iChanCenterRotatedY = iChanCenter.X() * sinPhiMod + iChanCenter.Y() * cosPhiMod;
        ROOT::Math::XYVector iWorld(iModCenterX + iChanCenterRotatedX, iModCenterY + iChanCenterRotatedY);
        m_padWorldPositions.push_back(iWorld);
      }
    }
  }

◆ Print() [1/2]

void Print ( void ) const

Print some debug information.

Definition at line 178 of file ARICHBtestGeometryPar.cc.

  {
 
  }

◆ Print() [2/2]

void Print ( void ) const

Print some debug information.

Definition at line 303 of file ARICHGeometryPar.cc.

304 {

305 }

◆ QE()

double QE ( double e ) const

get photocathode quantum efficiency at energy e.

Definition at line 292 of file ARICHGeometryPar.cc.

  {
    if (e < 0.001) return 0;
    double dlam = 1240 / e - m_LambdaFirst;
    if (dlam < 0) return 0;
    int i = int(dlam / m_LambdaStep);
    if (i > m_NpointsQE - 2) return 0;
    return m_QE[i] + (m_QE[i + 1] - m_QE[i]) / m_LambdaStep * (dlam - i * m_LambdaStep);
  }

◆ read()

void read ( const GearDir & content )

gets geometry parameters from gearbox.

Definition at line 122 of file ARICHGeometryPar.cc.

  {
 
    //------------------------------
    // Get ARICH geometry parameters
    //------------------------------
    GearDir detParams(content, "Detector");
    m_modXSize = detParams.getLength("Module/moduleXSize");
    m_modZSize = detParams.getLength("Module/moduleZSize");
    m_winThick = detParams.getLength("Module/windowThickness");
    m_nPadX = detParams.getInt("Module/padXNum");
    m_nPads = m_nPadX * m_nPadX;
    m_padSize = detParams.getLength("Module/padSize");
    m_chipGap = detParams.getLength("Module/chipGap");
    m_detZpos = detParams.getLength("Plane/zPosition");
    m_qeScale = detParams.getDouble("Module/qeScale");
    m_windowAbsorbtion = detParams.getDouble("Module/windowAbsorbtion");
    m_chipNegativeCrosstalk = detParams.getDouble("Module/chipNegativeCrosstalk");
    string Type = content.getString("@type", "");
    if (Type == "beamtest") return;
    m_detInnerRadius = detParams.getLength("Plane/tubeInnerRadius");
    m_detOuterRadius = detParams.getLength("Plane/tubeOuterRadius");
 
    GearDir mirrParams(content, "Mirrors");
    if (mirrParams) {
      m_nMirrors = mirrParams.getInt("nMirrors");
      m_mirrorThickness =  mirrParams.getLength("mirrorThickness");
      m_mirrorOuterRad = mirrParams.getLength("outerRadius");
      m_mirrorStartAng = mirrParams.getAngle("startAngle");
      m_mirrorZPos = mirrParams.getLength("Zposition");
    }
 
    GearDir qeParams(content, "QE");
    m_ColEffi = qeParams.getDouble("ColEffi");
    m_LambdaFirst = qeParams.getLength("LambdaFirst") / Unit::nm;
    m_LambdaStep = qeParams.getLength("LambdaStep") / Unit::nm;
    m_NpointsQE = 0;
    for (int i = 0; i < MAXPTS_QE; i++) {
      int ii = i + 1;
      stringstream ss; string cc;
      ss << ii; ss >> cc;
      string path = "Qeffi[@component='point-" + cc + "']/";
      GearDir qe(qeParams, path);
      if (!qe) break;
      m_NpointsQE++;
      m_QE[i] = qe.getDouble("qe");
    }
 
    GearDir aerogel(content, "Aerogel");
    m_tileNr = aerogel.getInt("tileNr");
    m_tileGap = aerogel.getLength("tileGap") / Unit::cm;
    m_aeroRin = aerogel.getLength("tubeInnerRadius") / Unit::cm;
    m_aeroRout = aerogel.getLength("tubeOuterRadius") / Unit::cm;
 
    int i = 0;
    for (const GearDir& ring : aerogel.getNodes("tileNphi/Ring")) {
      m_tileNphi[i] = ring.getInt();
      i++;
    }
 
  }

◆ readMirrorAlignment()

void readMirrorAlignment ( const GearDir & content )

private

Reads mirror plates alignment parameters.

Definition at line 569 of file ARICHGeometryPar.cc.

  {
    GearDir modParams(content, "Mirrors/Alignment");
 
    BOOST_FOREACH(const GearDir & plate, modParams.getNodes("Plate")) {
      int id = atoi(plate.getString("@id").c_str());
      double dr = plate.getLength("dr");
      double dphi = plate.getAngle("dphi");
      double dtheta = plate.getAngle("dtheta");
      VectorUtil::setMag(m_mirrorpoint[id - 1], m_mirrorpoint[id - 1].R() + dr);
      VectorUtil::setTheta(m_mirrornorm[id - 1], m_mirrornorm[id - 1].Theta() + dtheta);
      VectorUtil::setPhi(m_mirrornorm[id - 1], m_mirrornorm[id - 1].Phi() + dphi);
    }
  }

◆ readModuleInfo()

void readModuleInfo ( const GearDir & content )

read parameters of each module from gearbox.

Definition at line 269 of file ARICHGeometryPar.cc.

  {
    istringstream chstream;
    int ch; double qq;
    GearDir modParams(content, "ModuleInfo");
    uint8_t defqe = uint8_t(modParams.getDouble("DefaultQE") * 100);
    m_ChannelQE.assign(m_nPads * m_fR.size(), defqe);
    BOOST_FOREACH(const GearDir & module, modParams.getNodes("Module")) {
      int modid = atoi(module.getString("@id", "").c_str());
      chstream.str(module.getString("ChannelsQE"));
      while (chstream >> ch >> qq) {
        int chid = (modid - 1) * m_nPads + ch;
        m_ChannelQE[chid] = uint8_t(qq * 100);
      }
      chstream.clear();
      chstream.str(module.getString("DeadChannels"));
      while (chstream >> ch) {
        setActive(modid, ch, false);
      }
      chstream.clear();
    }
  }

◆ setActive()

void setActive	(	int	module,
		int	channel,
		bool	val
	)

set the channel on/off

Definition at line 549 of file ARICHGeometryPar.cc.

  {
    int ch  = (module - 1) * m_nPads + channel;
    int bit = ch % 32;
    unsigned int idx = ch / 32;
    if (idx >= m_DetectorMask.size()) {
      B2WARNING(idx  << " Wrong detector mask size >= " << m_DetectorMask.size());
    }
    if (active) m_DetectorMask[idx] |= (1 << bit);
    else        m_DetectorMask[idx] &= ~(1 << bit);
  }

◆ setAerogelThickness()

void setAerogelThickness	(	int	ilayer,
		double	thick
	)

set thickness of "ilayer" aerogel layer

Definition at line 523 of file ARICHGeometryPar.cc.

  {
    m_nRad++;
    m_aeroThickness[layer] = thick;
  }

◆ setAerogelZPosition()

void setAerogelZPosition	(	int	ilayer,
		double	zPos
	)

set z position of "ilayer" aerogel layer

Definition at line 529 of file ARICHGeometryPar.cc.

  {
    m_aeroZPosition[layer] = zPos;
  }

◆ setAeroRefIndex()

void setAeroRefIndex	(	int	ilayer,
		double	n
	)

set refractive index of "ilayer" aerogel layer

Definition at line 518 of file ARICHGeometryPar.cc.

  {
    if (n) m_aeroRefIndex[layer] = n;
  }

◆ setAeroTransLength()

void setAeroTransLength	(	int	ilayer,
		double	trlen
	)

set transmission length of "ilayer" aerogel layer

Definition at line 513 of file ARICHGeometryPar.cc.

  {
    m_aeroTrLength[layer] = trlength;
  }

◆ setAverageAgel()

void setAverageAgel ( bool v )

Set the flag for the reconstruction by using the average aerogel refractive index.

Definition at line 144 of file ARICHBtestGeometryPar.cc.

  {
    m_averageAgel = v;
  }

◆ setFrameRotation()

void setFrameRotation ( double v )

Set the rotation angle of the Aerogel RICH frame.

Definition at line 132 of file ARICHBtestGeometryPar.cc.

  {
    m_frameRotation = ROOT::Math::Rotation3D() * ROOT::Math::RotationY(v);
    B2INFO("ARICHBtestGeometryPar::setFrameRotation=" << v);
 
  }

◆ setMwpc()

void setMwpc ( ARICHTracking * m_mwpc )

Set the pointer of the tracking MWPCs.

Definition at line 72 of file ARICHBtestGeometryPar.cc.

  {
    m_mwpc = m;
  }

◆ setOffset()

void setOffset ( const ROOT::Math::XYZVector & v )

Set of the setup global offset.

Definition at line 165 of file ARICHBtestGeometryPar.cc.

  {
    m_offset = v;
    B2INFO("ARICHBtestGeometryPar::setOffset=" << m_offset.X() << ","  << m_offset.Y() << ","  <<  m_offset.Z());
 
  }

◆ setRotationCenter()

void setRotationCenter ( const ROOT::Math::XYZVector & v )

Set the rotation center of the Aerogel RICH frame.

Definition at line 119 of file ARICHBtestGeometryPar.cc.

  {
    m_rotationCenter = v;
    B2INFO("ARICHBtestGeometryPar::setRotationCenter=" << m_rotationCenter.X() << ","  << m_rotationCenter.Y() << ","  <<
           m_rotationCenter.Z());
 
  }

◆ setTrackingShift()

void setTrackingShift ( const ROOT::Math::XYZVector & v )

Set the tracking shift.

Definition at line 105 of file ARICHBtestGeometryPar.cc.

  {
    m_trackingShift = v;
    B2INFO("ARICHBtestGeometryPar::setTrackingShift=" << m_trackingShift.X() << ","  << m_trackingShift.Y() << ","  <<
           m_trackingShift.Z());
 
  }

◆ setWindowRefIndex()

void setWindowRefIndex ( double refInd )

set detector module window refractive index

Definition at line 534 of file ARICHGeometryPar.cc.

  {
    m_winRefInd = refInd;
  }

◆ ~ARICHBtestGeometryPar()

~ARICHBtestGeometryPar ( )

virtual

Destructor.

Definition at line 67 of file ARICHBtestGeometryPar.cc.

68 {

69 }

◆ ~ARICHGeometryPar()

~ARICHGeometryPar ( )

virtual

Destructor.

Definition at line 46 of file ARICHGeometryPar.cc.

47 {

48 }

Variable Documentation

◆ p_B4ARICHBtestGeometryParDB

ARICHBtestGeometryPar * p_B4ARICHBtestGeometryParDB = 0

staticprivate

Pointer that saves the instance of this class.

Definition at line 196 of file ARICHBtestGeometryPar.h.

◆ p_B4ARICHGeometryParDB

ARICHGeometryPar * p_B4ARICHGeometryParDB = 0

staticprivate

Pointer that saves the instance of this class.

Definition at line 248 of file ARICHGeometryPar.h.

Modules

Classes

Functions

Variables

Detailed Description

Function Documentation

◆ AddHapdChannelPositionPair()

◆ AddHapdElectronicMapPair()

◆ ARICHBtestGeometryPar() [1/2]

◆ ARICHBtestGeometryPar() [2/2]

◆ ARICHGeometryPar()

◆ chipLocPosition()

◆ clear() [1/2]

◆ clear() [2/2]

◆ doBeamBackgroundStudy()

◆ frontEndMapping()

◆ getAerogelRefIndex()

◆ getAerogelThickness()

◆ getAerogelTileID()

◆ getAerogelTransmissionLength()

◆ getAerogelZPosition()

◆ getAeroTubeRin()

◆ getAeroTubeRout()

◆ getAverageAgel()

◆ getBoardFromMerger()

◆ getChannelCenterGlob()

◆ getChannelCenterLoc()

◆ getChannelID()

◆ getChannelQE()

◆ getChipID()

◆ getChipLocPos()

◆ getChipNegativeCrosstalk()

◆ getColEffi()

◆ getCopyNo()

◆ getDetectorPadSize()

◆ getDetectorWindowRefIndex()

◆ getDetectorWindowThickness()

◆ getDetectorXPadNumber()

◆ getDetectorZPosition() [1/2]

◆ getDetectorZPosition() [2/2]

◆ getFrameRotation()

◆ GetHapdChannelPosition()

◆ GetHapdElectronicMap()

◆ getMergerFromCooper()

◆ getMirrorNormal()

◆ getMirrorPoint()

◆ getMirrorsStartAngle()

◆ getMirrorsZPosition()

◆ getModAngle()

◆ getMwpc()

◆ getNBoardsOnMerger()

◆ getNMCopies()

◆ getNMirrors()

◆ getNphiTiles()

◆ getNrTiles()

◆ getNumberOfAerogelRadiators()

◆ getOffset()

◆ getOrigin()

◆ getOriginG4()

◆ getQEScaling()

◆ getRotationCenter()

◆ getSensitiveSurfaceSize()

◆ getTileGap()

◆ getTrackingShift()

◆ getWindowAbsorbtion()

◆ initDetectorMask()

◆ Initialize() [1/2]

◆ Initialize() [2/2]

◆ Instance() [1/2]

◆ Instance() [2/2]

◆ isActive()

◆ isInit()

◆ isSimple()

◆ mirrorPositions()

◆ mirrorPositionSimple()

◆ modulesPosition()

◆ modulesPositionSimple()

◆ padPositions()

◆ Print() [1/2]

◆ Print() [2/2]