TgcDigitMaker Class Reference

#include <TgcDigitMaker.h>

Inheritance diagram for TgcDigitMaker:

Collaboration diagram for TgcDigitMaker:

Public Member Functions
	TgcDigitMaker (const TgcHitIdHelper *hitIdHelper, const MuonGM::MuonDetectorManager *mdManager, unsigned int runperiod, const bool doFourBunch)
virtual	~TgcDigitMaker ()
StatusCode	initialize ()
	Initializes TgcHitIdHelper, TgcIdHelper and random number of a stream for the digitization. More...
TgcDigitCollection *	executeDigi (const TGCSimHit *hit, const double globalHitTime, const TgcDigitASDposData *ASDpos, const TgcDigitTimeOffsetData *TOffset, const TgcDigitCrosstalkData *Crosstalk, CLHEP::HepRandomEngine *rndmEngine)
	A single hit can be digitized in the two directions independently: R and phi directions. More...
bool	msgLvl (const MSG::Level lvl) const
	Test the output level. More...
MsgStream &	msg () const
	The standard message stream. More...
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream. More...
void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

Private Types
enum	NumberOfDimensions {   N_STATIONNAME = 8, OFFSET_STATIONNAME = 41, N_STATIONETA = 11, OFFSET_STATIONETA = -5,   N_STATIONPHI = 48, OFFSET_STATIONPHI = 1, N_GASGAP = 3, OFFSET_GASGAP = 1,   N_ISSTRIP = 2, OFFSET_ISSTRIP = 0, N_ABSSTATIONETA = 5, OFFSET_ABSSTATIONETA = 1,   N_STRIPCHANNEL = 32, OFFSET_STRIPCHANNEL = 1 }
enum	TgcStation { kOUTER = 0, kINNER, N_STATION }
enum	TgcSensor { kWIRE = 0, kSTRIP, N_SENSOR }

Private Member Functions
StatusCode	readFileOfTimeJitter ()
	Reads parameters for intrinsic time response from timejitter.dat. More...
float	timeJitter (const Amg::Vector3D &, CLHEP::HepRandomEngine *rndmEngine) const
	Calculates intrinsic time response according to incident angle of a track based on time response parameters. More...
bool	efficiencyCheck (const TgcSensor sensor, CLHEP::HepRandomEngine *rndmEngine) const
	Determines whether a hit is detected or not. More...
bool	efficiencyCheck (const std::string &stationName, const int stationEta, const int stationPhi, const int gasGap, const TgcSensor sensor, const double energyDeposit) const
uint16_t	bcTagging (const double digittime, const double window, const double offset) const
StatusCode	readFileOfEnergyThreshold ()
	Read share/TGC_Digitization_energyThreshold.dat file. More...
StatusCode	readFileOfDeadChamber ()
	Read share/TGC_Digitization_deadChamber.dat file. More...
StatusCode	readFileOfStripPosition ()
	Read share/TGC_Digitization_StripPosition.dat file. More...
double	getEnergyThreshold (const std::string &stationName, int stationEta, int stationPhi, int gasGap, const TgcSensor sensor) const
	Get energy threshold value for each chamber. More...
void	randomCrossTalk (const TgcDigitCrosstalkData *crosstalk, const Identifier elemId, const int gasGap, const TgcSensor sensor, const int channel, const float posInStrip, const float digitTime, const float time_offset, CLHEP::HepRandomEngine *rndmEngine, TgcDigitCollection *digits) const
bool	isDeadChamber (const std::string &stationName, int stationEta, int stationPhi, int gasGap)
	Method to check a chamber is dead or active. More...
float	getStripPosition (const std::string &stationName, int stationEta, int channel) const
	Method to get position of Strip channel. More...
float	getDistanceToAsdFromSensor (const TgcDigitASDposData *readCdo, const int iStationName, const int stationEta, const int stationPhi, const TgcSensor sensor, const int channel, const float position) const
	Method to get propagation time to the ASD from the sensor. More...
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Static Private Member Functions
static void	addDigit (const Identifier id, const uint16_t bctag, TgcDigitCollection *digits)
static int	getIStationName (const std::string &staionName)
	Get stationName integer from stationName string. More...
static double	getSigPropTimeDelay (const float cableDistance)
	Method to get signal propagation time delay. More...
static float	timeDiffByCableRadiusOfInner (const int iStationName, const int stationPhi, const int channel)
	Method to get time difference by cable radius of inner. More...
static float	getTimeOffset (const TgcDigitTimeOffsetData *readCdo, const uint16_t station_num, const int station_eta, const TgcSensor sensor)
	Method to get time offset to absorb signal delay. More...
static float	getCrosstalkProbability (const TgcDigitCrosstalkData *readCdo, const uint16_t layer_id, const TgcSensor sensor, const unsigned int index_prob)
	Method to get the channel crosstalk probability. More...

Private Attributes
double	m_energyThreshold [N_STATIONNAME][N_STATIONETA][N_STATIONPHI][N_GASGAP][N_ISSTRIP] {}
	Energy threshold value for each chamber. More...
bool	m_isDeadChamber [N_STATIONNAME][N_STATIONETA][N_STATIONPHI][N_GASGAP] {}
	Dead chamber flag for each chamber. More...
float	m_StripPos [N_STATIONNAME][N_ABSSTATIONETA][N_STRIPCHANNEL] {}
	Position of Strip Channel (Longer base or Shorter base) More...
std::vector< std::vector< float > >	m_vecAngle_Time
const TgcHitIdHelper *	m_hitIdHelper
unsigned int	m_runperiod
const MuonGM::MuonDetectorManager *	m_mdManager
const TgcIdHelper *	m_idHelper
float	m_efficiency [N_SENSOR] {}
double	m_gateTimeWindow [N_STATION][N_SENSOR] {}
	define the time windows for signals from wiregangs and strips. More...
double	m_bunchCrossingTime
bool	m_doFourBunchDigitization {false}
	Activate four bunch digitization. More...
std::string	m_nm
	Message source name. More...
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels) More...
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer. More...
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level. More...
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging) More...

Detailed Description

for SCT_ChargeTrappingTool

All functionality of TGC digitization is implemented to this class.

Definition at line 40 of file TgcDigitMaker.h.

Member Enumeration Documentation

NumberOfDimensions

enum TgcDigitMaker::NumberOfDimensions

private

Enumerator
N_STATIONNAME
OFFSET_STATIONNAME
N_STATIONETA
OFFSET_STATIONETA
N_STATIONPHI
OFFSET_STATIONPHI
N_GASGAP
OFFSET_GASGAP
N_ISSTRIP
OFFSET_ISSTRIP
N_ABSSTATIONETA
OFFSET_ABSSTATIONETA
N_STRIPCHANNEL
OFFSET_STRIPCHANNEL

Definition at line 85 of file TgcDigitMaker.h.

                            {
        N_STATIONNAME = 8,
        OFFSET_STATIONNAME = 41,
        N_STATIONETA = 11,
        OFFSET_STATIONETA = -5,
        N_STATIONPHI = 48,
        OFFSET_STATIONPHI = 1,
        N_GASGAP = 3,
        OFFSET_GASGAP = 1,
        N_ISSTRIP = 2,
        OFFSET_ISSTRIP = 0,
        N_ABSSTATIONETA = 5,
        OFFSET_ABSSTATIONETA = 1,
        N_STRIPCHANNEL = 32,
        OFFSET_STRIPCHANNEL = 1
    };

TgcSensor

enum TgcDigitMaker::TgcSensor

private

Enumerator
kWIRE
kSTRIP
N_SENSOR

Definition at line 102 of file TgcDigitMaker.h.

{ kWIRE = 0, kSTRIP, N_SENSOR };

TgcStation

enum TgcDigitMaker::TgcStation

private

Enumerator
kOUTER
kINNER
N_STATION

Definition at line 101 of file TgcDigitMaker.h.

{ kOUTER = 0, kINNER, N_STATION };

Constructor & Destructor Documentation

TgcDigitMaker()

