nano_mu::DTTrigGeomUtils Class Reference

#include <MuNtupleUtils.h>

Classes
struct	chambCoord

Public Member Functions
	DTTrigGeomUtils (edm::ConsumesCollector &&collector, bool dirInDeg=true)
	Constructor. More...
void	getFromES (const edm::Run &run, const edm::EventSetup &environment)
	Update EventSetup information. More...
bool	hasPosRF (int wh, int sec)
	Checks id the chamber has positive RF;. More...
chambCoord	trigToReco (const L1MuDTChambPhDigi *trig)
	Return local position and direction in chamber RF - legacy. More...
chambCoord	trigToReco (const L1Phase2MuDTPhDigi *trig)
	Return local position and direction in chamber RF - analytical method. More...

Private Attributes
ESTokenHandle< DTGeometry, MuonGeometryRecord, edm::Transition::BeginRun >	m_dtGeom
std::array< double, 4 >	m_zcn
std::array< double, 4 >	m_zsl1
std::array< double, 4 >	m_zsl3

Static Private Attributes
static constexpr double	PH1_PHI_R = 4096.
static constexpr double	PH1_PHIB_R = 512.
static constexpr double	PH2_PHI_R = 65536. / 0.8
static constexpr double	PH2_PHIB_R = 2048. / 1.4

Detailed Description

Definition at line 86 of file MuNtupleUtils.h.

Constructor & Destructor Documentation

DTTrigGeomUtils()

DTTrigGeomUtils::DTTrigGeomUtils	(	edm::ConsumesCollector &&	collector,
		bool	dirInDeg = true
	)

Constructor.

Definition at line 19 of file MuNtupleUtils.cc.

References eostools::move().

    : m_dtGeom{std::move(collector), "idealForDigi"} {}

Member Function Documentation

getFromES()

void nano_mu::DTTrigGeomUtils::getFromES ( const edm::Run &  run, const edm::EventSetup &  environment  )

inline

Update EventSetup information.

Definition at line 106 of file MuNtupleUtils.h.

References DTGeometry::chamber(), reco_calib_source_client_cfg::environment, nano_mu::ESTokenHandle< T, R, TR >::getFromES(), m_dtGeom, m_zcn, m_zsl1, m_zsl3, GeomDet::position(), DTChamber::superLayer(), GeomDet::surface(), and GloballyPositioned< T >::toLocal().

Referenced by MuDTTPGPhiFlatTableProducer::getFromES().

                                                                        {
      m_dtGeom.getFromES(environment);
      for (int i_st = 0; i_st != 4; ++i_st) {
        const DTChamberId chId(-2, i_st + 1, 4);
        const DTChamber* chamb = m_dtGeom->chamber(chId);
        const DTSuperLayer* sl1 = chamb->superLayer(DTSuperLayerId(chId, 1));
        const DTSuperLayer* sl3 = chamb->superLayer(DTSuperLayerId(chId, 3));
        m_zsl1[i_st] = chamb->surface().toLocal(sl1->position()).z();
        m_zsl3[i_st] = chamb->surface().toLocal(sl3->position()).z();
        m_zcn[i_st] = 0.5 * (m_zsl1[i_st] + m_zsl3[i_st]);
      }
    };

hasPosRF()

bool nano_mu::DTTrigGeomUtils::hasPosRF ( int  wh, int  sec  )

inline

Checks id the chamber has positive RF;.

Definition at line 103 of file MuNtupleUtils.h.

References fileinputsource_cfi::sec.

{ return wh > 0 || (wh == 0 && sec % 4 > 1); };

trigToReco() [1/2]

nano_mu::DTTrigGeomUtils::chambCoord DTTrigGeomUtils::trigToReco

(

const L1MuDTChambPhDigi *

trig

)

Return local position and direction in chamber RF - legacy.

Definition at line 22 of file MuNtupleUtils.cc.

References funct::cos(), HLT_2024v13_cff::delta_phi, DeadROC_duringRun::dir, Geom::pi(), fileinputsource_cfi::sec, funct::sin(), funct::tan(), ClusterTask_cfi::trig, and x.

