MTDTimeCalib.cc

Go to the documentation of this file.

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

#include "DataFormats/ForwardDetId/interface/BTLDetId.h"
#include "DataFormats/ForwardDetId/interface/ETLDetId.h"

#include "Geometry/CommonDetUnit/interface/GeomDetType.h"
#include "Geometry/MTDGeometryBuilder/interface/ProxyMTDTopology.h"
#include "Geometry/MTDGeometryBuilder/interface/RectangularMTDTopology.h"
#include "Geometry/MTDNumberingBuilder/interface/GeometricTimingDet.h"

MTDTimeCalib::MTDTimeCalib(edm::ParameterSet const& conf, 
               const MTDGeometry* geom, const MTDTopology* topo):
  geom_(geom),
  topo_(topo),
  btlTimeOffset_( conf.getParameter<double>("BTLTimeOffset") ),
  etlTimeOffset_( conf.getParameter<double>("ETLTimeOffset") ),
  btlLightCollTime_( conf.getParameter<double>("BTLLightCollTime") ),
  btlLightCollSlope_( conf.getParameter<double>("BTLLightCollSlope") )
{
}
 

float MTDTimeCalib::getTimeCalib(const MTDDetId& id) const
{
  if (id.subDetector() != MTDDetId::FastTime)
    {
      throw cms::Exception("MTDTimeCalib") << "MTDDetId: " << std::hex
                              << id.rawId()
                              << " is invalid!" << std::dec
                              << std::endl;
    }

  float time_calib = 0.;

  if ( id.mtdSubDetector() == MTDDetId::BTL )
    {
      time_calib += btlTimeOffset_;
      BTLDetId hitId(id);
      DetId geoId = hitId.geographicalId( (BTLDetId::CrysLayout) topo_->getMTDTopologyMode() ); //for BTL topology gives different layout id
      const MTDGeomDet* thedet = geom_->idToDet(geoId);
      
      if( thedet == nullptr ) {
    throw cms::Exception("MTDTimeCalib") << "GeographicalID: " << std::hex
                          << geoId.rawId()
                        << " (" << id.rawId()<< ") is invalid!" << std::dec
                        << std::endl;
      }
      const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
      const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());    
      
      if ( topo_->getMTDTopologyMode() == (int) BTLDetId::CrysLayout::tile )
    {
      time_calib -= btlLightCollTime_; //simply remove the offset introduced at sim level
    }
      else if ( topo_->getMTDTopologyMode() == (int) BTLDetId::CrysLayout::bar ||
        topo_->getMTDTopologyMode() == (int) BTLDetId::CrysLayout::barphiflat 
        )
    {
      //for bars in phi
      time_calib -= 0.5*topo.pitch().first*btlLightCollSlope_; //time offset for bar time is L/2v 
    }
      else if ( topo_->getMTDTopologyMode() == (int) BTLDetId::CrysLayout::barzflat ) 
    {
      //for bars in z
      time_calib -= 0.5*topo.pitch().second*btlLightCollSlope_; //time offset for bar time is L/2v 
    }
    }
  else if ( id.mtdSubDetector() == MTDDetId::ETL )
    {
      time_calib += etlTimeOffset_;
    }
  else
    {
      throw cms::Exception("MTDTimeCalib") << "MTDDetId: " << std::hex
                              << id.rawId()
                              << " is invalid!" << std::dec
                              << std::endl;
    }

  return time_calib;
}

#include "FWCore/Utilities/interface/typelookup.h"

//--- Now use the Framework macros to set it all up:
TYPELOOKUP_DATA_REG(MTDTimeCalib);