LAT/DataSetInfo.hh at master · wisecg/LAT · GitHub

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
#ifndef __DataSetInfo_hh__
#define __DataSetInfo_hh__

#include <iostream>
#include <map>
#include <cmath>
#include "glob.h"

#include "GATDataSet.hh"

using namespace std;

// ======================================================================
// Contents:

// FindDataSet - Returns the DS number of a given run.
// GetLNRunCoverage - Given a run and DS number, verify that this run is covered by the LN fill tag

// LoadDataSet - Contains run ranges for all DS's, updates a GATDataSet.
// LoadBlindDataSet - Contains run ranges for all BLIND DS's, updates a GATDataSet.

// GetDataSetSequences - Returns a map of all the sub-ranges in a DS.
// GetBlindDataSetSequences - Returns a map of all the sub-ranges in a DS.

// GetDSStartTime - find the first unix timestamp of the first run in the DS.
// GetDSRunAndStartTimes - DEPRECATED.  Start time and run time of each DS (as of 0nbb PRL)

// LoadDetectorList - Gives a list of the detector names for each module.
// GetTotalActiveMass - Total active mass for each dataset.
//                      This could be calculated instead of hardcoded in the future.
// LoadActiveMasses - Returns a map of all active masses.
// LoadActiveMassUncertainties - Returns a map of all active mass uncertainties.

// LoadBadDetectorMap - Returns a map of bad (i.e. not biased, unusuable) detectors.
// LoadVetoDetectorMap - Returns a map of veto-only detectors.

// GetChannelSelectionPath - Returns a string with the path to the highest
//                          version of the channel selection files.
// LoadEnrNatMap - quick way to tell if a given detID is enriched (1) or natural (0).
// CheckModule - Given a detector ID, look up which module it lives in.
// GetVetoActiveMass - Modifies total mass to not include veto-only detectors.
// GetENFC - Parameters for corrected trapENFCal
// GetAvsE - AvsE parameters
// GetDCR* - DCR parameters
// LoadDS4MuonList - Static muon list for DS-4, calculated manually
//                   by $GATDIR/mjd-veto/skim-veto.cc

// LoadLNFillTimes1 - Returns a vector of M1 LN fills.
// LoadLNFillTimes2 - Returns a vector of M2 LN fills.
// ======================================================================

int FindDataSet(int run);

void LoadDataSet(GATDataSet& ds, int dsNum, int subNum=-1);
void LoadDataSet(GATDataSet& ds, int dsNum, int subNum, map<int,vector<int>>& ranges);

void LoadBlindDataSet(GATDataSet& ds, int dsNum, int subNum=-1);
void LoadBlindDataSet(GATDataSet& ds, int dsNum, int subNum, map<int,vector<int>>& ranges);

int GetDataSetSequences(int dsNum);
int GetBlindDataSetSequences(int dsNum);

double GetDSStartTime(int dsNum, bool blindMode=0);
void GetDSRunAndStartTimes(int dsNum, double &runTime_s, double &startTime0);

vector<int> LoadDetectorList(int module);
void GetTotalActiveMass(int dsNum, double& m1Total, double& m1Enr, double& m1Nat, double& m2Total, double& m2Enr, double& m2Nat);
map<int,double> LoadActiveMasses(int dsNum);
map<int,double> LoadActiveMassUncertainties(int dsNum);

map<int,bool> LoadBadDetectorMap(int dsNum);
map<int,bool> LoadVetoDetectorMap(int dsNum);

std::string GetChannelSelectionPath(int dsNum, int officialVersion = -1);
map<int, bool> LoadEnrNatMap();
int CheckModule(int detID);
void GetVetoActiveMass(map<int,double> actM4Det_g, map<int,bool> detIDIsVetoOnly,
  double &mVeto_M1Total_kg, double &mVeto_M2Total_kg);

double GetENFC(int chan, int dsNum, double trapENF, int run);
double GetAvsE(int chan, double TSCurrent50nsMax, double TSCurrent100nsMax, double TSCurrent200nsMax,
  double trapENF, double trapENFCal, int dsNum, int run);

double GetDCR90(int chan, double nlcblrwfSlope, double trapMax, int dsNum, int run);
double GetDCRCTC90(int chan, double nlcblrwfSlope, double trapE, double trapMax, int dsNum);
double GetDCR85(int chan, double nlcblrwfSlope, double trapMax, int dsNum, int run);
double GetDCR95(int chan, double nlcblrwfSlope, double trapMax, int dsNum, int run);
double GetDCR98(int chan, double nlcblrwfSlope, double trapMax, int dsNum, int run);
double GetDCR99(int chan, double nlcblrwfSlope, double trapMax, int dsNum, int run);
double GetDCR995(int chan, double nlcblrwfSlope, double trapMax, int dsNum, int run);
double GetDCR999(int chan, double nlcblrwfSlope, double trapMax, int dsNum, int run);

bool GetLNRunCoverage(int dsNum, int run);
void LoadDS4MuonList(vector<int> &muRuns, vector<double> &muRunTStarts,
  vector<double> &muTimes,vector<int> &muTypes, vector<double> &muUncert);
void LoadLNFillTimes1(vector<double>& lnFillTimes1, int dsNum);
void LoadLNFillTimes2(vector<double>& lnFillTimes2, int dsNum);

#endif