DTLinearDriftFromDBAlgo Class Reference

#include <DTLinearDriftFromDBAlgo.h>

Inheritance diagram for DTLinearDriftFromDBAlgo:

Public Member Functions
bool	compute (const DTLayer *layer, const DTDigi &digi, LocalPoint &leftPoint, LocalPoint &rightPoint, LocalError &error) const override
bool	compute (const DTLayer *layer, const DTRecHit1D &recHit1D, const float &angle, const GlobalPoint &globPos, DTRecHit1D &newHit1D) const override
bool	compute (const DTLayer *layer, const DTRecHit1D &recHit1D, const float &angle, DTRecHit1D &newHit1D) const override
	DTLinearDriftFromDBAlgo (const edm::ParameterSet &config)
	Constructor. More...
void	setES (const edm::EventSetup &setup) override
	Pass the Event Setup to the algo at each event. More...
	~DTLinearDriftFromDBAlgo () override
	Destructor. More...
Public Member Functions inherited from DTRecHitBaseAlgo
	DTRecHitBaseAlgo (const edm::ParameterSet &config)
	Constructor. More...
virtual edm::OwnVector< DTRecHit1DPair >	reconstruct (const DTLayer *layer, const DTLayerId &layerId, const DTDigiCollection::Range &digiRange)
	Build all hits in the range associated to the layerId, at the 1st step. More...
virtual	~DTRecHitBaseAlgo ()
	Destructor. More...

Private Member Functions
virtual bool	compute (const DTLayer *layer, const DTWireId &wireId, const float digiTime, const GlobalPoint &globPos, DTRecHit1D &newHit1D, int step) const
virtual bool	compute (const DTLayer *layer, const DTWireId &wireId, const float digiTime, const GlobalPoint &globPos, LocalPoint &leftPoint, LocalPoint &rightPoint, LocalError &error, int step) const

Private Attributes
const bool	debug
const bool	doVdriftCorr
const MagneticField *	field
const float	maxTime
const float	minTime
const DTMtime *	mTimeMap
int	nominalB
const bool	stepTwoFromDigi
const DTRecoConditions *	uncertMap
bool	useUncertDB

Additional Inherited Members
Protected Attributes inherited from DTRecHitBaseAlgo
std::unique_ptr< DTTTrigBaseSync >	theSync

Detailed Description

Concrete implementation of DTRecHitBaseAlgo. Compute drift distance using constant drift velocity read from database.

Author

S.Bolognesi - INFN Torino

Definition at line 18 of file DTLinearDriftFromDBAlgo.h.

Constructor & Destructor Documentation

DTLinearDriftFromDBAlgo()

DTLinearDriftFromDBAlgo::DTLinearDriftFromDBAlgo

(

const edm::ParameterSet &

config

)

Constructor.

Definition at line 26 of file DTLinearDriftFromDBAlgo.cc.

    : DTRecHitBaseAlgo(config),
      mTimeMap(nullptr),
      field(nullptr),
      nominalB(-1),
      minTime(config.getParameter<double>("minTime")),
      maxTime(config.getParameter<double>("maxTime")),
      doVdriftCorr(config.getParameter<bool>("doVdriftCorr")),
      // Option to force going back to digi time at Step 2
      stepTwoFromDigi(config.getParameter<bool>("stepTwoFromDigi")),
      useUncertDB(config.getParameter<bool>("useUncertDB")),
      // Set verbose output
      debug(config.getUntrackedParameter<bool>("debug")) {}

~DTLinearDriftFromDBAlgo()

DTLinearDriftFromDBAlgo::~DTLinearDriftFromDBAlgo ( )

override

Destructor.

Definition at line 40 of file DTLinearDriftFromDBAlgo.cc.

{}

Member Function Documentation

compute() [1/5]

bool DTLinearDriftFromDBAlgo::compute ( const DTLayer *  layer, const DTDigi &  digi, LocalPoint &  leftPoint, LocalPoint &  rightPoint, LocalError &  error  ) const

overridevirtual

First step in computation of Left/Right hits from a Digi. The results are the local position (in DTLayer frame) of the Left and Right hit, and the error (which is common). Returns false on failure. The hit is assumed to be at the wire center.

Implements DTRecHitBaseAlgo.

