PremModel.cxx
//......................................................................
///
/// Set the radius of the outermost layer of the model.
///
/// This usually corresponds to the atmosphere and is useful for
/// computing oscillations with variable neutrino production height.
///
/// @param thick - The thickness of the outer layer in km
///
Set the radius of the outermost layer of the model.

This usually corresponds to the atmosphere and is useful for
computing oscillations with variable neutrino production height.
Set the radius of the outermost layer of the model . This usually corresponds to the atmosphere and is useful for computing oscillations with variable neutrino production height .
void PremModel::SetTopLayerSize(double thick)
{
  if(thick <= 0){
    cout << "Layer thickness should be positive. Do nothing." << endl;
    return;
  }
  int nlayers = fPremLayers.size();
  if(nlayers < 1){
    cout << "PremModel has no layers. Do nothing." << endl;
    return;
  }
  double bottomRadius;
  if(nlayers>1) bottomRadius = fPremLayers[nlayers-2].radius;
  else          bottomRadius = 0;
  double topRadius = bottomRadius + thick;
  fPremLayers[nlayers-1].radius = topRadius;
}
void PremModel :: SetTopLayerSize ( double thick ) { if ( thick <= 0 ) { cout << " " << endl ; return ; } int nlayers = fPremLayers . size ( ) ; if ( nlayers < 1 ) { cout << " " << endl ; return ; } double bottomRadius ; if ( nlayers > 1 ) bottomRadius = fPremLayers [ nlayers - 2 ] . radius ; else bottomRadius = 0 ; double topRadius = bottomRadius + thick ; fPremLayers [ nlayers - 1 ] . radius = topRadius ; }
......................................................................
Set the radius of the outermost layer of the model.
...................................................................... Set the radius of the outermost layer of the model .