#include <DTTTrigSyncFromDB.h>

Inheritance diagram for DTTTrigSyncFromDB:

Public Member Functions
	DTTTrigSyncFromDB (const edm::ParameterSet &config)
	Constructor. More...

double	emulatorOffset (const DTWireId &wireId, double &tTrig, double &t0cell) const override

double	offset (const DTLayer *layer, const DTWireId &wireId, const GlobalPoint &globPos, double &tTrig, double &wirePropCorr, double &tofCorr) const override

double	offset (const DTWireId &wireId) const override

void	setES (const edm::EventSetup &setup) override
	Pass the Event Setup to the algo at each event. More...

	~DTTTrigSyncFromDB () override
	Destructor. More...

Public Member Functions inherited from DTTTrigBaseSync
	DTTTrigBaseSync ()
	Constructor. More...

virtual double	emulatorOffset (const DTWireId &wireId) const

double	offset (const DTLayer *layer, const DTWireId &wireId, const GlobalPoint &globalPos) const

virtual	~DTTTrigBaseSync ()
	Destructor. More...

Private Attributes
const bool	debug

bool	doT0Correction

bool	doTOFCorrection

bool	doWirePropCorrection

double	theBXspace

int	theTOFCorrType

std::string	thetTrigLabel

double	theVPropWire

int	theWirePropCorrType

const DTTtrig *	tTrigMap

const DTT0 *	tZeroMap

Detailed Description

Concrete implementation of a DTTTrigBaseSync. This class define the offset for RecHit building of data and simulation. The offset is computes as:
offset = t0 + tTrig + wirePropCorr - tofCorr
where:

t0 from test pulses (taken from DB, it is assumed to be in ns; can be switched off)
ttrig from the fit of time boxrising edge (taken from DB, it is assumed to be in ns) (At the moment a single value is read for ttrig offset but this may change in the future)
signal propagation along the wire (can be switched off): it is assumed the ttrig accounts on average for correction from the center of the wire to the frontend. Here we just have to correct for the distance of the hit from the wire center.
TOF correction (can be switched off for cosmics): the ttrig already accounts for average TOF correction, depending on the granularity used for the ttrig computation we just have to correct for the TOF from the center of the chamber, SL, layer or wire to the hit position. NOTE: particles are assumed as coming from the IP.

The emulatorOffset is computed as:
offset = int(ttrig/BXspace)*BXspace + t0
where:

t0 from test pulses (taken from DB, it is assumed to be in ns; can be switched off)
ttrig from the fit of time box rising edge (taken from DB, it is assumed to be in ns)
BXspace BX spacing (in ns). Can be configured.

NOTE: this should approximate what is seen online by the BTI

Author: G. Cerminara - INFN Torino

Definition at line 52 of file DTTTrigSyncFromDB.h.

Constructor & Destructor Documentation

DTTTrigSyncFromDB::DTTTrigSyncFromDB ( const edm::ParameterSet & config )

Constructor.

Definition at line 25 of file DTTTrigSyncFromDB.cc.

     : debug(config.getUntrackedParameter<bool>("debug")),
       // The velocity of signal propagation along the wire (cm/ns)
       theVPropWire(config.getParameter<double>("vPropWire")),
       // Switch on/off the T0 correction from pulses
       doT0Correction(config.getParameter<bool>("doT0Correction")),
       // Switch on/off the TOF correction for particles from IP
       doTOFCorrection(config.getParameter<bool>("doTOFCorrection")),
       theTOFCorrType(config.getParameter<int>("tofCorrType")),
       // Switch on/off the correction for the signal propagation along the wire
       doWirePropCorrection(config.getParameter<bool>("doWirePropCorrection")),
       theWirePropCorrType(config.getParameter<int>("wirePropCorrType")),
       // spacing of BX in ns
       theBXspace(config.getUntrackedParameter<double>("bxSpace", 25.)),
       thetTrigLabel(config.getParameter<string>("tTrigLabel")) {}

DTTTrigSyncFromDB::~DTTTrigSyncFromDB ( )

override

Destructor.

Definition at line 41 of file DTTTrigSyncFromDB.cc.

41 {}

Member Function Documentation

double DTTTrigSyncFromDB::emulatorOffset	(	const DTWireId &	wireId,
		double &	tTrig,
		double &	t0cell
	)		const

overridevirtual

Time (ns) to be subtracted to the digi time for emulation purposes It does not take into account TOF and signal propagation along the wire It also returns the different contributions separately:

tTrig is the offset (t_trig)
t0cell is the t0 from pulses

Implements DTTTrigBaseSync.

Definition at line 177 of file DTTTrigSyncFromDB.cc.

