ScintillatorSimulator Class Reference

Digitize EKLMSim2Hits to get EKLM StripHits. More...

#include <ScintillatorSimulator.h>

Classes
struct	Photoelectron
	Photoelectron data. More...

Public Member Functions
	ScintillatorSimulator (const KLMScintillatorDigitizationParameters *digPar, ScintillatorFirmware *fitter, double digitizationInitialTime, bool debug)
	Constructor.
	ScintillatorSimulator (const ScintillatorSimulator &)=delete
	Copy constructor (disabled).
ScintillatorSimulator &	operator= (const ScintillatorSimulator &)=delete
	Operator = (disabled).
	~ScintillatorSimulator ()
	Destructor.
void	simulate (const std::multimap< KLMChannelNumber, const KLMSimHit * >::iterator &firstHit, const std::multimap< KLMChannelNumber, const KLMSimHit * >::iterator &end)
	Simulate a strip.
KLMScintillatorFirmwareFitResult *	getFPGAFit ()
	Get fit data.
enum ScintillatorFirmwareFitStatus	getFitStatus () const
	Get fit status.
double	getNPhotoelectrons ()
	Get number of photoelectrons (fit result).
int	getNGeneratedPhotoelectrons ()
	Get generated number of photoelectrons.
double	getEnergy ()
	Get total energy deposited in the strip (sum over ssimulation hits).
void	setFEEData (const KLMScintillatorFEEData *FEEData)
	Set FEE data.
void	generatePhotoelectrons (double stripLen, double distSiPM, int nPhotons, double timeShift, bool isReflected)
	Generate photoelectrons.
void	fillSiPMOutput (float *hist, bool useDirect, bool useReflected)
	Fill SiPM output.
float	getMCTime () const
	Get MC time.
float	getSiPMMCTime () const
	Get SiPM MC time.

Private Member Functions
void	reallocPhotoElectronBuffers (int size)
	Reallocate photoelectron buffers.
void	prepareSimulation ()
	Prepare simulation.
void	performSimulation ()
	Perform common simulation stage.
int *	sortPhotoelectrons (int nPhotoelectrons)
	Sort photoelectrons.
void	addRandomSiPMNoise ()
	Add random noise to the signal (amplitude-dependend).
void	simulateADC ()
	Simulate ADC (create digital signal from analog),.
void	debugOutput ()
	Debug output (signal and fit result histograms).

Private Attributes
const KLMTime *	m_Time
	Time.
const KLMScintillatorDigitizationParameters *	m_DigPar
	Parameters.
ScintillatorFirmware *	m_fitter
	Fitter.
double	m_DigitizationInitialTime
	Initial digitization time.
bool	m_Debug
	Debug mode (generates additional output files with histograms).
double	m_histRange
	Time range, (number of digitizations) * (ADC sampling time).
float *	m_amplitudeDirect
	Analog amplitude (direct).
float *	m_amplitudeReflected
	Analog amplitude (reflected).
float *	m_amplitude
	Analog amplitude.
int *	m_ADCAmplitude
	Digital amplitude.
double *	m_SignalTimeDependence
	Buffer for signal time dependence calculation.
double *	m_SignalTimeDependenceDiff
	Buffer for signal time dependence calculation.
struct Photoelectron *	m_Photoelectrons
	Buffer for photoelectron data.
int *	m_PhotoelectronIndex
	Buffer for photoelectron indices.
int *	m_PhotoelectronIndex2
	Buffer for photoelectron indices.
int	m_PhotoelectronBufferSize
	Size of photoelectron data buffer.
enum ScintillatorFirmwareFitStatus	m_FPGAStat
	FPGA fit status.
KLMScintillatorFirmwareFitResult	m_FPGAFit
	FPGA fit data.
int	m_npe
	Number of photoelectrons (generated).
double	m_Energy
	Total energy deposited in the strip.
std::string	m_stripName
	Name of the strip.
double	m_Pedestal
	Pedestal.
double	m_PhotoelectronAmplitude
	Photoelectron amplitude.
int	m_Threshold
	Threshold.
float	m_MCTime
	MC time.
float	m_SiPMMCTime
	MC time at SiPM.

Detailed Description

Digitize EKLMSim2Hits to get EKLM StripHits.

Definition at line 30 of file ScintillatorSimulator.h.

Constructor & Destructor Documentation

ScintillatorSimulator()

