# aggregate = [[1], [2, 3, 4, 19, 21, 22], [5], [6, 7, 8], [9], [10, 11, 12], [14], [15], [20], [24, 25, 26], [29, 30], [32, 33, 34], [35], [36, 37, 38], [39], [40], [41], [42], [43], [44]]
## corr > .7, removing low points (13,16,17,18,23,27,28,29,30,31,42,43): 17 agg regions. agreement 91%. around 1.5s per point on commodore.
# aggregate = [[1], [2, 3, 4, 19, 21, 22], [5], [6, 7, 8], [9], [10, 11, 12], [14], [15], [20], [24, 25, 26], [32, 33, 34], [35], [36, 37, 38], [39], [40], [41], [44]]
## corr > .7, removing low points (13,16,17,18,23,24,25,26,27,28,29,30,31,42,43): 16 agg regions. agreement 91%. around 1.4s per point on commodore.
#aggregate = [[1], [2, 3, 4, 19, 21, 22], [5], [6, 7, 8], [9], [10, 11, 12], [14], [15], [20], [32, 33, 34], [35], [36, 37, 38], [39], [40], [41], [44]]
## corr > .7, moving low points (13,16,17,18,23,24,25,26,27,28,29,30,31,42,43) to separate region: 17 agg regions. agreement 94%. around 1.6s per point on commodore.
aggregate = [[1], [2, 3, 4, 19, 21, 22], [5], [6, 7, 8], [9], [10, 11, 12],
[14], [15], [20], [32, 33, 34], [35], [36, 37, 38], [39], [40],
[41], [44],
[13, 16, 17, 18, 23, 24, 25, 26, 27, 28, 29, 30, 31, 42, 43]]
info.createCovarianceMatrix ( "orig/CMS-PAS-SUS-16-052_Figure-aux_001.root", \
"c1/CovarianceSRs", addOrder=False, max_datasets = max_datasets,
aggregate = aggregate )
## fixme for now i use the root files of CMS16033, as CMS16050 currently
## only publishes the numbers in a table, not in a digitized format.
#creator = DatasetCreator ( "orig/CMS-SUS-16-033_Figure_009.root:DataObs", \
# "orig/PostFitHistograms.root:PostFitTotal", readDatasetNames=False )
creator = DatasetsFromLatex("orig/tables.tex",
max_datasets=max_datasets,
c_obs=6,
c_bg=5,
ds_name="#0",
aggregate=aggregate)
dsnames = [] ## collect the dataset names
for agg in creator.aggregate:
# aggregation = [[1, 6, 10, 11, 22, 38], [2], [3], [4], [5], [7, 12, 13, 16, 17, 18, 21, 28, 29, 30, 32, 33, 34, 36, 37, 58, 68, 69, 70, 82], [8], [9], [14], [15], [19], [20], [23], [24], [25], [26], [27], [31], [35], [39], [40], [41], [42], [43], [44], [45], [46], [47], [48], [49], [50], [51], [52], [53], [54], [55], [56], [57], [59], [60], [61], [62], [63], [64], [65], [66], [67], [71], [72], [73], [74], [75], [76], [77], [78], [79], [80], [81], [83], [84]]
## corr>.4: remove 0-points (23,25,41,77), 56 aggregate regions, agreements 98/97/97. Takes around 10s per point on higgs. overexcludes for one case (T2tt?)
aggregation = [[1, 6, 10, 11, 22, 38], [2], [3], [4], [5], [7, 12, 13, 16, 17, 18, 21, 28, 29, 30, 32, 33, 34, 36, 37, 58, 68, 69, 70, 82], [8], [9], [14], [15], [19], [20], [24], [26], [27], [31], [35], [39], [40], [42], [43], [44], [45], [46], [47], [48], [49], [50], [51], [52], [53], [54], [55], [56], [57], [59], [60], [61], [62], [63], [64], [65], [66], [67], [71], [72], [73], [74], [75], [76], [78], [79], [80], [81], [83], [84]]
## corr>.4: remove low-points (5,8,23,24,25,35,41,45,48,52,59,62,66,67,74,77,78,80,84), 41 aggregate regions, agreements 91/97/94. Takes around 5.5s per point on commodore.
# aggregation = [[1, 6, 10, 11, 22, 38], [2], [3], [4], [7, 12, 13, 16, 17, 18, 21, 28, 29, 30, 32, 33, 34, 36, 37, 58, 68, 69, 70, 82], [9], [14], [15], [19], [20], [26], [27], [31], [39], [40], [42], [43], [44], [46], [47], [49], [50], [51], [53], [54], [55], [56], [57], [60], [61], [63], [64], [65], [71], [72], [73], [75], [76], [79], [81], [83]]
## corr>.4: move low-points (5,8,23,24,25,35,41,45,48,52,59,62,66,67,74,77,78,80,84) to its own region, 42 aggregate regions, agreements 97/90/93. Takes around 7s per point on higgs.
# aggregation = [[1, 6, 10, 11, 22, 38], [2], [3], [4], [7, 12, 13, 16, 17, 18, 21, 28, 29, 30, 32, 33, 34, 36, 37, 58, 68, 69, 70, 82], [9], [14], [15], [19], [20], [26], [27], [31], [39], [40], [42], [43], [44], [46], [47], [49], [50], [51], [53], [54], [55], [56], [57], [60], [61], [63], [64], [65], [71], [72], [73], [75], [76], [79], [81], [83], [5,8,23,24,25,35,41,45,48,52,59,62,66,67,74,77,78,80,84] ]
## corr>.3:
# aggregation = [[1, 2, 6, 10, 11, 15, 22, 23, 38, 53, 54, 59, 63], [3, 4, 9, 51], [5, 7, 12, 13, 16, 17, 18, 21, 24, 28, 29, 30, 32, 33, 34, 35, 36, 37, 57, 58, 68, 69, 70, 82], [14], [19], [20], [26], [27], [31], [39], [40], [42], [43], [44], [46], [47], [49], [50], [55], [56], [60], [61], [64], [65], [71], [72], [73], [75], [76], [79], [81], [83]]
## low pointers [1, 2, 5, 6, 7, 8, 11, 12, 13, 23, 24, 25, 29, 33, 34, 35, 36, 37, 38, 41, 45, 48, 52, 59, 62, 66, 67, 69, 74, 77, 78, 80, 84]
info.createCovarianceMatrix ( "orig/CMS-SUS-16-050_Figure-aux_001.root", \
"total_covar", addOrder=True, max_datasets = max_datasets, aggregate=aggregation )
creator = DatasetsFromLatex ( "orig/tables.tex", max_datasets = max_datasets,
ds_name = "t#1b#2MT2#3MET#4", aggregate=aggregation )
for ctr,dataset in enumerate(creator):
#+++++++ next txName block ++++++++++++++
T2tt = dataset.addTxName('T2tt')
T2tt.checked =''
T2tt.constraint = "[[['t']],[['t']]]"
T2tt.conditionDescription = None
T2tt.condition =None
T2tt.source ='CMS'
#+++++++ next mass plane block ++++++++++++++
T2tt_1 = T2tt.addMassPlane([[x,y]]*2)
objNames = 'acc_bin%d' % ctr