LArWheelSliceSolidTests.cxx File Reference

#include <iostream>
#include <stdexcept>
#include <format>
#include <map>
#include "TRandom3.h"
#include "TF1.h"
#include "TNtupleD.h"
#include "TFile.h"
#include "TROOT.h"
#include "GeoSpecialShapes/LArWheelCalculator.h"
#include "LArWheelSliceSolid.h"
#include "CLHEP/Units/SystemOfUnits.h"
#include "CLHEP/Units/PhysicalConstants.h"
#include <stdlib.h>
#include "CxxUtils/checker_macros.h"

Go to the source code of this file.

Functions
static void	get_r (const G4VSolid *p, G4double z, G4double &rmin, G4double &rmax)
double	LArWheelSliceSolid_fcn_area (double *x, double *p)
double	LArWheelSliceSolid_fcn_vol (double *x, double *p)
double	LArWheelSliceSolid_fcn_area_on_pc (double *x, double *p)
double	LArWheelSliceSolid_get_dl (double *x, double *par, G4int side)
static double	fcn_length (double *x, double *p)
double	LArWheelSliceSolid_fcn_side_area (double *x, double *p)

Variables
	ATLAS_NO_CHECK_FILE_THREAD_SAFETY
static double	IntPrecision = 0.0001
static TRandom *	rnd = 0
static std::map< double, LArWheelSliceSolid * >	solid

Function Documentation

fcn_length()

double fcn_length ( double * x, double * p )

static

Definition at line 591 of file LArWheelSliceSolidTests.cxx.

{
  return LArWheelSliceSolid_get_dl(x, p, 1) + LArWheelSliceSolid_get_dl(x, p, -1);
}

get_r()

void get_r ( const G4VSolid * p, G4double z, G4double & rmin, G4double & rmax )

static

Definition at line 58 of file LArWheelSliceSolidTests.cxx.

{
    G4ThreeVector from(10.*CLHEP::m, 0., z);
    rmax = from[0] - p->DistanceToIn(from, G4ThreeVector(-1., 0., 0.));
    from[0] = 0.;
    rmin = p->DistanceToIn(from, G4ThreeVector(1., 0., 0.));
}

LArWheelSliceSolid_fcn_area()

double LArWheelSliceSolid_fcn_area	(	double *	x,
		double *	p )

Definition at line 524 of file LArWheelSliceSolidTests.cxx.

{
  const double &z = x[0];
  const double &r = p[0];
  const double &index = p[1];
 
  G4ThreeVector a(0., r, z);
  double b = solid[index]->GetCalculator()->AmplitudeOfSurface(a, -1, 121) // sagging ignored, effect should be negligible, use arbitrary fan number
    - solid[index]->GetCalculator()->AmplitudeOfSurface(a, 1, 121);
  return b;
}

LArWheelSliceSolid_fcn_area_on_pc()

double LArWheelSliceSolid_fcn_area_on_pc	(	double *	x,
		double *	p )

Definition at line 544 of file LArWheelSliceSolidTests.cxx.

{
  const double &z = x[0];
  const double &index = p[0];
 
  G4double rmin, rmax;
  get_r(solid[index]->m_BoundingShape, z, rmin, rmax);
 
  double result = 0.;
  G4ThreeVector a(0., rmin, z);
  result += solid[index]->GetCalculator()->AmplitudeOfSurface(a, -1, 232) // sagging ignored, effect should be negligible, use arbitrary fan number
    - solid[index]->GetCalculator()->AmplitudeOfSurface(a, 1, 232);
  a[1] = rmax;
  result += solid[index]->GetCalculator()->AmplitudeOfSurface(a, -1, 343)
    - solid[index]->GetCalculator()->AmplitudeOfSurface(a, 1, 343);
 
  return result;
}

LArWheelSliceSolid_fcn_side_area()

double LArWheelSliceSolid_fcn_side_area	(	double *	x,
		double *	p )

Definition at line 596 of file LArWheelSliceSolidTests.cxx.

{
  const double &r = x[0];
  const double &index = p[0];
 
  return solid[index]->get_length_at_r(r);
}

LArWheelSliceSolid_fcn_vol()

double LArWheelSliceSolid_fcn_vol	(	double *	x,
		double *	p )

Definition at line 536 of file LArWheelSliceSolidTests.cxx.

{
  const double &r = x[0];
  const double &index = p[0];
 
  return solid[index]->get_area_at_r(r);
}

LArWheelSliceSolid_get_dl()

double LArWheelSliceSolid_get_dl	(	double *	x,
		double *	par,
		G4int	side )

Definition at line 564 of file LArWheelSliceSolidTests.cxx.

{
  const double &z = x[0];
  const double &r = par[0];
  const double &index = par[1];
 
  const double h = 0.001;
 
  //check what happens if z+h > m_Zmax etc
  G4ThreeVector p(0., r, z + h);
  G4double D1 = solid[index]->GetCalculator()->AmplitudeOfSurface(p, side, 5665);  // sagging ignored, effect should be negligible, use arbitrary fan number
  p[2] = z - h;
  D1 -= solid[index]->GetCalculator()->AmplitudeOfSurface(p, side, 5665);
  D1 /= 2 * h;
 
  p[2] = z + h / 2;
  G4double D2 = solid[index]->GetCalculator()->AmplitudeOfSurface(p, side, 5665);
  p[2] = z - h / 2;
  D2 -= solid[index]->GetCalculator()->AmplitudeOfSurface(p, side, 5665);
  D2 /= h;
 
  G4double D = (D2 * 4 - D1) / 3.;
  G4double dl = sqrt(1 + D * D);
 
  return dl;
}

Variable Documentation

ATLAS_NO_CHECK_FILE_THREAD_SAFETY

ATLAS_NO_CHECK_FILE_THREAD_SAFETY

Definition at line 25 of file LArWheelSliceSolidTests.cxx.

IntPrecision

double IntPrecision = 0.0001

static

Definition at line 27 of file LArWheelSliceSolidTests.cxx.

rnd

TRandom* rnd = 0

static

Definition at line 68 of file LArWheelSliceSolidTests.cxx.

solid

std::map<double, LArWheelSliceSolid *> solid

static

Definition at line 522 of file LArWheelSliceSolidTests.cxx.