#include "Geant/Solid.h" #include "Vtk/HEP/Geant/vtkGeantEvent.h" #include "HEP/Geant/Scene.h" #include "HEP/Geant/EmitterPrimary.hh" #include "Vtk/Math/vtkContainerBox.h" #include "Vtk/Math/vtkDense.h" #include "Math/Units.h" #include "Vtk/uLibVtkViewer.h" #include "Vtk/HEP/Geant/vtkGeantEvent.h" #include "Vtk/HEP/Geant/vtkEmitterPrimary.h" #include "Vtk/Math/vtkContainerBox.h" #include "Vtk/HEP/Detectors/vtkDetectorChamber.h" #include "Vtk/HEP/Detectors/vtkDetectorChamber.h" #include #include #include #include #include #include #include #include #include #include using namespace uLib; struct AppState { Geant::Scene* scene; Vtk::Viewer* viewer; Vector results; }; void KeyPressCallbackFunction(vtkObject* caller, long unsigned int eventId, void* clientData, void* callData) { auto* interactor = static_cast(caller); auto* state = static_cast(clientData); std::string key = interactor->GetKeySym(); if (key == "Return") { std::cout << "--> Firing muon from current emitter position..." << std::endl; // Run one event state->scene->RunSimulation(1, state->results); if (!state->results.empty()) { // Get the last event Geant::GeantEvent* lastEvent = &state->results.back(); std::cout << " Collected event " << lastEvent->GetEventID() << " with " << lastEvent->Path().size() << " steps." << std::endl; // Wrap it for VTK Vtk::GeantEvent* vtkEvent = new Vtk::GeantEvent(lastEvent); state->viewer->AddPuppet(*vtkEvent); // Re-render state->viewer->GetRenderer()->Render(); state->viewer->GetRenderWindow()->Render(); } else { std::cout << " No event collected." << std::endl; } } } int main(int argc, char** argv) { // 1. Setup Geant4 Scene Geant::Scene scene; scene.ConstructWorldBox(Vector3f(30_m, 30_m, 30_m), "G4_AIR"); ContainerBox iron_box; iron_box.Scale(Vector3f(10_m, 10_m, 10_m)); iron_box.SetPosition(Vector3f(-5_m, -5_m, -5_m)); Geant::BoxSolid* iron_cube = new Geant::BoxSolid("IronCube", &iron_box); iron_cube->SetNistMaterial("G4_Fe"); iron_cube->Update(); scene.AddSolid(iron_cube); // Top Detector Chamber DetectorChamber* top_chamber_box = new DetectorChamber(); top_chamber_box->Scale(Vector3f(10_m, 10_m, 40_cm)); top_chamber_box->SetPosition(Vector3f(-5_m, -5_m, 7_m)); Geant::BoxSolid* top_chamber = new Geant::BoxSolid("TopChamber", top_chamber_box); top_chamber->SetNistMaterial("G4_AIR"); top_chamber->Update(); scene.AddSolid(top_chamber); // Bottom Detector Chamber DetectorChamber* bottom_chamber_box = new DetectorChamber(); bottom_chamber_box->Scale(Vector3f(10_m, 10_m, 40_cm)); bottom_chamber_box->SetPosition(Vector3f(-5_m, -5_m, -7.1_m)); Geant::BoxSolid* bottom_chamber = new Geant::BoxSolid("BottomChamber", bottom_chamber_box); bottom_chamber->SetNistMaterial("G4_AIR"); bottom_chamber->Update(); scene.AddSolid(bottom_chamber); Geant::EmitterPrimary* emitter = new Geant::EmitterPrimary(); emitter->SetPosition(Vector3f(0, 0, 14_m)); scene.SetEmitter(emitter); scene.Initialize(); // 2. Setup VTK Viewer Vtk::Viewer viewer; viewer.GetRenderer()->SetBackground(0.05, 0.05, 0.1); // Visualize world box Vtk::ContainerBox* vtkWorld = new Vtk::ContainerBox(scene.GetWorldBox()); vtkWorld->ShowScaleMeasures(true); vtkWorld->SetRepresentation(Vtk::Puppet::Wireframe); vtkWorld->SetSelectable(false); viewer.AddPuppet(*vtkWorld); // Visualize iron cube Vtk::ContainerBox* vtkIron = new Vtk::ContainerBox(&iron_box); vtkIron->SetOpacity(0.2); vtkIron->SetRepresentation(Vtk::Puppet::Surface); viewer.AddPuppet(*vtkIron); // Visualize Top Chamber Vtk::DetectorChamber* vtkTop = new Vtk::DetectorChamber(top_chamber_box); vtkTop->SetOpacity(0.5); vtkTop->SetColor(0.2, 0.8, 0.2); vtkTop->SetRepresentation(Vtk::Puppet::Surface); viewer.AddPuppet(*vtkTop); // Visualize Bottom Chamber Vtk::DetectorChamber* vtkBottom = new Vtk::DetectorChamber(bottom_chamber_box); vtkBottom->SetOpacity(0.5); vtkBottom->SetColor(0.2, 0.8, 0.2); vtkBottom->SetRepresentation(Vtk::Puppet::Surface); viewer.AddPuppet(*vtkBottom); // Visualize Emitter Vtk::EmitterPrimary* vtkEmitter = new Vtk::EmitterPrimary(*emitter); viewer.AddPuppet(*vtkEmitter); // 3. Event Handling AppState state = { &scene, &viewer, {} }; vtkSmartPointer keyCallback = vtkSmartPointer::New(); keyCallback->SetCallback(KeyPressCallbackFunction); keyCallback->SetClientData(&state); viewer.GetInteractor()->AddObserver(vtkCommand::KeyPressEvent, keyCallback); std::cout << "=================================================" << std::endl; std::cout << " Geant Muon Interactive Emitter Test" << std::endl; std::cout << " Use the Widget to move/rotate the Arrow (Emitter)" << std::endl; std::cout << " Press [ENTER] to fire a muon" << std::endl; std::cout << " Press [q] to exit" << std::endl; std::cout << "=================================================" << std::endl; viewer.ZoomAuto(); viewer.Start(); return 0; }