ScintillatorSimulator	(	const KLMScintillatorDigitizationParameters *	digPar,
		ScintillatorFirmware *	fitter,
		double	digitizationInitialTime,
		bool	debug
	)

Constructor.

Parameters

[in]	digPar	Digitization parameters.
[in]	digitizationInitialTime	Initial digitization time.
[in]	fitter	Fitter.
[in]	debug	Use debug mode.

Definition at line 62 of file ScintillatorSimulator.cc.

                                              :
  m_Time(&(KLMTime::Instance())),
  m_DigPar(digPar),
  m_fitter(fitter),
  m_DigitizationInitialTime(digitizationInitialTime),
  m_Debug(debug),
  m_FPGAStat(c_ScintillatorFirmwareNoSignal),
  m_npe(0),
  m_Energy(0),
  m_MCTime(-1),
  m_SiPMMCTime(-1)
{
  int i;
  const double samplingTime = m_DigPar->getADCSamplingTDCPeriods() *
                              m_Time->getTDCPeriod();
  /* cppcheck-suppress variableScope */
  double time, attenuationTime;
  m_histRange = m_DigPar->getNDigitizations() * samplingTime;
  m_Pedestal = m_DigPar->getADCPedestal();
  m_PhotoelectronAmplitude = m_DigPar->getADCPEAmplitude();
  m_Threshold = m_DigPar->getADCThreshold();
  /* Amplitude arrays. */
  m_amplitudeDirect = (float*)malloc(m_DigPar->getNDigitizations() *
                                     sizeof(float));
  if (m_amplitudeDirect == nullptr)
    B2FATAL(MemErr);
  m_amplitudeReflected = (float*)malloc(m_DigPar->getNDigitizations() *
                                        sizeof(float));
  if (m_amplitudeReflected == nullptr)
    B2FATAL(MemErr);
  m_amplitude = (float*)malloc(m_DigPar->getNDigitizations() * sizeof(float));
  if (m_amplitude == nullptr)
    B2FATAL(MemErr);
  m_ADCAmplitude = (int*)malloc(m_DigPar->getNDigitizations() * sizeof(int));
  if (m_ADCAmplitude == nullptr)
    B2FATAL(MemErr);
  m_SignalTimeDependence = (double*)malloc((m_DigPar->getNDigitizations() + 1) *
                                           sizeof(double));
  if (m_SignalTimeDependence == nullptr)
    B2FATAL(MemErr);
  m_SignalTimeDependenceDiff = (double*)malloc(m_DigPar->getNDigitizations() *
                                               sizeof(double));
  if (m_SignalTimeDependenceDiff == nullptr)
    B2FATAL(MemErr);
  attenuationTime = 1.0 / m_DigPar->getPEAttenuationFrequency();
  for (i = 0; i <= m_DigPar->getNDigitizations(); i++) {
    time = samplingTime * i;
    m_SignalTimeDependence[i] =
      exp(-digPar->getPEAttenuationFrequency() * time) * attenuationTime /
      samplingTime;
    if (i > 0) {
      m_SignalTimeDependenceDiff[i - 1] = m_SignalTimeDependence[i - 1] -
                                          m_SignalTimeDependence[i];
    }
  }
  m_Photoelectrons = nullptr;
  m_PhotoelectronIndex = nullptr;
  m_PhotoelectronIndex2 = nullptr;
  reallocPhotoElectronBuffers(100);
}

~ScintillatorSimulator()

~ScintillatorSimulator

(

)

Destructor.

Definition at line 127 of file ScintillatorSimulator.cc.

{
  free(m_amplitudeDirect);
  free(m_amplitudeReflected);
  free(m_amplitude);
  free(m_ADCAmplitude);
  free(m_SignalTimeDependence);
  free(m_SignalTimeDependenceDiff);
  free(m_Photoelectrons);
  free(m_PhotoelectronIndex);
  free(m_PhotoelectronIndex2);
}

Member Function Documentation

addRandomSiPMNoise()

void addRandomSiPMNoise ( )

private

Add random noise to the signal (amplitude-dependend).

Definition at line 252 of file ScintillatorSimulator.cc.

{
  int i;
  for (i = 0; i < m_DigPar->getNDigitizations(); i++)
    m_amplitude[i] = m_amplitude[i] +
                     gRandom->Poisson(m_DigPar->getMeanSiPMNoise());
}

debugOutput()

void debugOutput ( )

private

Debug output (signal and fit result histograms).

