Hello
I Want to use Magboltz to write custom gas files, or I want to use Magboltz to do the analysis of a proportional counter properly. I think while installing the Garfield ++ the Magboltz didn’t installed properly because i am not able to generate any gas file.
Can anyone help me rectifying this issue
I am using Ubuntu 22.04.5 LTS
Thanks
Hi,
can you give more details? What do you mean by “not able to generate any gas file”? Did you try running this example?
What output did you get from the program?
Hi
Thank you for replying, Actually i am trying to use the gas properties in my Geant4 setup. Now in the included example i wasn’t able to figure out where and how to use magboltz for the same. I am trying to implement electric field properties for which i needed Garfield ++ but i am stuck.
Can you suggest if there is any way to simulate energy deposition of charged particles in electric field.
Thanks
I’m afraid I’m not sure I understand your question. You mention a few different things in your post(s):
All of these are ingredients in a typical Garfield++ simulation, but they are distinct things. Which of them is your question about?
Hi
Actually, I have to do all three analyses that you have mentioned. After your suggestion I have used Magboltz perfectly, there were some library linking issues that I was facing earlier, and now I have rectified that. I have played with the transport properties one and now the gas file is generating. I will try Heed this week, and if i face any issue i will ask on this thread only.
Thanks
Hi again, I have successfully modified the construction file for the Geant4Garfield ++ interface code according to my requirements and ran it to see the visualisation. It was working and then I started to modify other files as well.
When I modified the GarfieldPhysics.cc file and then tried to run the executable I got the following error
=========================================================================
ujjwal@Jarvis:~/Geant4GarfieldInterface/build$ ./exampleGeant4Interface
Available UI session types: [ Qt, tcsh, csh ]
Warning: Ignoring XDG_SESSION_TYPE=wayland on Gnome. Use QT_QPA_PLATFORM=wayland to run on Wayland anyway.
RandomEngineRoot::Seed: 1
Geant4 version Name: geant4-11-02-patch-02 [MT] (21-June-2024)
Copyright : Geant4 Collaboration
References : NIM A 506 (2003), 250-303
: IEEE-TNS 53 (2006), 270-278
: NIM A 835 (2016), 186-225
WWW : http://geant4.org/
<<< Geant4 Physics List simulation engine: QGSP_BERT_HP
RegisterPhysics: G4EmStandard
RegisterPhysics: G4GammaLeptoNuclearPhys
RegisterPhysics: Decay
RegisterPhysics: G4RadioactiveDecay
RegisterPhysics: nElasticPhysics_HP
RegisterPhysics: hInelastic QGSP_BERT_HP
RegisterPhysics: stopping
RegisterPhysics: ionInelasticFTFP_BIC
RegisterPhysics: IonElasticPhysics
Using
Visualization Manager instantiating with verbosity “warnings (3)”…
Visualization Manager initialising…
Registering graphics systems…
You have successfully registered the following graphics systems.
Registered graphics systems are:
ASCIITree (ATree)
DAWNFILE (DAWNFILE)
G4HepRepFile (HepRepFile)
RayTracer (RayTracer)
VRML2FILE (VRML2FILE)
gMocrenFile (gMocrenFile)
TOOLSSG_OFFSCREEN (TSG_OFFSCREEN)
TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE)
OpenGLImmediateQt (OGLIQt, OGLI)
OpenGLStoredQt (OGLSQt, OGL, OGLS)
RayTracerX (RayTracerX)
TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG)
TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB)
Default graphics system is: OGL (based on build flags).
Default window size hint is: 600x600-0+0 (based on G4VisManager initialisation).
Note: Parameters specified on the command line will override these defaults.
Use “vis/open” without parameters to get these defaults.
Registering model factories…
You have successfully registered the following model factories.
Registered model factories:
generic
drawByAttribute
drawByCharge
drawByOriginVolume
drawByParticleID
drawByEncounteredVolume
Registered models:
None
Registered filter factories:
attributeFilter
chargeFilter
originVolumeFilter
particleFilter
encounteredVolumeFilter
Registered filters:
None
You have successfully registered the following user vis actions.
Run Duration User Vis Actions: none
End of Event User Vis Actions: none
End of Run User Vis Actions: none
Some /vis commands (optionally) take a string to specify colour.
“/vis/list” to see available colours.
/control/saveHistory
/run/verbose 2
#/run/numberOfThreads 12
/control/execute physics.mac
/exampleGarfield/physics/setIonizationModel PAI false
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy e- 0 1e+8 MeV
PAI is applicable for G4Particle e- between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy e+ 0 1e+8 MeV
PAI is applicable for G4Particle e+ between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy mu- 0 1e+8 MeV
PAI is applicable for G4Particle mu- between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy mu+ 0 1e+8 MeV
PAI is applicable for G4Particle mu+ between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy proton 0 1e+8 MeV
PAI is applicable for G4Particle proton between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy pi+ 0 1e+8 MeV
PAI is applicable for G4Particle pi+ between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy pi- 0 1e+8 MeV
PAI is applicable for G4Particle pi- between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy alpha 0 1e+8 MeV
PAI is applicable for G4Particle alpha between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy He3 0 1e+8 MeV
PAI is applicable for G4Particle He3 between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGeant4ParticleTypeAndEnergy GenericIon- 0 1e+8 MeV
PAI is applicable for G4Particle GenericIon- between 0 MeV and 1e+08 MeV
/exampleGarfield/physics/setGarfieldParticleTypeAndEnergy gamma 1e-6 1e+8 MeV
Garfield model (Heed) is applicable for G4Particle gamma between 1e-06 MeV and 1e+08 MeV
/exampleGarfield/physics/setGarfieldParticleTypeAndEnergy e- 1e-6 1e-3 MeV
Garfield model (Heed) is applicable for G4Particle e- between 1e-06 MeV and 0.001 MeV
/run/initialize
userDetector->Construct() start.
***** Table : Nb of materials = 5 *****
Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.095 m
Imean: 85.700 eV temperature: 293.15 K pressure: 1.00 atm
—> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole
—> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 %
—> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 %
ElmMassFraction: 0.01 % ElmAbundance 0.02 %
—> Element: N (N) Z = 7.0 N = 14 A = 14.007 g/mole
—> Isotope: N14 Z = 7 N = 14 A = 14.00 g/mole abundance: 99.632 %
—> Isotope: N15 Z = 7 N = 15 A = 15.00 g/mole abundance: 0.368 %
ElmMassFraction: 75.53 % ElmAbundance 78.44 %
—> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole
—> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
—> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
—> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
ElmMassFraction: 23.18 % ElmAbundance 21.07 %
—> Element: Ar (Ar) Z = 18.0 N = 40 A = 39.948 g/mole
—> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.337 %
—> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.063 %
—> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.600 %
ElmMassFraction: 1.28 % ElmAbundance 0.47 %
Material: G4_STAINLESS-STEEL density: 8.000 g/cm3 RadL: 1.738 cm Nucl.Int.Length: 16.678 cm
Imean: 282.977 eV temperature: 293.15 K pressure: 1.00 atm
—> Element: Fe (Fe) Z = 26.0 N = 56 A = 55.845 g/mole
—> Isotope: Fe54 Z = 26 N = 54 A = 53.94 g/mole abundance: 5.845 %
—> Isotope: Fe56 Z = 26 N = 56 A = 55.93 g/mole abundance: 91.754 %
—> Isotope: Fe57 Z = 26 N = 57 A = 56.94 g/mole abundance: 2.119 %
—> Isotope: Fe58 Z = 26 N = 58 A = 57.93 g/mole abundance: 0.282 %
ElmMassFraction: 74.62 % ElmAbundance 74.00 %
—> Element: Cr (Cr) Z = 24.0 N = 52 A = 51.996 g/mole
—> Isotope: Cr50 Z = 24 N = 50 A = 49.95 g/mole abundance: 4.345 %
—> Isotope: Cr52 Z = 24 N = 52 A = 51.94 g/mole abundance: 83.789 %
—> Isotope: Cr53 Z = 24 N = 53 A = 52.94 g/mole abundance: 9.501 %
—> Isotope: Cr54 Z = 24 N = 54 A = 53.94 g/mole abundance: 2.365 %
ElmMassFraction: 16.90 % ElmAbundance 18.00 %
—> Element: Ni (Ni) Z = 28.0 N = 59 A = 58.693 g/mole
—> Isotope: Ni58 Z = 28 N = 58 A = 57.94 g/mole abundance: 68.077 %
—> Isotope: Ni60 Z = 28 N = 60 A = 59.93 g/mole abundance: 26.223 %
—> Isotope: Ni61 Z = 28 N = 61 A = 60.93 g/mole abundance: 1.140 %
—> Isotope: Ni62 Z = 28 N = 62 A = 61.93 g/mole abundance: 3.635 %
—> Isotope: Ni64 Z = 28 N = 64 A = 63.93 g/mole abundance: 0.926 %
ElmMassFraction: 8.48 % ElmAbundance 8.00 %
Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.312 cm
Imean: 727.000 eV temperature: 293.15 K pressure: 1.00 atm
—> Element: W (W) Z = 74.0 N = 184 A = 183.842 g/mole
—> Isotope: W180 Z = 74 N = 180 A = 179.95 g/mole abundance: 0.120 %
—> Isotope: W182 Z = 74 N = 182 A = 181.95 g/mole abundance: 26.500 %
—> Isotope: W183 Z = 74 N = 183 A = 182.95 g/mole abundance: 14.310 %
—> Isotope: W184 Z = 74 N = 184 A = 183.95 g/mole abundance: 30.640 %
—> Isotope: W186 Z = 74 N = 186 A = 185.95 g/mole abundance: 28.430 %
ElmMassFraction: 100.00 % ElmAbundance 100.00 %
Material: CH4 density: 0.667 mg/cm3 RadL: 696.640 m Nucl.Int.Length: 907.478 m
Imean: 45.542 eV temperature: 293.15 K pressure: 1.00 atm
—> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole
—> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 %
—> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 %
ElmMassFraction: 74.87 % ElmAbundance 20.00 %
—> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
—> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
—> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
ElmMassFraction: 25.13 % ElmAbundance 80.00 %
Material: P10 density: 1.500 mg/cm3 RadL: 138.339 m Nucl.Int.Length: 726.668 m
Imean: 155.640 eV temperature: 293.15 K pressure: 1.00 atm
—> Element: Ar (Ar) Z = 18.0 N = 40 A = 39.948 g/mole
—> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.337 %
—> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.063 %
—> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.600 %
ElmMassFraction: 90.00 % ElmAbundance 41.96 %
—> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole
—> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 %
—> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 %
ElmMassFraction: 7.49 % ElmAbundance 11.61 %
—> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
—> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
—> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
ElmMassFraction: 2.51 % ElmAbundance 46.43 %
Checking overlaps for volume PRCOuter:0 (G4Box) … OK!
Checking overlaps for volume PRCInner:0 (G4Box) … OK!
Checking overlaps for volume AnodeWire:0 (G4Tubs) … OK!
MediumMagboltz::SetComposition: Ar/CH4 (90/10)
MediumMagboltz::Mixer:
4000 linear energy steps between 0 and 40 eV.
ARGON ANISOTROPIC 2014�� 3n�16I�
mass: 39.948 amu
ionisation threshold: 15.7596 eV
44 inelastic terms (44 excitations, 0 superelastic, 0 other)
CH4 2015 ANISOTROPIC �� 3n�16I�
mass: 16.0426 amu
ionisation threshold: 12.65 eV
34 inelastic terms (25 excitations, 2 superelastic, 7 other)
MediumMagboltz::Mixer:
Lowest ionisation threshold in the mixture: 12.65 eV (CH4)
MediumMagboltz::Mixer:
Energy [eV] Collision Rate [ns-1]
2.50 1122.23
7.50 3517.83
12.50 6805.63
17.50 10411.14
22.50 13030.24
27.50 14339.80
32.50 13896.79
37.50 13164.04
MediumMagboltz::ComputePhotonCollisionTable:
Energy [eV] Mean free path [um]
1.00 ------------
3.00 ------------
5.00 ------------
7.00 ------------
9.00 929.7192
11.00 212.0893
13.00 83.0227
15.00 77.8212
17.00 11.1085
19.00 10.4533
MediumMagboltz::EnablePenningTransfer:
Warning: present gas table has no ionisation rates.
Ignore this message if you are using microscopic tracking only.
MediumMagboltz::EnablePenningTransfer:
Penning transfer probability for 44 Ar excitation levels set to r = 0.57
MediumMagboltz::LoadGasFile:
Cannot open file p10_modified.gas.
I am attaching the construction as well as the physics list file please have a look
Thanks
GarfieldDetectorConstruction.cc (8.6 KB)
GarfieldPhysics.cc (14.0 KB)
Cannot open file p10_modified.gas.
Looks like this file is not in your working directory.
Thankyou for pointing this out, After adding this and linking it in CMakeLists.txt file I got the following error
=========================================================================
MediumMagboltz::EnablePenningTransfer:
Warning: present gas table has no ionisation rates.
Ignore this message if you are using microscopic tracking only.
MediumMagboltz::EnablePenningTransfer:
Penning transfer probability for 44 Ar excitation levels set to r = 0.57
MediumMagboltz::LoadGasFile:
Reading file p10_modified.gas.
Version 12.
Gas composition set to Ar/CH4 (90/10).
*** Break *** segmentation violation
=========================================================================
Can you highlight what modifications you made to the code?
Yeah sure
from (Orignal code)
void GarfieldPhysics::InitializePhysics() {
// Define the gas mixture.
fMediumMagboltz = new Garfield::MediumMagboltz();
fMediumMagboltz->SetComposition("ar", 70., "co2", 30.);
fMediumMagboltz->SetTemperature(293.15);
fMediumMagboltz->SetPressure(760.);
fMediumMagboltz->Initialise(true);
// Set the Penning transfer efficiency.
const double rPenning = 0.57;
const double lambdaPenning = 0.;
fMediumMagboltz->EnablePenningTransfer(rPenning, lambdaPenning, "ar");
fMediumMagboltz->LoadGasFile("ar_70_co2_30_1000mbar.gas");
fComponentAnalyticField = new Garfield::ComponentAnalyticField();
fComponentAnalyticField->SetMedium(fMediumMagboltz);
// Wire radius [cm]
constexpr double rWire = 25.e-4;
// Tube radius [cm]
constexpr double rTube = 1.451;
// Voltages
constexpr double vWire = 1000.;
constexpr double vTube = 0.;
// Add the wire in the center.
fComponentAnalyticField->AddWire(0., 0., 2 * rWire, vWire, "w");
// Add the tube.
fComponentAnalyticField->AddTube(rTube, vTube, 0, "t");
fSensor = new Garfield::Sensor();
fSensor->AddComponent(fComponentAnalyticField);
fTrackHeed = new Garfield::TrackHeed();
fTrackHeed->SetSensor(fSensor);
fTrackHeed->EnableDeltaElectronTransport();
}
to (Modified)
void GarfieldPhysics::InitializePhysics() {
// Define the gas mixture (P10: 90% argon, 10% methane)
fMediumMagboltz = new Garfield::MediumMagboltz();
fMediumMagboltz->SetComposition("ar", 90., "ch4", 10.);
fMediumMagboltz->SetTemperature(295.15); // Standard temperature in Kelvin
fMediumMagboltz->SetPressure(814.); // Standard pressure in Torr
fMediumMagboltz->Initialise(true);
// Set the Penning transfer efficiency
const double rPenning = 0.57;
const double lambdaPenning = 0.;
fMediumMagboltz->EnablePenningTransfer(rPenning, lambdaPenning, "ar");
fMediumMagboltz->LoadGasFile("p10_modified.gas");
fComponentAnalyticField = new Garfield::ComponentAnalyticField();
fComponentAnalyticField->SetMedium(fMediumMagboltz);
// Set the detector's wire and tube dimensions
constexpr double rWire = 50.e-6; // Tungsten anode wire radius in cm
constexpr double rTube = 4.77; // PRC tube radius in cm (adjusted for your setup)
constexpr double vWire = 3000.; // Voltage on the anode wire in V
constexpr double vTube = 0.; // Voltage on the tube (grounded)
// Add the wire at the center of the tube
fComponentAnalyticField->AddWire(0., 0., 2 * rWire, vWire, "w");
// Add the tube surrounding the wire
fComponentAnalyticField->AddTube(rTube, vTube, 0, "t");
// Create and configure the sensor
fSensor = new Garfield::Sensor();
fSensor->AddComponent(fComponentAnalyticField);
// Set up the track simulation with Heed
fTrackHeed = new Garfield::TrackHeed();
fTrackHeed->SetSensor(fSensor);
fTrackHeed->EnableDeltaElectronTransport();
}
Thanks
Hi, can anybody help me understand why the Garfield++ ad Geant4 code isn’t triggering an avalanche and showing gain ??
GarfieldDetectorConstruction.cc (7.0 KB)
GarfieldPhysics.cc (14.4 KB)
GarfieldPhysicsList.cc (4.7 KB)
========================================================================
The output is as follows
Event 99 starts.
Primary muon generated at (0, 20.05, 0) m with direction (0, -1, 0)
Step in PRCOuter (Shell), Edep = 2.1999 MeV
Step in PRCInner (Gas), Edep = 0.00569154 MeV
Step in PRCInner (Gas), Edep = 0.00359258 MeV
Step in PRCInner (Gas), Edep = 0.000840207 MeV
Step in PRCInner (Gas), Edep = 0.00163683 MeV
Step in PRCInner (Gas), Edep = 6.46213e-05 MeV
Step in PRCInner (Gas), Edep = 0.00602431 MeV
Step in PRCInner (Gas), Edep = 0.00164082 MeV
Step in PRCOuter (Shell), Edep = 2.52587 MeV
Step in PRCInner (Gas), Edep = 0.0009472 MeV
Step in PRCInner (Gas), Edep = 0.000371256 MeV
Step in PRCInner (Gas), Edep = 0.000164142 MeV
Step in PRCInner (Gas), Edep = 0.000185416 MeV
Step in PRCInner (Gas), Edep = 0.000176215 MeV
Step in PRCInner (Gas), Edep = 0.000370921 MeV
---> End of event: 99
Gas: total energy: 21.7061 keV Wire: total energy: 0 eV Shell: total energy: 4.72578 MeV
Total Detector Energy: 4.74748 MeV Avalanche size: 0 Gain: 0
... write file : Garfield.root - done
... close file : Garfield.root - done
GarfieldRunAction: Closed file Garfield.root
There are histograms that can be viewed with visualization:
7 h1 histograms(s)
List them with "/analysis/list".
View them immediately with "/vis/plot" or "/vis/reviewPlots".
But...there are no entries. To make your histograms available for
plotting in this UI session, use CloseFile(false) in your
EndOfRunAction and Reset() in your BeginOfRunAction.
Graphics systems deleted.
Visualization Manager deleting...
GarfieldPhysics: Instance destructed
GarfieldPhysics: Singleton disposed
This is the run.mac
/exampleGarfield/physics/setIonizationModel Heed true
/run/initialize
/gun/particle mu-
/gun/energy 100 GeV
/gun/position 0 1.05 0 m
/gun/direction 0 -1 0
/run/beamOn 100
Hi,
have you managed to make progress on understanding your issue in the mean time?
I just had a quick look at the code in your GarfieldPhysics.cc and I noticed that you have a line
fMediumMagboltz->SetMaxElectronEnergy(100000.);
This is an extremely (unphysically) high value.
Hi,
I have made significant progress since the last query. Now everything seems to work, but the processing is slow, i want to run 100000 events in this interface to see the energy deposition and then match it with the experimental data. But right now it takes about 2.5 hours just to compute 1000 events.
Thankyou for pointing it out, i will reduce it, i set it up as a placeholder value, Can you suggest what should be the ideal value so that i can check it by running ?
Thanks
This topic was automatically closed 14 days after the last reply. New replies are no longer allowed.