91 keys = geom0.dt.keys()
93 elif muSystem ==
"CSC":
94 keys = geom0.csc.keys()
98 nkeys, nkeysr, nkeyspass, nmoved, nnotmoved = 0,0,0,0,0
99 nfail_toideal, nfail_deltafinal, nfail_lowstat, nfail_nsigma = 0,0,0,0
100 nok_toideal, nok_deltafinal, nok_lowstat, nok_nsigma = 0,0,0,0
105 if len(key) != 3: is_ch =
False 109 ch_g1 = geom0.dt[ch_key]
110 ch_g2 = geomN.dt[ch_key]
112 if len(key) != 4: is_ch =
False 116 ch_g1 = geom0.csc[ch_key]
117 ch_g2 = geomN.csc[ch_key]
121 if ch_key
in movedChamberKeys:
124 if g1.relativeto != g2.relativeto:
125 print "%s %s relativeto=\"%s\" versus relativeto=\"%s\"" % (muSystem,
str(key), g1.relativeto, g2.relativeto)
130 if r.postal_address[1:] == ch_key:
135 if is_ch:
print muSystem,
str(key),
"not found in the report. Continue..." 139 if rep.status !=
"PASS":
140 if is_ch:
print muSystem,
str(key),
"status is not PASS: %s Continue..." % rep.status
143 if is_ch: nkeyspass+=1
150 if muSystem ==
"DT" and chWasntMoved:
153 if abs(ch_g2.x) > 2.
or abs(ch_g2.y) > 2.
or abs(ch_g2.phiy) > 0.02
or abs(ch_g2.phiz) > 0.02:
157 print "Warning!!!", muSystem,
str(key), \
158 "moved too much with respect to ideal: dx=%.2f mm dy=%.2f mm dphiy=%.2f mrad dphiz=%.2f mrad skipping..." % (ch_g2.x*10, ch_g2.y*10, ch_g2.phiy*1000, ch_g2.phiz*1000)
159 if is_ch
and ok: nok_toideal +=1
164 if abs(rep.deltax.value) > 0.03
or abs(rep.deltay.value) > 0.03
or abs(rep.deltaphiy.value) > 0.0003
or abs(rep.deltaphiz.value) > 0.0006:
167 nfail_deltafinal += 1
168 print "Warning!!!", muSystem,
str(key), \
169 "moved too much at final iteration: dx=%.2f mm dy=%.2f mm dphiy=%.2f mrad dphiz=%.2f mrad skipping..." % \
170 (rep.deltax.value*10, rep.deltay.value*10, rep.deltaphiy.value*1000, rep.deltaphiz.value*1000)
172 if abs(rep.deltax.value) > 0.03
or abs(rep.deltaphiy.value) > 0.0003
or abs(rep.deltaphiz.value) > 0.0006:
175 nfail_deltafinal += 1
176 print "Warning!!!", muSystem,
str(key), \
177 "moved too much at final iteration: dx=%.2f mm dphiy=%.2f mrad dphiz=%.2f mrad skipping..." % \
178 (rep.deltax.value*10, rep.deltaphiy.value*1000, rep.deltaphiz.value*1000)
179 if is_ch
and ok: nok_deltafinal +=1
182 if rep.deltax.error > 0.5:
186 print "Warning!!!", muSystem,
str(key),
"low statistics chamber with too big dx.error = %.2f mm skipping..." % (rep.deltax.error*10.)
187 if is_ch
and ok: nok_lowstat +=1
192 if abs(ch_g1.x - ch_g2.x) < theNSigma * math.sqrt ( rep.deltax.error*rep.deltax.error + 0.02*0.02 ) :
196 print muSystem,
str(key),
"not moved: xN-x0 = %.3f - %.3f = %.3f < %.3f mm" % \
197 ( ch_g2.x*10., ch_g1.x*10., (ch_g2.x-ch_g1.x)*10., theNSigma * math.sqrt ( rep.deltax.error*rep.deltax.error + 0.02*0.02 )*10.)
199 if ok: chWasntMoved =
False 201 if not ok
or chWasntMoved:
continue 206 movedChamberKeys.append(ch_key)
207 print muSystem,
str(key),
"moved!" 211 geom0.dt[key] = copy.copy(geomN.dt[key])
213 geom0.csc[key] = copy.copy(geomN.csc[key])
215 nnotmoved = nkeys - nmoved
216 nsig =
int(theNSigma)
218 %(muSystem)s REPORT for %(nsig)d sigma policy: 220 %(nkeys)d\t chambers in geometry 221 %(nkeysr)d\t chambers in report 222 %(nkeyspass)d\t have PASS status 223 %(nok_toideal)d\t pass big shift with respect to ideal check 224 %(nok_deltafinal)d\t pass big deltas during final iteration 225 %(nok_lowstat)d\t pass low statistics (or big deltax.error) check 228 %(nnotmoved)d\t chambers not moved 229 Numbers of chambers were not moved due to: 230 %(nfail_toideal)d\t big shift with respect to ideal 231 %(nfail_deltafinal)d\t big deltas during final iteration 232 %(nfail_lowstat)d\t low statistics (or big deltax.error) 233 %(nfail_nsigma)d\t |x_final - x_initial| < nsigma
Abs< T >::type abs(const T &t)
1.8.11