Definition at line 513 of file ScintillatorSimulator.cc.

{
  int i;
  std::string str;
  StoreObjPtr<EventMetaData> event;
  TFile* hfile = nullptr;
  TH1D* histAmplitudeDirect = nullptr;
  TH1D* histAmplitudeReflected = nullptr;
  TH1D* histAmplitude = nullptr;
  TH1D* histADCAmplitude = nullptr;
  try {
    histAmplitudeDirect =
      new TH1D("histAmplitudeDirect", m_stripName.c_str(),
               m_DigPar->getNDigitizations(), 0, m_histRange);
    histAmplitudeReflected =
      new TH1D("histAmplitudeReflected", m_stripName.c_str(),
               m_DigPar->getNDigitizations(), 0, m_histRange);
    histAmplitude =
      new TH1D("histAmplitude", m_stripName.c_str(),
               m_DigPar->getNDigitizations(), 0, m_histRange);
    histADCAmplitude =
      new TH1D("histADCAmplitude", m_stripName.c_str(),
               m_DigPar->getNDigitizations(), 0, m_histRange);
  } catch (std::bad_alloc& ba) {
    B2FATAL(MemErr);
  }
  for (i = 0; i < m_DigPar->getNDigitizations(); i++) {
    histAmplitudeDirect->SetBinContent(i + 1, m_amplitudeDirect[i]);
    histAmplitudeReflected->SetBinContent(i + 1, m_amplitudeReflected[i]);
    histAmplitude->SetBinContent(i + 1, m_amplitude[i]);
    histADCAmplitude->SetBinContent(i + 1, m_ADCAmplitude[i]);
  }
  str = std::string("experiment_") + std::to_string(event->getExperiment()) +
        "_run_" + std::to_string(event->getRun()) + "_event_" +
        std::to_string(event->getEvent()) + "_" + m_stripName + ".root";
  try {
    hfile = new TFile(str.c_str(), "NEW");
  } catch (std::bad_alloc& ba) {
    B2FATAL(MemErr);
  }
  hfile->Append(histAmplitudeDirect);
  hfile->Append(histAmplitudeReflected);
  hfile->Append(histAmplitude);
  hfile->Append(histADCAmplitude);
  hfile->Write();
  hfile->Close();
}

fillSiPMOutput()

void fillSiPMOutput	(	float *	hist,
		bool	useDirect,
		bool	useReflected
	)

Fill SiPM output.

Parameters

[in,out]	hist	Output histogram (signal is added to it).
[in]	useDirect	Use direct photons.
[in]	useReflected	Use reflected photons.

Definition at line 418 of file ScintillatorSimulator.cc.

{
  /* cppcheck-suppress variableScope */
  int i, bin, maxBin;
  double attenuationTime, sig, expSum;
  /* cppcheck-suppress variableScope */
  int ind1, ind2, ind3;
  int* indexArray;
  if (m_npe == 0)
    return;
  attenuationTime = 1.0 / m_DigPar->getPEAttenuationFrequency();
  indexArray = sortPhotoelectrons(m_npe);
  ind1 = 0;
  expSum = 0;
  while (1) {
    ind2 = ind1;
    bin = m_Photoelectrons[indexArray[ind1]].bin;
    while (1) {
      ind2++;
      if (ind2 == m_npe)
        break;
      if (bin != m_Photoelectrons[indexArray[ind2]].bin)
        break;
    }
    /* Now ind1 .. ind2 - photoelectrons in current bin. */
    for (ind3 = ind1; ind3 != ind2; ind3++) {
      if (m_Photoelectrons[indexArray[ind3]].isReflected && !useReflected)
        continue;
      if (!m_Photoelectrons[indexArray[ind3]].isReflected && !useDirect)
        continue;
      if (bin >= 0) {
        sig = attenuationTime - m_Photoelectrons[indexArray[ind3]].expTime *
              m_SignalTimeDependence[bin + 1];
        hist[bin] = hist[bin] + sig;
      }
      expSum = expSum + m_Photoelectrons[indexArray[ind3]].expTime;
    }
    if (ind2 == m_npe)
      maxBin = m_DigPar->getNDigitizations() - 1;
    else
      maxBin = m_Photoelectrons[indexArray[ind2]].bin;
    for (i = bin + 1; i <= maxBin; i++) {
      sig = m_SignalTimeDependenceDiff[i] * expSum;
      hist[i] = hist[i] + sig;
    }
    if (ind2 == m_npe)
      break;
    ind1 = ind2;
  }
}

