#include <LUCID_DigitizationToolBox.h>

Inheritance diagram for LUCID_DigitizationToolBox:

Collaboration diagram for LUCID_DigitizationToolBox:

Public Member Functions
	LUCID_DigitizationToolBox ()
	LUCID_DigitizationToolBox (int numTubes, double qdcChannelsPerPE, double qdcPedestal, double qdcFedNoiseFactor, double tdcPmtNoiseFactor, double tdcFedNoiseFactor, double TotalPMTgain, double AmpFactor, double Q1bin, int NoiseCharge, int numDyinodes, double dynodeGammaFactor, std::vector< double > pmtSmearing, std::vector< double > pmtScaling, std::vector< double > gasScaling, std::vector< double > npeThreshold, bool fillRootTree)
StatusCode	recordContainers (const ServiceHandle< StoreGateSvc > &, const std::string &)
StatusCode	fillDigitContainer (TimedHitCollection< LUCID_SimHit > &, CLHEP::HepRandomEngine *)
StatusCode	fillDigitContainer (LUCID_SimHitCollection , CLHEP::HepRandomEngine )
StatusCode	setDebugTree ()

Protected Member Functions
TTree *	getDebugTree ()
StatusCode	createAndStoreDigit (unsigned short tubeID, CLHEP::HepRandomEngine *rndEngine)
double	DynodeChainSimulation (double npe, CLHEP::HepRandomEngine *rndEngine) const

Static Protected Member Functions
static unsigned int	roundoff (double x)
static double	DynodeGainSmearing (double npe, CLHEP::HepRandomEngine *rndEngine)

Protected Attributes
LUCID_DigitContainer *	m_digitContainer {}
int	m_numTubes {}
double	m_qdcChannelsPerPE {}
double	m_qdcPedestal {}
double	m_qdcFedNoiseFactor {}
double	m_tdcPmtNoiseFactor {}
double	m_tdcFedNoiseFactor {}
double	m_TotalPMTgain {}
double	m_AmpFactor {}
double	m_Q1bin {}
int	m_NoiseCharge {}
int	m_numDyinodes {}
double	m_dynodeGammaFactor {}
double	m_gainPerDynode {}
double	m_ChargeToQdcFactor {}
std::vector< double >	m_pmtSmearing {}
std::vector< double >	m_pmtScaling {}
std::vector< double >	m_gasScaling {}
std::vector< double >	m_npeThreshold {}
bool	m_fillRootTree {}
unsigned short	m_tubeID {}
double	m_npe {}
unsigned short	m_npeGas {}
unsigned short	m_npePmt {}
unsigned short	m_qdc {}
double	m_tdcTot {}
unsigned short	m_tdc {}
TTree *	m_tubeInfo {}

Detailed Description

Definition at line 18 of file LUCID_DigitizationToolBox.h.

Constructor & Destructor Documentation

◆ LUCID_DigitizationToolBox() [1/2]

LUCID_DigitizationToolBox::LUCID_DigitizationToolBox ( )

Definition at line 19 of file LUCID_DigitizationToolBox.cxx.

                                                    : 
  m_tubeID(0),
  m_npe   (0),
  m_npeGas(0), 
  m_npePmt(0),
  m_qdc   (0),
  m_tdcTot(0),
  m_tdc   (0)
{
  m_digitContainer = nullptr;
  m_tubeInfo       = nullptr;
  
  m_numTubes          = 0;  
  m_qdcChannelsPerPE  = 0;
  m_qdcPedestal       = 0;
  m_qdcFedNoiseFactor = 0;
  m_tdcPmtNoiseFactor = 0;
  m_tdcFedNoiseFactor = 0;
 
  m_TotalPMTgain      = 0;
  m_AmpFactor         = 0;
  m_Q1bin             = 0;
  m_NoiseCharge       = 0;
  m_numDyinodes       = 0;
  m_dynodeGammaFactor = 0;
  m_gainPerDynode     = 0;
  m_ChargeToQdcFactor = 0;
 
  m_fillRootTree = 0;
}

◆ LUCID_DigitizationToolBox() [2/2]

LUCID_DigitizationToolBox::LUCID_DigitizationToolBox	(	int	numTubes,
		double	qdcChannelsPerPE,
		double	qdcPedestal,
		double	qdcFedNoiseFactor,
		double	tdcPmtNoiseFactor,
		double	tdcFedNoiseFactor,
		double	TotalPMTgain,
		double	AmpFactor,
		double	Q1bin,
		int	NoiseCharge,
		int	numDyinodes,
		double	dynodeGammaFactor,
		std::vector< double >	pmtSmearing,
		std::vector< double >	pmtScaling,
		std::vector< double >	gasScaling,
		std::vector< double >	npeThreshold,
		bool	fillRootTree )

