Exposures.h
Go to the documentation of this file.
1 #pragma once
2 
3 // \date: 07/24/2015
4 
5 namespace ana
6 {
7  // A collection of Exposure-related constants
8 
9  // Total POT/Exposure/Livetime for combined pre-/post-shutdown and top-off
10  // FD NuMI data. Doc-db 13556 is the reference for all info herein for the First Analysis.
11  // Date format: [mm/dd/yyyy]
12 
13  // Pre-shutdown: run 12942 to 17132 (02/06/2014 to 09/05/2014)
14  // Post-shutdown: run 17893 through run 19096 (10/24/2014 to 03/14/2015)
15  // Top-off: run 19097 to run 19587 (03/14/2015 to 05/15/2015)
16 
17  // POT
18  const double kFirstAnaPrePOT = 1.20e+20;
19  const double kFirstAnaPostPOT = 1.34e+20;
20  const double kFirstAnaTopOffPOT = 0.91e+20;
21 
22  // Total POT used for First Analysis: Correct as of 07/24/2015
23  const double kFirstAnaPOT = kFirstAnaPrePOT +
24  kFirstAnaPostPOT +
26 
27  // Exposure [POT x kT]
28  const double kFirstAnaPreExposure = 5.7e+20;
29  const double kFirstAnaPostExposure = 12.6e+20;
30  const double kFirstAnaTopOffExposure = 9.0e+20;
31 
32  const double kFirstAnaExposure = kFirstAnaPreExposure +
33  kFirstAnaPostExposure +
35 
36  const double kFirstAnaFullDetEquivalentPOT = 2.74e+20;
37 
38  // Livetime [s]
39  const double kFirstAnaPreLivetime = 125.1;
40  const double kFirstAnaPostLivetime = 67.5;
41  const double kFirstAnaTopOffLivetime = 35.0;
42 
43  const double kFirstAnaLivetime = kFirstAnaPreLivetime +
44  kFirstAnaPostLivetime +
46 
47  // Effective fiducial mass [kT] for each diblock
48  const double kFirstAnaFiducialMass[15] = {0.,
49  0.,
50  0.,
51  0.,
52  2.3,
53  3.0,
54  3.8,
55  4.6,
56  5.4,
57  6.1,
58  6.9,
59  7.7,
60  8.5,
61  9.2,
62  10.0};
63 
64 ////////////////////////////////////Second Ana//////////////////////////////////
65 
66  // Second Analysis POT and Livetime accounting
67  // Current accounting from docdb-15390
68  // Period 1: 12942 - 17139
69  // Period 2: 17891 - 19746
70  // Epoch 3b: 20923 - 21230
71  // Epoch 3c: 21231 - 22019
72  // Top off Epoch 3d: 22020 - 22900
73  const double kSecondAnaPeriod1POT = 1.23e+20;
74  const double kSecondAnaPeriod2POT = 2.64e+20;
75  const double kSecondAnaEpoch3bPOT = 0.16e+20;
76  const double kSecondAnaEpoch3cPOT = 0.77e+20;
77  const double kSecondAnaEpoch3dPOT = 1.89e+20;
78  const double kSecondAnaPeriod3POT = kSecondAnaEpoch3bPOT + kSecondAnaEpoch3cPOT + kSecondAnaEpoch3dPOT;
79 
80  //Run matched MC exists for periods 1 to 3c. Keep track of POT for this
81  //period. In the end we will scale up epoch3c to cover to top-off
82  const double kSecondAnaPeriod1to3cPOT = kSecondAnaPeriod1POT +
83  kSecondAnaPeriod2POT +
84  kSecondAnaEpoch3bPOT +
86 
87  const double kSecondAnaPOT = kSecondAnaPeriod1POT +
88  kSecondAnaPeriod2POT +
90 
91  //Exposure numbers [kT x POT] will be added once cosmic rejection is set,
92  //which fixes fiducial mass
93 
94  // For current FD MC which covers Period1 to Epoch 3c this is the appropriate
95  // scale factor to be applied to the whole MC sample. docdb-15390
96  const double kSecondAnaRunMatchedMCPOTScale = 6.957e+20;
97 
98  //Livetime [s]
99  const double kSecondAnaPeriod1Livetime = 128.1;
100  const double kSecondAnaPeriod2Livetime = 119.1;
101  const double kSecondAnaEpoch3bLivetime = 8.7;
102  const double kSecondAnaEpoch3cLivetime = 37.0;
103  const double kSecondAnaEpoch3dLivetime = 60.7;
104  const double kSecondAnaPeriod3Livetime = kSecondAnaEpoch3bLivetime +
105  kSecondAnaEpoch3cLivetime +
107 
108  const double kSecondAnaPeriod1to3cLivetime = kSecondAnaPeriod1Livetime +
109  kSecondAnaPeriod2Livetime +
110  kSecondAnaEpoch3bLivetime +
112 
113  const double kSecondAnaLivetime = kSecondAnaPeriod1Livetime +
114  kSecondAnaPeriod2Livetime +
116 
117 
118  // Effective fiducial mass [kT] for each diblock
119  const double kSecondAnaFiducialMass[15] = {0.,
120  0.,
121  0.,
122  0.,
123  2.5,
124  3.3,
125  4.1,
126  4.9,
127  5.6,
128  6.4,
129  7.2,
130  8.0,
131  8.8,
132  9.5,
133  10.3};
134 
135  //Count only 14 diblock POT
136  const double kOnly14DBPeriod1POT = 0.02e+20;
137  const double kOnly14DBPeriod2POT = 2.41e+20;
138  const double kOnly14DBEpoch3bPOT = 0.15e+20;
139  const double kOnly14DBEpoch3cPOT = 0.76e+20;
140  const double kOnly14DBEpoch3dPOT = 1.88e+20;
141 
142  const double kOnly14DBPeriod1to3cPOT = kOnly14DBPeriod1POT +
143  kOnly14DBPeriod2POT +
144  kOnly14DBEpoch3bPOT +
146 
147  const double kOnly14DBPOT = kOnly14DBPeriod1to3cPOT + kOnly14DBEpoch3dPOT;
148 
149  const double kOnly14DBPeriod1Livetime = 2.3;
150  const double kOnly14DBPeriod2Livetime = 104.7;
151  const double kOnly14DBEpoch3bLivetime = 8.6;
152  const double kOnly14DBEpoch3cLivetime = 36.7;
153  const double kOnly14DBEpoch3dLivetime = 60.2;
154 
155  const double kOnly14DBPeriod1to3cLivetime = kOnly14DBPeriod1Livetime +
156  kOnly14DBPeriod2Livetime +
157  kOnly14DBEpoch3bLivetime +
159 
160  const double kOnly14DBLivetime = kOnly14DBPeriod1to3cLivetime + kOnly14DBEpoch3dLivetime;
161 
162  ////////////////////////////////////2017 Analysis//////////////////////////////////
163  // See Dan's talk - docdb 16455 for the idea behind the livetime estimate
164  const double kLivetimeSensitivity2017 = 440;
165  const double kPOTSensitivity2017 = 9e+20;
166 
167  //////////////////////////////////////////////////////////////////////////////
168  ////////// Updated 2017 Analysis Exposure
169  ////////// See Vlad's talk for details (DocDB 21140)
170  const double kAna2017Epoch3ePOT = 1.036e+20;
171  const double kAna2017Period5POT = 1.763e+20;
172  const double kAna2017Period3POT = kSecondAnaPeriod3POT + kAna2017Epoch3ePOT;
173 
174  const double kAna2017POT = kSecondAnaPeriod1POT +
175  kSecondAnaPeriod2POT +
176  kAna2017Period3POT +
178 
179  const double kAna2017FiducialMass[15] = {0.,
180  0.,
181  0.,
182  0.,
183  2.57,
184  3.36,
185  4.14,
186  4.93,
187  5.72,
188  6.50,
189  7.29,
190  8.07,
191  8.86,
192  9.64,
193  10.43};
194 
195  const double kAna2017FullDetEquivPOT = 8.85e+20;
196 
197  const double kAna2017Epoch3eLivetime = 32.6;
198  const double kAna2017Period5Livetime = 51.9;
199  const double kAna2017Period3Livetime = kSecondAnaPeriod3Livetime + kAna2017Epoch3eLivetime;
200  const double kAna2017Livetime = kSecondAnaPeriod1Livetime +
201  kSecondAnaPeriod2Livetime +
202  kAna2017Period3Livetime +
204 
205  ///////////////////////////////////////2018 Analysis/////////////////////////////////////
206 
207  const double kAna2018FHCPOT = 9.47964e20;
208  const double kAna2018RHCPOT = 6.90562e20;
209 
210  const double kAna2018SensitivityFHCNDPOT = 8.04e20;
211  const double kAna2018SensitivityRHCNDPOT = 3.10e20;
212 
213  const double kAna2018FHCLivetime = 438.1584;
214  const double kAna2018RHCLivetime = 179.6173;
215 
216  ///////////////////////////////////////2019 Analysis/////////////////////////////////////
217  // For future reference, if a sum is taken for the full FHC Livetime histogram one would
218  // find 438.1582 for the FHC livetime. This is more proper than the sum of individual
219  // periods with somewhat different decimal places reported. However, the change is
220  // very very tiny, so we are leaving it at the 2018 value and doing the proper thing
221  // for RHC. See doc-37260
222 
223  const double kAna2019FHCPOT = 9.47964e20;
224  const double kAna2019RHCPOT = 12.3319e20;
225 
226  const double kAna2019FHCLivetime = 438.1584;
227  const double kAna2019RHCLivetime = 317.0466;
228 
229  ///////////////////////////////////////2020 Analysis/////////////////////////////////////
230  // see doc-45139
231 
232 
233  const double kAna2020FHCPOT = 14.2283e20;
234  const double kAna2020FHCFullDetEquivPOT = 13.6e+20;
235  const double kAna2020RHCPOT = 12.5003e20;
236 
237  const double kAna2020FHCLivetime = 555.3415;
238  const double kAna2020RHCLivetime = 321.1179;
239 
240  /////////////////////////////////////// Future Analysis /////////////////////////////////
241  // ping me if you strongly against any of these values
242  // Karl: New estimates are 31.5e20, not 36e20. Scaling livetimes appropriately from
243  // 1600s to 1400s.
244  const double kFutureFHCPOT = 31.5e20;
245  const double kFutureRHCPOT = 31.5e20;
246 
247  const double kFutureFHCLivetime = 1400;
248  const double kFutureRHCLivetime = 1400;
249 
250 
251  // ===== FHC =====
252  // ----- POT -----
253  // period 1 1.22728E+20
254  // period 2 2.64473E+20
255  // period 3 3.84544E+20
256  // period 5 1.76219E+20
257  // ----- Livetime -----
258  // period 1 128.116
259  // period 2 1.19E+02
260  // period 3 138.936
261  // period 5 51.9274
262 
263 
264  // ===== RHC =====
265  // ----- POT -----
266  // period 4 6.55380E+19
267  // period 6 3.06065E+20
268  // epoch 7a 7.78044E+19
269  // epoch 7b 1.85763E+20
270  // epoch 7c 5.53911E+19
271  // ----- Livetime -----
272  // period 4 19.577
273  // period 6 76.2761
274  // epoch 7a 24.7618
275  // epoch 7b 45.7463
276  // epoch 7c 13.2561
277 
278 
279 }
const double kFirstAnaTopOffExposure
Definition: Exposures.h:30
const double kLivetimeSensitivity2017
Definition: Exposures.h:164
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
const double kFirstAnaPreLivetime
Definition: Exposures.h:39
const double kFirstAnaTopOffLivetime
Definition: Exposures.h:41
const double kSecondAnaPeriod2POT
Definition: Exposures.h:74
const double kAna2019RHCLivetime
Definition: Exposures.h:227
const double kOnly14DBPeriod1POT
Definition: Exposures.h:136
const double kSecondAnaPeriod1to3cPOT
Definition: Exposures.h:82
const double kSecondAnaEpoch3dPOT
Definition: Exposures.h:77
const double kAna2017Period5POT
Definition: Exposures.h:171
const double kFirstAnaFiducialMass[15]
Definition: Exposures.h:48
const double kFirstAnaLivetime
Definition: Exposures.h:43
const double kSecondAnaEpoch3bPOT
Definition: Exposures.h:75
const double kFutureRHCPOT
Definition: Exposures.h:245
const double kAna2017POT
Definition: Exposures.h:174
const double kFutureFHCPOT
Definition: Exposures.h:244
const double kAna2018SensitivityFHCNDPOT
Definition: Exposures.h:210
const double kSecondAnaPeriod1POT
Definition: Exposures.h:73
const double kOnly14DBPeriod2POT
Definition: Exposures.h:137
const double kFirstAnaPostExposure
Definition: Exposures.h:29
const double kAna2017Period3Livetime
Definition: Exposures.h:199
const double kAna2020FHCFullDetEquivPOT
Definition: Exposures.h:234
const double kSecondAnaPeriod1to3cLivetime
Definition: Exposures.h:108
const double kAna2018RHCPOT
Definition: Exposures.h:208
const double kOnly14DBEpoch3cPOT
Definition: Exposures.h:139
const double kSecondAnaPeriod3Livetime
Definition: Exposures.h:104
const double kFirstAnaFullDetEquivalentPOT
Definition: Exposures.h:36
const double kAna2017Livetime
Definition: Exposures.h:200
const double kAna2017Epoch3eLivetime
Definition: Exposures.h:197
const double kFutureFHCLivetime
Definition: Exposures.h:247
const double kFirstAnaPreExposure
Definition: Exposures.h:28
const double kOnly14DBEpoch3dPOT
Definition: Exposures.h:140
const double kOnly14DBPOT
Definition: Exposures.h:147
const double kAna2020FHCLivetime
Definition: Exposures.h:237
const double kFirstAnaTopOffPOT
Definition: Exposures.h:20
const double kOnly14DBEpoch3bLivetime
Definition: Exposures.h:151
const double kOnly14DBEpoch3dLivetime
Definition: Exposures.h:153
const double kOnly14DBPeriod2Livetime
Definition: Exposures.h:150
const double kPOTSensitivity2017
Definition: Exposures.h:165
const double kSecondAnaEpoch3cLivetime
Definition: Exposures.h:102
const double kAna2019FHCLivetime
Definition: Exposures.h:226
const double kAna2017Period5Livetime
Definition: Exposures.h:198
const double kAna2020FHCPOT
Definition: Exposures.h:233
const double kFirstAnaPostLivetime
Definition: Exposures.h:40
const double kFutureRHCLivetime
Definition: Exposures.h:248
const double kAna2020RHCPOT
Definition: Exposures.h:235
const double kAna2019RHCPOT
Definition: Exposures.h:224
const double kOnly14DBPeriod1to3cPOT
Definition: Exposures.h:142
const double kSecondAnaPeriod1Livetime
Definition: Exposures.h:99
const double kSecondAnaEpoch3dLivetime
Definition: Exposures.h:103
const double kAna2017FullDetEquivPOT
Definition: Exposures.h:195
const double kOnly14DBEpoch3cLivetime
Definition: Exposures.h:152
const double kFirstAnaPOT
Definition: Exposures.h:23
const double kAna2020RHCLivetime
Definition: Exposures.h:238
const double kAna2018SensitivityRHCNDPOT
Definition: Exposures.h:211
const double kOnly14DBPeriod1Livetime
Definition: Exposures.h:149
const double kSecondAnaRunMatchedMCPOTScale
Definition: Exposures.h:96
const double kAna2019FHCPOT
Definition: Exposures.h:223
const double kOnly14DBPeriod1to3cLivetime
Definition: Exposures.h:155
const double kOnly14DBLivetime
Definition: Exposures.h:160
const double kSecondAnaEpoch3cPOT
Definition: Exposures.h:76
const double kSecondAnaPeriod2Livetime
Definition: Exposures.h:100
const double kSecondAnaFiducialMass[15]
Definition: Exposures.h:119
const double kAna2018FHCPOT
Definition: Exposures.h:207
const double kOnly14DBEpoch3bPOT
Definition: Exposures.h:138
const double kSecondAnaEpoch3bLivetime
Definition: Exposures.h:101
Float_t e
Definition: plot.C:35
const double kSecondAnaPOT
Definition: Exposures.h:87
const double kAna2017Epoch3ePOT
Definition: Exposures.h:170
const double kFirstAnaPrePOT
Definition: Exposures.h:18
const double kAna2017Period3POT
Definition: Exposures.h:172
const double kFirstAnaPostPOT
Definition: Exposures.h:19
const double kSecondAnaPeriod3POT
Definition: Exposures.h:78
const double kAna2018FHCLivetime
Definition: Exposures.h:213
const double kAna2018RHCLivetime
Definition: Exposures.h:214
const double kAna2017FiducialMass[15]
Definition: Exposures.h:179
const double kSecondAnaLivetime
Definition: Exposures.h:113
const double kFirstAnaExposure
Definition: Exposures.h:32