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