Definition at line 51 of file LUCID_DigitizationToolBox.cxx.

                                                 : 
  m_numTubes            (numTubes),
  m_qdcChannelsPerPE    (qdcChannelsPerPE),
  m_qdcPedestal         (qdcPedestal),
  m_qdcFedNoiseFactor   (qdcFedNoiseFactor),
  m_tdcPmtNoiseFactor   (tdcPmtNoiseFactor),
  m_tdcFedNoiseFactor   (tdcFedNoiseFactor),
  m_TotalPMTgain        (TotalPMTgain),
  m_AmpFactor           (AmpFactor),
  m_Q1bin               (Q1bin),
  m_NoiseCharge         (NoiseCharge),
  m_numDyinodes         (numDyinodes),
  m_dynodeGammaFactor   (dynodeGammaFactor),
  m_pmtSmearing         (std::move(pmtSmearing)),
  m_pmtScaling          (std::move(pmtScaling)),
  m_gasScaling          (std::move(gasScaling)),
  m_npeThreshold        (std::move(npeThreshold)),
  m_fillRootTree        (fillRootTree),    
  m_tubeID              (0),
  m_npe                 (0),
  m_npeGas              (0), 
  m_npePmt              (0),
  m_qdc                 (0),
  m_tdcTot              (0),
  m_tdc                 (0)
{
  m_gainPerDynode     = std::pow(m_TotalPMTgain/m_dynodeGammaFactor, 1.0/static_cast<double>(m_numDyinodes)); 
  m_ChargeToQdcFactor = (double)1.6e-19 * m_AmpFactor/m_Q1bin;
 
  m_digitContainer = nullptr;
  m_tubeInfo       = nullptr;
}

Member Function Documentation

◆ createAndStoreDigit()

StatusCode LUCID_DigitizationToolBox::createAndStoreDigit	(	unsigned short	tubeID,
		CLHEP::HepRandomEngine *	rndEngine )

protected

Definition at line 133 of file LUCID_DigitizationToolBox.cxx.

                                                                                                                {
  
  double npePmt     = m_npePmt * m_pmtScaling[tubeID] * CLHEP::RandGaussQ::shoot(rndEngine, 1, m_pmtSmearing[tubeID]);
  double npeGas     = m_npeGas * m_gasScaling[tubeID];
  double npeTot     = npeGas + npePmt;
  double qdcSmeared = DynodeChainSimulation(npeTot, rndEngine);
    
  m_qdc = roundoff(qdcSmeared);
    
  m_npe = (m_qdc - m_qdcPedestal)/m_qdcChannelsPerPE;
 
  double tdcTotSmeared = 0;
    
  tdcTotSmeared = CLHEP::RandGaussQ::shoot(rndEngine, m_tdcTot     , m_tdcPmtNoiseFactor);
  tdcTotSmeared = CLHEP::RandGaussQ::shoot(rndEngine, tdcTotSmeared, m_tdcFedNoiseFactor);
  tdcTotSmeared = (npeTot) ? tdcTotSmeared/npeTot : tdcTotSmeared;
  tdcTotSmeared = (tdcTotSmeared < 0) ? 0 : tdcTotSmeared;
 
  m_tdc = static_cast<unsigned int>(tdcTotSmeared);
    
  if (m_fillRootTree) m_tubeInfo->Fill(); 
 
  m_digitContainer->push_back(new LUCID_Digit(m_tubeID, 
                          m_npe, 
                          m_npeGas, 
                          m_npePmt, 
                          m_qdc, 
                          m_tdc,
                          m_npe > m_npeThreshold[tubeID]));
  return StatusCode::SUCCESS;
}

◆ DynodeChainSimulation()

double LUCID_DigitizationToolBox::DynodeChainSimulation	(	double	npe,
		CLHEP::HepRandomEngine *	rndEngine ) const

protected

Definition at line 242 of file LUCID_DigitizationToolBox.cxx.

                                                                                                        {
  
  double smearedQDC = 0;
  
  if (npe == 0) return CLHEP::RandGaussQ::shoot(rndEngine, npe, m_ChargeToQdcFactor * m_NoiseCharge) + m_qdcPedestal;
  
  for (int i_pe=1; i_pe <= m_numDyinodes; i_pe++ ) {
    
    if (i_pe == 1) smearedQDC = DynodeGainSmearing(npe        * m_gainPerDynode * m_dynodeGammaFactor, rndEngine);  
    else           smearedQDC = DynodeGainSmearing(smearedQDC * m_gainPerDynode                      , rndEngine);  
  }
  
  smearedQDC = smearedQDC*m_ChargeToQdcFactor;
  
  return static_cast<double>(CLHEP::RandGaussQ::shoot(rndEngine, smearedQDC, m_NoiseCharge * m_ChargeToQdcFactor)) + m_qdcPedestal;  
}

