Gasfile issue during running

ahmer_k · March 22, 2021, 3:24pm

@hschindl i was working on it and producing different gas files. now u want to produce another gas file by using example gasfile in garfield. but it is giving an error and not running see. as i have put name and parameter in generate.cc. can you help what is issue?

ahmer_k · March 22, 2021, 3:47pm

dear another thing i’m replacing the gas file in the example geant4garfieldinterface in the file detectorconstruction.cc but its not compiling what is issue here i want you to help me in this regard too. please waiting for your answwer.

hschindl · March 22, 2021, 3:48pm

The parameters of the function SetComposition are (in this order)

name of the first component in the mixture (string),
fraction of the first component in the mixture (double),
name of the second component (string),
fraction of the second component (double),
and so on (up to six components).
So, in your case:

gas.SetComposition("C2H4F2", 70., "C4H10", 30.);

hschindl · March 22, 2021, 3:59pm

Hard to say without seeing your code, but it looks like C2F4H2 is a G4Element object while it should be a G4Material (first error) and you are using an object (Ar) that doesn’t exist.

ahmer_k · March 24, 2021, 3:49pm

oook dear i’ll show you more so you can explain if there is something else.

ahmer_k · March 25, 2021, 4:15pm

here the detector construction.cc with in the geant4garfieldinterface example of garfield. and i have tried to change and then error comes in make command which i have asked you previous day. i have made changes but why is errors?

ahmer_k · March 25, 2021, 4:24pm

in this above pic you can see i have put the changes with c2f4h2 100% and with its gas file but still issues in make command.

hschindl · March 25, 2021, 4:46pm

C2F4H2 is not an element, so I guess the line

nistManager->findOrBuildElement("C2F4H2", isotopes);

won’t work. I think you have to build a C2F4H2 G4Material the same way you do it for CO2.

Another issue, as I mentioned in my previous reply, is in the line

ArCO2_93_7->AddElement(Ar, fractionmass = 0.93);

It causes an error because there is no object called Ar in your function. This line (and the one before/after) is probably a left over from the original code (which did use Ar/CO2).

By the way, it would be great if (in the future) you could post your code as text instead of a picture.

ahmer_k · March 26, 2021, 5:38am

ok dear i’ll check then let you know. thanks a lot.

ahmer_k · March 29, 2021, 2:24pm

ahmer@ahmer-Dell-Precision-M3800 ~ cd Desktop ahmer@ahmer-Dell-Precision-M3800 ~/Desktop cd garfeidl
bash: cd: garfeidl: No such file or directory
ahmer@ahmer-Dell-Precision-M3800 ~/Desktop cd garfield ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield cd Examples/GasFile
ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile ls ar_80_co2_20_0T.gas ar_80_co2_20_2T.gas CMakeLists.txt generate.C merge.C read.C ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile mkdir build
ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile cd build ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile/build source /home/ahmer/Desktop/Geant4/geant4.10.06.p01-install/bin/geant4.sh
ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile/build source /home/ahmer/Desktop/root/test_build/bin/thisroot.sh ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile/build source /home/ahmer/Desktop/garfield/build/setup.sh
ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile/build cmake .. -- The C compiler identification is GNU 5.4.0 -- The CXX compiler identification is GNU 5.4.0 -- Check for working C compiler: /usr/bin/cc -- Check for working C compiler: /usr/bin/cc - works -- Detecting C compiler ABI info -- Detecting C compiler ABI info - done -- Detecting C compile features -- Detecting C compile features - done -- Check for working CXX compiler: /usr/bin/c++ -- Check for working CXX compiler: /usr/bin/c++ - works -- Detecting CXX compiler ABI info -- Detecting CXX compiler ABI info - done -- Detecting CXX compile features -- Detecting CXX compile features - done -- Configuring done -- Generating done -- Build files have been written to: /home/ahmer/Desktop/garfield/Examples/GasFile/build ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile/build make -j8
Scanning dependencies of target read
Scanning dependencies of target generate
[ 25%] Building CXX object CMakeFiles/generate.dir/generate.C.o
[ 50%] Building CXX object CMakeFiles/read.dir/read.C.o
[ 75%] Linking CXX executable generate
[100%] Linking CXX executable read
[100%] Built target generate
[100%] Built target read
ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile/build ls ar_80_co2_20_0T.gas ar_80_co2_20_2T.gas CMakeCache.txt CMakeFiles cmake_install.cmake generate Makefile read ahmer@ahmer-Dell-Precision-M3800 ~/Desktop/garfield/Examples/GasFile/build ./generate
Heed:
Database path: /home/ahmer/Desktop/garfield/Heed/heed++/database
MediumMagboltz::GetGasName:
Gas G4_FREON-12 is not recognized.
MediumMagboltz::SetComposition:
Error setting the composition. No valid components.
MediumMagboltz::GenerateGasTable: Found 0 excitations and 0 ionisations.
MediumMagboltz::GenerateGasTable: E = 100 V/cm, B = 0 T, angle: 1.5708 rad

 ERROR IN GAS INPUT : NGAS=    0