generatePhotoelectrons()

void generatePhotoelectrons	(	double	stripLen,
		double	distSiPM,
		int	nPhotons,
		double	timeShift,
		bool	isReflected
	)

Generate photoelectrons.

Parameters

[in]	stripLen	Strip length.
[in]	distSiPM	Distance from hit to SiPM.
[in]	nPhotons	Number of photons to be simulated.
[in]	timeShift	Time of hit.
[in]	isReflected	Whether the hits are reflected or not.

Definition at line 308 of file ScintillatorSimulator.cc.

{
  const double samplingTime = m_DigPar->getADCSamplingTDCPeriods() *
                              m_Time->getTDCPeriod();
  const double maxHitTime = m_DigPar->getNDigitizations() * samplingTime;
  int i;
  /* cppcheck-suppress variableScope */
  double hitTime, deExcitationTime, cosTheta, hitDist, selection;
  double inverseLightSpeed, inverseAttenuationLength;
  inverseLightSpeed = 1.0 / m_DigPar->getFiberLightSpeed();
  inverseAttenuationLength = 1.0 / m_DigPar->getAttenuationLength();
  /* Generation of photoelectrons. */
  for (i = 0; i < nPhotons; i++) {
    if (m_npe >= m_PhotoelectronBufferSize)
      reallocPhotoElectronBuffers(m_npe + 100);
    cosTheta = gRandom->Uniform(m_DigPar->getMinCosTheta(), 1);
    if (!isReflected)
      hitDist = distSiPM / cosTheta;
    else
      hitDist = (2.0 * stripLen - distSiPM) / cosTheta;
    /* Fiber absorption. */
    selection = gRandom->Uniform();
    if (selection > exp(-hitDist * inverseAttenuationLength))
      continue;
    /* Account for mirror reflective index. */
    if (isReflected) {
      selection = gRandom->Uniform();
      if (selection > m_DigPar->getMirrorReflectiveIndex())
        continue;
    }
    deExcitationTime =
      gRandom->Exp(m_DigPar->getScintillatorDeExcitationTime()) +
      gRandom->Exp(m_DigPar->getFiberDeExcitationTime());
    hitTime = hitDist * inverseLightSpeed + deExcitationTime +
              timeShift - m_DigitizationInitialTime;
    if (hitTime >= maxHitTime)
      continue;
    if (hitTime >= 0)
      m_Photoelectrons[m_npe].bin = floor(hitTime / samplingTime);
    else
      m_Photoelectrons[m_npe].bin = -1;
    m_Photoelectrons[m_npe].expTime =
      exp(m_DigPar->getPEAttenuationFrequency() * hitTime);
    m_Photoelectrons[m_npe].isReflected = isReflected;
    m_PhotoelectronIndex[m_npe] = m_npe;
    m_npe++;
  }
}

getEnergy()

double getEnergy

(

)

Get total energy deposited in the strip (sum over ssimulation hits).

Definition at line 508 of file ScintillatorSimulator.cc.

{
  return m_Energy;
}

getFitStatus()

enum KLM::ScintillatorFirmwareFitStatus getFitStatus

(

)

const

Get fit status.

Returns

Status of the fit.

Definition at line 490 of file ScintillatorSimulator.cc.

{
  return m_FPGAStat;
}

getFPGAFit()

KLMScintillatorFirmwareFitResult * getFPGAFit

(

)

Get fit data.

Definition at line 485 of file ScintillatorSimulator.cc.

{
  return &m_FPGAFit;
}

getMCTime()

float getMCTime ( ) const

inline

Get MC time.

Returns

MC time.

Definition at line 132 of file ScintillatorSimulator.h.

      {
        return m_MCTime;
      }

getNGeneratedPhotoelectrons()

int getNGeneratedPhotoelectrons

(

)

Get generated number of photoelectrons.

Definition at line 503 of file ScintillatorSimulator.cc.

{
  return m_npe;
}

getNPhotoelectrons()

double getNPhotoelectrons

(

)

Get number of photoelectrons (fit result).

Definition at line 495 of file ScintillatorSimulator.cc.

{
  double intg;
  intg = m_FPGAFit.getAmplitude();
  return intg * m_DigPar->getPEAttenuationFrequency() /
         m_PhotoelectronAmplitude;
}

getSiPMMCTime()

float getSiPMMCTime ( ) const

