//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
//
// $Id: ReactorRunAction.cc,v 1.9 2006/06/29 17:48:16 gunter Exp $
// GEANT4 tag $Name: geant4-09-01-patch-02 $
// 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

#include "ReactorRunAction.hh"
#include "ReactorRun.hh"
//-- In order to obtain detector information.
#include "G4RunManager.hh"
#include "ReactorConstruction.hh"
#include "G4UnitsTable.hh"
#include "G4THitsMap.hh"

#include "G4Run.hh"



//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

ReactorRunAction::ReactorRunAction()
{
  // - Prepare data member for ReactorRun.
  //   vector represents a list of MultiFunctionalDetector names.
  theSDName.push_back(G4String("Al-SD"));
  theSDName.push_back(G4String("Poly-SD"));
  G4cout<<"0"<< G4endl;
  theSDName.push_back(G4String("GammaSD"));
  
 // theSDName.push_back(G4String("airSD"));
//  theSDName.push_back(G4String("AlAir2SD"));
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

ReactorRunAction::~ReactorRunAction()
{
  theSDName.clear();
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
G4Run* ReactorRunAction::GenerateRun()
{
  G4cout << " In ReactorRunAction.cc GenerateRun " <<G4endl;
  // Generate new RUN object, which is specially
  // dedicated for MultiFunctionalDetector scheme.
  //  Detail description can be found in RE02Run.hh/cc.
  return new ReactorRun(theSDName);
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
void ReactorRunAction::BeginOfRunAction(const G4Run* aRun)
{
  G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void ReactorRunAction::EndOfRunAction(const G4Run* aRun)
{ 
  G4cout << "### Run " << aRun->GetRunID() << " end." << G4endl;

    //- ReactorRun object.
  ReactorRun* gammaRun = (ReactorRun*)aRun;
  //--- Dump all socred quantities involved in RE02Run.
  gammaRun->DumpAllScorer();
  //---
  
 //---------------------------------------------
  // Dump accumulated quantities for this RUN.
  //---------------------------------------------

 G4THitsMap<G4double>* PtotED = gammaRun->GetHitsMap("Poly-SD/totalEDep");
 G4THitsMap<G4double>* totED = gammaRun->GetHitsMap("Al-SD/totalEDep"); 

 
 //------------------------------------------------------------------------------------
 
 G4cout << "1" <<G4endl;
   // G4THitsMap<G4double>* gEcurrentTherm = gammaRun->GetHitsMap("GammaSD/gammaSurfCurrTherm");
  G4THitsMap<G4double>* gEcurrent0 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr000");
 G4THitsMap<G4double>* gEcurrent1 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr001");
 G4THitsMap<G4double>* gEcurrent2 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr002");
 G4THitsMap<G4double>* gEcurrent3 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr003");
 G4THitsMap<G4double>* gEcurrent4 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr004");
 G4THitsMap<G4double>* gEcurrent5 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr005");
 G4THitsMap<G4double>* gEcurrent6 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr006");
 G4THitsMap<G4double>* gEcurrent7 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr007");
 G4THitsMap<G4double>* gEcurrent8 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr008");
 G4THitsMap<G4double>* gEcurrent9 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr009");
 G4THitsMap<G4double>* gEcurrent10 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr010");
 G4THitsMap<G4double>* gEcurrent11 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr011");
 G4THitsMap<G4double>* gEcurrent12 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr012");
 G4THitsMap<G4double>* gEcurrent13 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr013");
 G4THitsMap<G4double>* gEcurrent14 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr014");
 G4THitsMap<G4double>* gEcurrent15 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr015");
 G4THitsMap<G4double>* gEcurrent16 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr016");
 G4THitsMap<G4double>* gEcurrent17 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr017");
 G4THitsMap<G4double>* gEcurrent18 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr018");
 G4THitsMap<G4double>* gEcurrent19 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr019");
 G4THitsMap<G4double>* gEcurrent20 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr020");
 G4THitsMap<G4double>* gEcurrent21 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr021");
 G4THitsMap<G4double>* gEcurrent22 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr022");
 G4THitsMap<G4double>* gEcurrent23 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr023");
 G4THitsMap<G4double>* gEcurrent24 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr024");
 G4THitsMap<G4double>* gEcurrent25 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr025");
 G4THitsMap<G4double>* gEcurrent26 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr026");
 G4THitsMap<G4double>* gEcurrent27 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr027");
 G4THitsMap<G4double>* gEcurrent28 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr028");
 G4THitsMap<G4double>* gEcurrent29 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr029");
 G4THitsMap<G4double>* gEcurrent30 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr030");
 G4THitsMap<G4double>* gEcurrent31 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr031");
 G4THitsMap<G4double>* gEcurrent32 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr032");
 G4THitsMap<G4double>* gEcurrent33 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr033");
 G4THitsMap<G4double>* gEcurrent34 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr034");
 G4THitsMap<G4double>* gEcurrent35 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr035");
 G4THitsMap<G4double>* gEcurrent36 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr036");
 G4THitsMap<G4double>* gEcurrent37 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr037");
 G4THitsMap<G4double>* gEcurrent38 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr038");
 G4THitsMap<G4double>* gEcurrent39 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr039");
 G4THitsMap<G4double>* gEcurrent40 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr040");
 G4THitsMap<G4double>* gEcurrent41 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr041");
 G4THitsMap<G4double>* gEcurrent42 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr042");
 G4THitsMap<G4double>* gEcurrent43 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr043");
 G4THitsMap<G4double>* gEcurrent44 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr044");
 G4THitsMap<G4double>* gEcurrent45 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr045");
 G4THitsMap<G4double>* gEcurrent46 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr046");
 G4THitsMap<G4double>* gEcurrent47 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr047");
 G4THitsMap<G4double>* gEcurrent48 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr048");
 G4THitsMap<G4double>* gEcurrent49 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr049");
 G4THitsMap<G4double>* gEcurrent50 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr050");
 G4THitsMap<G4double>* gEcurrent51 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr051");
 G4THitsMap<G4double>* gEcurrent52 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr052");
 G4THitsMap<G4double>* gEcurrent53 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr053");
 G4THitsMap<G4double>* gEcurrent54 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr054");
 G4THitsMap<G4double>* gEcurrent55 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr055");
 G4THitsMap<G4double>* gEcurrent56 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr056");
 G4THitsMap<G4double>* gEcurrent57 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr057");
 G4THitsMap<G4double>* gEcurrent58 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr058");
 G4THitsMap<G4double>* gEcurrent59 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr059");
 G4THitsMap<G4double>* gEcurrent60 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr060");
 G4THitsMap<G4double>* gEcurrent61 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr061");
 G4THitsMap<G4double>* gEcurrent62 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr062");
 G4THitsMap<G4double>* gEcurrent63 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr063");
 G4THitsMap<G4double>* gEcurrent64 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr064");
 G4THitsMap<G4double>* gEcurrent65 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr065");
 G4THitsMap<G4double>* gEcurrent66 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr066");
 G4THitsMap<G4double>* gEcurrent67 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr067");
 G4THitsMap<G4double>* gEcurrent68 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr068");
 G4THitsMap<G4double>* gEcurrent69 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr069");
 G4THitsMap<G4double>* gEcurrent70 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr070");
 G4THitsMap<G4double>* gEcurrent71 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr071");
 G4THitsMap<G4double>* gEcurrent72 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr072");
 G4THitsMap<G4double>* gEcurrent73 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr073");
 G4THitsMap<G4double>* gEcurrent74 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr074");
 G4THitsMap<G4double>* gEcurrent75 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr075");
 G4THitsMap<G4double>* gEcurrent76 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr076");
 G4THitsMap<G4double>* gEcurrent77 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr077");
 G4THitsMap<G4double>* gEcurrent78 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr078");
 G4THitsMap<G4double>* gEcurrent79 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr079");
 G4THitsMap<G4double>* gEcurrent80 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr080");
 G4THitsMap<G4double>* gEcurrent81 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr081");
 G4THitsMap<G4double>* gEcurrent82 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr082");
 G4THitsMap<G4double>* gEcurrent83 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr083");
 G4THitsMap<G4double>* gEcurrent84 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr084");
 G4THitsMap<G4double>* gEcurrent85 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr085");
 G4THitsMap<G4double>* gEcurrent86 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr086");
 G4THitsMap<G4double>* gEcurrent87 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr087");
 G4THitsMap<G4double>* gEcurrent88 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr088");
 G4THitsMap<G4double>* gEcurrent89 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr089");
 G4THitsMap<G4double>* gEcurrent90 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr090");
 G4THitsMap<G4double>* gEcurrent91 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr091");
 G4THitsMap<G4double>* gEcurrent92 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr092");
 G4THitsMap<G4double>* gEcurrent93 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr093");
 G4THitsMap<G4double>* gEcurrent94 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr094");
 G4THitsMap<G4double>* gEcurrent95 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr095");
 G4THitsMap<G4double>* gEcurrent96 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr096");
 G4THitsMap<G4double>* gEcurrent97 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr097");
 G4THitsMap<G4double>* gEcurrent98 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr098");
 G4THitsMap<G4double>* gEcurrent99 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr099");
 G4THitsMap<G4double>* gEcurrent100 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr100");
 G4THitsMap<G4double>* gEcurrent101 = gammaRun->GetHitsMap("GammaSD/gammaSurfCurr101");
 G4cout << "2" <<G4endl;
 
 //------------------------------------------------------------------------------------
 
 
 G4double nEvt = (G4double)(aRun->GetNumberOfEvent());
  G4cout << "=============================================================" <<G4endl;
	G4cout << " Number of event processed : "<< nEvt << G4endl; 
 G4cout << "=============================================================" <<G4endl;
  	       G4double* AltotED = (*totED)[0,0,0];
		       if ( !AltotED ) AltotED = new G4double(0.0);
  	       G4double* ACtotED = (*PtotED)[0,0,0];
		       if ( !ACtotED ) ACtotED = new G4double(0.0);

			G4double* gEcurr[100];
			G4double nAl = 0.0;
		   gEcurr[1] = (*gEcurrent1)[0,0,0];
  	       gEcurr[2] = (*gEcurrent2)[0,0,0];
  	       gEcurr[3] = (*gEcurrent3)[0,0,0];
  	       gEcurr[4] = (*gEcurrent4)[0,0,0];
	       gEcurr[5] = (*gEcurrent5)[0,0,0];
	       gEcurr[6] = (*gEcurrent6)[0,0,0];
	       gEcurr[7] = (*gEcurrent7)[0,0,0];
	       gEcurr[8] = (*gEcurrent8)[0,0,0];
	       gEcurr[9] = (*gEcurrent9)[0,0,0];
	       gEcurr[10] = (*gEcurrent10)[0,0,0];
	       gEcurr[11] = (*gEcurrent11)[0,0,0];
	       gEcurr[12] = (*gEcurrent12)[0,0,0];
	       gEcurr[13] = (*gEcurrent13)[0,0,0];
	       gEcurr[14] = (*gEcurrent14)[0,0,0];
	       gEcurr[15] = (*gEcurrent15)[0,0,0];
	       gEcurr[16] = (*gEcurrent16)[0,0,0];
	       gEcurr[17] = (*gEcurrent17)[0,0,0];
	       gEcurr[18] = (*gEcurrent18)[0,0,0];
	       gEcurr[19] = (*gEcurrent19)[0,0,0];
	       gEcurr[20] = (*gEcurrent20)[0,0,0];
	       gEcurr[21] = (*gEcurrent21)[0,0,0];
	       gEcurr[22] = (*gEcurrent22)[0,0,0];
	       gEcurr[23] = (*gEcurrent23)[0,0,0];
	       gEcurr[24] = (*gEcurrent24)[0,0,0];
	       gEcurr[25] = (*gEcurrent25)[0,0,0];
	       gEcurr[26] = (*gEcurrent26)[0,0,0];
	       gEcurr[27] = (*gEcurrent27)[0,0,0];
	       gEcurr[28] = (*gEcurrent28)[0,0,0];
	       gEcurr[29] = (*gEcurrent29)[0,0,0];
	       gEcurr[30] = (*gEcurrent30)[0,0,0];
	       gEcurr[31] = (*gEcurrent31)[0,0,0];
	       gEcurr[32] = (*gEcurrent32)[0,0,0];
	       gEcurr[33] = (*gEcurrent33)[0,0,0];
	       gEcurr[34] = (*gEcurrent34)[0,0,0];
	       gEcurr[35] = (*gEcurrent35)[0,0,0];
	       gEcurr[36] = (*gEcurrent36)[0,0,0];
	       gEcurr[37] = (*gEcurrent37)[0,0,0];
	       gEcurr[38] = (*gEcurrent38)[0,0,0];
	       gEcurr[39] = (*gEcurrent39)[0,0,0];
	       gEcurr[40] = (*gEcurrent40)[0,0,0];
	       gEcurr[41] = (*gEcurrent41)[0,0,0];
	       gEcurr[42] = (*gEcurrent42)[0,0,0];
	       gEcurr[43] = (*gEcurrent43)[0,0,0];
	       gEcurr[44] = (*gEcurrent44)[0,0,0];
	       gEcurr[45] = (*gEcurrent45)[0,0,0];
	       gEcurr[46] = (*gEcurrent46)[0,0,0];
	       gEcurr[47] = (*gEcurrent47)[0,0,0];
	       gEcurr[48] = (*gEcurrent48)[0,0,0];
	       gEcurr[49] = (*gEcurrent49)[0,0,0];
	       gEcurr[50] = (*gEcurrent50)[0,0,0];
	       gEcurr[51] = (*gEcurrent51)[0,0,0];
	       gEcurr[52] = (*gEcurrent52)[0,0,0];
	       gEcurr[53] = (*gEcurrent53)[0,0,0];
	       gEcurr[54] = (*gEcurrent54)[0,0,0];
	       gEcurr[55] = (*gEcurrent55)[0,0,0];
	       gEcurr[56] = (*gEcurrent56)[0,0,0];
	       gEcurr[57] = (*gEcurrent57)[0,0,0];
	       gEcurr[58] = (*gEcurrent58)[0,0,0];
	       gEcurr[59] = (*gEcurrent59)[0,0,0];
	       gEcurr[60] = (*gEcurrent60)[0,0,0];
	       gEcurr[61] = (*gEcurrent61)[0,0,0];
	       gEcurr[62] = (*gEcurrent62)[0,0,0];
	       gEcurr[63] = (*gEcurrent63)[0,0,0];
	       gEcurr[64] = (*gEcurrent64)[0,0,0];
	       gEcurr[65] = (*gEcurrent65)[0,0,0];
	       gEcurr[66] = (*gEcurrent66)[0,0,0];
	       gEcurr[67] = (*gEcurrent67)[0,0,0];
	       gEcurr[68] = (*gEcurrent68)[0,0,0];
	       gEcurr[69] = (*gEcurrent69)[0,0,0];
	       gEcurr[70] = (*gEcurrent70)[0,0,0];
	       gEcurr[71] = (*gEcurrent71)[0,0,0];
	       gEcurr[72] = (*gEcurrent72)[0,0,0];
	       gEcurr[73] = (*gEcurrent73)[0,0,0];
	       gEcurr[74] = (*gEcurrent74)[0,0,0];
	       gEcurr[75] = (*gEcurrent75)[0,0,0];
	       gEcurr[76] = (*gEcurrent76)[0,0,0];
	       gEcurr[77] = (*gEcurrent77)[0,0,0];
	       gEcurr[78] = (*gEcurrent78)[0,0,0];
	       gEcurr[79] = (*gEcurrent79)[0,0,0];
	       gEcurr[80] = (*gEcurrent80)[0,0,0];
	       gEcurr[81] = (*gEcurrent81)[0,0,0];
	       gEcurr[82] = (*gEcurrent82)[0,0,0];
	       gEcurr[83] = (*gEcurrent83)[0,0,0];
	       gEcurr[84] = (*gEcurrent84)[0,0,0];
	       gEcurr[85] = (*gEcurrent85)[0,0,0];
	       gEcurr[86] = (*gEcurrent86)[0,0,0];
	       gEcurr[87] = (*gEcurrent87)[0,0,0];
	       gEcurr[88] = (*gEcurrent88)[0,0,0];
	       gEcurr[89] = (*gEcurrent89)[0,0,0];
	       gEcurr[90] = (*gEcurrent90)[0,0,0];
	       gEcurr[91] = (*gEcurrent91)[0,0,0];
	       gEcurr[92] = (*gEcurrent92)[0,0,0];
	       gEcurr[93] = (*gEcurrent93)[0,0,0];
	       gEcurr[94] = (*gEcurrent94)[0,0,0];
	       gEcurr[95] = (*gEcurrent95)[0,0,0];
	       gEcurr[96] = (*gEcurrent96)[0,0,0];
	       gEcurr[97] = (*gEcurrent97)[0,0,0];
	       gEcurr[98] = (*gEcurrent98)[0,0,0];
	       gEcurr[99] = (*gEcurrent99)[0,0,0];
	       gEcurr[100] = (*gEcurrent100)[0,0,0];
			  for ( G4int i = 1; i < 101; i++){

		       if ( !gEcurr[i] ) gEcurr[i] = new G4double(0.0);
			    nAl = nAl + (*gEcurr[i]);
}	   
		G4cout << "  total # through Al " << nAl  <<"  " << G4endl;
	    G4cout << "Tot Eng dep. in 3cm poly moderator" << std::setw(12) << G4BestUnit(*ACtotED,"Energy") << G4endl; 
	    G4cout << "Tot Eng dep. in 1um Al" << std::setw(12) << G4BestUnit(*AltotED,"Energy") << G4endl; 
	    G4cout << "Tot Eng dep. in 1um Al per transmitted e-" << std::setw(12) << G4BestUnit(*AltotED/nAl,"Energy") << G4endl; 
	   
  	       
	  G4cout << "=============================================================" <<G4endl;

	  G4cout << "   " <<G4endl;
	G4cout << "=============================================================" <<G4endl;
	G4cout << "  Gammas heading inward from active core " << G4endl;
//********Electrons************	  
//	G4cout << "  Electrons transmitted through Poly Moderator " << G4endl;

		G4cout << "=============================================================" <<G4endl;
	
		  //***********Al****************	  

//********Gammas************	  
	G4cout << "  Gammas leaving the Poly Moderator and entering the Al Moderator  " << G4endl;
//********Electrons************	  
//	G4cout << "  Electrons transmitted through Al Moderator " << G4endl;
//**************************	  

		G4cout << "=============================================================" <<G4endl;
	//********Gammas************	  
	  G4cout << "  Gammas leaving the micrometer Al tube  and entering the Gamdet in nn-cavity" << G4endl;
//********Electrons************	  
//	  G4cout << "  Electrons through Gamdet in nn-cavity " << G4endl;
//***************************	  

		  for ( G4int i = 1; i < 101; i++){

//		  G4double kmax = std::pow(10.,(G4int)(i/50)*3)*(i % 50+2)*1000/50*eV;

	G4double ebin = std::pow(10.,(G4int)(i/50)*3)*(i % 50+2)*1000/50*eV  ;
       G4cout <<   std::setw(5)  << G4BestUnit(ebin,"Energy") << "  " <<  std::setw(12) << (*gEcurr[i])  << G4endl;
}



		G4cout << "=============================================================" <<G4endl;


  
  G4cout << "============================================="<<G4endl;
}  
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......