I have modified a Geant Monte Carlo program I set up a year or so ago. Now it includes the Birks' saturation law, light attenuation in green and clear fibers and PMT response (12 stages are simulated).

There are suspicions that GMINOS does not include properly the PMT response and the latter can affect noticably the detector's spectra shapes. Hopefully it will be interesting and important to compare this MC with the real CalDet data and what GMINOS produces.

The output of the MC produces a tree similar to our mapped rooted files with real data (well, no faf errors etc. here of course). So each entry basically has an adc count associated with a particular plane, strip, side etc.

To convert adc to mips I simulated muons. Their pulseheight distributions are shown in this plot.
 


 

I fitted the distributions ("green",  "clear" and "summed" separately) with landau function for one strip  and used the  conversion factor for all the strips.

Then I simulated beam electrons of 5 different energies 1, 1.5, 2, 2.5 and 3 GeV. Again the output is in adc counts and with a small analysis code I extract the total visible energy distribution and use the muon calibration discussed above to convert everything to mip.

This is 1 GeV electrons distiribution:


 

And 2 GeV......


 
 

And 3 GeV.....
 


 
 
 

The detector response versus beam energy for "green", "clear" and "summed" readout separately:


 

The next step is to compare this MC with GMINOS and data. The plot below shows the total visible energy distribution
for 2 GeV electron sample obtained with my MC and GMINOS (provided by Matthew). As GMINOS file was available only in npes I just scaled my plot.


 

And here's the plot comparing this MC with data and GMINOS for 2 GeV positrons beam.  Again I had to scale the x-axis to match the means so "MIP" is arbitrary in this plot. The positron data showed below has been provided by Tricia.  She is using two sorts of cuts, detector cuts (dead chips, cerenkov, plane trigger, no faf and sparsification errors etc) and computational cuts(first plane cut,
number of hit planes, total mip in the first plane etc). I'm comparing here the data with the detector cuts only (NO computational cuts).

For other beam energies I don't have GMINOS files at the moment but I do the data files analysed by Tricia (only detector cuts again).