def AnalysisFunction(index): minlap = 10**9 mint1 = 500.0 mint2 = 500.0 for t1 in range(-10, 10, 2): for t2 in range(20, 200, 2): Histogram1 = Utilities.ProduceRatioHistogram( Utilities.Particles[0] + ".root", t1, t2) Histogram2 = Utilities.ProduceRatioHistogram( Utilities.Particles[1] + ".root", t1, t2) overlap = Utilities.BhattacharyaCoeff(Histogram1, Histogram2) del Histogram1 del Histogram2 if overlap < minlap: minlap = overlap mint1 = t1 mint2 = t2 print "\n" print "Please place the text below into the database file: CLASSIFIER_ALPHA_SEEDED.ratdb located in rat/data, replacing any existing entry with the same index." print "\n" print "{" print "type: \"CLASSIFIER_ALPHA_SEEDED\"," print "version: 1," print "index: \"" + index + "\"," print "run_range: [0, 0]," print "pass: 0," print "production: false," print "timestamp: \"\"," print "comment: \"\"," print "\n", print "t1: " + str(mint1) + "," print "t2: " + str(mint2) + "," print "tmax: " + str(Utilities.maxTimeResid) + "," print "}" print "\n"
def AnalysisFunction(index): minlap = 10**9 minFraction = 1000.0 for tailFraction in Utilities.TailFracs: Histogram1 = Utilities.ProduceRatioHistogram(Utilities.Particles[0] + ".root", tailFraction) Histogram2 = Utilities.ProduceRatioHistogram(Utilities.Particles[1] + ".root", tailFraction) overlap = Utilities.BhattacharyaCoeff(Histogram1, Histogram2) del Histogram1 del Histogram2 if overlap < minlap: minlap = overlap minFraction = tailFraction print "\n" print "Please place the text below into the database file: CLASSIFIER_ALPHA_UNSEEDED.ratdb located in rat/data, replacing any existing entry with the same index." print "\n" print "{" print "type: \"CLASSIFIER_ALPHA_UNSEEDED\"," print "version: 1," print "index: \"" + index + "\"," print "run_range: [0, 0]," print "pass : 0," print "production: false," print "timestamp: \"\"," print "comment: \"\"," print "\n", print "tail_fraction: " + str(minFraction) + "," print "pmt_Min: " + str(Utilities.minTimeResid) + "," print "pmt_Max: " + str(Utilities.maxTimeResid) + "," print "}" print "\n"