diff --git a/doc/musrSim.pdf b/doc/musrSim.pdf
index f83e101..f17dd13 100644
Binary files a/doc/musrSim.pdf and b/doc/musrSim.pdf differ
diff --git a/doc/musrSim.tex b/doc/musrSim.tex
index 80e2413..f1498c1 100644
--- a/doc/musrSim.tex
+++ b/doc/musrSim.tex
@@ -566,6 +566,11 @@ care.
 	See the file musrPhysicsList.cc for the list of defined processes (e.g. G4eMultipleScattering,
 	G4eIonisation, ...) and Geant4 manual for the detail description of the processes.
 
+\item{\bf /musr/command process SetLossFluctuations\_OFF \emph{particle} \emph{process}}\\
+        Switches off the statistical fluctuations in the ionisation losses.  Can be usefull
+	for some technical studies.  Presently works only for \emph{G4eIonisation} and
+	for \emph{G4MuIonisation}.
+
 %	There is one special process, combined from G4MultipleScattering and G4CoulombScattering,
 %	defined by the following command:\\
 %{\bf /musr/command process addProcess  \emph{particle} MultipleAndCoulombScattering \emph{ordAtRestDoIt} \emph{ordAlongSteptDoIt} \emph{ordPostStepDoIt} \emph{G4Region1}  [\emph{G4Region2}]  [\emph{G4Region3}]}\\
diff --git a/src/musrPhysicsList.cc b/src/musrPhysicsList.cc
index b54154b..f7b98bb 100644
--- a/src/musrPhysicsList.cc
+++ b/src/musrPhysicsList.cc
@@ -296,7 +296,8 @@ void musrPhysicsList::ConstructEM()
       if (strcmp(tmpString1,"process")!=0)  continue; 
       
       float tmpCutValue;
-      if ((strcmp(tmpString2,"addProcess")==0)||(strcmp(tmpString2,"addDiscreteProcess")==0)||(strcmp(tmpString2,"addModel")==0)) {
+      if ((strcmp(tmpString2,"addProcess")==0)||(strcmp(tmpString2,"addDiscreteProcess")==0)||
+	  (strcmp(tmpString2,"addModel")==0)||(strcmp(tmpString2,"SetLossFluctuations_OFF")==0)) {
 	char charParticleName[100], charProcessName[100];
 	sscanf(&line[0],"%*s %*s %s %s %s",tmpString2,charParticleName,charProcessName);
 	G4cout<<"musrPhysicsList: Defining process "<<charProcessName<<" for "<<charParticleName<<G4endl;
@@ -452,8 +453,28 @@ void musrPhysicsList::ConstructEM()
 
 	}
 
+	else if (strcmp(tmpString2,"SetLossFluctuations_OFF")==0) {
+	  G4ProcessTable* processTable = G4ProcessTable::GetProcessTable();
+	  if (strcmp(charProcessName,"G4eIonisation")==0) {
+	    G4cout<<"musrPhysicsList.cc:  Switching off loss fluctuations for "<<charParticleName<<" in "<<charProcessName<<G4endl;
+	    G4eIonisation* mmm = (G4eIonisation*) processTable->FindProcess("eIoni",particleDefinition);
+	    mmm->SetLossFluctuations(false);
+	  }
+	  else if (strcmp(charProcessName,"G4MuIonisation")==0) {
+	    G4cout<<"musrPhysicsList.cc:  Switching off loss fluctuations for "<<charParticleName<<" in "<<charProcessName<<G4endl;
+	    G4MuIonisation* mmm = (G4MuIonisation*) processTable->FindProcess("muIoni",particleDefinition);
+	    mmm->SetLossFluctuations(false);
+	  }
+	  else {
+	    G4cout<<"musrPhysicsList.cc:  Switching off loss fluctuations for "
+		  <<charParticleName<<" in "<<charProcessName<<" NOT SUPPORTED YET"<<G4endl;
+	  }
+	}
+
+
       }
 
+
       else if (strcmp(tmpString2,"cutForGamma")==0) {
 	sscanf(&line[0],"%*s %*s %*s %g",&tmpCutValue);
 	cutForGamma = tmpCutValue;