DTTTrigSyncFromDB Class Reference

Inheritance diagram for DTTTrigSyncFromDB:

Public Member Functions
	DTTTrigSyncFromDB (const edm::ParameterSet &config, edm::ConsumesCollector)
	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
edm::ESGetToken< DTT0, DTT0Rcd >	t0Token_
double	theBXspace
int	theTOFCorrType
double	theVPropWire
int	theWirePropCorrType
const DTTtrig *	tTrigMap
const edm::ESGetToken< DTTtrig, DTTtrigRcd >	ttrigToken_
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 57 of file DTTTrigSyncFromDB.cc.

Constructor & Destructor Documentation

DTTTrigSyncFromDB()

DTTTrigSyncFromDB::DTTTrigSyncFromDB	(	const edm::ParameterSet &	config,
		edm::ConsumesCollector	cc
	)

Constructor.

Definition at line 116 of file DTTTrigSyncFromDB.cc.

References gpuPixelDoublets::cc, doT0Correction, and t0Token_.

    : ttrigToken_(cc.esConsumes(edm::ESInputTag("", config.getParameter<string>("tTrigLabel")))),
      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.)) {
  if (doT0Correction) {
    t0Token_ = cc.esConsumes(edm::ESInputTag("", config.getParameter<string>("t0Label")));
  }
}

~DTTTrigSyncFromDB()

DTTTrigSyncFromDB::~DTTTrigSyncFromDB ( )

override

Destructor.

Definition at line 136 of file DTTTrigSyncFromDB.cc.

{}

Member Function Documentation

emulatorOffset()

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 268 of file DTTTrigSyncFromDB.cc.

References gather_cfg::cout, doT0Correction, DTT0::get(), DTTtrig::get(), createfilelist::int, dttriganalyzer_cfi::kFactor, DTTimeUnits::ns, DTLayerId::superlayerId(), FrontierCondition_GT_autoExpress_cfi::t0, theBXspace, dttriganalyzer_cfi::tTrig, 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;
}

offset() [1/2]

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 154 of file DTTTrigSyncFromDB.cc.

References gather_cfg::cout, debug, doTOFCorrection, doWirePropCorrection, Exception, nano_mu_digi_cff::layer, PV3DBase< T, PVType, FrameType >::mag(), mag(), GloballyPositioned< T >::position(), GeomDet::surface(), theTOFCorrType, theVPropWire, theWirePropCorrType, dttriganalyzer_cfi::tTrig, and DTWireId::wire().

                                                        {
  // 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;
}

offset() [2/2]

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 247 of file DTTTrigSyncFromDB.cc.

References gather_cfg::cout, doT0Correction, DTT0::get(), DTTtrig::get(), dttriganalyzer_cfi::kFactor, DTTimeUnits::ns, DTLayerId::superlayerId(), FrontierCondition_GT_autoExpress_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;
}

setES()

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

overridevirtual

Pass the Event Setup to the algo at each event.

Implements DTTTrigBaseSync.

Definition at line 138 of file DTTTrigSyncFromDB.cc.

References gather_cfg::cout, debug, doT0Correction, singleTopDQM_cfi::setup, t0Token_, tTrigMap, ttrigToken_, tZeroMap, DTT0::version(), and DTTtrig::version().

                                                     {
  if (doT0Correction) {
    // Get the map of t0 from pulses from the Setup
    tZeroMap = &setup.getData(t0Token_);
    if (debug) {
      cout << "[DTTTrigSyncFromDB] t0 version: " << tZeroMap->version() << endl;
    }
  }

  // Get the map of ttrig from the Setup
  tTrigMap = &setup.getData(ttrigToken_);
  if (debug) {
    cout << "[DTTTrigSyncFromDB] ttrig version: " << tTrigMap->version() << endl;
  }
}

Member Data Documentation

debug

const bool DTTTrigSyncFromDB::debug

private

Definition at line 98 of file DTTTrigSyncFromDB.cc.

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().

doT0Correction

bool DTTTrigSyncFromDB::doT0Correction

private

Definition at line 102 of file DTTTrigSyncFromDB.cc.

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

doTOFCorrection

bool DTTTrigSyncFromDB::doTOFCorrection

private

Definition at line 104 of file DTTTrigSyncFromDB.cc.

Referenced by offset().

doWirePropCorrection

bool DTTTrigSyncFromDB::doWirePropCorrection

private

Definition at line 107 of file DTTTrigSyncFromDB.cc.

Referenced by offset().

t0Token_

edm::ESGetToken<DTT0, DTT0Rcd> DTTTrigSyncFromDB::t0Token_

private

Definition at line 93 of file DTTTrigSyncFromDB.cc.

Referenced by DTTTrigSyncFromDB(), and setES().

theBXspace

double DTTTrigSyncFromDB::theBXspace

private

Definition at line 110 of file DTTTrigSyncFromDB.cc.

Referenced by emulatorOffset().

theTOFCorrType

int DTTTrigSyncFromDB::theTOFCorrType

private

Definition at line 105 of file DTTTrigSyncFromDB.cc.

Referenced by offset().

theVPropWire

double DTTTrigSyncFromDB::theVPropWire

private

Definition at line 100 of file DTTTrigSyncFromDB.cc.

Referenced by offset().

theWirePropCorrType

int DTTTrigSyncFromDB::theWirePropCorrType

private

Definition at line 108 of file DTTTrigSyncFromDB.cc.

Referenced by offset().

tTrigMap

const DTTtrig* DTTTrigSyncFromDB::tTrigMap

private

Definition at line 96 of file DTTTrigSyncFromDB.cc.

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

ttrigToken_

const edm::ESGetToken<DTTtrig, DTTtrigRcd> DTTTrigSyncFromDB::ttrigToken_

private

Definition at line 94 of file DTTTrigSyncFromDB.cc.

Referenced by setES().

tZeroMap

const DTT0* DTTTrigSyncFromDB::tZeroMap

private

Definition at line 95 of file DTTTrigSyncFromDB.cc.

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