Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
hcalControls.h
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file hcalControls.h
1 
2 #ifndef _HCALCONTROLS_H_
3 #define _HCALCONTROLS_H_
4 #include <Riostream.h>
5 #include <TString.h>
6 
7 
8 enum OBJECTKINDS {
9  CALOR,
10  COUNTER
11 };
12 enum GAINS {
13  HIGH,
14  LOW
15 };
16 enum STACKS{
17  EMC, // EMCal
18  HINNER, // Inner HCal
19  HOUTER, // Outer HCal
20  HODO, // Scintillating Hodoscope
21  SCINT, // Scintillators (Counters along the beam line)+Cherenkovs
22  CHER // CHerenkov counters
23 };
24 
25 // **************************************************************************
26 
27 // CURRENTLY OPERATING DAQ
28 static const int READOUTCHANNELS = 192; // four ADC stations (EMC/HCInner/HCOuter/Service)
29 
30 static const int TRIGGERCHANNELS = 0; // 4;
31 
32 // defaults
33 static const int chinspected = 32; // high and low gains in HCal Inner
34 static const int triggerinsp = 0;
35 
36 // HCAL 2016 EXTERNAL COUNTERS AND HODOSCOPES
37 
38 // T1044 CONFIGURATION
39 static const int CALSTACKS = 6;
40 static const int T1044TRIGGERCH = 14;
41 static const int T1044HODO = 16;
42 static const int NSAMPLES = 24;
43 static const Double_t PEDESTAL = 2048.;
44 // thresholds below are in units ADC in selected TRIGGER gain range
45 // all thresholds are in 8 adc counts steps
46 // mixed amplitudes are stored in LG units, so the actual trigger thresholds is computed as
47 // 8./(TRGAINRANGE==0? HLGRATIO : 1)*TWRAMPTHR
48 // 8./(TRGAINRANGE==0? HLGRATIO : 1)*STTOTAMPTHR
49 // the 8 corresponds to lowerst 3 bits suppressed by Chi before thresholding
50 static Int_t TRGAINRANGE = 0; // gain range used to trigger
51 static const Int_t STTOTAMPTHR = 10; // total sum threshold in daq.c
52 static const Int_t TWRAMPTHR = 1; // individual tower threshold in daq.c
53 static const Int_t STHITMIN = 3; // minimal number of hits above threshold in stack in daq.c
54 static const Int_t STHITMAX = 14;
55 // next two thresholds are computed based upon pedestal widths in towers and stacks
56 // Towers(RMS) ~ 1.2; Stacks (RMS) ~4.7 (sum, cleaned), 7.1 (Fit) (in units of LG)
57 static const Double_t TWRZEROSUPTHR = 1.; // tower zero suppression threshold (units of LG)
58 static const Double_t STZEROSUPTHR = 5.;
59 static const Double_t HLGRATIO = 16.;
60 static const Float_t hlgratios[] = { 1., 32., 16., 1., 1., 1.};
61 static Int_t emcGainSelection = 1;
62 // HCAL PROTOTYPE
63 static const int HCALTOWERS = 16;
64 static const int HCALROWS = 4;
65 static const int HCALCOLUMNS = 4;
66 static const int HCALGAINS = 2;
67 
68 // EMC 2016 PROTOTYPE
69 static const int EMCTOWERS = 64;
70 static const int EMCROWS = 8;
71 static const int EMCCOLUMNS = 8;
72 static const int EMCGAINS = 1;
73 
74 static int CHTOTAL = 0; // all channels stored into ROOT file
75 static int ACTIVECHANNELS = 0; // total number of mapped channels
76 static int detchannels = 0;
77 static int hgDetChannels = 0;
78 static int lgDetChannels = 0;
79 static const int ovrflow = 1990; // 2048
80 static const int undflow = -1990; //-2048
81 // CALORIMETER CALIBRATION
82 // GeV/ADC count (LG)
83 static const Double_t stECalib[] = {0.008, 0.010, 0.005, 1., 1., 1.};
84 // HINNER has HGain dropped by x2.
85 // This kind of Scaling to be used only to get uncalibrated sums of amplitudes
86 static Double_t stAScale[] = {1., 2., 1., 1., 1., 1.};
87 
88 // TILE
89 static const int TILEFIBERS = 4;
90 static const int TILECHANNELS = 2*TILEFIBERS;
91 
92 // TILE CONTROLS
93 //static const int NSAMPLES = 24;
94 static const int TILETRIGGERCH = 8;
95 static const int TRIGGERRES = 3;
96 static const int CHANNELTHRESHOLDS = 25;
97 static const int HITMULTTHRESHOLDS = 8;
98 static const int RISETIME = 4;
99 static const int FALLTIME = 5;
100 static const int NPARAMETERS = 6;
101 static const int DISPLAYX = 4; // maximum number of channels per display panel
102 // number of channels in detector (set by Helper, if left zero on the first call to hLabHelper
103 // - it will be set to ACTIVECHANNELS-TRIGGERCHANNELS
104 
105 // **************************************************************************
106 // all informative ADC (TILE) channels (stored into .root file)
107 static const int chInUse[] = { 115, 113, 119, 117, 123, 121, 127, 125 };
108 
109 // **************************************************************************
110 // DEFAULTS FOR INSPECTION AND PLOTTING
111 //static const int feechinsp[] = {64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,
112 // 80,81,82,83,84,85,86,87,88,89,90,91.92,93.94,95};
113 // Mapping for HINNER 3/2/2016
114 static const int feechinsp[] = {64,72,80,88,66,74,82,90,68,76,84,92,70,78,86,94,
115  65,73,81,89,67,75,83,91,69,77,85,93,71,79,87,95};
116 // **************************************************************************
117 // 2016 FNAL Beam Test Prototype
118 static const int hcalInnerCh[] = { 64, 72, 80, 88, 66, 74, 82, 90, 68, 76, 84, 92, 70, 78, 86, 94,
119  65, 73, 81, 89, 67, 75, 83, 91, 69, 77, 85, 93, 71, 79, 87, 95};
120 static const int hcalOuterCh[] = {112,120,128,136,114,122,130,138,116,124,132,140,118,126,134,142,
121  113,121,129,137,115,123,131,139,117,125,133,141,119,127,135,143};
122 // static const int hcalOuterCh[] = {16,24,32,40,18,26,34,42,20,28,36,44,22,30,38,46,
123 // 17,25,33,41,19,27,35,43,21,29,37,45,23,31,39,47};
124 
125 
126 
127 
128 
129 static const int emcCh[] = { 5, 4, 7, 6, 1, 0, 3, 2, 13, 12, 15, 14, 9, 8, 11, 10,
130  21, 20, 23, 22, 17, 16, 19, 18, 29, 28, 31, 30, 25, 24, 27, 26,
131  37, 36, 39, 38, 33, 32, 35, 34, 45, 44, 47, 46, 41, 40, 43, 42,
132  53, 52, 55, 54, 49, 48, 51, 50, 61, 60, 63, 62, 57, 56, 59, 58};
133 static const int hodoCh[] = { 96, 97, 98, 99,100,101,102,103,104,105,106,107,108,109,110,111};
134 static const int counters[] = { 144,145,146,147,148,149,150,
135  160,161,162,163,164,165,166}; // including cherenkov's
136 
137 // **************************************************************************
138 // SMART TILE VERSION 0
139 
140 static const Double_t tileSizeX = 25.;
141 static const Double_t tileSizeY = 15.;
142 static const Int_t minProjEntries = 50; // min # of hits to try fitting the light yield
143 // TILECHANNELS
144 static const int feech1[] = { 115, 119, 123, 127, 113, 117, 121, 125 }; // HG channels only. Valid for earlier runs
145 static const int feech2[] = { 115, 113, 119, 117, 123, 121, 127, 125 }; // HG channels only. Valid starting Run 803
146 // TRIGGERCHANNELS
147 //static const int trch[] = { 134, 135, 128, 129, 130, 131, 132, 133 };
148 
149 static const int trch1125[] = { 128, 129, 130, 131 }; // two of the four channels are used by trigger tags
150 static const int trch1152[] = { 98, 99, 129, 131 }; // first are two gain ranges for Tower then sc trigger tags
151 // Calibration data at different times
152 static const Double_t sc_779[] = { 3.6, 5.2, 5.2, 4.5, 6.6, 5.3, 7.9, 4.2 }; // single cell calibration - cosmic run 779
153 static const Double_t mu_779[] = { 72.7, 130.6, 141.5, 132.1, 132., 138., 135.5, 68.4 }; // cosmic calibration Run 779
154 
155 static const Double_t sc_900[] = { 6.6, 12.3, 12.9, 12.3, 14.6, 14.2, 12.25, 6.42 }; // cosmic run, low amp. signals
156 
157 static const Double_t sc_1061[] = { 6.71, 12.74, 13.13, 12.76, 14.30, 14.35, 12.70, 6.56 }; // cosmic run, low amp. signals, mode 1
158 //const Double_t sc_1061[] = { 6.65, 12.68, 12.86, 12.55, 14.25 , 14.10, 12.23, 6.59 }; // cosmic run, low amp. signals, mode 0
159 
160 static const Double_t sc_1123[] = { 6.61, 11.83, 12.77, 12.41, 13.75, 13.55, 12.31, 6.44 }; // cosmic run, low amp. signals, mode 0
161 // this is just an educated guess for the time being. Estimate includes Gain difference at the same bias.
162 
163 // **************************************************************************
164 // SMART TILE DATA
165 // data and ROOT files on RCF
166 //static const TString RCFdataDir("/direct/phenix+user04/sphenix/bnl/buffer/");
167 //static const TString RCFrootDir("/direct/phenix+user04/sphenix/bnl/buffer/rootf/");
168 // data and ROOT files on sphnxdaq
169 //static const TString HLABdataDir("/scratch/buffer/bnl/");
170 //static const TString HLABrootDir("/scratch/buffer/bnl/rootf/");
171 // **************************************************************************
172 
173 // PROTOTYPE 2016 COSMIC & LED DATA
174 // data and ROOT files on RCF
175 // static const TString RCFdataDir("/gpfs/mnt/gpfs02/sphenix/sim/sim01/user04/sphenix/bnl/buffer/");
176 // before April 1st 2016
177 // static const TString RCFdataDir("/gpfs/mnt/gpfs02/sphenix/data/data01/t1044-2016a/cosmics/");
178 // static const TString RCFrootDir("/gpfs/mnt/gpfs02/sphenix/sim/sim01/user04phenix/bnl/buffer/tbsest/");
179 // Beginning April 1st, 2016
180 //static const TString RCFdataDir("/sphenix/data/data01/t1044-2016a/fnal/cosmics/");
181 static const TString RCFdataDir("/sphenix/data/data01/t1044-2016a/fnal/");
182 static const TString RCFrootDir("/sphenix/data/data01/kistenev/rootf/");
183 
184 // data and ROOT files on sphnxdaq
185 static const TString HLABdataDir("/scratch/buffer/bnl/");
186 static const TString HLABrootDir("/scratch/buffer/bnl/rootf/");
187 // **************************************************************************
188 
189 // USEFUL CONSTANTS
190 #define piMass 0.1356
191 
192 
193 
194 #endif // _HCALCONTROLS_H_