Dear Collaborators,

As promised herewith I am sending my report on the alignment work we
have done so far.

The attached ps file has information on the slides from my 'talk'.

The slides can be found at

www.physics.valpo.edu/Restricted/e614/align091401/

Good luck.

Shirvel

%!PS-Adobe-3.0
%%Title: (Microsoft Word - Alignment_Report_13sept01)
%%Creator: (Microsoft Word: LaserWriter 8 8.6)
%%CreationDate: (8:20 PM Thursday, September 13, 2001)
%%For: (tdss)
%%Routing: (mailto:\000email@domain.com)
%%Pages: 3
%%DocumentFonts: Palatino-Bold Palatino-Roman Palatino-Italic
%%DocumentNeededFonts: Palatino-Bold Palatino-Roman Palatino-Italic
%%DocumentSuppliedFonts:
%%DocumentData: Clean7Bit
%%PageOrder: Ascend
%%Orientation: Portrait
%%DocumentMedia: Default 612 792 0 () ()
%RBINumCopies: 1
%RBINupNess: 1 1
%RBIPCFileName: (HP4050_5.PPD)
%RBIPPDFileVersion: (1.0.1)
%ADO_ImageableArea: 30 33.12 582 761.04
%%EndComments
%%BeginDefaults
%%ViewingOrientation: 1 0 0 1
%%EndDefaults
userdict/dscInfo 5 dict dup begin
/Title(Microsoft Word - Alignment_Report_13sept01)def
/Creator(Microsoft Word: LaserWriter 8 8.6)def
/CreationDate(8:20 PM Thursday, September 13, 2001)def
/For(tdss)def
/Pages 3 def
end put
%%BeginProlog
%%BeginFeature: *JobPatchFile 1

