# 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