LUT Davis, Odd behavior

Starting off,

Geant4 Version: 11.1.0
OS: Ubuntu 20.04.6
cmake version: 3.16.3

My goal has been to duplicate the simulation results of the LUT Davis benchmark/implementation paper, “An integrated model of scintillator-reflector properties for advanced simulations of optical transport - PubMed” referenced in the manual. My issue is that I am unable to duplicate the results, I think the issue is due to my definition of the LUT Davis optical surface.

The results from the paper I am comparing against are the polished no reflector case, shown in black in Fig. 12., I digitized the results and plotted against my own.


In green is the digitized data from the paper, and the other two curves are from when I used the Unified model and the LUT Davis model. I included the UNIFIED model case because there should be little difference between the bare polished UNIFIED model and the bare polished LUT Davis model.

Because my results when using the UNIFIED model are close to the paper result, I believe I can rule out the possibility that the problem is with my model geometry or physics. Hence, I am inclined to lay the blame on my definition of the LUT Davis surface.

So my questions are:

  1. Is there anything wrong with my wrong with my definition of the LUT_Davis surface?
  2. If there is nothing wrong with my surface definition what else could it be?
  3. Exhausting all other possibilities, could there be an issue with how newer versions of Geant4 implements the LUT Davis surfaces?

‘’’
G4OpticalSurface* LUT_Surface = new G4OpticalSurface(“LUT_Surface”);
LUT_Surface->SetType(dielectric_LUTDAVIS);
LUT_Surface->SetFinish(Polished_LUT);
LUT_Surface->SetModel(DAVIS);

G4LogicalBorderSurface *CrystaltoWorld  = new G4LogicalBorderSurface("CrystaltoWorld", 
                                                                     physCrystal, 
                                                                     physWorld, 
                                                                     LUT_Surface);

‘’’

Construction.cc (11.9 KB)
Construction.hh (959 Bytes)