/limitinterpolate { dup type /dicttype eq { dup /Interpolate known { dup /Interpolate get {
	1 1 idtransform 2 index /ImageMatrix get dtransform exch abs 1 exch div exch abs 1 exch div 2 copy lt 
	{ exch }if pop 4 lt { dup /Interpolate false put } if } if } if } if
} bind def
/image { limitinterpolate systemdict /image get exec } bind def
/imagemask { limitinterpolate systemdict /imagemask get exec } bind def
%%EndFeature
/md 183 dict def md begin/currentpacking where {pop /sc_oldpacking currentpacking def true setpacking}if
%%BeginFile: lw8_feature-1.01
%%Copyright: Copyright 1990-1998 Adobe Systems Incorporated and Apple Computer Incorporated. All Rights Reserved.
/bd{bind def}bind def
/ld{load def}bd
/xs{exch store}bd
/Z{0 def}bd
/T true def
/F false def
/level2
/languagelevel where
{
pop languagelevel 2 ge
}{
F
}ifelse
def
/odictstk Z
/oopstk Z
/fcl
{
count oopstk sub dup 0 gt
{
{pop}repeat
}{
pop
}ifelse
countdictstack odictstk sub dup 0 gt
{
{end}repeat
}{
pop
}ifelse
}bd
/sfcl2
{
/odictstk countdictstack store
count/oopstk xs
}bd
/efcl2
{
stopped{$error/newerror F put}if
fcl
}bd
/noload Z
/startnoload
{
{/noload save store}if
}bd
/endnoload
{
{noload restore}if
}bd
/setcopies{
level2
{
1 dict begin/NumCopies exch def currentdict end setpagedevice
}{
userdict/#copies 3 -1 roll put
}ifelse
}def
level2 startnoload
/ststpgdev{}def
/dopgdev{}def
/stpgdev{}def
/buf Z
/didstop T def
/sfcl
{
/didstop T store
/odictstk countdictstack store
count/oopstk xs
currentfile cvx stopped
{
$error/newerror F put
didstop
{
save/didstop xs
/buf vmstatus exch sub exch pop dup 0 lt{pop 0}if
dup 64000 gt{pop 64000}if string store
{
currentfile buf readline
{
(}efcl)eq{exit}if
}{
/UnexpectedEOF errordict/rangecheck get exec
}ifelse
}loop
didstop restore
}if
}if
fcl
}bd
/efcl
{
/didstop F store
exec
stop
}bd
level2 endnoload level2 not startnoload
/setpagedevice where{pop/realstpgdev/setpagedevice ld}if
/SC_topddict Z
/SC_spdict Z
/$spusrdict F def
/dopgdev
{
userdict/setpagedevice undef
$spusrdict
{
userdict/setpagedevice/realstpgdev load put
/$spusrdict F store
}if
SC_topddict realstpgdev
}bd
/stpgdev
{
SC_topddict dup 3 -1 roll
{
SC_spdict 2 index known
{
SC_spdict 2 index get
dup 3 -1 roll
{
put dup
}forall
pop put dup
}{
put dup
}ifelse
}forall
pop pop
}bd
/ststpgdev
{
/setpagedevice where
{
userdict eq
{
/$spusrdict T store
}if
}if
userdict/setpagedevice/stpgdev load put
/SC_topddict 0 dict store
/SC_spdict 3 dict begin
/InputAttributes 0 dict def
/Policies 0 dict def
/OutputAttributes 0 dict def
currentdict
end
store
}def
/sfcl/sfcl2 ld
/efcl/efcl2 ld
level2 not endnoload
%%EndFile
%%BeginFile: lw8_basic-3.0
/xdf{exch def}bd
/:L/lineto
/lw/setlinewidth
/:M/moveto
/rl/rlineto
/rm/rmoveto
/:C/curveto
/:T/translate
/:K/closepath
/:mf/makefont
/gS/gsave
/gR/grestore
/np/newpath
12{ld}repeat
/framewidth -1 def
/QDframwid -1 def
/numframes Z
/mTS matrix def
/$m matrix def
/av 86 def
/por T def
/normland F def
/psb-nosave{}def
/pse-nosave{}def
/us Z
/psb{/us save store}bd
/pse{us restore}bd
/level3
/languagelevel where
{
pop languagelevel 3 ge
}{
F
}ifelse
def
level2 startnoload
/setjob
{
statusdict/jobname 3 -1 roll put
}bd
/devg/DeviceGray def
/devr/DeviceRGB def
/devc/DeviceCMYK def
level2 endnoload level2 not startnoload
/setjob
{
1 dict begin/JobName xdf currentdict end setuserparams
}bd
/devg[/DeviceGray]def
/devr[/DeviceRGB]def
/devc[/DeviceCMYK]def
level2 not endnoload
/pm Z
/mT Z
/sD Z
/initializepage
{
mT $m currentmatrix mTS concatmatrix pop
/pm save store mT concat
}bd
/endp
{
pm restore
}bd
/adjRect
{
dup 2 mul 6 2 roll
4 index sub exch 5 -1 roll sub exch
4 2 roll
4 index add exch 5 -1 roll add exch
4 2 roll
}bd
/frame1up
{
gS
mTS setmatrix
QDframwid lw
/setstrokeadjust where{pop T setstrokeadjust}if
clippath pathbbox
2 index sub exch
3 index sub exch
currentlinewidth framewidth mul
adjRect
numframes dup 0 lt{pop 0}if
{
4 copy
rS
currentlinewidth framewidth
mul 4 mul
adjRect
}repeat
pop pop pop pop
gR
}bd
/$c devr def
/rectclip where
{
pop/rC/rectclip ld
}{
/rC
{
np 4 2 roll
:M
1 index 0 rl
0 exch rl
neg 0 rl
:K
clip np
}bd
}ifelse
/rectfill where
{
pop/rF/rectfill ld
}{
/rF
{
gS
np
4 2 roll
:M
1 index 0 rl
0 exch rl
neg 0 rl
fill
gR
}bd
}ifelse
/rectstroke where
{
pop/rS/rectstroke ld
}{
/rS
{
gS
np
4 2 roll
:M
1 index 0 rl
0 exch rl
neg 0 rl
:K
stroke
gR
}bd
}ifelse
%%EndFile
level3 startnoload
%%BeginFile: lw8_safeclipL12-1.0					
/rectclip where
{
pop/rCa/rectclip ld
}{
/rCa
{
np 0
1 index length 4 idiv
{
2 copy
4 getinterval
aload pop
4 2 roll
:M
1 index 0 rl
0 exch rl
neg 0 rl
:K
4 add
}repeat
clip np
pop pop
}bd
}ifelse
/savedstack Z
/subsavedstack Z
/execstring Z
/saferCa
{
/execstring xs
/odictstk countdictstack store
/oopstk 0 store
count 0 ne
{
savedstack 0 eq{
count 100 gt{count}{100}ifelse
array/savedstack xs
}{
count savedstack length gt
{
count array/savedstack xs
}if
}ifelse
count savedstack 0 3 -1 roll
getinterval
astore/subsavedstack xs
}{
/subsavedstack 0 store
}ifelse
execstring cvx stopped
{
$error/newerror F put
newpath
}if
fcl
subsavedstack 0 ne{
subsavedstack aload pop
}if
}bd
%%EndFile
level3 endnoload level3 not startnoload
%%BeginFile: lw8_safeclipL3-1.0	
/saferCa/pop ld
%%EndFile
level3 not endnoload
%%BeginFile: lw8_level1_colorspace-2.0
/G/setgray ld
/:F1/setgray ld
/:F/setrgbcolor ld
/:F4/setcmykcolor where
{
pop
/setcmykcolor ld
}{
{
3
{
dup
3 -1 roll add
dup 1 gt{pop 1}if
1 exch sub
4 1 roll
}repeat
pop
setrgbcolor
}bd
}ifelse
/:Fx
{
counttomark
{0{G}0{:F}{:F4}}
exch get
exec
pop
}bd
/$cs Z
/:rg{devr :ss}bd
/:sc{$cs :ss}bd
/:dc
{
dup type/arraytype eq{0 get}if
dup/DeviceCMYK eq
{
pop devc
}{
/DeviceGray eq
{
devg
}{
devr
}ifelse
}ifelse
/$cs xdf
}bd
/:sgl{}def
/:dr{}bd
/:fCRD{pop}bd
/:ckcs{}bd
/:ss{/$c xdf}bd
%%EndFile
%%BeginFile: lw8_basic_text-2.1
/S/show ld
/A{
0.0 exch ashow
}bd
/R{
0.0 exch 32 exch widthshow
}bd
/W{
0.0 3 1 roll widthshow
}bd
/J{
0.0 32 4 2 roll 0.0 exch awidthshow
}bd
/V{
0.0 4 1 roll 0.0 exch awidthshow
}bd
/fcflg T def
/fc{
fcflg{
vmstatus exch sub 50000 lt{
(%%[ Warning: Running out of memory ]%%\r)print flush/fcflg F store
}if pop
}if
}bd
/$f[1 0 0 -1 0 0]def
/:ff{$f :mf}bd
/MacEncoding StandardEncoding 256 array copy def
MacEncoding dup 9/space put dup 39/quotesingle put 96/grave put
/Adieresis/Aring/Ccedilla/Eacute/Ntilde/Odieresis/Udieresis/aacute
/agrave/acircumflex/adieresis/atilde/aring/ccedilla/eacute/egrave
/ecircumflex/edieresis/iacute/igrave/icircumflex/idieresis/ntilde/oacute
/ograve/ocircumflex/odieresis/otilde/uacute/ugrave/ucircumflex/udieresis
/dagger/degree/cent/sterling/section/bullet/paragraph/germandbls
/registered/copyright/trademark/acute/dieresis/notequal/AE/Oslash
/infinity/plusminus/lessequal/greaterequal/yen/mu/partialdiff/summation
/product/pi/integral/ordfeminine/ordmasculine/Omega/ae/oslash
/questiondown/exclamdown/logicalnot/radical/florin/approxequal/Delta/guillemotleft
/guillemotright/ellipsis/space/Agrave/Atilde/Otilde/OE/oe
/endash/emdash/quotedblleft/quotedblright/quoteleft/quoteright/divide/lozenge
/ydieresis/Ydieresis/fraction/Euro/guilsinglleft/guilsinglright/fi/fl
/daggerdbl/periodcentered/quotesinglbase/quotedblbase/perthousand
/Acircumflex/Ecircumflex/Aacute/Edieresis/Egrave/Iacute/Icircumflex/Idieresis/Igrave
/Oacute/Ocircumflex/apple/Ograve/Uacute/Ucircumflex/Ugrave/dotlessi/circumflex/tilde
/macron/breve/dotaccent/ring/cedilla/hungarumlaut/ogonek/caron
MacEncoding 128 128 getinterval astore pop
level2 startnoload
/copyfontdict
{
findfont dup length dict
begin
{
1 index/FID ne{def}{pop pop}ifelse
}forall
}bd
/$ckeyd md def
/:skey
{
1 index maxlength 2 index
length sub
ge
{
begin
/$mkeys 20 dict def
/$mkeys load
end
dup/$ckeyd xs
}if
3 1 roll put
}bd
/fD1pass
{
40
$ckeyd
//md
ne
{
pop 1
}if
$ckeyd exch
:skey
}bd
/:searchdict Z
/:searchdict
{
exch 2 copy
known
{
get
}{
exch/$mkeys
get
:searchdict
}ifelse
}bd
/lU
{
//md exch 2 copy
known
{
get
}{
exch
/$mkeys get
:searchdict
}ifelse
}bd
/:xs{
exch 0 1 2 index length 1 sub{
3 copy
gsave 1 getinterval show grestore
exch
get 0.0
rmoveto
}for pop pop
}bd
/:ys{
exch 0 1 2 index length 1 sub{
3 copy
gsave 1 getinterval show grestore
exch
get 0.0 exch
rmoveto
}for pop pop
}bd
/:xys{
exch 0 1 2 index length 1 sub{
3 copy
gsave 1 getinterval show grestore
exch 2 mul 2 copy
1 add get 3 -2 roll get exch
rmoveto
}for pop pop
}bd
level2 endnoload level2 not startnoload
/copyfontdict
{
findfont dup length dict
copy
begin
}bd
/fD1pass/def ld
/lU/load ld
level2{
/:xs/xshow ld
/:ys/yshow ld
/:xys/xyshow ld
}if
level2 not endnoload
/fD Z
/sf Z
/scf Z
/sf1pass
{
lU setfont
}bd
/scf1pass
{
exch lU exch scalefont fD
}bd
/scf2pass
{
scalefont fD
}bd
md/fontname known not{
/fontname/customfont def
}if
/Encoding Z
/:mre
{
copyfontdict
/Encoding MacEncoding def
fontname currentdict
end
definefont :ff fD
}bd
/:bsr
{
copyfontdict
/Encoding Encoding 256 array copy def
Encoding dup
}bd
/pd{put dup}bd
/:esr
{
pop pop
fontname currentdict
end
definefont :ff fD
}bd
/ps Z
/fz{/ps xs}bd
/cF/currentfont ld
/mbf
{
/makeblendedfont where
{
pop
makeblendedfont
/ABlend exch definefont
}{
pop
}ifelse
fD
}def
%%EndFile
/currentpacking where {pop sc_oldpacking setpacking}if end
%%EndProlog
%%BeginSetup
md begin
/fD/def ld/sf/setfont ld /scf/scf2pass ld
 600/languagelevel where{pop languagelevel 2 ge}{false}ifelse{1 dict dup/WaitTimeout 4 -1 roll put setuserparams}{statusdict/waittimeout 3 -1 roll put}ifelse