◆ DynodeGainSmearing()

double LUCID_DigitizationToolBox::DynodeGainSmearing	(	double	npe,
		CLHEP::HepRandomEngine *	rndEngine )

staticprotected

Definition at line 260 of file LUCID_DigitizationToolBox.cxx.

                                                                                               {
  
  if (npe < 20) return CLHEP::RandPoissonQ::shoot(rndEngine, npe);
  else          return CLHEP::RandGaussQ::shoot  (rndEngine, npe, std::sqrt(npe));
}

◆ fillDigitContainer() [1/2]

StatusCode LUCID_DigitizationToolBox::fillDigitContainer	(	LUCID_SimHitCollection *	thpclucid,
		CLHEP::HepRandomEngine *	rndEngine )

Definition at line 199 of file LUCID_DigitizationToolBox.cxx.

                                                                                                                           {
      
  for (m_tubeID = 0; m_tubeID < m_numTubes; m_tubeID++) {
    
    if (m_tubeID >= 16 &&  m_tubeID <= 19) continue; // skip fiber tubes on sideA
    if (m_tubeID >= 36 &&  m_tubeID <= 39) continue; // skip fiber tubes on sideC
    
    m_npe = m_npeGas = m_npePmt = m_qdc = m_tdc = m_tdcTot = 0;
        
    LUCID_SimHitCollection::const_iterator hitIt  = thpclucid->begin();
    LUCID_SimHitCollection::const_iterator hitItE = thpclucid->end();
 
    for (; hitIt != hitItE; ++hitIt) {
      
      if (m_tubeID != (*hitIt).GetTubeID()) continue;
      
      if (!(*hitIt).isDetected(rndEngine)) continue; // applying quantum efficiency
 
      if ((*hitIt).GetGenVolume() == 1) m_npeGas++;
      if ((*hitIt).GetGenVolume() == 2) m_npeGas++;
      if ((*hitIt).GetGenVolume() == 3) m_npePmt++;
      
      m_tdcTot += (*hitIt).GetPreStepTime();
    }
 
    if (createAndStoreDigit(m_tubeID, rndEngine).isFailure()) return StatusCode::FAILURE;
  }
  
  return StatusCode::SUCCESS;
}

◆ fillDigitContainer() [2/2]

StatusCode LUCID_DigitizationToolBox::fillDigitContainer	(	TimedHitCollection< LUCID_SimHit > &	thpclucid,
		CLHEP::HepRandomEngine *	rndEngine )

Definition at line 166 of file LUCID_DigitizationToolBox.cxx.

                                                                                                                                     {
  
  for (m_tubeID = 0; m_tubeID < m_numTubes; m_tubeID++) {
    
    if (m_tubeID >= 16 &&  m_tubeID <= 19) continue; // skip fiber tubes on sideA
    if (m_tubeID >= 36 &&  m_tubeID <= 39) continue; // skip fiber tubes on sideC
    
    m_npe = m_npeGas = m_npePmt = m_qdc = m_tdc = m_tdcTot = 0;
    
    TimedHitCollection<LUCID_SimHit>                 thpc = thpclucid;
    TimedHitCollection<LUCID_SimHit>::const_iterator i, e, hitIt;
    
    while (thpc.nextDetectorElement(i, e)) for (hitIt = i; hitIt != e; ++hitIt) {
      
      if (m_tubeID != (*hitIt)->GetTubeID()) continue;
 
      if (!(*hitIt)->isDetected(rndEngine)) continue; // applying quantum efficiency
 
      if ((*hitIt)->GetGenVolume() == 1) m_npeGas++;
      if ((*hitIt)->GetGenVolume() == 2) m_npeGas++;
      if ((*hitIt)->GetGenVolume() == 3) m_npePmt++;
      
      m_tdcTot += (*hitIt)->GetPreStepTime();
    }
    
    if (createAndStoreDigit(m_tubeID, rndEngine).isFailure()) return StatusCode::FAILURE;
  }
  
  return StatusCode::SUCCESS;
 
}

◆ getDebugTree()

TTree * LUCID_DigitizationToolBox::getDebugTree ( )

inlineprotected

