Translator mirroring the realistic Digitization. More...

#include <RealisticTDCCountTranslator.h>

Inheritance diagram for RealisticTDCCountTranslator:

Public Member Functions
	RealisticTDCCountTranslator (bool useInWirePropagationDelay=false)
	Constructor, with the additional information, if propagation in the wire shall be considered.

	~RealisticTDCCountTranslator ()
	Destructor.

double	getDriftLength (unsigned short tdcCount, const WireID &wireID=WireID(), double timeOfFlightEstimator=0, bool leftRight=false, double z=0, double alpha=0, double theta=static_cast< double >(TMath::Pi()/2.), unsigned short adcCount=0) override
	Get Drift length.

double	getDriftTime (unsigned short tdcCount, const WireID &wireID, double timeOfFlightEstimator, double z, unsigned short adcCount) override
	Get Drift time.

double	getDriftLengthResolution (double driftLength, const WireID &wireID=WireID(), bool leftRight=false, double z=0, double alpha=0, double=static_cast< double >(TMath::Pi()/2.)) override
	Get position resolution^2 corresponding to the drift length from getDriftLength of this class.

Private Attributes
bool	m_useInWirePropagationDelay
	Flag to activate the propagation delay of the sense wire.

bool	m_realData = false
	Flag to distinguish betw.

StoreObjPtr< EventT0 >	m_eventTimeStoreObject
	Event timing.

const CDCGeoControlPar &	m_gcp
	Cached reference to CDC GeoControlPar object.

const CDCSimControlPar &	m_scp
	Cached reference to CDC SimControlPar object.

const CDCGeometryPar &	m_cdcp
	Cached reference to CDC GeometryPar object.

const double	m_tdcBinWidth
	Cached TDC bin width (ns).

double	m_fudgeFactor = 1.
	Cached fudge factor for space resol.

Detailed Description

Translator mirroring the realistic Digitization.

Definition at line 25 of file RealisticTDCCountTranslator.h.

Constructor & Destructor Documentation

◆ RealisticTDCCountTranslator()

RealisticTDCCountTranslator ( bool useInWirePropagationDelay = false )

explicit

Constructor, with the additional information, if propagation in the wire shall be considered.

Definition at line 17 of file RealisticTDCCountTranslator.cc.

                                                                                       :
  m_useInWirePropagationDelay(useInWirePropagationDelay), m_gcp(CDCGeoControlPar::getInstance()),
  m_scp(CDCSimControlPar::getInstance()), m_cdcp(CDCGeometryPar::Instance()),
  m_tdcBinWidth(m_cdcp.getTdcBinWidth())
{
  m_realData = !Environment::Instance().isMC();
  if (m_realData) {
    m_fudgeFactor = m_cdcp.getFudgeFactorForSigma(0);
  } else {
    m_fudgeFactor = m_cdcp.getFudgeFactorForSigma(1);
  }
  //  B2INFO("RealisticTDCCountTranslator:: m_fudgeFactor= " << m_fudgeFactor);
 
#if defined(CDC_DEBUG)
  cout << " " << endl;
  cout << "RealisticTDCCountTranslator constructor" << endl;
  cout << "m_cdcp=" << &m_cdcp << endl;
  cout << "m_tdcBinWidth=" << m_tdcBinWidth << endl;
#endif
}

◆ ~RealisticTDCCountTranslator()

~RealisticTDCCountTranslator ( )

inline

Destructor.

Definition at line 32 of file RealisticTDCCountTranslator.h.

32{};

Member Function Documentation

◆ getDriftLength()

double getDriftLength	(	unsigned short	tdcCount,
		const WireID &	wireID = `WireID()`,
		double	timeOfFlightEstimator = `0`,
		bool	leftRight = `false`,
		double	z = `0`,
		double	alpha = `0`,
		double	theta = `static_cast<double>(TMath::Pi() / 2.)`,
		unsigned short	adcCount = `0`
	)

overridevirtual

Get Drift length.

Parameters

tdcCount	TDC count (ns).
wireID	Encoded sense wire ID.
timeOfFlightEstimator	Time of flight (ns).
leftRight	left/right flag.
z	z-position on the wire (cm).
alpha	Track incident angle in r-phi plane (rad).
theta	Track incident angle in s-z plane (=polar angle) (rad).
adcCount	ADC count.

Returns: Drift length (cm) if drift time >= 0; v*(drift time) (i.e. negative value) otherwise, where v is set to one (cm/ns) for simplicity.

Implements TDCCountTranslatorBase.

Definition at line 92 of file RealisticTDCCountTranslator.cc.

{
  const double driftTime = getDriftTime(tdcCount, wireID, timeOfFlightEstimator, z, adcCount);
 
  unsigned short layer = wireID.getICLayer();
 
  //Now we have an estimate for the time it took from the ionisation to the hitting of the wire.
  //Need to reverse calculate the relation between drift lenght and drift time.
  //  double driftL = std::copysign(m_cdcp.getDriftLength(fabs(driftTime), layer, leftRight, alpha, theta), driftTime);
  //Note: The above treatment for negative drifttime is now done in m_cdcp.getDriftLength, so the line is commented out
  double driftL = m_cdcp.getDriftLength(driftTime, layer, leftRight, alpha, theta);
 
#if defined(CDC_DEBUG)
  cout << " " << endl;
  cout << "RealisticTDCCountTranslator::getDriftLength" << endl;
  cout << "driftTime=" << driftTime << endl;
  cout << "layer=" << layer << endl;
  cout << "leftright=" << leftRight << endl;
  cout << "driftL= " << driftL << endl;
#endif
 
  return driftL;
}