sfcl{
%%BeginFeature: *Resolution 600x600x2dpi

	<> setpagedevice
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPHalftone PrinterDefault

%%EndFeature

}efcl

sfcl{
%%BeginFeature: *EconoMode PrinterDefault

%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPPaperPolicy PromptUser

   <> setpagedevice
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *Collate True

    <> setpagedevice
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPJobRetention HPJobRetentionOff

	<< /CollateDetails 
	<< /Type 8 /Hold 0 >> >> setpagedevice

%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPUserName FileSharingName

	<< /CollateDetails 
		<< /Type 8
			/UserName /dscInfo where 
				{ /dscInfo get dup /For known
					{/For get}
					{pop (No User Name)} ifelse}
				{(No User Name)}ifelse
			dup length 80 gt { 0 80 getinterval } if
			dup true exch { 32 eq not { false exch pop } if } forall
			{ pop (No User Name) } if
			0 1 2 index length 1 sub
			{ dup 2 index exch get dup 97 ge 1 index 122 le and
				{ 32 sub 2 index 3 1 roll put }
				{ pop pop } ifelse 
			} for
		>>
	>> setpagedevice

%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPJobName DocName

	<< /CollateDetails
		<< /Type 8 
			/JobMgrName /dscInfo where 
				{ /dscInfo get dup /Title known
					{/Title get}
					{pop ()} ifelse}
				{()}ifelse
			/dscInfo where 
				{ /dscInfo get dup /Creator known
					{
						/Creator get dup 0 exch
						{ dup 32 eq exch dup 64 gt exch 122 le and or { 1 add } { exit } ifelse } forall
						0 exch getinterval anchorsearch { pop } if
					}
					{pop} ifelse
				} if
			{ (:) search
				{ pop pop }
				{ exit }
				ifelse
			} loop
			dup 0 exch { false ( -) { 2 index eq or } forall exch pop { 1 add } { exit } ifelse } forall
			dup 0 eq { pop } { dup 2 index length exch sub getinterval } ifelse
			dup length 80 gt { 0 80 getinterval } if
			dup true exch { 32 eq not { false exch pop } if } forall
			{ pop () } if
			0 1 2 index length 1 sub
			{ dup 2 index exch get dup 97 ge 1 index 122 le and
				{ 32 sub 2 index 3 1 roll put }
				{ pop pop } ifelse 
			} for
		>>
	>> setpagedevice

%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPJobHoldKey HoldKeyNone
<< /CollateDetails << /Type 8 /HoldType 0 >> >> setpagedevice
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *Smoothing PrinterDefault

%%EndFeature

}efcl