TgcDigitMaker::TgcDigitMaker	(	const TgcHitIdHelper *	hitIdHelper,
		const MuonGM::MuonDetectorManager *	mdManager,
		unsigned int	runperiod,
		const bool	doFourBunch
	)

Definition at line 31 of file TgcDigitMaker.cxx.

    : AthMessaging("TgcDigitMaker"), m_doFourBunchDigitization(doFourBunch) {
    m_hitIdHelper = hitIdHelper;
    m_mdManager = mdManager;
    m_runperiod = runperiod;
    m_idHelper = nullptr;
    m_efficiency[kWIRE] = m_efficiency[kSTRIP] =
        1.000;  // 100% efficiency for TGCSimHit_p1
    m_gateTimeWindow[kOUTER][kWIRE] =
        29.32;  // 29.32ns = 26ns + 4  * 0.83ns(outer station)
    m_gateTimeWindow[kOUTER][kSTRIP] =
        40.94;  // 40.94ns = 26ns + 18 * 0.83ns(outer station)
    m_gateTimeWindow[kINNER][kWIRE] =
        33.47;  // 33.47ns = 26ns + 9  * 0.83ns(inner station)
    m_gateTimeWindow[kINNER][kSTRIP] =
        45.09;                    // 45.09ns = 26ns + 23 * 0.83ns(inner station)
    m_bunchCrossingTime = 24.95;  // 24.95 ns =(40.08 MHz)^(-1)
}

~TgcDigitMaker()

TgcDigitMaker::~TgcDigitMaker ( )

virtualdefault

Member Function Documentation

addDigit()

void TgcDigitMaker::addDigit ( const Identifier  id, const uint16_t  bctag, TgcDigitCollection *  digits  )

staticprivate

Definition at line 672 of file TgcDigitMaker.cxx.

                                                         {
    for (int bc = TgcDigit::BC_PREVIOUS; bc <= TgcDigit::BC_NEXTNEXT; bc++) {
        if ((bctag >> (bc - 1)) & 0x1) {
            bool duplicate = false;
            for (const auto digit : *digits) {
                if (id == digit->identify() && digit->bcTag() == bc) {
                    duplicate = true;
                    break;
                }
            }
            if (!duplicate) {
                std::unique_ptr<TgcDigit> multihitDigit =
                    std::make_unique<TgcDigit>(id, bc);
                digits->push_back(multihitDigit.release());
            }
        }
    }
    }

bcTagging()

uint16_t TgcDigitMaker::bcTagging ( const double  digittime, const double  window, const double  offset  ) const

private

Definition at line 648 of file TgcDigitMaker.cxx.

                                                             {
    const double calc_coll_time =
        digitTime - offset;  // calculated collision time
    uint16_t bctag = 0;
    if (-m_bunchCrossingTime < calc_coll_time &&
        calc_coll_time < window - m_bunchCrossingTime) {
        bctag |= 0x1;
    }
    if (0. < calc_coll_time && calc_coll_time < window) {
        bctag |= 0x2;
    }
    if (m_bunchCrossingTime < calc_coll_time &&
        calc_coll_time < window + m_bunchCrossingTime) {
        bctag |= 0x4;
    }
    if (2. * m_bunchCrossingTime < calc_coll_time &&
        calc_coll_time < window + 2. * m_bunchCrossingTime &&
        m_doFourBunchDigitization) {
        bctag |= 0x8;
    }
    return bctag;
}

efficiencyCheck() [1/2]

bool TgcDigitMaker::efficiencyCheck ( const std::string &  stationName, const int  stationEta, const int  stationPhi, const int  gasGap, const TgcSensor  sensor, const double  energyDeposit  ) const

private

Definition at line 638 of file TgcDigitMaker.cxx.

                                                                      {
    // If the energy deposit is equal to or greater than the threshold value of
    // the chamber, return true.
    return (energyDeposit >= getEnergyThreshold(stationName, stationEta,
                                                stationPhi, gasGap, sensor));
}

efficiencyCheck() [2/2]

bool TgcDigitMaker::efficiencyCheck ( const TgcSensor  sensor, CLHEP::HepRandomEngine *  rndmEngine  ) const

private

Determines whether a hit is detected or not.

Definition at line 629 of file TgcDigitMaker.cxx.

                                                                            {
    if (CLHEP::RandFlat::shoot(rndmEngine, 0.0, 1.0) < m_efficiency[sensor])
        return true;
    ATH_MSG_DEBUG("efficiencyCheck(): Hit removed for sensor= "
                  << sensor << "(0=WIRE,1=STRIP)");
    return false;
}

executeDigi()

TgcDigitCollection * TgcDigitMaker::executeDigi	(	const TGCSimHit *	hit,
		const double	globalHitTime,
		const TgcDigitASDposData *	ASDpos,
		const TgcDigitTimeOffsetData *	TOffset,
		const TgcDigitCrosstalkData *	Crosstalk,
		CLHEP::HepRandomEngine *	rndmEngine
	)

A single hit can be digitized in the two directions independently: R and phi directions.

The information is restored by MuonGeoModel and accessible from the class of TGCReadOutElement. As described the preceding section, R direction is digitized based on the wire ganging information in the database. Digits in phi direction are calculated based on the formula in which the structure of strips in TGC is well expressed for amdb_simrec.P.03. For amdb_simrec.Q or later, phi direction is digitized based on the parameter in the database as the same manner of those in R direction. The method determines the response time for digits which is commonly used for signals from wire gangs and strips, and signal propagation time along wires and strips. In case that response time is outside of the time window to be set, that hit is removed. This method also removes some hits based on the detection efficiency to be set.

