lem::DistanceMap Class Reference

Calculate and cache electrostatic potential between cells. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-21/LEM/func/DistanceMap.h"

Public Member Functions
double	InvDist (int aplane, int bplane, int acell, int bcell) const
double	InvDist (const LiteHit &a, const LiteHit &b) const
TH2 *	MakeHist () const

Static Public Member Functions
static const DistanceMap &	Instance ()
	Singleton. More...
static void	SetPlaneScale (double ps)
	Must call before first call to Instance() More...
static void	SetCellScale (double cs)
static void	SetDecayPower (double dp)
static void	SetExpMode (bool em=true)

Static Protected Member Functions
static double	DistFunc (double dp, double dc)

Protected Attributes
double	fInvDistMap [kNumPlanes][kNumCells]

Static Protected Attributes
static double	fgPlaneScale = 3.5
static double	fgCellScale = 10.5
static double	fgDecayPower = -0.25
static bool	fgExpMode = false

Private Member Functions
	DistanceMap ()

Detailed Description

Calculate and cache electrostatic potential between cells.

Definition at line 20 of file DistanceMap.h.

Constructor & Destructor Documentation

lem::DistanceMap::DistanceMap ( )

private

Definition at line 33 of file DistanceMap.cxx.

References plot_validation_datamc::c, confusionMatrixTree::count, DistFunc(), E, fInvDistMap, lem::kNumCells, lem::kNumPlanes, and util::sqr().

  {
    for(int p = 0; p < kNumPlanes; ++p){
      for(int c = 0; c < kNumCells; ++c){
        // For long enough distances don't bother to do the integral over cells
        if(p > 5 || c > 5){
          fInvDistMap[p][c] = DistFunc(p, c);
          continue;
        }

        double E = 0;

        const int kNumSteps = 10;
        for(int dp1 = 0; dp1 < kNumSteps; ++dp1){
          for(int dp2 = 0; dp2 < kNumSteps; ++dp2){
            const double dp = (dp1-dp2)/double(kNumSteps);
            for(int dc1 = 0; dc1 < kNumSteps; ++dc1){
              for(int dc2 = 0; dc2 < kNumSteps; ++dc2){
                // Avoid division by zero
                if(p == 0 && c == 0 && dp1 == dp2 && dc1 == dc2) continue;

                const double dc = (dc1-dc2)/double(kNumSteps);

                E += DistFunc(p+dp, c+dc);
              }
            }
          }
        } // end for dp1

        // Don't count the terms that we continue'd out of
        int count = util::sqr(util::sqr(kNumSteps));
        if(p == 0 && c == 0) count -= util::sqr(kNumSteps);

        fInvDistMap[p][c] = E/count;
      } // end for c
    } // end for p
  }

Member Function Documentation

double lem::DistanceMap::DistFunc ( double  dp, double  dc  )

staticprotected

Definition at line 86 of file DistanceMap.cxx.

References stan::math::exp(), fgCellScale, fgDecayPower, fgExpMode, fgPlaneScale, cet::pow(), and util::pythag().

Referenced by DistanceMap(), and SetExpMode().

  {
    const double dr = util::pythag(fgPlaneScale*dp, fgCellScale*dc);

    if(fgExpMode){
      // Can control the coefficient with planeScale and cellScale
      return exp(-dr);
    }

    return pow(dr, fgDecayPower);
  }

const DistanceMap & lem::DistanceMap::Instance ( )

static

Singleton.

Definition at line 26 of file DistanceMap.cxx.

Referenced by lem::FillPotential(), and lem::MatchableEvent::MatchableEvent().

  {
    static DistanceMap instance;
    return instance;
  }

double lem::DistanceMap::InvDist	(	int	aplane,
		int	bplane,
		int	acell,
		int	bcell
	)		const

Definition at line 99 of file DistanceMap.cxx.

References abs(), ana::assert(), fInvDistMap, lem::kNumCells, and lem::kNumPlanes.