Definition at line 72 of file DTLinearDriftFromDBAlgo.cc.

                                                                                                                      {
  // Get the wireId
  DTLayerId layerId = layer->id();
  const DTWireId wireId(layerId, digi.wire());
 
  // Get Wire position
  if (!layer->specificTopology().isWireValid(digi.wire()))
    return false;
  LocalPoint locWirePos(layer->specificTopology().wirePosition(digi.wire()), 0, 0);
  const GlobalPoint globWirePos = layer->toGlobal(locWirePos);
 
  return compute(layer, wireId, digi.time(), globWirePos, leftPoint, rightPoint, error, 1);
}

References relativeConstraints::error, DTLayer::id(), DTTopology::isWireValid(), DTLayer::specificTopology(), DTDigi::time(), GeomDet::toGlobal(), DTDigi::wire(), and DTTopology::wirePosition().

Referenced by compute().

compute() [2/5]

bool DTLinearDriftFromDBAlgo::compute ( const DTLayer *  layer, const DTRecHit1D &  recHit1D, const float &  angle, const GlobalPoint &  globPos, DTRecHit1D &  newHit1D  ) const

overridevirtual

Third (and final) step in hits position computation. Also the hit position along the wire is available and can be used to correct the drift time for particle TOF and propagation of signal along the wire. NOTE: Only position and error of the new hit are modified

Implements DTRecHitBaseAlgo.

Definition at line 109 of file DTLinearDriftFromDBAlgo.cc.

                                                                  {
  return compute(layer, recHit1D.wireId(), recHit1D.digiTime(), globPos, newHit1D, 3);
}

References compute(), DTRecHit1D::digiTime(), and DTRecHit1D::wireId().

compute() [3/5]

bool DTLinearDriftFromDBAlgo::compute ( const DTLayer *  layer, const DTRecHit1D &  recHit1D, const float &  angle, DTRecHit1D &  newHit1D  ) const

overridevirtual

Second step in hit position computation. It is the same as first step since the angular information is not used NOTE: Only position and error of the new hit are modified

Implements DTRecHitBaseAlgo.

Definition at line 88 of file DTLinearDriftFromDBAlgo.cc.

                                                                  {
  if (!stepTwoFromDigi) {
    newHit1D.setPositionAndError(recHit1D.localPosition(), recHit1D.localPositionError());
    return true;
  }
 
  const DTWireId wireId = recHit1D.wireId();
 
  // Get Wire position
  if (!layer->specificTopology().isWireValid(wireId.wire()))
    return false;
  LocalPoint locWirePos(layer->specificTopology().wirePosition(wireId.wire()), 0, 0);
  const GlobalPoint globWirePos = layer->toGlobal(locWirePos);
 
  return compute(layer, wireId, recHit1D.digiTime(), globWirePos, newHit1D, 2);
}

References compute(), DTRecHit1D::digiTime(), DTTopology::isWireValid(), DTRecHit1D::localPosition(), DTRecHit1D::localPositionError(), DTRecHit1D::setPositionAndError(), DTLayer::specificTopology(), stepTwoFromDigi, GeomDet::toGlobal(), DTWireId::wire(), DTRecHit1D::wireId(), and DTTopology::wirePosition().

compute() [4/5]

bool DTLinearDriftFromDBAlgo::compute ( const DTLayer *  layer, const DTWireId &  wireId, const float  digiTime, const GlobalPoint &  globPos, DTRecHit1D &  newHit1D, int  step  ) const

privatevirtual

Definition at line 200 of file DTLinearDriftFromDBAlgo.cc.

                                                      {
  LocalPoint leftPoint;
  LocalPoint rightPoint;
  LocalError error;
 
  if (compute(layer, wireId, digiTime, globPos, leftPoint, rightPoint, error, step)) {
    // Set the position and the error of the rechit which is being updated
    switch (newHit1D.lrSide()) {
      case DTEnums::Left: {
        // Keep the original y position of newHit1D: for step==3, it's the
        // position along the wire. Needed for rotation alignment
        LocalPoint leftPoint3D(leftPoint.x(), newHit1D.localPosition().y(), leftPoint.z());
        newHit1D.setPositionAndError(leftPoint3D, error);
        break;
      }
 
      case DTEnums::Right: {
        // as above: 3d position
        LocalPoint rightPoint3D(rightPoint.x(), newHit1D.localPosition().y(), rightPoint.z());
        newHit1D.setPositionAndError(rightPoint3D, error);
        break;
      }
 
      default:
        throw cms::Exception("InvalidDTCellSide") << "[DTLinearDriftFromDBAlgo] Compute at Step " << step
                                                  << ", Hit side " << newHit1D.lrSide() << " is invalid!" << endl;
        return false;
    }
 
    return true;
  } else {
    return false;
  }
}

