CaloCellHelpers.h File Reference

#include <cmath>
#include "CaloEvent/CaloCell.h"
#include "CaloGeoHelpers/CaloPhiRange.h"

Include dependency graph for CaloCellHelpers.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
double	phiMean (double a, double b)
bool	crossedPhi (const CaloCell &cell, double phi_entrance, double phi_exit)
	Return true if the cell crossed was crossed by the track in phi. More...
double	getPathLengthInEta (const CaloCell &cell, double eta_entrance, double eta_exit)
	Return the % of path length crossed by the track inside a cell in eta. More...
double	getPathLengthInZ (double zMin, double zMax, double z_entrance, double z_exit)
	Return the % of path length crossed by the track inside a cell in Z. More...
double	getPathLengthInZ (const CaloCell &cell, double z_entrance, double z_exit)
	Return the % of path length crossed by the track inside a cell in Z. More...
double	pathInsideCell (const CaloCell &cell, const Amg::Vector3D &entry, const Amg::Vector3D &exit)
	Return the % of the path crossed inside the cell, given the parameters for the extrapolation at entrance and exit of the layer. More...

Function Documentation

crossedPhi()

bool crossedPhi ( const CaloCell &  cell, double  phi_entrance, double  phi_exit  )

inline

Return true if the cell crossed was crossed by the track in phi.

Definition at line 16 of file CaloCellHelpers.h.

                                                                                   {
    double mean_phi = phiMean(phi_entrance, phi_exit);
    double dphi = fabs(CaloPhiRange::diff(phi_entrance, phi_exit));
    double phi_min = mean_phi - dphi, phi_max = mean_phi + dphi;
 
    return (CaloPhiRange::diff(cell.phi() + cell.caloDDE()->dphi() / 2., phi_min) > 0 &&
            CaloPhiRange::diff(phi_max, cell.phi() - cell.caloDDE()->dphi() / 2.) > 0);
}

getPathLengthInEta()

double getPathLengthInEta ( const CaloCell &  cell, double  eta_entrance, double  eta_exit  )

inline

Return the % of path length crossed by the track inside a cell in eta.

Definition at line 26 of file CaloCellHelpers.h.

                                                                                             {
    double etaMin = cell.eta() - 0.5 * cell.caloDDE()->deta();
    double etaMax = cell.eta() + 0.5 * cell.caloDDE()->deta();
    if (fabs(eta_entrance - eta_exit) < 1e-6)  // to avoid FPE
        return eta_entrance > etaMin && eta_entrance < etaMax;
 
    double etaMinTrack = std::min(eta_entrance, eta_exit);
    double etaMaxTrack = std::max(eta_entrance, eta_exit);
    return (std::min(etaMax, etaMaxTrack) - std::max(etaMin, etaMinTrack)) / (etaMaxTrack - etaMinTrack);
}

getPathLengthInZ() [1/2]

double getPathLengthInZ ( const CaloCell &  cell, double  z_entrance, double  z_exit  )

inline

Return the % of path length crossed by the track inside a cell in Z.

Definition at line 48 of file CaloCellHelpers.h.

                                                                                       {
    return getPathLengthInZ(cell.z() - 0.5 * cell.caloDDE()->dz(), cell.z() + 0.5 * cell.caloDDE()->dz(), z_entrance, z_exit);
}

getPathLengthInZ() [2/2]

double getPathLengthInZ ( double  zMin, double  zMax, double  z_entrance, double  z_exit  )

inline

Return the % of path length crossed by the track inside a cell in Z.

Definition at line 38 of file CaloCellHelpers.h.

                                                                                           {
    if (fabs(z_entrance - z_exit) < 1e-6)  // to avoid FPE
        return z_entrance > zMin && z_entrance < zMax;
 
    double zMinTrack = std::min(z_entrance, z_exit);
    double zMaxTrack = std::max(z_entrance, z_exit);
    return (std::min(zMax, zMaxTrack) - std::max(zMin, zMinTrack)) / (zMaxTrack - zMinTrack);
}

pathInsideCell()

double pathInsideCell ( const CaloCell &  cell, const Amg::Vector3D &  entry, const Amg::Vector3D &  exit  )

inline

Return the % of the path crossed inside the cell, given the parameters for the extrapolation at entrance and exit of the layer.

Definition at line 84 of file CaloCellHelpers.h.

                                                                                                      {
    if (!crossedPhi(cell, entry.phi(), exit.phi())) return 0;
    double pathCrossed = 0;
    if (cell.caloDDE()->getSubCalo() != CaloCell_ID::TILE)
        pathCrossed = getPathLengthInEta(cell, entry.eta(), exit.eta());
    else if (cell.caloDDE()->getSampling() == CaloCell_ID::TileBar1)  // ladder shape cells, divide the problem in 2
        pathCrossed = 0;                                              // getPathLengthInZTileBar1(cell, entry.z(), exit.z());
    else
        pathCrossed = getPathLengthInZ(cell, entry.z(), exit.z());
    if (pathCrossed <= 0) return 0;
    return pathCrossed;
}

phiMean()

double phiMean ( double  a, double  b  )

inline

Definition at line 13 of file CaloCellHelpers.h.

{ return 0.5 * (a + b) + (a * b < 0) * M_PI; }