CMS 3D CMS Logo
CSCDriftParamLowB.cc
Go to the documentation of this file.
1 #include "SimMuon/CSCDigitizer/src/CSCDriftSim.h"
2 #include <cmath>
3 #include <iostream>
4 
5 // the following functions were all generated by f2c.
6 
7 double CSCDriftSim::avgPathLengthLowB() {
8 
9  static const double coeff[7] = { .28872209,.21705601,.063908389,-.012224924,
10  .012901814,-.0089355058,-.015769145 };
11 
12  // System generated locals
13  float ret_val;
14 
15  // Local variables
16  double x10, x11, x12, x13, x20, x21, x22;
17 
18 /* ! Parameterization of drift path length - low field chambers */
19 /* ***********************************************************************
20  */
21 /* DOC MC_BHGH_SLEN *
22  */
23 /* *
24 */
25 /* DOC Function : Parameterization of the drift path length *
26  */
27 /* DOC in the muon endcap CSCs. *
28  */
29 /* *
30 */
31 /* DOC References : None *
32  */
33 /* *
34 */
35 /* DOC Arguments : YCELL - distance from the anode wire in the *
36  */
37 /* anode-cathode coordinate plane *
38 */
39 /* DOC ZCELL - distance from the anode wire in the wire *
40  */
41 /* plane coordinate *
42 */
43 /* DOC Errors : None *
44  */
45 /* *
46 */
47 /* DOC Returns : Drift path length for low field CSC chambers *
48 */
49 /* *
50 */
51 /* DOC Created : 15-OCT-1996 Author : Jeff Rowe *
52  */
53 /* ***********************************************************************
54  */
55 
56 
57  x10 = 1.;
58  x11 = fabs(ycell) * 2. - 1.;
59  x12 = x11 * 2. * x11 - x10;
60  x13 = x11 * 2. * x12 - x11;
61  x20 = 1.;
62  x21 = fabs(zcell) * 2. - 1.;
63  x22 = x21 * 2. * x21 - x20;
64  ret_val = coeff[0] + coeff[1] * x11 + coeff[2] * x21 + coeff[3] * x12 *
65  x22 + coeff[4] * x22 + coeff[5] * x13 + coeff[6] * x11 * x21;
66  return ret_val;
67 }
68 
69 
70 double CSCDriftSim::pathSigmaLowB() {
71  /* Initialized data */
72 
73  static const double coeff[12] = { .002854441,8.701339e-4,.0053064193,
74  .0012356508,9.8627318e-4,.0013802449,-5.4633755e-4,.0026078648,
75  .0044171026,6.2029063e-4,.0011392474,-.0056275595 };
76 
77  /* System generated locals */
78  float ret_val;
79 
80  /* Local variables */
81  double /*x10,*/ x11, x12, x13, x14, x15, x16, /*x20,*/ x21, x22, x23,
82  x24, x25, x26, x27, x28, x29, x210, x211, x212, x213;
83 
84 /* ! Parameterization of path length dispersion- low field chambers */
85 /* ***********************************************************************
86  */
87 /* DOC MC_BHGH_SSIG *
88  */
89 /* *
90 */
91 /* DOC Function : Parameterization of the drift path length *
92  */
93 /* DOC dispersion in the muon endcap CSCs. *
94  */
95 /* *
96 */
97 /* DOC References : None *
98  */
99 /* *
100 */
101 /* DOC Arguments : YCELL - distance from the anode wire in the *
102  */
103 /* anode-cathode coordinate plane *
104 */
105 /* DOC ZCELL - distance from the anode wire in the wire *
106  */
107 /* plane coordinate *
108 */
109 /* **NOTE** Both distances normalize to cell dim=1x1 *
110 */
111 /* DOC Errors : None *
112  */
113 /* *
114 */
115 /* DOC Returns : Path length dispersion for low field CSC chambers *
116 */
117 /* *
118 */
119 /* DOC Created : 15-OCT-1996 Author : Jeff Rowe *
120  */
121 /* ***********************************************************************
122  */
123 
124 
125  //x10 = 1.; //not used later
126  x11 = fabs(ycell) * 2. - 1.;
127  x12 = x11 * x11;
128  x13 = x11 * x12;
129  x14 = x11 * x13;
130  x15 = x11 * x14;
131  x16 = x11 * x15;
132  //x20 = 1.; //not used later
133  x21 = fabs(zcell) * 2. - 1.;
134  x22 = x21 * x21;
135  x23 = x21 * x22;
136  x24 = x21 * x23;
137  x25 = x21 * x24;
138  x26 = x21 * x25;
139  x27 = x21 * x26;
140  x28 = x21 * x27;
141  x29 = x21 * x28;
142  x210 = x21 * x29;
143  x211 = x21 * x210;
144  x212 = x21 * x211;
145  x213 = x21 * x212;
146  ret_val = coeff[0] + coeff[1] * x23 + coeff[2] * x24 + coeff[3] * x13 +
147  coeff[4] * x11 * x21 + coeff[5] * x11 * x22 + coeff[6] * x16 +
148  coeff[7] * x213 + coeff[8] * x212 + coeff[9] * x21 + coeff[10] *
149  x11 * x23 + coeff[11] * x26;
150 
151  return ret_val;
152 }
153 
154 
155 double CSCDriftSim::avgDriftTimeLowB() {
156  /* Initialized data */
157 
158  static const double coeff[20] = { 42.981588,25.732805,26.539129,16.719016,
159  10.862044,7.4859085,5.0353142,3.3620548,1.9057762,2.2207695,
160  -2.6787582,1.2977292,-.8358091,1.2452612,.74813469,-.57581194,
161  .32705275,-.85426489,-.55688158,-.38384903 };
162 
163  /* System generated locals */
164  float ret_val;
165 
166  /* Local variables */
167  double x10, x11, x12, x13, x14, x20, x21, x22, x23, x24, x25,
168  x26, x27, x28, x29;
169 
170 /* ! Parameterization of drift time - low field chambers */
171 /* ***********************************************************************
172  */
173 /* DOC MC_BLOW_TIME *
174  */
175 /* *
176 */
177 /* DOC Function : Parameterization of the drift time *
178  */
179 /* DOC in the muon endcap CSCs. *
180  */
181 /* *
182 */
183 /* DOC References : None *
184  */
185 /* *
186 */
187 /* DOC Arguments : YCELL - distance from the anode wire in the *
188  */
189 /* anode-cathode coordinate plane *
190 */
191 /* (ycell=1 > d_acat) *
192 */
193 /* DOC ZCELL - distance from the anode wire in the wire *
194  */
195 /* plane coordinate (zcell=1 > d_anod/2.) *
196 */
197 /* DOC Errors : None *
198  */
199 /* *
200 */
201 /* DOC Returns : Drift time for low field CSC chambers *
202  */
203 /* *
204 */
205 /* DOC Created : 15-OCT-1996 Author : Jeff Rowe *
206  */
207 /* ***********************************************************************
208  */
209 
210 
211 
212 
213  x10 = 1.;
214  x11 = fabs(ycell) * 2. - 1.;
215  x12 = x11 * 2. * x11 - x10;
216  x13 = x11 * 2. * x12 - x11;
217  x14 = x11 * 2. * x13 - x12;
218  x20 = 1.;
219  x21 = fabs(zcell) * 2. - 1.;
220  x22 = x21 * 2. * x21 - x20;
221  x23 = x21 * 2. * x22 - x21;
222  x24 = x21 * 2. * x23 - x22;
223  x25 = x21 * 2. * x24 - x23;
224  x26 = x21 * 2. * x25 - x24;
225  x27 = x21 * 2. * x26 - x25;
226  x28 = x21 * 2. * x27 - x26;
227  x29 = x21 * 2. * x28 - x27;
228 
229  ret_val = coeff[0] + coeff[1] * x11 + coeff[2] * x21 + coeff[3] * x22 +
230  coeff[4] * x23 + coeff[5] * x24 + coeff[6] * x25 + coeff[7] * x26
231  + coeff[8] * x12 + coeff[9] * x27 + coeff[10] * x11 * x21 + coeff[
232  11] * x28 + coeff[12] * x13 + coeff[13] * x12 * x21 + coeff[14] *
233  x29 + coeff[15] * x13 * x21 + coeff[16] * x14 + coeff[17] * x11 *
234  x22 + coeff[18] * x11 * x23 + coeff[19] * x11 * x24;
235  return ret_val;
236 }
237 
238 
239 double CSCDriftSim::driftTimeSigmaLowB() {
240  /* Initialized data */
241 
242  static const double coeff[19] = { 6.2681223,3.5916437,5.5425809,4.6974052,
243  3.8472392,3.1019155,2.4323913,1.8695623,1.3909068,.99056625,
244  .52066638,.64671229,.68023389,.40512251,.40202738,.23908531,
245  -.41245784,-.32196924,-.29890696 };
246 
247  /* Local variables */
248  double /*x10,*/ x11, x20, x21, x22, x23, x24, x25, x26, x27, x28,
249  x29, x210, x211, x212;
250 
251 /* ! Parameterization of drift time dispersion- low field chambers */
252 /* ***********************************************************************
253  */
254 /* DOC MC_BHGH_TSIG *
255  */
256 /* *
257 */
258 /* DOC Function : Parameterization of the drift time dispersion *
259  */
260 /* DOC in the muon endcap CSCs. *
261  */
262 /* *
263 */
264 /* DOC References : None *
265  */
266 /* *
267 */
268 /* DOC Arguments : YCELL - distance from the anode wire in the *
269  */
270 /* anode-cathode coordinate plane *
271 */
272 /* (ycell=1 > d_acat) *
273 */
274 /* DOC ZCELL - distance from the anode wire in the wire *
275  */
276 /* plane coordinate (zcell=1 > d_anod/2.) *
277 */
278 /* DOC Errors : None *
279  */
280 /* *
281 */
282 /* DOC Returns : Drift time dispersion for low field CSC chambers *
283  */
284 /* *
285 */
286 /* DOC Created : 15-OCT-1996 Author : Jeff Rowe *
287  */
288 /* ***********************************************************************
289  */
290 
291 
292 
293  //x10 = 1.; //not used
294  x11 = fabs(ycell) * 2. - 1.;
295  x20 = 1.;
296  x21 = fabs(zcell) * 2. - 1.;
297  x22 = x21 * 2. * x21 - x20;
298  x23 = x21 * 2. * x22 - x21;
299  x24 = x21 * 2. * x23 - x22;
300  x25 = x21 * 2. * x24 - x23;
301  x26 = x21 * 2. * x25 - x24;
302  x27 = x21 * 2. * x26 - x25;
303  x28 = x21 * 2. * x27 - x26;
304  x29 = x21 * 2. * x28 - x27;
305  x210 = x21 * 2. * x29 - x28;
306  x211 = x21 * 2. * x210 - x29;
307  x212 = x21 * 2. * x211 - x210;
308 
309  return coeff[0] * x21 + coeff[1] + coeff[2] * x22 + coeff[3] * x23 +
310  coeff[4] * x24 + coeff[5] * x25 + coeff[6] * x26 + coeff[7] * x27
311  + coeff[8] * x28 + coeff[9] * x29 + coeff[10] * x11 + coeff[11] *
312  x11 * x21 + coeff[12] * x210 + coeff[13] * x11 * x22 + coeff[14] *
313  x211 + coeff[15] * x212 + coeff[16] * x11 * x25 + coeff[17] *
314  x11 * x27 + coeff[18] * x11 * x26;
315 }
316 
CSCDriftSim::zcell
double zcell
Definition: CSCDriftSim.h:66
CSCDriftSim::driftTimeSigmaLowB
double driftTimeSigmaLowB()
Definition: CSCDriftParamLowB.cc:239
CSCDriftSim::avgDriftTimeLowB
double avgDriftTimeLowB()
Definition: CSCDriftParamLowB.cc:155
CSCDriftSim::pathSigmaLowB
double pathSigmaLowB()
Definition: CSCDriftParamLowB.cc:70
CSCDriftSim::ycell
double ycell
Definition: CSCDriftSim.h:66
CSCDriftSim::avgPathLengthLowB
double avgPathLengthLowB()
Definition: CSCDriftParamLowB.cc:7
Generated for CMSSW Reference Manual by doxygen 1.8.11