Changeset 1340 for trunk/source/digits_hits/scorer/src/G4PSSphereSurfaceCurrent.cc

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/source/digits_hits/scorer/src/G4PSSphereSurfaceCurrent.cc (modified) (5 diffs)

trunk/source/digits_hits/scorer/src/G4PSSphereSurfaceCurrent.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4PSSphereSurfaceCurrent.cc,v 1.4 2009/11/14 00:01:13 asaim Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: G4PSSphereSurfaceCurrent.cc,v 1.6 2010/07/23 04:35:38 taso Exp $
+// GEANT4 tag $Name:  $
 //
 // G4PSSphereSurfaceCurrent
@@ -49 +49 @@
 //
 // Created: 2005-11-14  Tsukasa ASO, Akinori Kimura.
+// 2010-07-22   Introduce Unit specification.
 // 
 ///////////////////////////////////////////////////////////////////////////////
@@ -56 +57 @@
     :G4VPrimitiveScorer(name,depth),HCID(-1),fDirection(direction),
      weighted(true),divideByArea(true)
-{;}
+{
+    DefineUnitAndCategory();
+    SetUnit("percm2");
+}
+
+G4PSSphereSurfaceCurrent::G4PSSphereSurfaceCurrent(G4String name, 
+                                                   G4int direction, 
+                                                   const G4String& unit,
+                                                   G4int depth)
+    :G4VPrimitiveScorer(name,depth),HCID(-1),fDirection(direction),
+     weighted(true),divideByArea(true)
+{
+    DefineUnitAndCategory();
+    SetUnit(unit);
+}
 
 G4PSSphereSurfaceCurrent::~G4PSSphereSurfaceCurrent()
@@ -79 +94 @@
   }
 
-//  if( solid->GetEntityType() != "G4Sphere" ){
-//    G4Exception("G4PSSphereSurfaceCurrentScorer. - Solid type is not supported.");
-//    return FALSE;
-//  }
   G4Sphere* sphereSolid = (G4Sphere*)(solid);
 
@@ -177 +188 @@
   for(; itr != EvtMap->GetMap()->end(); itr++) {
     G4cout << "  copy no.: " << itr->first  << "  current  : " ;
-    if ( divideByArea ) G4cout << *(itr->second)*cm*cm << " [/cm2]" ;
-    else G4cout << *(itr->second)*cm*cm << " [track]" ;
+    if ( divideByArea ) {
+        G4cout << *(itr->second)/GetUnitValue() 
+               << " [" <<GetUnit()<<"]";
+    }else {
+        G4cout << *(itr->second) << " [tracks]" ;
+    }
     G4cout  << G4endl;
   }
 }
 
+
+void G4PSSphereSurfaceCurrent::SetUnit(const G4String& unit)
+{
+    if ( divideByArea ) {
+        CheckAndSetUnit(unit,"Per Unit Surface");
+    } else {
+        if (unit == "" ){
+            unitName = unit;
+            unitValue = 1.0;
+        }else{
+            G4String msg = "Invalid unit ["+unit+"] (Current  unit is [" +GetUnit()+"] )";
+            G4Exception(GetName(),"DetScorer0000",JustWarning,msg);
+        }
+    }
+}
+
+void G4PSSphereSurfaceCurrent::DefineUnitAndCategory(){
+   // Per Unit Surface
+   new G4UnitDefinition("percentimeter2","percm2","Per Unit Surface",(1./cm2));
+   new G4UnitDefinition("permillimeter2","permm2","Per Unit Surface",(1./mm2));
+   new G4UnitDefinition("permeter2","perm2","Per Unit Surface",(1./m2));
+}