G4PSSphereSurfaceCurrent3D.cc

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 //
00027 // $Id$
00028 //
00029 // G4PSSphereSurfaceCurrent3D
00030 #include "G4PSSphereSurfaceCurrent3D.hh"
00031 
00033 // (Description)
00034 //   This is a primitive scorer class for scoring Surface Current.
00035 //  Current version assumes only for G4Sphere shape, and the surface
00036 //  is defined at the inside of the sphere.
00037 //   The current is given in the unit of area. 
00038 //    e.g.  (Number of tracks)/mm2.
00039 //
00040 // Surface is defined  at the inside of sphere.
00041 // Direction                  -Rmin   +Rmax
00042 //   0  IN || OUT            ->|<-     |      fCurrent_InOut
00043 //   1  IN                   ->|       |      fCurrent_In
00044 //   2  OUT                    |<-     |      fCurrent_Out
00045 //
00046 // Created: 2007-08-14  Tsukasa ASO
00047 // 2010-07-22   Introduce Unit specification.
00048 // 
00050 
00051 G4PSSphereSurfaceCurrent3D::G4PSSphereSurfaceCurrent3D(G4String name, G4int direction,
00052                                G4int ni, G4int nj, G4int nk,
00053                                G4int depi, G4int depj, G4int depk)
00054     :G4PSSphereSurfaceCurrent(name,direction),
00055      fDepthi(depi),fDepthj(depj),fDepthk(depk)
00056 {
00057   fNi=ni;
00058   fNj=nj;
00059   fNk=nk;
00060 }
00061 
00062 G4PSSphereSurfaceCurrent3D::G4PSSphereSurfaceCurrent3D(G4String name, 
00063                                                        G4int direction,
00064                                                        const G4String& unit,
00065                                G4int ni, G4int nj, G4int nk,
00066                                G4int depi, G4int depj, G4int depk)
00067     :G4PSSphereSurfaceCurrent(name,direction),
00068      fDepthi(depi),fDepthj(depj),fDepthk(depk)
00069 {
00070   fNi=ni;
00071   fNj=nj;
00072   fNk=nk;
00073   SetUnit(unit);
00074 }
00075 
00076 G4PSSphereSurfaceCurrent3D::~G4PSSphereSurfaceCurrent3D()
00077 {;}
00078 
00079 G4int G4PSSphereSurfaceCurrent3D::GetIndex(G4Step* aStep)
00080 {
00081   const G4VTouchable* touchable = aStep->GetPreStepPoint()->GetTouchable();
00082   G4int i = touchable->GetReplicaNumber(fDepthi);
00083   G4int j = touchable->GetReplicaNumber(fDepthj);
00084   G4int k = touchable->GetReplicaNumber(fDepthk);
00085   
00086   return i*fNj*fNk+j*fNk+k;
00087 }

Generated on Mon May 27 17:49:29 2013 for Geant4 by  doxygen 1.4.7