-
Notifications
You must be signed in to change notification settings - Fork 184
/
Copy path_morph_functions_8cc_source.html
818 lines (816 loc) · 125 KB
/
_morph_functions_8cc_source.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
<!-- HTML header for doxygen 1.8.6-->
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.9.1"/>
<title>CombineHarvester: CombinePdfs/src/MorphFunctions.cc Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="resize.js"></script>
<script type="text/javascript" src="navtreedata.js"></script>
<script type="text/javascript" src="navtree.js"></script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/searchdata.js"></script>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
$(document).ready(function() { init_search(); });
/* @license-end */
</script>
<script type="text/x-mathjax-config">
MathJax.Hub.Config({
extensions: ["tex2jax.js"],
jax: ["input/TeX","output/HTML-CSS"],
});
</script>
<script type="text/javascript" async="async" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
<link href="customdoxygen.css" rel="stylesheet" type="text/css"/>
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 70px;">
<td style="padding-left: 0.5em;">
<div id="projectname">CombineHarvester
</div>
</td>
<td> <div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<img id="MSearchSelect" src="search/mag_sel.svg"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
alt=""/>
<input type="text" id="MSearchField" value="Search" accesskey="S"
onfocus="searchBox.OnSearchFieldFocus(true)"
onblur="searchBox.OnSearchFieldFocus(false)"
onkeyup="searchBox.OnSearchFieldChange(event)"/>
</span><span class="right">
<a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.svg" alt=""/></a>
</span>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.9.1 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
var searchBox = new SearchBox("searchBox", "search",false,'Search','.html');
/* @license-end */
</script>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
<div id="nav-tree">
<div id="nav-tree-contents">
<div id="nav-sync" class="sync"></div>
</div>
</div>
<div id="splitbar" style="-moz-user-select:none;"
class="ui-resizable-handle">
</div>
</div>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
$(document).ready(function(){initNavTree('_morph_functions_8cc_source.html',''); initResizable(); });
/* @license-end */
</script>
<div id="doc-content">
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div class="header">
<div class="headertitle">
<div class="title">MorphFunctions.cc</div> </div>
</div><!--header-->
<div class="contents">
<a href="_morph_functions_8cc.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="preprocessor">#include "<a class="code" href="_morph_functions_8h.html">CombineHarvester/CombinePdfs/interface/MorphFunctions.h</a>"</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="preprocessor">#include <iostream></span></div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="preprocessor">#include <set></span></div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="preprocessor">#include <vector></span></div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="preprocessor">#include <string></span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="preprocessor">#include "boost/lexical_cast.hpp"</span></div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="preprocessor">#include "boost/format.hpp"</span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="preprocessor">#include "boost/multi_array.hpp"</span></div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="preprocessor">#include "TGraphErrors.h"</span></div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="preprocessor">#include "RooFitResult.h"</span></div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="preprocessor">#include "RooRealVar.h"</span></div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="preprocessor">#include "RooDataHist.h"</span></div>
<div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="preprocessor">#include "RooProduct.h"</span></div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="preprocessor">#include "RooConstVar.h"</span></div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="preprocessor">#include "<a class="code" href="_logging_8h.html">CombineHarvester/CombineTools/interface/Logging.h</a>"</span></div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span>  </div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="keyword">namespace </span><a class="code" href="namespacech.html">ch</a> {</div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span>  </div>
<div class="line"><a name="l00020"></a><span class="lineno"><a class="line" href="namespacech.html#a011bcc8eb52bfa118fd44ade245e9e45"> 20</a></span> std::string <a class="code" href="namespacech.html#a011bcc8eb52bfa118fd44ade245e9e45">BuildRooMorphing</a>(RooWorkspace& ws, <a class="code" href="classch_1_1_combine_harvester.html">CombineHarvester</a>& cb,</div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span>  std::string <span class="keyword">const</span>& bin, std::string <span class="keyword">const</span>& <a class="code" href="structch_1_1syst_1_1process.html">process</a>,</div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span>  RooAbsReal& mass_var, std::string norm_postfix,</div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span>  <span class="keywordtype">bool</span> allow_morph, <span class="keywordtype">bool</span> verbose, <span class="keywordtype">bool</span> force_template_limit, TFile * file) {</div>
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span>  <span class="comment">// To keep the code concise we'll make some using-declarations here</span></div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span>  <span class="keyword">using</span> std::set;</div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span>  <span class="keyword">using</span> std::vector;</div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span>  <span class="keyword">using</span> std::string;</div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span>  <span class="keyword">using</span> boost::lexical_cast;</div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span>  <span class="keyword">using</span> boost::multi_array;</div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span>  <span class="keyword">using</span> boost::extents;</div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span>  </div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span>  <span class="comment">// RooFit can be rather noisy on stdout, unless verbose = true we'll kill</span></div>
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span>  <span class="comment">// everything that's not an error, and restore the user's MsgLevel at the end</span></div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span>  RooFit::MsgLevel backup_msg_level =</div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span>  RooMsgService::instance().globalKillBelow();</div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span>  <span class="keywordflow">if</span> (!verbose) RooMsgService::instance().setGlobalKillBelow(RooFit::WARNING);</div>
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  </div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span>  <span class="comment">// Everything here works on the assumption that the bin,process,mass</span></div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span>  <span class="comment">// combinations we find in this CH instance refer to single processes. So bin</span></div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span>  <span class="comment">// should uniquely label an event category, and mass should be defined and</span></div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span>  <span class="comment">// convertible from string to float for every process entry.</span></div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  <span class="keywordflow">if</span> (verbose) {</div>
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span>  std::cout << <span class="stringliteral">">> Bin: "</span> << bin << <span class="stringliteral">" Process: "</span> << <a class="code" href="structch_1_1syst_1_1process.html">process</a> << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  }</div>
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  TString key = bin + <span class="stringliteral">"_"</span> + <a class="code" href="structch_1_1syst_1_1process.html">process</a>;</div>
<div class="line"><a name="l00047"></a><span class="lineno"> 47</span>  </div>
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span>  <span class="comment">// For the sake of efficiency, take a filtered shallow-copy of the input CH</span></div>
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  <span class="comment">// instance that just contains the bin,process combination we've been asked to</span></div>
<div class="line"><a name="l00050"></a><span class="lineno"> 50</span>  <span class="comment">// turn into a RooMorphingPdf</span></div>
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span>  <a class="code" href="classch_1_1_combine_harvester.html">CombineHarvester</a> cb_bp = cb.<a class="code" href="classch_1_1_combine_harvester.html#ade0f7c161d078a189404a274b93336d2">cp</a>().<a class="code" href="classch_1_1_combine_harvester.html#a83dd50569a1f1a4b5a46f2ad2b5c1cdb">bin</a>({bin}).<a class="code" href="structch_1_1syst_1_1process.html">process</a>({<a class="code" href="structch_1_1syst_1_1process.html">process</a>});</div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span>  </div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span>  <span class="comment">// Get a vector of the mass values</span></div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  vector<string> m_str_vec = <a class="code" href="namespacech.html#a5e23234ab65cfa48c99b83d94e291125">Set2Vec</a>(cb_bp.<a class="code" href="classch_1_1_combine_harvester.html#a7a0c6913b055e11798269034ae324298">SetFromProcs</a>(</div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  std::mem_fn(&<a class="code" href="classch_1_1_object.html#ae2f9e80af391e105a6b6bc27107e3c9e">ch::Process::mass</a>)));</div>
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span>  </div>
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span>  <span class="comment">// We now sort these based on numerical value, hence the lexical_cast to double</span></div>
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  <span class="comment">// At the moment this also serves as our check that all mass values are</span></div>
<div class="line"><a name="l00060"></a><span class="lineno"> 60</span>  <span class="comment">// float-convertible, as an exception will be thrown by lexical_cast if not.</span></div>
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  std::sort(m_str_vec.begin(), m_str_vec.end(),</div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  [](<span class="keywordtype">string</span> <span class="keyword">const</span>& s1, <span class="keywordtype">string</span> <span class="keyword">const</span>& s2) {</div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  return lexical_cast<double>(s1) < lexical_cast<double>(s2);</div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  });</div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  </div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  <span class="comment">// Convert the sorted vector of mass strings to an actual vector of doubles</span></div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span>  vector<double> m_vec;</div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span>  <span class="keywordflow">for</span> (<span class="keyword">auto</span> <span class="keyword">const</span>& s : m_str_vec) {</div>
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  <span class="keywordflow">if</span> (verbose) std::cout << <span class="stringliteral">">>>> Mass point: "</span> << s << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00070"></a><span class="lineno"> 70</span>  m_vec.push_back(lexical_cast<double>(s));</div>
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span>  }</div>
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  <span class="comment">// So, we have m mass points to consider</span></div>
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  <span class="keywordtype">unsigned</span> m = m_vec.size();</div>
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  </div>
<div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  <span class="comment">// ss = "shape systematic"</span></div>
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span>  <span class="comment">// Make a list of the names of shape systematics affecting this process</span></div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  vector<string> ss_vec =</div>
<div class="line"><a name="l00079"></a><span class="lineno"> 79</span>  <a class="code" href="namespacech.html#a5e23234ab65cfa48c99b83d94e291125">Set2Vec</a>(cb_bp.<a class="code" href="classch_1_1_combine_harvester.html#ade0f7c161d078a189404a274b93336d2">cp</a>().<a class="code" href="classch_1_1_combine_harvester.html#af13ff7dd8c281cd9e9cfd6b11fc05b95">syst_type</a>({<span class="stringliteral">"shape"</span>}).<a class="code" href="classch_1_1_combine_harvester.html#ac3ef7424cb86d12fd294d9c6de98d533">syst_name_set</a>());</div>
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  <span class="comment">// Now check if all shape systematics are present for all mass points</span></div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  <span class="keywordflow">for</span> (<span class="keyword">auto</span> <span class="keyword">const</span>& s : m_str_vec) {</div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  <span class="keywordflow">if</span> (cb_bp.<a class="code" href="classch_1_1_combine_harvester.html#ade0f7c161d078a189404a274b93336d2">cp</a>().<a class="code" href="classch_1_1_combine_harvester.html#af13ff7dd8c281cd9e9cfd6b11fc05b95">syst_type</a>({<span class="stringliteral">"shape"</span>}).<a class="code" href="classch_1_1_combine_harvester.html#a25c231ed8b1614c5d50f192210374474">mass</a>({s}).<a class="code" href="classch_1_1_combine_harvester.html#ac3ef7424cb86d12fd294d9c6de98d533">syst_name_set</a>().size() !=</div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  ss_vec.size()) {</div>
<div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  <span class="keywordflow">throw</span> std::runtime_error(<a class="code" href="_logging_8h.html#ab140ae8826e1f99c5aa0e920a9c13112">FNERROR</a>(</div>
<div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  <span class="stringliteral">"Some mass points do not have the full set of shape systematics, "</span></div>
<div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  <span class="stringliteral">"this is currently unsupported"</span>));</div>
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  }</div>
<div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  }</div>
<div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  </div>
<div class="line"><a name="l00090"></a><span class="lineno"> 90</span>  <span class="keywordtype">unsigned</span> ss = ss_vec.size(); <span class="comment">// number of shape systematics</span></div>
<div class="line"><a name="l00091"></a><span class="lineno"> 91</span>  </div>
<div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  <span class="comment">// ls = "lnN systematic"</span></div>
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  <span class="comment">// Make a list of the regular lnN normalisation systematics affecting this</span></div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  <span class="comment">// process</span></div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  vector<string> ls_vec =</div>
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  <a class="code" href="namespacech.html#a5e23234ab65cfa48c99b83d94e291125">Set2Vec</a>(cb_bp.<a class="code" href="classch_1_1_combine_harvester.html#ade0f7c161d078a189404a274b93336d2">cp</a>().<a class="code" href="classch_1_1_combine_harvester.html#af13ff7dd8c281cd9e9cfd6b11fc05b95">syst_type</a>({<span class="stringliteral">"lnN"</span>}).<a class="code" href="classch_1_1_combine_harvester.html#ac3ef7424cb86d12fd294d9c6de98d533">syst_name_set</a>());</div>
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  <span class="keywordtype">unsigned</span> ls = ls_vec.size(); <span class="comment">// number of lnN systematics</span></div>
<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  </div>
<div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  <span class="comment">// Create a bunch of empty arrays to store the information we need in a more</span></div>
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  <span class="comment">// convenient format. We use the boost multi_array class because it's a nice</span></div>
<div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  <span class="comment">// multi-dimensional array implementation, and safer to use than standard</span></div>
<div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  <span class="comment">// C-arrays</span></div>
<div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  </div>
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  <span class="comment">// Store pointers to each ch::Process (one per mass point) in the CH instance</span></div>
<div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  multi_array<ch::Process *, 1> pr_arr(extents[m]);</div>
<div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  <span class="comment">// Store pointers pointers to each ch::Systematic for each mass point (hence</span></div>
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  <span class="comment">// an ss * m 2D array)</span></div>
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  multi_array<ch::Systematic *, 2> ss_arr(extents[ss][m]);</div>
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span>  <span class="comment">// With shape systematics we have to support cases where the value in the</span></div>
<div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  <span class="comment">// datacard is != 1.0, i.e. we are scaling the parameter that goes into the</span></div>
<div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  <span class="comment">// Gaussian constraint PDF - we have to pass this factor on when we build the</span></div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  <span class="comment">// normal vertical-interpolation PDF for each mass point. We will store these</span></div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  <span class="comment">// scale factors in this array, one per shape systematic.</span></div>
<div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  multi_array<double, 1> ss_scale_arr(extents[ss]);</div>
<div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  <span class="comment">// Really just for book-keeping, we'll set this flag to true when the shape</span></div>
<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  <span class="comment">// systematic scale factor != 1</span></div>
<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  multi_array<bool, 1> ss_must_scale_arr(extents[ss]);</div>
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  <span class="comment">// Similar to the ss_arr above, store the lnN ch::Systematic objects. Note</span></div>
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  <span class="comment">// implicit in all this is the assumption that the processes at each mass</span></div>
<div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  <span class="comment">// point have exactly the same list of systematic uncertainties.</span></div>
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  multi_array<ch::Systematic *, 2> ls_arr(extents[ls][m]);</div>
<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  </div>
<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  <span class="comment">// This array holds pointers to the RooRealVar objects that will become our</span></div>
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span>  <span class="comment">// shape nuisance parameters, e.g. "CMS_scale_t_mutau_8TeV"</span></div>
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  multi_array<std::shared_ptr<RooRealVar>, 1> ss_scale_var_arr(extents[ss]);</div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  <span class="comment">// And this array holds the constant scale factors that we'll build from</span></div>
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  <span class="comment">// ss_scale_arr above</span></div>
<div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  multi_array<std::shared_ptr<RooConstVar>, 1> ss_scale_fac_arr(extents[ss]);</div>
<div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  <span class="comment">// Finally this array will contain the scale_var * scale_fac product where we</span></div>
<div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  <span class="comment">// need to provide a scaled value of the nuisance parameter instead of the</span></div>
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  <span class="comment">// parameter itself in building the vertical-interp. PDF</span></div>
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  multi_array<std::shared_ptr<RooProduct>, 1> ss_scale_prod_arr(extents[ss]);</div>
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  </div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  </div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  <span class="comment">// Now let's fill some of these arrays...</span></div>
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  <span class="comment">// The ch::Process pointers</span></div>
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  cb_bp.<a class="code" href="classch_1_1_combine_harvester.html#ade0f7c161d078a189404a274b93336d2">cp</a>().<a class="code" href="classch_1_1_combine_harvester.html#a25c231ed8b1614c5d50f192210374474">mass</a>({m_str_vec[mi]}).ForEachProc([&](<a class="code" href="classch_1_1_process.html">ch::Process</a> *p) {</div>
<div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  pr_arr[mi] = p;</div>
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  });</div>
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> ssi = 0; ssi < ss; ++ssi) {</div>
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  <span class="comment">// The ch::Systematic pointers for shape systematics</span></div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  cb_bp.<a class="code" href="classch_1_1_combine_harvester.html#ade0f7c161d078a189404a274b93336d2">cp</a>().<a class="code" href="classch_1_1_combine_harvester.html#a25c231ed8b1614c5d50f192210374474">mass</a>({m_str_vec[mi]}).syst_name({ss_vec[ssi]})</div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  .ForEachSyst([&](<a class="code" href="classch_1_1_systematic.html">ch::Systematic</a> *n) {</div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  ss_arr[ssi][mi] = n;</div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  });</div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  }</div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> lsi = 0; lsi < ls; ++lsi) {</div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  <span class="comment">// The ch::Systematic pointers for lnN systematics</span></div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  cb_bp.<a class="code" href="classch_1_1_combine_harvester.html#ade0f7c161d078a189404a274b93336d2">cp</a>().<a class="code" href="classch_1_1_combine_harvester.html#a25c231ed8b1614c5d50f192210374474">mass</a>({m_str_vec[mi]}).syst_name({ls_vec[lsi]})</div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  .ForEachSyst([&](<a class="code" href="classch_1_1_systematic.html">ch::Systematic</a> *n) {</div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span>  ls_arr[lsi][mi] = n;</div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  });</div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  }</div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  }</div>
<div class="line"><a name="l00156"></a><span class="lineno"> 156</span>  </div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  <span class="comment">// We need to build a RooArgList of the vertical morphing parameters for the</span></div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  <span class="comment">// vertical-interpolation pdf - this will be the same for each mass point so</span></div>
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  <span class="comment">// we only build it once</span></div>
<div class="line"><a name="l00161"></a><span class="lineno"> 161</span>  RooArgList ss_list;</div>
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> ssi = 0; ssi < ss; ++ssi) {</div>
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  <span class="comment">// Make the nuisance parameter var. We use shared_ptrs here as a convenience</span></div>
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  <span class="comment">// because they will take care of automatically cleaning up the memory at</span></div>
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  <span class="comment">// the end</span></div>
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  ss_scale_var_arr[ssi] =</div>
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  std::make_shared<RooRealVar>(ss_vec[ssi].c_str(), <span class="stringliteral">""</span>, 0);</div>
<div class="line"><a name="l00168"></a><span class="lineno"> 168</span>  <span class="comment">// We'll make a quick check that the scale factor for this systematic is the</span></div>
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span>  <span class="comment">// same for all mass points. We could do a separate scaling at each mass</span></div>
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span>  <span class="comment">// point but this would create a lot of complications</span></div>
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span>  set<double> scales;</div>
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span>  <span class="comment">// Insert the scale from each mass point into the set, if it has a size</span></div>
<div class="line"><a name="l00173"></a><span class="lineno"> 173</span>  <span class="comment">// larger than one at the end then we have a problem!</span></div>
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  scales.insert(ss_arr[ssi][mi]->scale());</div>
<div class="line"><a name="l00176"></a><span class="lineno"> 176</span>  }</div>
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  <span class="keywordflow">if</span> (scales.size() > 1) {</div>
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  <span class="comment">// Don't let the user proceed, we can't build the model they want</span></div>
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  std::runtime_error(<a class="code" href="_logging_8h.html#ab140ae8826e1f99c5aa0e920a9c13112">FNERROR</a>(</div>
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  <span class="stringliteral">"Shape morphing parameters that vary with mass are not allowed"</span>));</div>
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  } <span class="keywordflow">else</span> {</div>
<div class="line"><a name="l00182"></a><span class="lineno"> 182</span>  <span class="comment">// Everything ok, set the scale value in its array</span></div>
<div class="line"><a name="l00183"></a><span class="lineno"> 183</span>  ss_scale_arr[ssi] = *(scales.begin());</div>
<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  <span class="comment">// Handle the case where the scale factor is != 1</span></div>
<div class="line"><a name="l00185"></a><span class="lineno"> 185</span>  <span class="keywordflow">if</span> (std::fabs(ss_scale_arr[ssi] - 1.0) > 1E-6) {</div>
<div class="line"><a name="l00186"></a><span class="lineno"> 186</span>  ss_must_scale_arr[ssi] = <span class="keyword">true</span>;</div>
<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  <span class="comment">// Build the RooConstVar with the value of the scale factor</span></div>
<div class="line"><a name="l00188"></a><span class="lineno"> 188</span>  ss_scale_fac_arr[ssi] = std::make_shared<RooConstVar>(</div>
<div class="line"><a name="l00189"></a><span class="lineno"> 189</span>  TString::Format(<span class="stringliteral">"%g"</span>, ss_scale_arr[ssi]), <span class="stringliteral">""</span>,</div>
<div class="line"><a name="l00190"></a><span class="lineno"> 190</span>  ss_scale_arr[ssi]);</div>
<div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  <span class="comment">// Create the product of the floating nuisance parameter and the</span></div>
<div class="line"><a name="l00192"></a><span class="lineno"> 192</span>  <span class="comment">// constant scale factor</span></div>
<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  ss_scale_prod_arr[ssi] = std::make_shared<RooProduct>(</div>
<div class="line"><a name="l00194"></a><span class="lineno"> 194</span>  ss_vec[ssi] + <span class="stringliteral">"_scaled_"</span> + key, <span class="stringliteral">""</span>,</div>
<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  RooArgList(*(ss_scale_var_arr[ssi]), *(ss_scale_fac_arr[ssi])));</div>
<div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  <span class="comment">// Add this to the list</span></div>
<div class="line"><a name="l00197"></a><span class="lineno"> 197</span>  ss_list.add(*(ss_scale_prod_arr[ssi]));</div>
<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  } <span class="keywordflow">else</span> {</div>
<div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  <span class="comment">// If the scale factor is 1.0 then we just add the nuisance parameter</span></div>
<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  <span class="comment">// directly to our list</span></div>
<div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  ss_list.add(*(ss_scale_var_arr[ssi]));</div>
<div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  }</div>
<div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  }</div>
<div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  }</div>
<div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  </div>
<div class="line"><a name="l00207"></a><span class="lineno"> 207</span>  <span class="comment">// Summarise the info on the shape systematics and scale factors</span></div>
<div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  <span class="keywordflow">if</span> (verbose) {</div>
<div class="line"><a name="l00209"></a><span class="lineno"> 209</span>  std::cout << <span class="stringliteral">">> Shape systematics: "</span> << ss << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00210"></a><span class="lineno"> 210</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> ssi = 0; ssi < ss; ++ssi) {</div>
<div class="line"><a name="l00211"></a><span class="lineno"> 211</span>  std::cout << boost::format(<span class="stringliteral">"%-50s %-5i %-8.3g\n"</span>)</div>
<div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  % ss_vec[ssi] % ss_must_scale_arr[ssi] % ss_scale_arr[ssi];</div>
<div class="line"><a name="l00213"></a><span class="lineno"> 213</span>  }</div>
<div class="line"><a name="l00214"></a><span class="lineno"> 214</span>  }</div>
<div class="line"><a name="l00215"></a><span class="lineno"> 215</span>  </div>
<div class="line"><a name="l00216"></a><span class="lineno"> 216</span>  </div>
<div class="line"><a name="l00217"></a><span class="lineno"> 217</span>  <span class="comment">// lms = "lnN morphing systematic"</span></div>
<div class="line"><a name="l00218"></a><span class="lineno"> 218</span>  <span class="comment">// Now we have some work to do with the lnN systematics. We can consider two cases:</span></div>
<div class="line"><a name="l00219"></a><span class="lineno"> 219</span>  <span class="comment">// a) The uncertainty is the same for each mass point => we can leave it in</span></div>
<div class="line"><a name="l00220"></a><span class="lineno"> 220</span>  <span class="comment">// the datacard as is and let text2workspace do its normal thing</span></div>
<div class="line"><a name="l00221"></a><span class="lineno"> 221</span>  <span class="comment">// b) The uncertainty varies between mass points => we can't capture this</span></div>
<div class="line"><a name="l00222"></a><span class="lineno"> 222</span>  <span class="comment">// information in the text datacard in the usual way, so we'll build a RooFit</span></div>
<div class="line"><a name="l00223"></a><span class="lineno"> 223</span>  <span class="comment">// object that effectively makes the lnN uncertainty a function of the mass</span></div>
<div class="line"><a name="l00224"></a><span class="lineno"> 224</span>  <span class="comment">// variable</span></div>
<div class="line"><a name="l00225"></a><span class="lineno"> 225</span>  <span class="comment">// We'll use "lms" to refer to case b), which we'll try to figure out now...</span></div>
<div class="line"><a name="l00226"></a><span class="lineno"> 226</span>  vector<string> lms_vec; <span class="comment">// vec of systematic names</span></div>
<div class="line"><a name="l00227"></a><span class="lineno"> 227</span>  set<string> lms_set; <span class="comment">// set of systematic names</span></div>
<div class="line"><a name="l00228"></a><span class="lineno"> 228</span>  <span class="comment">// index positions in our full ls_arr array for the lms systematics</span></div>
<div class="line"><a name="l00229"></a><span class="lineno"> 229</span>  vector<unsigned > lms_vec_idx; </div>
<div class="line"><a name="l00230"></a><span class="lineno"> 230</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> lsi = 0; lsi < ls; ++lsi) {</div>
<div class="line"><a name="l00231"></a><span class="lineno"> 231</span>  <span class="comment">// Extra complication is that the user might have been evil and mixed</span></div>
<div class="line"><a name="l00232"></a><span class="lineno"> 232</span>  <span class="comment">// symmetric and asymmetric lnN values, we'll try and detect changes in</span></div>
<div class="line"><a name="l00233"></a><span class="lineno"> 233</span>  <span class="comment">// either</span></div>
<div class="line"><a name="l00234"></a><span class="lineno"> 234</span>  set<double> k_hi;</div>
<div class="line"><a name="l00235"></a><span class="lineno"> 235</span>  set<double> k_lo;</div>
<div class="line"><a name="l00236"></a><span class="lineno"> 236</span>  <span class="comment">// Go through each mass point for this systematic and add the uncertainty</span></div>
<div class="line"><a name="l00237"></a><span class="lineno"> 237</span>  <span class="comment">// values (so-called "kappa" values)</span></div>
<div class="line"><a name="l00238"></a><span class="lineno"> 238</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00239"></a><span class="lineno"> 239</span>  <a class="code" href="classch_1_1_systematic.html">Systematic</a> *n = ls_arr[lsi][mi];</div>
<div class="line"><a name="l00240"></a><span class="lineno"> 240</span>  k_hi.insert(n-><a class="code" href="classch_1_1_systematic.html#a1b40f502a9e5343cc1947bf452cd80d5">value_u</a>());</div>
<div class="line"><a name="l00241"></a><span class="lineno"> 241</span>  <span class="keywordflow">if</span> (n-><a class="code" href="classch_1_1_systematic.html#a5764d6db4954b8aac554776c149d3201">asymm</a>()) {</div>
<div class="line"><a name="l00242"></a><span class="lineno"> 242</span>  k_lo.insert(n-><a class="code" href="classch_1_1_systematic.html#ac0438bd51a527d0b2c1a70dd303ff4c3">value_d</a>());</div>
<div class="line"><a name="l00243"></a><span class="lineno"> 243</span>  }</div>
<div class="line"><a name="l00244"></a><span class="lineno"> 244</span>  }</div>
<div class="line"><a name="l00245"></a><span class="lineno"> 245</span>  <span class="comment">// If either of these sets has more than one entry then this is a lms case</span></div>
<div class="line"><a name="l00246"></a><span class="lineno"> 246</span>  <span class="keywordflow">if</span> (k_hi.size() > 1 || k_lo.size() > 1) {</div>
<div class="line"><a name="l00247"></a><span class="lineno"> 247</span>  lms_vec.push_back(ls_vec[lsi]);</div>
<div class="line"><a name="l00248"></a><span class="lineno"> 248</span>  lms_set.insert(ls_vec[lsi]);</div>
<div class="line"><a name="l00249"></a><span class="lineno"> 249</span>  lms_vec_idx.push_back(lsi);</div>
<div class="line"><a name="l00250"></a><span class="lineno"> 250</span>  }</div>
<div class="line"><a name="l00251"></a><span class="lineno"> 251</span>  }</div>
<div class="line"><a name="l00252"></a><span class="lineno"> 252</span>  <span class="keywordtype">unsigned</span> lms = lms_vec.size();</div>
<div class="line"><a name="l00253"></a><span class="lineno"> 253</span>  <span class="comment">// New array for the pointers to the lms Systematic objects</span></div>
<div class="line"><a name="l00254"></a><span class="lineno"> 254</span>  multi_array<ch::Systematic *, 2> lms_arr(extents[lms][m]);</div>
<div class="line"><a name="l00255"></a><span class="lineno"> 255</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> lmsi = 0; lmsi < lms; ++lmsi) {</div>
<div class="line"><a name="l00256"></a><span class="lineno"> 256</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00257"></a><span class="lineno"> 257</span>  lms_arr[lmsi][mi] = ls_arr[lms_vec_idx[lmsi]][mi];</div>
<div class="line"><a name="l00258"></a><span class="lineno"> 258</span>  }</div>
<div class="line"><a name="l00259"></a><span class="lineno"> 259</span>  }</div>
<div class="line"><a name="l00260"></a><span class="lineno"> 260</span>  </div>
<div class="line"><a name="l00261"></a><span class="lineno"> 261</span>  <span class="comment">// We will need to create the nuisance parameters for these now</span></div>
<div class="line"><a name="l00262"></a><span class="lineno"> 262</span>  multi_array<std::shared_ptr<RooRealVar>, 1> lms_var_arr(extents[lms]);</div>
<div class="line"><a name="l00263"></a><span class="lineno"> 263</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> lmsi = 0; lmsi < lms; ++lmsi) {</div>
<div class="line"><a name="l00264"></a><span class="lineno"> 264</span>  lms_var_arr[lmsi] =</div>
<div class="line"><a name="l00265"></a><span class="lineno"> 265</span>  std::make_shared<RooRealVar>(lms_vec[lmsi].c_str(), <span class="stringliteral">""</span>, 0);</div>
<div class="line"><a name="l00266"></a><span class="lineno"> 266</span>  }</div>
<div class="line"><a name="l00267"></a><span class="lineno"> 267</span>  </div>
<div class="line"><a name="l00268"></a><span class="lineno"> 268</span>  <span class="comment">// Give a summary of the lms systematics to the user</span></div>
<div class="line"><a name="l00269"></a><span class="lineno"> 269</span>  <span class="keywordflow">if</span> (verbose) {</div>
<div class="line"><a name="l00270"></a><span class="lineno"> 270</span>  std::cout << <span class="stringliteral">">> lnN morphing systematics: "</span> << lms << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00271"></a><span class="lineno"> 271</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> lmsi = 0; lmsi < lms; ++lmsi) {</div>
<div class="line"><a name="l00272"></a><span class="lineno"> 272</span>  std::cout << <span class="stringliteral">">>>> "</span> << lms_vec[lmsi] << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00273"></a><span class="lineno"> 273</span>  }</div>
<div class="line"><a name="l00274"></a><span class="lineno"> 274</span>  }</div>
<div class="line"><a name="l00275"></a><span class="lineno"> 275</span>  </div>
<div class="line"><a name="l00276"></a><span class="lineno"> 276</span>  <span class="comment">// Now we move into a phase of building all the objects we need:</span></div>
<div class="line"><a name="l00277"></a><span class="lineno"> 277</span>  </div>
<div class="line"><a name="l00278"></a><span class="lineno"> 278</span>  <span class="comment">// 2D array of all input histograms, size is (mass points * (nominal +</span></div>
<div class="line"><a name="l00279"></a><span class="lineno"> 279</span>  <span class="comment">// 2*shape-systs)). The factor of 2 needed for the Up and Down shapes</span></div>
<div class="line"><a name="l00280"></a><span class="lineno"> 280</span>  multi_array<std::shared_ptr<TH1F>, 2> hist_arr(extents[m][1+ss*2]);</div>
<div class="line"><a name="l00281"></a><span class="lineno"> 281</span>  <span class="comment">// We also need the array of process yields vs mass, because this will have to</span></div>
<div class="line"><a name="l00282"></a><span class="lineno"> 282</span>  <span class="comment">// be interpolated too</span></div>
<div class="line"><a name="l00283"></a><span class="lineno"> 283</span>  multi_array<double, 1> rate_arr(extents[m]);</div>
<div class="line"><a name="l00284"></a><span class="lineno"> 284</span>  <span class="comment">// Also store the yield uncertainty - we don't actually need this for the signal</span></div>
<div class="line"><a name="l00285"></a><span class="lineno"> 285</span>  <span class="comment">// model, but will include it in the debug TGraph</span></div>
<div class="line"><a name="l00286"></a><span class="lineno"> 286</span>  multi_array<double, 1> rate_unc_arr(extents[m]);</div>
<div class="line"><a name="l00287"></a><span class="lineno"> 287</span>  <span class="comment">// The vertical-interpolation PDF needs the TH1 inputs in the format of a TList</span></div>
<div class="line"><a name="l00288"></a><span class="lineno"> 288</span>  multi_array<std::shared_ptr<TList>, 1> list_arr(extents[m]);</div>
<div class="line"><a name="l00289"></a><span class="lineno"> 289</span>  <span class="comment">// Combine always treats the normalisation part of shape systematics as</span></div>
<div class="line"><a name="l00290"></a><span class="lineno"> 290</span>  <span class="comment">// distinct from the actual shape morphing. Essentially the norm part is</span></div>
<div class="line"><a name="l00291"></a><span class="lineno"> 291</span>  <span class="comment">// treated as an asymmetric lnN. We have to make the kappa_hi and kappa_lo a</span></div>
<div class="line"><a name="l00292"></a><span class="lineno"> 292</span>  <span class="comment">// function of the mass parameter too, so we need two more arrays in (ss * m)</span></div>
<div class="line"><a name="l00293"></a><span class="lineno"> 293</span>  multi_array<double, 2> ss_k_hi_arr(extents[ss][m]);</div>
<div class="line"><a name="l00294"></a><span class="lineno"> 294</span>  multi_array<double, 2> ss_k_lo_arr(extents[ss][m]);</div>
<div class="line"><a name="l00295"></a><span class="lineno"> 295</span>  <span class="comment">// For each shape systematic we will build a RooSpline1D, configured to</span></div>
<div class="line"><a name="l00296"></a><span class="lineno"> 296</span>  <span class="comment">// interpolate linearly between the kappa values</span></div>
<div class="line"><a name="l00297"></a><span class="lineno"> 297</span>  multi_array<std::shared_ptr<RooSpline1D>, 1> ss_spl_hi_arr(extents[ss]);</div>
<div class="line"><a name="l00298"></a><span class="lineno"> 298</span>  multi_array<std::shared_ptr<RooSpline1D>, 1> ss_spl_lo_arr(extents[ss]);</div>
<div class="line"><a name="l00299"></a><span class="lineno"> 299</span>  <span class="comment">// To define the actually process scaling as a function of these kappa values,</span></div>
<div class="line"><a name="l00300"></a><span class="lineno"> 300</span>  <span class="comment">// we need an AsymPow object per mass point</span></div>
<div class="line"><a name="l00301"></a><span class="lineno"> 301</span>  multi_array<std::shared_ptr<AsymPow>, 1> ss_asy_arr(extents[ss]);</div>
<div class="line"><a name="l00302"></a><span class="lineno"> 302</span>  </div>
<div class="line"><a name="l00303"></a><span class="lineno"> 303</span>  <span class="comment">// Similar set of objects needed for the lms normalisation systematics</span></div>
<div class="line"><a name="l00304"></a><span class="lineno"> 304</span>  multi_array<double, 2> lms_k_hi_arr(extents[lms][m]);</div>
<div class="line"><a name="l00305"></a><span class="lineno"> 305</span>  multi_array<double, 2> lms_k_lo_arr(extents[lms][m]);</div>
<div class="line"><a name="l00306"></a><span class="lineno"> 306</span>  multi_array<std::shared_ptr<RooSpline1D>, 1> lms_spl_hi_arr(extents[lms]);</div>
<div class="line"><a name="l00307"></a><span class="lineno"> 307</span>  multi_array<std::shared_ptr<RooSpline1D>, 1> lms_spl_lo_arr(extents[lms]);</div>
<div class="line"><a name="l00308"></a><span class="lineno"> 308</span>  multi_array<std::shared_ptr<AsymPow>, 1> lms_asy_arr(extents[lms]);</div>
<div class="line"><a name="l00309"></a><span class="lineno"> 309</span>  </div>
<div class="line"><a name="l00310"></a><span class="lineno"> 310</span>  <span class="comment">// Now we'll fill these objects..</span></div>
<div class="line"><a name="l00311"></a><span class="lineno"> 311</span>  </div>
<div class="line"><a name="l00312"></a><span class="lineno"> 312</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00313"></a><span class="lineno"> 313</span>  <span class="comment">// Grab the nominal process histograms. We also have to convert every</span></div>
<div class="line"><a name="l00314"></a><span class="lineno"> 314</span>  <span class="comment">// histogram to a uniform integer binning, because this is what</span></div>
<div class="line"><a name="l00315"></a><span class="lineno"> 315</span>  <span class="comment">// text2workspace will do for all the non-morphed processes in our datacard,</span></div>
<div class="line"><a name="l00316"></a><span class="lineno"> 316</span>  <span class="comment">// and we need the binning of these to be in sync.</span></div>
<div class="line"><a name="l00317"></a><span class="lineno"> 317</span>  hist_arr[mi][0] =</div>
<div class="line"><a name="l00318"></a><span class="lineno"> 318</span>  std::make_shared<TH1F>(<a class="code" href="namespacech.html#a2cbcfd049543bfe9f63958bc671deeae">RebinHist</a>(<a class="code" href="namespacech.html#af9b462e52384bf17acf24c57bee7a5d3">AsTH1F</a>(pr_arr[mi]->shape())));</div>
<div class="line"><a name="l00319"></a><span class="lineno"> 319</span>  <span class="comment">// If the user supplied a TFile pointer we'll dump a bunch of info into it</span></div>
<div class="line"><a name="l00320"></a><span class="lineno"> 320</span>  <span class="comment">// for debugging</span></div>
<div class="line"><a name="l00321"></a><span class="lineno"> 321</span>  <span class="keywordflow">if</span> (file) {</div>
<div class="line"><a name="l00322"></a><span class="lineno"> 322</span>  file->WriteTObject(pr_arr[mi]->shape(), key + <span class="stringliteral">"_"</span> + m_str_vec[mi]);</div>
<div class="line"><a name="l00323"></a><span class="lineno"> 323</span>  }</div>
<div class="line"><a name="l00324"></a><span class="lineno"> 324</span>  <span class="comment">// Store the process rate</span></div>
<div class="line"><a name="l00325"></a><span class="lineno"> 325</span>  rate_arr[mi] = pr_arr[mi]->rate();</div>
<div class="line"><a name="l00326"></a><span class="lineno"> 326</span>  <span class="keyword">auto</span> proc_hist = pr_arr[mi]->ClonedScaledShape();</div>
<div class="line"><a name="l00327"></a><span class="lineno"> 327</span>  proc_hist->IntegralAndError(1, proc_hist->GetNbinsX(), rate_unc_arr[mi]);</div>
<div class="line"><a name="l00328"></a><span class="lineno"> 328</span>  <span class="comment">// Do the same for the Up and Down shapes</span></div>
<div class="line"><a name="l00329"></a><span class="lineno"> 329</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> ssi = 0; ssi < ss; ++ssi) {</div>
<div class="line"><a name="l00330"></a><span class="lineno"> 330</span>  hist_arr[mi][1 + 2 * ssi] =</div>
<div class="line"><a name="l00331"></a><span class="lineno"> 331</span>  std::make_shared<TH1F>(<a class="code" href="namespacech.html#a2cbcfd049543bfe9f63958bc671deeae">RebinHist</a>(<a class="code" href="namespacech.html#af9b462e52384bf17acf24c57bee7a5d3">AsTH1F</a>(ss_arr[ssi][mi]->shape_u())));</div>
<div class="line"><a name="l00332"></a><span class="lineno"> 332</span>  hist_arr[mi][2 + 2 * ssi] =</div>
<div class="line"><a name="l00333"></a><span class="lineno"> 333</span>  std::make_shared<TH1F>(<a class="code" href="namespacech.html#a2cbcfd049543bfe9f63958bc671deeae">RebinHist</a>(<a class="code" href="namespacech.html#af9b462e52384bf17acf24c57bee7a5d3">AsTH1F</a>(ss_arr[ssi][mi]->shape_d())));</div>
<div class="line"><a name="l00334"></a><span class="lineno"> 334</span>  <span class="keywordflow">if</span> (file) {</div>
<div class="line"><a name="l00335"></a><span class="lineno"> 335</span>  file->WriteTObject(ss_arr[ssi][mi]->shape_u(),</div>
<div class="line"><a name="l00336"></a><span class="lineno"> 336</span>  key + <span class="stringliteral">"_"</span> + m_str_vec[mi] + <span class="stringliteral">"_"</span> + ss_vec[ssi] + <span class="stringliteral">"Up"</span>);</div>
<div class="line"><a name="l00337"></a><span class="lineno"> 337</span>  file->WriteTObject(ss_arr[ssi][mi]->shape_d(),</div>
<div class="line"><a name="l00338"></a><span class="lineno"> 338</span>  key + <span class="stringliteral">"_"</span> + m_str_vec[mi] + <span class="stringliteral">"_"</span> + ss_vec[ssi] + <span class="stringliteral">"Down"</span>);</div>
<div class="line"><a name="l00339"></a><span class="lineno"> 339</span>  }</div>
<div class="line"><a name="l00340"></a><span class="lineno"> 340</span>  <span class="comment">// Store the uncertainty ("kappa") values for the shape systematics</span></div>
<div class="line"><a name="l00341"></a><span class="lineno"> 341</span>  ss_k_hi_arr[ssi][mi] = ss_arr[ssi][mi]->value_u();</div>
<div class="line"><a name="l00342"></a><span class="lineno"> 342</span>  ss_k_lo_arr[ssi][mi] = ss_arr[ssi][mi]->value_d();</div>
<div class="line"><a name="l00343"></a><span class="lineno"> 343</span>  <span class="comment">// For the normalisation we scale the kappa instead of putting the scaling</span></div>
<div class="line"><a name="l00344"></a><span class="lineno"> 344</span>  <span class="comment">// parameter as the variable</span></div>
<div class="line"><a name="l00345"></a><span class="lineno"> 345</span>  <span class="keywordflow">if</span> (std::fabs(ss_scale_arr[ssi] - 1.0) > 1E-6) {</div>
<div class="line"><a name="l00346"></a><span class="lineno"> 346</span>  ss_k_hi_arr[ssi][mi] = std::pow(ss_arr[ssi][mi]->value_u(), ss_scale_arr[ssi]);</div>
<div class="line"><a name="l00347"></a><span class="lineno"> 347</span>  ss_k_lo_arr[ssi][mi] = std::pow(ss_arr[ssi][mi]->value_d(), ss_scale_arr[ssi]);</div>
<div class="line"><a name="l00348"></a><span class="lineno"> 348</span>  }</div>
<div class="line"><a name="l00349"></a><span class="lineno"> 349</span>  }</div>
<div class="line"><a name="l00350"></a><span class="lineno"> 350</span>  <span class="comment">// And now the uncertainty values for the lnN systematics that vary with mass</span></div>
<div class="line"><a name="l00351"></a><span class="lineno"> 351</span>  <span class="comment">// We'll force these to be asymmetric even if they're not</span></div>
<div class="line"><a name="l00352"></a><span class="lineno"> 352</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> lmsi = 0; lmsi < lms; ++lmsi) {</div>
<div class="line"><a name="l00353"></a><span class="lineno"> 353</span>  lms_k_hi_arr[lmsi][mi] = lms_arr[lmsi][mi]->value_u();</div>
<div class="line"><a name="l00354"></a><span class="lineno"> 354</span>  <span class="keywordflow">if</span> (lms_arr[lmsi][mi]->asymm()) {</div>
<div class="line"><a name="l00355"></a><span class="lineno"> 355</span>  lms_k_lo_arr[lmsi][mi] = lms_arr[lmsi][mi]->value_d();</div>
<div class="line"><a name="l00356"></a><span class="lineno"> 356</span>  } <span class="keywordflow">else</span> {</div>
<div class="line"><a name="l00357"></a><span class="lineno"> 357</span>  lms_k_lo_arr[lmsi][mi] = 1. / lms_arr[lmsi][mi]->value_u();</div>
<div class="line"><a name="l00358"></a><span class="lineno"> 358</span>  }</div>
<div class="line"><a name="l00359"></a><span class="lineno"> 359</span>  }</div>
<div class="line"><a name="l00360"></a><span class="lineno"> 360</span>  }</div>
<div class="line"><a name="l00361"></a><span class="lineno"> 361</span>  </div>
<div class="line"><a name="l00362"></a><span class="lineno"> 362</span>  <span class="comment">// Now we've made all our histograms, we'll put them in the TList format</span></div>
<div class="line"><a name="l00363"></a><span class="lineno"> 363</span>  <span class="comment">// of [nominal, syst_1_Up, syst_1_Down, ... , syst_N_Up, syst_N_Down]</span></div>
<div class="line"><a name="l00364"></a><span class="lineno"> 364</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00365"></a><span class="lineno"> 365</span>  list_arr[mi] = std::make_shared<TList>();</div>
<div class="line"><a name="l00366"></a><span class="lineno"> 366</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> xi = 0; xi < (1 + ss * 2); ++xi) {</div>
<div class="line"><a name="l00367"></a><span class="lineno"> 367</span>  list_arr[mi]->Add(hist_arr[mi][xi].get());</div>
<div class="line"><a name="l00368"></a><span class="lineno"> 368</span>  }</div>
<div class="line"><a name="l00369"></a><span class="lineno"> 369</span>  }</div>
<div class="line"><a name="l00370"></a><span class="lineno"> 370</span>  </div>
<div class="line"><a name="l00371"></a><span class="lineno"> 371</span>  <span class="comment">// Print the values of the yields and kappa factors that will be inputs</span></div>
<div class="line"><a name="l00372"></a><span class="lineno"> 372</span>  <span class="comment">// to our spline interpolation objects</span></div>
<div class="line"><a name="l00373"></a><span class="lineno"> 373</span>  <span class="keywordflow">if</span> (verbose) {</div>
<div class="line"><a name="l00374"></a><span class="lineno"> 374</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00375"></a><span class="lineno"> 375</span>  std::cout << boost::format(<span class="stringliteral">"%-10s"</span>) % m_str_vec[mi];</div>
<div class="line"><a name="l00376"></a><span class="lineno"> 376</span>  }</div>
<div class="line"><a name="l00377"></a><span class="lineno"> 377</span>  std::cout << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00378"></a><span class="lineno"> 378</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00379"></a><span class="lineno"> 379</span>  std::cout << boost::format(<span class="stringliteral">"%-10.5g"</span>) % rate_arr[mi];</div>
<div class="line"><a name="l00380"></a><span class="lineno"> 380</span>  }</div>
<div class="line"><a name="l00381"></a><span class="lineno"> 381</span>  std::cout << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00382"></a><span class="lineno"> 382</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> ssi = 0; ssi < ss; ++ssi) {</div>
<div class="line"><a name="l00383"></a><span class="lineno"> 383</span>  std::cout << ss_vec[ssi] << <span class="stringliteral">" Up"</span> << std::endl;</div>
<div class="line"><a name="l00384"></a><span class="lineno"> 384</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00385"></a><span class="lineno"> 385</span>  std::cout << boost::format(<span class="stringliteral">"%-10.5g"</span>) % ss_k_hi_arr[ssi][mi];</div>
<div class="line"><a name="l00386"></a><span class="lineno"> 386</span>  }</div>
<div class="line"><a name="l00387"></a><span class="lineno"> 387</span>  std::cout << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00388"></a><span class="lineno"> 388</span>  std::cout << ss_vec[ssi] << <span class="stringliteral">" Down"</span> << std::endl;</div>
<div class="line"><a name="l00389"></a><span class="lineno"> 389</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00390"></a><span class="lineno"> 390</span>  std::cout << boost::format(<span class="stringliteral">"%-10.5g"</span>) % ss_k_lo_arr[ssi][mi];</div>
<div class="line"><a name="l00391"></a><span class="lineno"> 391</span>  }</div>
<div class="line"><a name="l00392"></a><span class="lineno"> 392</span>  std::cout << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00393"></a><span class="lineno"> 393</span>  }</div>
<div class="line"><a name="l00394"></a><span class="lineno"> 394</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> lmsi = 0; lmsi < lms; ++lmsi) {</div>
<div class="line"><a name="l00395"></a><span class="lineno"> 395</span>  std::cout << lms_vec[lmsi] << <span class="stringliteral">" Up"</span> << std::endl;</div>
<div class="line"><a name="l00396"></a><span class="lineno"> 396</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00397"></a><span class="lineno"> 397</span>  std::cout << boost::format(<span class="stringliteral">"%-10.5g"</span>) % lms_k_hi_arr[lmsi][mi];</div>
<div class="line"><a name="l00398"></a><span class="lineno"> 398</span>  }</div>
<div class="line"><a name="l00399"></a><span class="lineno"> 399</span>  std::cout << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00400"></a><span class="lineno"> 400</span>  std::cout << lms_vec[lmsi] << <span class="stringliteral">" Down"</span> << std::endl;</div>
<div class="line"><a name="l00401"></a><span class="lineno"> 401</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00402"></a><span class="lineno"> 402</span>  std::cout << boost::format(<span class="stringliteral">"%-10.5g"</span>) % lms_k_lo_arr[lmsi][mi];</div>
<div class="line"><a name="l00403"></a><span class="lineno"> 403</span>  }</div>
<div class="line"><a name="l00404"></a><span class="lineno"> 404</span>  std::cout << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00405"></a><span class="lineno"> 405</span>  }</div>
<div class="line"><a name="l00406"></a><span class="lineno"> 406</span>  }</div>
<div class="line"><a name="l00407"></a><span class="lineno"> 407</span>  </div>
<div class="line"><a name="l00408"></a><span class="lineno"> 408</span>  <span class="comment">// Can do more sophistical spline interpolation if we want, but let's use</span></div>
<div class="line"><a name="l00409"></a><span class="lineno"> 409</span>  <span class="comment">// simple LINEAR interpolation for now</span></div>
<div class="line"><a name="l00410"></a><span class="lineno"> 410</span>  TString interp = <span class="stringliteral">"LINEAR"</span>;</div>
<div class="line"><a name="l00411"></a><span class="lineno"> 411</span>  </div>
<div class="line"><a name="l00412"></a><span class="lineno"> 412</span>  <span class="comment">// Here when the force_template_limit option is requested </span></div>
<div class="line"><a name="l00413"></a><span class="lineno"> 413</span>  <span class="comment">// we have to add some extra terms to the vector of masses to ensure that </span></div>
<div class="line"><a name="l00414"></a><span class="lineno"> 414</span>  <span class="comment">// the signal pdf goes to 0 outside of the MC template range.</span></div>
<div class="line"><a name="l00415"></a><span class="lineno"> 415</span>  </div>
<div class="line"><a name="l00416"></a><span class="lineno"> 416</span>  vector<double> new_m_vec(m_vec);</div>
<div class="line"><a name="l00417"></a><span class="lineno"> 417</span>  <span class="comment">// Insert an entry at either end of the vector for a mass just slightly</span></div>
<div class="line"><a name="l00418"></a><span class="lineno"> 418</span>  <span class="comment">// outside of the range</span></div>
<div class="line"><a name="l00419"></a><span class="lineno"> 419</span>  new_m_vec.insert(new_m_vec.begin(),m_vec[0]-1E-6);</div>
<div class="line"><a name="l00420"></a><span class="lineno"> 420</span>  new_m_vec.push_back(m_vec[m-1]+1E-6);</div>
<div class="line"><a name="l00421"></a><span class="lineno"> 421</span>  <span class="comment">// Create a corresponding rate array with 0 entries for these new masses</span></div>
<div class="line"><a name="l00422"></a><span class="lineno"> 422</span>  multi_array<double, 1> new_rate_arr(extents[m+2]);</div>
<div class="line"><a name="l00423"></a><span class="lineno"> 423</span>  new_rate_arr[0] = 0.0;</div>
<div class="line"><a name="l00424"></a><span class="lineno"> 424</span>  <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> i = 0; i < m; ++i) new_rate_arr[i+1] = rate_arr[i] ;</div>
<div class="line"><a name="l00425"></a><span class="lineno"> 425</span>  new_rate_arr[m+1] = 0.0;</div>
<div class="line"><a name="l00426"></a><span class="lineno"> 426</span>  </div>
<div class="line"><a name="l00427"></a><span class="lineno"> 427</span>  <span class="keywordflow">if</span> (verbose && force_template_limit) {</div>
<div class="line"><a name="l00428"></a><span class="lineno"> 428</span>  std::cout << <span class="stringliteral">">>>> Forcing rate to 0 outside of template range:"</span> << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00429"></a><span class="lineno"> 429</span>  <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> mi = 0; mi < m+2; ++mi) {</div>
<div class="line"><a name="l00430"></a><span class="lineno"> 430</span>  std::cout << boost::format(<span class="stringliteral">"%-10.5g"</span>) % new_m_vec[mi];</div>
<div class="line"><a name="l00431"></a><span class="lineno"> 431</span>  }</div>
<div class="line"><a name="l00432"></a><span class="lineno"> 432</span>  std::cout << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00433"></a><span class="lineno"> 433</span>  <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> mi = 0; mi < m+2; ++mi) {</div>
<div class="line"><a name="l00434"></a><span class="lineno"> 434</span>  std::cout << boost::format(<span class="stringliteral">"%-10.5g"</span>) % new_rate_arr[mi];</div>
<div class="line"><a name="l00435"></a><span class="lineno"> 435</span>  }</div>
<div class="line"><a name="l00436"></a><span class="lineno"> 436</span>  std::cout << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00437"></a><span class="lineno"> 437</span>  }</div>
<div class="line"><a name="l00438"></a><span class="lineno"> 438</span>  <span class="comment">// Create the 1D spline directly from the rate array</span></div>
<div class="line"><a name="l00440"></a><span class="lineno"> 440</span> <span class="comment"></span> RooSpline1D rate_spline(<span class="stringliteral">"interp_rate_"</span>+key, <span class="stringliteral">""</span>, mass_var, </div>
<div class="line"><a name="l00441"></a><span class="lineno"> 441</span>  force_template_limit ? m+2 : m, </div>
<div class="line"><a name="l00442"></a><span class="lineno"> 442</span>  force_template_limit ? new_m_vec.data() : m_vec.data(),</div>
<div class="line"><a name="l00443"></a><span class="lineno"> 443</span>  force_template_limit ? new_rate_arr.data() : rate_arr.data(),</div>
<div class="line"><a name="l00444"></a><span class="lineno"> 444</span>  interp);</div>
<div class="line"><a name="l00446"></a><span class="lineno"> 446</span>  </div>
<div class="line"><a name="l00447"></a><span class="lineno"> 447</span>  <span class="keywordflow">if</span> (file) {</div>
<div class="line"><a name="l00448"></a><span class="lineno"> 448</span>  TGraphErrors tmp(m, m_vec.data(), rate_arr.data(), <span class="keyword">nullptr</span>, rate_unc_arr.data());</div>
<div class="line"><a name="l00449"></a><span class="lineno"> 449</span>  file->WriteTObject(&tmp, <span class="stringliteral">"interp_rate_"</span>+key);</div>
<div class="line"><a name="l00450"></a><span class="lineno"> 450</span>  }</div>
<div class="line"><a name="l00451"></a><span class="lineno"> 451</span>  <span class="comment">// Collect all terms that will go into the total normalisation:</span></div>
<div class="line"><a name="l00452"></a><span class="lineno"> 452</span>  <span class="comment">// nominal * systeff_1 * systeff_2 * ... * systeff_N</span></div>
<div class="line"><a name="l00453"></a><span class="lineno"> 453</span>  RooArgList rate_prod(rate_spline);</div>
<div class="line"><a name="l00454"></a><span class="lineno"> 454</span>  <span class="comment">// For each shape systematic build a 1D spline for kappa_hi and kappa_lo</span></div>
<div class="line"><a name="l00455"></a><span class="lineno"> 455</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> ssi = 0; ssi < ss; ++ssi) {</div>
<div class="line"><a name="l00456"></a><span class="lineno"> 456</span>  ss_spl_hi_arr[ssi] = std::make_shared<RooSpline1D>(<span class="stringliteral">"spline_hi_"</span> +</div>
<div class="line"><a name="l00457"></a><span class="lineno"> 457</span>  key + <span class="stringliteral">"_"</span> + ss_vec[ssi], <span class="stringliteral">""</span>, mass_var, m, m_vec.data(),</div>
<div class="line"><a name="l00458"></a><span class="lineno"> 458</span>  ss_k_hi_arr[ssi].origin(), interp);</div>
<div class="line"><a name="l00459"></a><span class="lineno"> 459</span>  ss_spl_lo_arr[ssi] = std::make_shared<RooSpline1D>(<span class="stringliteral">"spline_lo_"</span> +</div>
<div class="line"><a name="l00460"></a><span class="lineno"> 460</span>  key + <span class="stringliteral">"_"</span> + ss_vec[ssi], <span class="stringliteral">""</span>, mass_var, m, m_vec.data(),</div>
<div class="line"><a name="l00461"></a><span class="lineno"> 461</span>  ss_k_lo_arr[ssi].origin(), interp);</div>
<div class="line"><a name="l00462"></a><span class="lineno"> 462</span>  <span class="keywordflow">if</span> (file) {</div>
<div class="line"><a name="l00463"></a><span class="lineno"> 463</span>  TGraph tmp_hi(m, m_vec.data(), ss_k_hi_arr[ssi].origin());</div>
<div class="line"><a name="l00464"></a><span class="lineno"> 464</span>  file->WriteTObject(&tmp_hi, <span class="stringliteral">"spline_hi_"</span> + key + <span class="stringliteral">"_"</span> + ss_vec[ssi]);</div>
<div class="line"><a name="l00465"></a><span class="lineno"> 465</span>  TGraph tmp_lo(m, m_vec.data(), ss_k_lo_arr[ssi].origin());</div>
<div class="line"><a name="l00466"></a><span class="lineno"> 466</span>  file->WriteTObject(&tmp_lo, <span class="stringliteral">"spline_lo_"</span> + key + <span class="stringliteral">"_"</span> + ss_vec[ssi]);</div>
<div class="line"><a name="l00467"></a><span class="lineno"> 467</span>  }</div>
<div class="line"><a name="l00468"></a><span class="lineno"> 468</span>  <span class="comment">// Then build the AsymPow object for each systematic as a function of the</span></div>
<div class="line"><a name="l00469"></a><span class="lineno"> 469</span>  <span class="comment">// kappas and the nuisance parameter</span></div>
<div class="line"><a name="l00470"></a><span class="lineno"> 470</span>  ss_asy_arr[ssi] = std::make_shared<AsymPow>(<span class="stringliteral">"systeff_"</span> +</div>
<div class="line"><a name="l00471"></a><span class="lineno"> 471</span>  key + <span class="stringliteral">"_"</span> + ss_vec[ssi], <span class="stringliteral">""</span>,</div>
<div class="line"><a name="l00472"></a><span class="lineno"> 472</span>  *(ss_spl_lo_arr[ssi]), *(ss_spl_hi_arr[ssi]),</div>
<div class="line"><a name="l00473"></a><span class="lineno"> 473</span>  *(ss_scale_var_arr[ssi]));</div>
<div class="line"><a name="l00474"></a><span class="lineno"> 474</span>  rate_prod.add(*(ss_asy_arr[ssi]));</div>
<div class="line"><a name="l00475"></a><span class="lineno"> 475</span>  }</div>
<div class="line"><a name="l00476"></a><span class="lineno"> 476</span>  <span class="comment">// Same procedure for the lms normalisation systematics: build the splines</span></div>
<div class="line"><a name="l00477"></a><span class="lineno"> 477</span>  <span class="comment">// then the AsymPows and add to the rate_prod list</span></div>
<div class="line"><a name="l00478"></a><span class="lineno"> 478</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> lmsi = 0; lmsi < lms; ++lmsi) {</div>
<div class="line"><a name="l00479"></a><span class="lineno"> 479</span>  lms_spl_hi_arr[lmsi] = std::make_shared<RooSpline1D>(<span class="stringliteral">"spline_hi_"</span> +</div>
<div class="line"><a name="l00480"></a><span class="lineno"> 480</span>  key + <span class="stringliteral">"_"</span> + lms_vec[lmsi], <span class="stringliteral">""</span>, mass_var, m, m_vec.data(),</div>
<div class="line"><a name="l00481"></a><span class="lineno"> 481</span>  lms_k_hi_arr[lmsi].origin(), interp);</div>
<div class="line"><a name="l00482"></a><span class="lineno"> 482</span>  lms_spl_lo_arr[lmsi] = std::make_shared<RooSpline1D>(<span class="stringliteral">"spline_lo_"</span> +</div>
<div class="line"><a name="l00483"></a><span class="lineno"> 483</span>  key + <span class="stringliteral">"_"</span> + lms_vec[lmsi], <span class="stringliteral">""</span>, mass_var, m, m_vec.data(),</div>
<div class="line"><a name="l00484"></a><span class="lineno"> 484</span>  lms_k_lo_arr[lmsi].origin(), interp);</div>
<div class="line"><a name="l00485"></a><span class="lineno"> 485</span>  lms_asy_arr[lmsi] = std::make_shared<AsymPow>(<span class="stringliteral">"systeff_"</span> +</div>
<div class="line"><a name="l00486"></a><span class="lineno"> 486</span>  key + <span class="stringliteral">"_"</span> + lms_vec[lmsi], <span class="stringliteral">""</span>, *(lms_spl_lo_arr[lmsi]),</div>
<div class="line"><a name="l00487"></a><span class="lineno"> 487</span>  *(lms_spl_hi_arr[lmsi]), *(lms_var_arr[lmsi]));</div>
<div class="line"><a name="l00488"></a><span class="lineno"> 488</span>  rate_prod.add(*(lms_asy_arr[lmsi]));</div>
<div class="line"><a name="l00489"></a><span class="lineno"> 489</span>  }</div>
<div class="line"><a name="l00490"></a><span class="lineno"> 490</span>  <span class="comment">// We'll come back to this rate_prod a bit later.</span></div>
<div class="line"><a name="l00491"></a><span class="lineno"> 491</span>  </div>
<div class="line"><a name="l00492"></a><span class="lineno"> 492</span>  <span class="comment">// Now some fun with binning. We anticipate that the process histograms we</span></div>
<div class="line"><a name="l00493"></a><span class="lineno"> 493</span>  <span class="comment">// have been supplied could have a finer binning than is actually wanted for</span></div>
<div class="line"><a name="l00494"></a><span class="lineno"> 494</span>  <span class="comment">// the analysis (and the fit), in order to avoid known problems with the RMS</span></div>
<div class="line"><a name="l00495"></a><span class="lineno"> 495</span>  <span class="comment">// of a peaking distribution not being morphed smoothly from mass point to</span></div>
<div class="line"><a name="l00496"></a><span class="lineno"> 496</span>  <span class="comment">// mass point if the binning is too wide. The RooMorphingPdf will handle</span></div>
<div class="line"><a name="l00497"></a><span class="lineno"> 497</span>  <span class="comment">// re-binning on the fly, but we have to tell it how to rebin. To do this we</span></div>
<div class="line"><a name="l00498"></a><span class="lineno"> 498</span>  <span class="comment">// assume the observed data histogram has the target binning.</span></div>
<div class="line"><a name="l00499"></a><span class="lineno"> 499</span>  TH1F data_hist = cb_bp.<a class="code" href="classch_1_1_combine_harvester.html#a0d3eca727963c0c14a3c04bbdce2ca5e">GetObservedShape</a>();</div>
<div class="line"><a name="l00500"></a><span class="lineno"> 500</span>  TH1F proc_hist = cb_bp.<a class="code" href="classch_1_1_combine_harvester.html#a3d265eb50ea8ec70192788d985733d9d">GetShape</a>();</div>
<div class="line"><a name="l00501"></a><span class="lineno"> 501</span>  <span class="comment">// The x-axis variable has to be called "CMS_th1x", as this is what</span></div>
<div class="line"><a name="l00502"></a><span class="lineno"> 502</span>  <span class="comment">// text2workspace will use for all the normal processes</span></div>
<div class="line"><a name="l00503"></a><span class="lineno"> 503</span>  RooRealVar xvar(<span class="stringliteral">"CMS_th1x"</span>, <span class="stringliteral">"CMS_th1x"</span>, 0,</div>
<div class="line"><a name="l00504"></a><span class="lineno"> 504</span>  <span class="keyword">static_cast<</span><span class="keywordtype">float</span><span class="keyword">></span>(data_hist.GetNbinsX()));</div>
<div class="line"><a name="l00505"></a><span class="lineno"> 505</span>  xvar.setBins(data_hist.GetNbinsX());</div>
<div class="line"><a name="l00506"></a><span class="lineno"> 506</span>  </div>
<div class="line"><a name="l00507"></a><span class="lineno"> 507</span>  <span class="comment">// Create a second x-axis variable, named specific to the bin, that will be</span></div>
<div class="line"><a name="l00508"></a><span class="lineno"> 508</span>  <span class="comment">// for the finer-binned input</span></div>
<div class="line"><a name="l00509"></a><span class="lineno"> 509</span>  RooRealVar morph_xvar((<span class="stringliteral">"CMS_th1x_"</span>+bin).c_str(), <span class="stringliteral">""</span>, 0,</div>
<div class="line"><a name="l00510"></a><span class="lineno"> 510</span>  <span class="keyword">static_cast<</span><span class="keywordtype">float</span><span class="keyword">></span>(proc_hist.GetNbinsX()));</div>
<div class="line"><a name="l00511"></a><span class="lineno"> 511</span>  <span class="comment">// We're not going to need roofit to evaluate anything as a function of this</span></div>
<div class="line"><a name="l00512"></a><span class="lineno"> 512</span>  <span class="comment">// morphing x-axis variable, so we set it constant</span></div>
<div class="line"><a name="l00513"></a><span class="lineno"> 513</span>  morph_xvar.setConstant();</div>
<div class="line"><a name="l00514"></a><span class="lineno"> 514</span>  morph_xvar.setBins(proc_hist.GetNbinsX());</div>
<div class="line"><a name="l00515"></a><span class="lineno"> 515</span>  <span class="keywordflow">if</span> (verbose) {</div>
<div class="line"><a name="l00516"></a><span class="lineno"> 516</span>  xvar.Print();</div>
<div class="line"><a name="l00517"></a><span class="lineno"> 517</span>  morph_xvar.Print();</div>
<div class="line"><a name="l00518"></a><span class="lineno"> 518</span>  }</div>
<div class="line"><a name="l00519"></a><span class="lineno"> 519</span>  </div>
<div class="line"><a name="l00520"></a><span class="lineno"> 520</span>  <span class="comment">// Now we can build the array of vertical-interpolation pdfs (the same that</span></div>
<div class="line"><a name="l00521"></a><span class="lineno"> 521</span>  <span class="comment">// text2workspace builds for every process), one per mass point</span></div>
<div class="line"><a name="l00522"></a><span class="lineno"> 522</span>  multi_array<std::shared_ptr<FastVerticalInterpHistPdf2>, 1> vpdf_arr(</div>
<div class="line"><a name="l00523"></a><span class="lineno"> 523</span>  extents[m]);</div>
<div class="line"><a name="l00524"></a><span class="lineno"> 524</span>  RooArgList vpdf_list;</div>
<div class="line"><a name="l00525"></a><span class="lineno"> 525</span>  </div>
<div class="line"><a name="l00526"></a><span class="lineno"> 526</span>  TString vert_name = key + <span class="stringliteral">"_"</span>;</div>
<div class="line"><a name="l00527"></a><span class="lineno"> 527</span>  </div>
<div class="line"><a name="l00528"></a><span class="lineno"> 528</span>  <span class="comment">// Follow what ShapeTools.py does and set the smoothing region</span></div>
<div class="line"><a name="l00529"></a><span class="lineno"> 529</span>  <span class="comment">// to the minimum of all of the shape scales</span></div>
<div class="line"><a name="l00530"></a><span class="lineno"> 530</span>  <span class="keywordtype">double</span> qrange = 1.;</div>
<div class="line"><a name="l00531"></a><span class="lineno"> 531</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> ssi = 0; ssi < ss; ++ssi) {</div>
<div class="line"><a name="l00532"></a><span class="lineno"> 532</span>  <span class="keywordflow">if</span> (ss_scale_arr[ssi] < qrange) qrange = ss_scale_arr[ssi];</div>
<div class="line"><a name="l00533"></a><span class="lineno"> 533</span>  }</div>
<div class="line"><a name="l00534"></a><span class="lineno"> 534</span>  </div>
<div class="line"><a name="l00535"></a><span class="lineno"> 535</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> mi = 0; mi < m; ++mi) {</div>
<div class="line"><a name="l00536"></a><span class="lineno"> 536</span>  <span class="comment">// Construct it with the right binning, the right histograms and the right</span></div>
<div class="line"><a name="l00537"></a><span class="lineno"> 537</span>  <span class="comment">// scaling parameters</span></div>
<div class="line"><a name="l00538"></a><span class="lineno"> 538</span>  vpdf_arr[mi] = std::make_shared<FastVerticalInterpHistPdf2>(</div>
<div class="line"><a name="l00539"></a><span class="lineno"> 539</span>  vert_name + m_str_vec[mi] + <span class="stringliteral">"_vmorph"</span>, <span class="stringliteral">""</span>, morph_xvar, *(list_arr[mi]),</div>
<div class="line"><a name="l00540"></a><span class="lineno"> 540</span>  ss_list, qrange, 0);</div>
<div class="line"><a name="l00541"></a><span class="lineno"> 541</span>  <span class="comment">// Add it to a list that we'll supply to the RooMorphingPdf</span></div>
<div class="line"><a name="l00542"></a><span class="lineno"> 542</span>  vpdf_list.add(*(vpdf_arr[mi]));</div>
<div class="line"><a name="l00543"></a><span class="lineno"> 543</span>  }</div>
<div class="line"><a name="l00544"></a><span class="lineno"> 544</span>  TString morph_name = key + <span class="stringliteral">"_morph"</span>;</div>
<div class="line"><a name="l00545"></a><span class="lineno"> 545</span>  <span class="comment">// At long last, we can build our pdf, giving it:</span></div>
<div class="line"><a name="l00546"></a><span class="lineno"> 546</span>  <span class="comment">// xvar : the fixed "CMS_th1x" x-axis variable with uniform integer binning</span></div>
<div class="line"><a name="l00547"></a><span class="lineno"> 547</span>  <span class="comment">// mass_var: the floating mass value for the interpolation</span></div>
<div class="line"><a name="l00548"></a><span class="lineno"> 548</span>  <span class="comment">// vpdf_list: the list of vertical-interpolation pdfs</span></div>
<div class="line"><a name="l00549"></a><span class="lineno"> 549</span>  <span class="comment">// m_vec: the corresponding list of mass points</span></div>
<div class="line"><a name="l00550"></a><span class="lineno"> 550</span>  <span class="comment">// allow_morph: if false will just evaluate to the closest pdf in mass</span></div>
<div class="line"><a name="l00551"></a><span class="lineno"> 551</span>  <span class="comment">// data_hist.GetXaxis(): The original (non-uniform) target binning</span></div>
<div class="line"><a name="l00552"></a><span class="lineno"> 552</span>  <span class="comment">// proc_hist.GetXaxis(): The original (non-uniform) morphing binning</span></div>
<div class="line"><a name="l00554"></a><span class="lineno"> 554</span> <span class="comment"></span> RooMorphingPdf morph_pdf(morph_name, <span class="stringliteral">""</span>, xvar, mass_var, vpdf_list,</div>
<div class="line"><a name="l00555"></a><span class="lineno"> 555</span>  m_vec, allow_morph, *(data_hist.GetXaxis()),</div>
<div class="line"><a name="l00556"></a><span class="lineno"> 556</span>  *(proc_hist.GetXaxis()));</div>
<div class="line"><a name="l00558"></a><span class="lineno"> 558</span>  <span class="comment">// And we can make the final normalisation product</span></div>
<div class="line"><a name="l00559"></a><span class="lineno"> 559</span>  <span class="comment">// The value of norm_postfix is very important. text2workspace will only look</span></div>
<div class="line"><a name="l00560"></a><span class="lineno"> 560</span>  <span class="comment">// for for a term with the pdf name + "_norm". But it might be the user wants</span></div>
<div class="line"><a name="l00561"></a><span class="lineno"> 561</span>  <span class="comment">// to add even more terms to this total normalisation, so we give them the option</span></div>
<div class="line"><a name="l00562"></a><span class="lineno"> 562</span>  <span class="comment">// of using some other suffix.</span></div>
<div class="line"><a name="l00564"></a><span class="lineno"> 564</span> <span class="comment"></span> RooProduct morph_rate(morph_name + <span class="stringliteral">"_"</span> + TString(norm_postfix), <span class="stringliteral">""</span>,</div>
<div class="line"><a name="l00565"></a><span class="lineno"> 565</span>  rate_prod);</div>
<div class="line"><a name="l00567"></a><span class="lineno"> 567</span>  </div>
<div class="line"><a name="l00568"></a><span class="lineno"> 568</span>  <span class="comment">// Dump even more plots</span></div>
<div class="line"><a name="l00569"></a><span class="lineno"> 569</span>  <span class="keywordflow">if</span> (file) <a class="code" href="namespacech.html#a21ce0ae95c56bde7e1ecc1f21a839d2d">MakeMorphDebugPlots</a>(&morph_pdf, &mass_var, m_vec, file, &data_hist);</div>
<div class="line"><a name="l00570"></a><span class="lineno"> 570</span>  </div>
<div class="line"><a name="l00571"></a><span class="lineno"> 571</span>  <span class="comment">// Load our pdf and norm objects into the workspace</span></div>
<div class="line"><a name="l00572"></a><span class="lineno"> 572</span>  ws.import(morph_pdf, RooFit::RecycleConflictNodes());</div>
<div class="line"><a name="l00573"></a><span class="lineno"> 573</span>  ws.import(morph_rate, RooFit::RecycleConflictNodes());</div>
<div class="line"><a name="l00574"></a><span class="lineno"> 574</span>  </div>
<div class="line"><a name="l00575"></a><span class="lineno"> 575</span>  <span class="keywordflow">if</span> (!verbose) RooMsgService::instance().setGlobalKillBelow(backup_msg_level);</div>
<div class="line"><a name="l00576"></a><span class="lineno"> 576</span>  </div>
<div class="line"><a name="l00577"></a><span class="lineno"> 577</span>  <span class="comment">// Now we can clean up the CH instance a bit</span></div>
<div class="line"><a name="l00578"></a><span class="lineno"> 578</span>  <span class="comment">// We only need one entry for each Process or Systematic now, not one per mass point</span></div>
<div class="line"><a name="l00579"></a><span class="lineno"> 579</span>  <span class="comment">// We'll modify the first mass point to house our new pdf and norm, and drop</span></div>
<div class="line"><a name="l00580"></a><span class="lineno"> 580</span>  <span class="comment">// the rest.</span></div>
<div class="line"><a name="l00581"></a><span class="lineno"> 581</span>  std::string mass_min = m_str_vec.at(0);</div>
<div class="line"><a name="l00582"></a><span class="lineno"> 582</span>  </div>
<div class="line"><a name="l00583"></a><span class="lineno"> 583</span>  <span class="comment">// Dump Process entries for other mass points</span></div>
<div class="line"><a name="l00584"></a><span class="lineno"> 584</span>  cb.<a class="code" href="classch_1_1_combine_harvester.html#aab6a698c192aa091178bbd8dfbd0266f">FilterProcs</a>([&](<a class="code" href="classch_1_1_process.html">ch::Process</a> <span class="keyword">const</span>* p) {</div>
<div class="line"><a name="l00585"></a><span class="lineno"> 585</span>  <span class="keywordflow">return</span> p-><a class="code" href="classch_1_1_object.html#a7b588d78c22464ebe2b793cbaa2366d3">bin</a>() == bin && p-><a class="code" href="classch_1_1_object.html#a740d853956df951a7cb08094e21f0332">process</a>() == <a class="code" href="structch_1_1syst_1_1process.html">process</a> && p-><a class="code" href="classch_1_1_object.html#ae2f9e80af391e105a6b6bc27107e3c9e">mass</a>() != mass_min;</div>
<div class="line"><a name="l00586"></a><span class="lineno"> 586</span>  });</div>
<div class="line"><a name="l00587"></a><span class="lineno"> 587</span>  <span class="comment">// Dump Systematic entries for other mass points, but only if the type is shape</span></div>
<div class="line"><a name="l00588"></a><span class="lineno"> 588</span>  <span class="comment">// or a lnN that we found varied as a function of the mass. We've already built</span></div>
<div class="line"><a name="l00589"></a><span class="lineno"> 589</span>  <span class="comment">// these uncertainties into our normalisation term. Constant lnN systematics can</span></div>
<div class="line"><a name="l00590"></a><span class="lineno"> 590</span>  <span class="comment">// remain in the datacard and be added by text2workspace</span></div>
<div class="line"><a name="l00591"></a><span class="lineno"> 591</span>  cb.<a class="code" href="classch_1_1_combine_harvester.html#a91821b9464ebaacb5b92be2b3bc794f8">FilterSysts</a>([&](<a class="code" href="classch_1_1_systematic.html">ch::Systematic</a> <span class="keyword">const</span>* n) {</div>
<div class="line"><a name="l00592"></a><span class="lineno"> 592</span>  <span class="keywordflow">return</span> (n-><a class="code" href="classch_1_1_object.html#a7b588d78c22464ebe2b793cbaa2366d3">bin</a>() == bin && n-><a class="code" href="classch_1_1_object.html#a740d853956df951a7cb08094e21f0332">process</a>() == <a class="code" href="structch_1_1syst_1_1process.html">process</a>) &&</div>
<div class="line"><a name="l00593"></a><span class="lineno"> 593</span>  ((n-><a class="code" href="classch_1_1_object.html#ae2f9e80af391e105a6b6bc27107e3c9e">mass</a>() != mass_min) || (n-><a class="code" href="classch_1_1_systematic.html#a2c91f65744afb7460bfffdf2a4987196">type</a>() == <span class="stringliteral">"shape"</span>) ||</div>
<div class="line"><a name="l00594"></a><span class="lineno"> 594</span>  (lms_set.count(n-><a class="code" href="classch_1_1_systematic.html#a872469110fe349da45922431d6d24bc7">name</a>())));</div>
<div class="line"><a name="l00595"></a><span class="lineno"> 595</span>  });</div>
<div class="line"><a name="l00596"></a><span class="lineno"> 596</span>  <span class="comment">// With the remaining Process entry (should only be one if we did this right),</span></div>
<div class="line"><a name="l00597"></a><span class="lineno"> 597</span>  <span class="comment">// Make the mass generic ("*"), drop the TH1 and set the rate to 1.0, as this</span></div>
<div class="line"><a name="l00598"></a><span class="lineno"> 598</span>  <span class="comment">// will now be read from our norm object</span></div>
<div class="line"><a name="l00599"></a><span class="lineno"> 599</span>  cb.<a class="code" href="classch_1_1_combine_harvester.html#a44064524ffe47ccac63ef8025a6f46dc">ForEachProc</a>([&](<a class="code" href="classch_1_1_process.html">ch::Process</a> * p) {</div>
<div class="line"><a name="l00600"></a><span class="lineno"> 600</span>  <span class="keywordflow">if</span> (p-><a class="code" href="classch_1_1_object.html#a7b588d78c22464ebe2b793cbaa2366d3">bin</a>() == bin && p-><a class="code" href="classch_1_1_object.html#a740d853956df951a7cb08094e21f0332">process</a>() == <a class="code" href="structch_1_1syst_1_1process.html">process</a>) {</div>
<div class="line"><a name="l00601"></a><span class="lineno"> 601</span>  p->set_mass(<span class="stringliteral">"*"</span>);</div>
<div class="line"><a name="l00602"></a><span class="lineno"> 602</span>  p->set_shape(nullptr, false);</div>
<div class="line"><a name="l00603"></a><span class="lineno"> 603</span>  p->set_rate(1.0);</div>
<div class="line"><a name="l00604"></a><span class="lineno"> 604</span>  }</div>
<div class="line"><a name="l00605"></a><span class="lineno"> 605</span>  });</div>
<div class="line"><a name="l00606"></a><span class="lineno"> 606</span>  <span class="comment">// Just declare the mass to be generic for the remaining systematics</span></div>
<div class="line"><a name="l00607"></a><span class="lineno"> 607</span>  cb.<a class="code" href="classch_1_1_combine_harvester.html#a3b73ec07ca4de14bf76383b2921f344a">ForEachSyst</a>([&](<a class="code" href="classch_1_1_systematic.html">ch::Systematic</a> * n) {</div>
<div class="line"><a name="l00608"></a><span class="lineno"> 608</span>  <span class="keywordflow">if</span> (n-><a class="code" href="classch_1_1_object.html#a7b588d78c22464ebe2b793cbaa2366d3">bin</a>() == bin && n-><a class="code" href="classch_1_1_object.html#a740d853956df951a7cb08094e21f0332">process</a>() == <a class="code" href="structch_1_1syst_1_1process.html">process</a>) {</div>
<div class="line"><a name="l00609"></a><span class="lineno"> 609</span>  n->set_mass(<span class="stringliteral">"*"</span>);</div>
<div class="line"><a name="l00610"></a><span class="lineno"> 610</span>  }</div>
<div class="line"><a name="l00611"></a><span class="lineno"> 611</span>  });</div>
<div class="line"><a name="l00612"></a><span class="lineno"> 612</span>  </div>
<div class="line"><a name="l00613"></a><span class="lineno"> 613</span>  <span class="comment">// And we're done, but note that we haven't populated the Process entry with</span></div>
<div class="line"><a name="l00614"></a><span class="lineno"> 614</span>  <span class="comment">// the PDF or norm objects we created, as we assume the user has more work to</span></div>
<div class="line"><a name="l00615"></a><span class="lineno"> 615</span>  <span class="comment">// do on their RooWorkspace before copying into the CH instance. Once they've</span></div>
<div class="line"><a name="l00616"></a><span class="lineno"> 616</span>  <span class="comment">// imported the workspace:</span></div>
<div class="line"><a name="l00617"></a><span class="lineno"> 617</span>  <span class="comment">//</span></div>
<div class="line"><a name="l00618"></a><span class="lineno"> 618</span>  <span class="comment">// cb.AddWorkspace(ws);</span></div>
<div class="line"><a name="l00619"></a><span class="lineno"> 619</span>  <span class="comment">//</span></div>
<div class="line"><a name="l00620"></a><span class="lineno"> 620</span>  <span class="comment">// They can populate the Process entries, e.g. if all signals were replaced</span></div>
<div class="line"><a name="l00621"></a><span class="lineno"> 621</span>  <span class="comment">// with Morphing PDFs:</span></div>
<div class="line"><a name="l00622"></a><span class="lineno"> 622</span>  <span class="comment">// </span></div>
<div class="line"><a name="l00623"></a><span class="lineno"> 623</span>  <span class="comment">// cb.cp().signals().ExtractPdfs(cb, "htt", "$BIN_$PROCESS_morph");</span></div>
<div class="line"><a name="l00624"></a><span class="lineno"> 624</span>  <span class="keywordflow">return</span> std::string(morph_name);</div>
<div class="line"><a name="l00625"></a><span class="lineno"> 625</span> }</div>
<div class="line"><a name="l00626"></a><span class="lineno"> 626</span>  </div>
<div class="line"><a name="l00627"></a><span class="lineno"><a class="line" href="namespacech.html#a21ce0ae95c56bde7e1ecc1f21a839d2d"> 627</a></span> <span class="keywordtype">void</span> <a class="code" href="namespacech.html#a21ce0ae95c56bde7e1ecc1f21a839d2d">MakeMorphDebugPlots</a>(RooMorphingPdf* pdf, RooAbsReal* mass,</div>
<div class="line"><a name="l00628"></a><span class="lineno"> 628</span>  std::vector<double> <span class="keyword">const</span>& masses, TFile* f, TH1 *ref_bins) {</div>
<div class="line"><a name="l00629"></a><span class="lineno"> 629</span>  RooRealVar *rrv = <span class="keyword">dynamic_cast<</span>RooRealVar*<span class="keyword">></span>(mass);</div>
<div class="line"><a name="l00630"></a><span class="lineno"> 630</span>  <span class="keywordflow">if</span> (!rrv) <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00631"></a><span class="lineno"> 631</span>  f->cd();</div>
<div class="line"><a name="l00632"></a><span class="lineno"> 632</span>  f->mkdir(pdf->GetName());</div>
<div class="line"><a name="l00633"></a><span class="lineno"> 633</span>  gDirectory->cd(pdf->GetName());</div>
<div class="line"><a name="l00634"></a><span class="lineno"> 634</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> i = 0; i < masses.size(); ++i) {</div>
<div class="line"><a name="l00635"></a><span class="lineno"> 635</span>  rrv->setVal(masses[i]);</div>
<div class="line"><a name="l00636"></a><span class="lineno"> 636</span>  TH1 * h = pdf->createHistogram(<span class="stringliteral">"CMS_th1x"</span>);</div>
<div class="line"><a name="l00637"></a><span class="lineno"> 637</span>  h->AddDirectory(<span class="keyword">false</span>);</div>
<div class="line"><a name="l00638"></a><span class="lineno"> 638</span>  TH1 * h2 = <span class="keyword">nullptr</span>;</div>
<div class="line"><a name="l00639"></a><span class="lineno"> 639</span>  <span class="keywordflow">if</span> (ref_bins) { </div>
<div class="line"><a name="l00640"></a><span class="lineno"> 640</span>  h2 = (TH1*)ref_bins->Clone();</div>
<div class="line"><a name="l00641"></a><span class="lineno"> 641</span>  h2->Reset();</div>
<div class="line"><a name="l00642"></a><span class="lineno"> 642</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> x = 1; x <= h2->GetNbinsX(); ++x) {</div>
<div class="line"><a name="l00643"></a><span class="lineno"> 643</span>  h2->SetBinContent(x, h->GetBinContent(x));</div>
<div class="line"><a name="l00644"></a><span class="lineno"> 644</span>  }</div>
<div class="line"><a name="l00645"></a><span class="lineno"> 645</span>  } <span class="keywordflow">else</span> {</div>
<div class="line"><a name="l00646"></a><span class="lineno"> 646</span>  h2 = h;</div>
<div class="line"><a name="l00647"></a><span class="lineno"> 647</span>  h = <span class="keyword">nullptr</span>;</div>
<div class="line"><a name="l00648"></a><span class="lineno"> 648</span>  }</div>
<div class="line"><a name="l00649"></a><span class="lineno"> 649</span>  h2->AddDirectory(<span class="keyword">false</span>);</div>
<div class="line"><a name="l00650"></a><span class="lineno"> 650</span>  h2->SetName(TString::Format(<span class="stringliteral">"actual_point_%g"</span>, masses[i]));</div>
<div class="line"><a name="l00651"></a><span class="lineno"> 651</span>  gDirectory->WriteTObject(h2);</div>
<div class="line"><a name="l00652"></a><span class="lineno"> 652</span>  <span class="keywordflow">if</span> (h) <span class="keyword">delete</span> h;</div>
<div class="line"><a name="l00653"></a><span class="lineno"> 653</span>  <span class="keywordflow">if</span> (h2) <span class="keyword">delete</span> h2;</div>
<div class="line"><a name="l00654"></a><span class="lineno"> 654</span>  }</div>
<div class="line"><a name="l00655"></a><span class="lineno"> 655</span>  <span class="keywordtype">double</span> m = masses.front();</div>
<div class="line"><a name="l00656"></a><span class="lineno"> 656</span>  <span class="keywordtype">double</span> step = 1;</div>
<div class="line"><a name="l00657"></a><span class="lineno"> 657</span>  <span class="keywordflow">if</span> (((masses.back() - masses.front()) / step) > 100.) step = step * 10.;</div>
<div class="line"><a name="l00658"></a><span class="lineno"> 658</span>  <span class="keywordflow">if</span> (((masses.back() - masses.front()) / step) > 100.) step = step * 10.;</div>
<div class="line"><a name="l00659"></a><span class="lineno"> 659</span>  <span class="keywordflow">if</span> (((masses.back() - masses.front()) / step) < 10.) step = step/10.;</div>
<div class="line"><a name="l00660"></a><span class="lineno"> 660</span>  <span class="keywordflow">while</span> (m <= masses.back()) {</div>
<div class="line"><a name="l00661"></a><span class="lineno"> 661</span>  rrv->setVal(m);</div>
<div class="line"><a name="l00662"></a><span class="lineno"> 662</span>  TH1* hm = pdf->createHistogram(<span class="stringliteral">"CMS_th1x"</span>);</div>
<div class="line"><a name="l00663"></a><span class="lineno"> 663</span>  hm->AddDirectory(<span class="keyword">false</span>);</div>
<div class="line"><a name="l00664"></a><span class="lineno"> 664</span>  TH1 * hm2 = <span class="keyword">nullptr</span>;</div>
<div class="line"><a name="l00665"></a><span class="lineno"> 665</span>  <span class="keywordflow">if</span> (ref_bins) { </div>
<div class="line"><a name="l00666"></a><span class="lineno"> 666</span>  hm2 = (TH1*)ref_bins->Clone();</div>
<div class="line"><a name="l00667"></a><span class="lineno"> 667</span>  hm2->Reset();</div>
<div class="line"><a name="l00668"></a><span class="lineno"> 668</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> x = 1; x <= hm2->GetNbinsX(); ++x) {</div>
<div class="line"><a name="l00669"></a><span class="lineno"> 669</span>  hm2->SetBinContent(x, hm->GetBinContent(x));</div>
<div class="line"><a name="l00670"></a><span class="lineno"> 670</span>  }</div>
<div class="line"><a name="l00671"></a><span class="lineno"> 671</span>  } <span class="keywordflow">else</span> {</div>
<div class="line"><a name="l00672"></a><span class="lineno"> 672</span>  hm2 = hm;</div>
<div class="line"><a name="l00673"></a><span class="lineno"> 673</span>  hm = <span class="keyword">nullptr</span>;</div>
<div class="line"><a name="l00674"></a><span class="lineno"> 674</span>  }</div>
<div class="line"><a name="l00675"></a><span class="lineno"> 675</span>  hm2->AddDirectory(<span class="keyword">false</span>);</div>
<div class="line"><a name="l00676"></a><span class="lineno"> 676</span>  hm2->SetName(TString::Format(<span class="stringliteral">"morph_point_%g"</span>, m));</div>
<div class="line"><a name="l00677"></a><span class="lineno"> 677</span>  <span class="comment">//It struggles with m=0, instead adds really small number close to 0</span></div>
<div class="line"><a name="l00678"></a><span class="lineno"> 678</span>  <span class="keywordflow">if</span>(fabs(m) < 1E-5) hm2->SetName(TString::Format(<span class="stringliteral">"morph_point_0"</span>));</div>
<div class="line"><a name="l00679"></a><span class="lineno"> 679</span>  gDirectory->WriteTObject(hm2);</div>
<div class="line"><a name="l00680"></a><span class="lineno"> 680</span>  <span class="keywordflow">if</span> (hm) <span class="keyword">delete</span> hm;</div>
<div class="line"><a name="l00681"></a><span class="lineno"> 681</span>  <span class="keywordflow">if</span> (hm2) <span class="keyword">delete</span> hm2;</div>
<div class="line"><a name="l00682"></a><span class="lineno"> 682</span>  m += step;</div>
<div class="line"><a name="l00683"></a><span class="lineno"> 683</span>  }</div>
<div class="line"><a name="l00684"></a><span class="lineno"> 684</span>  f->cd();</div>
<div class="line"><a name="l00685"></a><span class="lineno"> 685</span> }</div>
<div class="line"><a name="l00686"></a><span class="lineno"> 686</span> }</div>
<div class="ttc" id="a_logging_8h_html"><div class="ttname"><a href="_logging_8h.html">Logging.h</a></div></div>
<div class="ttc" id="a_logging_8h_html_ab140ae8826e1f99c5aa0e920a9c13112"><div class="ttname"><a href="_logging_8h.html#ab140ae8826e1f99c5aa0e920a9c13112">FNERROR</a></div><div class="ttdeci">#define FNERROR(x)</div><div class="ttdef"><b>Definition:</b> <a href="_logging_8h_source.html#l00009">Logging.h:9</a></div></div>
<div class="ttc" id="a_morph_functions_8h_html"><div class="ttname"><a href="_morph_functions_8h.html">MorphFunctions.h</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html"><div class="ttname"><a href="classch_1_1_combine_harvester.html">ch::CombineHarvester</a></div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester_8h_source.html#l00030">CombineHarvester.h:30</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_a0d3eca727963c0c14a3c04bbdce2ca5e"><div class="ttname"><a href="classch_1_1_combine_harvester.html#a0d3eca727963c0c14a3c04bbdce2ca5e">ch::CombineHarvester::GetObservedShape</a></div><div class="ttdeci">TH1F GetObservedShape()</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester___evaluate_8cc_source.html#l00443">CombineHarvester_Evaluate.cc:443</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_a25c231ed8b1614c5d50f192210374474"><div class="ttname"><a href="classch_1_1_combine_harvester.html#a25c231ed8b1614c5d50f192210374474">ch::CombineHarvester::mass</a></div><div class="ttdeci">CombineHarvester & mass(std::vector< std::string > const &vec, bool cond=true)</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester___filters_8cc_source.html#l00097">CombineHarvester_Filters.cc:97</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_a3b73ec07ca4de14bf76383b2921f344a"><div class="ttname"><a href="classch_1_1_combine_harvester.html#a3b73ec07ca4de14bf76383b2921f344a">ch::CombineHarvester::ForEachSyst</a></div><div class="ttdeci">void ForEachSyst(Function func)</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester_8h_source.html#l00646">CombineHarvester.h:646</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_a3d265eb50ea8ec70192788d985733d9d"><div class="ttname"><a href="classch_1_1_combine_harvester.html#a3d265eb50ea8ec70192788d985733d9d">ch::CombineHarvester::GetShape</a></div><div class="ttdeci">TH1F GetShape()</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester___evaluate_8cc_source.html#l00287">CombineHarvester_Evaluate.cc:287</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_a44064524ffe47ccac63ef8025a6f46dc"><div class="ttname"><a href="classch_1_1_combine_harvester.html#a44064524ffe47ccac63ef8025a6f46dc">ch::CombineHarvester::ForEachProc</a></div><div class="ttdeci">void ForEachProc(Function func)</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester_8h_source.html#l00636">CombineHarvester.h:636</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_a7a0c6913b055e11798269034ae324298"><div class="ttname"><a href="classch_1_1_combine_harvester.html#a7a0c6913b055e11798269034ae324298">ch::CombineHarvester::SetFromProcs</a></div><div class="ttdeci">std::set< R > SetFromProcs(T func)</div><div class="ttdoc">Fill an std::set using only the Process entries.</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester_8h_source.html#l00615">CombineHarvester.h:615</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_a83dd50569a1f1a4b5a46f2ad2b5c1cdb"><div class="ttname"><a href="classch_1_1_combine_harvester.html#a83dd50569a1f1a4b5a46f2ad2b5c1cdb">ch::CombineHarvester::bin</a></div><div class="ttdeci">CombineHarvester & bin(std::vector< std::string > const &vec, bool cond=true)</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester___filters_8cc_source.html#l00013">CombineHarvester_Filters.cc:13</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_a91821b9464ebaacb5b92be2b3bc794f8"><div class="ttname"><a href="classch_1_1_combine_harvester.html#a91821b9464ebaacb5b92be2b3bc794f8">ch::CombineHarvester::FilterSysts</a></div><div class="ttdeci">CombineHarvester & FilterSysts(Function func)</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester_8h_source.html#l00674">CombineHarvester.h:674</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_aab6a698c192aa091178bbd8dfbd0266f"><div class="ttname"><a href="classch_1_1_combine_harvester.html#aab6a698c192aa091178bbd8dfbd0266f">ch::CombineHarvester::FilterProcs</a></div><div class="ttdeci">CombineHarvester & FilterProcs(Function func)</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester_8h_source.html#l00667">CombineHarvester.h:667</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_ac3ef7424cb86d12fd294d9c6de98d533"><div class="ttname"><a href="classch_1_1_combine_harvester.html#ac3ef7424cb86d12fd294d9c6de98d533">ch::CombineHarvester::syst_name_set</a></div><div class="ttdeci">std::set< std::string > syst_name_set()</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester___filters_8cc_source.html#l00271">CombineHarvester_Filters.cc:271</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_ade0f7c161d078a189404a274b93336d2"><div class="ttname"><a href="classch_1_1_combine_harvester.html#ade0f7c161d078a189404a274b93336d2">ch::CombineHarvester::cp</a></div><div class="ttdeci">CombineHarvester cp()</div><div class="ttdoc">Creates and returns a shallow copy of the CombineHarvester instance.</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester_8cc_source.html#l00220">CombineHarvester.cc:220</a></div></div>
<div class="ttc" id="aclassch_1_1_combine_harvester_html_af13ff7dd8c281cd9e9cfd6b11fc05b95"><div class="ttname"><a href="classch_1_1_combine_harvester.html#af13ff7dd8c281cd9e9cfd6b11fc05b95">ch::CombineHarvester::syst_type</a></div><div class="ttdeci">CombineHarvester & syst_type(std::vector< std::string > const &vec, bool cond=true)</div><div class="ttdef"><b>Definition:</b> <a href="_combine_harvester___filters_8cc_source.html#l00136">CombineHarvester_Filters.cc:136</a></div></div>
<div class="ttc" id="aclassch_1_1_object_html_a740d853956df951a7cb08094e21f0332"><div class="ttname"><a href="classch_1_1_object.html#a740d853956df951a7cb08094e21f0332">ch::Object::process</a></div><div class="ttdeci">virtual std::string const & process() const</div><div class="ttdef"><b>Definition:</b> <a href="_object_8h_source.html#l00020">Object.h:20</a></div></div>
<div class="ttc" id="aclassch_1_1_object_html_a7b588d78c22464ebe2b793cbaa2366d3"><div class="ttname"><a href="classch_1_1_object.html#a7b588d78c22464ebe2b793cbaa2366d3">ch::Object::bin</a></div><div class="ttdeci">virtual std::string const & bin() const</div><div class="ttdef"><b>Definition:</b> <a href="_object_8h_source.html#l00017">Object.h:17</a></div></div>
<div class="ttc" id="aclassch_1_1_object_html_ae2f9e80af391e105a6b6bc27107e3c9e"><div class="ttname"><a href="classch_1_1_object.html#ae2f9e80af391e105a6b6bc27107e3c9e">ch::Object::mass</a></div><div class="ttdeci">virtual std::string const & mass() const</div><div class="ttdef"><b>Definition:</b> <a href="_object_8h_source.html#l00038">Object.h:38</a></div></div>
<div class="ttc" id="aclassch_1_1_process_html"><div class="ttname"><a href="classch_1_1_process.html">ch::Process</a></div><div class="ttdef"><b>Definition:</b> <a href="_process_8h_source.html#l00015">Process.h:15</a></div></div>
<div class="ttc" id="aclassch_1_1_systematic_html"><div class="ttname"><a href="classch_1_1_systematic.html">ch::Systematic</a></div><div class="ttdef"><b>Definition:</b> <a href="_systematic_8h_source.html#l00013">Systematic.h:13</a></div></div>
<div class="ttc" id="aclassch_1_1_systematic_html_a1b40f502a9e5343cc1947bf452cd80d5"><div class="ttname"><a href="classch_1_1_systematic.html#a1b40f502a9e5343cc1947bf452cd80d5">ch::Systematic::value_u</a></div><div class="ttdeci">double value_u() const</div><div class="ttdef"><b>Definition:</b> <a href="_systematic_8h_source.html#l00028">Systematic.h:28</a></div></div>
<div class="ttc" id="aclassch_1_1_systematic_html_a2c91f65744afb7460bfffdf2a4987196"><div class="ttname"><a href="classch_1_1_systematic.html#a2c91f65744afb7460bfffdf2a4987196">ch::Systematic::type</a></div><div class="ttdeci">std::string const & type() const</div><div class="ttdef"><b>Definition:</b> <a href="_systematic_8h_source.html#l00025">Systematic.h:25</a></div></div>
<div class="ttc" id="aclassch_1_1_systematic_html_a5764d6db4954b8aac554776c149d3201"><div class="ttname"><a href="classch_1_1_systematic.html#a5764d6db4954b8aac554776c149d3201">ch::Systematic::asymm</a></div><div class="ttdeci">bool asymm() const</div><div class="ttdef"><b>Definition:</b> <a href="_systematic_8h_source.html#l00037">Systematic.h:37</a></div></div>
<div class="ttc" id="aclassch_1_1_systematic_html_a872469110fe349da45922431d6d24bc7"><div class="ttname"><a href="classch_1_1_systematic.html#a872469110fe349da45922431d6d24bc7">ch::Systematic::name</a></div><div class="ttdeci">std::string const & name() const</div><div class="ttdef"><b>Definition:</b> <a href="_systematic_8h_source.html#l00022">Systematic.h:22</a></div></div>
<div class="ttc" id="aclassch_1_1_systematic_html_ac0438bd51a527d0b2c1a70dd303ff4c3"><div class="ttname"><a href="classch_1_1_systematic.html#ac0438bd51a527d0b2c1a70dd303ff4c3">ch::Systematic::value_d</a></div><div class="ttdeci">double value_d() const</div><div class="ttdef"><b>Definition:</b> <a href="_systematic_8h_source.html#l00031">Systematic.h:31</a></div></div>
<div class="ttc" id="anamespacech_html"><div class="ttname"><a href="namespacech.html">ch</a></div><div class="ttdef"><b>Definition:</b> <a href="_algorithm_8h_source.html#l00010">Algorithm.h:10</a></div></div>
<div class="ttc" id="anamespacech_html_a011bcc8eb52bfa118fd44ade245e9e45"><div class="ttname"><a href="namespacech.html#a011bcc8eb52bfa118fd44ade245e9e45">ch::BuildRooMorphing</a></div><div class="ttdeci">std::string BuildRooMorphing(RooWorkspace &ws, CombineHarvester &cb, std::string const &bin, std::string const &process, RooAbsReal &mass_var, std::string norm_postfix, bool allow_morph, bool verbose, bool force_template_limit=false, TFile *file=nullptr)</div><div class="ttdoc">[part1]</div><div class="ttdef"><b>Definition:</b> <a href="_morph_functions_8cc_source.html#l00020">MorphFunctions.cc:20</a></div></div>
<div class="ttc" id="anamespacech_html_a21ce0ae95c56bde7e1ecc1f21a839d2d"><div class="ttname"><a href="namespacech.html#a21ce0ae95c56bde7e1ecc1f21a839d2d">ch::MakeMorphDebugPlots</a></div><div class="ttdeci">void MakeMorphDebugPlots(RooMorphingPdf *pdf, RooAbsReal *mass, std::vector< double > const &masses, TFile *f, TH1 *ref_bins=nullptr)</div><div class="ttdef"><b>Definition:</b> <a href="_morph_functions_8cc_source.html#l00627">MorphFunctions.cc:627</a></div></div>
<div class="ttc" id="anamespacech_html_a2cbcfd049543bfe9f63958bc671deeae"><div class="ttname"><a href="namespacech.html#a2cbcfd049543bfe9f63958bc671deeae">ch::RebinHist</a></div><div class="ttdeci">TH1F RebinHist(TH1F const &hist)</div><div class="ttdef"><b>Definition:</b> <a href="_utilities_8cc_source.html#l00170">Utilities.cc:170</a></div></div>
<div class="ttc" id="anamespacech_html_a5e23234ab65cfa48c99b83d94e291125"><div class="ttname"><a href="namespacech.html#a5e23234ab65cfa48c99b83d94e291125">ch::Set2Vec</a></div><div class="ttdeci">std::vector< T > Set2Vec(std::set< T > const &in)</div><div class="ttdef"><b>Definition:</b> <a href="_utilities_8h_source.html#l00101">Utilities.h:101</a></div></div>
<div class="ttc" id="anamespacech_html_af9b462e52384bf17acf24c57bee7a5d3"><div class="ttname"><a href="namespacech.html#af9b462e52384bf17acf24c57bee7a5d3">ch::AsTH1F</a></div><div class="ttdeci">TH1F AsTH1F(TH1 const *hist)</div><div class="ttdef"><b>Definition:</b> <a href="_morph_functions_8h_source.html#l00020">MorphFunctions.h:20</a></div></div>
<div class="ttc" id="astructch_1_1syst_1_1process_html"><div class="ttname"><a href="structch_1_1syst_1_1process.html">ch::syst::process</a></div><div class="ttdef"><b>Definition:</b> <a href="_systematics_8h_source.html#l00036">Systematics.h:36</a></div></div>
</div><!-- fragment --></div><!-- contents -->
</div><!-- doc-content -->
<!-- HTML footer for doxygen 1.8.6-->
<!-- start footer part -->
</body>
</html>