References compute(), relativeConstraints::error, Exception, DTEnums::Left, DTRecHit1D::localPosition(), DTRecHit1D::lrSide(), DTEnums::Right, DTRecHit1D::setPositionAndError(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

compute() [5/5]

bool DTLinearDriftFromDBAlgo::compute ( const DTLayer *  layer, const DTWireId &  wireId, const float  digiTime, const GlobalPoint &  globPos, LocalPoint &  leftPoint, LocalPoint &  rightPoint, LocalError &  error, int  step  ) const

privatevirtual

cf. http://arxiv.org/PS_cache/arxiv/pdf/0911/0911.4895v2.pdf

Definition at line 118 of file DTLinearDriftFromDBAlgo.cc.

                                                      {
  // Subtract the offset to the digi time accordingly to the DTTTrigBaseSync concrete instance
  float driftTime = digiTime - theSync->offset(layer, wireId, globPos);
 
  // check for out-of-time
  if (driftTime < minTime || driftTime > maxTime) {
    if (debug)
      cout << "[DTLinearDriftFromDBAlgo]*** Drift time out of window for in-time hits " << driftTime << endl;
 
    if (step == 1) {  //FIXME: protection against failure at 2nd and 3rd steps, must be checked!!!
      // Hits are interpreted as coming from out-of-time pile-up and recHit
      // is ignored.
      return false;
    }
  }
 
  // Small negative times interpreted as hits close to the wire.
  if (driftTime < 0.)
    driftTime = 0;
 
  // Read the vDrift and reso for this wire
  float vDrift = 0;
  float hitResolution = 0;
  // vdrift is cm/ns , resolution is cm
  mTimeMap->get(wireId.superlayerId(),
                vDrift,
                hitResolution,  // Value from vdrift DB; replaced below if useUncertDB card is set
                DTVelocityUnits::cm_per_ns);
 
  if (useUncertDB) {
    // Read the uncertainty from the DB for the given channel and step
    double args[1] = {double(step - 1)};
    hitResolution = uncertMap->get(wireId, args);
  }
 
  //only in step 3
  if (doVdriftCorr && step == 3 && nominalB != 0) {
    if (abs(wireId.wheel()) == 2 && wireId.station() == 1 && wireId.superLayer() != 2) {
      // Variation of vdrift along Y due to B field,
      // vdrift is lower  a negative Y (lower global |Z|)
      const float k_param = 1.2e-04;
      LocalPoint local_pos = layer->toLocal(globPos);
      vDrift = vDrift * (1. - k_param * local_pos.y());
    }
  }
 
  // Compute the drift distance
  float drift = driftTime * vDrift;
 
  // Get Wire position
  if (!layer->specificTopology().isWireValid(wireId.wire()))
    return false;
  LocalPoint locWirePos(layer->specificTopology().wirePosition(wireId.wire()), 0, 0);
  //Build the two possible points and the error on the position
  leftPoint = LocalPoint(locWirePos.x() - drift, locWirePos.y(), locWirePos.z());
  rightPoint = LocalPoint(locWirePos.x() + drift, locWirePos.y(), locWirePos.z());
  error = LocalError(hitResolution * hitResolution, 0., 0.);
 
  if (debug) {
    cout << "[DTLinearDriftFromDBAlgo] Compute drift distance, for digi at wire: " << wireId << endl
         << "       Step:           " << step << endl
         << "       Digi time:      " << digiTime << endl
         << "       Drift time:     " << driftTime << endl
         << "       Drift distance: " << drift << endl
         << "       Hit Resolution: " << hitResolution << endl
         << "       Left point:     " << leftPoint << endl
         << "       Right point:    " << rightPoint << endl
         << "       Error:          " << error << endl;
  }
 
  return true;
}

References funct::abs(), writedatasetfile::args, DTVelocityUnits::cm_per_ns, gather_cfg::cout, debug, doVdriftCorr, shallow::drift(), relativeConstraints::error, DTRecoConditions::get(), DTMtime::get(), DTLinearDriftAlgo_cfi::hitResolution, DTTopology::isWireValid(), maxTime, mTimeMap, nominalB, DTLayer::specificTopology(), DTChamberId::station(), DTSuperLayerId::superLayer(), DTLayerId::superlayerId(), DTRecHitBaseAlgo::theSync, GeomDet::toLocal(), uncertMap, useUncertDB, DTChamberId::wheel(), DTWireId::wire(), DTTopology::wirePosition(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

setES()

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

overridevirtual

Pass the Event Setup to the algo at each event.

Implements DTRecHitBaseAlgo.

Definition at line 42 of file DTLinearDriftFromDBAlgo.cc.

                                                           {
  if (debug)
    cout << "[DTLinearDriftFromDBAlgo] setES called" << endl;
  theSync->setES(setup);
  // Get the map of ttrig from the Setup
  ESHandle<DTMtime> mTimeHandle;
  setup.get<DTMtimeRcd>().get(mTimeHandle);
  mTimeMap = &*mTimeHandle;
 
  ESHandle<MagneticField> magfield;
  setup.get<IdealMagneticFieldRecord>().get(magfield);
  field = &*magfield;
  nominalB = field->nominalValue();
 
  if (useUncertDB) {
    ESHandle<DTRecoConditions> uncerts;
    setup.get<DTRecoConditionsUncertRcd>().get(uncerts);
    uncertMap = &*uncerts;
    if (uncertMap->version() > 1)
      edm::LogError("NotImplemented") << "DT Uncertainty DB version unsupported: " << uncertMap->version();
  }
 
  if (debug) {
    cout << "[DTLinearDriftFromDBAlgo] meanTimer version: " << mTimeMap->version() << endl;
    if (useUncertDB)
      cout << "                          uncertDB  version: " << uncertMap->version() << endl;
  }
}

References gather_cfg::cout, debug, field, get, volumeBasedMagneticField_160812_cfi::magfield, mTimeMap, nominalB, MagneticField::nominalValue(), singleTopDQM_cfi::setup, DTRecHitBaseAlgo::theSync, uncertMap, useUncertDB, DTRecoConditions::version(), and DTMtime::version().

Member Data Documentation

debug

const bool DTLinearDriftFromDBAlgo::debug

private

Definition at line 107 of file DTLinearDriftFromDBAlgo.h.

Referenced by compute(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), runTauIdMVA.TauIDEmbedder::runTauID(), and setES().

doVdriftCorr

const bool DTLinearDriftFromDBAlgo::doVdriftCorr

private

Definition at line 96 of file DTLinearDriftFromDBAlgo.h.

Referenced by compute().

field

const MagneticField* DTLinearDriftFromDBAlgo::field

private

Definition at line 83 of file DTLinearDriftFromDBAlgo.h.

Referenced by setES().

maxTime

const float DTLinearDriftFromDBAlgo::maxTime

private

Definition at line 93 of file DTLinearDriftFromDBAlgo.h.

Referenced by compute().

minTime

const float DTLinearDriftFromDBAlgo::minTime

private

Definition at line 90 of file DTLinearDriftFromDBAlgo.h.

mTimeMap

const DTMtime* DTLinearDriftFromDBAlgo::mTimeMap

private

Definition at line 80 of file DTLinearDriftFromDBAlgo.h.

Referenced by compute(), and setES().

nominalB

int DTLinearDriftFromDBAlgo::nominalB

private

Definition at line 84 of file DTLinearDriftFromDBAlgo.h.

Referenced by compute(), and setES().

stepTwoFromDigi

const bool DTLinearDriftFromDBAlgo::stepTwoFromDigi

private

Definition at line 100 of file DTLinearDriftFromDBAlgo.h.

Referenced by compute().

uncertMap

const DTRecoConditions* DTLinearDriftFromDBAlgo::uncertMap

private

Definition at line 87 of file DTLinearDriftFromDBAlgo.h.

Referenced by compute(), and setES().

useUncertDB

bool DTLinearDriftFromDBAlgo::useUncertDB

private

Definition at line 104 of file DTLinearDriftFromDBAlgo.h.

Referenced by compute(), and setES().