EcalEBPhase2TimeReconstructor Class Reference

#include <EcalEBPhase2TimeReconstructor.h>

Public Member Functions
	EcalEBPhase2TimeReconstructor (bool debug)
virtual void	process (std::vector< int > &addout, std::vector< int > &ampRecoOutput, std::vector< int64_t > &output)
void	setParameters (uint32_t raw, const EcalEBPhase2TPGTimeWeightIdMap *ecaltpgTimeWeightMap, const EcalTPGWeightGroup *ecaltpgWeightGroup)
virtual	~EcalEBPhase2TimeReconstructor ()

Private Member Functions
void	process ()
int	setInput (int input)

Private Attributes
uint64_t	ampIn_ [maxSamplesUsed_]
int	buffer_ [maxSamplesUsed_]
bool	debug_
bool	extraShift_ [2] = {false, false}
int	inputsAlreadyIn_
uint64_t	invAmpAr_ [512]
int	processedOutput_
int	shift_
int	weights_ [maxSamplesUsed_]

Static Private Attributes
static const int	maxSamplesUsed_ = 12

Detailed Description

Author

L. Lutton, N. Marinelli - Univ. of Notre Dame Description: forPhase II Measures the timing of a xTal signal

Definition at line 16 of file EcalEBPhase2TimeReconstructor.h.

Constructor & Destructor Documentation

EcalEBPhase2TimeReconstructor()

EcalEBPhase2TimeReconstructor::EcalEBPhase2TimeReconstructor

(

bool

debug

)

Definition at line 9 of file EcalEBPhase2TimeReconstructor.cc.

    : debug_(debug), inputsAlreadyIn_(0), shift_(maxSamplesUsed_) {}

~EcalEBPhase2TimeReconstructor()

EcalEBPhase2TimeReconstructor::~EcalEBPhase2TimeReconstructor ( )

virtual

Definition at line 12 of file EcalEBPhase2TimeReconstructor.cc.

{}

Member Function Documentation

process() [1/2]

void EcalEBPhase2TimeReconstructor::process ( )

private

Definition at line 86 of file EcalEBPhase2TimeReconstructor.cc.

References ampIn_, buffer_, debug_, mps_fire::i, inputsAlreadyIn_, invAmpAr_, LogDebug, convertSQLitetoXML_cfg::output, processedOutput_, shift_, and weights_.

Referenced by process(), and EcalEBPhase2TrigPrimAlgo::run().

                                            {
  //UB FIXME: 5
  processedOutput_ = 0;
  if (inputsAlreadyIn_ < 12)
    return;
  int64_t output = 0;
  for (int i = 0; i < 12; i++) {
    output += (weights_[i] * buffer_[i]);
    if (debug_)
      LogDebug("") << " TimeFilter buffer " << buffer_[i] << " weight " << weights_[i] << " output " << output
                   << std::endl;
  }
  output = output >> shift_;
  if (debug_)
    LogDebug("") << " TimeFilter local  output " << output << std::endl;
  //Dividing output by the result of the amplitude reconstruction via an approximation using the invAmpAr lookup table
  int ampInd = 0;
  if (debug_)
    LogDebug("") << " inputsAlreadyIn_ " << inputsAlreadyIn_ << std::endl;
  if (inputsAlreadyIn_ > 12) {
    ampInd = 1;
  }

  if (debug_)
    LogDebug("") << " Begininning Final TimeFilter Calculation" << std::endl;

  int64_t tmpOutput = output * invAmpAr_[ampIn_[ampInd]];
  if (debug_)
    LogDebug("") << " output*tmpInvAmpAr " << tmpOutput << std::endl;

  output = tmpOutput >> 20;
  if (debug_)
    LogDebug("") << " output after bit shift " << output << std::endl;

  if (output < -1024)
    output = -1023;
  else if (output > 1024)
    output = 1023;
  if (debug_)
    LogDebug("") << " output after if/else " << output << std::endl;
  processedOutput_ = output;

  if (debug_)
    LogDebug("") << " TimeFilter final output " << processedOutput_ << std::endl;
}

process() [2/2]