inline

Get SiPM MC time.

Returns

SiPM MC yime.

Definition at line 141 of file ScintillatorSimulator.h.

      {
        return m_SiPMMCTime;
      }

performSimulation()

void performSimulation ( )

private

Perform common simulation stage.

Definition at line 230 of file ScintillatorSimulator.cc.

{
  if (m_Debug) {
    fillSiPMOutput(m_amplitudeDirect, true, false);
    fillSiPMOutput(m_amplitudeReflected, false, true);
    for (int i = 0; i < m_DigPar->getNDigitizations(); i++)
      m_amplitude[i] = m_amplitudeDirect[i] + m_amplitudeReflected[i];
  } else
    fillSiPMOutput(m_amplitude, true, true);
  /* SiPM noise and ADC. */
  if (m_DigPar->getMeanSiPMNoise() > 0)
    addRandomSiPMNoise();
  simulateADC();
  m_FPGAStat = m_fitter->fit(m_ADCAmplitude, m_Threshold, &m_FPGAFit);
  if (m_FPGAStat != c_ScintillatorFirmwareSuccessfulFit)
    return;
  if (m_Debug) {
    if (m_npe >= 10)
      debugOutput();
  }
}

prepareSimulation()

void prepareSimulation ( )

private

Prepare simulation.

Definition at line 148 of file ScintillatorSimulator.cc.

{
  m_MCTime = -1;
  m_SiPMMCTime = -1;
  m_npe = 0;
  m_Energy = 0;
  for (int i = 0; i < m_DigPar->getNDigitizations(); i++) {
    if (m_Debug) {
      m_amplitudeDirect[i] = 0;
      m_amplitudeReflected[i] = 0;
    } else
      m_amplitude[i] = 0;
  }
}

reallocPhotoElectronBuffers()

void reallocPhotoElectronBuffers ( int  size )

private

Reallocate photoelectron buffers.

Parameters

[in]

size

New size of buffers.

Definition at line 41 of file ScintillatorSimulator.cc.

{
  m_PhotoelectronBufferSize = size;
  /*
   * Here there is a memory leak in case of realloc() failure, but it does not
   * matter because a fatal error is issued in this case.
   */
  m_Photoelectrons =
    (struct Photoelectron*)realloc(m_Photoelectrons,
                                   size * sizeof(struct Photoelectron));
  m_PhotoelectronIndex = (int*)realloc(m_PhotoelectronIndex,
                                       size * sizeof(int));
  m_PhotoelectronIndex2 = (int*)realloc(m_PhotoelectronIndex2,
                                        size * sizeof(int));
  if (size != 0) {
    if (m_Photoelectrons == nullptr || m_PhotoelectronIndex == nullptr ||
        m_PhotoelectronIndex2 == nullptr)
      B2FATAL(MemErr);
  }
}

setFEEData()

void setFEEData

(

const KLMScintillatorFEEData *

FEEData

)

Set FEE data.

Definition at line 140 of file ScintillatorSimulator.cc.

{
  m_Pedestal = FEEData->getPedestal();
  m_PhotoelectronAmplitude = FEEData->getPhotoelectronAmplitude();
  m_Threshold = FEEData->getThreshold();
}

simulate()

void simulate	(	const std::multimap< KLMChannelNumber, const KLMSimHit * >::iterator &	firstHit,
		const std::multimap< KLMChannelNumber, const KLMSimHit * >::iterator &	end
	)

Simulate a strip.

Parameters

[in]	firstHit	First hit in this strip.
[in]	end	End of hit range.

Definition at line 163 of file ScintillatorSimulator.cc.

