LArWheelSliceSolidTests.cxx File Reference

#include <iostream>
#include <stdexcept>
#include "boost/io/ios_state.hpp"
#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
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)
double	LArWheelSliceSolid_fcn_side_area (double *x, double *p)

Variables
	ATLAS_NO_CHECK_FILE_THREAD_SAFETY

Function Documentation

LArWheelSliceSolid_fcn_area()

double LArWheelSliceSolid_fcn_area	(	double *	x,
		double *	p
	)

Definition at line 528 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 548 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 600 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 540 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 568 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.