Definition at line 50 of file LUCID_DigitizationToolBox.h.

50{ return m_tubeInfo; }

◆ recordContainers()

StatusCode LUCID_DigitizationToolBox::recordContainers	(	const ServiceHandle< StoreGateSvc > &	digitsStore,
		const std::string &	key_digitCnt )

Definition at line 232 of file LUCID_DigitizationToolBox.cxx.

                                                                                                                                    {
  
  m_digitContainer = new LUCID_DigitContainer();
 
  StatusCode sc = digitsStore->record(m_digitContainer, key_digitCnt);
  
  return sc;
}

◆ roundoff()

unsigned int LUCID_DigitizationToolBox::roundoff ( double x )

staticprotected

Definition at line 122 of file LUCID_DigitizationToolBox.cxx.

                                                         {
  
  assert(x >= INT_MIN-0.5);
  assert(x <= INT_MAX+0.5);
  
  if (x >= 0) return static_cast<int>(x+0.5);
  
  return static_cast<int>(x-0.5);
}

◆ setDebugTree()

StatusCode LUCID_DigitizationToolBox::setDebugTree ( )

Definition at line 101 of file LUCID_DigitizationToolBox.cxx.

                                                   {
  SmartIF<ITHistSvc> digitHistSvc{Gaudi::svcLocator()->service("THistSvc")};
  if (!digitHistSvc) { return StatusCode::FAILURE; }
 
  m_tubeInfo = new TTree("t", "LUCID_LUMI_SUMMARY");
  
  StatusCode sc = digitHistSvc->regTree("/AANT/LUCID_LUMI_SUMMARY", m_tubeInfo);
  
  if (sc.isFailure()) return StatusCode::FAILURE;
  
  m_tubeInfo->Branch("tubeID", &m_tubeID, "tubeID/s");
  m_tubeInfo->Branch("npe"   , &m_npe   , "npe/D");
  m_tubeInfo->Branch("npeGas", &m_npeGas, "npeGas/s");
  m_tubeInfo->Branch("npePmt", &m_npePmt, "npePmt/s");
  m_tubeInfo->Branch("qdc"   , &m_qdc   , "qdc/s");
  m_tubeInfo->Branch("tdc"   , &m_tdc   , "tdc/s");
  
  return StatusCode::SUCCESS;
}

Member Data Documentation

◆ m_AmpFactor

double LUCID_DigitizationToolBox::m_AmpFactor {}

protected

Definition at line 69 of file LUCID_DigitizationToolBox.h.

69{};

◆ m_ChargeToQdcFactor

double LUCID_DigitizationToolBox::m_ChargeToQdcFactor {}

protected

Definition at line 76 of file LUCID_DigitizationToolBox.h.

76{};

◆ m_digitContainer

LUCID_DigitContainer* LUCID_DigitizationToolBox::m_digitContainer {}

protected

Definition at line 59 of file LUCID_DigitizationToolBox.h.

59{};

◆ m_dynodeGammaFactor

double LUCID_DigitizationToolBox::m_dynodeGammaFactor {}

protected

Definition at line 73 of file LUCID_DigitizationToolBox.h.

73{};

◆ m_fillRootTree

bool LUCID_DigitizationToolBox::m_fillRootTree {}

protected

Definition at line 83 of file LUCID_DigitizationToolBox.h.

83{};

◆ m_gainPerDynode

double LUCID_DigitizationToolBox::m_gainPerDynode {}

protected

Definition at line 75 of file LUCID_DigitizationToolBox.h.

75{};

◆ m_gasScaling

std::vector<double> LUCID_DigitizationToolBox::m_gasScaling {}

protected

Definition at line 80 of file LUCID_DigitizationToolBox.h.

80{};

◆ m_NoiseCharge

int LUCID_DigitizationToolBox::m_NoiseCharge {}

protected

Definition at line 71 of file LUCID_DigitizationToolBox.h.

71{};

◆ m_npe

double LUCID_DigitizationToolBox::m_npe {}

protected

Definition at line 86 of file LUCID_DigitizationToolBox.h.

86{};

◆ m_npeGas

unsigned short LUCID_DigitizationToolBox::m_npeGas {}

protected

Definition at line 87 of file LUCID_DigitizationToolBox.h.

87{};

◆ m_npePmt

unsigned short LUCID_DigitizationToolBox::m_npePmt {}

protected

Definition at line 88 of file LUCID_DigitizationToolBox.h.

88{};

◆ m_npeThreshold

std::vector<double> LUCID_DigitizationToolBox::m_npeThreshold {}

protected

