BTLUncalibRecHitAlgo Class Reference

Inheritance diagram for BTLUncalibRecHitAlgo:

Public Member Functions
	BTLUncalibRecHitAlgo (const edm::ParameterSet &conf, edm::ConsumesCollector &sumes)
	Constructor. More...
void	getEvent (const edm::Event &) final
	get event and eventsetup information More...
void	getEventSetup (const edm::EventSetup &) final
FTLUncalibratedRecHit	makeRecHit (const BTLDataFrame &dataFrame) const final
	make the rec hit More...
	~BTLUncalibRecHitAlgo () override
	Destructor. More...
Public Member Functions inherited from MTDUncalibratedRecHitAlgoBase< DataFrame >
virtual FTLUncalibratedRecHit	makeRecHit (const DataFrame &dataFrame) const =0
	make the rec hit More...
	MTDUncalibratedRecHitAlgoBase (const edm::ParameterSet &conf, edm::ConsumesCollector &sumes)
	Constructor. More...
const std::string &	name () const
virtual	~MTDUncalibratedRecHitAlgoBase ()
	Destructor. More...

Private Attributes
const double	adcLSB_
const uint32_t	adcNBits_
const double	adcSaturation_
const double	c_LYSO_
const double	timeCorr_p0_
const double	timeCorr_p1_
const double	timeCorr_p2_
const reco::FormulaEvaluator	timeError_
const double	toaLSBToNS_

Detailed Description

Definition at line 8 of file BTLUncalibRecHitAlgo.cc.

Constructor & Destructor Documentation

BTLUncalibRecHitAlgo()

BTLUncalibRecHitAlgo::BTLUncalibRecHitAlgo ( const edm::ParameterSet &  conf, edm::ConsumesCollector &  sumes  )

inline

Constructor.

Definition at line 11 of file BTLUncalibRecHitAlgo.cc.

      : MTDUncalibratedRecHitAlgoBase<BTLDataFrame>(conf, sumes),
        adcNBits_(conf.getParameter<uint32_t>("adcNbits")),
        adcSaturation_(conf.getParameter<double>("adcSaturation")),
        adcLSB_(adcSaturation_ / (1 << adcNBits_)),
        toaLSBToNS_(conf.getParameter<double>("toaLSB_ns")),
        timeError_(conf.getParameter<std::string>("timeResolutionInNs")),
        timeCorr_p0_(conf.getParameter<double>("timeCorr_p0")),
        timeCorr_p1_(conf.getParameter<double>("timeCorr_p1")),
        timeCorr_p2_(conf.getParameter<double>("timeCorr_p2")),
        c_LYSO_(conf.getParameter<double>("c_LYSO")) {}

~BTLUncalibRecHitAlgo()

BTLUncalibRecHitAlgo::~BTLUncalibRecHitAlgo ( )

inlineoverride

Destructor.

Definition at line 24 of file BTLUncalibRecHitAlgo.cc.

{}

Member Function Documentation

getEvent()

void BTLUncalibRecHitAlgo::getEvent ( const edm::Event &  )

inlinefinalvirtual

get event and eventsetup information

Implements MTDUncalibratedRecHitAlgoBase< DataFrame >.

Definition at line 27 of file BTLUncalibRecHitAlgo.cc.

{}

getEventSetup()

void BTLUncalibRecHitAlgo::getEventSetup ( const edm::EventSetup &  )

inlinefinalvirtual

Implements MTDUncalibratedRecHitAlgoBase< DataFrame >.

Definition at line 28 of file BTLUncalibRecHitAlgo.cc.

{}

makeRecHit()

FTLUncalibratedRecHit BTLUncalibRecHitAlgo::makeRecHit ( const BTLDataFrame &  dataFrame ) const

final

make the rec hit

Definition at line 45 of file BTLUncalibRecHitAlgo.cc.