N= 1 NGAS= 0 FRAC= 0.000
N= 2 NGAS= 0 FRAC= 0.000
N= 3 NGAS= 0 FRAC= 0.000
N= 4 NGAS= 0 FRAC= 0.000
N= 5 NGAS= 0 FRAC= 0.000
N= 6 NGAS= 0 FRAC= 0.000
MediumMagboltz::RunMagboltz: Results:
Drift velocity along E: 0.00000000 cm/ns +/- 0.00%
Drift velocity along Bt: 0.00000000 cm/ns +/- 0.00%
Drift velocity along ExB: 0.00000000 cm/ns +/- 0.00%
Lorentz angle: 0.000 degree
Longitudinal diffusion: -nan cm1/2 +/- 0.00%
Transverse diffusion: -nan cm1/2 +/- 0.00%
Townsend coefficient: 0.0000 cm-1 +/- 0.00%
Attachment coefficient: 0.0000 cm-1 +/- 0.00%
MediumMagboltz::GenerateGasTable: E = 143.845 V/cm, B = 0 T, angle: 1.5708 rad

ahmer_k · March 29, 2021, 2:26pm

in this down pasted code i have tried to build the gas file for 100% freon.
but as mentioned name in the list of materials is G4_FREON-12 but the code is not picking it as you are seeing in the code the errors how to fix it?

hschindl · March 29, 2021, 7:11pm

Hi,
G4_FREON-12 is not a valid identifier for C₂H₄F₂. You can use one of the following:

"C2H2F4", "C2F4H2",
"FREON 134", "FREON 134A", "FREON-134", "FREON-134-A",
"TETRAFLUOROETHANE"

ahmer_k · March 30, 2021, 3:39pm

ok dear i’ll try with others. but G4_FREON-12 was in the /material/nist/list material. i have checked then why is it not working?

ahmer_k · March 30, 2021, 3:43pm

and dear can you send me the material list so i can also check material and its proper font.? my email. ahmer.khan713@gmail.com

hschindl · March 30, 2021, 4:09pm

That means that Geant4 has a material with that name in its list of materials, but that is completely separate from Garfield++/Magboltz. To set the gas mixture to C₂H₄F₂ in MediumMagboltz, you need to use one of the labels/names I mentioned in my previous message.

By the way, I suspect that G4_FREON-12 might correspond to Dichlorodifluoromethane - Wikipedia and not 1,1,1,2-Tetrafluoroethane - Wikipedia.

hschindl · March 30, 2021, 4:12pm

Do you mean a list of gases that are available in Magboltz?
See Appendix B of the Garfield++ user guide.

ahmer_k · March 30, 2021, 5:12pm

OK dear thanks i’ll see that too.

ahmer_k · March 30, 2021, 5:33pm

in the above picture now i’m basically introduce my tetrafluoroethane gas file in the geant4garfieldinterface example. this picture is basically the picture of the detectorconstruction.cc src file. and also i have copied the generated file of c2h2f4 in the geant4garfieldinterface folder no other file is it right? kindly if you can tell me how to do it?

ahmer_k · March 30, 2021, 5:39pm

and also as i have to place the density, ionizing potential, and n-components? how to find that?

hschindl · March 31, 2021, 9:59am

I’m not a Geant4 expert, but I think you could do something like this:

G4Element* H = nistManager->FindOrBuildelement("H", isotopes);
G4Element* C = nistManager->FindOrBuildelement("C", isotopes);
G4Element* F = nistManager->FindOrBuildelement("F", isotopes);
 
G4Material* C2H2F4 = new G4Material("C2H2F4", 
    density = ... * CLHEP::mg / CLHEP::cm3, ncomponents = 3, kStateGas);
 C2H2F4->AddElement(C, natoms = 2);
 C2H2F4->AddElement(H, natoms = 2);
 C2H2F4->AddElement(F, natoms = 4);

I don’t know the density off hand, so you’ll need to replace the three dots by the actual number.

Incidentally, I would recommend to tidy up the function DefineMaterials. You define materials (e. g. Ar/CO₂ 70:30, Ar/CO₂ 93:7) in this function which seem to be remnants of the original code and which I believe you don’t use.
Also note that this function defines only the materials on the Geant4 side of your program, it doesn’t do anything related to Garfield.

Regarding the density: you can assume that it is an ideal gas, so the number density at 0° C and atmospheric pressure is N = 2.6867774e19 cm^-3.
To get the mass density you need to multiply this number by the atomic weight A of the molecule and the atomic mass unit.

I think I asked this question before, but are you sure you need to interface Geant4 and Garfield++ for your study and cannot simply use Garfield++ only?