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 G4PSFlatSurfaceFlux_h 00031 #define G4PSFlatSurfaceFlux_h 1 00032 00033 #include "G4VPrimitiveScorer.hh" 00034 #include "G4THitsMap.hh" 00035 #include "G4Box.hh" 00036 #include "G4PSDirectionFlag.hh" 00038 // (Description) 00039 // This is a primitive scorer class for scoring Surface Flux. 00040 // Current version assumes only for G4Box shape, and the surface 00041 // is defined at the -Z plane of the box. 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 -Z surface. 00048 // Direction -Z +Z 00049 // 0 IN || OUT ->|<- | fFlux_InOut 00050 // 1 IN ->| | fFlux_In 00051 // 2 OUT |<- | fFlux_Out 00052 // 00053 // Created: 2005-11-14 Tsukasa ASO, Akinori Kimura. 00054 // 00055 // 18-Nov-2005 T.Aso, To use always positive value for anglefactor. 00056 // Bug fix. Area definition. 00057 // 29-Mar-2007 T.Aso, Bug fix for momentum direction at outgoing flux. 00058 // 2010-07-22 Introduce Unit specification. 00059 // 2010-07-22 Add weighted and divideByArea options 00061 00062 class G4PSFlatSurfaceFlux : public G4VPrimitiveScorer 00063 { 00064 public: // with description 00065 G4PSFlatSurfaceFlux(G4String name,G4int direction, G4int depth=0); 00066 G4PSFlatSurfaceFlux(G4String name,G4int direction, 00067 const G4String& unit, G4int depth=0); 00068 virtual ~G4PSFlatSurfaceFlux(); 00069 00070 inline void Weighted(G4bool flg=true) { weighted = flg; } 00071 // Multiply track weight 00072 00073 inline void DivideByArea(G4bool flg=true) { divideByArea = flg; } 00074 // Divided by Area. 00075 00076 protected: // with description 00077 virtual G4bool ProcessHits(G4Step*,G4TouchableHistory*); 00078 G4int IsSelectedSurface(G4Step*,G4Box*); 00079 00080 public: 00081 virtual void Initialize(G4HCofThisEvent*); 00082 virtual void EndOfEvent(G4HCofThisEvent*); 00083 virtual void clear(); 00084 virtual void DrawAll(); 00085 virtual void PrintAll(); 00086 00087 virtual void SetUnit(const G4String& unit); 00088 00089 protected: 00090 virtual void DefineUnitAndCategory(); 00091 00092 private: 00093 G4int HCID; 00094 G4int fDirection; 00095 G4THitsMap<G4double>* EvtMap; 00096 G4bool weighted; 00097 G4bool divideByArea; 00098 }; 00099 00100 #endif 00101