DriftTimeUtil.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#pragma once
 
#include <cdc/geometry/CDCGeometryPar.h>
#include <TRandom.h>
#include <cdc/geometry/CDCGeoControlPar.h>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    struct DriftTimeUtil {
 
      static double getDriftV(double driftTime,
                              unsigned short iCLayer,
                              unsigned short lr,
                              double alpha = 0.,
                              double theta = 0.5 * M_PI)
      {
        return CDC::CDCGeometryPar::Instance().getDriftV(driftTime, iCLayer, lr, alpha, theta);
      }
 
      static double getDriftLength(double driftTime,
                                   unsigned short iCLayer,
                                   unsigned short lr,
                                   double alpha = 0.,
                                   double theta = 0.5 * M_PI)
      {
        return CDC::CDCGeometryPar::Instance().getDriftLength(driftTime, iCLayer, lr, alpha, theta);
      }
 
      static double getDriftTime(double dist,
                                 unsigned short iCLayer,
                                 unsigned short lr,
                                 double alpha,
                                 double theta)
      {
        return CDC::CDCGeometryPar::Instance().getDriftTime(dist, iCLayer, lr, alpha, theta);
      }
 
      static double getPropTime(const WireID& wireID, double z)
      {
        unsigned short iCLayer = wireID.getICLayer();
        const CDC::CDCGeometryPar& geometryPar = CDC::CDCGeometryPar::Instance();
        const CDC::CDCGeoControlPar& controlPar = CDC::CDCGeoControlPar::getInstance();
        B2Vector3D backwardWirePos =
          geometryPar.wireBackwardPosition(wireID, CDC::CDCGeometryPar::c_Aligned);
        B2Vector3D forwardWirePos =
          geometryPar.wireForwardPosition(wireID, CDC::CDCGeometryPar::c_Aligned);
        double zDistance = z - backwardWirePos.Z();
 
        // Actually are very small correction
        double stereoFactor =
          (forwardWirePos - backwardWirePos).Mag() / (forwardWirePos.Z() - backwardWirePos.Z());
 
        double distance = zDistance * stereoFactor;
        if (controlPar.getSenseWireZposMode() == 1) {
          distance += geometryPar.getBwdDeltaZ(iCLayer);
        }
        return distance * geometryPar.getPropSpeedInv(iCLayer);
      }
 
      static double getTimeWalk(const WireID& wireID, unsigned short adcCount)
      {
        return CDC::CDCGeometryPar::Instance().getTimeWalk(wireID, adcCount);
      }
 
      static double getMeasuredTime(const WireID& wireID, unsigned short tdcCount, bool smear)
      {
        double smearing = 0;
        if (smear) {
          smearing = gRandom->Rndm() - 0.5;
        }
        double channelT0 = static_cast<double>(CDC::CDCGeometryPar::Instance().getT0(wireID));
        double measuredTime =
          channelT0 - (tdcCount + smearing) * CDC::CDCGeometryPar::Instance().getTdcBinWidth();
        if (measuredTime > 2000) {
          B2INFO("channelT0 " << channelT0);
          B2INFO("measuredTime " << measuredTime);
        }
        return measuredTime;
      }
 
    };
  }
}