Definition at line 87 of file TgcDigitMaker.cxx.

                                      {
    // timing constant parameters
    constexpr float sensor_propagation_time =
        3.3 * CLHEP::ns /
        CLHEP::m;  // Until MC10, 8.5*ns/m for wire, 8.7*ns/m for strip.
                   // Since MC11, 3.3*ns/m (the speed of light) is used
                   // from the Z->mumu data/MC comparison.
    constexpr float cable_propagation_time = 5.0 * CLHEP::ns / CLHEP::m;
 
    // position constant parameters
    constexpr float wire_pitch = 1.8 * CLHEP::mm;
    constexpr float zwidth_frame = 17. * CLHEP::mm;
 
    ATH_MSG_DEBUG("executeDigi() Got HIT Id.");
 
    int Id = hit->TGCid();
    std::string stationName = m_hitIdHelper->GetStationName(Id);
    int stationEta = m_hitIdHelper->GetStationEta(Id);
    int stationPhi = m_hitIdHelper->GetStationPhi(Id);
    int ilyr = m_hitIdHelper->GetGasGap(Id);
 
    // Check the chamber is dead or not.
    if (isDeadChamber(stationName, stationEta, stationPhi, ilyr))
        return nullptr;
 
    const Identifier elemId =
        m_idHelper->elementID(stationName, stationEta, stationPhi);
    ATH_MSG_DEBUG("executeDigi() - element identifier is: "
                  << m_idHelper->show_to_string(elemId));
 
    const MuonGM::TgcReadoutElement* tgcChamber =
        m_mdManager->getTgcReadoutElement(elemId);
    if (!tgcChamber) {
        ATH_MSG_WARNING("executeDigi() - no ReadoutElement found for "
                        << m_idHelper->show_to_string(elemId));
        return nullptr;
    }
 
    IdContext tgcContext = m_idHelper->module_context();
    IdentifierHash coll_hash;
    if (m_idHelper->get_hash(elemId, coll_hash, &tgcContext)) {
        ATH_MSG_WARNING("Unable to get TGC hash id from TGC Digit collection "
                        << "context begin_index = " << tgcContext.begin_index()
                        << " context end_index  = " << tgcContext.end_index()
                        << " the identifier is " << elemId);
        elemId.show();
    }
 
    std::unique_ptr<TgcDigitCollection> digits =
        std::make_unique<TgcDigitCollection>(elemId, coll_hash);
 
    const Amg::Vector3D centreChamber = tgcChamber->globalPosition();
    float height = tgcChamber->getRsize();
    float hmin = sqrt(pow(centreChamber.x(), 2) + pow(centreChamber.y(), 2)) -
                 height / 2.;
 
    // Direction cosine of incident angle of this track
    Amg::Vector3D direCos = hit->localDireCos();
 
    // local position
    Amg::Vector3D localPos = hit->localPosition();
 
    // Local z direction is global r direction.
    float distanceZ = 1.4 * CLHEP::mm / direCos[0] * direCos[2];
    float zLocal = localPos.z() + distanceZ;
 
    // Local y direction is global phi direction.
    float distanceY = 1.4 * CLHEP::mm / direCos[0] * direCos[1];
    // This ySign depends on the implementation of TGC geometry in G4 simulation
    // left-handed coordinate in A side(+z, stationEta>0)
    float ySign = (stationEta < 0) ? +1. : -1.;
    float yLocal = ySign * (localPos.y() + distanceY);
 
    // Time of flight correction for each chamber
    // the offset is set to 0 for ASD located at the postion where tof is
    // minimum in each chamber, i.e. the ASD at the smallest R in each chamber
    double tofCorrection =
        (sqrt(pow(hmin + zwidth_frame, 2) + pow(centreChamber.z(), 2)) /
         (299792458. * (CLHEP::m / CLHEP::s)));  // FIXME use CLHEP::c_light
 
    // bunch time
    float bunchTime = globalHitTime - tofCorrection;
 
    static const float jitterInitial = 9999.;
    float jitter = jitterInitial;  // calculated at wire group calculation and
                                   // used also for strip calculation
 
    int iStationName = getIStationName(stationName);
 
    double energyDeposit = hit->energyDeposit();  // Energy deposit in MeV
    // If TGCSimHit_p1 is used, energyDeposit is the default value of -9999.
    // If TGCSimHit_p2 is used, energyDeposit is equal to or greater than 0.
    // Therefore, the energyDeposit value can be used to distinguish
    // TGCSimHit_p1 and TGCSimHit_p2. For TGCSimHit_p1, old efficiency check
    // with only isStrip variable is used. For TGCSimHit_p2, new efficiency
    // check with chamber dependent energy threshold is used.
 
 
    if ((energyDeposit >= -1. &&
         efficiencyCheck(
             stationName, stationEta, stationPhi, ilyr, kWIRE,
             energyDeposit)) ||  // efficiency check for TGCSimHit_p2 at first,
        (energyDeposit < -1. &&
         efficiencyCheck(
             kWIRE, rndmEngine))) {  // Old efficiencyCheck for TGCSimHit_p1
 
        int iWireGroup[2];
        float posInWireGroup[2] = {0., 0.};
        for (int iPosition = 0; iPosition < 2; iPosition++) {
            int nWireOffset = (std::abs(stationEta) == 5 ||
                               stationName.compare(0, 2, "T4") == 0)
                                  ? 1
                                  : 0;
            // for chambers in which only the first wire  is not connected : 1
            // for chambers in which the first and last wires are not connected
            // OR all wires are connected : 0
 
            double zPosInSensArea =
                zLocal + static_cast<double>(tgcChamber->nWires(ilyr) -
                                             nWireOffset) *
                             wire_pitch / 2.;
 
            // check a hit in the sensitive area
            if (zPosInSensArea < 0. ||
                zPosInSensArea > tgcChamber->nWires(ilyr) * wire_pitch) {
                iWireGroup[iPosition] = 0;
                posInWireGroup[iPosition] = 0.;
                ATH_MSG_DEBUG(
                    "executeDigi(): Wire: Hit position located at outside of a "
                    "sensitive volume, "
                    << " id: " << stationName << "/" << stationEta << "/"
                    << stationPhi << "/" << ilyr
                    << " position: " << zPosInSensArea << " xlocal: " << zLocal
                    << " dir cosine: " << direCos[0] << "/" << direCos[1] << "/"
                    << direCos[2]
                    << " active region: " << height - zwidth_frame * 2.);
            } else {
                int igang = 1;
                int wire_index = 0;
                while (wire_pitch * (static_cast<float>(wire_index)) <
                           zPosInSensArea &&
                       igang <= tgcChamber->nWireGangs(ilyr)) {
                    wire_index += tgcChamber->nWires(ilyr, igang);
                    igang++;
                }
                posInWireGroup[iPosition] =
                    (zPosInSensArea / wire_pitch -
                     (static_cast<float>(wire_index))) /
                        (static_cast<float>(
                            tgcChamber->nWires(ilyr, igang - 1))) +
                    1.;
 
                iWireGroup[iPosition] = ((1 == igang) ? 1 : igang - 1);
            }
        }
 
        unsigned int jWG[2] = {
            (iWireGroup[0] <= iWireGroup[1]) ? (unsigned int)(0)
                                             : (unsigned int)(1),
            (iWireGroup[0] <= iWireGroup[1]) ? (unsigned int)(1)
                                             : (unsigned int)(0)};
        int iWG[2] = {iWireGroup[jWG[0]], iWireGroup[jWG[1]]};
        float posInWG[2] = {posInWireGroup[jWG[0]], posInWireGroup[jWG[1]]};
        if (iWG[0] != iWG[1]) {
            ATH_MSG_DEBUG("executeDigi(): Multihits found in WIRE GROUPs:"
                          << iWG[0] << " " << iWG[1]);
        }
 
        // === BC tagging from the hit timing ===
        for (int iwg = iWG[0]; iwg <= iWG[1]; iwg++) {
            if (1 <= iwg && iwg <= tgcChamber->nWireGangs(ilyr)) {
                // timing window offset
                float wire_timeOffset =
                    (TOffset != nullptr)
                        ? this->getTimeOffset(TOffset, iStationName, stationEta,
                                              kWIRE)
                        : 0.;
                // EI/FI has different offset between two ASDs.
                if (iStationName > 46)
                    wire_timeOffset +=
                        (iwg < 17) ? 0.5 * CLHEP::ns : -0.5 * CLHEP::ns;
 
                // TGC response time calculation
                float jit = timeJitter(direCos, rndmEngine);
                if (jit < jitter)
                    jitter = jit;
                float ySignPhi = (stationPhi % 2 == 1) ? -1. : +1.;
                // stationPhi%2==0 : +1. : ASD attached at the smaller phi side
                // of TGC stationPhi%2==1 : -1. : ASD attached at the larger phi
                // side of TGC
                float wTimeDiffByRadiusOfInner =
                    this->timeDiffByCableRadiusOfInner(iStationName, stationPhi,
                                                       iwg);
                double digit_time =
                    bunchTime + jitter + wTimeDiffByRadiusOfInner;
                float wASDDis{-1000.};
                if (ASDpos != nullptr) {
                    wASDDis = this->getDistanceToAsdFromSensor(
                        ASDpos, iStationName, abs(stationEta), stationPhi,
                        kWIRE, iwg, zLocal);
                    float wCableDis =
                        wASDDis + (ySignPhi * yLocal +
                                   tgcChamber->chamberWidth(zLocal) / 2.);
                    float wPropTime =
                        sensor_propagation_time *
                            (ySignPhi * yLocal +
                             tgcChamber->chamberWidth(zLocal) / 2.) +
                        cable_propagation_time * wASDDis;
                    digit_time +=
                        this->getSigPropTimeDelay(wCableDis) + wPropTime;
                }
 
                TgcStation station =
                    (m_idHelper->stationName(elemId) > 46) ? kINNER : kOUTER;
                uint16_t bctag =
                    bcTagging(digit_time, m_gateTimeWindow[station][kWIRE],
                              wire_timeOffset);
 
                if (bctag == 0x0) {
                    ATH_MSG_DEBUG(
                        "WireGroup: digitized time "
                        << digit_time << " ns is outside of time window from "
                        << wire_timeOffset << ". bunchTime: " << bunchTime
                        << " time jitter: " << jitter << " propagation time: "
                        << sensor_propagation_time *
                                   (ySignPhi * yLocal +
                                    tgcChamber->chamberWidth(zLocal) / 2.) +
                               cable_propagation_time * wASDDis
                        << " time difference by cable radius of inner station: "
                        << wTimeDiffByRadiusOfInner);
                } else {
                    Identifier newId = m_idHelper->channelID(
                        stationName, stationEta, stationPhi, ilyr, (int)kWIRE,
                        iwg);
                    addDigit(newId, bctag, digits.get());
 
                    if (iwg == iWG[0]) {
                        randomCrossTalk(Crosstalk, elemId, ilyr, kWIRE, iwg,
                                        posInWG[0], digit_time, wire_timeOffset,
                                        rndmEngine, digits.get());
                    }
 
                    ATH_MSG_DEBUG(
                        "WireGroup: newid breakdown digitTime x/y/z direcos "
                        "height_gang bctag: "
                        << newId << " " << stationName << "/" << stationEta
                        << "/" << stationPhi << "/" << ilyr << "/"
                        << "WIRE/" << iwg << " " << digit_time << " "
                        << localPos.x() << "/" << localPos.y() << "/"
                        << localPos.z() << " " << direCos.x() << "/"
                        << direCos.y() << "/" << direCos.z() << " " << height
                        << " " << tgcChamber->nWires(ilyr, iwg) << " "
                        << bctag);
                }
            } else {
                ATH_MSG_DEBUG(
                    "Wire Gang id is out of range. id, x/y/z, height: "
                    << iwg << " " << localPos.x() << "/" << localPos.y() << "/"
                    << localPos.z() << " " << height);
            }
        }
    }  // end of wire group calculation
 
    TgcSensor sensor = kSTRIP;
 
    if ((ilyr != 2 || (stationName.compare(0, 2, "T1") !=
                       0)) &&  // no stip in middle layers of T1*
        ((energyDeposit < -1. &&
          efficiencyCheck(
              sensor, rndmEngine)) ||  // Old efficiencyCheck for TGCSimHit_p1.
         (energyDeposit >= -1. &&
          efficiencyCheck(
              stationName, stationEta, stationPhi, ilyr, sensor,
              energyDeposit)))  // New efficiencyCheck for TGCSimHit_p2
    ) {
 
 
        int iStrip[2];
        float posInStrip[2] = {0., 0.};
        // Take into account of charge spread on cathod plane
        for (int iPosition = 0; iPosition < 2; iPosition++) {
 
            // check a hit in the sensitive area
            float zPos = zLocal + height / 2. - zwidth_frame;
            if (zPos < 0.) {
                iStrip[iPosition] = 0;
                posInStrip[iPosition] = 0.;
                ATH_MSG_DEBUG(
                    "Strip: Hit position located at outside of a sensitive "
                    "volume, id, position, xlocal0/1, dir cosine: "
                    << stationName << "/" << stationEta << "/" << stationPhi
                    << "/" << ilyr << zPos << " " << zLocal << " " << direCos[0]
                    << "/" << direCos[1] << "/" << direCos[2]);
            } else if (zPos > height - zwidth_frame * 2.) {
                iStrip[iPosition] = tgcChamber->nStrips(ilyr) + 1;
                posInStrip[iPosition] = 0.;
                ATH_MSG_DEBUG(
                    "Strip: Hit position located at outside of a sensitive "
                    "volume, id, position, active region: "
                    << stationName << "/" << stationEta << "/" << stationPhi
                    << "/" << ilyr << zPos << " "
                    << height - zwidth_frame * 2.);
            } else {  // sensitive area
 
                //
                // for layout P03
                //
                // number of strips in exclusive phi coverage of a chamber in
                // T[1-3] and T4
                const float dphi = tgcChamber->stripDeltaPhi();
                float phiLocal = atan2(yLocal, zLocal + height / 2. + hmin);
 
                ATH_MSG_DEBUG(
                    "dphi, phiLocal, yLocal, zLocall+ height/2.+hmin: "
                    << dphi << " " << phiLocal << " " << yLocal << " "
                    << zLocal + height / 2. + hmin);
 
                int istr = 0;
                if ((stationEta > 0 && ilyr == 1) ||
                    (stationEta < 0 && ilyr != 1)) {
                    istr = static_cast<int>(floor(phiLocal / dphi + 15.75)) + 1;
                    posInStrip[iPosition] =
                        phiLocal / dphi + 15.75 - static_cast<float>(istr - 1);
                    if (istr > 30) {  // treatment for two half strips
                        istr = static_cast<int>(floor(
                                   (phiLocal - dphi * 14.25) / (dphi / 2.))) +
                               31;
                        posInStrip[iPosition] =
                            (phiLocal - dphi * 14.25) / (dphi / 2.) -
                            static_cast<float>(istr - 31);
                    }
                } else {
                    istr = static_cast<int>(floor(phiLocal / dphi + 16.25)) + 1;
                    posInStrip[iPosition] =
                        phiLocal / dphi + 16.25 - static_cast<float>(istr - 1);
                    if (istr < 3) {  // treatment for two half strips
                        istr = static_cast<int>(floor(
                                   (phiLocal + dphi * 14.25) / (dphi / 2.))) +
                               3;
                        posInStrip[iPosition] =
                            (phiLocal + dphi * 14.25) / (dphi / 2.) -
                            static_cast<float>(istr - 3);
                    }
                }
                if (istr < 1) {
                    istr = 0;
                    posInStrip[iPosition] = 0.;
                } else if (32 < istr) {
                    istr = 33;
                    posInStrip[iPosition] = 0.;
                }
                iStrip[iPosition] = istr;
            }  // sensitive area
        }
 
        unsigned int jStr[2] = {
            (iStrip[0] <= iStrip[1]) ? (unsigned int)(0) : (unsigned int)(1),
            (iStrip[0] <= iStrip[1]) ? (unsigned int)(1) : (unsigned int)(0)};
        int iStr[2] = {iStrip[jStr[0]], iStrip[jStr[1]]};
        float posInStr[2] = {posInStrip[jStr[0]], posInStrip[jStr[1]]};
        if (iStr[0] != iStr[1]) {
            ATH_MSG_DEBUG("Multihits found in STRIPs:" << iStr[0] << " "
                                                       << iStr[1]);
        }
 
        // BC tagging from the timing
        float strip_timeOffset =
            (TOffset != nullptr)
                ? this->getTimeOffset(TOffset, iStationName, stationEta, kSTRIP)
                : 0.;
 
        for (int istr = iStr[0]; istr <= iStr[1]; istr++) {
            if (1 <= istr && istr <= 32) {
                // TGC response time calculation
                if (jitter > jitterInitial - 0.1) {
                    jitter = timeJitter(direCos, rndmEngine);
                }
                float zSignEta =
                    (abs(stationEta) % 2 == 1 && abs(stationEta) != 5) ? -1.
                                                                       : +1.;
                // if(abs(stationEta)%2 == 1 && abs(stationEta) != 5) : -1. :
                // ASD attached at the longer base of TGC else : +1. : ASD
                // attached at the shorter base of TGC
                float sTimeDiffByRadiusOfInner =
                    this->timeDiffByCableRadiusOfInner(iStationName, stationPhi,
                                                       istr);
                float sDigitTime =
                    bunchTime + jitter +
                    sensor_propagation_time *
                        (height / 2. + zwidth_frame + zSignEta * zLocal) +
                    sTimeDiffByRadiusOfInner;
                float sASDDis{-1000};
                if (ASDpos != nullptr) {
                    sASDDis = this->getDistanceToAsdFromSensor(
                        ASDpos, iStationName, abs(stationEta), stationPhi,
                        sensor, istr,
                        getStripPosition(stationName, abs(stationEta), istr));
                    float sCableDis = sASDDis + (height / 2. + zwidth_frame +
                                                 zSignEta * zLocal);
                    sDigitTime += this->getSigPropTimeDelay(sCableDis) +
                                  sASDDis * cable_propagation_time;
                }
 
                TgcStation station =
                    (m_idHelper->stationName(elemId) > 46) ? kINNER : kOUTER;
                uint16_t bctag =
                    bcTagging(sDigitTime, m_gateTimeWindow[station][sensor],
                              strip_timeOffset);
 
                if (bctag == 0x0) {
                    ATH_MSG_DEBUG(
                        "Strip: Digitized time is outside of time window. "
                        << sDigitTime << " bunchTime: " << bunchTime
                        << " time jitter: " << jitter << " propagation time: "
                        << sensor_propagation_time *
                               (height / 2. + zwidth_frame + zSignEta * zLocal)
                        << " time difference by cable radius of inner station: "
                        << sTimeDiffByRadiusOfInner);
                } else {
                    Identifier newId = m_idHelper->channelID(
                        stationName, stationEta, stationPhi, ilyr, (int)sensor,
                        istr);
                    addDigit(newId, bctag, digits.get());
 
                    if (istr == iStr[0]) {
                        randomCrossTalk(Crosstalk, elemId, ilyr, sensor,
                                        iStr[0], posInStr[0], sDigitTime,
                                        strip_timeOffset, rndmEngine,
                                        digits.get());
                    }
 
                    ATH_MSG_DEBUG(
                        "Strip: newid breakdown digitTime x/y/z direcos "
                        "r_center bctag: "
                        << newId << " " << stationName << "/" << stationEta
                        << "/" << stationPhi << "/" << ilyr << "/" << sensor
                        << "/" << istr << " " << sDigitTime << " "
                        << localPos.x() << "/" << localPos.y() << "/"
                        << localPos.z() << " " << direCos.x() << "/"
                        << direCos.y() << "/" << direCos.z() << " "
                        << height / 2. + hmin << " " << bctag);
                }
            } else {
                ATH_MSG_DEBUG("Strip id is out of range: gap, id: "
                              << ilyr << " " << istr);
            }
        }
    }  // end of strip number calculation
 
    return digits.release();
}

