Log File of my COCOA work (resumed in 2007): -------------------------------------------- 1/15/2007 --------- > source muon_lpc_setup.cs > cd /uscms/home/hohlmann/work/cocoa > cd /uscmst1/prod/sw/cmssw/Releases/CMSSW/CMSSW_0_7_2 > eval `scramv1 run -csh` > cd /uscms/home/hohlmann/work/cocoa > scramv1 project CMSSW CMSSW_0_7_2 works ok and sets up the area: ... Installation procedure complete. Developer area located at: /uscms/home/hohlmann/work/cocoa/CMSSW_0_7_2 > cd /uscms/home/hohlmann/work/cocoa/CMSSW_0_7_2/src > cmscvsroot CMSSW > cvs login (pwd 98passwd) > cvs co Alignment checks out a whole bunch of code > cvs co CondFormats ditto... > cd /uscms/home/hohlmann/work/cocoa/CMSSW_0_7_2/src/Alignment > ls CocoaAnalysis/ CommonAlignment/ LTDRootSource/ CocoaApplication/ CommonAlignmentAlgorithm/ MillePedeAlignmentAlgorithm/ CocoaDaq/ CommonAlignmentParametrization/ MuonAlignment/ CocoaDDLObjects/ CommonAlignmentProducer/ MuonStandaloneAlgorithm/ CocoaFit/ CSA06AlignmentAlgorithm/ OptAlignASCIISource/ CocoaModel/ CVS/ OptAlignCSV2DB/ CocoaToDDL/ doc/ OptAlignGeneratedSource/ CocoaUtilities/ KalmanAlignmentAlgorithm/ TrackerAlignment/ > cd CocoaApplication/bin > ls BuildFile CVS/ simple2D.txt simple2DWithMirror.txt cocoa.cpp simple2D.meas simple2DWithMirror.meas > eval `scramv1 run -csh` > scramv1 build compiles a whole bunch of files > which cocoa /uscmst1/prod/sw/cmssw/Releases/CMSSW/CMSSW_0_7_2/bin/slc3_ia32_gcc323/cocoa > cd /uscms/home/hohlmann/work/cocoa/CMSSW_0_7_2/bin/slc3_ia32_gcc323 > ls cocoa* > eval `scramv1 run -csh` > which cocoa /uscms/home/hohlmann/work/cocoa/CMSSW_0_7_2/bin/slc3_ia32_gcc323/cocoa > cd /uscms/home/hohlmann/work/cocoa/CMSSW_0_7_2 I HAVE A NAGGING SUSPICION THAT THIS IS AN OLD VERSION OF CMSSW !! In GB's user guide it lists "scram list" instead of "scramv1 list", so I only got the old versions of CMSSW! Do it over: > scramv1 list CMSSW (as opposed to scram list CMSSW) reveals much more up-to-date version: ... CMSSW CMSSW_1_2_0_pre5 --> /uscmst1/prod/sw/cms/slc3_ia32_gcc323/cms/cmssw/CMSSW_1_2_0_pre5 CMSSW CMSSW_1_2_0_pre6 --> /uscmst1/prod/sw/cms/slc3_ia32_gcc323/cms/cmssw/CMSSW_1_2_0_pre6 CMSSW CMSSW_1_1_2 --> /uscmst1/prod/sw/cms/slc3_ia32_gcc323/cms/cmssw/CMSSW_1_1_2 CMSSW CMSSW_1_2_0_pre8 --> /uscmst1/prod/sw/cms/slc3_ia32_gcc323/cms/cmssw/CMSSW_1_2_0_pre8 CMSSW CMSSW_1_2_0_pre9 --> /uscmst1/prod/sw/cms/slc3_ia32_gcc323/cms/cmssw/CMSSW_1_2_0_pre9 CMSSW CMSSW_1_2_0 --> /uscmst1/prod/sw/cms/slc3_ia32_gcc323/cms/cmssw/CMSSW_1_2_0 CMSSW CMSSW_1_2_0_g4_81 --> /uscmst1/prod/sw/cms/slc3_ia32_gcc323/cms/cmssw/CMSSW_1_2_0_g4_81 CMSSW CMSSW_0_5_0 ... So, let's start this over with the latest version which appears to be CMSSW_1_2_0 ! > cd /uscmst1/prod/sw/cms/slc3_ia32_gcc323/cms/cmssw/CMSSW_1_2_0 > eval `scramv1 run -csh` > cd /uscms/home/hohlmann/work/cocoa > scramv1 project CMSSW CMSSW_1_2_0 ... Installation procedure complete. Developer area located at: /uscms/home/hohlmann/work/cocoa/CMSSW_1_2_0 > cd /uscms/home/hohlmann/work/cocoa/CMSSW_1_2_0/src > project CMSSW > cmscvsroot CMSSW > cvs login (pwd 98passwd) > cvs co Alignment checks out a whole bunch of code > cvs co CondFormats ditto... > cd /uscms/home/hohlmann/work/cocoa/CMSSW_1_2_0/src/Alignment/CocoaApplication/bin > eval `scramv1 run -csh` > scramv1 build > cd /uscms/home/hohlmann/work/cocoa/CMSSW_1_2_0/bin/slc3_ia32_gcc323 > ls cocoa* > eval `scramv1 run -csh` > which cocoa /uscms/home/hohlmann/work/cocoa/CMSSW_1_2_0/bin/slc3_ia32_gcc323/cocoa > cd /uscms/home/hohlmann/work/cocoa/CMSSW_1_2_0 > mkdir Analysis > cd Analysis > mkdir OldAnalysis > cd /uscms/home/hohlmann/work/cocoa/COCOA_3_2_5/src/Analysis/MyAnalysis > cp * /uscms/home/hohlmann/work/cocoa/CMSSW_1_2_0/Analysis/OldAnalysis I copy Gyongyi's text based input files for COCOA to my U: drive at my.fit.edu, so I can download them from the cmsuaf: > cd /uscms/home/hohlmann/work/cocoa/CMSSW_1_2_0/Analysis > sftp hohlmann@my.fit.edu > cd Cocoa (sic!) > cd Gyongyis_work_2006 > ls ME2_SLM25_GB_2006.txt slm25_b0t_data.txt slm25_b4t_data.txt slm_simu.txt > mget * > exit Ok, so let's give this a try and see if we can actually run COCOA: > cocoa ME2_SLM25_GB_2006.txt ... FIT STATUS COCOA_FitOK 0x9f7de08$$$$$$$$$$$$$$$ moreData Sets 0 ............ program ended OK TIME:PROGRAM ENDED : 2390000 0.22 TIME:TOTAL PROGRAM : 2390000 1.74 Yeah, it actually works ! It produces 3 output files: report.out fittedEntries.out longFittedEntries.out The last two files appear to contain the fit values for the chamber centers and angles. Point the "work" alias to the new directory: > xemacs .cshrc & edit alias work ... accordingly > cd > source .cshrc to make new alias effective 1/17/07 ------- > work > ev > which cocoa /uscms/home/hohlmann/work/cocoa/CMSSW_1_2_0/bin/slc3_ia32_gcc323/cocoa I verify in the report.out file that I can reproduce Gyongyi's values for positions of chambers and apins with the nominal positions. Takes one iteration, FIT quality is 2E-12. For 0T real data: FIT quality is 135.7 after 9 iterations; however, at the end of Iteration 8 it shows a FIT QUALITY of 2.04 ! Interesting... The numbers for chamber positions seem to be close to Gyopngyi's, but not identical. For 4T: 5 iterations, FIT Q=0.04 (<0.1 so it stops) I copied the report.out files to appropriately named .out files. Ok, try to introduce some SHIMS ! This is how the DCOPS position is defined: //#14: short inner (1) chamber, Point5 Schamber me+21_14 ... COPS dcops_outer+21_14>in ENTRY { length rightCCDYtoDowel2 bi_cops_hCCDcal cal_error_YCCD cal length leftCCDYtoDowel2 bi_cops_hCCDcal cal_error_YCCD cal length rightCCDXtoDowel2 bi_rightCCDX_disp cal_error_XCCD cal length leftCCDXtoDowel2 bi_leftCCDX2_disp cal_error_XCCD cal length upCCDXtoDowel2 bi_cops_vCCDcal cal_error_XCCD cal length upCCDYtoDowel2 bi_upCCDY_disp cal_error_YCCD cal length downCCDXtoDowel2 bi_cops_vCCDcal cal_error_XCCD cal length downCCDYtoDowel2 bi_downCCDY_disp cal_error_YCCD cal } centre X dowel_out2_yref dowel_out2_yref_error cal Y dowel_out2_zref dowel_out2_zref_error cal Z dowel_out2_xref dowel_out2_xref_error cal This is in the CSC ref. frame. Note that for X, it uses ..Yref, etc. Here are the actual values from the definition section of the sdf: dowel_out2_xref 992.094 dowel_out2_yref 272.081 dowel_out2_zref 50.100 dowel_out2_zref_tower 0 So, centre X 272.081 dowel_out2_yref_error cal Y 50.100 dowel_out2_zref_error cal Z 992.094 dowel_out2_xref_error cal Note this is in mm. Based on Gyongyi's drawing in Fig.4 of her note, I would conclude that the shimming direction, i.e. perpendicular to the chamber surface is the Y_CSC direction. It's called dowel_out2_zref, which makes me think that the "z" in "zref" actually marks the CMS direction, because Y_CSC = Z_CMS (barring a minus sign). Let's look at the inner DCOPS to confirm this: COPS dcops_inner+21_14>in ENTRY { length rightCCDYtoDowel2 bi_cops_hCCDcal cal_error_YCCD cal length leftCCDYtoDowel2 bi_cops_hCCDcal cal_error_YCCD cal length rightCCDXtoDowel2 bi_rightCCDX_disp cal_error_XCCD cal length leftCCDXtoDowel2 bi_leftCCDX2_disp cal_error_XCCD cal length upCCDXtoDowel2 bi_cops_vCCDcal cal_error_XCCD cal length upCCDYtoDowel2 bi_upCCDY_disp cal_error_YCCD cal length downCCDXtoDowel2 bi_cops_vCCDcal cal_error_XCCD cal length downCCDYtoDowel2 bi_downCCDY_disp cal_error_YCCD cal } centre X -dowel_in2_yref dowel_in2_yref_error cal Y dowel_in2_zref dowel_in2_zref_error cal Z -dowel_in2_xref dowel_in2_xref_error cal angles X 0 prec_dowel_ang cal Y 10 prec_dowel_ang cal Z 0 prec_dowel_ang cal where the values are dowel_in2_xref 884.330 dowel_in2_yref 13.107 dowel_in2_zref 50.100 dowel_in2_zref_tower 0 So, for the inner DCOPS on this short 20-deg chamber: centre X - 13.107 cal Y 50.100 cal Z -884.330 cal From Gyongyi's Fig.4, both coordinates in the chamber plane are negative, so this confirms that Y=50.1mm is the perpendicular coordinate. Good! For the other short chamber ME2_1_04, I see the same thing. ok. What about the long 10-deg chambers ? For ME+2_2_27 we have for the 2 DCOPS: centre X dowel_out_yref dowel_out_yref_error cal Y dowel_out_zref dowel_out_zref_error cal Z dowel_out_xref dowel_out_xref_error cal and centre X dowel_in_yref dowel_in_yref_error cal Y dowel_in_zref dowel_in_zref_error cal Z -dowel_in_xref dowel_in_xref_error cal with dowel_out_zref 50.100 dowel_in_zref 50.100 Same values as for short chamber ! For ME+2_2_08 I see the same thing, so that is all consistent. So the perpendicular DCOPS distances are in sequence (radially out -> in) LONG outer DCOPS dowel_out_zref inner DCOPS dowel_in_zref SHORT outer DCOPS dowel_out2_zref inner DCOPS dowel_in2_zref and the same thing on the other side. CONCLUSION: These parameter would have to be changed according to whatever shim was applied. Here is the shimming information from Dave Eartly via a file provided by Oleg: Long Long Short Short Ref ME2-2 Outer ME2-2 Inner ME2-1 Outer ME2-1 Inner 0.5mm (up) 0mm -4mm (down) -5mm (down) -8mm (down) So, based on that our new values including shimming should be: Name Change Absolute (all values in mm) value dowel_out_zref unchanged 50.1 dowel_in_zref -4.0 46.1 dowel_out2_zref -5.0 45.1 dowel_in2_zref -8.0 42.1 I rerun B=0T data with these values adjusted for the shims. Program finished ok after 9 iteration. FIT Q=2342.5 "has worsened". I did not include the shim in the ref. sensor, yet. After iteration 8 it is FIT Q=219.9. Fitted Chamber positions after Iteration 8 are X_SLM Y_SLM Z_SLM ME+2_2_08 (long) -5225.8 -79.6 169.4 ME+2_1_04 (short) -2371.3 -84.4 133.5 ME+2_1_14 (short) 2371.3 -83.9 120.8 ME+2_2_27 (long) 5225.8 -92.1 163.1 Basically all chambers are shifted to higher (=less negative) Y_SLM numbers compared to Gyongyi's Fig.8. Long chambers shift by 2-3mm; short chambers by 5-6mm. Note, all these fits give us very large rotations (order of tens of degrees) for the chambers. That is unphysical and I think we need to restrict them, e.g. by setting the rotations to "calibrated". Tomorrow... 1/19/2007 --------- Still need to put in the shims for the 2 reference DCOPS. According to Dave, it's 0.5mm up (= off the transferplate). Note that the SLM reference DCOPS are defined on the transfer plate and in fact the DCOPS dowel pin is defined as the center of the TP system: COPS me+2_reference5>in ... centre X 0 prec_pos_tp cal Y 0 prec_pos_tp cal Z 0 prec_pos_tp cal angles X 0 prec_ang_tp cal Y 0 prec_ang_tp cal Z 0 prec_ang_tp cal Based on Fig.6 in Gyongyi's note, the coordinate off the TP, i.e. parallel to Y_SLM is Y_TP, so I change this to centre X 0 prec_pos_tp cal Y 0.5 prec_pos_tp cal Z 0 prec_pos_tp cal for both SLM ref. DCOPS Result for B=0T 5 iteration, final FQ=1257 has worsened from fit 4: FQ=417 Run on B=4T with shims: Get good fit quality !!! FQ=0.012 after 6 iterations. 1/31/07: -------- The file structure on cmsuaf has changed. Home directory is now /uscms/home/hohlmann. Old afs home directory is available under the link afshome (/afs/fnal.gov/files/home/room1/hohlmann). I copy all files from there to the new home at /uscms/home/hohlmann and will continue from there, including with this file itself. I get too much output; try > calcul_type 1 (only unk values) Crahes the program ! Set back to 0. Try > report_verbose 2 This doesn't seem to shorten the output. Oh well... Running with shims on B=4T data I get for chamber centers in SLM coordinates: %%%% Optical Object: OCMS/slm22/me+21_14 CAL: 80 OCMS/slm22/me+21_14 centre_X 2371.2849 +- 99.997889 2371.279 +- 100 Q1 UNK: 81 OCMS/slm22/me+21_14 centre_Y -91.66121 +- 419.45851 -93.35 +- 100 Q2 UNK: 82 OCMS/slm22/me+21_14 centre_Z 113.91896 +- 185.91443 118.16258 +- 100 Q2 UNK: 83 OCMS/slm22/me+21_14 angles_X 12.46115 +- 7427.2526 0 +- 100 Q2 UNK: 84 OCMS/slm22/me+21_14 angles_Y 79.752854 +- 303.61173 80 +- 100 Q2 UNK: 85 OCMS/slm22/me+21_14 angles_Z 12.165216 +- 7315.397 0 +- 100 Q2 %%%% Optical Object: OCMS/slm22/me+21_04 CAL: 114 OCMS/slm22/me+21_04 centre_X -2371.2748 +- 99.99878 -2371.279 +- 100 Q1 UNK: 115 OCMS/slm22/me+21_04 centre_Y -91.00122 +- 436.68727 -93.35 +- 100 Q2 UNK: 116 OCMS/slm22/me+21_04 centre_Z 118.04817 +- 219.06667 118.16258 +- 100 Q2 UNK: 117 OCMS/slm22/me+21_04 angles_X 5.8358397 +- 9586.7357 0 +- 100 Q2 UNK: 118 OCMS/slm22/me+21_04 angles_Y -79.906834 +- 189.77356 -80 +- 100 Q2 UNK: 119 OCMS/slm22/me+21_04 angles_Z -5.6627869 +- 9444.9601 0 +- 100 Q2 %%%% Optical Object: OCMS/slm22/me+22_27 CAL: 148 OCMS/slm22/me+22_27 centre_X 5225.8427 +- 99.99909 5225.848 +- 100 Q1 UNK: 149 OCMS/slm22/me+22_27 centre_Y -93.037968 +- 440.74824 -93.35 +- 100 Q2 UNK: 150 OCMS/slm22/me+22_27 centre_Z 159.79689 +- 181.15744 157.20274 +- 100 Q2 UNK: 151 OCMS/slm22/me+22_27 angles_X 10.713041 +- 13611.598 0 +- 100 Q2 UNK: 152 OCMS/slm22/me+22_27 angles_Y 84.90371 +- 236.00481 85 +- 100 Q2 UNK: 153 OCMS/slm22/me+22_27 angles_Z 10.639026 +- 13560.143 0 +- 100 Q2 %%%% Optical Object: OCMS/slm22/me+22_08 CAL: 182 OCMS/slm22/me+22_08 centre_X -5225.8492 +- 99.998992 -5225.848 +- 100 Q1 UNK: 183 OCMS/slm22/me+22_08 centre_Y -91.420375 +- 468.23978 -93.35 +- 100 Q2 UNK: 184 OCMS/slm22/me+22_08 centre_Z 154.35718 +- 221.54191 157.20274 +- 100 Q2 UNK: 185 OCMS/slm22/me+22_08 angles_X 10.422072 +- 17409.428 0 +- 100 Q2 UNK: 186 OCMS/slm22/me+22_08 angles_Y -84.802861 +- 286.04385 -85 +- 100 Q2 UNK: 187 OCMS/slm22/me+22_08 angles_Z -10.357462 +- 17344.302 0 +- 100 Q2 Note how the fit prefers very similar rotations around SLM_x and SLM_z for each chamber: 12deg, 5deg, 10deg, 10deg. I force the rotations around x and z to more or less 0 by setting them as CAL coordinates in the sdf: Schamber me+21_14 centre X 2371.279 prec_chamber_pos_det cal Y -93.350 prec_chamber_pos_det unk Z 118.16258 prec_chamber_pos_det unk //old 118.164 angles X 0 prec_chamber_ang_det cal Y 80 prec_chamber_ang_det unk Z 0 prec_chamber_ang_det cal COCOA still fits, but now with FQ=5.2. Now results are %%%% Optical Object: OCMS/slm22/me+21_14 CAL: 80 OCMS/slm22/me+21_14 centre_X 2371.2851 +- 99.997886 2371.279 +- 100 Q1 UNK: 81 OCMS/slm22/me+21_14 centre_Y -87.337883 +- 389.47274 -93.35 +- 100 Q2 UNK: 82 OCMS/slm22/me+21_14 centre_Z 117.52311 +- 132.46871 118.16258 +- 100 Q2 CAL: 83 OCMS/slm22/me+21_14 angles_X 0.0072327292 +- 97.291332 0 +- 100 Q1 UNK: 84 OCMS/slm22/me+21_14 angles_Y 79.981443 +- 114.06522 80 +- 100 Q2 CAL: 85 OCMS/slm22/me+21_14 angles_Z -0.0052327517 +- 97.210778 0 +- 100 Q1 whereas before it was %%%% Optical Object: OCMS/slm22/me+21_14 CAL: 80 OCMS/slm22/me+21_14 centre_X 2371.2849 +- 99.997889 2371.279 +- 100 Q1 UNK: 81 OCMS/slm22/me+21_14 centre_Y -91.66121 +- 419.45851 -93.35 +- 100 Q2 UNK: 82 OCMS/slm22/me+21_14 centre_Z 113.91896 +- 185.91443 118.16258 +- 100 Q2 UNK: 83 OCMS/slm22/me+21_14 angles_X 12.46115 +- 7427.2526 0 +- 100 Q2 UNK: 84 OCMS/slm22/me+21_14 angles_Y 79.752854 +- 303.61173 80 +- 100 Q2 UNK: 85 OCMS/slm22/me+21_14 angles_Z 12.165216 +- 7315.397 0 +- 100 Q2 Note that the center_Y has shifted by ~4.3mm and the center_Z has also shifted by ~3.5mm. The CMS coordinates of the center for that chamber have changed as follows: with Angles UNK: OCMS/slm22/me+21_14 CENTRE GLOBAL: (-0.21391873,-2.3976154,8.0896612) with Angles CAL: OCMS/slm22/me+21_14 CENTRE GLOBAL: (-0.21043743,-2.3985484,8.0853379) i.e. a shift by 3.5mm in CMS_x and by 0.9mm in CMS_y. Next, I restrict the chamber angles to CAL for all 4 chambers: Results: Now the fit converges after only 2 iterations (!) with FQ=0.029 UNK: 6 OCMS/slm22/me+2_tp5/las25 angleBetweenAxis 90.311325 +- 905.57065 90 +- 100 Q2 %%%% Optical Object: OCMS/slm22/me+21_14 CAL: 80 OCMS/slm22/me+21_14 centre_X 2371.2849 +- 99.997886 2371.279 +- 100 Q1 UNK: 81 OCMS/slm22/me+21_14 centre_Y -87.344768 +- 389.47193 -93.35 +- 100 Q2 UNK: 82 OCMS/slm22/me+21_14 centre_Z 117.52208 +- 132.46837 118.16258 +- 100 Q2 CAL: 83 OCMS/slm22/me+21_14 angles_X 0.0071313375 +- 97.291164 0 +- 100 Q1 UNK: 84 OCMS/slm22/me+21_14 angles_Y 79.981514 +- 114.06524 80 +- 100 Q2 CAL: 85 OCMS/slm22/me+21_14 angles_Z -0.0053361151 +- 97.210635 0 +- 100 Q1 %%%% Optical Object: OCMS/slm22/me+21_04 CAL: 114 OCMS/slm22/me+21_04 centre_X -2371.2749 +- 99.998778 -2371.279 +- 100 Q1 UNK: 115 OCMS/slm22/me+21_04 centre_Y -88.995279 +- 390.07122 -93.35 +- 100 Q2 UNK: 116 OCMS/slm22/me+21_04 centre_Z 119.87129 +- 137.83507 118.16258 +- 100 Q2 CAL: 117 OCMS/slm22/me+21_04 angles_X 0.004767152 +- 97.300542 0 +- 100 Q1 UNK: 118 OCMS/slm22/me+21_04 angles_Y -79.943583 +- 119.81002 -80 +- 100 Q2 CAL: 119 OCMS/slm22/me+21_04 angles_Z 0.004570281 +- 97.218527 0 +- 100 Q1 %%%% Optical Object: OCMS/slm22/me+22_27 CAL: 148 OCMS/slm22/me+22_27 centre_X 5225.8426 +- 99.999089 5225.848 +- 100 Q1 UNK: 149 OCMS/slm22/me+22_27 centre_Y -90.484374 +- 392.86516 -93.35 +- 100 Q2 UNK: 150 OCMS/slm22/me+22_27 centre_Z 161.33317 +- 139.86528 157.20274 +- 100 Q2 CAL: 151 OCMS/slm22/me+22_27 angles_X 0.004070609 +- 94.009761 0 +- 100 Q1 UNK: 152 OCMS/slm22/me+22_27 angles_Y 84.990932 +- 72.33144 85 +- 100 Q2 CAL: 153 OCMS/slm22/me+22_27 angles_Z -0.0037977256 +- 93.964569 0 +- 100 Q1 %%%% Optical Object: OCMS/slm22/me+22_08 CAL: 182 OCMS/slm22/me+22_08 centre_X -5225.8487 +- 99.998989 -5225.848 +- 100 Q1 UNK: 183 OCMS/slm22/me+22_08 centre_Y -88.910243 +- 395.88835 -93.35 +- 100 Q2 UNK: 184 OCMS/slm22/me+22_08 centre_Z 156.08711 +- 146.08221 157.20274 +- 100 Q2 CAL: 185 OCMS/slm22/me+22_08 angles_X 0.00016079018 +- 99.998378 0 +- 100 Q1 UNK: 186 OCMS/slm22/me+22_08 angles_Y -84.877652 +- 79.045325 -85 +- 100 Q2 UNK: 187 OCMS/slm22/me+22_08 angles_Z 0.021861869 +- 275.46317 0 +- 100 Q2 and in global coordinates: OCMS/slm22/me+21_14 CENTRE GLOBAL: (-0.21043838,-2.398548,8.0853448) OCMS/slm22/me+21_04 CENTRE GLOBAL: (1.0192956,2.181805,8.0869953) OCMS/slm22/me+22_27 CENTRE GLOBAL: (-0.90693402,-5.1671781,8.0884844) OCMS/slm22/me+22_08 CENTRE GLOBAL: (1.7930954,4.9297382,8.0869102) Let's apply this to B=0T data and see what we get: After 11 fits FQ=180 and cannot improve. At least that is an improvement over FQ=417 that we got before (see above). Results are: UNK: 6 OCMS/slm22/me+2_tp5/las25 angleBetweenAxis 91.291283 +- 1026.0625 90 +- 100 Q2 %%%% Optical Object: OCMS/slm22/me+21_14 CAL: 80 OCMS/slm22/me+21_14 centre_X 2371.2797 +- 99.999951 2371.279 +- 100 Q1 UNK: 81 OCMS/slm22/me+21_14 centre_Y -83.789752 +- 390.63542 -93.35 +- 100 Q2 UNK: 82 OCMS/slm22/me+21_14 centre_Z 120.93852 +- 141.71115 118.16258 +- 100 Q2 CAL: 83 OCMS/slm22/me+21_14 angles_X 0.0066877698 +- 97.31348 0 +- 100 Q1 UNK: 84 OCMS/slm22/me+21_14 angles_Y 79.909078 +- 123.24599 80 +- 100 Q2 CAL: 85 OCMS/slm22/me+21_14 angles_Z -0.0066970732 +- 97.223173 0 +- 100 Q1 %%% Optical Object: OCMS/slm22/me+21_04 CAL: 114 OCMS/slm22/me+21_04 centre_X -2371.2806 +- 99.999974 -2371.279 +- 100 Q1 UNK: 115 OCMS/slm22/me+21_04 centre_Y -91.1083 +- 392.44969 -93.35 +- 100 Q2 UNK: 116 OCMS/slm22/me+21_04 centre_Z 128.19445 +- 307.84167 118.16258 +- 100 Q2 CAL: 117 OCMS/slm22/me+21_04 angles_X 0.0152544 +- 97.335925 0 +- 100 Q1 UNK: 118 OCMS/slm22/me+21_04 angles_Y -79.684086 +- 334.7513 -80 +- 100 Q2 CAL: 119 OCMS/slm22/me+21_04 angles_Z 0.016829964 +- 97.24299 0 +- 100 Q1 %%%% Optical Object: OCMS/slm22/me+22_27 CAL: 148 OCMS/slm22/me+22_27 centre_X 5225.8463 +- 99.99996 5225.848 +- 100 Q1 UNK: 149 OCMS/slm22/me+22_27 centre_Y -89.153392 +- 392.9219 -93.35 +- 100 Q2 UNK: 150 OCMS/slm22/me+22_27 centre_Z 164.80692 +- 148.38542 157.20274 +- 100 Q2 CAL: 151 OCMS/slm22/me+22_27 angles_X 0.0099871235 +- 94.011842 0 +- 100 Q1 UNK: 152 OCMS/slm22/me+22_27 angles_Y 85.002765 +- 79.26029 85 +- 100 Q2 CAL: 153 OCMS/slm22/me+22_27 angles_Z -0.0095248623 +- 93.965283 0 +- 100 Q1 %%%% Optical Object: OCMS/slm22/me+22_08 CAL: 182 OCMS/slm22/me+22_08 centre_X -5225.848 +- 99.999998 -5225.848 +- 100 Q1 UNK: 183 OCMS/slm22/me+22_08 centre_Y -89.670911 +- 396.63008 -93.35 +- 100 Q2 UNK: 184 OCMS/slm22/me+22_08 centre_Z 164.00271 +- 298.83628 157.20274 +- 100 Q2 CAL: 185 OCMS/slm22/me+22_08 angles_X -0.0010002623 +- 99.998603 0 +- 100 Q1 UNK: 186 OCMS/slm22/me+22_08 angles_Y -84.569845 +- 213.62565 -85 +- 100 Q2 UNK: 187 OCMS/slm22/me+22_08 angles_Z -0.042410335 +- 275.54668 0 +- 100 Q2 In global CMS coordinates: OCMS/slm22/me+21_14 CENTRE GLOBAL: (-0.207137,-2.3994271,8.0817898) OCMS/slm22/me+21_14/apin_inner+21_14 CENTRE GLOBAL: (-0.12422504,-1.4340983,8.0815652) OCMS/slm22/me+21_14/apin_outer+21_14 CENTRE GLOBAL: (-0.29170088,-3.3839892,8.0820188) OCMS/slm22/me+21_04 CENTRE GLOBAL: (1.0273366,2.1796563,8.0891083) OCMS/slm22/me+21_04/apin_inner+21_04 CENTRE GLOBAL: (0.61303361,1.3038213,8.0885703) OCMS/slm22/me+21_04/apin_outer+21_04 CENTRE GLOBAL: (1.4498942,3.0729414,8.089657) OCMS/slm22/me+22_27 CENTRE GLOBAL: (-0.90357959,-5.1680808,8.0871534) OCMS/slm22/me+22_27/apin_inner+22_27 CENTRE GLOBAL: (-0.61608006,-3.5380498,8.0865908) OCMS/slm22/me+22_27/apin_outer+22_27 CENTRE GLOBAL: (-1.1910791,-6.7981117,8.087716) OCMS/slm22/me+22_08 CENTRE GLOBAL: (1.8007411,4.9276888,8.0876709) OCMS/slm22/me+22_08/apin_inner+22_08 CENTRE GLOBAL: (1.2229715,3.3766124,8.0889195) OCMS/slm22/me+22_08/apin_outer+22_08 CENTRE GLOBAL: (2.3785108,6.4787652,8.0864223) - Interestingly, if I run this on simulated data, it takes 12 iterations to get FQ=0.17 now. Hmmm.. Ah, but there shouldn't be any shims included when running on simulated data, so that's why this is now messed up and why the Y_SLM coordinates of the chamber centers do not come out dead-center at -93.35mm as we would expect for simulated data. Just to be sure I double check that by using GB's original file WITHOUT shimming and changing the chamber angles from UNK to CAL and run on simulated data. Lo and behold it converges after one fit with FQ=2.2E-12 and all chamber angles at 0deg and Y_SLMs at -93.35mm. All is well! - Let me check our chamber center defintions against Richard Breedon's chamber center definition in the DDD. Here is his mail: Hello Marcus, Here are the distances from the DDD center of each CSC type to the alignment pin on the narrow end along the center line that connects the 2 pins. The dimensions come from Oleg's engineering drawings. I do not have this information yet for ME1/1. (Note that the distance from the outside of the narrow end of the frame to the alignment pin is 34.9 mm for all types. I am trying to confirm this for ME1/1 as well.) ME1/1 possibly 810.0 - 34.9 = 775.1 mm but this has to be checked ME1/2 912.1 mm ME1/3 861.6 ME2/1 988.1 ME3/1 888.1 ME4/1 788.6 ME 234/2 1655.1 COCOA: DDD: ------------------------------------------------------------------------------------------------------------ For ME+2_1: APIN apin_inner+21_14 centre X 0 prec_apin_pos fix Y 0 prec_apin_pos fix Z -968.883 prec_apin_pos fix 988.1 APIN apin_outer+21_14 centre X 0 prec_apin_pos fix Y 0 prec_apin_pos fix Z 988.187 prec_apin_pos fix 2/1/2007: --------- I copy Jim Bellinger's large DCOPS data file over to my UNIX side at FIT and gunzip it. I call it MTCC_DCOPS_centers_Oct10-Nov6_2006_from_Jim_debug.v6 and it is big (500MB)! I can open it with xemacs and if I set the lines to "truncate" it is even readable. After reading Jim's data description web page more, I find that he has made subsets for ME+2 SLM, including SLM2 (or SLM2-5). I copy that over to MTCC_DCOPS_raw_data_ME+2_slm2-5_Oct-Nov2006_from_Jim.txt on my UNIX at FIT. 2/3/2007: --------- I try to write a little FORTRAN routine dcops_format.for on my local machine that I can run within PAW so I can convert Jim Bellinger's data files into text file format that COCOA can read. After messing around I notice that the FORTRAN interpreter doesn't like to open and read .txt files with CAPS. So, use all lower case in data file names: mtcc_dcops_raw_data_mep2_slm2-5_oct-nov2006_only_fit_data.txt COMIS also can only read max. 128 columns and Jim's lines are MUCH longer.. Looks like I need to compile. So, I switch to cmsuaf where g77 exists and amazingly compiles my dcops_format.for code without any problems. It also can read very long lines as long as one reads it with READ(1,*), i.e. without any specific formatting. I can read Jim's complete file - success ! Next, I need to write COCOA-style text output. 2/5/2007: --------- Finished writing the output part of the code. When running on Jim's file I notice that there is no "perfect" B=0T event in there, i.e. all events show DCOPS with bad fits according to the status parameter. However, this is to be expected as the laser beam cannot get all the way across the disk and the DCOPS on the opposite side show bad fits with large errors. Sometimes the mean values come even back as negative. I checked this in Jim's input file, and indeed there are negative numbers in there, so it's not an error in my formatted read-in. Amazingly, COCOA seems to run right away on the newly formatted data file slm25_b0t_data_ev21_03Nov06_00_01.txt ! It even fits: 42(!) iterations with FQ=0.02 at the end... - Need to double-check that we are reading in the correct DCOPS ids: Gyongyi's web documentation claims 44 35 ME+2_reference2>IN 45 36 DCOPS_outer+22_08>OUT 46 37 DCOPS_inner+22_08>OUT 47 38 DCOPS_outer+21_04>OUT 48 39 DCOPS_inner+21_04>OUT 54 40 DCOPS_inner+21_14>IN 53 41 DCOPS_outer+21_14>IN 52 42 DCOPS_inner+22_27>IN 51 43 DCOPS_outer+22_27>IN 50 44 ME+2_reference5>IN and indeed Jim's ids are 35-44 in that sequence in the input file. In the cocoa sdf this is reflected by the following structure COPS id1 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+2_tp2/me+2_reference2>in COPS id2 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+22_08/dcops_outer+22_08>out COPS id3 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+22_08/dcops_inner+22_08>out COPS id4 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+21_04/dcops_outer+21_04>out COPS id5 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+21_04/dcops_inner+21_04>out COPS id6 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+21_14/dcops_inner+21_14>in COPS id7 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+21_14/dcops_outer+21_14>in COPS id8 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+22_27/dcops_inner+22_27>in COPS id9 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+22_27/dcops_outer+22_27>in COPS id10 OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+2_tp5/me+2_reference5>in For the other laser direction: COPS id11 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+2_tp2/me+2_reference2>in COPS id12 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+22_08/dcops_outer+22_08>out COPS id13 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+22_08/dcops_inner+22_08>out COPS id14 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+21_04/dcops_outer+21_04>out COPS id15 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+21_04/dcops_inner+21_04>out COPS id16 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+21_14/dcops_inner+21_14>in COPS id17 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+21_14/dcops_outer+21_14>in COPS id18 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+22_27/dcops_inner+22_27>in COPS id19 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+22_27/dcops_outer+22_27>in COPS id20 OCMS/slm22/me+2_tp5/las25 & OCMS/slm22/me+2_tp5/me+2_reference5>in Same sequence as for "/las22" (=las +2 pt.2)! With dcops_format.for I write the id from Jim's file also as a comment into the COCOA data file so that we always have that cross-check: Ex.: DATE: 03/11/06 00:01:48 // file id: 35 COPS id1 U 8.686 4.081 D 8.334 4.190 L 18.340 2.407 R 18.190 2.428 .... // file id: 44 COPS id20 U 9.041 2.614 D 10.275 4.549 L 17.769 4.484 R 18.351 3.857 This is fine with COCOA; it doesn't complain about the "//" in the input file. Looking at that file, I find that the correspondence appears to be fine: id1-10 is DCOPS id 35-44 in the input file and the same for id11-20. So, all appears to be correct with the ids. - Next I want to fit only 7 DCOPS if we run over B=0T data since we know that we miss the last 3 dcops in each direction. Need to modify the sdf accordingly I comment out the line in the sdf file: // COPS id8 // OCMS/slm22/me+2_tp2/las22 & OCMS/slm22/me+22_27/dcops_inner+22_27>in //U simulated_value prec_CCD //D simulated_value prec_CCD //L simulated_value prec_CCD //R simulated_value prec_CCD and the corresponding line in the data file: // file id: 42 // COPS id8 // U 11.455 427.631 // D 11.872 18.993 // L 2.300 10188.947 // R 3.770 16.166 COCOA runs ! Next, I exclude 3 DCOPS on each side, i.e. COPS id8,9,10 for las22 and COPS id11,12,13 for las 25 by commenting out the corresponding lines. COCOA still runs. Finishes after 8 iterations with terrible FQ=1.1E+08. Hmmmm... 2/14/07: -------- Before pursuing the B=0T fits further, I'm going to introduce the new official sensor names from Samir's updated authoritative list into the COCOA text file. Changes made: slm22 -> slm_p2_2 me+2 -> MEp2 in all names las22 -> laser_MEp2_2 and similar replacements of 22->2_2 and 21->2_1 OCMS -> CMS (system name) Amazingly, the code still runs (well, I checked at every single step in the replacement procedure). 2/16/07: -------- I make a new sdf with only one laser (pt.2) and 8 DCOPS and a corresponding measurement file slm25_b0t_data_ev21_03Nov06_00_01_laser2_only_10DCOPS.txt. After a couple of iterations COCOA crashes with "solve.c", line 91: singular matrix in function Usolve() "lufactor.c", line 227: singular matrix in function m_inverse() I exclude DCOPS 9 (outer DCOPS on opposite side ME2/2 chamber), which has negative numbers for the means, as well and now COCOA fits. But the output and FQ look really strange: Fit iteration 10 ... 10 Fit quality is = 1.7339622e+14 Chi2= 3.0919726e+13 / 18 dof From measurements= 3.0919726e+13 from calibrated parameters= 3856.8242 Fit iteration 10 ... 10 Fit quality is = 3.1027194e+13 Chi2= 1.0960436e+13 / 18 dof From measurements= 1.0960436e+13 from calibrated parameters= 7550.3798 Fit iteration 10 ... 10 Fit quality is = 1.0892914e+13 Chi2= 9.1981012e+12 / 18 dof From measurements= 9.1981011e+12 from calibrated parameters= 9973.4181 Fit iteration 10 ... 10 Fit quality is = 4.4128526e+12 GOOD QUALITY OF THE FIT FOR ITERATION 10 = -3.9039463e+12 < 0.1 A negative number for FQ here is really worrisome ! The fit values for chambers are also ridiculous: %%%% Optical Object: CMS/slm_p2_2/MEp2_2_08 CAL: 142 CMS/slm_p2_2/MEp2_2_08 centre_X -5225.8434 +- 100 -5225.848 +- 100 Q1 UNK: 143 CMS/slm_p2_2/MEp2_2_08 centre_Y -630776.69 +- 45858331 -93.35 +- 100 Q2 UNK: 144 CMS/slm_p2_2/MEp2_2_08 centre_Z -41919.824 +- 3011829.5 157.20274 +- 100 Q2 CAL: 145 CMS/slm_p2_2/MEp2_2_08 angles_X -30.717007 +- 100.00023 0 +- 100 Q1 UNK: 146 CMS/slm_p2_2/MEp2_2_08 angles_Y -70.24583 +- 2315805 -85 +- 100 Q2 UNK: 147 CMS/slm_p2_2/MEp2_2_08 angles_Z -3.7190203 +- 11288846 0 +- 100 Q2 Now, I notice that the values for U (CCD1) and D (CCD3) jump from 8mm for DCOPS_MEp2_reference_2_IN to 20mm for the first DCOPS on a chamber: DCOPS_outer_MEp2_2_08_OUT. Note that these DCOPS are very close together along the SLM (not more than 20cm or so). How can the laser beam be so skewed by 12mm and still thread almost all the DCOPS ?? Now, note this (old names, sorry): 44 35 ME+2_reference2>IN 45 36 DCOPS_outer+22_08>OUT 46 37 DCOPS_inner+22_08>OUT 47 38 DCOPS_outer+21_04>OUT 48 39 DCOPS_inner+21_04>OUT 54 40 DCOPS_inner+21_14>IN 53 41 DCOPS_outer+21_14>IN 52 42 DCOPS_inner+22_27>IN These are the DCOPS I am trying to fit. Note the "IN" and "OUT"'s! This is DCOPS orientation with respect to the center of the SLM. The three DCOPS on the other side (40-42) are looking into the center and the 4 DCOPS on the near side (36-39) are looking out from the center, but since they are on the other side of the center, they all have the same orientation, EXCEPT for the reference sensor at point2. It looks back in to the center! (I think all reference sensor look into the center due to symmetry of the hardware design). The consequence of that is the direction for the U (CCD1) and D (CCD3) CCDs is flipped by 180 degrees. Note that this does not affect L and R which have pixel 1 at the bottom and count up no matter if the DCOPS is facing inward or outward. Ok, so the REAL distance that ME+2_reference2>IN measures is 2048 pixels * 14 um/pixel - mean of the beam profile peak because it is flipped.. Note the "minus" sign. For example for event 21 in Jim's MTCC2 data set, this then gives BEFORE: U 8.686 4.081 D 8.334 4.190 FLIPPED: 28.672 mm - 8.686 mm = 19.986 for U U 19.986 4.081 D 20.338 4.190 This is much closer to what DCOPS_outer_MEp2_2_08_OUT measures: U 20.109 2.083 D 20.313 2.551 - By the way, note the notation ME+2_reference5>IN for the other reference sensor! It points inward just like the DCOPS on its near side, so here this is not a problem. It is only a problem for the reference sensor at pt.2! I'm glad we've been keeping that "IN" and "OUT" notation. - So, I write the above numbers directly into that data file and try again. It still won't fit! Grrr... - Take out DCOPS 8 and the long chamber on the other side and fit only 7 DCOPS. Now it crashes again on a matrix inversion error. - I make the angle of the Xlaser UNK to give it more freedom; it still crashes. The same if I make the X and Y of the Xlaser UNK... Maybe I need to go back to the original data values for the Pt.2 ref DCOPS, but turn it 180 degrees in COCOA: COPS DCOPS_MEp2_reference_2_IN ENTRY { length rightCCDYtoDowel2 bi_cops_hCCDcal cal_error_YCCD cal length leftCCDYtoDowel2 bi_cops_hCCDcal cal_error_YCCD cal length rightCCDXtoDowel2 bi_rightCCDX_disp cal_error_XCCD cal length leftCCDXtoDowel2 bi_leftCCDX2_disp cal_error_XCCD cal length upCCDXtoDowel2 bi_cops_vCCDcal cal_error_XCCD cal length upCCDYtoDowel2 bi_upCCDY_disp cal_error_YCCD cal length downCCDXtoDowel2 bi_cops_vCCDcal cal_error_XCCD cal length downCCDYtoDowel2 bi_downCCDY_disp cal_error_YCCD cal } centre X 0 prec_pos_tp cal // Y 0 prec_pos_tp cal // original Y 0.5 prec_pos_tp cal // with MTCC shims (acc. to Dave) Z 0 prec_pos_tp cal angles X 0 prec_ang_tp cal // Y 0 prec_ang_tp cal //original Y 180 prec_ang_tp cal // flip it around since it is facing IN ?! Z 0 prec_ang_tp cal Still crashes. Add DCOPS 8 back in. Now at least it runs. If I set Xlaser laser_MEp2_2 ENTRY { angle angleBetweenAxis 90. prec_ang_las unk } center X -22.500 prec_pos_tp cal Y 43.25 prec_pos_tp cal Z 0. prec_pos_tp fix angles X 0 prec_ang_tp unk //was cal Y 0 prec_ang_tp cal //was cal Z 0 prec_ang_tp cal I get Fit iteration 3 ... 3 Fit quality is = 1.2273549e+17 GOOD QUALITY OF THE FIT FOR ITERATION 3 = -1.7150689e+16 < 0.1 FITTED VALUES nEnt_unk Optical Object Parameter value (+-error) orig.val (+-error) quality %%%% Optical Object: CMS/slm_p2_2/MEp2_tp2/laser_MEp2_2 CAL: 4 CMS/slm_p2_2/MEp2_tp2/laser_MEp2_2 centre_X -29.867866 +- 20 -22.5 +- 20 Q1 CAL: 5 CMS/slm_p2_2/MEp2_tp2/laser_MEp2_2 centre_Y -27.686378 +- 20 43.25 +- 20 Q1 FIX: -1 CMS/slm_p2_2/MEp2_tp2/laser_MEp2_2 centre_Z 3.437895e-14 +- 0 0 +- 20 Q0 UNK: 6 CMS/slm_p2_2/MEp2_tp2/laser_MEp2_2 angles_X -54.174642 +- 1031.9659 0 +- 10 Q2 CAL: 7 CMS/slm_p2_2/MEp2_tp2/laser_MEp2_2 angles_Y -19.391387 +- 9.9998918 0 +- 10 Q1 CAL: 8 CMS/slm_p2_2/MEp2_tp2/laser_MEp2_2 angles_Z -4.0386343 +- 9.9996948 0 +- 10 Q1 UNK: 3 CMS/slm_p2_2/MEp2_tp2/laser_MEp2_2 angleBetweenAxis 581.91108 +- 1285.7005 90 +- 100 Q2 %%%% Optical Object: CMS/slm_p2_2/MEp2_tp2/DCOPS_MEp2_reference_2_IN CAL: 17 CMS/slm_p2_2/MEp2_tp2/DCOPS_MEp2_reference_2_IN centre_X 3.0960301e-13 +- 20 0 +- 20 Q1 CAL: 18 CMS/slm_p2_2/MEp2_tp2/DCOPS_MEp2_reference_2_IN centre_Y 0.5 +- 20 0.5 +- 20 Q1 CAL: 19 CMS/slm_p2_2/MEp2_tp2/DCOPS_MEp2_reference_2_IN centre_Z 1.7788932e-14 +- 20 0 +- 20 Q1 CAL: 20 CMS/slm_p2_2/MEp2_tp2/DCOPS_MEp2_reference_2_IN angles_X 0 +- 10 0 +- 10 Q1 CAL: 21 CMS/slm_p2_2/MEp2_tp2/DCOPS_MEp2_reference_2_IN angles_Y 0 +- 10 0 +- 10 Q1 CAL: 22 CMS/slm_p2_2/MEp2_tp2/DCOPS_MEp2_reference_2_IN angles_Z 0 +- 10 0 +- 10 Q1 So it keeps the ref. sensor in place and moves the laser around like crazy. Note the negative (!) FQ value! It's huge and negative. Really screwed up. This is using the "corrected" reference data: // file id: 35 COPS id1 // U 8.686 4.081 // D 8.334 4.190 // Subtract the width of the full CCD for U and D because this ref. sensor is // flipped by 180 degrees ! U 19.986 4.081 D 20.338 4.190 L 18.340 2.407 R 18.190 2.428 How is this program possibly able to get a converging fit on B=4T data ? I should retry all this with point 5 where the ref. DCOPS automatically points in the right direction. Maybe this data is just bad, but it doesn't seem bad. The values that the DCOPS see follow a clear trend that I can understand even without a fit. Also I need to look closer at the coordinate system for the laser on the TP. 2/21/07: -------- Here is the B=0T data file that Gyongyi used: COPS id1 U 8.694 4.200 D 8.274 130.200 L 18.242 54.600 R 18.382 63.000 COPS id2 U 20.048 8.400 D 20.244 4.200 L 19.208 47.600 R 18.172 35.000 So that looks quite similar. And here is her B=4T file: COPS id1 U 8.680 4.200 D 8.344 1.400 L 18.340 1.400 R 18.508 2.800 COPS id2 U 20.118 5.600 D 20.300 4.200 L 19.250 4.200 R 18.214 5.600 Also quite similar... ----- I try fitting the 03/11/06 00:01:48 data using only laser 5. COCOA crashes with singular matrices for 10, 9, 8, 7, 3 DCOPS. Files used for this are: sdf - ME2_SLM25_WITH_SHIMS_AND_ROTATIONS_RESTRICTED_LAS5_7DCOPS_Official_names.txt data - slm25_b0t_data_ev21_03Nov06_00_01_laser5_only_10DCOPS.txt Ok, so I am sick of this. Let's go the other way: Define a simple system with 1 laser (pt5) on 1 TP and 1 Long chamber (MEp2_2_27) for a total of 3 DCOPS.: sdf - ME2_SLM25_WITH_SHIMS_AND_ROTATIONS_RESTRICTED_LAS5_1ch_3DCOPS_Official_names.txt Now the system tree is: object system slm_line //object slm_line transfer_plate 2 Schamber 2 Lchamber object slm_line transfer_plate 1 Lchamber object transfer_plate Xlaser 2 COPS //object Schamber 2 COPS 2 APIN object Lchamber 2 COPS 2 APIN Aha, now it doesn't crash and ends with FQ=6.8 after 23 (!) Iterations. Here are the fit results: FITTED VALUES nEnt_unk Optical Object Parameter value (+-error) orig.val (+-error) quality %%%% Optical Object: CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 CAL: 4 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_X -22.5 +- 20 -22.5 +- 20 Q1 CAL: 5 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_Y 43.25 +- 20 43.25 +- 20 Q1 FIX: -1 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_Z -3.4482751e-13 +- 0 0 +- 20 Q0 UNK: 6 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_X -0.091599509 +- 1029.5787 0 +- 10 Q2 UNK: 7 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Y -3.4870994 +- 55630.464 0 +- 10 Q2 CAL: 8 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Z -1.1679575e-05 +- 9.9997625 0 +- 10 Q1 UNK: 3 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angleBetweenAxis 90.444881 +- 5921.1487 90 +- 100 Q2 %%%% Optical Object: CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN CAL: 17 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN centre_X 9.0212405e-14 +- 20 0 +- 20 Q1 CAL: 18 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN centre_Y 0.5 +- 20 0.5 +- 20 Q1 CAL: 19 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN centre_Z -2.5792469e-14 +- 20 0 +- 20 Q1 CAL: 20 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN angles_X 0 +- 10 0 +- 10 Q1 CAL: 21 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN angles_Y 0 +- 10 0 +- 10 Q1 CAL: 22 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN angles_Z 0 +- 10 0 +- 10 Q1 CAL: 9 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN rightCCDYtoDowel2 57.579 +- 65 57.579 +- 65 Q1 CAL: 10 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN leftCCDYtoDowel2 57.579 +- 65 57.579 +- 65 Q1 CAL: 11 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN rightCCDXtoDowel2 -8.171 +- 65 -8.171 +- 65 Q1 CAL: 12 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN leftCCDXtoDowel2 -45 +- 65 -45 +- 65 Q1 CAL: 13 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN upCCDXtoDowel2 -36.829 +- 65 -36.829 +- 65 Q1 CAL: 14 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN upCCDYtoDowel2 65.75 +- 65 65.75 +- 65 Q1 CAL: 15 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN downCCDXtoDowel2 -36.829 +- 65 -36.829 +- 65 Q1 CAL: 16 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN downCCDYtoDowel2 20.75 +- 65 20.75 +- 65 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT CAL: 31 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT centre_X -465 +- 20 -465 +- 20 Q1 CAL: 32 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT centre_Y 82.275 +- 20 82.275 +- 20 Q1 CAL: 33 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT centre_Z -31.902 +- 20 -31.902 +- 20 Q1 CAL: 34 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT angles_X 270 +- 10 270 +- 10 Q1 CAL: 35 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT angles_Y 180 +- 10 180 +- 10 Q1 CAL: 36 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT angles_Z 180 +- 10 180 +- 10 Q1 CAL: 23 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT rightCCDYtoDowel2 57.579 +- 65 57.579 +- 65 Q1 CAL: 24 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT leftCCDYtoDowel2 57.579 +- 65 57.579 +- 65 Q1 CAL: 25 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT rightCCDXtoDowel2 -8.171 +- 65 -8.171 +- 65 Q1 CAL: 26 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT leftCCDXtoDowel2 -45 +- 65 -45 +- 65 Q1 CAL: 27 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT upCCDXtoDowel2 -36.829 +- 65 -36.829 +- 65 Q1 CAL: 28 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT upCCDYtoDowel2 65.75 +- 65 65.75 +- 65 Q1 CAL: 29 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT downCCDXtoDowel2 -36.829 +- 65 -36.829 +- 65 Q1 CAL: 30 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT downCCDYtoDowel2 20.75 +- 65 20.75 +- 65 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_tp5 FIX: -1 CMS/slm_p2_2/MEp2_tp5 centre_X 7175.833 +- 0 7175.833 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_tp5 centre_Y -43.25 +- 0 -43.25 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_tp5 centre_Z -22.5 +- 0 -22.5 +- 100 Q0 CAL: 0 CMS/slm_p2_2/MEp2_tp5 angles_X 0 +- 697.64234 0 +- 698 Q1 CAL: 1 CMS/slm_p2_2/MEp2_tp5 angles_Y 89.999996 +- 698 90 +- 698 Q1 CAL: 2 CMS/slm_p2_2/MEp2_tp5 angles_Z 0 +- 697.64234 0 +- 698 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN CAL: 51 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN centre_X 280.19898 +- 144 280.199 +- 144 Q1 CAL: 52 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN centre_Y 50.100022 +- 538 50.1 +- 538 Q1 CAL: 53 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN centre_Z 1657.151 +- 100 1657.151 +- 100 Q1 CAL: 54 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN angles_X 0 +- 2 0 +- 2 Q1 CAL: 55 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN angles_Y -5 +- 2 185 +- 2 Q1 CAL: 56 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN angles_Z 5.2659502e-07 +- 1.9999991 0 +- 2 Q1 CAL: 43 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN rightCCDYtoDowel2 57.579489 +- 64.999812 57.579 +- 65 Q1 CAL: 44 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN leftCCDYtoDowel2 57.578513 +- 64.999814 57.579 +- 65 Q1 CAL: 45 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN rightCCDXtoDowel2 -8.170996 +- 65 -8.171 +- 65 Q1 CAL: 46 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN leftCCDXtoDowel2 -45.000004 +- 65 -45 +- 65 Q1 CAL: 47 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN upCCDXtoDowel2 -36.612009 +- 46.072392 -36.829 +- 65 Q1 CAL: 48 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN upCCDYtoDowel2 65.750002 +- 65 65.75 +- 65 Q1 CAL: 49 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN downCCDXtoDowel2 -37.045986 +- 46.073495 -36.829 +- 65 Q1 CAL: 50 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN downCCDYtoDowel2 20.749998 +- 65 20.75 +- 65 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN CAL: 65 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN centre_X 19.922023 +- 144 19.922 +- 144 Q1 CAL: 66 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN centre_Y 46.099978 +- 538 46.1 +- 538 Q1 CAL: 67 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN centre_Z -1317.833 +- 100 -1317.833 +- 100 Q1 CAL: 68 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN angles_X 0 +- 2 0 +- 2 Q1 CAL: 69 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN angles_Y -5 +- 2 185 +- 2 Q1 CAL: 70 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN angles_Z -6.0287156e-08 +- 1.999999 0 +- 2 Q1 CAL: 57 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN rightCCDYtoDowel2 57.578512 +- 64.999813 57.579 +- 65 Q1 CAL: 58 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN leftCCDYtoDowel2 57.579487 +- 64.999814 57.579 +- 65 Q1 CAL: 59 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN rightCCDXtoDowel2 -8.171004 +- 65 -8.171 +- 65 Q1 CAL: 60 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN leftCCDXtoDowel2 -44.999996 +- 65 -45 +- 65 Q1 CAL: 61 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN upCCDXtoDowel2 -36.854178 +- 46.004452 -36.829 +- 65 Q1 CAL: 62 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN upCCDYtoDowel2 65.75 +- 65 65.75 +- 65 Q1 CAL: 63 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN downCCDXtoDowel2 -36.803826 +- 46.003345 -36.829 +- 65 Q1 CAL: 64 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN downCCDYtoDowel2 20.75 +- 65 20.75 +- 65 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 centre_X 9.3284841e-16 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 centre_Y -5.9733559e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 centre_Z -1655.191 +- 0 -1655.191 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 angles_X 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 angles_Y 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 angles_Z 0 +- 0 0 +- 100 Q0 %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 centre_X 1.5978657e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 centre_Y 5.9733537e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 centre_Z 1655.191 +- 0 1655.191 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 angles_X 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 angles_Y 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 angles_Z 0 +- 0 0 +- 100 Q0 %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27 CAL: 37 CMS/slm_p2_2/MEp2_2_27 centre_X 5225.848 +- 100 5225.848 +- 100 Q1 UNK: 38 CMS/slm_p2_2/MEp2_2_27 centre_Y -91.34736 +- 899.22331 -93.35 +- 100 Q2 UNK: 39 CMS/slm_p2_2/MEp2_2_27 centre_Z 42.545318 +- 108352.96 157.20274 +- 100 Q2 CAL: 40 CMS/slm_p2_2/MEp2_2_27 angles_X 0.00010203609 +- 99.998187 0 +- 100 Q1 UNK: 41 CMS/slm_p2_2/MEp2_2_27 angles_Y 81.519244 +- 55625.95 85 +- 100 Q2 CAL: 42 CMS/slm_p2_2/MEp2_2_27 angles_Z 7.863816e-05 +- 99.998923 0 +- 100 Q1 %%%% Optical Object: CMS/slm_p2_2 FIX: -1 CMS/slm_p2_2 centre_X 289.7777 +- 0 289.7777 +- 1 Q0 FIX: -1 CMS/slm_p2_2 centre_Y -77.6457 +- 0 -77.6457 +- 1 Q0 FIX: -1 CMS/slm_p2_2 centre_Z 7998 +- 0 7998 +- 1 Q0 FIX: -1 CMS/slm_p2_2 angles_X -90 +- 0 -90 +- 1 Q0 FIX: -1 CMS/slm_p2_2 angles_Y 0 +- 0 0 +- 1 Q0 FIX: -1 CMS/slm_p2_2 angles_Z 75 +- 0 255 +- 1 Q0 and in global coordinates: OPTICALOBJECT list size 11 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 CENTRE GLOBAL: (-1.567466,-7.0089736,7.998) CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 ROT_MATRIX GLOBAL: xx=-0.97987978 xy=0.00028468413 xz=-0.19958843 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 ROT_MATRIX GLOBAL: yx=0.19958863 yy=0.0014389484 yz=-0.97987872 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 ROT_MATRIX GLOBAL: zx=8.2415366e-06 zy=-0.99999892 zz=-0.0014668162 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN CENTRE GLOBAL: (-1.5891979,-7.0031448,8.04075) CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN ROT_MATRIX GLOBAL: xx=-0.96592581 xy=-3.4203819e-05 xz=-0.25881911 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN ROT_MATRIX GLOBAL: yx=0.25881911 yy=-0.00012765037 yz=-0.9659258 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN ROT_MATRIX GLOBAL: zx=-4.4405097e-12 zy=-0.99999999 zz=0.00013215339 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT CENTRE GLOBAL: (-1.1317883,-7.0926911,7.9589708) CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT ROT_MATRIX GLOBAL: xx=-0.96592581 xy=-0.25881911 xz=3.4203819e-05 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT ROT_MATRIX GLOBAL: yx=0.25881911 yy=-0.9659258 yz=0.00012765037 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT ROT_MATRIX GLOBAL: zx=-4.4403872e-12 zy=0.00013215339 zz=0.99999999 CMS/slm_p2_2/MEp2_tp5 CENTRE GLOBAL: (-1.5891979,-7.0031447,8.04125) CMS/slm_p2_2/MEp2_tp5 ROT_MATRIX GLOBAL: xx=-0.96592581 xy=-3.4203819e-05 xz=-0.25881911 CMS/slm_p2_2/MEp2_tp5 ROT_MATRIX GLOBAL: yx=0.25881911 yy=-0.00012765037 yz=-0.9659258 CMS/slm_p2_2/MEp2_tp5 ROT_MATRIX GLOBAL: zx=-4.4405097e-12 zy=-0.99999999 zz=0.00013215339 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN CENTRE GLOBAL: (-1.4882106,-6.7510614,8.039248) CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN ROT_MATRIX GLOBAL: xx=0.97985769 xy=1.6216517e-07 xz=0.19969705 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN ROT_MATRIX GLOBAL: yx=-0.19969705 yy=-4.4466564e-07 yz=0.97985769 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN ROT_MATRIX GLOBAL: zx=2.476972e-07 zy=-1 zz=-4.0332514e-07 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN CENTRE GLOBAL: (-0.89184616,-3.8248653,8.0432468) CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN ROT_MATRIX GLOBAL: xx=0.97985769 xy=1.5198838e-07 xz=0.19969705 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN ROT_MATRIX GLOBAL: yx=-0.19969705 yy=-4.4350267e-07 yz=0.97985769 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN ROT_MATRIX GLOBAL: zx=2.3749315e-07 zy=-1 zz=-4.0421787e-07 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 CENTRE GLOBAL: (-0.83375038,-3.4919506,8.0893467) CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 ROT_MATRIX GLOBAL: xx=-0.99353378 xy=1.5303378e-07 xz=-0.11353692 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 ROT_MATRIX GLOBAL: yx=0.11353692 yy=-4.4362214e-07 yz=-0.99353378 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 ROT_MATRIX GLOBAL: zx=-2.0241172e-07 zy=-1 zz=4.2337859e-07 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 CENTRE GLOBAL: (-1.209601,-6.780927,8.0893481) CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 ROT_MATRIX GLOBAL: xx=-0.99353378 xy=1.5303378e-07 xz=-0.11353692 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 ROT_MATRIX GLOBAL: yx=0.11353692 yy=-4.4362214e-07 yz=-0.99353378 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 ROT_MATRIX GLOBAL: zx=-2.0241172e-07 zy=-1 zz=4.2337859e-07 CMS/slm_p2_2/MEp2_2_27 CENTRE GLOBAL: (-1.0216757,-5.1364388,8.0893474) CMS/slm_p2_2/MEp2_2_27 ROT_MATRIX GLOBAL: xx=-0.99353378 xy=1.5303378e-07 xz=-0.11353692 CMS/slm_p2_2/MEp2_2_27 ROT_MATRIX GLOBAL: yx=0.11353692 yy=-4.4362214e-07 yz=-0.99353378 CMS/slm_p2_2/MEp2_2_27 ROT_MATRIX GLOBAL: zx=-2.0241172e-07 zy=-1 zz=4.2337859e-07 CMS/slm_p2_2 CENTRE GLOBAL: (0.2897777,-0.0776457,7.998) CMS/slm_p2_2 ROT_MATRIX GLOBAL: xx=-0.25881905 xy=5.9143945e-17 xz=0.96592583 CMS/slm_p2_2 ROT_MATRIX GLOBAL: yx=-0.96592583 yy=-1.5847572e-17 yz=-0.25881905 CMS/slm_p2_2 ROT_MATRIX GLOBAL: zx=0 zy=-1 zz=6.1230318e-17 CMS CENTRE GLOBAL: (0,0,0) CMS ROT_MATRIX GLOBAL: xx=1 xy=0 xz=0 CMS ROT_MATRIX GLOBAL: yx=0 yy=1 yz=0 CMS ROT_MATRIX GLOBAL: zx=0 zy=0 zz=1 Observations: - After fitting, the laser has been rotated by -3.5 deg around Y in TP coordinates (angles_Y) According to Gyongyi's fig.6 in her note, this is the same coordinate as Y_SLM or Z_CMS. This would imply that we are seeing a rotation that would be picked up by the U and D CCDs. Note sure if this is really bourne out by the data: // file id: 42 Chamber Inner DCOPS COPS id18 U 10.341 1.744 D 10.291 1.780 L 17.038 2.017 R 17.337 1.965 // file id: 43 Chamber Outer DCOPS COPS id19 U 10.428 3.909 D 10.864 4.066 L 16.628 3.037 R 16.942 4.233 // file id: 44 Pt.5 ref COPS id20 U 9.041 2.614 D 10.275 4.549 L 17.769 4.484 R 18.351 3.857 Ref. DCOPS has an average of 9.5mm and the chamber DCOPS have 10.2mm, so a change of 0.7mm. Over a : distance of say 3m, this is 0.7/3000=0.23 mrad = 0.013 deg. Seems too small or rather 3.5 deg seems too big. - A rotation due to dipping the laser up or down would be a rotation around X_TP according to Gyongyi's fig.6.. The fit comes back with only 0.1 deg here. - Note that laser rotations around X_TP and Y_TP were set to UNK in this run, but rotations around Z_TP were set to CAL. Out of curiosity, let's set that to UNK and see what we get. Here is the answer: After 19 fits FQ=38577401 (bad!): CAL: 4 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_X -22.5 +- 20 -22.5 +- 20 Q1 CAL: 5 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_Y 43.25 +- 20 43.25 +- 20 Q1 FIX: -1 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_Z -3.8478638e-13 +- 0 0 +- 20 Q0 UNK: 6 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_X 14.765887 +- 1499.7587 0 +- 10 Q2 UNK: 7 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Y 7.9025596 +- 1336.0971 0 +- 10 Q2 UNK: 8 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Z -0.31609732 +- 1836.5306 0 +- 10 Q2 UNK: 3 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angleBetweenAxis 89.235694 +- 3509.062 90 +- 10 So, setting this to CAL appears to be a good thing. - The angle between the 2 crosshair lasers gets fitted to 90.4 deg -> seems ok - All the DCOPS values come back as nominal which is what we expect since we expect the chamber to be moved and not the DCOPS on the chamber. - TP is held in place by making all pos. FIX and rots CAL - The one chamber gets fitted to %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27 CAL: 37 CMS/slm_p2_2/MEp2_2_27 centre_X 5225.848 +- 100 5225.848 +- 100 Q1 UNK: 38 CMS/slm_p2_2/MEp2_2_27 centre_Y -91.34736 +- 899.22331 -93.35 +- 100 Q2 UNK: 39 CMS/slm_p2_2/MEp2_2_27 centre_Z 42.545318 +- 108352.96 157.20274 +- 100 Q2 CAL: 40 CMS/slm_p2_2/MEp2_2_27 angles_X 0.00010203609 +- 99.998187 0 +- 100 Q1 UNK: 41 CMS/slm_p2_2/MEp2_2_27 angles_Y 81.519244 +- 55625.95 85 +- 100 Q2 CAL: 42 CMS/slm_p2_2/MEp2_2_27 angles_Z 7.863816e-05 +- 99.998923 0 +- 100 Q1 That means it shifts in Y_SLM=Z_CMS by 2mm and in Z_SLM=R-phi_CMS by a whopping 115 mm !!! Wow. There are too many darn variables in this! Laser rotations plus chamber rots etc. Ok, seems like I haven't restricted chamber rots properly: UNK: 41 CMS/slm_p2_2/MEp2_2_27 angles_Y 81.519244 +- 55625.95 85 +- 100 Q2 Maybe this should be CAL, too ! I set it to CAL and rerun: AHA ! After 2 iterations I get FQ=0.00017. Cool. 2 Fit quality is = 0.00017173357 GOOD QUALITY OF THE FIT FOR ITERATION 2 = 0.00017173357 < 0.1 FITTED VALUES nEnt_unk Optical Object Parameter value (+-error) orig.val (+-error) quality %%%% Optical Object: CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 CAL: 4 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_X -22.5 +- 20 -22.5 +- 20 Q1 CAL: 5 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_Y 43.25 +- 20 43.25 +- 20 Q1 FIX: -1 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_Z -1.4850706e-14 +- 0 0 +- 20 Q0 UNK: 6 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_X -0.084324177 +- 1029.7356 0 +- 10 Q2 UNK: 7 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Y -0.0063403921 +- 708.73883 0 +- 10 Q2 CAL: 8 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Z -1.1727648e-05 +- 9.9997609 0 +- 10 Q1 UNK: 3 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angleBetweenAxis 90.476827 +- 1766.7487 90 +- 100 Q2 .. %%%% Optical Object: CMS/slm_p2_2/MEp2_tp5 FIX: -1 CMS/slm_p2_2/MEp2_tp5 centre_X 7175.833 +- 0 7175.833 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_tp5 centre_Y -43.25 +- 0 -43.25 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_tp5 centre_Z -22.5 +- 0 -22.5 +- 100 Q0 CAL: 0 CMS/slm_p2_2/MEp2_tp5 angles_X 6.220897e-06 +- 698 0 +- 698 Q1 CAL: 1 CMS/slm_p2_2/MEp2_tp5 angles_Y 90 +- 698 90 +- 698 Q1 CAL: 2 CMS/slm_p2_2/MEp2_tp5 angles_Z -6.3045603e-06 +- 698 0 +- 698 Q1 ... %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27 CAL: 37 CMS/slm_p2_2/MEp2_2_27 centre_X 5225.848 +- 100 5225.848 +- 100 Q1 UNK: 38 CMS/slm_p2_2/MEp2_2_27 centre_Y -91.291959 +- 652.97361 -93.35 +- 100 Q2 UNK: 39 CMS/slm_p2_2/MEp2_2_27 centre_Z 160.85508 +- 256.46633 157.20274 +- 100 Q2 CAL: 40 CMS/slm_p2_2/MEp2_2_27 angles_X 0.00010133113 +- 99.998124 0 +- 100 Q1 CAL: 41 CMS/slm_p2_2/MEp2_2_27 angles_Y 84.999998 +- 99.999686 85 +- 100 Q1 CAL: 42 CMS/slm_p2_2/MEp2_2_27 angles_Z -1.5271724e-06 +- 99.999685 0 +- 100 Q1 %%%% Optical Object: CMS/slm_p2_2 FIX: -1 CMS/slm_p2_2 centre_X 289.7777 +- 0 289.7777 +- 1 Q0 FIX: -1 CMS/slm_p2_2 centre_Y -77.6457 +- 0 -77.6457 +- 1 Q0 FIX: -1 CMS/slm_p2_2 centre_Z 7998 +- 0 7998 +- 1 Q0 FIX: -1 CMS/slm_p2_2 angles_X -90 +- 0 -90 +- 1 Q0 FIX: -1 CMS/slm_p2_2 angles_Y 0 +- 0 0 +- 1 Q0 FIX: -1 CMS/slm_p2_2 angles_Z 75 +- 0 255 +- 1 Q0 Now lookey here!! - laser dip angles come back with really small values and much more consistent with what I expect from my calculation above UNK: 6 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_X -0.084324177 +- 1029.7356 0 +- 10 Q2 UNK: 7 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Y -0.0063403921 +- 708.73883 0 +- 10 Q2 - chamber is barely rotated: CAL: 40 CMS/slm_p2_2/MEp2_2_27 angles_X 0.00010133113 +- 99.998124 0 +- 100 Q1 CAL: 41 CMS/slm_p2_2/MEp2_2_27 angles_Y 84.999998 +- 99.999686 85 +- 100 Q1 CAL: 42 CMS/slm_p2_2/MEp2_2_27 angles_Z -1.5271724e-06 +- 99.999685 0 +- 100 Q1 but shifted by a few mm: UNK: 38 CMS/slm_p2_2/MEp2_2_27 centre_Y -91.291959 +- 652.97361 -93.35 +- 100 Q2 UNK: 39 CMS/slm_p2_2/MEp2_2_27 centre_Z 160.85508 +- 256.46633 157.20274 +- 100 Q2 That's more like it !! So the lesson is that one should not only restrict X and Z rotation, but also the Y rotation that give the 5 and 10 deg offset. Ok, now let's try to build this up systematically: Add the next chamber in the SLM, the near short chamber: object system slm_line object slm_line transfer_plate 1 Schamber 1 Lchamber object transfer_plate Xlaser 2 COPS object Schamber 2 COPS 2 APIN object Lchamber 2 COPS 2 APIN Now it crashes again: [hohlmann@cmswn083 Analysis]$ cocoa ME2_SLM25_WITH_SHIMS_AND_ROTATIONS_RESTRICTED_LAS5_2ch_5DCOPS_Official_names.txt WARNING: in line 331 number of repeated OpticalObjects = 1. Please avoid the numbering WARNING: in line 331 number of repeated OpticalObjects = 1. Please avoid the numbering TIME:ENDED_READING : 800000 0.02 0 0x87be2e0$$$$$$$$$$$$$$$ moreData Sets 1 4976.7 adding chi2meas 2.47675e+07 id16: 0 (mm)R: 9.61 S: 0 Diff= 9.61 3785.17 adding chi2meas 3.9095e+07 id16: 1 (mm)R: 10.36 S: 0 Diff= 10.36 503.709 adding chi2meas 3.93487e+07 id16: 2 (mm)R: 18.81 S: 0 Diff= 18.81 9780.06 adding chi2meas 1.34998e+08 id16: 3 (mm)R: 19.032 S: 0 Diff= 19.032 6338.32 adding chi2meas 1.75173e+08 id17: 0 (mm)R: 11.447 S: 0 Diff= 11.447 3905.75 adding chi2meas 1.90427e+08 id17: 1 (mm)R: 11.272 S: 0 Diff= 11.272 6626.21 adding chi2meas 2.34334e+08 id17: 2 (mm)R: 17.142 S: 0 Diff= 17.142 4237.94 adding chi2meas 2.52294e+08 id17: 3 (mm)R: 17.134 S: 0 Diff= 17.134 -2286.79 adding chi2meas 2.57524e+08 id18: 0 (mm)R: 10.341 S: 14.3292 Diff= -3.98816 -2268.63 adding chi2meas 2.6267e+08 id18: 1 (mm)R: 10.291 S: 14.3292 Diff= -4.03816 3326.23 adding chi2meas 2.73734e+08 id18: 2 (mm)R: 17.038 S: 10.329 Diff= 6.709 3566.41 adding chi2meas 2.86453e+08 id18: 3 (mm)R: 17.337 S: 10.329 Diff= 7.008 -997.9 adding chi2meas 2.87449e+08 id19: 0 (mm)R: 10.428 S: 14.3288 Diff= -3.90079 -852.138 adding chi2meas 2.88175e+08 id19: 1 (mm)R: 10.864 S: 14.3288 Diff= -3.46479 756.997 adding chi2meas 2.88748e+08 id19: 2 (mm)R: 16.628 S: 14.329 Diff= 2.299 617.293 adding chi2meas 2.89129e+08 id19: 3 (mm)R: 16.942 S: 14.329 Diff= 2.613 -2022.95 adding chi2meas 2.93222e+08 id20: 0 (mm)R: 9.041 S: 14.329 Diff= -5.288 -891.185 adding chi2meas 2.94016e+08 id20: 1 (mm)R: 10.275 S: 14.329 Diff= -4.054 655.665 adding chi2meas 2.94446e+08 id20: 2 (mm)R: 17.769 S: 14.829 Diff= 2.94 913.145 adding chi2meas 2.9528e+08 id20: 3 (mm)R: 18.351 S: 14.829 Diff= 3.522 TIME:MAT_MEAS_FILLED: 840000 0.03 TIME:MAT_CAL_FILLED : 840000 0 TIME:BEFORE_INVERSE : 850000 0.01 "solve.c", line 91: singular matrix in function Usolve() "lufactor.c", line 227: singular matrix in function m_inverse() Note the different S values for short and long chamber. Especially that it S=0 for the inner chamber. Why? - I get rid of the WARNING: WARNING: in line 331 number of repeated OpticalObjects = 1. Please avoid the numbering by defining: object slm_line transfer_plate Schamber Lchamber instead of object slm_line transfer_plate 1 Schamber 1 Lchamber But it still crashes. 2/23/2007 --------- - I reverse the order of long and short chamber in the sdf: object system slm_line //object slm_line transfer_plate 2 Schamber 2 Lchamber //object slm_line transfer_plate Schamber Lchamber object slm_line transfer_plate Lchamber Schamber That doesn't make a difference. Still crashes. - Change the last laser angle to unk. Xlaser laser_MEp2_5 ENTRY { angle angleBetweenAxis 90. prec_ang_las unk } center angles Z 0 prec_ang_tp unk //was cal Still crashes. - I add the short chamber on the far end back in: Schamber MEp2_1_04 Now it reads in values: -44.0379 adding chi2meas 1939.33 id14: 0 (mm)R: 9.896 S: 14.3386 Diff= -4.44258 -14.3892 adding chi2meas 2146.38 id14: 1 (mm)R: 10.544 S: 14.3386 Diff= -3.79458 -533.495 adding chi2meas 286763 id14: 2 (mm)R: 6.032 S: 9.329 Diff= -3.297 -1429.97 adding chi2meas 2.33158e+06 id14: 3 (mm)R: 6.552 S: 9.329 Diff= -2.777 5884.42 adding chi2meas 3.69579e+07 id15: 0 (mm)R: 10.233 S: 0 Diff= 10.233 6030.7 adding chi2meas 7.33272e+07 id15: 1 (mm)R: 8.841 S: 0 Diff= 8.841 8037.39 adding chi2meas 1.37927e+08 id15: 2 (mm)R: 13.326 S: 0 Diff= 13.326 7465.1 adding chi2meas 1.93655e+08 id15: 3 (mm)R: 13.691 S: 0 Diff= 13.691 4976.7 adding chi2meas 2.18422e+08 id16: 0 (mm)R: 9.61 S: 0 Diff= 9.61 3785.17 adding chi2meas 2.3275e+08 id16: 1 (mm)R: 10.36 S: 0 Diff= 10.36 503.709 adding chi2meas 2.33003e+08 id16: 2 (mm)R: 18.81 S: 0 Diff= 18.81 9780.06 adding chi2meas 3.28653e+08 id16: 3 (mm)R: 19.032 S: 0 Diff= 19.032 -1601.1 adding chi2meas 3.31217e+08 id17: 0 (mm)R: 11.447 S: 14.3386 Diff= -2.89158 -1062.57 adding chi2meas 3.32346e+08 id17: 1 (mm)R: 11.272 S: 14.3386 Diff= -3.06658 3020.1 adding chi2meas 3.41467e+08 id17: 2 (mm)R: 17.142 S: 9.329 Diff= 7.813 1930.5 adding chi2meas 3.45193e+08 id17: 3 (mm)R: 17.134 S: 9.329 Diff= 7.805 -2286.79 adding chi2meas 3.50423e+08 id18: 0 (mm)R: 10.341 S: 14.3292 Diff= -3.98816 -2268.63 adding chi2meas 3.5557e+08 id18: 1 (mm)R: 10.291 S: 14.3292 Diff= -4.03816 3326.23 adding chi2meas 3.66633e+08 id18: 2 (mm)R: 17.038 S: 10.329 Diff= 6.709 3566.41 adding chi2meas 3.79353e+08 id18: 3 (mm)R: 17.337 S: 10.329 Diff= 7.008 -997.9 adding chi2meas 3.80348e+08 id19: 0 (mm)R: 10.428 S: 14.3288 Diff= -3.90079 -852.138 adding chi2meas 3.81075e+08 id19: 1 (mm)R: 10.864 S: 14.3288 Diff= -3.46479 756.997 adding chi2meas 3.81648e+08 id19: 2 (mm)R: 16.628 S: 14.329 Diff= 2.299 617.293 adding chi2meas 3.82029e+08 id19: 3 (mm)R: 16.942 S: 14.329 Diff= 2.613 -2022.95 adding chi2meas 3.86121e+08 id20: 0 (mm)R: 9.041 S: 14.329 Diff= -5.288 -891.185 adding chi2meas 3.86915e+08 id20: 1 (mm)R: 10.275 S: 14.329 Diff= -4.054 655.665 adding chi2meas 3.87345e+08 id20: 2 (mm)R: 17.769 S: 14.829 Diff= 2.94 913.145 adding chi2meas 3.88179e+08 id20: 3 (mm)R: 18.351 S: 14.829 Diff= 3.522 and begins fitting Fit iteration 0 ... 0 Fit quality is = 3.82764e+11 0 Fit quality predicted improvement in distance to minimum is = 6.50709e+07 0 Fit quality derivative at 0 = -7.65528e+11 Note that now the weird S values are in DCOPS id15 and id16 whereas before with only one short chamber they were in id16 and id17. - When I set Schamber MEp2_1_04 angles X 0 prec_chamber_ang_det cal Y -80 prec_chamber_ang_det unk Z 0 prec_chamber_ang_det cal it crashes after about one iteration, BUT with angles X 0 prec_chamber_ang_det cal Y -80 prec_chamber_ang_det cal Z 0 prec_chamber_ang_det cal it fits and doesn't crash: Results: Fit iteration 3 ... 3 Fit quality is = 0.03481277 GOOD QUALITY OF THE FIT FOR ITERATION 3 = 0.03481277 < 0.1 Not bad at all, so this fits 3 chambers now! FITTED VALUES nEnt_unk Optical Object Parameter value (+-error) orig.val (+-error) quality %%%% Optical Object: CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 CAL: 4 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_X -23.137983 +- 18.579452 -22.5 +- 20 Q1 CAL: 5 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_Y 42.810046 +- 18.584221 43.25 +- 20 Q1 FIX: -1 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_Z 1.9329808e-14 +- 0 0 +- 20 Q0 UNK: 6 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_X -0.084340938 +- 1028.8778 0 +- 10 Q2 UNK: 7 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Y -0.0059342079 +- 708.73878 0 +- 10 Q2 CAL: 8 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Z 4.2152483e-05 +- 9.99951 0 +- 10 Q1 UNK: 3 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angleBetweenAxis 90.390277 +- 1116.0697 90 +- 100 Q2 %%%% Optical Object: CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN CAL: 17 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN centre_X 0.63798338 +- 18.579452 0 +- 20 Q1 CAL: 18 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN centre_Y 0.93960279 +- 18.586565 0.5 +- 20 Q1 CAL: 19 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN centre_Z -0.00058256026 +- 19.999997 0 +- 20 Q1 CAL: 20 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN angles_X -6.7930022e-07 +- 10 0 +- 10 Q1 CAL: 21 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN angles_Y 2.7026413e-07 +- 10 0 +- 10 Q1 CAL: 22 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN angles_Z -0.0005798926 +- 9.999382 0 +- 10 Q1 CAL: 9 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN rightCCDYtoDowel2 60.317901 +- 30.378372 57.579 +- 65 Q1 CAL: 10 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN leftCCDYtoDowel2 59.485243 +- 33.3521 57.579 +- 65 Q1 CAL: 11 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN rightCCDXtoDowel2 -8.1894389 +- 64.998849 -8.171 +- 65 Q1 CAL: 12 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN leftCCDXtoDowel2 -45.012833 +- 64.998915 -45 +- 65 Q1 CAL: 13 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN upCCDXtoDowel2 -32.824123 +- 24.208748 -36.829 +- 65 Q1 CAL: 14 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN upCCDYtoDowel2 65.750043 +- 65 65.75 +- 65 Q1 CAL: 15 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN downCCDXtoDowel2 -34.064703 +- 24.408386 -36.829 +- 65 Q1 CAL: 16 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN downCCDYtoDowel2 20.75003 +- 65 20.75 +- 65 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_tp5 FIX: -1 CMS/slm_p2_2/MEp2_tp5 centre_X 7175.833 +- 0 7175.833 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_tp5 centre_Y -43.25 +- 0 -43.25 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_tp5 centre_Z -22.5 +- 0 -22.5 +- 100 Q0 CAL: 0 CMS/slm_p2_2/MEp2_tp5 angles_X 0 +- 697.99999 0 +- 698 Q1 CAL: 1 CMS/slm_p2_2/MEp2_tp5 angles_Y 89.999607 +- 697.99991 90 +- 698 Q1 CAL: 2 CMS/slm_p2_2/MEp2_tp5 angles_Z 0 +- 698 0 +- 698 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 centre_X 3.2859879e-14 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 centre_Y -4.039846e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 centre_Z -1655.191 +- 0 -1655.191 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 angles_X 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 angles_Y 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 angles_Z 0 +- 0 0 +- 100 Q0 %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 centre_X -1.0997376e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 centre_Y 4.039852e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 centre_Z 1655.191 +- 0 1655.191 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 angles_X 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 angles_Y 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 angles_Z 0 +- 0 0 +- 100 Q0 %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27 CAL: 37 CMS/slm_p2_2/MEp2_2_27 centre_X 5225.848 +- 100 5225.848 +- 100 Q1 UNK: 38 CMS/slm_p2_2/MEp2_2_27 centre_Y -91.733375 +- 652.2979 -93.35 +- 100 Q2 UNK: 39 CMS/slm_p2_2/MEp2_2_27 centre_Z 161.49349 +- 256.35984 157.20274 +- 100 Q2 CAL: 40 CMS/slm_p2_2/MEp2_2_27 angles_X 6.8986745e-05 +- 99.996742 0 +- 100 Q1 CAL: 41 CMS/slm_p2_2/MEp2_2_27 angles_Y 85.000007 +- 99.999972 85 +- 100 Q1 CAL: 42 CMS/slm_p2_2/MEp2_2_27 angles_Z -9.8815737e-06 +- 99.999717 0 +- 100 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 centre_X 1.4906743e-15 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 centre_Y 7.3143362e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 centre_Z -968.883 +- 0 -968.883 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 angles_X 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 angles_Y 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 angles_Z 0 +- 0 0 +- 100 Q0 %%%% Optical Object: CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 centre_X 2.5728663e-14 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 centre_Y 5.1693203e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 centre_Z 988.187 +- 0 988.187 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 angles_X 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 angles_Y 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 angles_Z 0 +- 0 0 +- 100 Q0 %%%% Optical Object: CMS/slm_p2_2/MEp2_1_14 CAL: 71 CMS/slm_p2_2/MEp2_1_14 centre_X 2371.279 +- 100 2371.279 +- 100 Q1 UNK: 72 CMS/slm_p2_2/MEp2_1_14 centre_Y -91.537789 +- 1236.663 -93.35 +- 100 Q2 UNK: 73 CMS/slm_p2_2/MEp2_1_14 centre_Z 121.36203 +- 499.25084 118.16258 +- 100 Q2 CAL: 74 CMS/slm_p2_2/MEp2_1_14 angles_X 1.77014e-05 +- 99.994517 0 +- 100 Q1 CAL: 75 CMS/slm_p2_2/MEp2_1_14 angles_Y 79.999995 +- 99.999998 80 +- 100 Q1 CAL: 76 CMS/slm_p2_2/MEp2_1_14 angles_Z -3.2147925e-06 +- 99.999998 0 +- 100 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 centre_X -6.8126437e-14 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 centre_Y 5.3885711e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 centre_Z -968.883 +- 0 -968.883 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 angles_X 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 angles_Y 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 angles_Z 0 +- 0 0 +- 100 Q0 %%%% Optical Object: CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 centre_X -4.4979915e-14 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 centre_Y 7.1823513e-13 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 centre_Z 988.187 +- 0 988.187 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 angles_X 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 angles_Y 0 +- 0 0 +- 100 Q0 FIX: -1 CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 angles_Z 0 +- 0 0 +- 100 Q0 %%%% Optical Object: CMS/slm_p2_2/MEp2_1_04 CAL: 105 CMS/slm_p2_2/MEp2_1_04 centre_X -2371.279 +- 100 -2371.279 +- 100 Q1 UNK: 106 CMS/slm_p2_2/MEp2_1_04 centre_Y -112.39127 +- 3138.1779 -93.35 +- 100 Q2 UNK: 107 CMS/slm_p2_2/MEp2_1_04 centre_Z 121.96971 +- 1315.8838 118.16258 +- 100 Q2 CAL: 108 CMS/slm_p2_2/MEp2_1_04 angles_X -7.672071e-05 +- 99.997934 0 +- 100 Q1 CAL: 109 CMS/slm_p2_2/MEp2_1_04 angles_Y -79.999992 +- 99.999991 -80 +- 100 Q1 CAL: 110 CMS/slm_p2_2/MEp2_1_04 angles_Z 2.1077554e-06 +- 99.999999 0 +- 100 Q1 %%%% Optical Object: CMS/slm_p2_2 FIX: -1 CMS/slm_p2_2 centre_X 289.7777 +- 0 289.7777 +- 1 Q0 FIX: -1 CMS/slm_p2_2 centre_Y -77.6457 +- 0 -77.6457 +- 1 Q0 FIX: -1 CMS/slm_p2_2 centre_Z 7998 +- 0 7998 +- 1 Q0 FIX: -1 CMS/slm_p2_2 angles_X -90 +- 0 -90 +- 1 Q0 FIX: -1 CMS/slm_p2_2 angles_Y 0 +- 0 0 +- 1 Q0 FIX: -1 CMS/slm_p2_2 angles_Z 75 +- 0 255 +- 1 Q0 In Global Coordinates: OPTICALOBJECT list size 21 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 CENTRE GLOBAL: (-1.5668483,-7.0091331,7.99844) CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 ROT_MATRIX GLOBAL: xx=-0.9659544 xy=0.00038155989 xz=-0.25871209 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 ROT_MATRIX GLOBAL: yx=0.25871237 yy=0.0014217865 yz=-0.96595336 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 ROT_MATRIX GLOBAL: zx=-7.3570829e-07 zy=-0.99999892 zz=-0.0014720952 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN CENTRE GLOBAL: (-1.589814,-7.002979,8.0403104) CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN ROT_MATRIX GLOBAL: xx=-0.9659276 xy=-9.7553327e-06 xz=-0.25881242 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN ROT_MATRIX GLOBAL: yx=0.25881242 yy=2.6972759e-06 yz=-0.9659276 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN ROT_MATRIX GLOBAL: zx=1.0121034e-05 zy=-1 zz=-8.0571949e-08 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT CENTRE GLOBAL: (-1.1317849,-7.0926774,7.958975) CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT ROT_MATRIX GLOBAL: xx=-0.9659276 xy=-0.25881241 xz=-1.7786055e-08 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT ROT_MATRIX GLOBAL: yx=0.25881241 yy=-0.9659276 yz=-6.6374249e-08 CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_transfer_5_OUT ROT_MATRIX GLOBAL: zx=-1.5618631e-12 zy=-6.8715971e-08 zz=1 CMS/slm_p2_2/MEp2_tp5 CENTRE GLOBAL: (-1.5891979,-7.0031447,8.04125) CMS/slm_p2_2/MEp2_tp5 ROT_MATRIX GLOBAL: xx=-0.9659276 xy=1.7786055e-08 xz=-0.25881241 CMS/slm_p2_2/MEp2_tp5 ROT_MATRIX GLOBAL: yx=0.25881241 yy=6.6374249e-08 yz=-0.9659276 CMS/slm_p2_2/MEp2_tp5 ROT_MATRIX GLOBAL: zx=-1.5619856e-12 zy=-1 zz=-6.8715971e-08 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN CENTRE GLOBAL: (-1.4704842,-6.7505439,8.0396357) CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN ROT_MATRIX GLOBAL: xx=0.96592579 xy=-2.4571727e-07 xz=0.25881917 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN ROT_MATRIX GLOBAL: yx=-0.25881917 yy=-1.3537539e-06 yz=0.96592579 CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN ROT_MATRIX GLOBAL: zx=1.130328e-07 zy=-1 zz=-1.3712221e-06 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN CENTRE GLOBAL: (-0.69756037,-3.8659534,8.0436316) CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN ROT_MATRIX GLOBAL: xx=0.96592579 xy=-2.5558327e-07 xz=0.25881917 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN ROT_MATRIX GLOBAL: yx=-0.25881917 yy=-1.3520142e-06 yz=0.96592579 CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN ROT_MATRIX GLOBAL: zx=1.0305272e-07 zy=-1 zz=-1.3720953e-06 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 CENTRE GLOBAL: (-0.61935945,-3.5371799,8.0897311) CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 ROT_MATRIX GLOBAL: xx=-0.98480773 xy=-2.53718e-07 xz=-0.1736483 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 ROT_MATRIX GLOBAL: yx=0.1736483 yy=-1.352343e-06 yz=-0.98480773 CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 ROT_MATRIX GLOBAL: zx=1.503138e-08 zy=-1 zz=1.3758556e-06 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 CENTRE GLOBAL: (-1.1942017,-6.7972697,8.0897357) CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 ROT_MATRIX GLOBAL: xx=-0.98480773 xy=-2.53718e-07 xz=-0.1736483 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 ROT_MATRIX GLOBAL: yx=0.1736483 yy=-1.352343e-06 yz=-0.98480773 CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 ROT_MATRIX GLOBAL: zx=1.503138e-08 zy=-1 zz=1.3758556e-06 CMS/slm_p2_2/MEp2_2_27 CENTRE GLOBAL: (-0.90678056,-5.1672248,8.0897334) CMS/slm_p2_2/MEp2_2_27 ROT_MATRIX GLOBAL: xx=-0.98480773 xy=-2.53718e-07 xz=-0.1736483 CMS/slm_p2_2/MEp2_2_27 ROT_MATRIX GLOBAL: yx=0.1736483 yy=-1.352343e-06 yz=-0.98480773 CMS/slm_p2_2/MEp2_2_27 ROT_MATRIX GLOBAL: zx=1.503138e-08 zy=-1 zz=1.3758556e-06 CMS/slm_p2_2/MEp2_1_14/DCOPS_outer_MEp2_1_14_IN CENTRE GLOBAL: (-0.56424002,-3.3641415,8.0444382) CMS/slm_p2_2/MEp2_1_14/DCOPS_outer_MEp2_1_14_IN ROT_MATRIX GLOBAL: xx=0.96592585 xy=-4.0685804e-08 xz=0.25881897 CMS/slm_p2_2/MEp2_1_14/DCOPS_outer_MEp2_1_14_IN ROT_MATRIX GLOBAL: yx=-0.25881897 yy=-3.6203515e-07 yz=0.96592585 CMS/slm_p2_2/MEp2_1_14/DCOPS_outer_MEp2_1_14_IN ROT_MATRIX GLOBAL: zx=5.4402094e-08 zy=-1 zz=-3.6022937e-07 CMS/slm_p2_2/MEp2_1_14/DCOPS_inner_MEp2_1_14_IN CENTRE GLOBAL: (-0.11659626,-1.5197136,8.0474375) CMS/slm_p2_2/MEp2_1_14/DCOPS_inner_MEp2_1_14_IN ROT_MATRIX GLOBAL: xx=-0.96592585 xy=-4.1448547e-08 xz=-0.25881897 CMS/slm_p2_2/MEp2_1_14/DCOPS_inner_MEp2_1_14_IN ROT_MATRIX GLOBAL: yx=0.25881897 yy=-3.6196893e-07 yz=-0.96592585 CMS/slm_p2_2/MEp2_1_14/DCOPS_inner_MEp2_1_14_IN ROT_MATRIX GLOBAL: zx=-5.3648203e-08 zy=-1 zz=3.6036282e-07 CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 CENTRE GLOBAL: (-0.12228411,-1.43434,8.0895374) CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 ROT_MATRIX GLOBAL: xx=-0.99619471 xy=-4.1448547e-08 xz=-0.087155663 CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 ROT_MATRIX GLOBAL: yx=0.087155663 yy=-3.6196893e-07 yz=-0.99619471 CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 ROT_MATRIX GLOBAL: zx=9.7431803e-09 zy=-1 zz=3.6420401e-07 CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 CENTRE GLOBAL: (-0.29285384,-3.3839628,8.0895381) CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 ROT_MATRIX GLOBAL: xx=-0.99619471 xy=-4.1448547e-08 xz=-0.087155663 CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 ROT_MATRIX GLOBAL: yx=0.087155663 yy=-3.6196893e-07 yz=-0.99619471 CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 ROT_MATRIX GLOBAL: zx=9.7431803e-09 zy=-1 zz=3.6420401e-07 CMS/slm_p2_2/MEp2_1_14 CENTRE GLOBAL: (-0.20672775,-2.3995361,8.0895378) CMS/slm_p2_2/MEp2_1_14 ROT_MATRIX GLOBAL: xx=-0.99619471 xy=-4.1448547e-08 xz=-0.087155663 CMS/slm_p2_2/MEp2_1_14 ROT_MATRIX GLOBAL: yx=0.087155663 yy=-3.6196893e-07 yz=-0.99619471 CMS/slm_p2_2/MEp2_1_14 ROT_MATRIX GLOBAL: zx=9.7431803e-09 zy=-1 zz=3.6420401e-07 CMS/slm_p2_2/MEp2_1_04/DCOPS_inner_MEp2_1_04_OUT CENTRE GLOBAL: (0.65946839,1.3742514,8.0682924) CMS/slm_p2_2/MEp2_1_04/DCOPS_inner_MEp2_1_04_OUT ROT_MATRIX GLOBAL: xx=-0.96592579 xy=-5.5637703e-07 xz=-0.25881918 CMS/slm_p2_2/MEp2_1_04/DCOPS_inner_MEp2_1_04_OUT ROT_MATRIX GLOBAL: yx=0.25881918 yy=-1.1780385e-06 yz=-0.96592579 CMS/slm_p2_2/MEp2_1_04/DCOPS_inner_MEp2_1_04_OUT ROT_MATRIX GLOBAL: zx=2.3251995e-07 zy=-1 zz=1.2818989e-06 CMS/slm_p2_2/MEp2_1_04/DCOPS_outer_MEp2_1_04_OUT CENTRE GLOBAL: (1.1940116,3.1953947,8.06529) CMS/slm_p2_2/MEp2_1_04/DCOPS_outer_MEp2_1_04_OUT ROT_MATRIX GLOBAL: xx=0.96592579 xy=-5.593711e-07 xz=0.25881918 CMS/slm_p2_2/MEp2_1_04/DCOPS_outer_MEp2_1_04_OUT ROT_MATRIX GLOBAL: yx=-0.25881918 yy=-1.176643e-06 yz=0.96592579 CMS/slm_p2_2/MEp2_1_04/DCOPS_outer_MEp2_1_04_OUT ROT_MATRIX GLOBAL: zx=-2.357732e-07 zy=-1 zz=-1.2813258e-06 CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 CENTRE GLOBAL: (0.61185579,1.3031597,8.1103925) CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 ROT_MATRIX GLOBAL: xx=0.90630773 xy=-5.5637703e-07 xz=0.42261839 CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 ROT_MATRIX GLOBAL: yx=-0.42261839 yy=-1.1780385e-06 yz=0.90630773 CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 ROT_MATRIX GLOBAL: zx=-6.3880476e-09 zy=-1 zz=-1.3028006e-06 CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 CENTRE GLOBAL: (1.4389496,3.0768674,8.11039) CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 ROT_MATRIX GLOBAL: xx=0.90630773 xy=-5.5637703e-07 xz=0.42261839 CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 ROT_MATRIX GLOBAL: yx=-0.42261839 yy=-1.1780385e-06 yz=0.90630773 CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 ROT_MATRIX GLOBAL: zx=-6.3880476e-09 zy=-1 zz=-1.3028006e-06 CMS/slm_p2_2/MEp2_1_04 CENTRE GLOBAL: (1.0213236,2.1812658,8.1103913) CMS/slm_p2_2/MEp2_1_04 ROT_MATRIX GLOBAL: xx=0.90630773 xy=-5.5637703e-07 xz=0.42261839 CMS/slm_p2_2/MEp2_1_04 ROT_MATRIX GLOBAL: yx=-0.42261839 yy=-1.1780385e-06 yz=0.90630773 CMS/slm_p2_2/MEp2_1_04 ROT_MATRIX GLOBAL: zx=-6.3880476e-09 zy=-1 zz=-1.3028006e-06 CMS/slm_p2_2 CENTRE GLOBAL: (0.2897777,-0.0776457,7.998) CMS/slm_p2_2 ROT_MATRIX GLOBAL: xx=-0.25881905 xy=5.9143945e-17 xz=0.96592583 CMS/slm_p2_2 ROT_MATRIX GLOBAL: yx=-0.96592583 yy=-1.5847572e-17 yz=-0.25881905 CMS/slm_p2_2 ROT_MATRIX GLOBAL: zx=0 zy=-1 zz=6.1230318e-17 CMS CENTRE GLOBAL: (0,0,0) Some Observations: %%%% Optical Object: CMS/slm_p2_2/MEp2_2_27 CAL: 37 CMS/slm_p2_2/MEp2_2_27 centre_X 5225.848 +- 100 5225.848 +- 100 Q1 UNK: 38 CMS/slm_p2_2/MEp2_2_27 centre_Y -91.733375 +- 652.2979 -93.35 +- 100 Q2 UNK: 39 CMS/slm_p2_2/MEp2_2_27 centre_Z 161.49349 +- 256.35984 157.20274 +- 100 Q2 CAL: 40 CMS/slm_p2_2/MEp2_2_27 angles_X 6.8986745e-05 +- 99.996742 0 +- 100 Q1 CAL: 41 CMS/slm_p2_2/MEp2_2_27 angles_Y 85.000007 +- 99.999972 85 +- 100 Q1 CAL: 42 CMS/slm_p2_2/MEp2_2_27 angles_Z -9.8815737e-06 +- 99.999717 0 +- 100 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_1_14 CAL: 71 CMS/slm_p2_2/MEp2_1_14 centre_X 2371.279 +- 100 2371.279 +- 100 Q1 UNK: 72 CMS/slm_p2_2/MEp2_1_14 centre_Y -91.537789 +- 1236.663 -93.35 +- 100 Q2 UNK: 73 CMS/slm_p2_2/MEp2_1_14 centre_Z 121.36203 +- 499.25084 118.16258 +- 100 Q2 CAL: 74 CMS/slm_p2_2/MEp2_1_14 angles_X 1.77014e-05 +- 99.994517 0 +- 100 Q1 CAL: 75 CMS/slm_p2_2/MEp2_1_14 angles_Y 79.999995 +- 99.999998 80 +- 100 Q1 CAL: 76 CMS/slm_p2_2/MEp2_1_14 angles_Z -3.2147925e-06 +- 99.999998 0 +- 100 Q1 %%%% Optical Object: CMS/slm_p2_2/MEp2_1_04 CAL: 105 CMS/slm_p2_2/MEp2_1_04 centre_X -2371.279 +- 100 -2371.279 +- 100 Q1 UNK: 106 CMS/slm_p2_2/MEp2_1_04 centre_Y -112.39127 +- 3138.1779 -93.35 +- 100 Q2 UNK: 107 CMS/slm_p2_2/MEp2_1_04 centre_Z 121.96971 +- 1315.8838 118.16258 +- 100 Q2 CAL: 108 CMS/slm_p2_2/MEp2_1_04 angles_X -7.672071e-05 +- 99.997934 0 +- 100 Q1 CAL: 109 CMS/slm_p2_2/MEp2_1_04 angles_Y -79.999992 +- 99.999991 -80 +- 100 Q1 CAL: 110 CMS/slm_p2_2/MEp2_1_04 angles_Z 2.1077554e-06 +- 99.999999 0 +- 100 Q1 - For all 3 chambers, the chamber angles are constrained to CAL and come out to be basically as their nominal values (0's,+-80,85). Good. - Y_SLM=Z_CMS values are Ch. 2_27 1_14 1_04 Y_SLM -91.7 -91.5 -112.4 [mm] So 1_04 is much lower in Y_SLM by 21mm, but that presumably reflects the data: COPS id14 U 9.896 100.881 D 10.544 263.711 L 6.032 6.180 R 6.552 1.942 // file id: 39 COPS id15 U 10.233 1.739 D 8.841 1.466 L 13.326 1.658 R 13.691 1.834 whereas for the other two chambers L and R measure 17-19 mm. That's a shift of about 18-(13+6)/2 = 10mm. Let's predict where the laser SHOULD hit if the chamber were at nominal height: From near short inner DCOPS to far short inner DCOPS is about 3m (pos. along SLM about +- 1.5m). With a fitted dip angle of 1.4 mrad (see next item), that would be a displacement of 3*1.4 = 4mm down or about 14mm. Data are 13mm. Not too bad. Note that these DCOPS on near and far have the same shimming, so should not be affected by shimming. From far inner to far outer DCOPS is 2mm, so expect 2*1.4 = 3mm plus shim difference of another 3mm, so we expect to go from 14 to 14-6=8mm. Data are 6mm. Again, not too far off! So, why does the fit put the far short chamber so much lower in Y_SLM? 20mm is a lot. - Note that the laser dip angle was not constrained here: UNK: 6 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_X -0.084340938 +- 1028.8778 0 +- 10 Q2 UNK: 7 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Y -0.0059342079 +- 708.73878 0 +- 10 Q2 CAL: 8 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 angles_Z 4.2152483e-05 +- 9.99951 0 +- 10 Q1 Note that the dip angle should correspond to a rotation around X_TP. The corresponding angles_X = -0.084 deg = 1.4+-1 mrad is not too unreasonable. If we assume that the laser was tuned so that it hit the 8th DCOPS at 6mm, i.e. 8mm below center, that would produce an angle of (8mm below + 4mm shim down)/11m = 1.1 mrad If they picked a usual right-handed coordinate system, then a negative angle would mean a rotation around X with the y-axis moving away from the z-axis. This is indeed a "dip down" motion. Ok, seems consistent. - But the laser position was constrained: CAL: 4 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_X -23.137983 +- 18.579452 -22.5 +- 20 Q1 CAL: 5 CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 centre_Y 42.810046 +- 18.584221 43.25 +- 20 Q1 Note that it actually gets moved appreciably in its TP ref. system: by 0.6mm in X and by 0.4mm in Y - Now, in SLM_Z or CMS_rphi, both short chambers come out very similar UNK: 73 CMS/slm_p2_2/MEp2_1_14 centre_Z 121.36203 +- 499.25084 118.16258 +- 100 Q2 UNK:107 CMS/slm_p2_2/MEp2_1_04 centre_Z 121.96971 +- 1315.8838 118.16258 +- 100 Q2 with a shift of about 3mm relative to nominal. This is the direction perp. to the SLM but in the station plane. Note the errors for these two and for MEP2_2_27, though: Ch. 2_27 1_14 1_04 Sigma_z 256 500 1315 [micron] The farther away from the laser, the worse the error. Makes sense ! - I can't get the DCOPS Y_SLM position from the output since the DCOPS are defined in the chamber system, BUT I can get it for all objects using global coordinates: CMS/slm_p2_2/MEp2_tp5/laser_MEp2_5 CENTRE GLOBAL: (-1.5668483,-7.0091331,7.99844) CMS/slm_p2_2/MEp2_tp5/DCOPS_MEp2_reference_5_IN CENTRE GLOBAL: (-1.589814,-7.002979,8.0403104) CMS/slm_p2_2/MEp2_tp5 CENTRE GLOBAL: (-1.5891979,-7.0031447,8.04125) CMS/slm_p2_2/MEp2_2_27/DCOPS_outer_MEp2_2_27_IN CENTRE GLOBAL: (-1.4704842,-6.7505439,8.0396357) CMS/slm_p2_2/MEp2_2_27/DCOPS_inner_MEp2_2_27_IN CENTRE GLOBAL: (-0.69756037,-3.8659534,8.0436316) CMS/slm_p2_2/MEp2_2_27/apin_inner_MEp2_2_27 CENTRE GLOBAL: (-0.61935945,-3.5371799,8.0897311) CMS/slm_p2_2/MEp2_2_27/apin_outer_MEp2_2_27 CENTRE GLOBAL: (-1.1942017,-6.7972697,8.0897357) CMS/slm_p2_2/MEp2_2_27 CENTRE GLOBAL: (-0.90678056,-5.1672248,8.0897334) CMS/slm_p2_2/MEp2_1_14/DCOPS_outer_MEp2_1_14_IN CENTRE GLOBAL: (-0.56424002,-3.3641415,8.0444382) CMS/slm_p2_2/MEp2_1_14/DCOPS_inner_MEp2_1_14_IN CENTRE GLOBAL: (-0.11659626,-1.5197136,8.0474375) CMS/slm_p2_2/MEp2_1_14/apin_inner_MEp2_1_14 CENTRE GLOBAL: (-0.12228411,-1.43434,8.0895374) CMS/slm_p2_2/MEp2_1_14/apin_outer_MEp2_1_14 CENTRE GLOBAL: (-0.29285384,-3.3839628,8.0895381) CMS/slm_p2_2/MEp2_1_14 CENTRE GLOBAL: (-0.20672775,-2.3995361,8.0895378) CMS/slm_p2_2/MEp2_1_04/DCOPS_inner_MEp2_1_04_OUT CENTRE GLOBAL: (0.65946839,1.3742514,8.0682924) CMS/slm_p2_2/MEp2_1_04/DCOPS_outer_MEp2_1_04_OUT CENTRE GLOBAL: (1.1940116,3.1953947,8.06529) CMS/slm_p2_2/MEp2_1_04/apin_inner_MEp2_1_04 CENTRE GLOBAL: (0.61185579,1.3031597,8.1103925) CMS/slm_p2_2/MEp2_1_04/apin_outer_MEp2_1_04 CENTRE GLOBAL: (1.4389496,3.0768674,8.11039) CMS/slm_p2_2/MEp2_1_04 CENTRE GLOBAL: (1.0213236,2.1812658,8.1103913) The reference z-position should be the z of the TP; i.e. 8.04125 m. I think what we get here are DCOPS dowel pin positions, which are going to higher z's. Ref : to lower z by 0.9mm long ch. outer: to lower z by 1.2mm inner: to higher z by 2.4mm near short out: to higher z by 3.6mm inn: to higher z by 6.2mm far short out: to higher z by 27mm !! inn: to higher z by 24mm !! The values for near side seem almost consistent with shimming (+0.5, 0, -4, -5,-8), but the values for the chamber on the far side seem again very high! Is the shimming for that chamber implemented in the wrong direction in the sdf ? Hmmm... It would actually help to see these global values for simulation so I get an idea what they "should" be; especially chambers and DCOPS relative to TP z-position. - Note that PG actually places the weird chamber MEp2_1_04 at lower CMS_z by 25cm relative to MEp2_2_27!! That's bizarre. See GB's Table 5! I verify GB's numbers against the PG report on YE+2. They look right. (There is an offset of 8843 mm on Z as the PG group gives Z relative to the center of the YE+2 disk itself and GB gives it relative to global CMS (=IP). CMS = PG + 8843. Here PG values are negative as the face of the YE+2 disk is towards neg.Z relative to the center of the YE+2 disk.) But wait... MEp2_2_27 is a back chamber, so it is at higher CMS Z! 25 cm could be the chamber thickness. - So, what's next ? Try to do the same for the laser from pt.2 plus 3 chambers ? I also should collate this comparison against PG from a converging B=0T fit into an Excel file a la GB. 2/25/2007 --------- - Set any of the coordinates of the center of the SLM itself to unk -> Cocoa crashes with the usual matrix inversion error. - Ok, let's try fitting Las2 for 7 DCOPS; this time with all chamber angles constrained to CAL and ref. DCOPS data values at original value. -> crashes - So, let's try to be systematic: laser 2 plus only near long chamber (3 DCOPS only) Fits, but: Too many iterations 49 and fit DOES NOT CONVERGE FIT STATUS COCOA_FitCannotImprovem FQ=17771139 -1382.75 adding chi2meas 1.912e+06 id1: 0 (mm)R: 8.686 S: 14.329 Diff= -5.643 -1430.79 adding chi2meas 3.95915e+06 id1: 1 (mm)R: 8.334 S: 14.329 Diff= -5.995 1458.66 adding chi2meas 6.08684e+06 id1: 2 (mm)R: 18.34 S: 14.829 Diff= 3.511 1384.27 adding chi2meas 8.00304e+06 id1: 3 (mm)R: 18.19 S: 14.829 Diff= 3.361 2774.94 adding chi2meas 1.57034e+07 id2: 0 (mm)R: 20.109 S: 14.3288 Diff= 5.78021 2345.83 adding chi2meas 2.12063e+07 id2: 1 (mm)R: 20.313 S: 14.3288 Diff= 5.98421 317.56 adding chi2meas 2.13071e+07 id2: 2 (mm)R: 18.134 S: 14.329 Diff= 3.805 2633.79 adding chi2meas 2.82439e+07 id2: 3 (mm)R: 19.162 S: 14.329 Diff= 4.833 11744.7 adding chi2meas 1.66181e+08 id3: 0 (mm)R: 18.674 S: 0 Diff= 18.674 6887.21 adding chi2meas 2.13614e+08 id3: 1 (mm)R: 18.196 S: 0 Diff= 18.196 6826.46 adding chi2meas 2.60215e+08 id3: 2 (mm)R: 18.213 S: 0 Diff= 18.213 9937.11 adding chi2meas 3.58961e+08 id3: 3 (mm)R: 18.96 S: 0 Diff= 18.96 Note again the weird S=0 values. %%%% Optical Object: CMS/slm_p2_2/MEp2_2_08 CAL: 37 CMS/slm_p2_2/MEp2_2_08 centre_X -5225.848 +- 100 -5225.848 +- 100 Q1 UNK: 38 CMS/slm_p2_2/MEp2_2_08 centre_Y -87.936391 +- 1225.4716 -93.35 +- 100 Q2 UNK: 39 CMS/slm_p2_2/MEp2_2_08 centre_Z 3356.2644 +- 136313.68 157.20274 +- 100 Q2 CAL: 40 CMS/slm_p2_2/MEp2_2_08 angles_X 0.0007552241 +- 99.999482 0 +- 100 Q1 CAL: 41 CMS/slm_p2_2/MEp2_2_08 angles_Y -84.999994 +- 99.999997 -85 +- 100 Q1 CAL: 42 CMS/slm_p2_2/MEp2_2_08 angles_Z 0.00063171218 +- 99.999944 0 +- 100 Q1 Note center_z shifts by 3.2m !! - Let's try the corrected U,D values, i.e. 2048*14 - measured data value (see above) Same answer: Too many iterations 49 and fit DOES NOT CONVERGE FIT STATUS COCOA_FitCannotImprove Fit iteration 49 ... 49 Fit quality is = 17980572 %%%% Optical Object: CMS/slm_p2_2/MEp2_2_08 CAL: 37 CMS/slm_p2_2/MEp2_2_08 centre_X -5225.848 +- 100 -5225.848 +- 100 Q1 UNK: 38 CMS/slm_p2_2/MEp2_2_08 centre_Y -89.921597 +- 1291.7304 -93.35 +- 100 Q2 UNK: 39 CMS/slm_p2_2/MEp2_2_08 centre_Z 3674.3222 +- 132394.04 157.20274 +- 100 Q2 CAL: 40 CMS/slm_p2_2/MEp2_2_08 angles_X 0.00050524342 +- 99.999509 0 +- 100 Q1 CAL: 41 CMS/slm_p2_2/MEp2_2_08 angles_Y -84.999995 +- 99.999997 -85 +- 100 Q1 CAL: 42 CMS/slm_p2_2/MEp2_2_08 angles_Z 0.00037758573 +- 99.999981 0 +- 100 Q1 Note center_z shifts by 3.5m !! 2/27/2007 --------- Login to cmsuaf.fnal.gov with cryptocard doesn't work. Call help desk at FNAL. They suggest to login directly to a worker node: ssh cmswn051.fnal.gov. This works with the cryptocard. 3/2/2007 -------- I get some interesting feedback from Celso Martinez (IFCA), who is doing COCOA work for the link. Here is the discussion: --- Hi Celso, well, that is already quite useful info for me! I have to admit that I probably haven't paid enough attention to these numbers because from the COCOA output it wasn't quite clear to me what they mean. So, you are saying that in the line below the meaning is 44.0379 adding chi2meas 1939.33 id14: 0 (mm)R: 9.896 S: 14.3386 Diff=-4.44258 Diff n | Diff |2 ---- Sum | ---- | sigma i=1 | sigma| for the n-th meas. (here n=1) Is that correct ? Now, the sigmas in the data file come from the uncertainty on the mean in the beam profile fits of the DCOPS and if it is a good profile, these can be fractions of a single pixel, so small values like 2 micron etc. are actually not unreasonable. So, then the problem is that the initial residual "Diff" is way too big. Now the measurement value (R) is what it is. That implies that the problem seems to be the S: values! These currently reflect the nominal positions where the laser hits the DCOPS. However, those are obviously inappropriate because in the nominal configuration the laser goes through the exact center of the DCOPS, whereas in the real B=0T data it is dipped down, which causes it to hit many mm off-center. How can I correct that ? How do I set the S: values to something more appropriate for the real world configuration so that the initial chi**2 is not so ridiculously large ? Also, notice that in some instances the S: is not the center value of the DCOPS (14.32mm), but 0: 44.0379 adding chi2meas 1939.33 id14: 0 (mm)R: 9.896 S: 14.3386 Diff= -4.44258 -14.3892 adding chi2meas 2146.38 id14: 1 (mm)R: 10.544 S: 14.3386 Diff= -3.79458 -533.495 adding chi2meas 286763 id14: 2 (mm)R: 6.032 S: 9.329 Diff= -3.297 -1429.97 adding chi2meas 2.33158e+06 id14: 3 (mm)R: 6.552 S: 9.329 Diff= -2.777 5884.42 adding chi2meas 3.69579e+07 id15: 0 (mm)R: 10.233 S: 0 Diff= 10.233 6030.7 adding chi2meas 7.33272e+07 id15: 1 (mm)R: 8.841 S: 0 Diff= 8.841 8037.39 adding chi2meas 1.37927e+08 id15: 2 (mm)R: 13.326 S: 0 Diff= 13.326 7465.1 adding chi2meas 1.93655e+08 id15: 3 (mm)R: 13.691 S: 0 Diff= 13.691 4976.7 adding chi2meas 2.18422e+08 id16: 0 (mm)R: 9.61 S: 0 Diff= 9.61 3785.17 adding chi2meas 2.3275e+08 id16: 1 (mm)R: 10.36 S: 0 Diff= 10.36 503.709 adding chi2meas 2.33003e+08 id16: 2 (mm)R: 18.81 S: 0 Diff= 18.81 9780.06 adding chi2meas 3.28653e+08 id16: 3 (mm)R: 19.032 S: 0 Diff= 19.032 That doesn't make sense. As I say in the log file, I had also noticed that which DCOP sensor shows such 0's depends on how many DCOPS I am trying to fit. That seems pretty strange. Any ideas on that? Anyway, this is really good feedback. Let's keep this kind of a discussion going ! Thanks! Marcus > Hi Marcus, > > I had a quick look just to the chi2 of each measurement and indeed they are huge! (look for > chi2meas in your log). In the left you see the number of sigmas and in the right you see the > residual (diff). normally the residual should be within some sigmas of the measurements, but you > have some with more than 1000 sigmas! they way to arrange this is to change the precision of each > measurement. > > For example in the first line where you have a Diff = -4.44258 the number to the left is -44.0370 > so I presume that your precision is about 100 microns for that measurement; but further down you > have a Diff=10.233 with a chi of 5884 that meas a precision of about 2 microns... > > All this numbers as I tell you should be kept low...at least that�s my experience... Greetings, > > Celso --- So, there is a lot of good info in the screen output that just whizzes by.(Why is this not included in the report.out?) So, let's write it to a file: $ cocoa ME2_SLM25_WITH_SHIMS_AND_ROTATIONS_RESTRICTED_LAS5_3ch_7DCOPS_CONVERGING_Official_names.txt > screen.out Here is what I get after 3 iterations: -0.741478 adding chi2meas 0.54979 id14: 0 (mm)R: 9.896 S: 9.9708 Diff= -0.0748011 1.93822 adding chi2meas 4.3065 id14: 1 (mm)R: 10.544 S: 10.0329 Diff= 0.511131 -0.197805 adding chi2meas 4.34563 id14: 2 (mm)R: 6.032 S: 6.03322 Diff= -0.00122243 0.0628385 adding chi2meas 4.34958 id14: 3 (mm)R: 6.552 S: 6.55188 Diff= 0.000122032 5884.42 adding chi2meas 3.46264e+07 id15: 0 (mm)R: 10.233 S: 0 Diff= 10.233 6030.7 adding chi2meas 7.09957e+07 id15: 1 (mm)R: 8.841 S: 0 Diff= 8.841 8037.39 adding chi2meas 1.35595e+08 id15: 2 (mm)R: 13.326 S: 0 Diff= 13.326 7465.1 adding chi2meas 1.91323e+08 id15: 3 (mm)R: 13.691 S: 0 Diff= 13.691 4976.7 adding chi2meas 2.16091e+08 id16: 0 (mm)R: 9.61 S: 0 Diff= 9.61 3785.17 adding chi2meas 2.30418e+08 id16: 1 (mm)R: 10.36 S: 0 Diff= 10.36 503.709 adding chi2meas 2.30672e+08 id16: 2 (mm)R: 18.81 S: 0 Diff= 18.81 9780.06 adding chi2meas 3.26321e+08 id16: 3 (mm)R: 19.032 S: 0 Diff= 19.032 0.0373458 adding chi2meas 3.26321e+08 id17: 0 (mm)R: 11.447 S: 11.4469 Diff= 6.74465e-05 -0.059678 adding chi2meas 3.26321e+08 id17: 1 (mm)R: 11.272 S: 11.2722 Diff= -0.000172231 0.211985 adding chi2meas 3.26321e+08 id17: 2 (mm)R: 17.142 S: 17.1415 Diff= 0.000548405 -0.0602411 adding chi2meas 3.26321e+08 id17: 3 (mm)R: 17.134 S: 17.1342 Diff= -0.000243555 0.0103033 adding chi2meas 3.26321e+08 id18: 0 (mm)R: 10.341 S: 10.341 Diff= 1.7969e-05 -0.0105247 adding chi2meas 3.26321e+08 id18: 1 (mm)R: 10.291 S: 10.291 Diff= -1.87339e-05 -0.0289562 adding chi2meas 3.26321e+08 id18: 2 (mm)R: 17.038 S: 17.0381 Diff= -5.84047e-05 -0.0128892 adding chi2meas 3.26321e+08 id18: 3 (mm)R: 17.337 S: 17.337 Diff= -2.53272e-05 -0.200938 adding chi2meas 3.26321e+08 id19: 0 (mm)R: 10.428 S: 10.4288 Diff= -0.000785468 0.20903 adding chi2meas 3.26321e+08 id19: 1 (mm)R: 10.864 S: 10.8632 Diff= 0.000849915 -0.0261488 adding chi2meas 3.26321e+08 id19: 2 (mm)R: 16.628 S: 16.6281 Diff= -7.9414e-05 0.032077 adding chi2meas 3.26321e+08 id19: 3 (mm)R: 16.942 S: 16.9419 Diff= 0.000135782 -2.47755 adding chi2meas 3.26321e+08 id20: 0 (mm)R: 9.041 S: 9.04748 Diff= -0.00647631 -2.97582 adding chi2meas 3.26321e+08 id20: 1 (mm)R: 10.275 S: 10.2885 Diff= -0.013537 2.01075 adding chi2meas 3.26321e+08 id20: 2 (mm)R: 17.769 S: 17.76 Diff= 0.0090162 2.48365 adding chi2meas 3.26321e+08 id20: 3 (mm)R: 18.351 S: 18.3414 Diff= 0.00957945 TIME:MAT_MEAS_FILLED: 1340000 0.06 TIME:MAT_CAL_FILLED : 1350000 0.01 TIME:BEFORE_INVERSE : 1360000 0.01 TIME:AFTER_INVERSE : 1370000 0.01 TIME:MAT_MULTIPLIED : 1370000 0 Fit iteration 3 ... 3 Fit quality is = 0.0348128 3 Fit quality predicted improvement in distance to minimum is = 0.0348128 3 Fit quality derivative at 0 = -0.0696255 TIME:QUALITY_CHECKED: 1420000 0.05 SRPARPOS Optical Object Parameter Fit.Value Orig.Value FITTED VALUES CORR:CU (0) (3) 2.30062e-05 .. FIT STATUS COCOA_FitOK 0x9d92700$$$$$$$$$$$$$$$ moreData Sets 0 ............ program ended OK TIME:PROGRAM ENDED : 1490000 0.07 TIME:TOTAL PROGRAM : 1490000 0.88 - Observations: 1) Residuals in terms of sigma are really small for most measurement, that's great ! What kills the Chi**2 is the two DCOPS id15 and id16 that have S:0 values throughout the whole fitting process !!! So that needs fixing. How? Gyongyi tells me that she had the same problem a while back even with nominal data! 2) Even though the initial S: values are far off, after the first iteration R: and S: values are already very close. Also, notice that the files fittedEntries.out longFittedEntries.out indeed contain chamber positions in text form. I could try to read those in with PAW to make plots. Here is a "COCOA To Do" list as I see it at this moment: -------------------------------------------------------- 1) (Most important): Solve the S: 0 problem 2) Write all COCOA output into root file for analysis 3) Get input service going for COCOA so that it can read Samir's root files 4) Put chi**2 information in the report.out so that all important information is in a single file 5) Provide more error diagnostic in the case that COCOA crashes with matrix inversion errors 6) Give errors on all fitted quantities in the CMS coordinate system 7) Update COCOA documentation