发布时间 : 星期四 文章Geant4 程序编写中的常用代码更新完毕开始阅读caec6b4067ec102de3bd8920
Geant4 程序编写中的常用代码 2010-06-01 20:51:44| 分类： 学习与工作 | 标签：geant4 |字号 订阅 系个人原创，未经同意请勿转载。cnscott的博客 1. 各种常用源粒子的产生 cnscott的博客 G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable(); ? 光子：
G4ParticleDefinition* particle
= particleTable->FindParticle(\? 中子：
G4ParticleDefinition* particle
= particleTable->FindParticle(\? 电子：
G4ParticleDefinition* particle
= particleTable->FindParticle(\? Alpha：
G4ParticleDefinition* particle
= particleTable->FindParticle(\? 质子：
G4ParticleDefinition* particle
= particleTable->FindParticle(\? 离子：
G4ParticleDefinition* particle
= particleTable->GetIon(6, 12, 0.0); //原子序数，质量数和激发能量(基态12C）
源粒子的能量、位置、方向、电荷等相关信息要在ParticleGun中设置。
2. 获取日期，记录程序开始运行的时刻 cnscott的博客 #include \ #include \ struct tm *ptr; time_t lt; lt =time(NULL); ptr=gmtime(<);
printf(asctime(ptr)); //格林威治时间 printf(ctime(<)); //本地时间
3. 记录程序运行所用时间 cnscott的博客 #include \ G4Timer myTimer; myTimer.Start(); myTimer.Stop();
G4cout << myTimer << G4endl;
4. 保存文件 cnscott的博客 #include
fstream dataFile(FileName, ios::out); //文件输出流 G4String mydata=”Save my data.”;
dataFile << mydata << G4endl; //数据可以是整数、浮点数或是字符串等。
5. 检查几何体是否重叠 cnscott的博客
G4PVPlacement和G4PVParameterised提供了成员函数: G4bool CheckOverlaps(G4int res=1000, G4double tol=0., G4bool verbose=true)
//当返回值为true时，表明几何存在重叠。
6. 获取相关类的指针 cnscott的博客
G4RunManager* runManager= G4RunManager::GetRunManager(); (MyRunAction *)runManager->GetUserRunAction(); (MyEventAction *)runManager->GetUserEventAction(); (MyStackingAction *)runManager->GetUserStackingAction(); (MySteppingAction *)runManager->GetUserSteppingAction(); (MyTrackingAction *)runManager->GetUserTrackingAction(); G4Run* run = G4RunManager::GetRunManager()->GetCurrentRun(); G4Event* event = G4RunManager::GetRunManager()->GetCurrentEvent(); //版本号
G4String version = G4RunManager::GetRunManager()->GetVersionString()
7. 随机数操作 cnscott的博客
Geant4的随机数产生主要是使用CLHEP函数库。
http://proj-clhep.web.cern.ch/proj-clhep/manual/UserGuide/Random/Random.html
CLHEP::HepRandom::setTheEngine(HepRandomEngine *theNewEngine); //设置随机数引擎，默认引擎为CLHEP::HepJamesRandom。其他可选的引擎包括： CLHEP::DRand48Engine, CLHEP::DualRand, CLHEP::Hurd160Engine, CLHEP::Hurd288Engine, CLHEP::MTwistEngine,
CLHEP::NonRandomEngine, CLHEP::RandEngine, CLHEP::RanecuEngine,
CLHEP::Ranlux64Engine, CLHEP::RanluxEngine, CLHEP::RanshiEngine, and CLHEP::TripleRand.
CLHEP::HepRandom::setTheSeed(long seed, int lux=3); //为默认引擎设置随机数种子，其他引擎可能使用不同函数
G4UniformRand();
//产生均匀分布随机数，等同于CLHEP::HepRandom::getTheEngine()->flat();
G4RandGauss::shoot(double mean, double stdDev);
//产生高斯分布随机数，等同于CLHEP::RandGaussQ::shoot(double mean, double stdDev);
8. 从G4Step中读取有用的信息 cnscott的博客
void MySteppingAction::UserSteppingAction(const G4Step* stp) {
G4Track* theTrack = stp->GetTrack(); //获取G4Track指针
theTrack->GetVolume()->GetName(); //获取当前几何体（PhysicalVolume）名称
G4StepPoint* preStepPoint = stp->GetPreStepPoint();
G4StepPoint* postStepPoint = stp->GetPostStepPoint(); //获取两个StepPoint的指针
if (preStepPoint->GetStepStatus() == fGeomBoundary) //判断是否位于几何体边界
std::cout << \
if (postStepPoint->GetStepStatus() == fGeomBoundary) //判断是否位于几何体边界
std::cout << \
if (theTrack->GetTrackID() == 1) //判断是否是初级粒子（初级粒子TrackID=1） std::cout << \；