--- /dev/null
+LHC_ENTRIES 64\r
+0 dip/acc/LHC/Beam/Energy.Energy i:1 10 3752 1258370110.680 3751 1258370111.940 3752 1258370112.200 3751 1258370113.460 3752 1258370114.980 3751 1005.240 3752 1006.500 3751 1007.760 3752 1008.280 3751 1009.540\r
+1 dip/acc/LHC/Beam/Intensity/Beam1.totalIntensity f:1 5 56037552334.431664 1258370110.0 56392545816.963661 1258370111.0 55815681407.849152 1258370112.0 56870307878.871315 1258370113.0 57936471638.075775 1258370114.0 \r
+2 dip/acc/LHC/Beam/Intensity/Beam2.totalIntensity f:1 5 71143942900.913422 1258370110.0 71092334378.715363 1258370111.0 70901421649.982712 1258370112.0 72992148850.006042 1258370113.0 72566669567.373245 1258370114.0\r
+3 dip/acc/LHC/Beam/LuminosityAverage/BRANB.4L2.meanLuminosity f:1 4 0.069952 1258370110.1 0.069253 1258370111.1 0.068560 1258370112.1 0.067874 1258370113.1 \r
+4 dip/acc/LHC/Beam/LuminosityAverage/BRANB.4R2.meanLuminosity f:1 4 0.253554 1258370110.1 0.251018 1258370111.1 0.248508 1258370112.1 0.246023 1258370113.1 \r
+5 dip/acc/LHC/RunControl/MachineMode.value s:1 2 BEAM SETUP 1260525283.679 CICCIO 126055300.000 \r
+6 dip/acc/LHC/RunControl/BeamMode.value s:1 1 STABLE BEAMS 1258370110.703 \r
+7 dip/acc/LHC/RunControl/BeamType/Beam1.payload i:1 1 0 1258370110.78 \r
+8 dip/acc/LHC/RunControl/BeamType/Beam2.payload i:1 1 0 1258370110.78 \r
+9 dip/acc/LHC/Machine/CollimatorPositions/TCLIA.4R2.lvdt_gap_downstream f:1 4 39.938952 1258370110.429 39.938831 1258370111.429 39.938908 1258370112.429 39.938830 1258370113.429\r
+10 dip/acc/LHC/Machine/CollimatorPositions/TCLIA.4R2.lvdt_gap_upstream f:1 4 39.839512 1258370110.429 39.839831 1258370111.429 39.839972 1258370112.429 39.839581 1258370113.429\r
+11 dip/acc/LHC/Machine/CollimatorPositions/TCLIA.4R2.lvdt_left_downstream f:1 4 19.989592 1258370110.429 19.989519 1258370111.429 19.989489 1258370112.429 19.989249 1258370113.429\r
+12 dip/acc/LHC/Machine/CollimatorPositions/TCLIA.4R2.lvdt_left_upstream f:1 4 19.920967 1258370110.429 19.920670 1258370111.429 19.920538 1258370112.429 19.920503 1258370113.429\r
+13 dip/acc/LHC/Machine/CollimatorPositions/TCLIA.4R2.lvdt_right_downstream f:1 4 -19.956480 1258370110.429 -19.956582 1258370111.429 -19.956627 1258370112.429 -19.956438 1258370113.429\r
+14 dip/acc/LHC/Machine/CollimatorPositions/TCLIA.4R2.lvdt_right_upstream f:1 4 -19.995873 1258370110.429 -19.995757 1258370111.429 -19.995164 1258370112.429 -19.994854 1258370113.429\r
+15 dip/acc/LHC/Beam/Intensity/Beam1.acqTimestamp i:1 5 326585588 1258370110.0 1326581913 1258370111.0 -1968387108 1258370112.0 -968384533 1258370113.0 31613692 1258370114.0\r
+16 dip/acc/LHC/Beam/Intensity/Beam1.primitiveLifetime f:1 5 0.014667 1258370110.0 0.011406 1258370111.0 -0.035169 1258370112.0 -0.006725 1258370113.0 0.009240 1258370114.0\r
+17 dip/acc/LHC/Beam/Intensity/Beam2.acqTimestamp i:1 5 327768588 1258370110.0 1327761613 1258370111.0 -1967204583 1258370112.0 -967202833 1258370113.0 32797142 1258370114.0 \r
+18 dip/acc/LHC/Beam/Intensity/Beam2.primitiveLifetime f:1 5 -0.007079 1258370110.0 0.018377 1258370111.0 -0.008751 1258370112.0 0.032130 1258370113.0 0.011080 1258370114.0\r
+19 dip/acc/LHC/Timing/PostMortem.oCounter i:1 1 389 1258370110.868 \r
+20 dip/acc/LHC/Timing/PostMortem.payload i:1 1 0 1258370110.868 \r
+21 dip/acc/LHC/Beam/Size/Beam1.acqCounter i:1 5 453238 1258370110.0 453239 1258370111.0 453240 1258370112.0 453241 1258370113.0 453242 1258370114.0\r
+22 dip/acc/LHC/Beam/Size/Beam1.acqState i:1 1 0 1258370110.0 \r
+23 dip/acc/LHC/Beam/Size/Beam1.amplitudeSet1 f:1 1 5330.270927 1258370110.0 \r
+24 dip/acc/LHC/Beam/Size/Beam1.amplitudeSet2 f:1 1 9958.275904 1258370110.0 \r
+25 dip/acc/LHC/Beam/Size/Beam1.nbDataSets i:1 1 2 1258370110.0 \r
+26 dip/acc/LHC/Beam/Size/Beam1.planeSet1 i:1 1 1 1258370110.0 \r
+27 dip/acc/LHC/Beam/Size/Beam1.planeSet2 i:1 1 2 1258370110.0 \r
+28 dip/acc/LHC/Beam/Size/Beam1.positionSet1 f:1 1 128.000000 1258370110.0 \r
+29 dip/acc/LHC/Beam/Size/Beam1.positionSet2 f:1 1 126.000000 1258370110.0 \r
+30 dip/acc/LHC/Beam/Size/Beam1.sigmaSet1 f:1 1 1.276596 1258370110.0 \r
+31 dip/acc/LHC/Beam/Size/Beam1.sigmaSet2 f:1 1 4.680851 1258370110.0 \r
+32 dip/acc/LHC/Beam/Size/Beam1.acqMsg s:1 5 OK 1258370110.0 OK 1258370111.0 OK 1258370112.0 OK 1258370113.0 OK 1258370114.0\r
+33 dip/acc/LHC/Beam/Size/Beam1.acqTime s:1 5 2009/12/12 20:42:56.000238 1258370110.0 2009/12/12 20:42:57.000238 1258370111.0 2009/12/12 20:42:58.000238 1258370112.0 2009/12/12 20:42:59.000238 1258370113.0 2009/12/12 20:43:00.000238 1258370114.0\r
+34 dip/acc/LHC/Beam/Size/Beam1.cycleName s:1 5 LHC.USER.LHC 1258370110.0 LHC.USER.LHC 1258370111.0 LHC.USER.LHC 1258370112.0 LHC.USER.LHC 1258370113.0 LHC.USER.LHC 1258370114.0\r
+35 dip/acc/LHC/Beam/Size/Beam1.cycleTime s:1 5 2009/12/12 20:42:57.001757 1258370110.0 2009/12/12 20:42:58.001761 1258370111.0 2009/12/12 20:42:59.001754 1258370112.0 2009/12/12 20:43:00.001779 1258370113.0 2009/12/12 20:43:01.001762 1258370114.0\r
+36 dip/acc/LHC/Beam/Size/Beam1.deviceName s:1 5 LHC.BSRTS.5R4.B1 1258370110.0 LHC.BSRTS.5R4.B1 1258370111.0 LHC.BSRTS.5R4.B1 1258370112.0 LHC.BSRTS.5R4.B1 1258370113.0 LHC.BSRTS.5R4.B1 1258370114.0\r
+37 dip/acc/LHC/Beam/Size/Beam2.acqCounter i:1 5 549085 1258370110.0 549086 1258370111.0 549087 1258370112.0 549088 1258370113.0 549089 1258370114.00\r
+38 dip/acc/LHC/Beam/Size/Beam2.acqState i:1 1 0 1258370110.0 \r
+39 dip/acc/LHC/Beam/Size/Beam2.amplitudeSet1 f:1 5 155157.385887 1258370110.0 159907.032000 1258370111.0 157872.861437 1258370112.0 156018.146620 1258370113.0 155337.050676 1258370114.0\r
+40 dip/acc/LHC/Beam/Size/Beam2.amplitudeSet2 f:1 5 114027.398515 1258370110.0 113327.081202 1258370111.0 113809.867556 1258370112.0 112864.620849 1258370113.0 113857.283803 1258370114.0\r
+41 dip/acc/LHC/Beam/Size/Beam2.nbDataSets i:1 1 2 1258370110.0 \r
+42 dip/acc/LHC/Beam/Size/Beam2.planeSet1 i:1 1 1 1258370110.0 \r
+43 dip/acc/LHC/Beam/Size/Beam2.planeSet2 i:1 1 2 1258370110.0 \r
+44 dip/acc/LHC/Beam/Size/Beam2.positionSet1 f:1 4 23.278800 1258370110.0 23.324400 1258370111.0 23.347200 1258370112.0 23.210400 1258370113.0\r
+45 dip/acc/LHC/Beam/Size/Beam2.positionSet2 f:1 5 18.959617 1258370110.0 18.966229 1258370111.0 18.967650 1258370112.0 18.948230 1258370113.0 18.933193 1258370114.0\r
+46 dip/acc/LHC/Beam/Size/Beam2.sigmaSet1 f:1 4 1.387404 1258370110.0 1.300085 1258370111.0 1.338894 1258370112.0 1.397106 1258370113.0\r
+47 dip/acc/LHC/Beam/Size/Beam2.sigmaSet2 f:1 5 2.475035 1258370110.0 2.523874 1258370111.0 2.482082 1258370112.0 2.541485 1258370113.0 2.519686 1258370114.0\r
+48 dip/acc/LHC/Beam/Size/Beam2.acqMsg s:1 5 OK 1258370110.0 OK 1258370111.0 OK 1258370112.0 OK 1258370113.0 OK 1258370114.0\r
+49 dip/acc/LHC/Beam/Size/Beam2.acqTime s:1 5 2009/12/12 20:42:56.000238 1258370110.0 2009/12/12 20:42:57.000238 1258370111.0 2009/12/12 20:42:58.000238 1258370112.0 2009/12/12 20:42:59.000238 1258370113.0 2009/12/12 20:43:00.000238 1258370114.0\r
+50 dip/acc/LHC/Beam/Size/Beam2.cycleName s:1 5 LHC.USER.LHC 1258370110.0 LHC.USER.LHC 1258370111.0 LHC.USER.LHC 1258370112.0 LHC.USER.LHC 1258370113.0 LHC.USER.LHC 1258370114.0\r
+51 dip/acc/LHC/Beam/Size/Beam2.cycleTime s:1 5 2009/12/12 20:42:57.001732 1258370110.0 2009/12/12 20:42:58.001739 1258370111.0 2009/12/12 20:42:59.001743 1258370112.0 2009/12/12 20:43:00.001736 1258370113.0 2009/12/12 20:43:01.002996 1258370114.0\r
+52 dip/acc/LHC/Beam/Size/Beam2.deviceName s:1 5 LHC.BSRTS.5L4.B2 1258370110.0 LHC.BSRTS.5L4.B2 1258370111.0 LHC.BSRTS.5L4.B2 1258370112.0 LHC.BSRTS.5L4.B2 1258370113.0 LHC.BSRTS.5L4.B2 1258370114.0\r
+53 dip/acc/LHC/Safety/MaxSPSProbeBeamA.probeBeamFlagAct i:1 1 0 1258370110.574 \r
+54 dip/acc/LHC/Safety/MaxSPSProbeBeamB.probeBeamFlagAct i:1 1 0 1258370111.651 \r
+55 dip/acc/LHC/Beam/LuminosityAverage/BRANB.4L2.acqMode i:1 3 5112 1258370110.1 5368 1258370111.1 5112 1258370112.1\r
+56 dip/acc/LHC/Beam/LuminosityAverage/BRANB.4L2.meanCrossingAngle f:1 1 0.000000 1258370110.1 \r
+57 dip/acc/LHC/Beam/LuminosityAverage/BRANB.4R2.acqMode i:1 3 5128 1258370110.1 4872 1258370111.1 5128 1258370112.1\r
+58 dip/acc/LHC/Beam/LuminosityAverage/BRANB.4R2.meanCrossingAngle f:1 1 0.000000 1258370110.1 \r
+59 dip/acc/LHC/RunControl/SafeBeam/Beam1.payload i:1 1 15 1258370110.782 \r
+60 dip/acc/LHC/RunControl/SafeBeam/Beam2.payload i:1 1 15 1258370110.468 \r
+61 dip/acc/LHC/RunControl/CirculatingBunchConfig/Beam1.value i:2808 1 1 1001 4231 22081 27731 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1258370111.33 \r
+62 dip/acc/LHC/RunControl/CirculatingBunchConfig/Beam2.value i:2808 1 1 1258370111 13141 18791 27731 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1258370112.41 \r
+63 dip/acc/LHC/RunControl/FillNumber.payload i:1 1 912 1258370111.60 \r
+==== END_OF_BLOCK ====\r
#include "AliCDBMetaData.h"
#include "AliLog.h"
#include "AliESDVertex.h"
+#include "AliLHCReader.h"
+#include "AliLHCData.h"
+#include "AliDCSArray.h"
class AliDCSValue;
class AliShuttleInterface;
const Int_t AliGRPPreprocessor::fgknDAQLbPar = 8; // num parameters in the logbook for PHYSICS runs, when beamType from DAQ logbook == NULL
const Int_t AliGRPPreprocessor::fgknDAQLbParReduced = 7; // num parameters in the logbook for the other cases
- const Int_t AliGRPPreprocessor::fgknDCSDP = 51; // number of dcs dps
+ const Int_t AliGRPPreprocessor::fgknDCSDP = 48; // number of dcs dps
const Int_t AliGRPPreprocessor::fgknDCSDPHallProbes = 40; // number of dcs dps
+ const Int_t AliGRPPreprocessor::fgknLHCDP = 5; // number of dcs dps from LHC data
const char* AliGRPPreprocessor::fgkDCSDataPoints[AliGRPPreprocessor::fgknDCSDP] = {
- "LHCState", // missing in DCS
"L3Polarity",
"DipolePolarity",
- "LHCLuminosity", // missing in DCS
- "BeamIntensity", // missing in DCS
"L3Current",
"DipoleCurrent",
"L3_BSF17_H1",
"Dipole_Outside_Temperature"
};
- const Short_t kSensors = 48; // start index position of sensor in DCS DPs
+ const Short_t kSensors = 45; // start index position of sensor in DCS DPs
const Short_t kNumSensors = 3; // Number of sensors in DCS DPs (CavernAtmosPressure, SurfaceAtmosPressure, CavernAtmosPressure2)
- const char* AliGRPPreprocessor::fgkLHCState[20] = {
- "P", "PREPARE",
- "J", "PREINJECTION",
- "I", "INJECTION",
- "F", "FILLING",
- "A", "ADJUST",
- "U", "UNSTABLE BEAMS",
- "S", "STABLE BEAMS",
- "D", "BEAM DUMP",
- "R", "RECOVER",
- "C", "PRECYCLE"
- };
-
+ const char* AliGRPPreprocessor::fgkLHCDataPoints[AliGRPPreprocessor::fgknLHCDP] = {
+ "dip/acc/LHC/Beam/Energy.Energy",
+ "dip/acc/LHC/RunControl/MachineMode.value",
+ "dip/acc/LHC/RunControl/BeamMode.value",
+ "dip/acc/LHC/RunControl/BeamType/Beam1.payload",
+ "dip/acc/LHC/RunControl/BeamType/Beam2.payload"
+ };
const char* kppError[] = {
"",
"(DAQ logbook ERROR)",
"(Trigger Configuration ERROR)",
"(DAQ logbook ERROR determining partition of the run)",
"(CTP timing ERROR)",
- "(SPD Mean Vertex ERROR)"
+ "(SPD Mean Vertex ERROR)",
+ "(DCS FXS Error for LHC Data)",
+ "(LHC Data Error)"
};
//_______________________________________________________________
AliGRPObject *grpobj = new AliGRPObject(); // object to store data
grpobj->SetBeamEnergyIsSqrtSHalfGeV(); // new format
+
//=================//
// DAQ logbook //
//=================//
//=================//
Log(Form("Starting DCS Query at %d and finishing at %d",GetStartTimeDCSQuery(),GetEndTimeDCSQuery()));
Int_t entries = ProcessDcsDPs( valueMap, grpobj );
- Log(Form("entries found = %d (should be %d)",entries, fgknDCSDP-4));
+ Log(Form("entries found = %d (should be %d)",entries, fgknDCSDP-1));
if (fdaqStartEndTimeOk){
- if( entries < fgknDCSDP-4 ) { // FIXME (!= ) LHState, LHCLuminosity, BeamIntensity, L3_BSF4_H3 are not working yet...
+ if( entries < fgknDCSDP-1 ) { // FIXME (!= ) L3_BSF4_H3 are not working yet...
Log(Form("Problem with the DCS data points!!! Only %d/%d entries found",entries,fgknDCSDP-4));
Log(Form("The DPs giving problems were:"));
for (Int_t iDP = 0; iDP < fgknDCSDP; iDP++){
TObjString *dpString = (TObjString*)ffailedDPs->At(iDP);
if (dpString){
TString name = dpString->String();
- if (name != "LHCState" && name != "LHCLuminosity" && name != "BeamIntensity" && name != "L3_BSF4_H3"){
+ if (name != "L3_BSF4_H3"){
Log(Form("******** %s ******** not present, but foreseen --> causing an ERROR",name.Data()));
}
else {
Log(Form("Incorrect field in DAQ logbook for partition = %s and detector = %s, going into error without trigger timing parameters...",partition.Data(),detector.Data()));
error |= 32;
}
- // storing AliGRPObject in OCDB
- AliCDBMetaData md;
- md.SetResponsible("Chiara Zampolli");
- md.SetComment("Output parameters from the GRP preprocessor.");
-
- Bool_t result = kTRUE;
- result = Store("GRP", "Data", grpobj, &md);
- delete grpobj;
+
+ //=================//
+ // LHC Data //
+ //=================//
+
+ UInt_t iLHCData = ProcessLHCData(grpobj);
+ if( iLHCData == 0 ) {
+ Log(Form("LHC Data from DCS FXS, successful!"));
+ } else if (iLHCData == 1) {
+ Log(Form("LHC Data, problems with DCS FXS!"));
+ error |= 256;
+ } else if (iLHCData ==3){
+ Log(Form("Problems in storing LHC Data - but not going into Error"));
+ } else{
+ Log(Form("LHC Data problems"));
+ error |= 512;
+ }
//==================//
// SPD Mean Vertex //
Log("SPD Mean Vertex not processed since runType != PHYSICS");
}
+ // storing AliGRPObject in OCDB
+
+ AliCDBMetaData md;
+ md.SetResponsible("Chiara Zampolli");
+ md.SetComment("Output parameters from the GRP preprocessor.");
+
+ Bool_t result = kTRUE;
+ result = Store("GRP", "Data", grpobj, &md);
+ delete grpobj;
+
if (result && !error ) {
Log("GRP Preprocessor Success");
return 0;
} else {
- Log( Form("GRP Preprocessor FAILS!!! %s%s%s%s%s%s%s%s",
+ Log( Form("GRP Preprocessor FAILS!!! %s%s%s%s%s%s%s%s%s%s",
kppError[(error&1)?1:0],
kppError[(error&2)?2:0],
kppError[(error&4)?3:0],
kppError[(error&16)?5:0],
kppError[(error&32)?6:0],
kppError[(error&64)?7:0],
- kppError[(error&128)?8:0]
+ kppError[(error&128)?8:0],
+ kppError[(error&256)?9:0],
+ kppError[(error&512)?10:0]
));
return error;
}
+
+
+}
+
+//_______________________________________________________________
+
+UInt_t AliGRPPreprocessor::ProcessLHCData(AliGRPObject *grpobj)
+{
+ //
+ //Getting the LHC Data from DCS FXS
+ //
+
+ TString timeStartString = (TString)GetRunParameter("DAQ_time_start");
+ TString timeEndString = (TString)GetRunParameter("DAQ_time_end");
+ if (timeStartString.IsNull() || timeStartString.IsNull()){
+ if (timeStartString.IsNull()){
+ AliError("DAQ_time_start not set in logbook! Setting statistical values for current DP to invalid");
+ }
+ else if (timeStartString.IsNull()){
+ AliError("DAQ_time_end not set in logbook! Setting statistical values for current DP to invalid");
+ }
+ return 2;
+ }
+
+ Double_t timeStart = timeStartString.Atof();
+ Double_t timeEnd = timeEndString.Atof();
+
+ // timeStart = 1260646960;
+ //timeEnd = 1260652740;
+
+ TString fileName = GetFile(kDCS, "LHCData","");
+ if (fileName.Length()>0){
+ AliInfo("Got The LHC Data file");
+ AliLHCReader* lhcReader = new AliLHCReader();
+ TMap* lhcMap = (TMap*)lhcReader->ReadLHCDP(fileName.Data());
+ if (lhcMap) {
+ Log(Form("LHCData map entries = %d",lhcMap->GetEntries()));
+
+ // Processing data to be put in AliGRPObject
+ TObjArray* energyArray = (TObjArray*)lhcMap->GetValue(fgkLHCDataPoints[0]);
+ Float_t energy = ProcessEnergy(energyArray,timeStart,timeEnd);
+ if (energy != -1) {
+ grpobj->SetBeamEnergy(energy);
+ grpobj->SetBeamEnergyIsSqrtSHalfGeV(kTRUE);
+ }
+
+ TObjArray* machineModeArray = (TObjArray*)lhcMap->GetValue(fgkLHCDataPoints[1]);
+ TObjArray* beamModeArray = (TObjArray*)lhcMap->GetValue(fgkLHCDataPoints[2]);
+ if (beamModeArray->GetEntries()!=1){
+ AliWarning("The beam mode changed! Setting it to UNKNOWN and storing array of strings");
+ grpobj->SetLHCState("UNKNOWN");
+ grpobj->SetLHCStateArray(beamModeArray);
+ }
+ else{
+ AliDCSArray* beamMode = (AliDCSArray*)beamModeArray->At(0);
+ TObjString* beamModeString = beamMode->GetStringArray(0);
+ AliInfo(Form("LHC State (corresponding to BeamMode) = %s",(beamModeString->String()).Data()));
+ grpobj->SetLHCState(beamModeString->String());
+ }
+
+ if (machineModeArray->GetEntries()!=1){
+ AliWarning("The machine mode changed! Setting it to UNKNOWN and storing array of strings");
+ grpobj->SetMachineMode("UNKNOWN");
+ grpobj->SetMachineModeArray(machineModeArray);
+ }
+ else{
+ AliDCSArray* machineMode = (AliDCSArray*)machineModeArray->At(0);
+ TObjString* machineModeString = machineMode->GetStringArray(0);
+ AliInfo(Form("Machine Mode = %s",(machineModeString->String()).Data()));
+ grpobj->SetMachineMode(machineModeString->String());
+ }
+ TObjArray* beam1Array = (TObjArray*)lhcMap->GetValue(fgkLHCDataPoints[3]);
+ TObjArray* beam2Array = (TObjArray*)lhcMap->GetValue(fgkLHCDataPoints[4]);
+
+ // Processing data to go to AliLHCData object
+ AliLHCData* dt = new AliLHCData(lhcMap,timeStart,timeEnd);
+
+ // storing AliLHCData in OCDB
+
+ AliCDBMetaData md;
+ md.SetResponsible("Ruben Shahoyan");
+ md.SetComment("LHC data from the GRP preprocessor.");
+
+ Bool_t result = kTRUE;
+ result = Store("GRP", "LHCData", dt, &md);
+ delete dt;
+ if (result) return 0;
+ else return 3;
+ }
+ else {
+ AliError("Cannot read correctly LHCData file");
+ return 2;
+ }
+ }
+
+ else {
+ AliError("No LHCData file found in DCS FXS");
+ return 1;
+ }
+
}
//_______________________________________________________________
UInt_t AliGRPPreprocessor::ProcessSPDMeanVertex()
{
+ //
//Getting the SPD Mean Vertex
+ //
+
TList* list = GetForeignFileSources("SPD", kDAQ, "VertexDiamond");
Bool_t storeResult = kTRUE;
if (list !=0x0 && list->GetEntries()!=0)
//
Int_t entries = 0; // counting the entries that are in the DCS DB, not taking care whether they have values or not
- Int_t nLHCEntries = 0;
Int_t nL3Entries = 0;
Int_t nDipoleEntries = 0;
Int_t nEnvEntries = 0;
Int_t nHallProbesEntries = 0;
- nLHCEntries = ProcessLHCDPs(valueMap, grpObj);
nL3Entries = ProcessL3DPs(valueMap, grpObj);
nDipoleEntries = ProcessDipoleDPs(valueMap, grpObj);
nEnvEntries = ProcessEnvDPs(valueMap, grpObj);
nHallProbesEntries = ProcessHPDPs(valueMap, grpObj);
grpObj->SetPolarityConventionLHC(); // after the dipole cables swap we comply with LHC convention
- Log(Form("nLHCEntries = %d, L3Entries = %d, nDipoleEntries =%d, nEnvEntries = %d, nHallProbesEntries = %d", nLHCEntries, nL3Entries, nDipoleEntries, nEnvEntries, nHallProbesEntries));
- entries = nLHCEntries + nL3Entries + nDipoleEntries + nEnvEntries + nHallProbesEntries;
+ Log(Form("L3Entries = %d, nDipoleEntries =%d, nEnvEntries = %d, nHallProbesEntries = %d", nL3Entries, nDipoleEntries, nEnvEntries, nHallProbesEntries));
+ entries = nL3Entries + nDipoleEntries + nEnvEntries + nHallProbesEntries;
return entries;
}
return nHPEntries;
}
-//_______________________________________________________________
-
-Int_t AliGRPPreprocessor::ProcessLHCDPs(const TMap* valueMap, AliGRPObject* grpObj)
-{
-
- //
- // processing of LHC related DCS DPs, i.e.:
- // LHCState
- // LHCLuminosity
- // BeamIntensity
- //
-
- Int_t nLHCEntries = 0;
- TObjArray *array = 0x0;
- Int_t indexDP = -1;
-
- AliInfo(Form("==========LHCState==========="));
- indexDP = kLHCState;
- array = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[indexDP]);
- if(!array) {
- Log(Form("%s not found in the map!!!",fgkDCSDataPoints[indexDP]));
- }
- else {
- if (array->GetEntries() == 0){
- AliError(Form("No entries found in array! setting %s to invalid...",fgkDCSDataPoints[indexDP]));
- }
- else {
- TString stringDCS = ProcessChar(array);
- if (stringDCS.Length()!=0) {
- Bool_t found = kFALSE;
- for( Int_t i=0; i<20; i+=2 ) {
- if( stringDCS.CompareTo(fgkLHCState[i]) == 0 ) {
- stringDCS = fgkLHCState[i+1];
- found = kTRUE;
- break;
- }
- }
- if (found){
- Log(Form("<%s> for run %d: %s",fgkDCSDataPoints[indexDP],fRun, stringDCS.Data()));
- grpObj->SetLHCState(stringDCS);
- }
- else{
- Log(Form("%s values found not valid!",fgkDCSDataPoints[indexDP]));
- grpObj->SetLHCState(AliGRPObject::GetInvalidString());
- }
- }
- else {
- Log(Form("%s not valid (null length), string set as invalid!",fgkDCSDataPoints[indexDP]));
- grpObj->SetLHCState(AliGRPObject::GetInvalidString());
- }
- }
- ffailedDPs->RemoveAt(indexDP);
- nLHCEntries++;
- }
-
- if (array) array = 0x0;
-
- AliInfo(Form("==========LHCLuminosity==========="));
- Bool_t outOfRange = kFALSE; // flag to monitor if any value collected by DCS is out of range
- indexDP = kLHCLuminosity;
- array = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[indexDP]);
- if(!array) {
- Log(Form("%s not found in the map!!!",fgkDCSDataPoints[indexDP]));
- }
- else {
- if (array->GetEntries() == 0){
- AliError(Form("No entries found in array! setting %s and its Spline Fit to invalid...",fgkDCSDataPoints[indexDP]));
- }
- else {
- Float_t *floatDCS = ProcessFloatAll(array);
- if (floatDCS != NULL){
- grpObj->SetLHCLuminosity(floatDCS);
- AliSplineFit* splfit = GetSplineFit(array,fgkDCSDataPoints[indexDP]);
- grpObj->SetLHCLuminositySplineFit(splfit);
- // delete splfit;
- }
- else {
- outOfRange = kTRUE;
- }
- if (floatDCS){
- delete[] floatDCS;
- floatDCS = 0x0;
- }
- }
- if (!outOfRange) {
- ffailedDPs->RemoveAt(indexDP);
- nLHCEntries++;
- }
- }
-
- if (array) array = 0x0;
-
- AliInfo(Form("==========BeamIntensity==========="));
- if (outOfRange) outOfRange = kFALSE; // resetting outOfRange if needed
- indexDP = kBeamIntensity;
- array = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[indexDP]);
- if(!array) {
- Log(Form("%s not found in the map!!!",fgkDCSDataPoints[indexDP]));
- }
- else {
- if (array->GetEntries() == 0){
- AliError(Form("No entries found in array! setting %s and its Spline Fit to invalid...",fgkDCSDataPoints[indexDP]));
- }
- else {
- Float_t *floatDCS = ProcessFloatAll(array);
- if (floatDCS != NULL){
- grpObj->SetBeamIntensity(floatDCS);
- AliSplineFit* splfit1 = GetSplineFit(array,fgkDCSDataPoints[indexDP]);
- grpObj->SetBeamIntensitySplineFit(splfit1);
- //delete splfit;
- }
- else{
- outOfRange = kTRUE;
- }
- if (floatDCS){
- delete[] floatDCS;
- floatDCS = 0x0;
- }
- }
- if (!outOfRange) {
- nLHCEntries++;
- ffailedDPs->RemoveAt(indexDP);
- }
- }
-
- return nLHCEntries;
-}
//_________________________________________________________________________
AliSplineFit* AliGRPPreprocessor::GetSplineFit(const TObjArray *array, const TString& stringID){
return lastRun;
}
+//-----------------------------------------------------------------
+Double_t AliGRPPreprocessor::CalculateMean(TObjArray* array){
+
+ //
+ // Calculating mean over TObjArray from LHC Data
+ //
+
+ TString timeStartString = (TString)GetRunParameter("DAQ_time_start");
+ TString timeEndString = (TString)GetRunParameter("DAQ_time_end");
+ if (timeStartString.IsNull() || timeStartString.IsNull()){
+ if (timeStartString.IsNull()){
+ AliError("DAQ_time_start not set in logbook! Setting statistical values for current DP to invalid");
+ }
+ else if (timeStartString.IsNull()){
+ AliError("DAQ_time_end not set in logbook! Setting statistical values for current DP to invalid");
+ }
+ return 0;
+ }
+
+ Int_t timeStart = (Int_t)(timeStartString.Atoi());
+ Int_t timeEnd = (Int_t)(timeEndString.Atoi());
+ timeStart = 1260646960;
+ timeEnd = 1260652740;
+ Double_t* parameters = new Double_t[5];
+ parameters[0] = -1.;
+ parameters[1] = -1.;
+ parameters[2] = -1.;
+ parameters[3] = -1.;
+ parameters[4] = -1.;
+ Int_t iCounts = 0;
+ Int_t iCountsRun = 0;
+ Int_t nCounts = array->GetEntries();
+ printf("ncounts = %d\n",nCounts);
+ Double_t valueBeforeSOR = 0;
+ Double_t valueAfterEOR = 0;
+ Double_t timestampBeforeSOR = -1.;
+ Double_t timestampAfterEOR = -1.;
+ Int_t ientrySOR = -1;
+ Int_t ientryEOR = -1;
+ Double_t* arrayValues = 0x0;
+ Double_t* arrayWeights = 0x0;
+ Bool_t truncMeanFlag = kTRUE; // flag to indicate whether Truncated Mean should be calculated or not
+ Bool_t sdFlag = kTRUE; // flag to indicate whether SD (wrt Mean/Median) should be calculated or not
+
+ for(Int_t i = 0; i < nCounts; i++) {
+ AliDCSArray *dcs = (AliDCSArray*)array->At(i);
+ if((dcs->GetTimeStamp() >= timeStart) &&(dcs->GetTimeStamp() <= timeEnd)) {
+ AliDebug(2,Form("%d-th entry = %f at timestamp %f\n",i,(Double_t)(dcs->GetInt(0)),dcs->GetTimeStamp()));
+ iCounts += 1;
+ // look for the last value before SOR and the first value before EOR
+ if ((dcs->GetTimeStamp() >= timeStart) && (dcs->GetTimeStamp() < timeStart)) {
+ timestampBeforeSOR = dcs->GetTimeStamp();
+ AliDebug(2,Form("timestamp of last value before SOR = %f, with DAQ_time_start = %d\n",timestampBeforeSOR,timeStart));
+ valueBeforeSOR = (Double_t)(dcs->GetInt(0));
+ }
+ else if ((dcs->GetTimeStamp() <= timeEnd) && (dcs->GetTimeStamp() > timeEnd) && timestampAfterEOR == -1){
+ timestampAfterEOR = dcs->GetTimeStamp();
+ valueAfterEOR = (Double_t)(dcs->GetInt(0));
+ AliDebug(2,Form("timestamp of first value after EOR = %f, with DAQ_time_end = %d\n",timestampAfterEOR,timeEnd));
+ }
+ // check if there are DPs between DAQ_time_start and DAQ_time_end
+ if((dcs->GetTimeStamp() >= timeStart) &&(dcs->GetTimeStamp() <= timeEnd)) {
+ if (ientrySOR == -1) ientrySOR = i; // first entry after SOR
+ if (ientryEOR < i) ientryEOR = i; // last entry before EOR
+ AliDebug(2,Form("entry between SOR and EOR\n"));
+ iCountsRun += 1;
+ }
+ }
+ else {
+ printf("DCS values for the parameter outside the queried interval: timestamp = %f\n",dcs->GetTimeStamp());
+ }
+ }
+
+ if (timestampBeforeSOR == -1.){
+ printf("No value found before SOR\n");
+ }
+ if (timestampAfterEOR == -1.){
+ printf("No value found after EOR\n");
+ }
+
+ printf("Number of valid entries (within DCS query interval) = %i, from a total amount of %i entries\n",iCounts,nCounts);
+ printf("Last value before DAQ_time_start (SOR) = %f at timestamp = %f\n",valueBeforeSOR,timestampBeforeSOR);
+ printf("First value after DAQ_time_end (EOR) = %f at timestamp = %f\n",valueAfterEOR,timestampAfterEOR);
+ printf("Found %d entries between DAQ_time_start (SOR) and DAQ_time_end (EOR)\n",iCountsRun);
+ printf("Index of first entry after DAQ_time_start (SOR) = %d\n ",ientrySOR);
+ printf("Index of first entry before DAQ_time_end (EOR) = %d\n ",ientryEOR);
+
+ Int_t nentriesUsed = 0;
+ if (iCountsRun > 1){
+ printf("Using entries between DAQ_time_start (SOR) and DAQ_time_end (EOR)\n");
+ printf("Calculating (weighted) Mean and Median\n" );
+ arrayValues = new Double_t[iCountsRun];
+ arrayWeights = new Double_t[iCountsRun];
+ nentriesUsed = iCountsRun;
+ for (Int_t i = ientrySOR; i <= ientryEOR; i++){
+ AliDCSArray *dcs = (AliDCSArray *)array->At(i);
+ Double_t timestamp2 = 0;
+ if (i < ientryEOR){
+ AliDCSArray *dcs1 = (AliDCSArray *)array->At(i+1);
+ timestamp2 = dcs1->GetTimeStamp();
+ }
+ else {
+ timestamp2 = (Double_t)timeEnd+1;
+ }
+ arrayWeights[i-ientrySOR] = (Double_t)((Double_t)timestamp2 - dcs->GetTimeStamp());
+ arrayValues[i-ientrySOR] = (Double_t)(dcs->GetInt(0));
+ printf("Entry %d: value = %f, weight = %f\n",i-ientrySOR,arrayValues[i-ientrySOR],arrayWeights[i-ientrySOR]);
+ }
+ parameters[0] = TMath::Mean(iCountsRun,arrayValues,arrayWeights);
+ parameters[2] = TMath::Median(iCountsRun,arrayValues,arrayWeights);
+ }
+ else if (iCountsRun == 1){
+ AliDCSArray* dcs = (AliDCSArray *)array->At(ientrySOR);
+ nentriesUsed = 2;
+ if (timestampBeforeSOR != -1 && timestampBeforeSOR != (Int_t)dcs->GetTimeStamp()){
+ printf("Using single entry between DAQ_time_start (SOR) and DAQ_time_end (EOR) and last entry before SOR. Truncated mean won't be calculated.\n");
+ arrayValues = new Double_t[2];
+ arrayWeights = new Double_t[2];
+ arrayValues[0] = valueBeforeSOR;
+ arrayWeights[0] = (Double_t)(dcs->GetTimeStamp()-(Double_t)timestampBeforeSOR);
+ arrayValues[1] = (Double_t)(dcs->GetInt(0));
+ arrayWeights[1] = (Double_t)((Double_t)timeEnd+1-dcs->GetTimeStamp());
+ printf("value0 = %f, with weight = %f\n",arrayValues[0],arrayWeights[0]);
+ printf("value1 = %f, with weight = %f\n",arrayValues[1],arrayWeights[1]);
+ parameters[0] = TMath::Mean(2,arrayValues,arrayWeights);
+ parameters[2] = TMath::Median(2,arrayValues,arrayWeights);
+ truncMeanFlag = kFALSE;
+ }
+ else{
+ printf("Cannot calculate mean, truncated mean, median, SD wrt mean, SD wrt median for current DP - only one value collected during the run, but no value before with which to calculate the statistical quantities\n");
+ parameters[0] = -1;
+ parameters[1] = -1;
+ parameters[2] = -1;
+ parameters[3] = -1;
+ parameters[4] = -1;
+ return parameters[0];
+ }
+ }
+ else { // iCountsRun == 0, using only the point immediately before SOR
+ if (timestampBeforeSOR == -1.){
+ printf("Cannot set mean, truncated mean, median, SD wrt mean, SD wrt median for current DP - no points during the run collected, and point before SOR missing\n");
+ parameters[0] = -1;
+ parameters[1] = -1;
+ parameters[2] = -1;
+ parameters[3] = -1;
+ parameters[4] = -1;
+ return parameters[0];
+ }
+ else {
+ printf("Using only last entry before SOR. Truncated mean and Standard deviations (wrt mean/median) won't be calculated.\n");
+ printf("value = %f\n",valueBeforeSOR);
+ parameters[0] = valueBeforeSOR;
+ parameters[2] = valueBeforeSOR;
+ truncMeanFlag = kFALSE;
+ sdFlag = kFALSE;
+ }
+ }
+
+ Double_t temp = 0;
+ Double_t temp1 = 0;
+ Double_t sumweights = 0;
+ Int_t entriesTruncMean = 0;
+ Double_t* arrayValuesTruncMean = new Double_t[nentriesUsed];
+ Double_t* arrayWeightsTruncMean = new Double_t[nentriesUsed];
+
+ // calculating SD wrt Mean and Median
+ printf("Calculating SD wrt Mean and SD wrt Median\n");
+ if (sdFlag){
+ for (Int_t i =0; i< nentriesUsed; i++){
+ //printf("Entry %d: value = %f, weight = %f\n",i,arrayValues[i],arrayWeights[i]);
+ temp += (arrayValues[i]-parameters[2])*(arrayValues[i]-parameters[2]);
+ temp1 += arrayWeights[i]*(arrayValues[i]-parameters[0])*(arrayValues[i]-parameters[0]);
+ sumweights += arrayWeights[i];
+ }
+ // setting SD wrt Mean
+ if (sumweights != 0 ){
+ parameters[3] = TMath::Sqrt(temp1/sumweights);
+ }
+ else {
+ printf("Sum of weights to calculate Standard Deviation (wrt mean) <= 0, setting the SD to invalid\n");
+ parameters[3] = -1;
+ }
+ // setting SD wrt Median
+ if (nentriesUsed != 0){
+ parameters[4] = TMath::Sqrt(temp/nentriesUsed);
+ }
+ else{
+ printf("Number of entries used to calculate Standard Deviation (wrt median) <= 0, setting the SD to invalid\n");
+ parameters[4] = -1;
+ }
+ }
+ else {
+ parameters[3] = -1;
+ parameters[4] = -1;
+ }
+
+ // calculating truncated mean (this comes afterwards since you need the SD wrt Mean)
+ if (truncMeanFlag){
+ printf("Calculating Truncated Mean\n");
+ for (Int_t i =0; i< nentriesUsed; i++){
+ //printf("Entry %d: value = %f, weight = %f\n",i,arrayValues[i],arrayWeights[i]);
+ if ((arrayValues[i]<=parameters[0]+3*parameters[3]) && (arrayValues[i]>=parameters[0]-3*parameters[3])){
+ arrayValuesTruncMean[entriesTruncMean]=arrayValues[i];
+ arrayWeightsTruncMean[entriesTruncMean]=arrayWeights[i];
+ printf("For Truncated Mean: Entry %d: value = %f, weight = %f\n",entriesTruncMean,arrayValuesTruncMean[entriesTruncMean],arrayWeightsTruncMean[entriesTruncMean]);
+ entriesTruncMean++;
+ }
+ else{
+ printf("Discarding entry\n");
+ }
+ }
+ // setting truncated mean
+ if (entriesTruncMean >1){
+ printf("%d entries used for truncated mean\n",entriesTruncMean);
+ parameters[1] = TMath::Mean(entriesTruncMean,arrayValuesTruncMean,arrayWeightsTruncMean);
+ }
+ else{
+ printf("Too few entries (%d) to calculate truncated mean\n",entriesTruncMean);
+ parameters[1] = -1;
+ }
+ }
+ else{
+ parameters[1] = -1;
+ }
+
+ printf("(weighted) mean = %f \n",parameters[0]);
+ printf("(weighted) truncated mean = %f \n",parameters[1]);
+ printf("median = %f \n",parameters[2]);
+ printf("(weighted) standard deviation with (weighted) mean = %f \n",parameters[3]);
+ printf("standard deviation with median = %f \n",parameters[4]);
+
+ return (parameters[0]);
+}
+//------------------------------------------------------------------------------------------------------
+Float_t AliGRPPreprocessor::ProcessEnergy(TObjArray* array, Double_t timeStart, Double_t timeEnd){
+
+ //
+ // Method to processo LHC Energy information
+ // Only the first value is returned, provided that it is withing DAQ_time_start and DAQ_time_end
+ //
+
+ Int_t nCounts = array->GetEntries();
+ Float_t energy = -1;
+ AliDebug(2,Form("Energy measurements = %d\n",nCounts));
+ for(Int_t i = 0; i < nCounts; i++) {
+ AliDCSArray *dcs = (AliDCSArray*)array->At(i);
+ if((dcs->GetTimeStamp() >= timeStart) &&(dcs->GetTimeStamp() <= timeEnd)) {
+ energy = (Float_t)(TMath::Nint(((Double_t)(dcs->GetInt(0)))*120/1000)); // sqrt(s)/2 energy in GeV
+ AliInfo(Form("Energy value found = %d, converting --> sqrt(s)/2 = %f (GeV)", dcs->GetInt(0),energy));
+ break;
+ }
+ else {
+ AliError("No energy values found between DAQ_time_start and DAQ_time_end - energy will remain invalid!");
+ }
+ }
+
+ return energy;
+}