void EcalEBPhase2TimeReconstructor::process ( std::vector< int > &  addout, std::vector< int > &  ampRecoOutput, std::vector< int64_t > &  output  )

virtual

Definition at line 37 of file EcalEBPhase2TimeReconstructor.cc.

References ampIn_, buffer_, debug_, mps_fire::i, inputsAlreadyIn_, dqmiolumiharvest::j, LogDebug, maxSamplesUsed_, process(), processedOutput_, ecalPh2::sampleSize, and setInput().

                                                                        {
  inputsAlreadyIn_ = 0;
  for (unsigned int i = 0; i < maxSamplesUsed_; i++) {
    buffer_[i] = 0;
  }

  //Taking in the results of the amplitude reconstruction
  //Bit shifting them for use as index of invAmpAr_ lookup table
  // move input amplitude (13 bits) to 9 bits to use as array index

  ampIn_[0] = ampRecoOutput[0] >> 4;
  ampIn_[1] = ampRecoOutput[1] >> 4;

  for (unsigned int i = 0; i < addout.size(); i++) {
    setInput(addout[i]);

    if (debug_) {
      LogDebug("") << "  EcalEBPhase2TimeReconstructor::process(std::vector<int> buffer_ " << std::endl;

      for (unsigned int j = 0; j < maxSamplesUsed_; j++) {
        LogDebug("") << " buffer_ " << buffer_[j];
      }
      LogDebug("") << "  " << std::endl;
    }

    if (i == (maxSamplesUsed_ - 1)) {
      if (debug_)
        LogDebug("") << "  EcalEBPhase2TimeReconstructor::process(std::vector<int>)    i = 11 " << std::endl;
      process();
      if (debug_)
        LogDebug("") << "  EcalEBPhase2TimeReconstructor::process(std::vector<int>)    after process() "
                     << processedOutput_ << std::endl;
      output[0] = processedOutput_;
      if (debug_)
        LogDebug("") << "  EcalEBPhase2TimeReconstructor::process(std::vector<int>)    after setting the output "
                     << output[0] << std::endl;
    } else if (i == (ecalPh2::sampleSize - 1)) {
      if (debug_)
        LogDebug("") << "  EcalEBPhase2TimeReconstructor::process(std::vector<int>)    i = 15 " << std::endl;
      process();
      output[1] = processedOutput_;
    }
  }

  return;
}

setInput()

int EcalEBPhase2TimeReconstructor::setInput ( int  input )

private

Definition at line 14 of file EcalEBPhase2TimeReconstructor.cc.

References buffer_, debug_, mps_fire::i, input, inputsAlreadyIn_, LogDebug, and maxSamplesUsed_.

Referenced by process().

                                                     {
  if (input > 0X7FFF) {
    edm::LogError("EcalEBPhase2TimeReconstructor::setInput") << "ERROR IN INPUT OF TIME FILTER" << std::endl;
    return -1;
  }
  if (inputsAlreadyIn_ < maxSamplesUsed_) {
    if (debug_)
      LogDebug("") << " EcalEBPhase2TimeReconstructor::setInput inputsAlreadyIn_<5 input " << input << std::endl;

    buffer_[inputsAlreadyIn_] = input;
    inputsAlreadyIn_++;
  } else {
    for (int i = 0; i < (maxSamplesUsed_ - 1); i++) {
      buffer_[i] = buffer_[i + 1];
      if (debug_)
        LogDebug("") << " EcalEBPhase2TimeReconstructor::setInput inputsAlreadyIn buffer " << buffer_[i] << std::endl;
    }
    buffer_[maxSamplesUsed_ - 1] = input;
    inputsAlreadyIn_++;
  }
  return 1;
}

setParameters()

void EcalEBPhase2TimeReconstructor::setParameters	(	uint32_t	raw,
		const EcalEBPhase2TPGTimeWeightIdMap *	ecaltpgTimeWeightMap,
		const EcalTPGWeightGroup *	ecaltpgWeightGroup
	)

Definition at line 132 of file EcalEBPhase2TimeReconstructor.cc.

References debug_, EcalEBPhase2TPGTimeWeightIdMap::getMap(), EcalTPGGroups::getMap(), mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, LogDebug, and weights_.

