Pro tvGeant4Image,date=date,SingleImage=SingleImage

If (not keyword_set(SingleImage)) then SingleImage = 1
If (not keyword_set(Pitch)) then Pitch = 15
If (not keyword_set(Thickness)) then Thickness = 25
If (not keyword_set(NEvents)) then NEvents = 400000
If (not keyword_set(PType)) then Ptype = 'Muons'
If (not keyword_set(OpticsFlag)) then OpticsFlag = 'OpticsOff'

;PRIMITIVE SUBTRACTION AND FLAT FIELD
GEANT4Dir='/nfs/slac/g/ki/ki08/lsst/lsstgeant4/'
FitsFileName=GEANT4Dir+strtrim(Pitch,2)+'x'+strtrim(Pitch,2)+'x'+$
	     strtrim(Thickness,2)+'Microns'+strtrim(NEvents,2)+$
	     PType+OpticsFlag+'.fits'
fits_read,FitsFileName,Im1,header1
Im1Hist=histogram(Im1,Nbins=100,max=10000,locations=locIm1)
window,xsize=768,ysize=768
tvscl,congrid(Im1,768,768,/center)
stop

end
