mt01runaction.cc~

利用linux环境下蒙特卡罗计算工具geant4实现伽玛射线与探测器乘积能量的模拟过程

//
// ********************************************************************
// * DISCLAIMER                                                       *
// *                                                                  *
// * The following disclaimer summarizes all the specific disclaimers *
// * of contributors to this software. The specific disclaimers,which *
// * govern, are listed with their locations in:                      *
// *   http://cern.ch/geant4/license                                  *
// *                                                                  *
// * 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.                                                             *
// *                                                                  *
// * This  code  implementation is the  intellectual property  of the *
// * GEANT4 collaboration.                                            *
// * By copying,  distributing  or modifying the Program (or any work *
// * based  on  the Program)  you indicate  your  acceptance of  this *
// * statement, and all its terms.                                    *
// ********************************************************************
//
//
// $Id: ExN02RunAction.cc,v 1.8 2004/03/08 15:14:49 maire Exp $
// GEANT4 tag $Name: geant4-07-00-patch-01 $

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

#include "mt01RunAction.hh"


#include "G4Run.hh"

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

mt01RunAction::mt01RunAction()
{}

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

mt01RunAction::~mt01RunAction()
{}

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

void mt01RunAction::BeginOfRunAction(const G4Run* aRun)
{

  G4cout <<"mataoflag"<<G4endl;


}

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

void mt01RunAction::EndOfRunAction(const G4Run*)
{
  
  G4int i=0;
   for(i=0;i<10001;i++)
    {
      switch(i)
{
 case 112:
   // G4cout<<i<<G4endl;
   G4cout<<((energycount[i])/96385.)<<G4endl;
      break;
 case 200:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 300:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 400:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 500:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 662:
   //G4cout<<i<<G4endl;
      G4cout<<((energycount[i])/78095.)<<G4endl;
      break;
  case 700:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 800:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 900:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 1130:
   //G4cout<<i<<G4endl;
      G4cout<<((energycount[i])/59550.)<<G4endl;
      break;
 case 1460:
   //G4cout<<i<<G4endl;
      G4cout<<((energycount[i])/51365.)<<G4endl;
      break; 
case 2615:
  //G4cout<<i<<G4endl;
      G4cout<<((energycount[i])/37237.)<<G4endl;
      break;
 case 4000:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break; 
case 5000:
  //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break; 
case 6000:
  //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 7000:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 8000:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break; 
case 9000:
  //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break;
 case 10000:
   //G4cout<<i<<G4endl;
      G4cout<<(energycount[i])<<G4endl;
      break; 

}
    }   
  G4cout<<"endflag"<<G4endl;
 }

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