getCrosstalkProbability()

float TgcDigitMaker::getCrosstalkProbability ( const TgcDigitCrosstalkData *  readCdo, const uint16_t  layer_id, const TgcSensor  sensor, const unsigned int  index_prob  )

staticprivate

Method to get the channel crosstalk probability.

Definition at line 1183 of file TgcDigitMaker.cxx.

                                                           {
    if (readCdo == nullptr)
        return 0.;  // no crosstalk
    return ((sensor == TgcSensor::kSTRIP)
                ? readCdo->getStripProbability(layer_id, index_prob)
                : readCdo->getWireProbability(layer_id, index_prob));
}

getDistanceToAsdFromSensor()

float TgcDigitMaker::getDistanceToAsdFromSensor ( const TgcDigitASDposData *  readCdo, const int  iStationName, const int  stationEta, const int  stationPhi, const TgcSensor  sensor, const int  channel, const float  position  ) const

private

Method to get propagation time to the ASD from the sensor.

Definition at line 1144 of file TgcDigitMaker.cxx.

                                                   {
    // EIFI has different parameter for phi, BW is same parameter for phi (i.e.
    // -99 in DB).
    int phiId = (iStationName >= 47) ? stationPhi : 0x1f;
    uint16_t chamberId = (iStationName << 8) + (stationEta << 5) + phiId;
 
    unsigned int asdNo = static_cast<unsigned int>(channel - 1) /
                         TgcDigitASDposData::N_CHANNELINPUT_TOASD;
 
    float asd_position = 0.;
    const auto& map = (sensor == kSTRIP) ? readCdo->stripAsdPos : readCdo->wireAsdPos;
    auto it = map.find(chamberId);
 
    if (it != map.end()) {
        asd_position = it->second[asdNo] * CLHEP::m;
    } else {
        ATH_MSG_WARNING("Undefined chamberID is provided! station="
                        << iStationName << ", eta=" << stationEta
                        << ", phi=" << phiId);
    }
 
    float distance = fabs(position - asd_position);
    return distance;
}