Referenced by EcalEBPhase2TrigPrimAlgo::run().

                                                                                                {
  uint32_t params_[12];
  const EcalTPGGroups::EcalTPGGroupsMap &groupmap = ecaltpgWeightGroup->getMap();
  if (debug_)
    LogDebug("") << " EcalEBPhase2TimeReconstructor::setParameters groupmap size " << groupmap.size() << std::endl;
  EcalTPGGroups::EcalTPGGroupsMapItr it = groupmap.find(raw);
  if (it != groupmap.end()) {
    uint32_t weightid = (*it).second;
    const EcalEBPhase2TPGTimeWeightIdMap::EcalEBPhase2TPGTimeWeightMap &weightmap = ecaltpgWeightMap->getMap();
    EcalEBPhase2TPGTimeWeightIdMap::EcalEBPhase2TPGTimeWeightMapItr itw = weightmap.find(weightid);

    (*itw).second.getValues(params_[0],
                            params_[1],
                            params_[2],
                            params_[3],
                            params_[4],
                            params_[5],
                            params_[6],
                            params_[7],
                            params_[8],
                            params_[9],
                            params_[10],
                            params_[11]);

    if (debug_)
      LogDebug("") << " EcalEBPhase2TimeReconstructor::setParameters time weights after the map  " << params_[0] << " "
                   << params_[1] << " " << params_[2] << " " << params_[3] << " " << params_[4] << " " << params_[5]
                   << " " << params_[6] << " " << params_[7] << " " << params_[8] << " " << params_[9] << " "
                   << params_[10] << " " << params_[11] << std::endl;

    // we have to transform negative coded in 16 bits into negative coded in 32 bits
    // maybe this should go into the getValue method??

    for (int i = 0; i < 12; ++i) {
      weights_[i] = (params_[i] & 0x8000) ? (int)(params_[i] | 0xffff8000) : (int)(params_[i]);
    }

  } else
    edm::LogWarning("EcalTPG")
        << " EcalEBPhase2TimeReconstructor::setParameters could not find EcalTPGGroupsMap entry for " << raw;
}

Member Data Documentation

ampIn_

uint64_t EcalEBPhase2TimeReconstructor::ampIn_[maxSamplesUsed_]

private

Definition at line 23 of file EcalEBPhase2TimeReconstructor.h.

Referenced by process().

buffer_

int EcalEBPhase2TimeReconstructor::buffer_[maxSamplesUsed_]

private

Definition at line 21 of file EcalEBPhase2TimeReconstructor.h.

Referenced by process(), and setInput().

debug_

bool EcalEBPhase2TimeReconstructor::debug_

private

Definition at line 19 of file EcalEBPhase2TimeReconstructor.h.

Referenced by process(), setInput(), and setParameters().

extraShift_

bool EcalEBPhase2TimeReconstructor::extraShift_[2] = {false, false}

private

Definition at line 25 of file EcalEBPhase2TimeReconstructor.h.

inputsAlreadyIn_

int EcalEBPhase2TimeReconstructor::inputsAlreadyIn_

private

Definition at line 20 of file EcalEBPhase2TimeReconstructor.h.

Referenced by process(), and setInput().

invAmpAr_

uint64_t EcalEBPhase2TimeReconstructor::invAmpAr_[512]

private

Definition at line 61 of file EcalEBPhase2TimeReconstructor.h.

Referenced by process().

maxSamplesUsed_

const int EcalEBPhase2TimeReconstructor::maxSamplesUsed_ = 12

staticprivate

Definition at line 18 of file EcalEBPhase2TimeReconstructor.h.

Referenced by process(), and setInput().

processedOutput_

int EcalEBPhase2TimeReconstructor::processedOutput_

private

Definition at line 28 of file EcalEBPhase2TimeReconstructor.h.

Referenced by process().

shift_

int EcalEBPhase2TimeReconstructor::shift_

private

Definition at line 24 of file EcalEBPhase2TimeReconstructor.h.

Referenced by process().

weights_

int EcalEBPhase2TimeReconstructor::weights_[maxSamplesUsed_]

private

Definition at line 22 of file EcalEBPhase2TimeReconstructor.h.

Referenced by process(), and setParameters().