24.1.2011 Kamil Sedlak

Corrected small bug in musrPrimaryGeneratorAction (definition of default
  values of the xDirection=0, yDirection=0 and zDirection=1)

GNUmakefile

@@ -3,7 +3,7 @@
 # GNUmakefile for examples module.  Gabriele Cosmo, 06/04/98.
 # --------------------------------------------------------------
 
-name := musrSim
+name := musrSimTest
 G4TARGET := $(name)
 G4EXLIB := true
 ##LDFLAGS := $(shell root-config --glibs)

doc/musrSim.tex

@@ -683,6 +683,40 @@ Three special volumes ``Target, M0, M1 and M2''.
 
 \end{description}
 
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+\section{Optical photons}
+Normaly the simulation of ``detected signal'' stops at the level of the deposited energy in 
+a sensitive volume (e.g.\ in a scintilator tile).  However, in some special cases, one would
+like to know how the light is propagated through the scintillators.  In such case simulation
+of optical photons is possible.  It will, however, significantly (in some cases dramatically)
+reduce the speed of the simulation.  The output of the optical photon simulation is stored
+in variables starting with the string ``odet\_''.  The user has to specify several parameters
+in order to simulate optical photons:
+%
+\begin{description}
+\item{\bf /musr/command G4OpticalPhotons \emph{true}}\\
+        Specifies to musrSim whether optical photons should be simulated or not.
+	If this parameter is set to \emph{false}, anything connected with optical photons
+        will be ignored internally in musrSim, and the user therefore does not have to
+        comment out other parameters connected with optical photons in the macro file.
+
+ \item{\bf /musr/command materialPropertiesTable \emph{MPT\_name} \emph{property} \emph{n} \emph{val(1)} ... \emph{val(n)} \emph{E(1)} ... \emph{E(n)} \\
+        Defines some optical \emph{property} of a given material (e.g.\ absorption lenght, 
+	refractive index, scintillation yield, ...} to a material property table.
+        \emph{MPT\_name} stands for the material property table name.  The table has
+	\emph{n} different values \emph{val(1)} ... \emph{val(n)} for the same number
+	of optical photon energies \emph{E(1)} ... \emph{E(n)} expressed in MeV.
+	If \emph{n}=0, the \emph{property} is called ``constant property''.
+	Possible \emph{property} keywords are: ABSLENGTH, RINDEX, FASTCOMPONENT, SLOWCOMPONENT,
+	SCINTILLATIONYIELD, and constant \emph{properties} are SCINTILLATIONYIELD, RESOLUTIONSCALE, 
+	FASTTIMECONSTANT, SLOWTIMECONSTANT, and YIELDRATIO.
+
+ \item {\bf /musr/command setMaterialPropertiesTable \emph{MPT\_name} \emph{material\_name}} \\
+        Assigns a material property table defined by ``/musr/command materialPropertiesTable''
+	to a given material.
+	
+\end{description}
+
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 \section{Some other parameters}
 %

src/musrPrimaryGeneratorAction.cc

@@ -65,6 +65,7 @@ musrPrimaryGeneratorAction::musrPrimaryGeneratorAction(
    p0(0), pSigma(0), pMinAllowed(0), pMaxAllowed(1e10*mm),
    xangle0(0), yangle0(0), xangleSigma(0), yangleSigma(0), pitch(0),
    UnpolarisedMuonBeam(false), TransversalyUnpolarisedMuonBeam(false), xPolarisIni(1.), yPolarisIni(0.), zPolarisIni(0.),
+   xDirection(0), yDirection(0), zDirection(1.),
    polarisFraction(1.),
    muonDecayTimeMin(-1), muonDecayTimeMax(-1), muonMeanLife(2197.03*ns),
    takeMuonsFromTurtleFile(false), z0_InitialTurtle(0),