93 keys = geom0.dt.keys()
95 elif muSystem ==
"CSC":
96 keys = geom0.csc.keys()
100 nkeys, nkeysr, nkeyspass, nmoved, nnotmoved = 0,0,0,0,0
101 nfail_toideal, nfail_deltafinal, nfail_lowstat, nfail_nsigma = 0,0,0,0
102 nok_toideal, nok_deltafinal, nok_lowstat, nok_nsigma = 0,0,0,0
107 if len(key) != 3: is_ch =
False 111 ch_g1 = geom0.dt[ch_key]
112 ch_g2 = geomN.dt[ch_key]
114 if len(key) != 4: is_ch =
False 118 ch_g1 = geom0.csc[ch_key]
119 ch_g2 = geomN.csc[ch_key]
123 if ch_key
in movedChamberKeys:
126 if g1.relativeto != g2.relativeto:
127 print(
"%s %s relativeto=\"%s\" versus relativeto=\"%s\"" % (muSystem,
str(key), g1.relativeto, g2.relativeto))
132 if r.postal_address[1:] == ch_key:
137 if is_ch:
print(muSystem,
str(key),
"not found in the report. Continue...")
141 if rep.status !=
"PASS":
142 if is_ch:
print(muSystem,
str(key),
"status is not PASS: %s Continue..." % rep.status)
145 if is_ch: nkeyspass+=1
152 if muSystem ==
"DT" and chWasntMoved:
155 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:
159 print(
"Warning!!!", muSystem,
str(key), \
160 "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))
161 if is_ch
and ok: nok_toideal +=1
166 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:
169 nfail_deltafinal += 1
170 print(
"Warning!!!", muSystem,
str(key), \
171 "moved too much at final iteration: dx=%.2f mm dy=%.2f mm dphiy=%.2f mrad dphiz=%.2f mrad skipping..." % \
172 (rep.deltax.value*10, rep.deltay.value*10, rep.deltaphiy.value*1000, rep.deltaphiz.value*1000))
174 if abs(rep.deltax.value) > 0.03
or abs(rep.deltaphiy.value) > 0.0003
or abs(rep.deltaphiz.value) > 0.0006:
177 nfail_deltafinal += 1
178 print(
"Warning!!!", muSystem,
str(key), \
179 "moved too much at final iteration: dx=%.2f mm dphiy=%.2f mrad dphiz=%.2f mrad skipping..." % \
180 (rep.deltax.value*10, rep.deltaphiy.value*1000, rep.deltaphiz.value*1000))
181 if is_ch
and ok: nok_deltafinal +=1
184 if rep.deltax.error > 0.5:
188 print(
"Warning!!!", muSystem,
str(key),
"low statistics chamber with too big dx.error = %.2f mm skipping..." % (rep.deltax.error*10.))
189 if is_ch
and ok: nok_lowstat +=1
194 if abs(ch_g1.x - ch_g2.x) < theNSigma * math.sqrt ( rep.deltax.error*rep.deltax.error + 0.02*0.02 ) :
198 print(muSystem,
str(key),
"not moved: xN-x0 = %.3f - %.3f = %.3f < %.3f mm" % \
199 ( 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.))
201 if ok: chWasntMoved =
False 203 if not ok
or chWasntMoved:
continue 208 movedChamberKeys.append(ch_key)
213 geom0.dt[key] = copy.copy(geomN.dt[key])
215 geom0.csc[key] = copy.copy(geomN.csc[key])
217 nnotmoved = nkeys - nmoved
218 nsig =
int(theNSigma)
220 %(muSystem)s REPORT for %(nsig)d sigma policy: 222 %(nkeys)d\t chambers in geometry 223 %(nkeysr)d\t chambers in report 224 %(nkeyspass)d\t have PASS status 225 %(nok_toideal)d\t pass big shift with respect to ideal check 226 %(nok_deltafinal)d\t pass big deltas during final iteration 227 %(nok_lowstat)d\t pass low statistics (or big deltax.error) check 230 %(nnotmoved)d\t chambers not moved 231 Numbers of chambers were not moved due to: 232 %(nfail_toideal)d\t big shift with respect to ideal 233 %(nfail_deltafinal)d\t big deltas during final iteration 234 %(nfail_lowstat)d\t low statistics (or big deltax.error) 235 %(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