References adcLSB_, CustomPhysics_cfi::amplitude, c_LYSO_, FTLDataFrameT< D, S, DECODE >::column(), reco::FormulaEvaluator::evaluate(), RemoveAddSevLevel::flag, nano_mu_digi_cff::float, FTLDataFrameT< D, S, DECODE >::id(), LogDebug, nHits, position, BTLRecHitsErrorEstimatorIM::positionError(), conifer::pow(), FTLDataFrameT< D, S, DECODE >::row(), FTLDataFrameT< D, S, DECODE >::sample(), hcalRecHitTable_cff::time, timeCorr_p0_, timeCorr_p1_, timeCorr_p2_, timeError_, and toaLSBToNS_.

                                                                                          {
  // The reconstructed amplitudes and times are saved in a std::pair
  //    BTL tile geometry (1 SiPM): only the first value of the amplitude
  //                                and time pairs is used.
  //    BTL bar geometry (2 SiPMs): both values of the amplitude and
  //                                time pairs are filled.

  std::pair<float, float> amplitude(0., 0.);
  std::pair<float, float> time(0., 0.);

  unsigned char flag = 0;

  const auto& sampleRight = dataFrame.sample(0);
  const auto& sampleLeft = dataFrame.sample(1);

  double nHits = 0.;

  LogDebug("BTLUncalibRecHit") << "Original input time t1,t2 " << float(sampleRight.toa()) * toaLSBToNS_ << ", "
                               << float(sampleLeft.toa()) * toaLSBToNS_ << std::endl;

  if (sampleRight.data() > 0) {
    amplitude.first = float(sampleRight.data()) * adcLSB_;
    time.first = float(sampleRight.toa()) * toaLSBToNS_;

    nHits += 1.;

    // Correct the time of the left SiPM for the time-walk
    time.first -= timeCorr_p0_ * pow(amplitude.first, timeCorr_p1_) + timeCorr_p2_;
    flag |= 0x1;
  }

  // --- If available, reconstruct the amplitude and time of the second SiPM
  if (sampleLeft.data() > 0) {
    amplitude.second = float(sampleLeft.data()) * adcLSB_;
    time.second = float(sampleLeft.toa()) * toaLSBToNS_;

    nHits += 1.;

    // Correct the time of the right SiPM for the time-walk
    time.second -= timeCorr_p0_ * pow(amplitude.second, timeCorr_p1_) + timeCorr_p2_;
    flag |= (0x1 << 1);
  }

  // --- Calculate the error on the hit time using the provided parameterization

  const std::array<double, 1> amplitudeV = {{(amplitude.first + amplitude.second) / nHits}};
  const std::array<double, 1> emptyV = {{0.}};

  double timeError = (nHits > 0. ? timeError_.evaluate(amplitudeV, emptyV) : -1.);

  // Calculate the position
  // Distance from center of bar to hit

  float position = 0.5f * (c_LYSO_ * (time.second - time.first));
  float positionError = BTLRecHitsErrorEstimatorIM::positionError();

  LogDebug("BTLUncalibRecHit") << "DetId: " << dataFrame.id().rawId() << " x position = " << position << " +/- "
                               << positionError;
  LogDebug("BTLUncalibRecHit") << "ADC+: set the charge to: (" << amplitude.first << ", " << amplitude.second << ")  ("
                               << sampleRight.data() << ", " << sampleLeft.data() << ") " << adcLSB_ << ' '
                               << std::endl;
  LogDebug("BTLUncalibRecHit") << "TDC+: set the time to: (" << time.first << ", " << time.second << ")  ("
                               << sampleRight.toa() << ", " << sampleLeft.toa() << ") " << toaLSBToNS_ << ' '
                               << std::endl;

  return FTLUncalibratedRecHit(
      dataFrame.id(), dataFrame.row(), dataFrame.column(), amplitude, time, timeError, position, positionError, flag);
}

Member Data Documentation

adcLSB_

const double BTLUncalibRecHitAlgo::adcLSB_

private

Definition at line 36 of file BTLUncalibRecHitAlgo.cc.

Referenced by makeRecHit().

adcNBits_

const uint32_t BTLUncalibRecHitAlgo::adcNBits_

private

Definition at line 34 of file BTLUncalibRecHitAlgo.cc.

adcSaturation_

const double BTLUncalibRecHitAlgo::adcSaturation_

private

Definition at line 35 of file BTLUncalibRecHitAlgo.cc.

c_LYSO_

const double BTLUncalibRecHitAlgo::c_LYSO_

private

Definition at line 42 of file BTLUncalibRecHitAlgo.cc.

Referenced by makeRecHit().

timeCorr_p0_

const double BTLUncalibRecHitAlgo::timeCorr_p0_

private

Definition at line 39 of file BTLUncalibRecHitAlgo.cc.

Referenced by makeRecHit().

timeCorr_p1_

const double BTLUncalibRecHitAlgo::timeCorr_p1_

private

Definition at line 40 of file BTLUncalibRecHitAlgo.cc.

Referenced by makeRecHit().

timeCorr_p2_

const double BTLUncalibRecHitAlgo::timeCorr_p2_

private

Definition at line 41 of file BTLUncalibRecHitAlgo.cc.

Referenced by makeRecHit().

timeError_

const reco::FormulaEvaluator BTLUncalibRecHitAlgo::timeError_

private

Definition at line 38 of file BTLUncalibRecHitAlgo.cc.

Referenced by makeRecHit().

toaLSBToNS_

const double BTLUncalibRecHitAlgo::toaLSBToNS_

private

Definition at line 37 of file BTLUncalibRecHitAlgo.cc.

Referenced by makeRecHit().