Referenced by MuDTTPGPhiFlatTableProducer::fillTable().

                                                                                                 {
  auto wh{trig->whNum()};
  auto sec{trig->scNum() + 1};
  auto st{trig->stNum()};
  auto phi{trig->phi()};
  auto phib{trig->phiB()};

  auto recoChamb = [&]() {
    if (st != 4) {
      return DTChamberId(wh, st, sec);
    }
    int reco_sec{(sec == 4 && phi > 0) ? 13 : (sec == 10 && phi > 0) ? 14 : sec};
    return DTChamberId(wh, st, reco_sec);
  };

  auto gpos{m_dtGeom->chamber(recoChamb())->position()};
  auto r{gpos.perp()};

  auto delta_phi{gpos.phi() - (sec - 1) * Geom::pi() / 6};

  // zcn is in local coordinates -> z invreases approching to vertex
  // LG: zcn offset was removed <- CB do we need to fix this?
  float x = r * tan((phi - (phi < 0 ? 1 : 0)) / PH1_PHI_R) * cos(delta_phi) - r * sin(delta_phi);
  float dir = (phib / PH1_PHIB_R + phi / PH1_PHI_R);

  // change sign in case of positive wheels
  if (hasPosRF(wh, sec)) {
    x = -x;
  } else {
    dir = -dir;
  }

  return {x, dir};
}

trigToReco() [2/2]

nano_mu::DTTrigGeomUtils::chambCoord DTTrigGeomUtils::trigToReco

(

const L1Phase2MuDTPhDigi *

trig

)

Return local position and direction in chamber RF - analytical method.

Definition at line 57 of file MuNtupleUtils.cc.

References funct::cos(), HLT_2024v13_cff::delta_phi, DeadROC_duringRun::dir, Geom::pi(), quality, fileinputsource_cfi::sec, funct::sin(), funct::tan(), ClusterTask_cfi::trig, and x.

                                                                                                  {
  auto wh{trig->whNum()};
  auto sec{trig->scNum() + 1};
  auto st{trig->stNum()};
  auto phi{trig->phi()};
  auto phib{trig->phiBend()};
  auto quality{trig->quality()};
  auto sl{trig->slNum()};

  auto recoChamb = [&]() {
    if (st != 4) {
      return DTChamberId(wh, st, sec);
    }
    int reco_sec{(sec == 4 && phi > 0) ? 13 : (sec == 10 && phi > 0) ? 14 : sec};
    return DTChamberId(wh, st, reco_sec);
  };

  auto gpos{m_dtGeom->chamber(recoChamb())->position()};
  auto r{gpos.perp()};

  auto delta_phi{gpos.phi() - (sec - 1) * Geom::pi() / 6};

  // CB to be potentially updated based on Silvia's results
  double zRF = 0;
  if (quality >= 6 && quality != 7)
    zRF = m_zcn[st - 1];
  if ((quality < 6 || quality == 7) && sl == 1)
    zRF = m_zsl1[st - 1];
  if ((quality < 6 || quality == 7) && sl == 3)
    zRF = m_zsl3[st - 1];

  // zcn is in local coordinates -> z invreases approching to vertex
  // LG: zcn offset was removed <- CB Mist confirm it is tryly accurate?
  float x = r * tan((phi - (phi < 0 ? 1 : 0)) / PH1_PHI_R) * (cos(delta_phi) - zRF) - r * sin(delta_phi);
  float dir = (phib / PH2_PHIB_R + phi / PH2_PHI_R);

  // change sign in case of positive wheels
  if (hasPosRF(wh, sec)) {
    x = -x;
  } else {
    dir = -dir;
  }

  return {x, dir};
}

Member Data Documentation

m_dtGeom

ESTokenHandle<DTGeometry, MuonGeometryRecord, edm::Transition::BeginRun> nano_mu::DTTrigGeomUtils::m_dtGeom

private

Definition at line 117 of file MuNtupleUtils.h.

Referenced by getFromES().

m_zcn

std::array<double, 4> nano_mu::DTTrigGeomUtils::m_zcn

private

Definition at line 122 of file MuNtupleUtils.h.

Referenced by getFromES().

m_zsl1

std::array<double, 4> nano_mu::DTTrigGeomUtils::m_zsl1

private

Definition at line 123 of file MuNtupleUtils.h.

Referenced by getFromES().

m_zsl3

std::array<double, 4> nano_mu::DTTrigGeomUtils::m_zsl3

private

Definition at line 124 of file MuNtupleUtils.h.

Referenced by getFromES().

PH1_PHI_R

constexpr double nano_mu::DTTrigGeomUtils::PH1_PHI_R = 4096.

staticprivate

Definition at line 126 of file MuNtupleUtils.h.

PH1_PHIB_R

constexpr double nano_mu::DTTrigGeomUtils::PH1_PHIB_R = 512.

staticprivate

Definition at line 127 of file MuNtupleUtils.h.

PH2_PHI_R

constexpr double nano_mu::DTTrigGeomUtils::PH2_PHI_R = 65536. / 0.8

staticprivate

Definition at line 129 of file MuNtupleUtils.h.

PH2_PHIB_R

constexpr double nano_mu::DTTrigGeomUtils::PH2_PHIB_R = 2048. / 1.4

staticprivate

Definition at line 130 of file MuNtupleUtils.h.