{
  m_stripName = "strip_" + std::to_string(firstHit->first);
  prepareSimulation();
  bklm::GeometryPar* geoPar = bklm::GeometryPar::instance();
  const KLMSimHit* hit = firstHit->second;
  double stripLength;
  bool isBKLM = (hit->getSubdetector() == KLMElementNumbers::c_BKLM);
  if (isBKLM) {
    const bklm::Module* module =
      geoPar->findModule(hit->getSection(), hit->getSector(), hit->getLayer());
    stripLength =
      2.0 * (hit->isPhiReadout() ?
             module->getPhiScintHalfLength(hit->getStrip()) :
             module->getZScintHalfLength(hit->getStrip()));
  } else {
    stripLength = EKLM::GeometryData::Instance().getStripLength(
                    hit->getStrip()) / CLHEP::mm * Unit::mm;
  }
  std::vector<const KLMSimHit*> hits;
  for (std::multimap<KLMChannelNumber, const KLMSimHit*>::iterator it = firstHit;
       it != end; ++it)
    hits.push_back(it->second);
  std::sort(hits.begin(), hits.end(), compareKLMSimHits);
  for (std::vector<const KLMSimHit*>::iterator it = hits.begin();
       it != hits.end(); ++it) {
    hit = *it;
    m_Energy = m_Energy + hit->getEnergyDeposit();
    /* Poisson mean for number of photons. */
    double nPhotons = hit->getEnergyDeposit() * m_DigPar->getNPEperMeV();
    /* Fill histograms. */
    double sipmDistance, time;
    if (isBKLM) {
      sipmDistance = hit->getPropagationTime() *
                     m_DigPar->getFiberLightSpeed();
      time = hit->getTime() + hit->getPropagationTime();
      if (m_MCTime < 0) {
        m_MCTime = hit->getTime();
        m_SiPMMCTime = time;
      } else {
        if (hit->getTime() < m_MCTime)
          m_MCTime = hit->getTime();
        if (time < m_SiPMMCTime)
          m_SiPMMCTime = time;
      }
    } else {
      sipmDistance = 0.5 * stripLength - hit->getLocalPosition().X();
      time = hit->getTime() + sipmDistance / m_DigPar->getFiberLightSpeed();
      if (m_MCTime < 0)
        m_MCTime = time;
      else
        m_MCTime = time < m_MCTime ? time : m_MCTime;
    }
    int generatedPhotons = gRandom->Poisson(nPhotons);
    generatePhotoelectrons(stripLength, sipmDistance, generatedPhotons,
                           hit->getTime(), false);
    if (m_DigPar->getMirrorReflectiveIndex() > 0) {
      generatedPhotons = gRandom->Poisson(nPhotons);
      generatePhotoelectrons(stripLength, sipmDistance, generatedPhotons,
                             hit->getTime(), true);
    }
  }
  performSimulation();
}

simulateADC()

void simulateADC ( )

private

Simulate ADC (create digital signal from analog),.

Definition at line 470 of file ScintillatorSimulator.cc.

{
  int i;
  /* cppcheck-suppress variableScope */
  double amp;
  if (m_Pedestal == 0 || m_PhotoelectronAmplitude == 0)
    B2FATAL("Incorrect EKLM ADC simulation parameters.");
  for (i = 0; i < m_DigPar->getNDigitizations(); i++) {
    amp = m_Pedestal - m_PhotoelectronAmplitude * m_amplitude[i];
    if (amp < m_DigPar->getADCSaturation())
      amp = m_DigPar->getADCSaturation();
    m_ADCAmplitude[i] = floor(amp);
  }
}

sortPhotoelectrons()

int * sortPhotoelectrons ( int  nPhotoelectrons )

private

Sort photoelectrons.

Parameters

[in]

nPhotoelectrons

Number of photoelectrons.

Returns

Pointer to index array.

Definition at line 260 of file ScintillatorSimulator.cc.

{
  int* currentIndexArray, *newIndexArray, *tmpIndexArray;
  int i, i1, i2, i1Max, i2Max, j, mergeSize;
  currentIndexArray = m_PhotoelectronIndex;
  newIndexArray = m_PhotoelectronIndex2;
  mergeSize = 1;
  while (mergeSize < nPhotoelectrons) {
    for (i = 0; i < nPhotoelectrons; i = i + 2 * mergeSize) {
      i1 = i;
      j = i;
      i2 = i + mergeSize;
      if (i2 > nPhotoelectrons)
        i2 = nPhotoelectrons;
      i1Max = i2;
      i2Max = i2 + mergeSize;
      if (i2Max > nPhotoelectrons)
        i2Max = nPhotoelectrons;
      while (i1 < i1Max || i2 < i2Max) {
        if (i1 < i1Max) {
          if (i2 < i2Max) {
            if (m_Photoelectrons[currentIndexArray[i1]].bin <
                m_Photoelectrons[currentIndexArray[i2]].bin) {
              newIndexArray[j] = currentIndexArray[i1];
              i1++;
            } else {
              newIndexArray[j] = currentIndexArray[i2];
              i2++;
            }
          } else {
            newIndexArray[j] = currentIndexArray[i1];
            i1++;
          }
        } else {
          newIndexArray[j] = currentIndexArray[i2];
          i2++;
        }
        j++;
      }
    }
    tmpIndexArray = currentIndexArray;
    currentIndexArray = newIndexArray;
    newIndexArray = tmpIndexArray;
    mergeSize = mergeSize * 2;
  }
  return currentIndexArray;
}

