Basic Information About VEMCalAnalyzer.cxxΒΆ

Produces Muon Spectrum for HG DOMs, fits the spectrum, produces the VEM value for HG DOMs Compares HG and LG DOMs in same tank, calibrates LG DOM wrt HG using the linear overlap region.

Listed below is each component of the VEMCalAnalyzer.cxx script, with a short explanation

hglgDiff(x_axis,profile_hist_params)

HG/LG charge ratio function

muonSpec(x_axis,profile_hist_params)

Muon spectrum function

VEMCalAnalyzer::VEMCalAnalyzer()

VEMCalAnalyzer::OpenFiles(filelist)

VEMCalAnalyzer::AddMuonHist(omkey_hg,roothist_1D)

VEM histograms for plotting

VEMCalAnalyzer::AddHGLGHist(omkey_hg,roothist_2D)

HGLG correlation for plotting

VEMCalAnalyzer::Analyze(outdir_for_rootfiles)

FitMuonSpectrum FitHGLGCorrelation

VEMCalAnalyzer::Smooth(roothist_1D,bin_range)

Removes negative values, and smooths histogram

VEMCalAnalyzer::ExpStartValues(roothist_1D,x1,x2,norm,slope,offset)

Calc start for histograms

VEMCalAnalyzer::FindPeakBin(roothist_1D,xmin,xmax,minimun_value,nsigma)

VEMCalAnalyzer::FitMuonSpectrum(roothist_1D,peak,width,,ratio_signal_bkg,Chi2)

Smooth, FindPeakBin Except values of <25pe, fit with gaussian Find peak position, width and amplitude ExpStartValues Add fit functions to histograms 1 VEM = 95% of peak value position GetChi2

VEMCalAnalyzer::ProfileX(roothist_2D)

Take 2D histogram ad make a 1D weighted histogam

VEMCalAnalyzer::HGLGStartValues(roothist_1D,x1,x2,x0,y0,p2,p3)

Calc mean start time from profile of 2D root histogram

VEMCalAnalyzer::FitHGLGCorrelation(roothist_2D,lgCorr,hgCorr,Chi2)

Take HGLG histogram Find difference weighted with errors Look for best fit with chi2 values

VEMCalAnalyzer::Summarize(rootfile,VEMCalResultMap)

All chi2, VEM, crossover printed to frame Check goodness fits

VEMCalAnalyzer::FillVEMHistograms(VEMCalResultMap,omkey,hg_minVEM,omkey,hg_maxVEM,omkey,hg_minVEM,omkey,hg_maxVEM)

VEMCalAnalyzer::FillCrossOverHistogram(VEMCalResultMap,omkey,min_crossover,omkey,max_crossover)

VEMCalAnalyzer::FillRChi2Histograms(VEMCalResultMap,omkey,min_chi2_muon,omkey,max_chi2_muon,omkey,min_chi2_hglg,omkey,max_chi2_hglg)

VEMCalAnalyzer::WriteXML(VEMCalResultMap,true_if_good)