Changeset 1340 for trunk/source/digits_hits/utils/src/G4VScoreWriter.cc

Timestamp:

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

Author:

garnier

Message:

update ti head

File:

• trunk/source/digits_hits/utils/src/G4VScoreWriter.cc (modified) (5 diffs)

trunk/source/digits_hits/utils/src/G4VScoreWriter.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4VScoreWriter.cc,v 1.5 2008/03/04 23:19:09 taso Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: G4VScoreWriter.cc,v 1.10 2010/07/27 01:44:54 akimura Exp $
+// GEANT4 tag $Name:  $
 //
 
@@ -87 +87 @@
     return;
   }
+
+
   std::map<G4int, G4double*> * score = msMapItr->second->GetMap();
-  ofile << "# primitive scorer name: " << msMapItr->first << G4endl;
+  ofile << "# primitive scorer name: " << msMapItr->first << std::endl;
+
+
+  G4double unitValue = fScoringMesh->GetPSUnitValue(psName);
+  G4String unit = fScoringMesh->GetPSUnit(psName);
+  G4String divisionAxisNames[3];
+  fScoringMesh->GetDivisionAxisNames(divisionAxisNames);
+  // index order
+  ofile << "# i" << divisionAxisNames[0]
+        << ", i" << divisionAxisNames[1]
+        << ", i" << divisionAxisNames[2];
+  // unit of scored value
+  ofile << ", value ";
+  if(unit.size() > 0) ofile << "[" << unit << "]";
+  ofile << G4endl;
 
   // "sequence" option: write header info 
@@ -111 +127 @@
           ofile << 0.;
         } else {
-          ofile << *(value->second);
+          ofile << *(value->second)/unitValue;
         }
 
@@ -160 +176 @@
   std::map<G4int, G4double*> * score;
   for(; msMapItr != fSMap.end(); msMapItr++) {
+
+    G4String psname = msMapItr->first;
+
     score = msMapItr->second->GetMap();
-    ofile << "# primitive scorer name: " << msMapItr->first << G4endl;
+    ofile << "# primitive scorer name: " << msMapItr->first << std::endl;
+
+    G4double unitValue = fScoringMesh->GetPSUnitValue(psname);
+    G4String unit = fScoringMesh->GetPSUnit(psname);
+    G4String divisionAxisNames[3];
+    fScoringMesh->GetDivisionAxisNames(divisionAxisNames);
+    // index order
+    ofile << "# i" << divisionAxisNames[0]
+          << ", i" << divisionAxisNames[1]
+          << ", i" << divisionAxisNames[2];
+    // unit of scored value
+    ofile << ", value ";
+    if(unit.size() > 0) ofile << "[" << unit << "]";
+    ofile << G4endl;
+
 
     // "sequence" option: write header info 
@@ -184 +217 @@
             ofile << 0.;
           } else {
-            ofile << *(value->second);
+            ofile << *(value->second)/unitValue;
           }