sfcl{
%%BeginFeature: *ManualFeed False

	<> setpagedevice
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *PageSize LetterSmall

	<> setpagedevice
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *Duplex None

    <> setpagedevice
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPwmBrightness Medium
userdict /HPwmSaturation { [ currentcmykcolor ] /DeviceCMYK setcolorspace { .55 mul } forall setcolor } put
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPwmText Draft
userdict /HPwmText (Draft) put
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPwmFont HelveticaB

  /Helvetica-Bold findfont dup length dict begin
    {1 index /FID ne {def} {pop pop} ifelse} forall
    /MacEncoding where
    { pop /Encoding MacEncoding def }
    { /Encoding ISOLatin1Encoding def } ifelse
    currentdict 
  end
  /HPwmFont exch definefont pop
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPwmFontSize pt48
userdict /HPwmSize 48 put
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPwmTextAngle Deg45
userdict /HPwmAngle 45 put
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPwmTextStyle Medium
userdict /HPwmStyle { HPwmText false charpath HPwmColor .48 setlinewidth stroke } bind put
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPwmPages AllPages
userdict /HPwmLocation true put
%%EndFeature

}efcl

sfcl{
%%BeginFeature: *HPwmSwitch Off

%%EndFeature

}efcl

(tdss)setjob
sfcl{/currentdistillerparams where{pop 1 dict dup/Orientation 2 put setpagedevice}if
}efcl

