92 keys = geom0.dt.keys()
94 elif muSystem ==
"CSC":
95 keys = geom0.csc.keys()
99 nkeys, nkeysr, nkeyspass, nmoved, nnotmoved = 0,0,0,0,0
100 nfail_toideal, nfail_deltafinal, nfail_lowstat, nfail_nsigma = 0,0,0,0
101 nok_toideal, nok_deltafinal, nok_lowstat, nok_nsigma = 0,0,0,0
106 if len(key) != 3: is_ch =
False 110 ch_g1 = geom0.dt[ch_key]
111 ch_g2 = geomN.dt[ch_key]
113 if len(key) != 4: is_ch =
False 117 ch_g1 = geom0.csc[ch_key]
118 ch_g2 = geomN.csc[ch_key]
122 if ch_key
in movedChamberKeys:
125 if g1.relativeto != g2.relativeto:
126 print(
"%s %s relativeto=\"%s\" versus relativeto=\"%s\"" % (muSystem,
str(key), g1.relativeto, g2.relativeto))
131 if r.postal_address[1:] == ch_key:
136 if is_ch:
print(muSystem,
str(key),
"not found in the report. Continue...")
140 if rep.status !=
"PASS":
141 if is_ch:
print(muSystem,
str(key),
"status is not PASS: %s Continue..." % rep.status)
144 if is_ch: nkeyspass+=1
151 if muSystem ==
"DT" and chWasntMoved:
154 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:
158 print(
"Warning!!!", muSystem,
str(key), \
159 "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))
160 if is_ch
and ok: nok_toideal +=1
165 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:
168 nfail_deltafinal += 1
169 print(
"Warning!!!", muSystem,
str(key), \
170 "moved too much at final iteration: dx=%.2f mm dy=%.2f mm dphiy=%.2f mrad dphiz=%.2f mrad skipping..." % \
171 (rep.deltax.value*10, rep.deltay.value*10, rep.deltaphiy.value*1000, rep.deltaphiz.value*1000))
173 if abs(rep.deltax.value) > 0.03
or abs(rep.deltaphiy.value) > 0.0003
or abs(rep.deltaphiz.value) > 0.0006:
176 nfail_deltafinal += 1
177 print(
"Warning!!!", muSystem,
str(key), \
178 "moved too much at final iteration: dx=%.2f mm dphiy=%.2f mrad dphiz=%.2f mrad skipping..." % \
179 (rep.deltax.value*10, rep.deltaphiy.value*1000, rep.deltaphiz.value*1000))
180 if is_ch
and ok: nok_deltafinal +=1
183 if rep.deltax.error > 0.5:
187 print(
"Warning!!!", muSystem,
str(key),
"low statistics chamber with too big dx.error = %.2f mm skipping..." % (rep.deltax.error*10.))
188 if is_ch
and ok: nok_lowstat +=1
193 if abs(ch_g1.x - ch_g2.x) < theNSigma * math.sqrt ( rep.deltax.error*rep.deltax.error + 0.02*0.02 ) :
197 print(muSystem,
str(key),
"not moved: xN-x0 = %.3f - %.3f = %.3f < %.3f mm" % \
198 ( 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.))
200 if ok: chWasntMoved =
False 202 if not ok
or chWasntMoved:
continue 207 movedChamberKeys.append(ch_key)
212 geom0.dt[key] = copy.copy(geomN.dt[key])
214 geom0.csc[key] = copy.copy(geomN.csc[key])
216 nnotmoved = nkeys - nmoved
217 nsig =
int(theNSigma)
219 %(muSystem)s REPORT for %(nsig)d sigma policy: 221 %(nkeys)d\t chambers in geometry 222 %(nkeysr)d\t chambers in report 223 %(nkeyspass)d\t have PASS status 224 %(nok_toideal)d\t pass big shift with respect to ideal check 225 %(nok_deltafinal)d\t pass big deltas during final iteration 226 %(nok_lowstat)d\t pass low statistics (or big deltax.error) check 229 %(nnotmoved)d\t chambers not moved 230 Numbers of chambers were not moved due to: 231 %(nfail_toideal)d\t big shift with respect to ideal 232 %(nfail_deltafinal)d\t big deltas during final iteration 233 %(nfail_lowstat)d\t low statistics (or big deltax.error) 234 %(nfail_nsigma)d\t |x_final - x_initial| < nsigma S & print(S &os, JobReport::InputFile const &f)
Abs< T >::type abs(const T &t)
1.8.11
