Fixed bad values for gd-scint in RAT

From: Christopher Tunnell (tunnell@mail.utexas.edu)
Date: Sat Dec 03 2005 - 15:50:37 CST


During the phone meeting, it was suggested by the KSU group that the
commit I made that placed all material properties in RATDB changed the
physics of the gd-loaded scintillator. After reviewing my work, it
seems that they were correct and I had made a mistake when converting
data to RATDB.

I have submitted a new commit which should fix this. Below I document
how I got the values I did for percentages of the components of the
gadolinium loaded scintillator. Also, attached is an IBD reconstructed
energy plot which should show things work with the new commit.

Thanks to the KSU-team for spotting this!

Just a note though, there is documentation on how to make changes to the
materials information in the RAT wiki.

Chris

Old BWDetectorConstruction:

  // Define the Gd-doped scintillator (copied from GLG4DetectorConstruction
  // with the CHOOZ numbers added)
  G4double density = 0.846*g/cm3;
  G4int nelements = 4;
  GdLoadedScint = new G4Material(name="GdLoadedScint",density,nelements);

  G4double PPO_fraction= 1.5*g/(1e3*cm3*density); // 1.5 g/l
  GdLoadedScint->AddMaterial(_mineralOil, 0.80 / (1.0+PPO_fraction) );
  GdLoadedScint->AddMaterial(Pseudocumene, 0.198 / (1.0+PPO_fraction) );
  GdLoadedScint->AddMaterial(_materialGd, 0.002 / (1.0+PPO_fraction) );
  GdLoadedScint->AddMaterial(PPO, PPO_fraction / (1.0+PPO_fraction) );

Check code I added to RAT to see what the proper weightings were:

  G4double density = 0.846*g/cm3;
  G4double pf= 1.5*g/(1e3*cm3*density);
  G4cout<<"CATCHME "<<pf<<" "<<0.8/(1+pf)<<" "<<0.198/(1+pf)<<"
"<<0.002/(1+pf)<<" "<<pf/(1+pf)<<G4endl;

Output on console:
CATCHME 0.00177305 0.798584 0.19765 0.00199646 0.00176991



c1.jpg

This archive was generated by hypermail 2.1.6 : Wed Dec 14 2005 - 00:01:02 CST