Referenced by lem::FillPotential(), InvDist(), and lem::MatchableEvent::MatchableEvent().

  {
    const int dp = abs(aplane-bplane);
    const int dc = abs(acell-bcell);

    // Very far away
    if(dp < 0) return 0;
    if(dp >= kNumPlanes) return 0;
    if(dc < 0) return 0;
    if(dc >= kNumCells) return 0;

    // Comparison between different views shouldn't happen
    assert(dp%2 == 0);

    return fInvDistMap[dp][dc];
  }

double lem::DistanceMap::InvDist	(	const LiteHit &	a,
		const LiteHit &	b
	)		const

Definition at line 117 of file DistanceMap.cxx.

References lem::LiteHit::Cell(), InvDist(), and lem::LiteHit::Plane().

  {
    return InvDist(a.Plane(), b.Plane(), a.Cell(), b.Cell());
  }

TH2 * lem::DistanceMap::MakeHist

(

)

const

Definition at line 72 of file DistanceMap.cxx.

References plot_validation_datamc::c, fInvDistMap, lem::kNumCells, and lem::kNumPlanes.

  {
    TH2* hDistMap = new TH2F("", ";#Deltaplanes;#Deltacells", 5, 0, 5, 5, 0, 5);

    for(int p = 0; p < kNumPlanes; ++p){
      for(int c = 0; c < kNumCells; ++c){
        hDistMap->Fill(p, c, 1/fInvDistMap[p][c]);
      } // end for c
    } // end for p

    return hDistMap;
  }

static void lem::DistanceMap::SetCellScale ( double  cs )

inlinestatic

Definition at line 33 of file DistanceMap.h.

References lem_server::cs, and fgCellScale.

Referenced by lem::FindMatchesTranspose::reconfigure(), lem::FindLEMMatches::reconfigure(), and lem::LEM::reconfigure().

{fgCellScale = cs;}

static void lem::DistanceMap::SetDecayPower ( double  dp )

inlinestatic

Definition at line 34 of file DistanceMap.h.

References fgDecayPower.

Referenced by lem::FindMatchesTranspose::reconfigure(), lem::FindLEMMatches::reconfigure(), and lem::LEM::reconfigure().

{fgDecayPower = dp;}

static void lem::DistanceMap::SetExpMode ( bool  em = true )

inlinestatic

Definition at line 35 of file DistanceMap.h.

References DistFunc(), and fgExpMode.

Referenced by lem::FindMatchesTranspose::reconfigure(), lem::FindLEMMatches::reconfigure(), and lem::LEM::reconfigure().

{fgExpMode = em;}

static void lem::DistanceMap::SetPlaneScale ( double  ps )

inlinestatic

Must call before first call to Instance()

Definition at line 32 of file DistanceMap.h.

References fgPlaneScale, and nd_projection_maker::ps.

Referenced by lem::FindMatchesTranspose::reconfigure(), lem::FindLEMMatches::reconfigure(), and lem::LEM::reconfigure().

{fgPlaneScale = ps;}

Member Data Documentation

double lem::DistanceMap::fgCellScale = 10.5

staticprotected

Definition at line 41 of file DistanceMap.h.

Referenced by DistFunc(), and SetCellScale().

double lem::DistanceMap::fgDecayPower = -0.25

staticprotected

Definition at line 42 of file DistanceMap.h.

Referenced by DistFunc(), and SetDecayPower().

bool lem::DistanceMap::fgExpMode = false

staticprotected

Definition at line 43 of file DistanceMap.h.

Referenced by DistFunc(), and SetExpMode().

double lem::DistanceMap::fgPlaneScale = 3.5

staticprotected

Definition at line 40 of file DistanceMap.h.

Referenced by DistFunc(), and SetPlaneScale().

double lem::DistanceMap::fInvDistMap[kNumPlanes][kNumCells]

protected

Definition at line 45 of file DistanceMap.h.

Referenced by DistanceMap(), InvDist(), and MakeHist().

---

The documentation for this class was generated from the following files:

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-21/LEM/func/DistanceMap.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-21/LEM/func/DistanceMap.cxx