README explaining the analysis chain
python script to produce a chi2 plot as a function of the varied parameter
KryptonCheck is used to get an idea about the systematics of the Krypton data curve, by using MaximumBin, Mean or a Fit. In the end a fit was used, and the mean and mean uncertainty were used further in analyser.exe to calculate the chi2 between data and MC.
updated KryptonCheck.C, added Fit of variable as a function of z -> copied it then over to analyser.C/.h
developed the code further into a class. Store now in an ntuple and produce plots from that output ntuple. Development continues.
Paolo's original script, from which I started developing. The modification of the values is a bit odd.
refill the TopBottom 2D map from top/bottom top/(top+bottom) and changed the axes to better match what I need to plot.
changes in the eventcounter and the log (from trace to routine)
- considerable reorganization of the code
- now it works to reproduce the plots of Paolo and Davide :)
- includes now also uncertainties -> precondition for chi2 calculation.
- added also Refill_MCTuning_Data.C to refill the Krypton-data for TopBottom in data, rather than doing this in the analyser.C/.h. Includes shifting of the z-axis and mirroring at the y-axis.
- separated the loading of the MC tree, the calculation of the CollEff and the TopBottom ratio into different functions.
- introduced also MC_CollEff_Optical() next to MC_CollEff() respectively MC_TopBottom_Optical() next to MC_TopBottom(), since the optical photon production (where one photon is simulated per event will be considerably different than e.g. a Krypton MC as far as e.g. the top/bottom asymmetry is concerned.
comment on refilling, which is done in Refill_MCTuning_Data.C
comment on optical MC, which is treated differently than e.g. Kr85. New default statistics: 10^7 photons.
changed from 10^6 to 10^7 optical photons per default. Changed from FusedSilicaUVAbs to GridSteelRindScale. These are minor changes.
Now a preliminary chi2 calculation is done. Including errors and everything.
member variables to hold chi2 and ndf for easier output. It works :)
Merge branch 'master' of https://github.com/berndreinhold/DS_Hawaii