Dear Paola,
thanks for your prompt reply. Indeed, this was my understanding too.
However, I found similar issue in defining a compound material like BGO
either
a = 16.00*g/mole;
G4Element* elO = new G4Element(name="Oxygen" ,symbol="O" , z= 8., a);
a = 71.64*g/mole;
G4Element* elGe = new G4Element(name="Germanium" ,symbol="Ge" , z=
32., a);
a = 208.98040*g/mole;
G4Element* elBi = new G4Element(name="Bismuth" ,symbol="Bi" , z=
83., a);
density = 7.31*g/cm3;
G4Material* BGO = new G4Material("BGO", density=density, ncomponents=3);
BGO->AddElement(elO , natoms=12);
BGO->AddElement(elGe, natoms= 3);
BGO->AddElement(elBi, natoms= 4);
or
a = 16.00*g/mole;
density=0.00133151*g/cm3;
G4Material *matO = new G4Material(name="matO" , z= 8., a, density);
a = 71.64*g/mole;
density=5.323*g/cm3;
G4Material *matGe = new G4Material(name="matGe" , z= 32., a, density);
a = 208.98040*g/mole;
density=9.747*g/cm3;
G4Material *matBi = new G4Material(name="matBi" , z= 83., a, density);
density = 7.31*g/cm3;
G4Material* myBGO = new G4Material("myBGO", density=density,
ncomponents=3);
myBGO->AddMaterial(matO , fractionmass=0.154485);
myBGO->AddMaterial(matGe, fractionmass=0.172926);
myBGO->AddMaterial(matBi, fractionmass=0.672589);
I have just included the above lines in the MyDetectorConstruction.cc
for the AlAuAl example.
In the first case the material output files is alaual_flukaMat_run.inp,
while in the second case it is alaual_flukaMat_run1.inp (the files are
attached).
Do you have any suggestion in this particular case?
My best,
Nicola
Il 02/10/2013 12:39, paola sala ha scritto:
> Dear Nicola,
> please do not use the Nist Manager to build materials, use G4Material
> instead
> Paola
> On 10/01/2013 10:01 AM, M. Nicola Mazziotta wrote:
>> Dear all,
>>
>> I am using the FLUGG tool to interface the GEANT4 geometry to FLUKA
>> code. However, I found an issue using the G4 material due to "NON
>> SEQUENTIAL MATERIAL NUMBER" message error in running FLUKA.
>>
>> To easily see this message is enough to made some changes in the
>> AlAuAl example changing the Al and Au material definitions in the
>> MyDetectorConstruction.cc as:
>>
>> ...
>> #include "G4NistManager.hh"
>> ...
>> ...
>>
>> density = 2.700*g/cm3;
>> a = 26.98*g/mole;
>> //G4Material* Al = new G4Material(name="Aluminum", z=13., a, density);
>>
>> G4NistManager* man = G4NistManager::Instance();
>> G4Material* Al = man->FindOrBuildMaterial("G4_Al");
>>
>> density = 19.300*g/cm3;
>> a = 196.97*g/mole;
>> //G4Material* Au = new G4Material(name="Gold", z=79., a, density);
>> G4Material* Au = man->FindOrBuildMaterial("G4_Au");
>>
>> I am using the latest version of FLUKA 2011.2b.4, the latest version
>> of FLUGG flugg_2009_4 and the Geant4 version 4.9.6.p02.
>>
>> Any suggestions would be greatly appreciated.
>>
>> Best regards,
>> Nicola
>>
>
Received on Wed Oct 02 2013 - 17:13:42 CEST