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"
