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 00030 #ifndef G4PSCylinderSurfaceFlux_h 00031 #define G4PSCylinderSurfaceFlux_h 1 00032 00033 #include "G4VPrimitiveScorer.hh" 00034 #include "G4THitsMap.hh" 00035 #include "G4Tubs.hh" 00036 #include "G4PSDirectionFlag.hh" 00038 // (Description) 00039 // This is a primitive scorer class for scoring Surface Flux. 00040 // Current version assumes only for G4Tubs shape, and the surface 00041 // is defined at the inner plane of the tube. 00042 // The surface flux is given in the unit of area. 00043 // e.g. sum of 1/cos(T)/mm2, where T is a incident angle of the 00044 // track on the surface. 00045 // 00046 // 00047 // Surface is defined at the inner surface of the tubs. 00048 // Direction R R+dR 00049 // 0 IN || OUT ->|<- | fFlux_InOut 00050 // 1 IN ->| | fFlux_In 00051 // 2 OUT |<- | fFlux_Out 00052 // 00053 // Created: 2007-03-29 Tsukasa ASO 00054 // 2010-07-22 Introduce Unit specification. 00055 // 2010-07-22 Add weighted and divideByArea options 00056 // 00058 00059 class G4PSCylinderSurfaceFlux : public G4VPrimitiveScorer 00060 { 00061 public: // with description 00062 G4PSCylinderSurfaceFlux(G4String name,G4int direction, G4int depth=0); 00063 G4PSCylinderSurfaceFlux(G4String name,G4int direction, 00064 const G4String& unit, G4int depth=0); 00065 virtual ~G4PSCylinderSurfaceFlux(); 00066 00067 inline void Weighted(G4bool flg=true) { weighted = flg; } 00068 // Multiply track weight 00069 00070 inline void DivideByArea(G4bool flg=true) { divideByArea = flg; } 00071 // Divided by Area. 00072 00073 00074 protected: // with description 00075 virtual G4bool ProcessHits(G4Step*,G4TouchableHistory*); 00076 G4int IsSelectedSurface(G4Step*,G4Tubs*); 00077 00078 public: 00079 virtual void Initialize(G4HCofThisEvent*); 00080 virtual void EndOfEvent(G4HCofThisEvent*); 00081 virtual void clear(); 00082 virtual void DrawAll(); 00083 virtual void PrintAll(); 00084 00085 virtual void SetUnit(const G4String& unit); 00086 00087 protected: 00088 virtual void DefineUnitAndCategory(); 00089 00090 private: 00091 G4int HCID; 00092 G4int fDirection; 00093 G4THitsMap<G4double>* EvtMap; 00094 G4bool weighted; 00095 G4bool divideByArea; 00096 }; 00097 00098 #endif 00099