getEnergyThreshold()

double TgcDigitMaker::getEnergyThreshold ( const std::string &  stationName, int  stationEta, int  stationPhi, int  gasGap, const TgcSensor  sensor  ) const

private

Get energy threshold value for each chamber.

Definition at line 925 of file TgcDigitMaker.cxx.

                                                                       {
    // Convert std::string stationName to int iStationName from 41 to 48
    int iStationName = getIStationName(stationName);
 
    // Subtract offsets to use these as the indices of the energyThreshold array
    iStationName -= OFFSET_STATIONNAME;
    stationEta -= OFFSET_STATIONETA;
    stationPhi -= OFFSET_STATIONPHI;
    gasGap -= OFFSET_GASGAP;
 
    double energyThreshold = +999999.;
 
    // If the indices are valid, the energyThreshold array is fetched.
    if ((iStationName >= 0 && iStationName < N_STATIONNAME) &&
        (stationEta >= 0 && stationEta < N_STATIONETA) &&
        (stationPhi >= 0 && stationPhi < N_STATIONPHI) &&
        (gasGap >= 0 && gasGap < N_GASGAP)) {
        energyThreshold = m_energyThreshold[iStationName][stationEta]
                                           [stationPhi][gasGap][sensor];
    }
 
    // Show the energy threshold value
    ATH_MSG_VERBOSE("getEnergyThreshold"
                    << " stationName= " << iStationName + OFFSET_STATIONNAME
                    << " stationEta= " << stationEta + OFFSET_STATIONETA
                    << " stationPhi= " << stationPhi + OFFSET_STATIONPHI
                    << " gasGap= " << gasGap + OFFSET_GASGAP << " sensor= "
                    << sensor << " energyThreshold(MeV)= " << energyThreshold);
 
    return energyThreshold;
}

getIStationName()

int TgcDigitMaker::getIStationName ( const std::string &  staionName )

staticprivate

Get stationName integer from stationName string.

Definition at line 1081 of file TgcDigitMaker.cxx.

                                                               {
    int iStationName = 0;
    if (stationName == "T1F")
        iStationName = 41;
    else if (stationName == "T1E")
        iStationName = 42;
    else if (stationName == "T2F")
        iStationName = 43;
    else if (stationName == "T2E")
        iStationName = 44;
    else if (stationName == "T3F")
        iStationName = 45;
    else if (stationName == "T3E")
        iStationName = 46;
    else if (stationName == "T4F")
        iStationName = 47;
    else if (stationName == "T4E")
        iStationName = 48;
 
    return iStationName;
}

getSigPropTimeDelay()

double TgcDigitMaker::getSigPropTimeDelay ( const float  cableDistance )

staticprivate

Method to get signal propagation time delay.

Definition at line 1138 of file TgcDigitMaker.cxx.

                                                                   {
    constexpr std::array<double, 2> par{0.0049 * CLHEP::ns, 0.0002 * CLHEP::ns};
    return par[0] * std::pow(cableDistance / CLHEP::m, 2) +
           par[1] * cableDistance / CLHEP::m;
}

getStripPosition()

float TgcDigitMaker::getStripPosition ( const std::string &  stationName, int  stationEta, int  channel  ) const

private

Method to get position of Strip channel.

Definition at line 1103 of file TgcDigitMaker.cxx.

                                                                         {
    // Convert std::string stationName to int iStationName from 41 to 48
    int iStationName = getIStationName(stationName);
 
    // Subtract offsets to use these as the indices of the energyThreshold array
    iStationName -= OFFSET_STATIONNAME;
    stationEta -= OFFSET_ABSSTATIONETA;
    channel -= OFFSET_STRIPCHANNEL;
 
    // Check the indices are valid
    if (iStationName < 0 || iStationName >= N_STATIONNAME)
        return 0.;
    if (stationEta < 0 || stationEta >= N_ABSSTATIONETA)
        return 0.;
    if (channel < 0 || channel >= N_STRIPCHANNEL)
        return 0.;
 
    return m_StripPos[iStationName][stationEta][channel];
}

getTimeOffset()

float TgcDigitMaker::getTimeOffset ( const TgcDigitTimeOffsetData *  readCdo, const uint16_t  station_num, const int  station_eta, const TgcSensor  sensor  )

staticprivate

Method to get time offset to absorb signal delay.

