optpilot-week.md

@@ -152,6 +152,56 @@ RRI2 probe _simulation_ histogram (for 5000 particles) with found peak locations
 ![RRI2-0038Y17](/uploads/b38cd4518844f0308ec0b500043ca647/RRI2-0038Y17.png)
 RRI2 probe _measurement_ with peak locations.
 
+Test 1
+=======
+First trial of optimizing three objectives:
+
+* difference between simulation and measurement of peak 1
+* difference between simulation and measurement of peak 2
+* difference between simulation and measurement of peak 3 up to 5
+
+Settings
+------------
+* 1000 generations,
+* 32 cores
+* 500 particles
+* no space charge
+* 7 turns
+
+```
+IF ( _optpilot_ ) {
+    // ---- OPTIMIZER SECTION -------
+    //dphiinit: DVAR, VARIABLE="phiinit", LOWERBOUND="106", UPPERBOUND="114";
+    //drinit: DVAR, VARIABLE="rinit", LOWERBOUND="2000", UPPERBOUND="2060";
+    //dprinit: DVAR, VARIABLE="prinit", LOWERBOUND="-0.02", UPPERBOUND="-0.01";
+    //dbenergy: DVAR, VARIABLE="benergy", LOWERBOUND="0.071", UPPERBOUND="0.072";
+    //dvars: DVARS=(dphiinit,drinit,dprinit, dbenergy);
+ 
+    //dpeak1:OBJECTIVE,EXPR="fabs(radialPeak("RRI2.peaks", 1) - radialPeak("RawData/RRI2-0038Y17.peaks", 1))";
+    //dpeak2:OBJECTIVE,EXPR="fabs(radialPeak("RRI2.peaks", 2) - radialPeak("RawData/RRI2-0038Y17.peaks", 2))";
+    //dpeak3_5:OBJECTIVE,EXPR="sumErrSqRadialPeak("RRI2.peaks", "RawData/RRI2-0038Y17.peaks", 3, 5)";
+ 
+    //objs: OBJECTIVES=(dpeak1,dpeak2,dpeak3_5);
+    //constrs: CONSTRAINTS = ();
+    //opt: OPTIMIZE, OBJECTIVES = objs, DVARS = dvars, CONSTRAINTS = constrs;
+ }
+```
+
+![dpeak1](/uploads/386181ad6663686058332c4e3f9382d4/dpeak1.png)
+
+Sum of all differences (simulation - measurement) of peak 1 per generation.
+
+![dpeak2](/uploads/93ae4469e41783713792cdc108629af2/dpeak2.png)
+
+Sum of all differences (simulation - measurement) of peak 2 per generation.
+
+![dpeak3_5](/uploads/4aa392ae2cef54cf6ad66d1dce16c623/dpeak3_5.png)
+
+Sum of all differences (simulation - measurement) of objective 3 per generation.
+
+We see a nice decreasing curve for peak 1. Also the curve for peak 2 is mainly decreasing. However, the third objective is also increasing over the generations.
+
+
 **Superconducting Optimisation**
 ===============
 Simple setup one dipole. Shift initial distribution in x-y plane and search for minimal emittance. Run 1000 generations.