MTDRecHitAlgo.cc

Go to the documentation of this file.

#include "RecoLocalFastTime/FTLCommonAlgos/interface/MTDRecHitAlgoBase.h"

#include "RecoLocalFastTime/Records/interface/MTDTimeCalibRecord.h"
#include "RecoLocalFastTime/FTLCommonAlgos/interface/MTDTimeCalib.h"

class MTDRecHitAlgo : public MTDRecHitAlgoBase {
public:
  MTDRecHitAlgo(const edm::ParameterSet& conf, edm::ConsumesCollector& sumes)
      : MTDRecHitAlgoBase(conf, sumes),
        thresholdToKeep_(conf.getParameter<double>("thresholdToKeep")),
        calibration_(conf.getParameter<double>("calibrationConstant")) {}

  ~MTDRecHitAlgo() override {}

  void getEvent(const edm::Event&) final {}
  void getEventSetup(const edm::EventSetup&) final;

  FTLRecHit makeRecHit(const FTLUncalibratedRecHit& uRecHit, uint32_t& flags) const final;

private:
  double thresholdToKeep_, calibration_;
  const MTDTimeCalib* time_calib_;
};

void MTDRecHitAlgo::getEventSetup(const edm::EventSetup& es) {
  edm::ESHandle<MTDTimeCalib> pTC;
  es.get<MTDTimeCalibRecord>().get("MTDTimeCalib", pTC);
  time_calib_ = pTC.product();
}

FTLRecHit MTDRecHitAlgo::makeRecHit(const FTLUncalibratedRecHit& uRecHit, uint32_t& flags) const {
  unsigned char flagsWord = uRecHit.flags();
  float timeError = uRecHit.timeError();

  float energy = 0.;
  float time = 0.;

  switch (flagsWord) {
    // BTL bar geometry with only the right SiPM information available
    case 0x2: {
      energy = uRecHit.amplitude().second;
      time = uRecHit.time().second;

      break;
    }
    // BTL bar geometry with left and right SiPMs information available
    case 0x3: {
      energy = 0.5 * (uRecHit.amplitude().first + uRecHit.amplitude().second);
      time = 0.5 * (uRecHit.time().first + uRecHit.time().second);

      break;
    }
    // ETL, BTL tile geometry, BTL bar geometry with only the left SiPM information available
    default: {
      energy = uRecHit.amplitude().first;
      time = uRecHit.time().first;

      break;
    }
  }

  // --- Energy calibration: for the time being this is just a conversion pC --> MeV
  energy *= calibration_;

  // --- Time calibration: for the time being just removes a time offset in BTL
  time += time_calib_->getTimeCalib(uRecHit.id());

  FTLRecHit rh(uRecHit.id(), uRecHit.row(), uRecHit.column(), energy, time, timeError);

  // Now fill flags
  // all rechits from the digitizer are "good" at present
  if (energy > thresholdToKeep_) {
    flags = FTLRecHit::kGood;
    rh.setFlag(flags);
  } else {
    flags = FTLRecHit::kKilled;
    rh.setFlag(flags);
  }

  return rh;
}

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_EDM_PLUGIN(MTDRecHitAlgoFactory, MTDRecHitAlgo, "MTDRecHitAlgo");