References gather_cfg::cout, doT0Correction, DTT0::get(), DTTtrig::get(), createfilelist::int, DTAnalyzerDetailed_cfi::kFactor, DTTimeUnits::ns, DTLayerId::superlayerId(), cscNeutronWriter_cfi::t0, theBXspace, tTrigMap, and tZeroMap.

                                                                                                     {
   float t0 = 0;
   float t0rms = 0;
   if (doT0Correction) {
     // Read the t0 from pulses for this wire (ns)
     tZeroMap->get(wireId, t0, t0rms, DTTimeUnits::ns);
   }
 
   // Read the ttrig for this wire
   float ttrigMean = 0;
   float ttrigSigma = 0;
   float kFactor = 0;
   // FIXME: should check the return value of the DTTtrigRcd::get(..) method
   if (tTrigMap->get(wireId.superlayerId(), ttrigMean, ttrigSigma, kFactor, DTTimeUnits::ns) != 0) {
     cout << "[DTTTrigSyncFromDB]*Error: ttrig not found for SL: " << wireId.superlayerId() << endl;
     //     FIXME: LogError.....
   }
 
   tTrig = ttrigMean + kFactor * ttrigSigma;
   t0cell = t0;
 
   return int(tTrig / theBXspace) * theBXspace + t0cell;
 }

double DTTTrigSyncFromDB::offset	(	const DTLayer *	layer,
		const DTWireId &	wireId,
		const GlobalPoint &	globPos,
		double &	tTrig,
		double &	wirePropCorr,
		double &	tofCorr
	)		const

overridevirtual

Time (ns) to be subtracted to the digi time, Parameters are the layer and the wireId to which the digi is referred and the estimation of the 3D hit position (globPos)

Implements DTTTrigBaseSync.

Definition at line 63 of file DTTTrigSyncFromDB.cc.

References DTTopology::cellLenght(), gather_cfg::cout, debug, doTOFCorrection, doWirePropCorrection, Exception, PV3DBase< T, PVType, FrameType >::mag(), mag(), GloballyPositioned< T >::position(), DTLayer::specificTopology(), DTLayer::superLayer(), GeomDet::surface(), theTOFCorrType, theVPropWire, theWirePropCorrType, GeomDet::toGlobal(), GeomDet::toLocal(), DTWireId::wire(), DTTopology::wirePosition(), and PV3DBase< T, PVType, FrameType >::y().

                                                         {
   // Correction for the float to int conversion while writeing the ttrig in ns into an int variable
   // (half a bin on average)
   // FIXME: this should disappear as soon as the ttrig object will become a float
   //   static const float f2i_convCorr = (25./64.); // ns //FIXME: check how the conversion is performed
 
   tTrig = offset(wireId);
 
   // Compute the time spent in signal propagation along wire.
   // NOTE: the FE is always at y>0
   wirePropCorr = 0;
   if (doWirePropCorrection) {
     switch (theWirePropCorrType) {
         // The ttrig computed from the timebox accounts on average for the signal propagation time
         // from the center of the wire to the frontend. Here we just have to correct for
         // the distance of the hit from the wire center.
       case 0: {
         float wireCoord = layer->toLocal(globPos).y();
         wirePropCorr = -wireCoord / theVPropWire;
         break;
         // FIXME: What if hits used for the time box are not distributed uniformly along the wire?
       }
       //On simulated data you need to subtract the total propagation time
       case 1: {
         float halfL = layer->specificTopology().cellLenght() / 2;
         float wireCoord = layer->toLocal(globPos).y();
         float propgL = halfL - wireCoord;
         wirePropCorr = propgL / theVPropWire;
         break;
       }
       default: {
         throw cms::Exception("[DTTTrigSyncFromDB]")
             << " Invalid parameter: wirePropCorrType = " << theWirePropCorrType << std::endl;
         break;
       }
     }
   }
 
   // Compute TOF correction:
   tofCorr = 0.;
   // TOF Correction can be switched off with appropriate parameter
   if (doTOFCorrection) {
     float flightToHit = globPos.mag();
     static const float cSpeed = 29.9792458;  // cm/ns
     switch (theTOFCorrType) {
       case 0: {
         // The ttrig computed from the real data accounts on average for the TOF correction
         // Depending on the granularity used for the ttrig computation we just have to correct for the
         // TOF from the center of the chamber, SL, layer or wire to the hit position.
         // At the moment only SL granularity is considered
         // Correction for TOF from the center of the SL to hit position
         const DTSuperLayer* sl = layer->superLayer();
         double flightToSL = sl->surface().position().mag();
         tofCorr = (flightToSL - flightToHit) / cSpeed;
         break;
       }
       case 1: {
         // On simulated data you need to consider only the TOF from 3D center of the wire to hit position
         // (because the TOF from the IP to the wire has been already subtracted in the digitization:
         // SimMuon/DTDigitizer/DTDigiSyncTOFCorr.cc corrType=2)
         float flightToWire =
             layer->toGlobal(LocalPoint(layer->specificTopology().wirePosition(wireId.wire()), 0., 0.)).mag();
         tofCorr = (flightToWire - flightToHit) / cSpeed;
         break;
       }
       default: {
         throw cms::Exception("[DTTTrigSyncFromDB]")
             << " Invalid parameter: tofCorrType = " << theTOFCorrType << std::endl;
         break;
       }
     }
   }
 
   if (debug) {
     cout << "[DTTTrigSyncFromDB] Channel: " << wireId << endl
          << "      Offset (ns): " << tTrig + wirePropCorr - tofCorr << endl
          << "      various contributions are: " << endl
          << "      tTrig + t0 (ns):   " << tTrig
          << endl
          //<< "      tZero (ns):   " << t0 << endl
          << "      Propagation along wire delay (ns): " << wirePropCorr << endl
          << "      TOF correction (ns): " << tofCorr << endl
          << endl;
   }
   //The global offset is the sum of various contributions
   return tTrig + wirePropCorr - tofCorr;
 }