Member Data Documentation

m_ADCAmplitude

int* m_ADCAmplitude

private

Digital amplitude.

Definition at line 214 of file ScintillatorSimulator.h.

m_amplitude

float* m_amplitude

private

Analog amplitude.

Definition at line 211 of file ScintillatorSimulator.h.

m_amplitudeDirect

float* m_amplitudeDirect

private

Analog amplitude (direct).

Definition at line 205 of file ScintillatorSimulator.h.

m_amplitudeReflected

float* m_amplitudeReflected

private

Analog amplitude (reflected).

Definition at line 208 of file ScintillatorSimulator.h.

m_Debug

bool m_Debug

private

Debug mode (generates additional output files with histograms).

Definition at line 199 of file ScintillatorSimulator.h.

m_DigitizationInitialTime

double m_DigitizationInitialTime

private

Initial digitization time.

Definition at line 196 of file ScintillatorSimulator.h.

m_DigPar

const KLMScintillatorDigitizationParameters* m_DigPar

private

Parameters.

Definition at line 190 of file ScintillatorSimulator.h.

m_Energy

double m_Energy

private

Total energy deposited in the strip.

Definition at line 244 of file ScintillatorSimulator.h.

m_fitter

ScintillatorFirmware* m_fitter

private

Fitter.

Definition at line 193 of file ScintillatorSimulator.h.

m_FPGAFit

KLMScintillatorFirmwareFitResult m_FPGAFit

private

FPGA fit data.

Definition at line 238 of file ScintillatorSimulator.h.

m_FPGAStat

enum ScintillatorFirmwareFitStatus m_FPGAStat

private

FPGA fit status.

Definition at line 235 of file ScintillatorSimulator.h.

m_histRange

double m_histRange

private

Time range, (number of digitizations) * (ADC sampling time).

Definition at line 202 of file ScintillatorSimulator.h.

m_MCTime

float m_MCTime

private

MC time.

Definition at line 259 of file ScintillatorSimulator.h.

m_npe

int m_npe

private

Number of photoelectrons (generated).

Definition at line 241 of file ScintillatorSimulator.h.

m_Pedestal

double m_Pedestal

private

Pedestal.

Definition at line 250 of file ScintillatorSimulator.h.

m_PhotoelectronAmplitude

double m_PhotoelectronAmplitude

private

Photoelectron amplitude.

Definition at line 253 of file ScintillatorSimulator.h.

m_PhotoelectronBufferSize

int m_PhotoelectronBufferSize

private

Size of photoelectron data buffer.

Definition at line 232 of file ScintillatorSimulator.h.

m_PhotoelectronIndex

int* m_PhotoelectronIndex

private

Buffer for photoelectron indices.

Definition at line 226 of file ScintillatorSimulator.h.

m_PhotoelectronIndex2

int* m_PhotoelectronIndex2

private

Buffer for photoelectron indices.

Definition at line 229 of file ScintillatorSimulator.h.

m_Photoelectrons

struct Photoelectron* m_Photoelectrons

private

Buffer for photoelectron data.

Definition at line 223 of file ScintillatorSimulator.h.

m_SignalTimeDependence

double* m_SignalTimeDependence

private

Buffer for signal time dependence calculation.

Definition at line 217 of file ScintillatorSimulator.h.

m_SignalTimeDependenceDiff

double* m_SignalTimeDependenceDiff

private

Buffer for signal time dependence calculation.

Definition at line 220 of file ScintillatorSimulator.h.

m_SiPMMCTime

float m_SiPMMCTime

private

MC time at SiPM.

Definition at line 262 of file ScintillatorSimulator.h.

m_stripName

std::string m_stripName

private

Name of the strip.

Definition at line 247 of file ScintillatorSimulator.h.

m_Threshold

int m_Threshold

private

Threshold.

Definition at line 256 of file ScintillatorSimulator.h.

m_Time

const KLMTime* m_Time

private

Time.

Definition at line 187 of file ScintillatorSimulator.h.

---

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

• klm/simulation/include/ScintillatorSimulator.h
• klm/simulation/src/ScintillatorSimulator.cc