◆ getDriftLengthResolution()

double getDriftLengthResolution	(	double	driftLength,
		const WireID &	wireID = `WireID()`,
		bool	leftRight = `false`,
		double	z = `0`,
		double	alpha = `0`,
		double	theta = `static_cast<double>(TMath::Pi() / 2.)`
	)

overridevirtual

Get position resolution^2 corresponding to the drift length from getDriftLength of this class.

this function returns the variance that is used as the CDC measurment resolution in track fitting

Parameters

driftLength	Drift length (cm).
wireID	Encoded sense wire ID.
leftRight	Left/right flag.
z	z-position on the wire (cm).
alpha	Track incident angle in r-phi plane (rad).
theta	Track incident angle in s-z plane (=polar angle) (rad).

Returns: Uncertainty^2 (cm^2) on the drift length.

Implements TDCCountTranslatorBase.

Definition at line 126 of file RealisticTDCCountTranslator.cc.

{
  static_cast<void>(z); //just to suppress warning of unused
  double resol = m_fudgeFactor * m_cdcp.getSigma(driftLength, wireID.getICLayer(), leftRight, alpha, theta);
  //  B2DEBUG(29, "fudgeFactor in TDCTranslator= " << m_fudgeFactor);
 
#if defined(CDC_DEBUG)
  cout << " " << endl;
  cout << "RealisticTDCCountTranslator::getDriftLengthResolution" << endl;
  cout << "spaceResol= " << resol << endl;
#endif
 
  return resol * resol;;
}

◆ getDriftTime()

double getDriftTime	(	unsigned short	tdcCount,
		const WireID &	wireID,
		double	timeOfFlightEstimator,
		double	z,
		unsigned short	adcCount
	)

overridevirtual

Get Drift time.

Parameters

tdcCount	TDC count (ns).
wireID	Encoded sense wire ID.
timeOfFlightEstimator	Time of flight (ns).
z	z-position on the wire (cm).
adcCount	ADC count.

Returns: Drift time (ns)

Implements TDCCountTranslatorBase.

Definition at line 39 of file RealisticTDCCountTranslator.cc.

{
  // translate TDC Count into time information:
  // N.B. No correction (+ or -0.5 count) is needed in the translation since no bias is in the real TDC count on average (info. from KEK electronics division).
  double driftTime = m_cdcp.getT0(wireID) - tdcCount * m_tdcBinWidth;
  //  std::cout << "t0,tdcbinw= " << m_cdcp.getT0(wireID) <<" "<< m_tdcBinWidth << std::endl;
 
  unsigned short layer = wireID.getICLayer();
 
  // Need to undo everything the digitization does in reverse order.
  // First: Undo propagation in wire, if it was used:
  if (m_useInWirePropagationDelay) {
    const B2Vector3D& backWirePos = m_cdcp.wireBackwardPosition(wireID, CDCGeometryPar::c_Aligned);
    const B2Vector3D& diffWirePos = m_cdcp.wireForwardPosition(wireID, CDCGeometryPar::c_Aligned) - backWirePos;
    //subtract distance divided by speed of electric signal in the wire from the drift time.
    //    std::cout << layer <<" "<< diffWirePos.Z() <<" "<< stereoAngleFactor << std::endl;
    double propLength = z - backWirePos.Z();
    double dZ = diffWirePos.Z();
    if (dZ != 0.) {
      propLength *= diffWirePos.Mag() / dZ;
    }
    if (m_gcp.getSenseWireZposMode() == 1) {
      //      std::cout <<"layer,zb,dzb,zf,dzf= "<< layer <<" "<< zb <<" "<< m_cdcp.getBwdDeltaZ(layer) <<" "<< m_cdcp.wireForwardPosition(wireID, CDCGeometryPar::c_Aligned).Z() <<" "<< m_cdcp.getFwdDeltaZ(layer) << std::endl;
      propLength += m_cdcp.getBwdDeltaZ(layer);
    }
    driftTime -= propLength * m_cdcp.getPropSpeedInv(layer);
  }
 
  // Second: correct for event time. If this wasn't simulated, m_eventTime can just be set to 0.
  if (m_eventTimeStoreObject.isValid() && m_eventTimeStoreObject->hasEventT0()) {
    driftTime -= m_eventTimeStoreObject->getEventT0();
  }
 
  //Third: If time of flight was simulated, this has to be undone, too. If it wasn't timeOfFlightEstimator should be taken as 0.
  driftTime -= timeOfFlightEstimator;
 
  //Forth: Time-walk correction
  if (m_realData) { //for data, always correct
    driftTime -= m_cdcp.getTimeWalk(wireID, adcCount);
    //    B2INFO("RealisticTDCCountTranslator:: time-walk corr. done.");
  } else if (m_scp.getTimeWalk()) { //for MC, ccorrect if the flag is on
    driftTime -= m_cdcp.getTimeWalk(wireID, adcCount);
    //    B2INFO("RealisticTDCCountTranslator:: time-walk corr. done for MC.");
  }
 
  return driftTime;
}