Handle layered index of refraction

[mjtueros]

• this requires:
• function to compute refraction index
• logic to read from zhaires and save into rawroot
• logic to read from rawroot and save with sim2root