Project

General

Profile

Task #14592

Milestone #14454: Refactoring LArG4

4. provide scripts to provide proper input of optical properties (Rayleigh scattering length, refraction index, absorption length) to Geant4.

Added by Katherine Lato over 3 years ago. Updated almost 3 years ago.

Status:
Closed
Priority:
Normal
Category:
Simulation
Target version:
-
Start date:
11/22/2016
Due date:
11/30/2016
% Done:

100%

Estimated time:
5.00 h
Spent time:
Duration: 9

Description

4. provide scripts to provide proper input of optical properties (Rayleigh scattering
length, refraction index, absorption length) to Geant4. [DONE]
Needs to be placed in LArSoft.

History

#1 Updated by Katherine Lato over 3 years ago

  • Estimated time set to 5.00 h

#2 Updated by Hans-Joachim Wenzel over 3 years ago

  • Status changed from Assigned to Resolved
  • % Done changed from 90 to 100

modified scripts to work with root 6. Then created scripts directory in LArG4 and put LAr.C there:

//-----------------------------------------------------------------------------------------------------------------
// Formulas and numbers are based on
// arXiv:1502.04213
// and private communication with one of the authors Emily Grace:
//
//
// to run the script from the root command line do:
// root [0] .L LAr.C++
// root [1] init(); // initialize
// to make some plots do:
// root [2] sellmeierLAr(); // plot refraction index between 110 and 700 nm
// root [3] sellmeierLAr(125.,130.); // plot refraction index between 125 and 130 nm
// root [4] rayleigh(); // plot rayleigh length between 110 and 400 nm
// root [5] rayleigh(125,130); // plot rayleigh length between 110 and 400 nm
//
// to print out table that can be cut and paste into the Geant4 gdml description:
// root [6] rindextable(); // refraction index between 110 and 700 nm in nsteps=100 for T=83.81 K
// root [7] rindextable(200,300,10,3); // refraction index between 200 and 300 nm in nsteps=10 for T=90 K
// root [8] rayleightable(); // rayleigh length between 110 and 700 nm in nsteps=100 for T=83.81 K
// root [9] rayleightable(200,300,10,3); // rayleigh length between 200 and 300 nm in nsteps=10 for T=90 K
//------------------------------------------------------------------------------------------------------------------

#3 Updated by Hans-Joachim Wenzel over 3 years ago

Besides the scripts that fill the arrays we will also add the functions that calculate various optical properties of liquid argon. Some are taken from:https://indico.cern.ch/event/44566/contributions/1101918/attachments/943057/1337650/dipompeo.pdf

#4 Updated by Katherine Lato almost 3 years ago

  • Status changed from Resolved to Closed


Also available in: Atom PDF