Definition at line 1172 of file TgcDigitMaker.cxx.

                                                           {
    uint16_t chamberId =
        (station_num << 3) + static_cast<uint16_t>(std::abs(station_eta));
    return ((sensor == TgcSensor::kSTRIP)
                ? readCdo->stripOffset.find(chamberId)->second
                : readCdo->wireOffset.find(chamberId)->second);
}

initialize()

StatusCode TgcDigitMaker::initialize

(

)

Initializes TgcHitIdHelper, TgcIdHelper and random number of a stream for the digitization.

In this method, readFileOfTimeJitter() reads and sets the parameters for intrinsic time response of TGC from timejitter.dat.

Definition at line 58 of file TgcDigitMaker.cxx.

                                     {
    // Initialize TgcIdHelper
    if (!m_hitIdHelper) {
        m_hitIdHelper = TgcHitIdHelper::GetHelper();
    }
 
    // initialize the TGC identifier helper
    m_idHelper = m_mdManager->tgcIdHelper();
 
    ATH_CHECK(readFileOfTimeJitter());
 
    // Read share/TGC_Digitization_energyThreshold.dat file and store values in
    // m_energyThreshold.
    ATH_CHECK(readFileOfEnergyThreshold());
 
    // Read share/TGC_Digitization_deadChamber.dat file and store values in
    // m_isDeadChamber.
    ATH_CHECK(readFileOfDeadChamber());
 
    // Read share/TGC_Digitization_StripPosition.dat file and store values in
    // m_StripPosition.
    ATH_CHECK(readFileOfStripPosition());
 
    return StatusCode::SUCCESS;
}

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  m_lvl = m_imsg ?
    static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
    MSG::INFO;
}

isDeadChamber()

bool TgcDigitMaker::isDeadChamber ( const std::string &  stationName, int  stationEta, int  stationPhi, int  gasGap  )

private

Method to check a chamber is dead or active.

Definition at line 1048 of file TgcDigitMaker.cxx.

                                                                              {
    bool v_isDeadChamber = true;
 
    // Convert std::string stationName to int iStationName from 41 to 48
    int iStationName = getIStationName(stationName);
 
    // Subtract offsets to use these as the indices of the energyThreshold array
    iStationName -= OFFSET_STATIONNAME;
    stationEta -= OFFSET_STATIONETA;
    stationPhi -= OFFSET_STATIONPHI;
    gasGap -= OFFSET_GASGAP;
 
    // If the indices are valid, the energyThreshold array is fetched.
    if ((iStationName >= 0 && iStationName < N_STATIONNAME) &&
        (stationEta >= 0 && stationEta < N_STATIONETA) &&
        (stationPhi >= 0 && stationPhi < N_STATIONPHI) &&
        (gasGap >= 0 && gasGap < N_GASGAP)) {
        v_isDeadChamber =
            m_isDeadChamber[iStationName][stationEta][stationPhi][gasGap];
    }
 
    // Show the energy threshold value
    ATH_MSG_VERBOSE("TgcDigitMaker::getEnergyThreshold"
                    << " stationName= " << iStationName + OFFSET_STATIONNAME
                    << " stationEta= " << stationEta + OFFSET_STATIONETA
                    << " stationPhi= " << stationPhi + OFFSET_STATIONPHI
                    << " gasGap= " << gasGap + OFFSET_GASGAP
                    << " isDeadChamber= " << v_isDeadChamber);
 
    return v_isDeadChamber;
}

msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 164 of file AthMessaging.h.

{
  MsgStream* ms = m_msg_tls.get();
  if (!ms) {
    if (!m_initialized.test_and_set()) initMessaging();
    ms = new MsgStream(m_imsg,m_nm);
    m_msg_tls.reset( ms );
  }
 
  ms->setLevel (m_lvl);
  return *ms;
}

msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level  lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 179 of file AthMessaging.h.

{ return msg() << lvl; }

msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level  lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