double DTTTrigSyncFromDB::offset ( const DTWireId & wireId ) const

overridevirtual

Time (ns) to be subtracted to the digi time. It does not take into account TOF and signal propagation along the wire

Implements DTTTrigBaseSync.

Definition at line 156 of file DTTTrigSyncFromDB.cc.

References gather_cfg::cout, doT0Correction, DTT0::get(), DTTtrig::get(), DTAnalyzerDetailed_cfi::kFactor, DTTimeUnits::ns, DTLayerId::superlayerId(), cscNeutronWriter_cfi::t0, tTrigMap, and tZeroMap.

                                                              {
   float t0 = 0;
   float t0rms = 0;
   if (doT0Correction) {
     // Read the t0 from pulses for this wire (ns)
     tZeroMap->get(wireId, t0, t0rms, DTTimeUnits::ns);
   }
 
   // Read the ttrig for this wire
   float ttrigMean = 0;
   float ttrigSigma = 0;
   float kFactor = 0;
   // FIXME: should check the return value of the DTTtrigRcd::get(..) method
   if (tTrigMap->get(wireId.superlayerId(), ttrigMean, ttrigSigma, kFactor, DTTimeUnits::ns) != 0) {
     cout << "[DTTTrigSyncFromDB]*Error: ttrig not found for SL: " << wireId.superlayerId() << endl;
     //     FIXME: LogError.....
   }
 
   return t0 + ttrigMean + kFactor * ttrigSigma;
 }

void DTTTrigSyncFromDB::setES ( const edm::EventSetup & setup )

overridevirtual

Pass the Event Setup to the algo at each event.

Implements DTTTrigBaseSync.

Definition at line 43 of file DTTTrigSyncFromDB.cc.

References gather_cfg::cout, debug, doT0Correction, edm::EventSetup::get(), thetTrigLabel, tTrigMap, tZeroMap, DTT0::version(), and DTTtrig::version().

                                                      {
   if (doT0Correction) {
     // Get the map of t0 from pulses from the Setup
     ESHandle<DTT0> t0Handle;
     setup.get<DTT0Rcd>().get(t0Handle);
     tZeroMap = &*t0Handle;
     if (debug) {
       cout << "[DTTTrigSyncFromDB] t0 version: " << tZeroMap->version() << endl;
     }
   }
 
   // Get the map of ttrig from the Setup
   ESHandle<DTTtrig> ttrigHandle;
   setup.get<DTTtrigRcd>().get(thetTrigLabel, ttrigHandle);
   tTrigMap = &*ttrigHandle;
   if (debug) {
     cout << "[DTTTrigSyncFromDB] ttrig version: " << tTrigMap->version() << endl;
   }
 }

Member Data Documentation

const bool DTTTrigSyncFromDB::debug

private

Definition at line 91 of file DTTTrigSyncFromDB.h.

Referenced by util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), util.rrapi.RRApi::get(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), offset(), runTauIdMVA.TauIDEmbedder::runTauID(), setES(), and pkg.AbstractPkg::write().

bool DTTTrigSyncFromDB::doT0Correction

private

Definition at line 95 of file DTTTrigSyncFromDB.h.

Referenced by emulatorOffset(), offset(), and setES().

bool DTTTrigSyncFromDB::doTOFCorrection

private

Definition at line 97 of file DTTTrigSyncFromDB.h.

Referenced by offset().

bool DTTTrigSyncFromDB::doWirePropCorrection

private

Definition at line 100 of file DTTTrigSyncFromDB.h.

Referenced by offset().

double DTTTrigSyncFromDB::theBXspace

private

Definition at line 103 of file DTTTrigSyncFromDB.h.

Referenced by emulatorOffset().

int DTTTrigSyncFromDB::theTOFCorrType

private

Definition at line 98 of file DTTTrigSyncFromDB.h.

Referenced by offset().

std::string DTTTrigSyncFromDB::thetTrigLabel

private

Definition at line 105 of file DTTTrigSyncFromDB.h.

Referenced by setES().

double DTTTrigSyncFromDB::theVPropWire

private

Definition at line 93 of file DTTTrigSyncFromDB.h.

Referenced by offset().

int DTTTrigSyncFromDB::theWirePropCorrType

private

Definition at line 101 of file DTTTrigSyncFromDB.h.

Referenced by offset().

const DTTtrig* DTTTrigSyncFromDB::tTrigMap

private

Definition at line 89 of file DTTTrigSyncFromDB.h.

Referenced by emulatorOffset(), offset(), and setES().

const DTT0* DTTTrigSyncFromDB::tZeroMap

private

Definition at line 88 of file DTTTrigSyncFromDB.h.

Referenced by emulatorOffset(), offset(), and setES().

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation