CMSSW Reference Manual

00001 
00009 #include "Alignment/LaserAlignment/interface/LaserAlignmentPosTEC.h"
00010 #include "FWCore/MessageLogger/interface/MessageLogger.h"
00011 
00012 LaserAlignmentPosTEC::LaserAlignmentPosTEC()
00013 {
00014 }
00015 
00016 LaserAlignmentPosTEC::~LaserAlignmentPosTEC()
00017 {
00018 }
00019 
00020 void LaserAlignmentPosTEC::alignment(edm::ParameterSet const & theConf, 
00021                                      AlignableTracker * theAlignableTracker,
00022                                      int theNumberOfIterations, 
00023                                      int theAlignmentIteration,
00024                                      std::vector<double>& theLaserPhi, 
00025                                      std::vector<double>& theLaserPhiError)
00026 {
00027   edm::LogInfo("LaserAlignmentPosTEC") << " ***************************************************** "
00028                                  << "\n *                Alignment of TEC+                  * "
00029                                  << "\n ***************************************************** ";
00030   
00031   int theMaxIteration = theNumberOfIterations - 1;
00032 
00033   // create an alignment Object for the positive TEC
00034   theLaserAlignmentTrackerPosTEC = new LaserAlignmentAlgorithmPosTEC(theConf, theAlignmentIteration);
00035       
00036   // do the iterations for the local fits
00037   for (int theIteration = 0; theIteration < theNumberOfIterations; theIteration++)
00038     {
00039       // fill fitted Phi position of the Beams and the Errors on Phi into 
00040       // a map. Afterwards add the beams to the LaserAlignmentTracker ...
00041           
00042       LogDebug("LaserAlignmentPosTEC") << "  AC1CMS: Total number of Iterations = " << theMaxIteration
00043                                  << "\n  AC1CMS: Current Iteration = " << theIteration
00044                                  << "\n  AC1CMS: Current Alignment Iteration = "  << theAlignmentIteration;
00045       
00046       // string is the name (i.e. Beam 0 in TEC+) and std::vector contains
00047       // 9 values for phi and the error on phi (on each Disc!) for the beam
00048       std::map<std::string, std::vector<double> > theBeamPosition;
00049       
00050       // fill now the map 
00051       // therefore loop over theLaserPhi and fill
00052       // Phi + Error in the map for each beam
00053       
00054       // Beams in TEC+
00055       std::vector<double> Beam0R4PosTEC;
00056       std::vector<double> Beam1R4PosTEC;
00057       std::vector<double> Beam2R4PosTEC;
00058       std::vector<double> Beam3R4PosTEC;
00059       std::vector<double> Beam4R4PosTEC;
00060       std::vector<double> Beam5R4PosTEC;
00061       std::vector<double> Beam6R4PosTEC;
00062       std::vector<double> Beam7R4PosTEC;
00063       
00064       std::vector<double> Beam0R6PosTEC;
00065       std::vector<double> Beam1R6PosTEC;
00066       std::vector<double> Beam2R6PosTEC;
00067       std::vector<double> Beam3R6PosTEC;
00068       std::vector<double> Beam4R6PosTEC;
00069       std::vector<double> Beam5R6PosTEC;
00070       std::vector<double> Beam6R6PosTEC;
00071       std::vector<double> Beam7R6PosTEC;
00072       
00073       // Beams in TEC+
00074       Beam0R4PosTEC.push_back(theLaserPhi.at(0 + theIteration * 434));
00075       Beam0R4PosTEC.push_back(theLaserPhiError.at(0 + theIteration * 434));
00076       Beam0R4PosTEC.push_back(theLaserPhi.at(1 + theIteration * 434));
00077       Beam0R4PosTEC.push_back(theLaserPhiError.at(1 + theIteration * 434));
00078       Beam0R4PosTEC.push_back(theLaserPhi.at(2 + theIteration * 434));
00079       Beam0R4PosTEC.push_back(theLaserPhiError.at(2 + theIteration * 434));
00080       Beam0R4PosTEC.push_back(theLaserPhi.at(3 + theIteration * 434));
00081       Beam0R4PosTEC.push_back(theLaserPhiError.at(3 + theIteration * 434));
00082       Beam0R4PosTEC.push_back(theLaserPhi.at(4 + theIteration * 434));
00083       Beam0R4PosTEC.push_back(theLaserPhiError.at(4 + theIteration * 434));
00084       Beam0R4PosTEC.push_back(theLaserPhi.at(5 + theIteration * 434));
00085       Beam0R4PosTEC.push_back(theLaserPhiError.at(5 + theIteration * 434));
00086       Beam0R4PosTEC.push_back(theLaserPhi.at(6 + theIteration * 434));
00087       Beam0R4PosTEC.push_back(theLaserPhiError.at(6 + theIteration * 434));
00088       Beam0R4PosTEC.push_back(theLaserPhi.at(7 + theIteration * 434));
00089       Beam0R4PosTEC.push_back(theLaserPhiError.at(7 + theIteration * 434));
00090       Beam0R4PosTEC.push_back(theLaserPhi.at(8 + theIteration * 434));
00091       Beam0R4PosTEC.push_back(theLaserPhiError.at(8 + theIteration * 434));
00092       
00093       Beam1R4PosTEC.push_back(theLaserPhi.at(9 + theIteration * 434));
00094       Beam1R4PosTEC.push_back(theLaserPhiError.at(9 + theIteration * 434));
00095       Beam1R4PosTEC.push_back(theLaserPhi.at(10 + theIteration * 434));
00096       Beam1R4PosTEC.push_back(theLaserPhiError.at(10 + theIteration * 434));
00097       Beam1R4PosTEC.push_back(theLaserPhi.at(11 + theIteration * 434));
00098       Beam1R4PosTEC.push_back(theLaserPhiError.at(11 + theIteration * 434));
00099       Beam1R4PosTEC.push_back(theLaserPhi.at(12 + theIteration * 434));
00100       Beam1R4PosTEC.push_back(theLaserPhiError.at(12 + theIteration * 434));
00101       Beam1R4PosTEC.push_back(theLaserPhi.at(13 + theIteration * 434));
00102       Beam1R4PosTEC.push_back(theLaserPhiError.at(13 + theIteration * 434));
00103       Beam1R4PosTEC.push_back(theLaserPhi.at(14 + theIteration * 434));
00104       Beam1R4PosTEC.push_back(theLaserPhiError.at(14 + theIteration * 434));
00105       Beam1R4PosTEC.push_back(theLaserPhi.at(15 + theIteration * 434));
00106       Beam1R4PosTEC.push_back(theLaserPhiError.at(15 + theIteration * 434));
00107       Beam1R4PosTEC.push_back(theLaserPhi.at(16 + theIteration * 434));
00108       Beam1R4PosTEC.push_back(theLaserPhiError.at(16 + theIteration * 434));
00109       Beam1R4PosTEC.push_back(theLaserPhi.at(17 + theIteration * 434));
00110       Beam1R4PosTEC.push_back(theLaserPhiError.at(17 + theIteration * 434));
00111           
00112           
00113       Beam2R4PosTEC.push_back(theLaserPhi.at(23 + theIteration * 434));
00114       Beam2R4PosTEC.push_back(theLaserPhiError.at(23 + theIteration * 434));
00115       Beam2R4PosTEC.push_back(theLaserPhi.at(24 + theIteration * 434));
00116       Beam2R4PosTEC.push_back(theLaserPhiError.at(24 + theIteration * 434));
00117       Beam2R4PosTEC.push_back(theLaserPhi.at(25 + theIteration * 434));
00118       Beam2R4PosTEC.push_back(theLaserPhiError.at(25 + theIteration * 434));
00119       Beam2R4PosTEC.push_back(theLaserPhi.at(26 + theIteration * 434));
00120       Beam2R4PosTEC.push_back(theLaserPhiError.at(26 + theIteration * 434));
00121       Beam2R4PosTEC.push_back(theLaserPhi.at(27 + theIteration * 434));
00122       Beam2R4PosTEC.push_back(theLaserPhiError.at(27 + theIteration * 434));
00123       Beam2R4PosTEC.push_back(theLaserPhi.at(28 + theIteration * 434));
00124       Beam2R4PosTEC.push_back(theLaserPhiError.at(28 + theIteration * 434));
00125       Beam2R4PosTEC.push_back(theLaserPhi.at(29 + theIteration * 434));
00126       Beam2R4PosTEC.push_back(theLaserPhiError.at(29 + theIteration * 434));
00127       Beam2R4PosTEC.push_back(theLaserPhi.at(30 + theIteration * 434));
00128       Beam2R4PosTEC.push_back(theLaserPhiError.at(30 + theIteration * 434));
00129       Beam2R4PosTEC.push_back(theLaserPhi.at(31 + theIteration * 434));
00130       Beam2R4PosTEC.push_back(theLaserPhiError.at(31 + theIteration * 434));
00131           
00132           
00133       Beam3R4PosTEC.push_back(theLaserPhi.at(37 + theIteration * 434));
00134       Beam3R4PosTEC.push_back(theLaserPhiError.at(37 + theIteration * 434));
00135       Beam3R4PosTEC.push_back(theLaserPhi.at(38 + theIteration * 434));
00136       Beam3R4PosTEC.push_back(theLaserPhiError.at(38 + theIteration * 434));
00137       Beam3R4PosTEC.push_back(theLaserPhi.at(39 + theIteration * 434));
00138       Beam3R4PosTEC.push_back(theLaserPhiError.at(39 + theIteration * 434));
00139       Beam3R4PosTEC.push_back(theLaserPhi.at(40 + theIteration * 434));
00140       Beam3R4PosTEC.push_back(theLaserPhiError.at(40 + theIteration * 434));
00141       Beam3R4PosTEC.push_back(theLaserPhi.at(41 + theIteration * 434));
00142       Beam3R4PosTEC.push_back(theLaserPhiError.at(41 + theIteration * 434));
00143       Beam3R4PosTEC.push_back(theLaserPhi.at(42 + theIteration * 434));
00144       Beam3R4PosTEC.push_back(theLaserPhiError.at(42 + theIteration * 434));
00145       Beam3R4PosTEC.push_back(theLaserPhi.at(43 + theIteration * 434));
00146       Beam3R4PosTEC.push_back(theLaserPhiError.at(43 + theIteration * 434));
00147       Beam3R4PosTEC.push_back(theLaserPhi.at(44 + theIteration * 434));
00148       Beam3R4PosTEC.push_back(theLaserPhiError.at(44 + theIteration * 434));
00149       Beam3R4PosTEC.push_back(theLaserPhi.at(45 + theIteration * 434));
00150       Beam3R4PosTEC.push_back(theLaserPhiError.at(45 + theIteration * 434));
00151           
00152       Beam4R4PosTEC.push_back(theLaserPhi.at(46 + theIteration * 434));
00153       Beam4R4PosTEC.push_back(theLaserPhiError.at(46 + theIteration * 434));
00154       Beam4R4PosTEC.push_back(theLaserPhi.at(47 + theIteration * 434));
00155       Beam4R4PosTEC.push_back(theLaserPhiError.at(47 + theIteration * 434));
00156       Beam4R4PosTEC.push_back(theLaserPhi.at(48 + theIteration * 434));
00157       Beam4R4PosTEC.push_back(theLaserPhiError.at(48 + theIteration * 434));
00158       Beam4R4PosTEC.push_back(theLaserPhi.at(49 + theIteration * 434));
00159       Beam4R4PosTEC.push_back(theLaserPhiError.at(49 + theIteration * 434));
00160       Beam4R4PosTEC.push_back(theLaserPhi.at(50 + theIteration * 434));
00161       Beam4R4PosTEC.push_back(theLaserPhiError.at(50 + theIteration * 434));
00162       Beam4R4PosTEC.push_back(theLaserPhi.at(51 + theIteration * 434));
00163       Beam4R4PosTEC.push_back(theLaserPhiError.at(51 + theIteration * 434));
00164       Beam4R4PosTEC.push_back(theLaserPhi.at(52 + theIteration * 434));
00165       Beam4R4PosTEC.push_back(theLaserPhiError.at(52 + theIteration * 434));
00166       Beam4R4PosTEC.push_back(theLaserPhi.at(53 + theIteration * 434));
00167       Beam4R4PosTEC.push_back(theLaserPhiError.at(53 + theIteration * 434));
00168       Beam4R4PosTEC.push_back(theLaserPhi.at(54 + theIteration * 434));
00169       Beam4R4PosTEC.push_back(theLaserPhiError.at(54 + theIteration * 434));
00170           
00171           
00172       Beam5R4PosTEC.push_back(theLaserPhi.at(60 + theIteration * 434));
00173       Beam5R4PosTEC.push_back(theLaserPhiError.at(60 + theIteration * 434));
00174       Beam5R4PosTEC.push_back(theLaserPhi.at(61 + theIteration * 434));
00175       Beam5R4PosTEC.push_back(theLaserPhiError.at(61 + theIteration * 434));
00176       Beam5R4PosTEC.push_back(theLaserPhi.at(62 + theIteration * 434));
00177       Beam5R4PosTEC.push_back(theLaserPhiError.at(62 + theIteration * 434));
00178       Beam5R4PosTEC.push_back(theLaserPhi.at(63 + theIteration * 434));
00179       Beam5R4PosTEC.push_back(theLaserPhiError.at(63 + theIteration * 434));
00180       Beam5R4PosTEC.push_back(theLaserPhi.at(64 + theIteration * 434));
00181       Beam5R4PosTEC.push_back(theLaserPhiError.at(64 + theIteration * 434));
00182       Beam5R4PosTEC.push_back(theLaserPhi.at(65 + theIteration * 434));
00183       Beam5R4PosTEC.push_back(theLaserPhiError.at(65 + theIteration * 434));
00184       Beam5R4PosTEC.push_back(theLaserPhi.at(66 + theIteration * 434));
00185       Beam5R4PosTEC.push_back(theLaserPhiError.at(66 + theIteration * 434));
00186       Beam5R4PosTEC.push_back(theLaserPhi.at(67 + theIteration * 434));
00187       Beam5R4PosTEC.push_back(theLaserPhiError.at(67 + theIteration * 434));
00188       Beam5R4PosTEC.push_back(theLaserPhi.at(68 + theIteration * 434));
00189       Beam5R4PosTEC.push_back(theLaserPhiError.at(68 + theIteration * 434));
00190           
00191       Beam6R4PosTEC.push_back(theLaserPhi.at(69 + theIteration * 434));
00192       Beam6R4PosTEC.push_back(theLaserPhiError.at(69 + theIteration * 434));
00193       Beam6R4PosTEC.push_back(theLaserPhi.at(70 + theIteration * 434));
00194       Beam6R4PosTEC.push_back(theLaserPhiError.at(70 + theIteration * 434));
00195       Beam6R4PosTEC.push_back(theLaserPhi.at(71 + theIteration * 434));
00196       Beam6R4PosTEC.push_back(theLaserPhiError.at(71 + theIteration * 434));
00197       Beam6R4PosTEC.push_back(theLaserPhi.at(72 + theIteration * 434));
00198       Beam6R4PosTEC.push_back(theLaserPhiError.at(72 + theIteration * 434));
00199       Beam6R4PosTEC.push_back(theLaserPhi.at(73 + theIteration * 434));
00200       Beam6R4PosTEC.push_back(theLaserPhiError.at(73 + theIteration * 434));
00201       Beam6R4PosTEC.push_back(theLaserPhi.at(74 + theIteration * 434));
00202       Beam6R4PosTEC.push_back(theLaserPhiError.at(74 + theIteration * 434));
00203       Beam6R4PosTEC.push_back(theLaserPhi.at(75 + theIteration * 434));
00204       Beam6R4PosTEC.push_back(theLaserPhiError.at(75 + theIteration * 434));
00205       Beam6R4PosTEC.push_back(theLaserPhi.at(76 + theIteration * 434));
00206       Beam6R4PosTEC.push_back(theLaserPhiError.at(76 + theIteration * 434));
00207       Beam6R4PosTEC.push_back(theLaserPhi.at(77 + theIteration * 434));
00208       Beam6R4PosTEC.push_back(theLaserPhiError.at(77 + theIteration * 434));
00209           
00210           
00211       Beam7R4PosTEC.push_back(theLaserPhi.at(83 + theIteration * 434));
00212       Beam7R4PosTEC.push_back(theLaserPhiError.at(83 + theIteration * 434));
00213       Beam7R4PosTEC.push_back(theLaserPhi.at(84 + theIteration * 434));
00214       Beam7R4PosTEC.push_back(theLaserPhiError.at(84 + theIteration * 434));
00215       Beam7R4PosTEC.push_back(theLaserPhi.at(85 + theIteration * 434));
00216       Beam7R4PosTEC.push_back(theLaserPhiError.at(85 + theIteration * 434));
00217       Beam7R4PosTEC.push_back(theLaserPhi.at(86 + theIteration * 434));
00218       Beam7R4PosTEC.push_back(theLaserPhiError.at(86 + theIteration * 434));
00219       Beam7R4PosTEC.push_back(theLaserPhi.at(87 + theIteration * 434));
00220       Beam7R4PosTEC.push_back(theLaserPhiError.at(87 + theIteration * 434));
00221       Beam7R4PosTEC.push_back(theLaserPhi.at(88 + theIteration * 434));
00222       Beam7R4PosTEC.push_back(theLaserPhiError.at(88 + theIteration * 434));
00223       Beam7R4PosTEC.push_back(theLaserPhi.at(89 + theIteration * 434));
00224       Beam7R4PosTEC.push_back(theLaserPhiError.at(89 + theIteration * 434));
00225       Beam7R4PosTEC.push_back(theLaserPhi.at(90 + theIteration * 434));
00226       Beam7R4PosTEC.push_back(theLaserPhiError.at(90 + theIteration * 434));
00227       Beam7R4PosTEC.push_back(theLaserPhi.at(91 + theIteration * 434));
00228       Beam7R4PosTEC.push_back(theLaserPhiError.at(91 + theIteration * 434));
00229           
00230           
00231       Beam0R6PosTEC.push_back(theLaserPhi.at(97 + theIteration * 434));
00232       Beam0R6PosTEC.push_back(theLaserPhiError.at(97 + theIteration * 434));
00233       Beam0R6PosTEC.push_back(theLaserPhi.at(98 + theIteration * 434));
00234       Beam0R6PosTEC.push_back(theLaserPhiError.at(98 + theIteration * 434));
00235       Beam0R6PosTEC.push_back(theLaserPhi.at(99 + theIteration * 434));
00236       Beam0R6PosTEC.push_back(theLaserPhiError.at(99 + theIteration * 434));
00237       Beam0R6PosTEC.push_back(theLaserPhi.at(100 + theIteration * 434));
00238       Beam0R6PosTEC.push_back(theLaserPhiError.at(100 + theIteration * 434));
00239       Beam0R6PosTEC.push_back(theLaserPhi.at(101 + theIteration * 434));
00240       Beam0R6PosTEC.push_back(theLaserPhiError.at(101 + theIteration * 434));
00241       Beam0R6PosTEC.push_back(theLaserPhi.at(102 + theIteration * 434));
00242       Beam0R6PosTEC.push_back(theLaserPhiError.at(102 + theIteration * 434));
00243       Beam0R6PosTEC.push_back(theLaserPhi.at(103 + theIteration * 434));
00244       Beam0R6PosTEC.push_back(theLaserPhiError.at(103 + theIteration * 434));
00245       Beam0R6PosTEC.push_back(theLaserPhi.at(104 + theIteration * 434));
00246       Beam0R6PosTEC.push_back(theLaserPhiError.at(104 + theIteration * 434));
00247       Beam0R6PosTEC.push_back(theLaserPhi.at(105 + theIteration * 434));
00248       Beam0R6PosTEC.push_back(theLaserPhiError.at(105 + theIteration * 434));
00249           
00250       Beam1R6PosTEC.push_back(theLaserPhi.at(106 + theIteration * 434));
00251       Beam1R6PosTEC.push_back(theLaserPhiError.at(106 + theIteration * 434));
00252       Beam1R6PosTEC.push_back(theLaserPhi.at(107 + theIteration * 434));
00253       Beam1R6PosTEC.push_back(theLaserPhiError.at(107 + theIteration * 434));
00254       Beam1R6PosTEC.push_back(theLaserPhi.at(108 + theIteration * 434));
00255       Beam1R6PosTEC.push_back(theLaserPhiError.at(108 + theIteration * 434));
00256       Beam1R6PosTEC.push_back(theLaserPhi.at(109 + theIteration * 434));
00257       Beam1R6PosTEC.push_back(theLaserPhiError.at(109 + theIteration * 434));
00258       Beam1R6PosTEC.push_back(theLaserPhi.at(110 + theIteration * 434));
00259       Beam1R6PosTEC.push_back(theLaserPhiError.at(110 + theIteration * 434));
00260       Beam1R6PosTEC.push_back(theLaserPhi.at(111 + theIteration * 434));
00261       Beam1R6PosTEC.push_back(theLaserPhiError.at(111 + theIteration * 434));
00262       Beam1R6PosTEC.push_back(theLaserPhi.at(112 + theIteration * 434));
00263       Beam1R6PosTEC.push_back(theLaserPhiError.at(112 + theIteration * 434));
00264       Beam1R6PosTEC.push_back(theLaserPhi.at(113 + theIteration * 434));
00265       Beam1R6PosTEC.push_back(theLaserPhiError.at(113 + theIteration * 434));
00266       Beam1R6PosTEC.push_back(theLaserPhi.at(114 + theIteration * 434));
00267       Beam1R6PosTEC.push_back(theLaserPhiError.at(114 + theIteration * 434));
00268           
00269       Beam2R6PosTEC.push_back(theLaserPhi.at(115 + theIteration * 434));
00270       Beam2R6PosTEC.push_back(theLaserPhiError.at(115 + theIteration * 434));
00271       Beam2R6PosTEC.push_back(theLaserPhi.at(116 + theIteration * 434));
00272       Beam2R6PosTEC.push_back(theLaserPhiError.at(116 + theIteration * 434));
00273       Beam2R6PosTEC.push_back(theLaserPhi.at(117 + theIteration * 434));
00274       Beam2R6PosTEC.push_back(theLaserPhiError.at(117 + theIteration * 434));
00275       Beam2R6PosTEC.push_back(theLaserPhi.at(118 + theIteration * 434));
00276       Beam2R6PosTEC.push_back(theLaserPhiError.at(118 + theIteration * 434));
00277       Beam2R6PosTEC.push_back(theLaserPhi.at(119 + theIteration * 434));
00278       Beam2R6PosTEC.push_back(theLaserPhiError.at(119 + theIteration * 434));
00279       Beam2R6PosTEC.push_back(theLaserPhi.at(120 + theIteration * 434));
00280       Beam2R6PosTEC.push_back(theLaserPhiError.at(120 + theIteration * 434));
00281       Beam2R6PosTEC.push_back(theLaserPhi.at(121 + theIteration * 434));
00282       Beam2R6PosTEC.push_back(theLaserPhiError.at(121 + theIteration * 434));
00283       Beam2R6PosTEC.push_back(theLaserPhi.at(122 + theIteration * 434));
00284       Beam2R6PosTEC.push_back(theLaserPhiError.at(122 + theIteration * 434));
00285       Beam2R6PosTEC.push_back(theLaserPhi.at(123 + theIteration * 434));
00286       Beam2R6PosTEC.push_back(theLaserPhiError.at(123 + theIteration * 434));
00287           
00288       Beam3R6PosTEC.push_back(theLaserPhi.at(124 + theIteration * 434));
00289       Beam3R6PosTEC.push_back(theLaserPhiError.at(124 + theIteration * 434));
00290       Beam3R6PosTEC.push_back(theLaserPhi.at(125 + theIteration * 434));
00291       Beam3R6PosTEC.push_back(theLaserPhiError.at(125 + theIteration * 434));
00292       Beam3R6PosTEC.push_back(theLaserPhi.at(126 + theIteration * 434));
00293       Beam3R6PosTEC.push_back(theLaserPhiError.at(126 + theIteration * 434));
00294       Beam3R6PosTEC.push_back(theLaserPhi.at(127 + theIteration * 434));
00295       Beam3R6PosTEC.push_back(theLaserPhiError.at(127 + theIteration * 434));
00296       Beam3R6PosTEC.push_back(theLaserPhi.at(128 + theIteration * 434));
00297       Beam3R6PosTEC.push_back(theLaserPhiError.at(128 + theIteration * 434));
00298       Beam3R6PosTEC.push_back(theLaserPhi.at(129 + theIteration * 434));
00299       Beam3R6PosTEC.push_back(theLaserPhiError.at(129 + theIteration * 434));
00300       Beam3R6PosTEC.push_back(theLaserPhi.at(130 + theIteration * 434));
00301       Beam3R6PosTEC.push_back(theLaserPhiError.at(130 + theIteration * 434));
00302       Beam3R6PosTEC.push_back(theLaserPhi.at(131 + theIteration * 434));
00303       Beam3R6PosTEC.push_back(theLaserPhiError.at(131 + theIteration * 434));
00304       Beam3R6PosTEC.push_back(theLaserPhi.at(132 + theIteration * 434));
00305       Beam3R6PosTEC.push_back(theLaserPhiError.at(132 + theIteration * 434));
00306           
00307       Beam4R6PosTEC.push_back(theLaserPhi.at(133 + theIteration * 434));
00308       Beam4R6PosTEC.push_back(theLaserPhiError.at(133 + theIteration * 434));
00309       Beam4R6PosTEC.push_back(theLaserPhi.at(134 + theIteration * 434));
00310       Beam4R6PosTEC.push_back(theLaserPhiError.at(134 + theIteration * 434));
00311       Beam4R6PosTEC.push_back(theLaserPhi.at(135 + theIteration * 434));
00312       Beam4R6PosTEC.push_back(theLaserPhiError.at(135 + theIteration * 434));
00313       Beam4R6PosTEC.push_back(theLaserPhi.at(136 + theIteration * 434));
00314       Beam4R6PosTEC.push_back(theLaserPhiError.at(136 + theIteration * 434));
00315       Beam4R6PosTEC.push_back(theLaserPhi.at(137 + theIteration * 434));
00316       Beam4R6PosTEC.push_back(theLaserPhiError.at(137 + theIteration * 434));
00317       Beam4R6PosTEC.push_back(theLaserPhi.at(138 + theIteration * 434));
00318       Beam4R6PosTEC.push_back(theLaserPhiError.at(138 + theIteration * 434));
00319       Beam4R6PosTEC.push_back(theLaserPhi.at(139 + theIteration * 434));
00320       Beam4R6PosTEC.push_back(theLaserPhiError.at(139 + theIteration * 434));
00321       Beam4R6PosTEC.push_back(theLaserPhi.at(140 + theIteration * 434));
00322       Beam4R6PosTEC.push_back(theLaserPhiError.at(140 + theIteration * 434));
00323       Beam4R6PosTEC.push_back(theLaserPhi.at(141 + theIteration * 434));
00324       Beam4R6PosTEC.push_back(theLaserPhiError.at(141 + theIteration * 434));
00325           
00326       Beam5R6PosTEC.push_back(theLaserPhi.at(142 + theIteration * 434));
00327       Beam5R6PosTEC.push_back(theLaserPhiError.at(142 + theIteration * 434));
00328       Beam5R6PosTEC.push_back(theLaserPhi.at(143 + theIteration * 434));
00329       Beam5R6PosTEC.push_back(theLaserPhiError.at(143 + theIteration * 434));
00330       Beam5R6PosTEC.push_back(theLaserPhi.at(144 + theIteration * 434));
00331       Beam5R6PosTEC.push_back(theLaserPhiError.at(144 + theIteration * 434));
00332       Beam5R6PosTEC.push_back(theLaserPhi.at(145 + theIteration * 434));
00333       Beam5R6PosTEC.push_back(theLaserPhiError.at(145 + theIteration * 434));
00334       Beam5R6PosTEC.push_back(theLaserPhi.at(146 + theIteration * 434));
00335       Beam5R6PosTEC.push_back(theLaserPhiError.at(146 + theIteration * 434));
00336       Beam5R6PosTEC.push_back(theLaserPhi.at(147 + theIteration * 434));
00337       Beam5R6PosTEC.push_back(theLaserPhiError.at(147 + theIteration * 434));
00338       Beam5R6PosTEC.push_back(theLaserPhi.at(148 + theIteration * 434));
00339       Beam5R6PosTEC.push_back(theLaserPhiError.at(148 + theIteration * 434));
00340       Beam5R6PosTEC.push_back(theLaserPhi.at(149 + theIteration * 434));
00341       Beam5R6PosTEC.push_back(theLaserPhiError.at(149 + theIteration * 434));
00342       Beam5R6PosTEC.push_back(theLaserPhi.at(150 + theIteration * 434));
00343       Beam5R6PosTEC.push_back(theLaserPhiError.at(150 + theIteration * 434));
00344           
00345       Beam6R6PosTEC.push_back(theLaserPhi.at(151 + theIteration * 434));
00346       Beam6R6PosTEC.push_back(theLaserPhiError.at(151 + theIteration * 434));
00347       Beam6R6PosTEC.push_back(theLaserPhi.at(152 + theIteration * 434));
00348       Beam6R6PosTEC.push_back(theLaserPhiError.at(152 + theIteration * 434));
00349       Beam6R6PosTEC.push_back(theLaserPhi.at(153 + theIteration * 434));
00350       Beam6R6PosTEC.push_back(theLaserPhiError.at(153 + theIteration * 434));
00351       Beam6R6PosTEC.push_back(theLaserPhi.at(154 + theIteration * 434));
00352       Beam6R6PosTEC.push_back(theLaserPhiError.at(154 + theIteration * 434));
00353       Beam6R6PosTEC.push_back(theLaserPhi.at(155 + theIteration * 434));
00354       Beam6R6PosTEC.push_back(theLaserPhiError.at(155 + theIteration * 434));
00355       Beam6R6PosTEC.push_back(theLaserPhi.at(156 + theIteration * 434));
00356       Beam6R6PosTEC.push_back(theLaserPhiError.at(156 + theIteration * 434));
00357       Beam6R6PosTEC.push_back(theLaserPhi.at(157 + theIteration * 434));
00358       Beam6R6PosTEC.push_back(theLaserPhiError.at(157 + theIteration * 434));
00359       Beam6R6PosTEC.push_back(theLaserPhi.at(158 + theIteration * 434));
00360       Beam6R6PosTEC.push_back(theLaserPhiError.at(158 + theIteration * 434));
00361       Beam6R6PosTEC.push_back(theLaserPhi.at(159 + theIteration * 434));
00362       Beam6R6PosTEC.push_back(theLaserPhiError.at(159 + theIteration * 434));
00363           
00364       Beam7R6PosTEC.push_back(theLaserPhi.at(160 + theIteration * 434));
00365       Beam7R6PosTEC.push_back(theLaserPhiError.at(160 + theIteration * 434));
00366       Beam7R6PosTEC.push_back(theLaserPhi.at(161 + theIteration * 434));
00367       Beam7R6PosTEC.push_back(theLaserPhiError.at(161 + theIteration * 434));
00368       Beam7R6PosTEC.push_back(theLaserPhi.at(162 + theIteration * 434));
00369       Beam7R6PosTEC.push_back(theLaserPhiError.at(162 + theIteration * 434));
00370       Beam7R6PosTEC.push_back(theLaserPhi.at(163 + theIteration * 434));
00371       Beam7R6PosTEC.push_back(theLaserPhiError.at(163 + theIteration * 434));
00372       Beam7R6PosTEC.push_back(theLaserPhi.at(164 + theIteration * 434));
00373       Beam7R6PosTEC.push_back(theLaserPhiError.at(164 + theIteration * 434));
00374       Beam7R6PosTEC.push_back(theLaserPhi.at(165 + theIteration * 434));
00375       Beam7R6PosTEC.push_back(theLaserPhiError.at(165 + theIteration * 434));
00376       Beam7R6PosTEC.push_back(theLaserPhi.at(166 + theIteration * 434));
00377       Beam7R6PosTEC.push_back(theLaserPhiError.at(166 + theIteration * 434));
00378       Beam7R6PosTEC.push_back(theLaserPhi.at(167 + theIteration * 434));
00379       Beam7R6PosTEC.push_back(theLaserPhiError.at(167 + theIteration * 434));
00380       Beam7R6PosTEC.push_back(theLaserPhi.at(168 + theIteration * 434));
00381       Beam7R6PosTEC.push_back(theLaserPhiError.at(168 + theIteration * 434));
00382 
00383       /* *************************************** */
00384 
00385       // create entry in the map for Beam 0 in Ring 4 of TEC+
00386       theBeamPosition["Beam 0 in Ring 4 in TEC+"] = Beam0R4PosTEC;
00387       // create entry in the map for Beam 1 in Ring 4 of TEC+
00388       theBeamPosition["Beam 1 in Ring 4 in TEC+"] = Beam1R4PosTEC;
00389       // create entry in the map for Beam 2 in Ring 4 of TEC+
00390       theBeamPosition["Beam 2 in Ring 4 in TEC+"] = Beam2R4PosTEC;
00391       // create entry in the map for Beam 3 in Ring 4 of TEC+
00392       theBeamPosition["Beam 3 in Ring 4 in TEC+"] = Beam3R4PosTEC;
00393       // create entry in the map for Beam 4 in Ring 4 of TEC+
00394       theBeamPosition["Beam 4 in Ring 4 in TEC+"] = Beam4R4PosTEC;
00395       // create entry in the map for Beam 5 in Ring 4 of TEC+
00396       theBeamPosition["Beam 5 in Ring 4 in TEC+"] = Beam5R4PosTEC;
00397       // create entry in the map for Beam 6 in Ring 4 of TEC+
00398       theBeamPosition["Beam 6 in Ring 4 in TEC+"] = Beam6R4PosTEC;
00399       // create entry in the map for Beam 7 in Ring 4 of TEC+
00400       theBeamPosition["Beam 7 in Ring 4 in TEC+"] = Beam7R4PosTEC;
00401           
00402       // create entry in the map for Beam 0 in Ring 6 of TEC+
00403       theBeamPosition["Beam 0 in Ring 6 in TEC+"] = Beam0R6PosTEC;
00404       // create entry in the map for Beam 1 in Ring 6 of TEC+
00405       theBeamPosition["Beam 1 in Ring 6 in TEC+"] = Beam1R6PosTEC;
00406       // create entry in the map for Beam 2 in Ring 6 of TEC+
00407       theBeamPosition["Beam 2 in Ring 6 in TEC+"] = Beam2R6PosTEC;
00408       // create entry in the map for Beam 3 in Ring 6 of TEC+
00409       theBeamPosition["Beam 3 in Ring 6 in TEC+"] = Beam3R6PosTEC;
00410       // create entry in the map for Beam 4 in Ring 6 of TEC+
00411       theBeamPosition["Beam 4 in Ring 6 in TEC+"] = Beam4R6PosTEC;
00412       // create entry in the map for Beam 5 in Ring 6 of TEC+
00413       theBeamPosition["Beam 5 in Ring 6 in TEC+"] = Beam5R6PosTEC;
00414       // create entry in the map for Beam 6 in Ring 6 of TEC+
00415       theBeamPosition["Beam 6 in Ring 6 in TEC+"] = Beam6R6PosTEC;
00416       // create entry in the map for Beam 7 in Ring 6 of TEC+
00417       theBeamPosition["Beam 7 in Ring 6 in TEC+"] = Beam7R6PosTEC;
00418           
00419       // *******************************************************
00420 
00421       // add the beams to millipede
00422       // Ring 4 of TEC+
00423       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 0 in Ring 4 in TEC+"], 0, 4);
00424       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 1 in Ring 4 in TEC+"], 1, 4);
00425       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 2 in Ring 4 in TEC+"], 2, 4);
00426       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 3 in Ring 4 in TEC+"], 3, 4);
00427       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 4 in Ring 4 in TEC+"], 4, 4);
00428       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 5 in Ring 4 in TEC+"], 5, 4);
00429       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 6 in Ring 4 in TEC+"], 6, 4);
00430       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 7 in Ring 4 in TEC+"], 7, 4);
00431       // Ring 6 of TEC+
00432       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 0 in Ring 6 in TEC+"], 0, 6);
00433       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 1 in Ring 6 in TEC+"], 1, 6);
00434       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 2 in Ring 6 in TEC+"], 2, 6);
00435       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 3 in Ring 6 in TEC+"], 3, 6);
00436       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 4 in Ring 6 in TEC+"], 4, 6);
00437       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 5 in Ring 6 in TEC+"], 5, 6);
00438       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 6 in Ring 6 in TEC+"], 6, 6);
00439       theLaserAlignmentTrackerPosTEC->addLaserBeam(theBeamPosition["Beam 7 in Ring 6 in TEC+"], 7, 6);
00440     }
00441           
00442   // finally do the fit of the global parameters when we have reached the last iteration
00443   edm::LogInfo("LASAlingPosTEC") << "<LaserAlignmentPosTEC::alignment()>: doing the global fit ... ";
00444               
00445   theLaserAlignmentTrackerPosTEC->doGlobalFit(theAlignableTracker);
00446   
00447   // reset Millepede
00448 //   theLaserAlignmentTrackerPosTEC->resetMillepede(theAlignmentIteration);
00449 
00450   // delete the alignment Object for the positive TEC to avoid problems
00451   // with the alignment of the negative TEC
00452   if (theLaserAlignmentTrackerPosTEC != 0) { delete theLaserAlignmentTrackerPosTEC; }
00453 }
LaserAlignmentPosTEC.cc
