DTNoDriftAlgo Class Reference

#include <DTNoDriftAlgo.h>

Inheritance diagram for DTNoDriftAlgo:

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
	DTNoDriftAlgo (const edm::ParameterSet &config)
	Constructor. More...
edm::OwnVector< DTRecHit1DPair >	reconstruct (const DTLayer *layer, const DTLayerId &layerId, const DTDigiCollection::Range &digiRange) override
void	setES (const edm::EventSetup &setup) override
	Pass the Event Setup to the algo at each event. More...
	~DTNoDriftAlgo () override
	Destructor. More...
Public Member Functions inherited from DTRecHitBaseAlgo
	DTRecHitBaseAlgo (const edm::ParameterSet &config)
	Constructor. 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 float	fixedDrift
const float	hitResolution
const float	maxTime
const float	minTime

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

Detailed Description

Concrete implementation of DTRecHitBaseAlgo. Create pair of RecHits at fixed distance from the wire.

Author

Martijn Mulders - CERN (marti.nosp@m.jn.m.nosp@m.ulder.nosp@m.s@ce.nosp@m.rn.ch) based on DTLinearDriftAlgo

Definition at line 15 of file DTNoDriftAlgo.h.

Constructor & Destructor Documentation

DTNoDriftAlgo()

DTNoDriftAlgo::DTNoDriftAlgo

(

const edm::ParameterSet &

config

)

Constructor.

Definition at line 18 of file DTNoDriftAlgo.cc.

    : DTRecHitBaseAlgo(config),
      fixedDrift(config.getParameter<double>("fixedDrift")),
      hitResolution(config.getParameter<double>("hitResolution")),  // Set to size of (half)cell
      minTime(config.getParameter<double>("minTime")),
      maxTime(config.getParameter<double>("maxTime")),
      debug(config.getUntrackedParameter<bool>("debug"))  // Set verbose output
{}

~DTNoDriftAlgo()

DTNoDriftAlgo::~DTNoDriftAlgo ( )

override

Destructor.

Definition at line 27 of file DTNoDriftAlgo.cc.

{}

Member Function Documentation

compute() [1/5]

bool DTNoDriftAlgo::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 78 of file DTNoDriftAlgo.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(), and reconstruct().

compute() [2/5]

bool DTNoDriftAlgo::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 103 of file DTNoDriftAlgo.cc.

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

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

compute() [3/5]

bool DTNoDriftAlgo::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 94 of file DTNoDriftAlgo.cc.

                                                        {
  newHit1D.setPositionAndError(recHit1D.localPosition(), recHit1D.localPositionError());
  return true;
}

References DTRecHit1D::localPosition(), DTRecHit1D::localPositionError(), and DTRecHit1D::setPositionAndError().

compute() [4/5]

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

privatevirtual

Definition at line 166 of file DTNoDriftAlgo.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:
        newHit1D.setPositionAndError(leftPoint, error);
        break;
 
      case DTEnums::Right:
        newHit1D.setPositionAndError(rightPoint, error);
        break;
 
      default:
        throw cms::Exception("InvalidDTCellSide") << "[DTNoDriftAlgo] 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::lrSide(), DTEnums::Right, and DTRecHit1D::setPositionAndError().

compute() [5/5]

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

privatevirtual

Definition at line 112 of file DTNoDriftAlgo.cc.

                                            {
  //}
 
  // Small negative times interpreted as hits close to the wire.
  //if (driftTime<0.) driftTime=0;
 
  // check for out-of-time
  if (digiTime < minTime || digiTime > maxTime) {
    if (debug)
      cout << "[DTNoDriftAlgo]*** Drift time out of window for in-time hits " << digiTime << 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;
    }
  }
 
  // Compute the drift distance
  float drift = fixedDrift;
 
  // 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 << "[DTNoDriftAlgo] Compute drift distance, for digi at wire: " << wireId << endl
         << "       Step:           " << step << endl
         << "       Digi time:      " << digiTime
         << endl
         //  << "       Drift time:     " << driftTime << endl
         << "       Fixed Drift distance: " << drift << endl
         << "       Hit Resolution: " << hitResolution << endl
         << "       Left point:     " << leftPoint << endl
         << "       Right point:    " << rightPoint << endl
         << "       Error:          " << error << endl;
  }
 
  return true;
}

References gather_cfg::cout, debug, shallow::drift(), relativeConstraints::error, fixedDrift, hitResolution, DTTopology::isWireValid(), maxTime, DTLayer::specificTopology(), DTWireId::wire(), DTTopology::wirePosition(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

reconstruct()

OwnVector< DTRecHit1DPair > DTNoDriftAlgo::reconstruct ( const DTLayer *  layer, const DTLayerId &  layerId, const DTDigiCollection::Range &  digiRange  )

overridevirtual

MM: Override virtual function from DTRecHitBaseAlgo--> for the NoDrift algorithm only a maximum of one hit per wire is allowed! Build all hits in the range associated to the layerId, at the 1st step.

Reimplemented from DTRecHitBaseAlgo.

Definition at line 34 of file DTNoDriftAlgo.cc.

                                                                                             {
  OwnVector<DTRecHit1DPair> result;
 
  // Loop over all digis in the given range
  for (DTDigiCollection::const_iterator digi = digiRange.first; digi != digiRange.second; digi++) {
    // Get the wireId
    DTWireId wireId(layerId, (*digi).wire());
 
    bool isDouble = false;
    for (OwnVector<DTRecHit1DPair>::const_iterator doubleWireCheck = result.begin(); doubleWireCheck != result.end();
         doubleWireCheck++) {
      if (wireId == (*doubleWireCheck).wireId()) {
        isDouble = true;
        //  std::cout << " Reject this hit with time " << (*digi).time() << std::endl;
        break;
      }
    }
 
    if (isDouble)
      continue;
 
    LocalError tmpErr;
    LocalPoint lpoint, rpoint;
    // Call the compute method
    bool OK = compute(layer, *digi, lpoint, rpoint, tmpErr);
 
    if (!OK)
      continue;
 
    // Build a new pair of 1D rechit
    DTRecHit1DPair* recHitPair = new DTRecHit1DPair(wireId, *digi);
 
    // Set the position and the error of the 1D rechits
    recHitPair->setPositionAndError(DTEnums::Left, lpoint, tmpErr);
    recHitPair->setPositionAndError(DTEnums::Right, rpoint, tmpErr);
 
    result.push_back(recHitPair);
  }
  return result;
}

References compute(), DTEnums::Left, mps_fire::result, DTEnums::Right, and DTRecHit1DPair::setPositionAndError().

setES()

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

overridevirtual

Pass the Event Setup to the algo at each event.

Implements DTRecHitBaseAlgo.

Definition at line 29 of file DTNoDriftAlgo.cc.

                                                 {
  //  theSync->setES(setup);
}

Member Data Documentation

debug

const bool DTNoDriftAlgo::debug

private

Definition at line 96 of file DTNoDriftAlgo.h.

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

fixedDrift

const float DTNoDriftAlgo::fixedDrift

private

Definition at line 84 of file DTNoDriftAlgo.h.

Referenced by compute().

hitResolution

const float DTNoDriftAlgo::hitResolution

private

Definition at line 87 of file DTNoDriftAlgo.h.

Referenced by compute().

maxTime

const float DTNoDriftAlgo::maxTime

private

Definition at line 93 of file DTNoDriftAlgo.h.

Referenced by compute().

minTime

const float DTNoDriftAlgo::minTime

private

Definition at line 90 of file DTNoDriftAlgo.h.