Changeset 607 in ETALON for SPESO


Ignore:
Timestamp:
May 10, 2016, 9:17:32 AM (8 years ago)
Author:
malovyts
Message:

Added new variables

Location:
SPESO/ana2015
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • SPESO/ana2015/data_get_vals.m

    r512 r607  
    3434for iloop=1:nloop
    3535switch vartype(iloop)
    36     case dvar.BEAMINFO
    37         vals(iloop,:)=allData.beamInfo(idx,varidx(iloop))'; %'
    38     case dvar.DATAINFO
    39             if (iloop==1)
    40               loopvar=var1;
    41             else
     36        case dvar.BEAMINFO
     37                vals(iloop,:)=allData.beamInfo(idx,varidx(iloop))'; %'
     38        case dvar.DATAINFO
     39                if (iloop==1)
     40                        loopvar=var1;
     41           else
    4242              loopvar=var2;
    43             end
    44         switch loopvar
    45             case dvar.DATA_VAL_LPM
    46             valsused(iloop,1)=1;
    47             valsused(iloop+2,1)=1;
    48                 vals(iloop,:)=allData.ch1signalValLPM(idx);
    49                 vals(iloop+2,:)=allData.ch2signalValLPM(idx);
    50             case dvar.DATA_PEAK_LPM
    51             valsused(iloop,1)=1;
    52             valsused(iloop+2,1)=1;
    53                 vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
    54                 vals(iloop+2,:)=allData.ch2signalNegPeakLPM(idx);
    55             case dvar.DATA_VAL_SPM
    56             valsused(iloop,1)=1;
    57             valsused(iloop+2,1)=1;
    58                 vals(iloop,:)=allData.ch1signalValSPM(idx);
    59                 vals(iloop+2,:)=allData.ch2signalValSPM(idx);
    60             case dvar.DATA_PEAK_SPM
    61             valsused(iloop,1)=1;
    62             valsused(iloop+2,1)=1;
    63                 vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);
    64                 vals(iloop+2,:)=allData.ch2signalNegPeakSPM(idx);
     43      end
     44      switch loopvar
     45                        case dvar.DATA_VAL_LPM
     46                                valsused(iloop,1)=1;
     47                                valsused(iloop+2,1)=1;
     48                                vals(iloop,:)=allData.ch1signalValLPM(idx);
     49                                vals(iloop+2,:)=allData.ch2signalValLPM(idx);
     50                        case dvar.DATA_PEAK_LPM
     51                                valsused(iloop,1)=1;
     52                                valsused(iloop+2,1)=1;
     53                                vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
     54                                vals(iloop+2,:)=allData.ch2signalNegPeakLPM(idx);
     55                        case dvar.DATA_VAL_SPM
     56                                valsused(iloop,1)=1;
     57                                valsused(iloop+2,1)=1;
     58                                vals(iloop,:)=allData.ch1signalValSPM(idx);
     59                                vals(iloop+2,:)=allData.ch2signalValSPM(idx);
     60                        case dvar.DATA_PEAK_SPM
     61                                valsused(iloop,1)=1;
     62                                valsused(iloop+2,1)=1;
     63                                vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);
     64                                vals(iloop+2,:)=allData.ch2signalNegPeakSPM(idx);
    6565
    66             case dvar.DATA1_VAL_LPM
    67                 valsused(iloop,1)=1;
    68                 vals(iloop,:)=allData.ch1signalValLPM(idx);
    69             case dvar.DATA1_PEAK_LPM
    70                 valsused(iloop,1)=1;
    71                 vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
    72             case dvar.DATA1_VAL_SPM
    73                 valsused(iloop,1)=1;
    74                 vals(iloop,:)=allData.ch1signalValSPM(idx);
    75             case dvar.DATA1_PEAK_SPM
    76                 valsused(iloop,1)=1;
    77                 vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);
     66                        case dvar.DATA1_VAL_LPM
     67                                valsused(iloop,1)=1;
     68                                vals(iloop,:)=allData.ch1signalValLPM(idx);
     69                        case dvar.DATA1_PEAK_LPM
     70                                valsused(iloop,1)=1;
     71                                vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
     72                        case dvar.DATA1_VAL_SPM
     73                                valsused(iloop,1)=1;
     74                                vals(iloop,:)=allData.ch1signalValSPM(idx);
     75                        case dvar.DATA1_PEAK_SPM
     76                                valsused(iloop,1)=1;
     77                                vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);
    7878
    79             case dvar.DATA2_VAL_LPM
    80                 valsused(iloop,1)=1;
    81                 vals(iloop,:)=allData.ch2signalValLPM(idx);
    82             case dvar.DATA2_PEAK_LPM
    83                 valsused(iloop,1)=1;
    84                 vals(iloop,:)=allData.ch2signalNegPeakLPM(idx);
    85             case dvar.DATA2_VAL_SPM
    86                 valsused(iloop,1)=1;
    87                 vals(iloop,:)=allData.ch2signalValSPM(idx);
    88             case dvar.DATA2_PEAK_SPM
    89                 valsused(iloop,1)=1;
    90                 vals(iloop,:)=allData.ch2signalNegPeakSPM(idx);
    91             % Select time in minutes   
    92             case dvar.DATA_TIME_MINUTES
    93                 valsused(iloop,1)=1;
    94                 %convert time HHMMSS to minutes
    95                 vals(iloop,:)=( fix(allData.hoursNum(idx)/10000)*60+...
     79                        case dvar.DATA2_VAL_LPM
     80                                valsused(iloop,1)=1;
     81                                vals(iloop,:)=allData.ch2signalValLPM(idx);
     82                        case dvar.DATA2_PEAK_LPM
     83                                valsused(iloop,1)=1;
     84                                vals(iloop,:)=allData.ch2signalNegPeakLPM(idx);
     85                        case dvar.DATA2_VAL_SPM
     86                                valsused(iloop,1)=1;
     87                                vals(iloop,:)=allData.ch2signalValSPM(idx);
     88                        case dvar.DATA2_PEAK_SPM
     89                                valsused(iloop,1)=1;
     90                                vals(iloop,:)=allData.ch2signalNegPeakSPM(idx);
     91                        % Select time in minutes       
     92                        case dvar.DATA_TIME_MINUTES
     93                                valsused(iloop,1)=1;
     94                                %convert time HHMMSS to minutes
     95                                vals(iloop,:)=( fix(allData.hoursNum(idx)/10000)*60+...
    9696                                fix(mod(allData.hoursNum(idx),10000)/100)+...
    9797                                mod(allData.hoursNum(idx),100)/60);
    98             case dvar.DATA_TIME_30
    99                 valsused(iloop,1)=1;
    100                 %convert time HHMMSS to minutes
    101                 vals(iloop,:)=mod(( fix(allData.hoursNum(idx)/10000)*60+...
     98                        case dvar.DATA_TIME_30
     99            valsused(iloop,1)=1;
     100            %convert time HHMMSS to minutes
     101            vals(iloop,:)=mod(( fix(allData.hoursNum(idx)/10000)*60+...
    102102                                fix(mod(allData.hoursNum(idx),10000)/100)+...
    103103                                mod(allData.hoursNum(idx),100)/60),30);
    104104           
    105             case dvar.DATA_RS1_OFFSET
    106                 RS_params=get_RS_params(allData.timestamp(1,idx));
    107                 vals(iloop,:)=allData.beamInfo(idx,5)'-RS_params.RS1_offset;
    108             case dvar.DATA_RS2_OFFSET
    109                 RS_params=get_RS_params(allData.timestamp(1,idx));
    110                 vals(iloop,:)=allData.beamInfo(idx,7)'-RS_params.RS2_offset;
    111             case dvar.DATA_RS3_OFFSET
    112                 RS_params=get_RS_params(allData.timestamp(1,idx));
    113                 vals(iloop,:)=allData.beamInfo(idx,9)'-RS_params.RS3_offset;
    114             case dvar.DATA_RS4_OFFSET
    115                 RS_params=get_RS_params(allData.timestamp(1,idx));
    116                 vals(iloop,:)=allData.beamInfo(idx,11)'-RS_params.RS4_offset;
    117             otherwise
    118                 warning_sp([ 'The varidx ' varidx ' is unknown.' ]);
    119                 error('end');
    120             end %switch
     105                        case dvar.DATA_RS1_OFFSET
     106                                RS_params=get_RS_params(allData.timestamp(1,idx));
     107                                vals(iloop,:)=allData.beamInfo(idx,5)'-RS_params.RS1_offset;
     108                        case dvar.DATA_RS2_OFFSET
     109                                RS_params=get_RS_params(allData.timestamp(1,idx));
     110                                vals(iloop,:)=allData.beamInfo(idx,7)'-RS_params.RS2_offset;
     111                        case dvar.DATA_RS3_OFFSET
     112                                RS_params=get_RS_params(allData.timestamp(1,idx));
     113                                vals(iloop,:)=allData.beamInfo(idx,9)'-RS_params.RS3_offset;
     114                        case dvar.DATA_RS4_OFFSET
     115                                RS_params=get_RS_params(allData.timestamp(1,idx));
     116                                vals(iloop,:)=allData.beamInfo(idx,11)'-RS_params.RS4_offset;
     117                       
     118                        case dvar.DATA_DET1_THETA
     119                                SPOS_offset=-645000; % position of the grating in TS
     120                                X0_pos=220000; % position of the grating in TX
     121                                raw_th=atan((X0_pos-allData.beamInfo(idx,15))./...
     122                                                                (allData.beamInfo(idx,13)-SPOS_offset)./...
     123                                                                0.6.*0.5)*180/pi;
     124                                % 0.6 - coeff between TS and mm
     125                                % 0.5 - coeff between TX and mm
     126                                vals(iloop,:)=(-raw_th+(1+sign(raw_th))/2*180);
     127                                 % conver from 0:-pi/2 - pi/2:0 to 0:pi
     128                        case dvar.DATA_DET2_THETA
     129                                SPOS_offset=-470000; % position of the grating in TS
     130                                X0_pos=220000; % position of the grating in TX
     131                                raw_th=atan((X0_pos-allData.beamInfo(idx,15))./...
     132                                                                (allData.beamInfo(idx,13)-SPOS_offset)./...
     133                                                                0.6.*0.5)*180/pi;
     134                                % 0.6 - coeff between TS and mm
     135                                % 0.5 - coeff between TX and mm
     136                                vals(iloop,:)=(-raw_th+(1+sign(raw_th))/2*180);
     137                                 % conver from 0:-pi/2 - pi/2:0 to 0:pi
     138                        case dvar.DATA_DET1_PHI
     139                                ZPOS_offset=-66000; % position of the grating in TS
     140                                X0_pos=220000; % position of the grating in TX
     141                                raw_th=atan((X0_pos-allData.beamInfo(idx,15))./...
     142                                                                (allData.beamInfo(idx,17)-ZPOS_offset)./...
     143                                                                2.48.*0.5)*180/pi;
     144                                % 2.48 - coeff between TZ and mm
     145                                % 0.5 - coeff between TX and mm
     146                                vals(iloop,:)=raw_th-sign(raw_th)*90;
     147                                 % conver from -0:pi/2 - pi/2:0 to -pi/2:pi/2
     148                        case dvar.DATA_DET2_PHI
     149                                ZPOS_offset=-66000; % position of the grating in TS
     150                                X0_pos=220000; % position of the grating in TX
     151                                raw_th=atan((X0_pos-allData.beamInfo(idx,15))./...
     152                                                                (allData.beamInfo(idx,17)-ZPOS_offset)./...
     153                                                                2.48.*0.5)*180/pi;
     154                                % 2.48 - coeff between TZ and mm
     155                                % 0.5 - coeff between TX and mm
     156                                vals(iloop,:)=raw_th-sign(raw_th)*90;
     157                                 % conver from 0:-pi/2 - pi/2:0 to -pi/2:pi/2
     158
     159                        otherwise
     160                                warning_sp([ 'The varidx ' varidx ' is unknown.' ]);
     161                                error('end');
     162                        end %switch
    121163    otherwise
    122         warning_sp([ 'The vartype ' vartype ' is unknown.' ]);
    123         error('end');
    124 end %switch
     164                        warning_sp([ 'The vartype ' vartype ' is unknown.' ]);
     165                        error('end');
     166        end %switch
    125167end %for
    126168

  • SPESO/ana2015/define_data_var.m

    r504 r607  
    4747dvar.DATA_RS4_OFFSET=dvar.DATAINFO+44;
    4848
     49dvar.DATA_DET1_THETA=dvar.DATAINFO+45;
     50dvar.DATA_DET2_THETA=dvar.DATAINFO+46;
     51dvar.DATA_DET1_PHI=dvar.DATAINFO+47;
     52dvar.DATA_DET2_PHI=dvar.DATAINFO+48;
     53
    4954dvar.TIMESTAMP=dvar.RAWINFO+1;
    5055dvar.RAW_HOUR=dvar.RAWINFO+2;

  • SPESO/ana2015/get_var_name.m

    r504 r607  
    55
    66switch(varidx)
    7     case dvar.BTIME
    8         str='Time (sec.)';       
    9     case dvar.RESEAU_POS
    10         str='Pos. reseau';       
    11     case dvar.DETECT_XPOS
    12         str='Pos X detector';
    13     case dvar.DETECT_ZPOS
    14         str='Pos Z detector';
    15     case dvar.DETECT_SPOS
    16         str='Pos S detector';
    17     case dvar.DETECT_RS1
    18         str='Rotation R1';
    19     case dvar.DETECT_RS2
    20         str='Rotation R2';
    21     case dvar.DETECT_RS3
    22         str='Rotation R3';
    23     case dvar.DETECT_RS4
    24         str='Rotation R4';
     7        case dvar.BTIME
     8                str='Time (sec.)';       
     9        case dvar.RESEAU_POS
     10                str='Pos. reseau';       
     11        case dvar.DETECT_XPOS
     12                str='Pos X detector';
     13        case dvar.DETECT_ZPOS
     14                str='Pos Z detector';
     15        case dvar.DETECT_SPOS
     16                str='Pos S detector';
     17        case dvar.DETECT_RS1
     18                str='Rotation R1';
     19        case dvar.DETECT_RS2
     20                str='Rotation R2';
     21        case dvar.DETECT_RS3
     22                str='Rotation R3';
     23        case dvar.DETECT_RS4
     24                str='Rotation R4';
    2525
    26     case dvar.XBPM
    27         str='BPM X position';
    28     case dvar.QICT1
    29         str='Q on ICT1';
    30     case dvar.QICT2
    31         str='Q on ICT2';
    32     case dvar.LPM
    33         str='LPM event';
     26        case dvar.XBPM
     27                str='BPM X position';
     28        case dvar.QICT1
     29                str='Q on ICT1';
     30        case dvar.QICT2
     31                str='Q on ICT2';
     32        case dvar.LPM
     33                str='LPM event';
    3434
    35     case dvar.DATA_VAL_LPM
    36         str='Signal value LPM';
    37     case dvar.DATA_PEAK_LPM
    38         str='Signal peak LPM';
    39     case dvar.DATA_VAL_SPM
    40         str='Signal value SPM';
    41     case dvar.DATA_PEAK_SPM
    42         str='Signal peak SPM';
     35        case dvar.DATA_VAL_LPM
     36                str='Signal value LPM';
     37        case dvar.DATA_PEAK_LPM
     38                str='Signal peak LPM';
     39        case dvar.DATA_VAL_SPM
     40                str='Signal value SPM';
     41        case dvar.DATA_PEAK_SPM
     42                str='Signal peak SPM';
    4343
    44     case dvar.DATA1_VAL_LPM
    45         str='Signal value LPM det. 1';
    46     case dvar.DATA1_PEAK_LPM
    47         str='Signal peak LPM det. 1';
    48     case dvar.DATA1_VAL_SPM
    49         str='Signal value SPM det. 1';
    50     case dvar.DATA1_PEAK_SPM
    51         str='Signal peak SPM det. 1';
     44        case dvar.DATA1_VAL_LPM
     45                str='Signal value LPM det. 1';
     46        case dvar.DATA1_PEAK_LPM
     47                str='Signal peak LPM det. 1';
     48        case dvar.DATA1_VAL_SPM
     49                str='Signal value SPM det. 1';
     50        case dvar.DATA1_PEAK_SPM
     51                str='Signal peak SPM det. 1';
    5252
    53     case dvar.DATA2_VAL_LPM
    54         str='Signal value LPM det. 2';
    55     case dvar.DATA2_PEAK_LPM
    56         str='Signal peak LPM det. 2';
    57     case dvar.DATA2_VAL_SPM
    58         str='Signal value SPM det. 2';
    59     case dvar.DATA2_PEAK_SPM
    60         str='Signal peak SPM det. 2';
     53        case dvar.DATA2_VAL_LPM
     54                str='Signal value LPM det. 2';
     55        case dvar.DATA2_PEAK_LPM
     56                str='Signal peak LPM det. 2';
     57        case dvar.DATA2_VAL_SPM
     58                str='Signal value SPM det. 2';
     59        case dvar.DATA2_PEAK_SPM
     60                str='Signal peak SPM det. 2';
    6161
    62     case dvar.DATA_TIME_MINUTES
    63         str='Time (minutes)';
    64     case dvar.DATA_TIME_30
    65         str='Time (minutes in half hour)';
     62        case dvar.DATA_TIME_MINUTES
     63                str='Time (minutes)';
     64        case dvar.DATA_TIME_30
     65                str='Time (minutes in half hour)';
    6666
    67     case dvar.DATA_RS1_OFFSET
    68         str='Theta, deg';
    69     case dvar.DATA_RS2_OFFSET
    70         str='Theta, deg';
    71     case dvar.DATA_RS3_OFFSET
    72         str='Phi, deg';
    73     case dvar.DATA_RS4_OFFSET
    74         str='Phi, deg';
     67        case dvar.DATA_RS1_OFFSET
     68                str='RS1 Theta, deg';
     69        case dvar.DATA_RS2_OFFSET
     70                str='RS2 Theta, deg';
     71        case dvar.DATA_RS3_OFFSET
     72                str='RS3 Phi, deg';
     73        case dvar.DATA_RS4_OFFSET
     74                str='RS4 Phi, deg';
    7575
    76     case dvar.TIMESTAMP
    77         str='Timestamp (day,sec)';       
    78     case dvar.RAW_HOUR
    79         str='Raw hour of Day (hhmmss)';       
     76        case dvar.DATA_DET1_THETA
     77                str='Theta, deg';
     78        case dvar.DATA_DET2_THETA
     79                str='Theta, deg';
     80        case dvar.DATA_DET1_PHI
     81                str='Phi, deg';
     82        case dvar.DATA_DET2_PHI
     83                str='Phi, deg';
    8084
    81     case dvar.DAYSAGO
    82         str='Days since now';
     85        case dvar.TIMESTAMP
     86                str='Timestamp (day,sec)';       
     87        case dvar.RAW_HOUR
     88                str='Raw hour of Day (hhmmss)';       
    8389
    84     otherwise
    85         str='Unknown value';
     90        case dvar.DAYSAGO
     91                str='Days since now';
     92
     93        otherwise
     94                str='Unknown value';
    8695end
    8796
Note: See TracChangeset for help on using the changeset viewer.