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" );