2 // ********************************************************************
 
    3 // * License and Disclaimer                                           *
 
    5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
 
    6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
 
    7 // * conditions of the Geant4 Software License,  included in the file *
 
    8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
 
    9 // * include a list of copyright holders.                             *
 
   11 // * Neither the authors of this software system, nor their employing *
 
   12 // * institutes,nor the agencies providing financial support for this *
 
   13 // * work  make  any representation or  warranty, express or implied, *
 
   14 // * regarding  this  software system or assume any liability for its *
 
   15 // * use.  Please see the license in the file  LICENSE  and URL above *
 
   16 // * for the full disclaimer and the limitation of liability.         *
 
   18 // * This  code  implementation is the result of  the  scientific and *
 
   19 // * technical work of the GEANT4 collaboration.                      *
 
   20 // * By using,  copying,  modifying or  distributing the software (or *
 
   21 // * any work based  on the software)  you  agree  to acknowledge its *
 
   22 // * use  in  resulting  scientific  publications,  and indicate your *
 
   23 // * acceptance of all terms of the Geant4 Software license.          *
 
   24 // ********************************************************************
 
   27 // $Id: G4EllipticalCone.icc 83572 2014-09-01 15:23:27Z gcosmo $
 
   30 // --------------------------------------------------------------------
 
   31 // GEANT 4 inline definitions file
 
   33 // G4EllipticalCone.icc
 
   35 // Implementation of inline methods of G4EllipticalCone
 
   36 // --------------------------------------------------------------------
 
   37 // Author: Dionysios Anninos
 
   38 // --------------------------------------------------------------------
 
   41 G4double G4EllipticalCone::GetSemiAxisMax () const
 
   43   return ySemiAxis > xSemiAxis ? ySemiAxis : xSemiAxis;
 
   47 G4double G4EllipticalCone::GetSemiAxisX () const
 
   53 G4double G4EllipticalCone::GetSemiAxisY () const
 
   59 G4double G4EllipticalCone::GetZMax() const
 
   65 G4double G4EllipticalCone::GetZTopCut() const
 
   71 void G4EllipticalCone::SetSemiAxis (G4double newxSemiAxis,
 
   72                                     G4double newySemiAxis,
 
   75   xSemiAxis = newxSemiAxis; 
 
   76   ySemiAxis = newySemiAxis; 
 
   78   semiAxisMax = xSemiAxis > ySemiAxis ? xSemiAxis : ySemiAxis;
 
   79   if (zTopCut > +zheight) { zTopCut = +zheight; }
 
   80   fRebuildPolyhedron = true;
 
   84 void G4EllipticalCone::SetZCut (G4double newzTopCut)
 
   86   if (newzTopCut > +zheight)
 
   87     { zTopCut = +zheight; }
 
   89     { zTopCut = newzTopCut; }
 
   90   fRebuildPolyhedron = true;
 
   94 G4double G4EllipticalCone::GetCubicVolume()
 
   96   if(fCubicVolume != 0 ) {;}
 
   99     if (zTopCut > +zheight )
 
  101       fCubicVolume = (8./3.)*CLHEP::pi*xSemiAxis*
 
  102                      ySemiAxis*zheight*zheight*zheight;
 
  106       fCubicVolume = CLHEP::pi*xSemiAxis*ySemiAxis*
 
  107                     (2./3.*std::pow(zTopCut,3.)+2.*sqr(zheight)*zTopCut);
 
  114 G4double G4EllipticalCone::GetSurfaceArea()
 
  116   if(fSurfaceArea != 0.) {;}
 
  117   else   { fSurfaceArea = G4VSolid::GetSurfaceArea(); }