ChinaNet

标题: 关于粒子收集的一些问题 [打印本页]

作者: 鸿飞冥冥 时间: 2015-4-8 09:55
标题: 关于粒子收集的一些问题
大家好：
初学编程对次级离子的收集不甚了解，还请高手指点。比如说，在许多的出射粒子中我只想要电子的信息，包括能谱，角分布。我该怎么实现呢，具体的该看哪些例子呢？
Geant4/geant4.9.0/examples/Novice/ExN07
用Particle filter这个method
比较简单的可以在steppingAction里面这么做
G4Track * theTrack = theStep->GetTrack();

// 判断是不是gamma光子
G4ParticleDefinition * particleType = theTrack->GetDefinition();
if(particleType==G4Gamma::GammaDefinition())
……

如果要判断电子改成 G4Electron::ElectronDefinition()就行了

Re:关于粒子收集的一些问题

给你一段我的代码作参考吧

void testSteppingAction::UserSteppingAction(const G4Step* theStep)
//这里是从G4自带的G4UserSteppingAction派生了一个类，
//这个函数就是替换G4UserSteppingAction中的UserSteppingAction(const G4Step* theStep)函数
{
G4Track * theTrack = theStep->GetTrack();

  // 判断是不是gamma光子
  G4ParticleDefinition * particleType = theTrack->GetDefinition();
  if(particleType==G4Gamma::GammaDefinition())
  {
  G4VPhysicalVolume * thePrePV = theTrack->GetVolume();
  G4VPhysicalVolume * theNextPV = theTrack->GetNextVolume();
  G4String thePrePVname = thePrePV->GetName();

  //判断光子是否进入统计位置
  if(theNextPV!=NULL)
  {
G4String theNextPVname = theNextPV->GetName();
//判断光子是否进入探测面
if((theNextPVname(0,10)=="Dec_Sphere")&&
   (thePrePVname(0,10)!="Dec_Sphere"))  //Dec_Sphere是我设定的探测面名字
{
   //G4cout<<"enter the Detector volume"<<G4endl; //for Debug only！
G4double KineEnergy=theTrack->GetKineticEnergy();  //获取当前step的光子动能
G4ThreeVector GammaPosition=theTrack->GetPosition();  //获取当前step位置
G4double KineTst=KineEnergy/MeV;  //将能量单位转换成MeV
//-============================================================================================//

//统计能谱
G4double gammaEnergy = KineEnergy;
EnergyID=int(gammaEnergy/keV); //将当前光子归入某个能量道
//G4cout<<"Energy is"<< EnergyID <<G4endl; //for Debug only！
theRun->CountID[EnergyID]+=1;  //该能量道的光子数+1，
                                                      //这里调用了RunAction中的数组
                                                      //CountID[MaxEnergyID]进行数据传递
//G4cout<<"Counts of"<< EnergyID <<"is"<< theRun->CountID[EnergyID] <<G4endl;
                                                                                                      //for Debug only！
//-============================================================================================//
//统计角分布
KineTst=KineTst*CalculateAirCoef(KineTst);
//上面的CalculateAirCoef()是自己定义的一个计算不同能量X射线空气质量能量吸收系数的函数
//上面的式子就是将能量转换成空气中X射线能量沉积
G4double zPosition=GammaPosition.z();  //获取位置的z方向（主束方向）坐标
zPosition=zPosition/mm; //坐标换成mm单位
G4double cosPhi=zPosition/1000; //计算极角余弦，
                                                      //这里1000是我设定的球状探测面的内径，外径是1000.1，
                                                      //单位mm
G4int Phi=int(acos(cosPhi)/deg); //计算出极角
//G4cout<<"Phi is "<< Phi <<G4endl; //for Debug only！

theRun->NumID[Phi]+=1;    //该极角内光子统计数+1
theRun->EDeposit[Phi]+=KineTst; //该极角内光子能量沉积+当前track到达探测面的光子能量沉积
G4double Area=cos(Phi*deg)-cos((Phi+1)*deg);  //计算每个极角占有的归一化面积
theRun->Dose[Phi]=theRun->EDeposit[Phi]/Area;  //将能量沉积按面积归一化，即转化成相对剂量
//G4cout <<"\t"<<Phi<<"\t"<<NumID[Phi]<<"\t\t"<<EDeposit[Phi]<<"\t"<<Dose[Phi]<<endl;
                                                                                                            //for Debug only！
  }
  }
}
}

作者: 鸿飞冥冥 时间: 2015-4-8 09:55
你好,我看了N07的粒子了,可是我要的信息,找不到 Particle filte
上面是你的问题
下面是N07中particle filter的段代码
可以只收集光子，电子等等

G4SDParticleFilter* gammaFilter = new G4SDParticleFilter

(filterName="gammaFilter",particleName="gamma");
  G4SDParticleFilter* electronFilter = new G4SDParticleFilter

(filterName="electronFilter",particleName="e-");
  G4SDParticleFilter* positronFilter = new G4SDParticleFilter

(filterName="positronFilter",particleName="e+");
  G4SDParticleFilter* epFilter = new G4SDParticleFilter

(filterName="epFilter");
  epFilter->add(particleName="e-");
  epFilter->add(particleName="e+");
G4VPrimitiveScorer* primitive;
primitive = new G4PSEnergyDeposit("eDep",j);
det->RegisterPrimitive(primitive);
primitive = new G4PSNofSecondary("nGamma",j);
primitive->SetFilter(gammaFilter);
det->RegisterPrimitive(primitive);
primitive = new G4PSNofSecondary("nElectron",j);
primitive->SetFilter(electronFilter);
det->RegisterPrimitive(primitive);
primitive = new G4PSNofSecondary("nPositron",j);
primitive->SetFilter(positronFilter);
det->RegisterPrimitive(primitive);
primitive = new G4PSMinKinEAtGeneration("minEkinGamma",j);
primitive->SetFilter(gammaFilter);
det->RegisterPrimitive(primitive);
primitive = new G4PSMinKinEAtGeneration("minEkinElectron",j);
primitive->SetFilter(electronFilter);
det->RegisterPrimitive(primitive);
primitive = new G4PSMinKinEAtGeneration("minEkinPositron",j);
primitive->SetFilter(positronFilter);
det->RegisterPrimitive(primitive);
primitive = new G4PSTrackLength("trackLength",j);
primitive->SetFilter(epFilter);
det->RegisterPrimitive(primitive);
primitive = new G4PSNofStep("nStep",j);
primitive->SetFilter(epFilter);
det->RegisterPrimitive(primitive);

作者: 鸿飞冥冥 时间: 2015-4-8 09:55
theRun->CountID[EnergyID]+=1;    //该能量道的光子数+1，
                                                      //这里调用了RunAction中的数组
                                                      //CountID[MaxEnergyID]进行数据传递

theRun没有定义阿，要怎么调用呢，能不能把调用讲讲～～
具体就是说SteppingAction和RunAction的头文件要怎么写呢

欢迎光临 ChinaNet (http://www.nuclear.net.cn/)