Geant4 案例分析（四）：B4

EthanXiao_2002

已于 2024-03-14 18:29:54 修改

阅读量665

点赞数

分类专栏： Geant4 文章标签： c++ 笔记

于 2023-10-18 11:39:06 首次发布

本文链接：https://blog.csdn.net/m0_51809467/article/details/133902115

版权

Geant4 专栏收录该内容

12 篇文章

订阅专栏

本文详细介绍了如何通过B4案例使用G4AnalysisManager生成并输出模拟数据，包括创建直方图（Histograms）、NTuple以及在SteppingAction中进行特定Volume的数据收集。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

B4 案例通过不同方式实现了模拟数据的输出，以B4a为例，通过G4AnalysisManager输出root文件

AnalysisManager的使用

//  RunAction
  auto analysisManager = G4AnalysisManager::Instance();
  // Creating histograms
  analysisManager->CreateH1("Eabs","Edep in absorber", 100, 0., 800*MeV);
  analysisManager->CreateH1("Egap","Edep in gap", 100, 0., 100*MeV);
  analysisManager->CreateH1("Labs","trackL in absorber", 100, 0., 1*m);
  analysisManager->CreateH1("Lgap","trackL in gap", 100, 0., 50*cm);

  // Creating ntuple
  //
  analysisManager->CreateNtuple("B4", "Edep and TrackL");
  analysisManager->CreateNtupleDColumn("Eabs");
  analysisManager->CreateNtupleDColumn("Egap");
  analysisManager->CreateNtupleDColumn("Labs");
  analysisManager->CreateNtupleDColumn("Lgap");
  analysisManager->FinishNtuple();
  
  
//  RunAction::BeginOfRunAction
  // Open an output file
  //
  G4String fileName = "B4.root";
  // Other supported output types:
  // G4String fileName = "B4.csv";
  // G4String fileName = "B4.hdf5";
  // G4String fileName = "B4.xml";
  analysisManager->OpenFile(fileName);
  G4cout << "Using " << analysisManager->GetType() << G4endl;
  
  
// EventAction::EndOfEventAction
  // fill histograms
  analysisManager->FillH1(0, fEnergyAbs);
  analysisManager->FillH1(1, fEnergyGap);
  analysisManager->FillH1(2, fTrackLAbs);
  analysisManager->FillH1(3, fTrackLGap);

  // fill ntuple
  analysisManager->FillNtupleDColumn(0, fEnergyAbs);
  analysisManager->FillNtupleDColumn(1, fEnergyGap);
  analysisManager->FillNtupleDColumn(2, fTrackLAbs);
  analysisManager->FillNtupleDColumn(3, fTrackLGap);
  analysisManager->AddNtupleRow();
  
  
// RunAction::EndOfRunAction

  auto analysisManager = G4AnalysisManager::Instance();
  analysisManager->Write();
  analysisManager->CloseFile();

补充：可视化

设置可视化，worldLV为隐形，并设置calorLV的颜色

  worldLV->SetVisAttributes (G4VisAttributes::GetInvisible());

  auto simpleBoxVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,1.0));
  simpleBoxVisAtt->SetVisibility(true);
  calorLV->SetVisAttributes(simpleBoxVisAtt);

补充：Stepping判断volume

获取当前粒子的volume与DetectorConstruction中预设的Volume，通过判断是否为同一对象来实现特定Volume中数据的输出。

// SteppingAction::SteppingAction
  fDetConstruction = detConstruction;
  fEventAction = eventAction;
  
// SteppingAction::UserSteppingAction
  // get volume of the current step
  auto volume = step->GetPreStepPoint()->GetTouchableHandle()->GetVolume();
  if ( volume == fDetConstruction->GetAbsorberPV() ) {
    fEventAction->AddAbs(edep,stepLength);
  }