Utils Namespace Reference

Functions
double	deltaPhi (double phi1, double phi2)
double	deltaR2 (double eta1, double eta2, double phi1, double phi2)
double	deltaR (double eta1, double eta2, double phi1, double phi2)
void	findNearestPoint (const Amg::Vector3D &inputPos, const Trk::CaloExtension *caloExtension, int &nearestIdx, Amg::Vector3D &nearestPos, Amg::Vector3D &nearestMom)
double	z0wrtPV (const xAOD::TrackParticle *trk, const xAOD::Vertex *vtx)
	Provide the trk DCA w.r.t. the PV. More...
const xAOD::Vertex *	selectPV (SG::ReadHandle< xAOD::VertexContainer > &container)
	Finds the Primary Vertex. More...
const xAOD::Jet *	findLeadingJet (SG::ReadHandle< xAOD::JetContainer > &container)
	Find leading jet. More...

Function Documentation

deltaPhi()

double Utils::deltaPhi	(	double	phi1,
		double	phi2
	)

Definition at line 11 of file CaloCellSelectorUtils.cxx.

                                              {
        double dPhi = std::fabs(phi1 - phi2);
        if (dPhi > M_PI) dPhi = 2 * M_PI - dPhi;
        return dPhi;
    }

deltaR()

double Utils::deltaR	(	double	eta1,
		double	eta2,
		double	phi1,
		double	phi2
	)

Definition at line 23 of file CaloCellSelectorUtils.cxx.

{ return std::sqrt(deltaR2(eta1, eta2, phi1, phi2)); }

deltaR2()

double Utils::deltaR2	(	double	eta1,
		double	eta2,
		double	phi1,
		double	phi2
	)

Definition at line 17 of file CaloCellSelectorUtils.cxx.

                                                                       {
        double dPhi = deltaPhi(phi1, phi2);
        double dEta = std::fabs(eta1 - eta2);
        return dEta * dEta + dPhi * dPhi;
    }

findLeadingJet()

const xAOD::Jet * Utils::findLeadingJet

(

SG::ReadHandle< xAOD::JetContainer > &

container

)

Find leading jet.

Definition at line 26 of file Trigger/TrigMonitoring/TrigMinBiasMonitoring/src/Utils.cxx.

                                                                         {
  auto compareJetPt = [](const auto& jet1, const auto& jet2) { return jet1->pt() < jet2->pt(); };
 
  const auto leadingJet = std::max_element(container->begin(), container->end(), compareJetPt);
  if (leadingJet == container->end()) {
    return nullptr;
  }
 
  return *leadingJet;
}

findNearestPoint()

void Utils::findNearestPoint	(	const Amg::Vector3D &	inputPos,
		const Trk::CaloExtension *	caloExtension,
		int &	nearestIdx,
		Amg::Vector3D &	nearestPos,
		Amg::Vector3D &	nearestMom
	)

Definition at line 25 of file CaloCellSelectorUtils.cxx.

                                                                            {
        const std::vector<Trk::CurvilinearParameters>& intersections = caloExtension->caloLayerIntersections();
        int nPts = intersections.size();
 
        int idxL, idxR, idxMid;
        Amg::Vector3D pos(0, 0, 0), mom(0, 0, 0);  // initialization should be(?) unnecessary, just to suppress scary warning msg...
        Amg::Vector3D posL(0, 0, 0), momL(0, 0, 0);
        Amg::Vector3D posR(0, 0, 0), momR(0, 0, 0);
 
        // find nearest crossing point
        idxL = 0;
        idxR = nPts - 1;
        while ((idxR - idxL) > 1) {
            idxMid = (idxL + idxR) / 2;
            pos = intersections[idxMid].position();
            mom = intersections[idxMid].momentum();
 
            if ((inputPos - pos).dot(mom) > 0) {
                idxL = idxMid;
                std::swap(posL, pos);
                std::swap(momL, mom);
            } else {
                idxR = idxMid;
                std::swap(posR, pos);
                std::swap(momR, mom);
            }
        }
 
        if (idxL == 0) {
            posL = intersections[0].position();
            momL = intersections[0].momentum();
        }
        if (idxR == (nPts - 1)) {
            posR = intersections[nPts - 1].position();
            momR = intersections[nPts - 1].momentum();
        }
 
        float mag2L = (inputPos - posL).mag2();
        float mag2R = (inputPos - posR).mag2();
        nearestIdx = (mag2L < mag2R) ? idxL : idxR;
        nearestPos = (mag2L < mag2R) ? posL : posR;
        nearestMom = (mag2L < mag2R) ? momL : momR;
    }

selectPV()

const xAOD::Vertex * Utils::selectPV

(

SG::ReadHandle< xAOD::VertexContainer > &

container

)

Finds the Primary Vertex.

Definition at line 12 of file Trigger/TrigMonitoring/TrigMinBiasMonitoring/src/Utils.cxx.

                                                                         {
  if (not container.isValid()) {
    return nullptr;
  }
 
  for (const xAOD::Vertex* vtx : *container) {
    if (vtx->vertexType() == xAOD::VxType::PriVtx) {
      return vtx;
    }
  }
 
  return nullptr;
}

z0wrtPV()

double Utils::z0wrtPV	(	const xAOD::TrackParticle *	trk,
		const xAOD::Vertex *	vtx
	)

Provide the trk DCA w.r.t. the PV.

Definition at line 8 of file Trigger/TrigMonitoring/TrigMinBiasMonitoring/src/Utils.cxx.

                                                                    {
  return (trk->z0() + trk->vz() - vtx->z()) * std::sin(trk->theta());
}