{
  if (!m_initialized.test_and_set()) initMessaging();
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

randomCrossTalk()

void TgcDigitMaker::randomCrossTalk ( const TgcDigitCrosstalkData *  crosstalk, const Identifier  elemId, const int  gasGap, const TgcSensor  sensor, const int  channel, const float  posInStrip, const float  digitTime, const float  time_offset, CLHEP::HepRandomEngine *  rndmEngine, TgcDigitCollection *  digits  ) const

private

Definition at line 961 of file TgcDigitMaker.cxx.

                                                                        {
    uint16_t station_number = m_idHelper->stationName(elemId);
    uint16_t station_eta = std::abs(m_idHelper->stationEta(elemId));
    uint16_t layer = gasGap;
    uint16_t layer_id = (station_number << 5) + (station_eta << 2) + layer;
 
    if (station_number < OFFSET_STATIONNAME ||
        station_number >= OFFSET_STATIONNAME + N_STATIONNAME ||
        station_eta <= 0 || station_eta > N_ABSSTATIONETA || layer <= 0 ||
        layer > N_GASGAP) {
        ATH_MSG_ERROR("Unexpected indices are provided!");
        return;
    }
 
    float prob1CrossTalk =
        this->getCrosstalkProbability(crosstalk, layer_id, sensor, 0);
    float prob11CrossTalk =
        this->getCrosstalkProbability(crosstalk, layer_id, sensor, 1);
    float prob20CrossTalk =
        this->getCrosstalkProbability(crosstalk, layer_id, sensor, 2);
    float prob21CrossTalk =
        this->getCrosstalkProbability(crosstalk, layer_id, sensor, 3);
 
    int nCrossTalks_neg = 0;
    int nCrossTalks_pos = 0;
 
    if (posInChan < prob1CrossTalk) {
        nCrossTalks_neg = 1;  // 1-0
    } else if (posInChan > 1. - prob1CrossTalk) {
        nCrossTalks_pos = 1;  // 0-1
    } else {
        double prob = CLHEP::RandFlat::shoot(rndmEngine, 0.0, 1.0);
        if (prob < prob11CrossTalk / (1. - 2. * prob1CrossTalk)) {
            nCrossTalks_neg = 1;
            nCrossTalks_pos = 1;  // 1-1
        } else if (prob < (prob20CrossTalk + prob11CrossTalk) /
                              (1. - 2. * prob1CrossTalk)) {
            if (posInChan < 0.5) {
                nCrossTalks_neg = 2;
            }  // 2-0
            else {
                nCrossTalks_pos = 2;
            }  // 0-2
        } else {
            if (prob <
                (prob20CrossTalk + prob11CrossTalk + 2. * prob21CrossTalk) /
                    (1. - 2. * prob1CrossTalk)) {
                if (posInChan < 0.5) {
                    nCrossTalks_neg = 2;
                    nCrossTalks_pos = 1;
                }  // 2-1
                else {
                    nCrossTalks_neg = 1;
                    nCrossTalks_pos = 2;
                }  // 1-2
            }
        }
    }
 
    if (nCrossTalks_neg == 0 && nCrossTalks_pos == 0)
        return;  // No cross-talk case
 
    // No time structure is implemented yet.
    float dt = digitTime;
    TgcStation station = (station_number > 46) ? kINNER : kOUTER;
    uint16_t bctag =
        bcTagging(dt, m_gateTimeWindow[station][sensor], time_offset);
    // obtain max channel number
    Identifier thisId =
        m_idHelper->channelID(elemId, gasGap, (int)sensor, channel);
    int maxChannelNumber = m_idHelper->channelMax(thisId);
 
    for (int jChan = channel - nCrossTalks_neg;
         jChan <= channel + nCrossTalks_pos; jChan++) {
        if (jChan == channel || jChan < 1 || jChan > maxChannelNumber)
            continue;
 
        Identifier newId =
            m_idHelper->channelID(elemId, gasGap, (int)sensor, jChan);
        addDigit(newId, bctag, digits);  // TgcDigit can be duplicated.
    }
}

readFileOfDeadChamber()

StatusCode TgcDigitMaker::readFileOfDeadChamber ( )

private

Read share/TGC_Digitization_deadChamber.dat file.

Definition at line 771 of file TgcDigitMaker.cxx.

                                                {
    // Indices to be used
    int iStationName, stationEta, stationPhi, gasGap;
 
    for (iStationName = 0; iStationName < N_STATIONNAME; iStationName++) {
        for (stationEta = 0; stationEta < N_STATIONETA; stationEta++) {
            for (stationPhi = 0; stationPhi < N_STATIONPHI; stationPhi++) {
                for (gasGap = 0; gasGap < N_GASGAP; gasGap++) {
                    m_isDeadChamber[iStationName][stationEta][stationPhi]
                                   [gasGap] = false;
                }
            }
        }
    }
 
    // Find path to the TGC_Digitization_deadChamber.dat file
    std::string fileName;
    if (m_runperiod == 1)
        fileName = "TGC_Digitization_deadChamber.dat";
    else if (m_runperiod == 2)
        fileName = "TGC_Digitization_2016deadChamber.dat";
    else if (m_runperiod == 3)
        fileName = "TGC_Digitization_NOdeadChamber.dat";
    else {
        ATH_MSG_ERROR("Run Period " << m_runperiod
                                    << " is unexpected in TgcDigitMaker - "
                                       "using NOdeadChamber configuration.");
        return StatusCode::FAILURE;
    }
    std::string fileWithPath = PathResolver::find_file(fileName, "DATAPATH");
    if (fileWithPath.empty()) {
        ATH_MSG_FATAL("readFileOfDeadChamber(): Could not find file "
                      << fileName);
        return StatusCode::FAILURE;
    }
 
    // Open the TGC_Digitization_deadChamber.dat file
    std::ifstream ifs;
    ifs.open(fileWithPath.c_str(), std::ios::in);
    if (ifs.bad()) {
        ATH_MSG_FATAL("readFileOfDeadChamber(): Could not open file "
                      << fileName);
        return StatusCode::FAILURE;
    }
 
    // Read the TGC_Digitization_deadChamber.dat file
    unsigned int nDeadChambers = 0;
    std::string comment;
    while (ifs.good()) {
        ifs >> iStationName >> stationEta >> stationPhi >> gasGap;
        bool valid = getline(ifs, comment).good();
        if (!valid)
            break;
 
        ATH_MSG_DEBUG("TgcDigitMaker::readFileOfDeadChamber"
                      << " stationName= " << iStationName << " stationEta= "
                      << stationEta << " stationPhi= " << stationPhi
                      << " gasGap= " << gasGap << " comment= " << comment);
 
        // Subtract offsets to use indices of isDeadChamber array
        iStationName -= OFFSET_STATIONNAME;
        stationEta -= OFFSET_STATIONETA;
        stationPhi -= OFFSET_STATIONPHI;
        gasGap -= OFFSET_GASGAP;
 
        // Check the indices are valid
        if (iStationName < 0 || iStationName >= N_STATIONNAME)
            continue;
        if (stationEta < 0 || stationEta >= N_STATIONETA)
            continue;
        if (stationPhi < 0 || stationPhi >= N_STATIONPHI)
            continue;
        if (gasGap < 0 || gasGap >= N_GASGAP)
            continue;
 
        m_isDeadChamber[iStationName][stationEta][stationPhi][gasGap] = true;
        nDeadChambers++;
 
        // If it is the end of the file, get out from while loop.
        if (ifs.eof())
            break;
    }
 
    // Close the TGC_Digitization_deadChamber.dat file
    ifs.close();
 
    ATH_MSG_INFO("readFileOfDeadChamber: the number of dead chambers = "
                 << nDeadChambers);
 
    return StatusCode::SUCCESS;
}

readFileOfEnergyThreshold()

StatusCode TgcDigitMaker::readFileOfEnergyThreshold ( )

private

Read share/TGC_Digitization_energyThreshold.dat file.

Definition at line 692 of file TgcDigitMaker.cxx.

                                                    {
    // Indices to be used
    int iStationName, stationEta, stationPhi, gasGap, isStrip;
 
    for (iStationName = 0; iStationName < N_STATIONNAME; iStationName++) {
        for (stationEta = 0; stationEta < N_STATIONETA; stationEta++) {
            for (stationPhi = 0; stationPhi < N_STATIONPHI; stationPhi++) {
                for (gasGap = 0; gasGap < N_GASGAP; gasGap++) {
                    for (isStrip = 0; isStrip < N_ISSTRIP; isStrip++) {
                        m_energyThreshold[iStationName][stationEta][stationPhi]
                                         [gasGap][isStrip] = -999999.;
                    }
                }
            }
        }
    }
 
    // Find path to the TGC_Digitization_energyThreshold.dat file
    const std::string fileName = "TGC_Digitization_energyThreshold.dat";
    std::string fileWithPath = PathResolver::find_file(fileName, "DATAPATH");
    if (fileWithPath.empty()) {
        ATH_MSG_FATAL("readFileOfEnergyThreshold(): Could not find file "
                      << fileName);
        return StatusCode::FAILURE;
    }
 
    // Open the TGC_Digitization_energyThreshold.dat file
    std::ifstream ifs;
    ifs.open(fileWithPath.c_str(), std::ios::in);
    if (ifs.bad()) {
        ATH_MSG_FATAL("readFileOfEnergyThreshold(): Could not open file "
                      << fileName);
        return StatusCode::FAILURE;
    }
 
    double energyThreshold;
    // Read the TGC_Digitization_energyThreshold.dat file
    while (ifs.good()) {
        ifs >> iStationName >> stationEta >> stationPhi >> gasGap >> isStrip >>
            energyThreshold;
        ATH_MSG_DEBUG("readFileOfEnergyThreshold"
                      << " stationName= " << iStationName << " stationEta= "
                      << stationEta << " stationPhi= " << stationPhi
                      << " gasGap= " << gasGap << " isStrip= " << isStrip
                      << " energyThreshold(MeV)= " << energyThreshold);
 
        // Subtract offsets to use indices of energyThreshold array
        iStationName -= OFFSET_STATIONNAME;
        stationEta -= OFFSET_STATIONETA;
        stationPhi -= OFFSET_STATIONPHI;
        gasGap -= OFFSET_GASGAP;
        isStrip -= OFFSET_ISSTRIP;
 
        // Check the indices are valid
        if (iStationName < 0 || iStationName >= N_STATIONNAME)
            continue;
        if (stationEta < 0 || stationEta >= N_STATIONETA)
            continue;
        if (stationPhi < 0 || stationPhi >= N_STATIONPHI)
            continue;
        if (gasGap < 0 || gasGap >= N_GASGAP)
            continue;
        if (isStrip < 0 || isStrip >= N_ISSTRIP)
            continue;
 
        m_energyThreshold[iStationName][stationEta][stationPhi][gasGap]
                         [isStrip] = energyThreshold;
 
        // If it is the end of the file, get out from while loop.
        if (ifs.eof())
            break;
    }
 
    // Close the TGC_Digitization_energyThreshold.dat file
    ifs.close();
 
    return StatusCode::SUCCESS;
}

readFileOfStripPosition()

StatusCode TgcDigitMaker::readFileOfStripPosition ( )

private

Read share/TGC_Digitization_StripPosition.dat file.

Definition at line 863 of file TgcDigitMaker.cxx.

                                                  {
    // Indices to be used
    int iStationName, stationEta, channel;
 
    for (iStationName = 0; iStationName < N_STATIONNAME; iStationName++) {
        for (stationEta = 0; stationEta < N_ABSSTATIONETA; stationEta++) {
            for (channel = 0; channel < N_STRIPCHANNEL; channel++) {
                m_StripPos[iStationName][stationEta][channel] = 0.;
            }
        }
    }
 
    // Find path to the TGC_Digitization_StripPosition.dat file
    const std::string fileName = "TGC_Digitization_StripPosition.dat";
    std::string fileWithPath = PathResolver::find_file(fileName, "DATAPATH");
    if (fileWithPath.empty()) {
        ATH_MSG_FATAL("readFileOfStripPosition(): Could not find file "
                      << fileName);
        return StatusCode::FAILURE;
    }
 
    // Open the TGC_Digitization_StripPosition.dat file
    std::ifstream ifs;
    ifs.open(fileWithPath.c_str(), std::ios::in);
    if (ifs.bad()) {
        ATH_MSG_FATAL("readFileOfStripPosition(): Could not open file "
                      << fileName);
        return StatusCode::FAILURE;
    }
 
    // Read the TGC_Digitization_StripPosition.dat file
    double strippos;
    while (ifs.good()) {
        ifs >> iStationName >> stationEta >> channel >> strippos;
        ATH_MSG_DEBUG("readFileOfStripPosition"
                      << " stationName= " << iStationName << " stationEta= "
                      << stationEta << " channel= " << channel
                      << " StripPosition= " << strippos);
 
        iStationName -= OFFSET_STATIONNAME;
        stationEta -= OFFSET_ABSSTATIONETA;
        channel -= OFFSET_STRIPCHANNEL;
 
        // Check the indices are valid
        if (iStationName < 0 || iStationName >= N_STATIONNAME)
            continue;
        if (stationEta < 0 || stationEta >= N_ABSSTATIONETA)
            continue;
        if (channel < 0 || channel >= N_STRIPCHANNEL)
            continue;
 
        m_StripPos[iStationName][stationEta][channel] = strippos;
        // If it is the end of the file, get out from while loop.
        if (ifs.eof())
            break;
    }
    // Close the TGC_Digitization_StripPosition.dat file
    ifs.close();
 
    return StatusCode::SUCCESS;
}

readFileOfTimeJitter()

StatusCode TgcDigitMaker::readFileOfTimeJitter ( )

private

Reads parameters for intrinsic time response from timejitter.dat.

Definition at line 547 of file TgcDigitMaker.cxx.

                                               {
    const char* const fileName = "TGC_Digitization_timejitter.dat";
    std::string fileWithPath = PathResolver::find_file(fileName, "DATAPATH");
 
    std::ifstream ifs;
    if (!fileWithPath.empty()) {
        ifs.open(fileWithPath.c_str(), std::ios::in);
    } else {
        ATH_MSG_FATAL("readFileOfTimeJitter(): Could not find file "
                      << fileName);
        return StatusCode::FAILURE;
    }
 
    if (ifs.bad()) {
        ATH_MSG_FATAL("readFileOfTimeJitter(): Could not open file "
                      << fileName);
        return StatusCode::FAILURE;
    }
 
    int angle = 0;
    int bins = 0;
    int i = 0;
    float prob = 0.;
    bool verbose = msgLvl(MSG::VERBOSE);
 
    while (ifs.good()) {
        ifs >> angle >> bins;
        if (ifs.eof())
            break;
        ATH_MSG_VERBOSE(
            "readFileOfTimeJitter(): Timejitter, angle, Number of bins, prob. "
            "dist.: "
            << angle << " " << bins << " ");
        m_vecAngle_Time.resize(i + 1);
        for (int j = 0; j < 41 /*bins*/; j++) {
            ifs >> prob;
            m_vecAngle_Time[i].push_back(prob);
            if (j == 0 && verbose)
                msg(MSG::VERBOSE) << "readFileOfTimeJitter(): ";
            if (verbose)
                msg(MSG::VERBOSE) << prob << " ";
        }
        if (verbose)
            msg(MSG::VERBOSE) << endmsg;
        i++;
    }
    ifs.close();
    return StatusCode::SUCCESS;
}

setLevel()

void AthMessaging::setLevel ( MSG::Level  lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

timeDiffByCableRadiusOfInner()

float TgcDigitMaker::timeDiffByCableRadiusOfInner ( const int  iStationName, const int  stationPhi, const int  channel  )

staticprivate

Method to get time difference by cable radius of inner.

Definition at line 1124 of file TgcDigitMaker.cxx.

                                                                     {
    float delay{0.};
    if (iStationName != 47 && iStationName != 48)
        return delay;  // Big Wheel has no delay for this.
 
    if (channel < 12 || (channel > 16 && channel < 27)) {
        int cablenum = (stationPhi >= 13) ? 25 - stationPhi : stationPhi;
        delay += 2.3 * CLHEP::ns - 0.06 * CLHEP::ns * float(cablenum);
    }
    return delay;
}

timeJitter()

float TgcDigitMaker::timeJitter ( const Amg::Vector3D &  direCosLocal, CLHEP::HepRandomEngine *  rndmEngine  ) const

private

Calculates intrinsic time response according to incident angle of a track based on time response parameters.

Definition at line 597 of file TgcDigitMaker.cxx.

                                                                        {
    float injectionAngle =
        atan2(fabs(direCosLocal[2]), fabs(direCosLocal[0])) / CLHEP::degree;
 
    int ithAngle = static_cast<int>(injectionAngle / 5.);
    float wAngle = injectionAngle / 5. - static_cast<float>(ithAngle);
    int jthAngle;
    if (ithAngle > 11) {
        ithAngle = 12;
        jthAngle = 12;
    } else {
        jthAngle = ithAngle + 1;
    }
 
    float jitter = 0.;
    float prob = 1.;
    float probRef = 0.;
 
    while (prob > probRef) {
        prob = CLHEP::RandFlat::shoot(rndmEngine, 0.0, 1.0);
        jitter = CLHEP::RandFlat::shoot(rndmEngine, 0.0, 1.0) * 40. *
                 CLHEP::ns;  // trial time jitter
        int ithJitter = static_cast<int>(jitter);
        // probability distribution calculated from weighted sum between
        // neighboring bins of angles
        probRef = (1. - wAngle) * m_vecAngle_Time[ithAngle][ithJitter] +
                  wAngle * m_vecAngle_Time[jthAngle][ithJitter];
    }
    return jitter;
}

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

m_bunchCrossingTime

double TgcDigitMaker::m_bunchCrossingTime

private

Definition at line 204 of file TgcDigitMaker.h.

m_doFourBunchDigitization

bool TgcDigitMaker::m_doFourBunchDigitization {false}

private

Activate four bunch digitization.

Definition at line 207 of file TgcDigitMaker.h.

m_efficiency

float TgcDigitMaker::m_efficiency[N_SENSOR] {}

private

Definition at line 195 of file TgcDigitMaker.h.

m_energyThreshold

double TgcDigitMaker::m_energyThreshold[N_STATIONNAME][N_STATIONETA][N_STATIONPHI][N_GASGAP][N_ISSTRIP] {}

private

Energy threshold value for each chamber.

Definition at line 180 of file TgcDigitMaker.h.

m_gateTimeWindow

double TgcDigitMaker::m_gateTimeWindow[N_STATION][N_SENSOR] {}

private

define the time windows for signals from wiregangs and strips.

The offsets are defined as relative time differences with respect to the time after TOF and cable length corrections. Bunch crossing time is specified.

Definition at line 203 of file TgcDigitMaker.h.

m_hitIdHelper

const TgcHitIdHelper* TgcDigitMaker::m_hitIdHelper

private

Definition at line 190 of file TgcDigitMaker.h.

m_idHelper

const TgcIdHelper* TgcDigitMaker::m_idHelper

private

Definition at line 194 of file TgcDigitMaker.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

m_isDeadChamber

bool TgcDigitMaker::m_isDeadChamber[N_STATIONNAME][N_STATIONETA][N_STATIONPHI][N_GASGAP] {}

private

Dead chamber flag for each chamber.

Definition at line 183 of file TgcDigitMaker.h.

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

m_mdManager

const MuonGM::MuonDetectorManager* TgcDigitMaker::m_mdManager

private

Definition at line 193 of file TgcDigitMaker.h.

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_runperiod

unsigned int TgcDigitMaker::m_runperiod

private

Definition at line 191 of file TgcDigitMaker.h.

m_StripPos

float TgcDigitMaker::m_StripPos[N_STATIONNAME][N_ABSSTATIONETA][N_STRIPCHANNEL] {}

private

Position of Strip Channel (Longer base or Shorter base)

Definition at line 186 of file TgcDigitMaker.h.

m_vecAngle_Time

std::vector<std::vector<float> > TgcDigitMaker::m_vecAngle_Time

private

Definition at line 188 of file TgcDigitMaker.h.

---

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

• TgcDigitMaker.h
• TgcDigitMaker.cxx