/mT[.24 0 0 -.24 30.001 761.057]def
/sD 16 dict def
%%IncludeFont: Palatino-Bold
/f44/Palatino-Bold
:mre
/f60 f44 58 scf
/f73 f44 50 scf
/f86 f44 41 scf
%%IncludeFont: Palatino-Roman
/f99/Palatino-Roman
:mre
/f115 f99 41 scf
%%IncludeFont: Palatino-Italic
/f128/Palatino-Italic
:mre
/f145 f128 41 scf
/Courier findfont[10 0 0 -10 0 0]:mf setfont
%%EndSetup
%%Page: 1 1
%%BeginPageSetup
initializepage
(tdss; page: 1 of 3)setjob
%%EndPageSetup
gS 0 0 2300 3033 rC
438 219 :M
f60 sf
.021 .002(Alignment of the TWIST Detector using the August data)J
826 276 :M
f73 sf
.022 .002(\(Notes to accompany the slides\))J
175 416 :M
f86 sf
.012 .001(Slides 1-6 - TheMethod)J
175 420 436 2 rF
175 507 :M
f115 sf
.023 .002(Contents in these slides should be pretty straight forward. The method used for the alignment is outlined.)J
175 553 :M
.023 .002(Some of these were presented at the last collaboration meeting.)J
175 645 :M
f86 sf
.016 .002(Slide 7 - A "complete iteration" Cycle)J
175 649 693 2 rF
175 737 :M
f115 sf
.014 .001(This slide shows a ")J
f145 sf
.027 .003(complete iteration)J
827 737 :M
f115 sf
.017 .002(" cycle in the alignment process.  Using a set of alignment offsets, a set of)J
175 782 :M
.016 .002(data runs are analyzed in )J
645 782 :M
.017 .002(Mofia and a file of hits along the fitted tracks is written out.  This hits file is read by)J
175 828 :M
.017 .002(the alignment code to do the )J
703 828 :M
.022 .002(translational and rotational alignments.  New offsets are determined.  These new)J
175 874 :M
.016 .002(offsets are used in Mofia and the data are analyzed again.)J
175 966 :M
f86 sf
.015 .001(Slides 8-9 - Tests with Monte Carlo)J
175 970 656 2 rF
175 1057 :M
f115 sf
.019 .002(These slides show how well the algorithm works when tested with Monte Carlo. The planes were given)J
175 1103 :M
.016 .002(translational offsets up to 4 mm, and rotational offsets up to 40 )J
1321 1103 :M
.01 .001(mRad.  The MC results )J
1749 1103 :M
.021 .002(show, that when both)J
175 1149 :M
.015 .002(types of offsets are present, we could align the detector to 0.5 micron and 0.02 mRad \(sigma\).  There is one)J
175 1195 :M
.023 .002(entry per plane on each plot.)J
175 1287 :M
f86 sf
.012 .001(Slides 10 - Data Analysis)J
175 1291 464 2 rF
175 1378 :M
f115 sf
.022 .002(This slide shows the current status of the convergence of the translational alignment. After 5 "complete)J
175 1424 :M
.022 .002(iterations", except for plane 22, others have converged. \(A "complete iteration" is defined on slide #7.\))J
175 1470 :M
.014 .001(It takes about one day to analyze the data in )J
984 1470 :M
.019 .002(Mofia and then obtain a new set of offsets.  Within the alignment)J
175 1516 :M
.021 .002(code it goes through about 50 iteration \(but converges after a few\).)J
175 1607 :M
.018 .002(There were a few setbacks at the beginning, when my code did not converge. First because I did not have the)J
175 1653 :M
.023 .002('best' cuts values, then because the Kalman filter did not have the rotations implemented.  Now everything)J
175 1699 :M
.018 .002(appears to work well. Please see the last page for the )J
1138 1699 :M
.019 .002(constants files and cuts used in the analysis.)J
175 1791 :M
f86 sf
(Slide 11)S
175 1795 149 2 rF
175 1882 :M
f115 sf
.018 .002(This slide shows the current values of the )J
936 1882 :M
.029 .003(translational offsets. Very little further change is expected in the)J
175 1928 :M
.019 .002(next few iterations. Only the offset of plane 22 is expected to change.  \(see a note on plane 22 later\). The)J
175 1974 :M
.014 .001(maximum offset is close to 600 microns.)J
175 2066 :M
f86 sf
(Slides 12-13)S
175 2070 222 2 rF
175 2158 :M
f115 sf
.019 .002(These slides show the convergence of rotational offsets. Again plane 22 has not converged. The maximum)J
175 2203 :M
.014 .001(offset is about 1 )J
470 2203 :M
.018 .002(mRad. The offsets shown in slide 13 may suggest that there is a twist in the detector \(i.e. the)J
175 2249 :M
.025 .003(fixed planes are not aligned\).)J
175 2341 :M
f86 sf
.018 .002(Slide 14 - Effects of Translational Alignment)J
175 2345 833 2 rF
175 2433 :M
f115 sf
.019 .002(Here you will see the residuals for plane 4 with and without the )J
1343 2433 :M
.023 .002(translational offsets. You can clearly see the)J
175 2478 :M
.02 .002(shifting of the peak. In this and subsequent plots, a gaussian function is fit to the peak and an exponential)J
175 2524 :M
.016 .002(term is added to fit the tails. When aligned, this function appears to do a good fit.)J
175 2616 :M
.02 .002(Note also the number of entries in each histogram. You can see very clearly how the alignment affects the)J
175 2662 :M
.02 .002(number of successfully fit tracks in the Kalman filter. More on this in a later slide.)J
175 2753 :M
f86 sf
.018 .002(Slide 15 - Evicence for Rotational Offsets)J
175 2757 765 2 rF
175 2845 :M
f115 sf
.021 .002(This slide shows the effects of rotational alignment. Because you expect to see rotational effects only at large u)J
175 2891 :M
.013 .001(or v  \(i.e. at the ends of the wires\) these two histograms show v-residuals of plane 12 when u > 5.0 cm and)J
endp
showpage
%%Page: 2 2
%%BeginPageSetup
initializepage
(tdss; page: 2 of 3)setjob
%%EndPageSetup
gS 0 0 2300 3033 rC
175 205 :M
f115 sf
.015 .002(when u < 5.0 cm.  If there are rotational offsets \(in additional to )J
1328 205 :M
.031 .003(translational offsets\) the distribution of the)J
175 251 :M
.013 .001(residuals is expected to shift in opposite directions when u > 0 or u < 0 \(or v >0 and v < 0\). This is )J
f145 sf
(clearly seen)S
175 297 :M
f115 sf
.02 .002(in these two plots \(one for the "upper half" and the other for the "lower half" of the plane\). )J
1811 297 :M
f145 sf
.016 .002(With no rotational)J
175 342 :M
.018 .002(alignment, the means are shifted by equal amounts in opposite directions)J
f115 sf
.015 .001(.  Note also the asymmetry of the curves -)J
175 388 :M
.02 .002(in opposite directions.  However, if you add both halves together, the effect will cancel out and the mean will)J
175 434 :M
.016 .002(be centered as shown in Slide #14.)J
175 526 :M
f86 sf
.021 .002(Slides 16 - Results for Rotational Alignment)J
175 530 824 2 rF
175 617 :M
f115 sf
.017 .002(These are the same two plots as on Slide #15, shown after the rotational alignment.  The means in the upper)J
175 663 :M
.02 .002(half and lower half are both centered at zero and the distributions are nice.)J
175 755 :M
f86 sf
(Slides 17-18)S
175 759 222 2 rF
175 847 :M
f115 sf
.018 .002(Slide 17 shows you how the magnitudes of the means of the residuals improved after each 'complete)J
175 892 :M
.02 .002(iteration'. Slide 18 shows the current means after the last iteration. It is clear that plane 22 has not converged.)J
175 938 :M
.017 .002(The means for all other planes are below 10 microns.)J
175 1030 :M
f86 sf
(Slide 19)S
175 1034 149 2 rF
175 1122 :M
f115 sf
.018 .002(This slide shows the sigmas of the residual distributions. You will see how the )J
1602 1122 :M
.017 .002(sigmas improved )J
1928 1122 :M
.028 .003(with each)J
175 1167 :M
.008(iteration.)A
175 1259 :M
f86 sf
(Slide 20)S
175 1263 149 2 rF
175 1351 :M
f115 sf
.025 .003(Here I have defined a term 'tracking efficiency' where)J
325 1442 :M
.022 .002(Tracking efficiency = number of 'good' event in a plane/total number of events in the run.)J
175 1534 :M
.025 .002(Here the 'good' refers to events passing the Kalman filter.)J
175 1626 :M
.02 .002(You can see how this efficiency is changing with each 'complete iteration' of alignment.  The lowest 3 are)J
175 1672 :M
.021 .002(planes 20, 21, and 22. When the alignment has fully converged, I expect the 'tracking efficiencies' of all the)J
175 1717 :M
.019 .002(planes to plateau out. It is still not flat.)J
175 1809 :M
f86 sf
(Slide 21)S
175 1813 149 2 rF
175 1901 :M
f115 sf
.017 .002(I have also tried to extract some drift chamber efficiencies as a )J
1307 1901 :M
.02 .002(by product of the alignment. For this I used)J
175 1947 :M
.018 .002(events that are close to the axis. These are expected to go through all the chambers. In my event selection I)J
175 1992 :M
.017 .002(have required planes 1, 2, 15 and 16. Then I calculated the efficiency by finding out how many )J
1888 1992 :M
.014 .001(times  each of)J
175 2038 :M
.02 .002(the other chambers have a hit. Hence these efficiencies are relative to the least efficient of planes \(1, 2, 15, or)J
175 2084 :M
.021 .002(16\). My results indicate that our chambers are below 90% efficient.)J
endp
showpage
%%Page: 3 3
%%BeginPageSetup
initializepage
(tdss; page: 3 of 3)setjob
%%EndPageSetup
gS 0 0 2300 3033 rC
175 212 :M
f73 sf
.012 .001(Files and cuts used in the present analysis)J
175 217 957 3 rF
175 306 :M
f115 sf
.016 .002(At present I am using the following cuts and files and cuts )J
1242 306 :M
.022 .002(\(these were suggested by Maher about 3 weeks)J
175 352 :M
.019 .002(ago. If these have changed please let me know\):)J
175 443 :M
(files)S
325 489 :M
(dc_t0.1790m)S
1225 489 :M
(pc_t0.1789)S
325 535 :M
(sc_t0.1789)S
1225 535 :M
(dc_res.00004)S
175 627 :M
(name sccuts)S
325 718 :M
(S1_MAX_NHITS   =    -1)S
1225 718 :M
(S2_MAX_NHITS   =    -1)S
325 764 :M
(S1_WIDTH_CUT =   -1.00)S
1225 764 :M
(S2_WIDTH_CUT =   -1.00)S
325 810 :M
(S1_MAX_TDC   =   -1.00)S
1225 810 :M
(S2_MAX_TDC   =   -1.00)S
325 856 :M
(S1_MIN_TDC   =   -1.00)S
1225 856 :M
(S2_MIN_TDC   =  )S
1557 856 :M
( -1.00)S
325 902 :M
(S1_MAX_ADC   =   -1.00)S
1225 902 :M
(S2_MAX_ADC   =   -1.00)S
325 947 :M
(S1_MIN_ADC   =   -1.00)S
1225 947 :M
(S2_MIN_ADC   =   -1.00)S
175 1039 :M
.029 .003(name rfcuts)J
325 1085 :M
(RF_MIN_TDC = -1.)S
1225 1085 :M
(RF_MAX_TDC = -1.)S
175 1177 :M
(name dccuts)S
325 1222 :M
(DC_MAXTDC_CUT = -1)S
1225 1222 :M
(DC_MINTDC_CUT = -1)S
325 1268 :M
(DC_MAXWTDC_CUT    = -1.)S
1225 1268 :M
(DC_MINWTDC_CUT    = -1.)S
325 1314 :M
(DC_MAX_HITS_IN_PLANE = 1)S
1225 1314 :M
(DC_MIN_PLANES      = 12)S
325 1360 :M
(DC_XTALK_WCUT    = 100.)S
1225 1360 :M
(DC_NOISE_WCUT    = 100.)S
175 1452 :M
(name pccuts)S
325 1497 :M
(PC_MAXTDC_CUT = 200.)S
1225 1497 :M
(PC_MINTDC_CUT = 0.)S
325 1543 :M
(PC_MAXWTDC_CUT    = -1.)S
1225 1543 :M
(PC_MINWTDC_CUT    = -1.)S
325 1589 :M
(PC_MAX_HITS_IN_PLANE = 3)S
1225 1589 :M
(PC_MIN_PLANES      = 3)S
325 1635 :M
(PC_XTALK_WCUT    = 5.)S
1225 1635 :M
(PC_NOISE_WCUT    = 5.)S
175 1727 :M
(name adccuts)S
325 1772 :M
(PC_Emin = 0.)S
1225 1772 :M
(PC_Emax = 1.e7)S
325 1818 :M
(SC_Emin = 0.)S
1225 1818 :M
(SC_Emax = 1.e7)S
175 1956 :M
.03 .003(name KalmanCuts)J
325 2002 :M
(ChiDiffCluster = 1.E-1)S
1225 2002 :M
.02 .002(ChiDiffTime = 1.E-1)J
325 2048 :M
(Chi2sCutCluster = 100000.)S
1225 2048 :M
.017 .002(Chi2sCutTime = 5.)J
325 2093 :M
(MaxIterateCluster = 5)S
1225 2093 :M
.021 .002(MaxIterateTime = 10)J
325 2139 :M
(enableKalman = .TRUE.)S
1225 2139 :M
.021 .002(NoiseExcludeMax = 2)J
175 2231 :M
.03 .003(name strset)J
325 2277 :M
(use)S
387 2277 :M
(_cos_increments = .FALSE.)S
1225 2277 :M
(str)S
1272 2277 :M
.018 .002(_angle_inc = 1.0)J
325 2323 :M
(str)S
372 2323 :M
.029 .003(_upper_angle_limit = 80.0)J
175 2414 :M
.028 .003(name TimeZero)J
325 2460 :M
(T0_TDC_MIN = 49995.)S
1225 2460 :M
(T0_TDC_MAX= 50005)S
325 2506 :M
(FitT0 )S
429 2506 :M
(= .FALSE.)S
1225 2506 :M
(TimeBackwards=)S
1543 2506 :M
(.TRUE)S
1663 2506 :M
(.)S
175 2598 :M
(name Rezs)S
325 2643 :M
(nEventsMax=2000)S
1225 2643 :M
.007(MinHistFitEntries=200)A
325 2689 :M
.004(calcResidualMethod=3)A
175 2781 :M
.037 .004(name Efficiency)J
325 2827 :M
(RadiusCut = 15.0)S
1225 2827 :M
.017 .002(CellCut = 2)J
325 2873 :M
(nHitPlanesCut = 3)S
1225 2873 :M
.017 .002(FindEff = .FALSE.)J
endp
showpage
%%Trailer
end
%%EOF

Created for the The Center for Subatomic Research E614 Project Projects Page.

Created by The CoCoBoard.