4 bool fixCDCWires =
true;
6 bool fixCDCGlobalX =
false;
7 bool fixCDCGlobalY =
false;
8 bool fixCDCGlobalZ =
true;
9 bool fixCDCGlobalPhiX =
true;
10 bool fixCDCGlobalPhiY =
true;
11 bool fixCDCGlobalPhiZ =
true;
13 bool fixFDCGains =
true;
14 bool fixFDCCathodeOffsets =
true;
15 bool fixFDCCathodeAngles =
true;
16 bool fixFDCCellOffsetsWires =
true;
17 bool fixFDCCellOffsetsCathodes =
true;
18 bool fixFDCWireRotationX =
true;
19 bool fixFDCWireRotationY =
true;
20 bool fixFDCWireRotationZ =
true;
22 bool fixFDCPitch =
true;
23 bool fixFDCGap =
true;
26 double translationPresigma = 0.0005;
27 double rotationPresigma = 0.0001;
28 double gainPresigma = 0.01;
29 double pitchPresigma = 0.001;
30 double t0Presigma=10.0;
33 outfile.open(
"Parameters.txt");
35 outfile <<
"Parameter" << endl;
38 if(fixCDCGlobalX) outfile <<
"1 0.0 -1" << endl;
39 else outfile <<
"1 0.0 " << translationPresigma << endl;
41 if(fixCDCGlobalY) outfile <<
"2 0.0 -1" << endl;
42 else outfile <<
"2 0.0 " << translationPresigma << endl;
44 if(fixCDCGlobalZ) outfile <<
"3 0.0 -1" << endl;
45 else outfile <<
"3 0.0 " << translationPresigma << endl;
47 if(fixCDCGlobalPhiX) outfile <<
"4 0.0 -1" << endl;
48 else outfile <<
"4 0.0 " << rotationPresigma << endl;
50 if(fixCDCGlobalPhiY) outfile <<
"5 0.0 -1" << endl;
51 else outfile <<
"5 0.0 " << rotationPresigma << endl;
53 if(fixCDCGlobalPhiZ) outfile <<
"6 0.0 -1" << endl;
54 else outfile <<
"6 0.0 " << rotationPresigma << endl;
56 for (
unsigned int i=16001; i <= 19522; i++){
57 if (fixCDCt0) outfile << i <<
" 0.0 -1.0" << endl;
58 else outfile << i <<
" 0.0 " << t0Presigma << endl;
61 for (
unsigned int i=1001; i <= 15088; i++){
62 if (fixCDCWires) outfile << i <<
" 0.0 -1.0" << endl;
63 else outfile << i <<
" 0.0 " << translationPresigma << endl;
95 for (
unsigned int i=1; i<=24; i++){
97 int indexOffset = 100000 + i*1000;
99 outfile<< indexOffset + 5 <<
" 0.0 -1.0" << endl;
102 outfile<< indexOffset + 5 <<
" 0.0 " << translationPresigma << endl;
104 if (fixFDCCellOffsetsWires){
105 outfile << indexOffset + 1 <<
" 0.0 -1.0" << endl;
108 outfile << indexOffset + 1 <<
" 0.0 " << translationPresigma << endl;
110 if (fixFDCCellOffsetsCathodes){
111 outfile << indexOffset + 100 <<
" 0.0 -1.0" << endl;
114 outfile << indexOffset + 100 <<
" 0.0 " << translationPresigma << endl;
117 if (fixFDCWireRotationX) outfile << indexOffset + 2 <<
" 0.0 -1.0" << endl;
118 else outfile << indexOffset + 2 <<
" 0.0 " << rotationPresigma << endl;
120 if (fixFDCWireRotationY) outfile << indexOffset + 3 <<
" 0.0 -1.0" << endl;
121 else outfile << indexOffset + 3 <<
" 0.0 " << rotationPresigma << endl;
123 if (fixFDCWireRotationZ) outfile << indexOffset + 4 <<
" 0.0 -1.0" << endl;
124 else outfile << indexOffset + 4 <<
" 0.0 " << rotationPresigma << endl;
126 if (fixFDCCathodeOffsets){
127 outfile << indexOffset + 101 <<
" 0.0 -1.0" << endl;
128 outfile << indexOffset + 102 <<
" 0.0 -1.0" << endl;
131 outfile << indexOffset + 101 <<
" 0.0 " << translationPresigma << endl;
132 outfile << indexOffset + 102 <<
" 0.0 " << translationPresigma << endl;
135 if (fixFDCCathodeAngles){
137 outfile << indexOffset + 103 <<
" 0.0 -1.0" << endl;
138 outfile << indexOffset + 104 <<
" 0.0 -1.0" << endl;
142 outfile << indexOffset + 103 <<
" 0.0 " << rotationPresigma << endl;
143 outfile << indexOffset + 104 <<
" 0.0 " << rotationPresigma << endl;
148 outfile << indexOffset + 200 <<
" 0.0 -1.0" << endl;
149 outfile << indexOffset + 202 <<
" 0.0 -1.0" << endl;
150 outfile << indexOffset + 204 <<
" 0.0 -1.0" << endl;
151 outfile << indexOffset + 205 <<
" 0.0 -1.0" << endl;
152 outfile << indexOffset + 207 <<
" 0.0 -1.0" << endl;
153 outfile << indexOffset + 209 <<
" 0.0 -1.0" << endl;
156 outfile << indexOffset + 200 <<
" 0.0 " << pitchPresigma << endl;
157 outfile << indexOffset + 202 <<
" 0.0 " << pitchPresigma << endl;
158 outfile << indexOffset + 204 <<
" 0.0 " << pitchPresigma << endl;
159 outfile << indexOffset + 205 <<
" 0.0 " << pitchPresigma << endl;
160 outfile << indexOffset + 207 <<
" 0.0 " << pitchPresigma << endl;
161 outfile << indexOffset + 209 <<
" 0.0 " << pitchPresigma << endl;
166 outfile << indexOffset + 201 <<
" 0.0 -1.0" << endl;
167 outfile << indexOffset + 203 <<
" 0.0 -1.0" << endl;
168 outfile << indexOffset + 206 <<
" 0.0 -1.0" << endl;
169 outfile << indexOffset + 208 <<
" 0.0 -1.0" << endl;
172 outfile << indexOffset + 201 <<
" 0.0 " << translationPresigma << endl;
173 outfile << indexOffset + 203 <<
" 0.0 " << translationPresigma << endl;
174 outfile << indexOffset + 206 <<
" 0.0 " << translationPresigma << endl;
175 outfile << indexOffset + 208 <<
" 0.0 " << translationPresigma << endl;
179 for (
unsigned int j = 301; j <= 517; j++){
180 if (fixFDCGains) outfile << indexOffset + j <<
" 0.0 -1.0" << endl;
181 else outfile << indexOffset + j <<
" 0.0 " << gainPresigma << endl;
184 for (
unsigned int j = 601; j <= 817; j++){
185 if (fixFDCGains) outfile << indexOffset + j <<
" 0.0 -1.0" << endl;
186 else outfile << indexOffset + j <<
" 0.0 " << gainPresigma << endl;
190 for (
unsigned int j = 901; j <= 996; j++){
191 if (fixFDCt0) outfile << indexOffset + j <<
" 0.0 -1.0" << endl;
192 else outfile << indexOffset + j <<
" 0.0 " << t0Presigma << endl;