db/d30/TPCRateData_8m_source.html

close all


clear all


%%


% n = 1000;

% full_rate = 100e3;

% trig_rate = 15e3;

% trigger_window = 18e-6;


% n = 10000000;

n = 100000;

% n = 500000;

full_rate = 100e3;

% full_rate = 50e3;

trig_rate = 15e3;

% trigger_window = 13e-6;

trigger_window = 17.5e-6;

% trigger_window = 35e-6;

BCO = 10e6;


% Active region readout only

% 0.69


trigger_window_in_BCO = int64(trigger_window * BCO);

event_time_gap = ...

    exprnd(1/full_rate,n,1);


DataLoadBCO = zeros(int64(n/full_rate*BCO*1.1), 1);

TriggerLoadBCO = zeros(int64(n/full_rate*BCO*1.1), 1);


trig_event  = binornd(1,trig_rate./full_rate, n, 1);


DataRateBeforeTriggerPerkHz = (1/100e3)* 966. *.5*.6;


%%


recorded_data = 0;

triggered_data = 0;


n_recorded = 0;

n_triggered = 0;

current_time = 0;


for i = 1:n


    start_time_bin = int64(current_time*BCO)+1;

    % add data

    DataLoadBCO(start_time_bin:(start_time_bin+trigger_window_in_BCO)) =...

        DataLoadBCO(start_time_bin:(start_time_bin+trigger_window_in_BCO)) + 1;


    % trigger

    if (trig_event(i) == 1)

        current_trigger_window = trigger_window;


        TriggerLoadBCO(start_time_bin:(start_time_bin+trigger_window_in_BCO)) =...

            TriggerLoadBCO(start_time_bin:(start_time_bin+trigger_window_in_BCO)) + 1;


        n_triggered = n_triggered + 1;


    end


    % time elaps

    current_time = current_time + event_time_gap(i);


end


%%


% fprintf('Throttled event / total = %.3f; Throttled data / total = %.3f; Triggered event / total = %.3f; Triggered data / total = = %.3f\n',...

%     n_recorded/n, recorded_data/n, n_triggered./n, triggered_data/n);

fprintf('------------\n');


fprintf('full_rate = %.0f kHz; trig_rate= %.0f kHz; trigger_window= %.1f us; \n',...

    full_rate/1e3,trig_rate/1e3, trigger_window*1e6 );

fprintf('Trigger rate*drift window = %.2f;Full rate*drift window= %.2f;Trigger rate/full rate= %.2f; \n',...

    trig_rate*trigger_window,full_rate*trigger_window,trig_rate/full_rate  );


fprintf('throttled data / total = %.3f; Triggered data / total = %.3f; throttled/trigger = %.3f \n',...

    sum((TriggerLoadBCO >0) .* DataLoadBCO)/ sum(DataLoadBCO), sum(TriggerLoadBCO .* DataLoadBCO)/ sum(DataLoadBCO), ...

     sum((TriggerLoadBCO >0) .* DataLoadBCO/ sum(TriggerLoadBCO .* DataLoadBCO)));


fprintf('throttled data rate = %.0f Gbps; Triggered data rate = %.0f \n',...

    sum((TriggerLoadBCO >0) .* DataLoadBCO)/ sum(DataLoadBCO) *full_rate*DataRateBeforeTriggerPerkHz ,...

    sum(TriggerLoadBCO .* DataLoadBCO)/ sum(DataLoadBCO)*full_rate*DataRateBeforeTriggerPerkHz);

fprintf('------------\n');


% fprintf('throttled event / total = %.3f; Triggered event / total = = %.3f;  throttled/trigger = %.3f\n',...

%     recorded_data/n, triggered_data/n, recorded_data/triggered_data );


%%


% recorded_data = 0;

% triggered_data = 0;

%

% n_recorded = 0;

% n_triggered = 0;

% current_trigger_window = -1;

% recorded_window = -1;

%

% for i = 1:n

%

%     % trigger

%     if (trig_event(i) == 1)

%          current_trigger_window = trigger_window;

%     end

%

%     % readout

%     if (current_trigger_window>0)

% %         n_recorded = n_recorded + 1;

% %         n_triggered = n_triggered + 1;

%

%         if (current_trigger_window>=trigger_window)

%

%             if (recorded_window>=trigger_window)

%                 recorded_data = recorded_data +0;

%             elseif (recorded_window>=0)

%                 recorded_data = recorded_data + (trigger_window - recorded_window)/trigger_window;

%                 recorded_window = trigger_window;

%             else

%                 recorded_data = recorded_data + 1;

%                 recorded_window = trigger_window;

%             end

%

%         else % if (current_trigger_window>=trigger_window)

%

%             if (recorded_window>=current_trigger_window)

%                 recorded_data = recorded_data +0;

%             elseif (recorded_window>=0)

%                 recorded_data = recorded_data + (current_trigger_window - recorded_window)/trigger_window;

%                 recorded_window = current_trigger_window;

%             else

%                 recorded_data = recorded_data + (current_trigger_window )/trigger_window;

%                 recorded_window = current_trigger_window;

%             end

%         end % if (current_trigger_window>=trigger_window)

%

%     end % (current_trigger_window>0)

%

%     % time elaps

%     current_trigger_window = current_trigger_window - event_time_gap(i);

%     recorded_window = recorded_window - event_time_gap(i);

%

% end