Definition at line 81 of file LUCID_DigitizationToolBox.h.

81{};

◆ m_numDyinodes

int LUCID_DigitizationToolBox::m_numDyinodes {}

protected

Definition at line 72 of file LUCID_DigitizationToolBox.h.

72{};

◆ m_numTubes

int LUCID_DigitizationToolBox::m_numTubes {}

protected

Definition at line 61 of file LUCID_DigitizationToolBox.h.

61{};

◆ m_pmtScaling

std::vector<double> LUCID_DigitizationToolBox::m_pmtScaling {}

protected

Definition at line 79 of file LUCID_DigitizationToolBox.h.

79{};

◆ m_pmtSmearing

std::vector<double> LUCID_DigitizationToolBox::m_pmtSmearing {}

protected

Definition at line 78 of file LUCID_DigitizationToolBox.h.

78{};

◆ m_Q1bin

double LUCID_DigitizationToolBox::m_Q1bin {}

protected

Definition at line 70 of file LUCID_DigitizationToolBox.h.

70{};

◆ m_qdc

unsigned short LUCID_DigitizationToolBox::m_qdc {}

protected

Definition at line 89 of file LUCID_DigitizationToolBox.h.

89{};

◆ m_qdcChannelsPerPE

double LUCID_DigitizationToolBox::m_qdcChannelsPerPE {}

protected

Definition at line 62 of file LUCID_DigitizationToolBox.h.

62{};

◆ m_qdcFedNoiseFactor

double LUCID_DigitizationToolBox::m_qdcFedNoiseFactor {}

protected

Definition at line 64 of file LUCID_DigitizationToolBox.h.

64{};

◆ m_qdcPedestal

double LUCID_DigitizationToolBox::m_qdcPedestal {}

protected

Definition at line 63 of file LUCID_DigitizationToolBox.h.

63{};

◆ m_tdc

unsigned short LUCID_DigitizationToolBox::m_tdc {}

protected

Definition at line 91 of file LUCID_DigitizationToolBox.h.

91{};

◆ m_tdcFedNoiseFactor

double LUCID_DigitizationToolBox::m_tdcFedNoiseFactor {}

protected

Definition at line 66 of file LUCID_DigitizationToolBox.h.

66{};

◆ m_tdcPmtNoiseFactor

double LUCID_DigitizationToolBox::m_tdcPmtNoiseFactor {}

protected

Definition at line 65 of file LUCID_DigitizationToolBox.h.

65{};

◆ m_tdcTot

double LUCID_DigitizationToolBox::m_tdcTot {}

protected

Definition at line 90 of file LUCID_DigitizationToolBox.h.

90{};

◆ m_TotalPMTgain

double LUCID_DigitizationToolBox::m_TotalPMTgain {}

protected

Definition at line 68 of file LUCID_DigitizationToolBox.h.

68{};

◆ m_tubeID

unsigned short LUCID_DigitizationToolBox::m_tubeID {}

protected

Definition at line 85 of file LUCID_DigitizationToolBox.h.

85{};

◆ m_tubeInfo

TTree* LUCID_DigitizationToolBox::m_tubeInfo {}

protected

Definition at line 93 of file LUCID_DigitizationToolBox.h.

93{};

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

Public Member Functions

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ LUCID_DigitizationToolBox() [1/2]

◆ LUCID_DigitizationToolBox() [2/2]

Member Function Documentation

◆ createAndStoreDigit()

◆ DynodeChainSimulation()

◆ DynodeGainSmearing()

◆ fillDigitContainer() [1/2]

◆ fillDigitContainer() [2/2]

◆ getDebugTree()

◆ recordContainers()

◆ roundoff()

◆ setDebugTree()

Member Data Documentation

◆ m_AmpFactor

◆ m_ChargeToQdcFactor

◆ m_digitContainer

◆ m_dynodeGammaFactor

◆ m_fillRootTree

◆ m_gainPerDynode

◆ m_gasScaling

◆ m_NoiseCharge

◆ m_npe

◆ m_npeGas

◆ m_npePmt

◆ m_npeThreshold

◆ m_numDyinodes

◆ m_numTubes

◆ m_pmtScaling

◆ m_pmtSmearing

◆ m_Q1bin

◆ m_qdc

◆ m_qdcChannelsPerPE

◆ m_qdcFedNoiseFactor

◆ m_qdcPedestal

◆ m_tdc

◆ m_tdcFedNoiseFactor

◆ m_tdcPmtNoiseFactor

◆ m_tdcTot

◆ m_TotalPMTgain

◆ m_tubeID

◆ m_tubeInfo