Search:
wiki:tuning_the_carls
Last modified 16 months ago Last modified on 05/22/12 12:48:45

    CARL*TRIG tuning after Mammoth                      5/20/99
    ------------------------------

    Many carlsubtrig triggers for the Mammoth sequence have been very
    long (sometimes maximum duration) and have contained multiple
    small events. Rick would've preferred that many these big triggers
    had been a bunch of single-event triggers.  He's given me a few
    examples of triggers that he thought CUSP would've split into more
    triggers. Let's see if we can teach carl*trig to make them into a
    series of smaller triggers.  We'll grab the waveforms and try replaying
    them with tankplayer!  I'll grab a chuck of time a bit longer than the
    original triggers (tack some on the beginning and end).

    Then we should run these data thru carl*trig with the same settings
    as ew03/ew04 and see what results we get.  Problem is, we have to adjust
    the timestamp on the data to match the computer's system clock so that
    carl*trig will work. Then we'll have to mentally adjust the times on the
    output so we can compare them with the original results.  We'll
    document the hypocenters and triggers from the original real-time runs.
    All files for these tests will live on campbell in ~dietz/carltrig.

    First some basics on how carlsubtrig works.
    Three carlsubtrig parameters control the duration of the network triggers.
    Here are the original ew03/ew04 settings for those parameters:

    NetTriggerDur    10   # Nominal duration (sec) for the base network trigger
    SubnetContrib    15   # Addition contribution (sec) to network trigger.
                          #   duration for each concurrently triggered subnet.
    MaxDuration     500   # Maximum duration (sec) allowed for network trigger.

    The current network trigger duration (let's call it NewTrig) is calculated
    once each second based on the number of currently triggered subnets (nSubTrig)
    using the following equation:

      NewTrig = NetTriggerDur + ( nSubTrig * SubnetContrib )

    The total network duration (TotalNetDuration) is a slightly different
    beast, defined as the time at which the first subnet(s) triggered until
    some time after the last subnet(s) untrigger, or until TotalNetDuration
    reaches MaxDuration.  The total network duration can reach MaxDuration
    in one of two ways:

    1. Many subnets are triggered at the same time. If

         nSubTrig >= ( MaxDuration - NetTriggerDur )/ SubnetContrib

       then the TotalNetDuration will be be limited to MaxDuration.

    2. One or more subnets re-trigger at frequent intervals.
       In the one-subnet case, TotalNetDuration will keep increasing if

          ReTriggerInterval <= NetTriggerDur + SubnetContrib

       In the more general case, if a current network trigger consists
       of nSubTrig triggered subnets, then TotalNetDuration will be
       extended if at least one subnet retriggers within

          ReTriggerInterval <= NetTriggerDur + nSubTrig * SubnetContrib

    The maximum-duration triggers from the Mammoth sequence have been caused
    by both of these methods.

    ===========================================================================
    Here's the first original trigger and time range in question:

    src     id      original trigger    dur     time range saved for testing
    ----   -----  ---------------------------   ----------------------------
    ew04   25308  19990515 23:26:27.00  500 s   start:  19990515 23:25:27.00
                                                  end:  19990515 23:36:27.00
                                                 file:  19990515.2325.tbuf

    Original Hypocenters:
    a 19990515 2326 24.25 37 32.30 118 48.71   6.55 D 2.14 24 110  6.  0.06  0.3  0.6 BW   51082739 1
    b 19990515 2327 46.85 37 27.91 118 53.61  22.70 D 1.35  9 169  5.  0.20  3.3  4.5 CW   51082740 1
    c 19990515 2329 48.57 37 29.55 118 48.66   6.90 D 1.08  6 195  7.  0.01  1.4  2.7 CW   51082741 1
    d 19990515 2330 26.34 37 32.17 118 48.18   9.35 D 1.44  8 178  6.  0.04  0.8  2.7 CW   51082742 1
    e 19990515 2330 58.42 37 29.59 118 50.17   8.70 D 1.81 16 143  5.  0.04  0.4  0.6 BW   51082743 1
    f 19990515 2333 03.78 37 30.47 118 49.37   2.43 D 1.05  7 170  5.  0.06  1.0  7.6 CW   51082744 1
    g 19990515 2333 39.33 37 32.29 118 49.76   0.04 D 2.35 15 117  5.  0.91  1.9  3.4 CW   51082745 1
    h 19990515 2334 36.70 37 29.20 118 49.79   4.56 D 3.24 24 121  5.  0.06  0.2  0.6 BW   51082746 1

    Original Triggers:             subnets                        contains
                        duration   tot/coinc                      hypocenters
    19990515 23:25:48.00   33 sec   1 / 1      EVENT ID: 025307
    19990515 23:26:27.00  500 sec   8 / 7      EVENT ID: 025308    a-g
    19990515 23:34:48.00  500 sec  44 /28      EVENT ID: 025309    h+

    Test1: carlsubtrig parameters same as ew04.
    -------------------------------------------
    NetTriggerDur    10   # Nominal duration (sec) for the base network trigger
    SubnetContrib    15   # Addition contribution (sec) to network trigger.
                          #   duration for each concurrently triggered subnet.
    MaxDuration     500   # Maximum duration (sec) allowed for network trigger.

    Test1 produced very similar results to the original! Eight hypocenters and
    three triggers, one short and two long.  The third trigger in Test1 is
    shorter than maximum duration (the original was maximum); this is probably
    because the waveform data ends while the third trigger was still active.
    Remember, the second trigger is really the one we're testing here.

    Test1 Hypocenters:
    a 19990521 1616 08.32 37 32.28 118 48.72   6.47 D 2.12 23 110  6.  0.06  0.3  0.6 BW   40099203 1
    b 19990521 1617 33.71 37 32.17 118 49.14   4.95 D 1.03  8 122  6.  0.02  0.5  1.3 BW   40099204 1
    c 19990521 1619 32.65 37 29.55 118 48.66   6.90 D 1.08  6 195  7.  0.01  1.4  2.7 CW   40099205 1
    d 19990521 1620 10.42 37 32.17 118 48.18   9.35 D 1.44  8 178  6.  0.04  0.8  2.7 CW   40099206 1
    e 19990521 1620 42.49 37 29.57 118 50.18   8.70 D 1.81 16 143  5.  0.04  0.4  0.6 BW   40099207 1
    f 19990521 1622 47.86 37 30.47 118 49.37   2.43 D 1.05  7 170  5.  0.06  1.0  7.6 CW   40099208 1
    g 19990521 1623 23.41 37 32.29 118 49.76   0.04 D 2.35 15 117  5.  0.91  1.9  3.4 CW   40099209 1
    h 19990521 1624 20.78 37 29.19 118 49.79   4.55 D 3.22 24 121  5.  0.05  0.2  0.6 BW   40099210 1

    Test1 Triggers:                   subnets                        contains
                          duration   tot/coinc                      hypocenters
    19990521 16:15:38.00   25 sec     1 /  1      EVENT ID: 027212
    19990521 16:16:12.00  500 sec     8 /  7      EVENT ID: 027213   a-g
    19990521 16:24:33.00  435 sec    41 / 26      EVENT ID: 027214   h+

    Test2: new carlsubtrig parameters (reduce SubnetContrib)
    --------------------------------------------------------
    NetTriggerDur    10   # Nominal duration (sec) for the base network trigger
    SubnetContrib    10   # Addition contribution (sec) to network trigger.
                          #   duration for each concurrently triggered subnet.
    MaxDuration     500   # Maximum duration (sec) allowed for network trigger.

    This will make total duration for triggers shorter by something
    on the order of nSubnet * 5 seconds.  I'm not sure if this will make
    one-subnet triggers too short or not.
    Note: Rick said that he wouldn't like to shorten the duration
          for the one-subnet earthquake, so on to Test3!

    Test2 Hypocenters:
    a 19990521 1818 51.53 37 32.28 118 48.72   6.47 D 2.12 23 110  6.  0.06  0.3  0.6 BW   40099212 1
    b 19990521 1820 16.92 37 32.17 118 49.14   4.95 D 1.03  8 122  6.  0.02  0.5  1.3 BW   40099213 1
    c 19990521 1822 15.86 37 29.55 118 48.66   6.90 D 1.08  6 195  7.  0.01  1.4  2.7 CW   40099214 1
    d 19990521 1822 53.63 37 32.17 118 48.18   9.35 D 1.44  8 178  6.  0.04  0.8  2.7 CW   40099215 1
    e 19990521 1823 25.70 37 29.57 118 50.18   8.70 D 1.81 16 143  5.  0.04  0.4  0.6 BW   40099216 1
    f 19990521 1825 31.07 37 30.47 118 49.37   2.43 D 1.05  7 170  5.  0.06  1.0  7.6 CW   40099217 1
    g 19990521 1826 06.62 37 32.29 118 49.76   0.04 D 2.35 15 117  5.  0.91  1.9  3.4 CW   40099218 1
    h 19990521 1827 03.99 37 29.19 118 49.79   4.55 D 3.22 24 121  5.  0.05  0.2  0.6 BW   40099219 1

    Test2 Triggers:                   subnets                       contains
                          duration   tot/coinc                      hypocenters
    19990521 18:18:21.00   20 sec     1 /  1     EVENT ID: 027215
    19990521 18:18:55.00   84 sec     8 /  7     EVENT ID: 027216     a
    19990521 18:20:20.00   69 sec     5 /  5     EVENT ID: 027217     b
    19990521 18:21:59.00   20 sec     1 /  1     EVENT ID: 027218
    19990521 18:22:20.00  500 sec    42 / 27     EVENT ID: 027219     c-h

    Test3: new carlsubtrig params (increase NetTriggerDur,reduce SubnetContrib)
    ---------------------------------------------------------------------------
    NetTriggerDur    15   # Nominal duration (sec) for the base network trigger
    SubnetContrib    10   # Addition contribution (sec) to network trigger.
                          #   duration for each concurrently triggered subnet.
    MaxDuration     500   # Maximum duration (sec) allowed for network trigger.

    This will make total duration for triggers shorter (than original) by
    something on the order of (nSubnet-1) * 5 seconds.  This will retain the
    original trigger length for one-subnet triggers, but make multiple-subnet
    triggers shorter.

    Test3 Hypocenters:
    a 19990521 1957 31.15 37 32.28 118 48.72   6.47 D 2.12 23 110  6.  0.06  0.3  0.6 BW   40099221 1
    b 19990521 1958 56.54 37 32.17 118 49.14   4.95 D 1.03  8 122  6.  0.02  0.5  1.3 BW   40099222 1
    c 19990521 2000 55.48 37 29.55 118 48.66   6.90 D 1.08  6 195  7.  0.01  1.4  2.7 CW   40099223 1
    d 19990521 2001 33.25 37 32.17 118 48.18   9.35 D 1.44  8 178  6.  0.04  0.8  2.7 CW   40099224 1
    e 19990521 2002 05.32 37 29.57 118 50.18   8.70 D 1.81 16 143  5.  0.04  0.4  0.6 BW   40099225 1
    f 19990521 2004 10.69 37 30.47 118 49.37   2.43 D 1.05  7 170  5.  0.06  1.0  7.6 CW   40099226 1
    g 19990521 2004 46.24 37 32.29 118 49.76   0.04 D 2.35 15 117  5.  0.91  1.9  3.4 CW   40099227 1
    h 19990521 2005 43.61 37 29.19 118 49.79   4.55 D 3.22 24 121  5.  0.05  0.2  0.6 BW   40099228 1

    Test3 Triggers:                   subnets                       contains
                          duration   tot/coinc                      hypocenters
    19990521 19:57:01.00   25 sec     1 /  1    EVENT ID: 027220
    19990521 19:57:34.00  160 sec     9 /  7    EVENT ID: 027221    a,b
    19990521 20:00:34.00  500 sec    44 / 29    EVENT ID: 027222    c-h

    Test4: new carlsubtrig params (increase NetTriggerDur,reduce SubnetContrib)
    ---------------------------------------------------------------------------
    NetTriggerDur    17   # Nominal duration (sec) for the base network trigger
    SubnetContrib     8   # Addition contribution (sec) to network trigger.
                          #   duration for each concurrently triggered subnet.
    MaxDuration     500   # Maximum duration (sec) allowed for network trigger.

    Maintain one-subnet duration!  Decrease multi-subnet durations.

    Test4 Hypocenters:
    a 19990521 2025 55.50 37 32.28 118 48.72   6.47 D 2.12 23 110  6.  0.06  0.3  0.6 BW   40099230 1
    b 19990521 2027 20.89 37 32.17 118 49.14   4.95 D 1.03  8 122  6.  0.02  0.5  1.3 BW   40099231 1
    c 19990521 2029 19.83 37 29.55 118 48.66   6.90 D 1.08  6 195  7.  0.01  1.4  2.7 CW   40099232 1
    d 19990521 2029 57.60 37 32.17 118 48.18   9.35 D 1.44  8 178  6.  0.04  0.8  2.7 CW   40099233 1
    e 19990521 2030 29.67 37 29.57 118 50.18   8.70 D 1.81 16 143  5.  0.04  0.4  0.6 BW   40099234 1
    f 19990521 2032 35.04 37 30.47 118 49.37   2.43 D 1.05  7 170  5.  0.06  1.0  7.6 CW   40099235 1
    g 19990521 2033 10.59 37 32.29 118 49.76   0.04 D 2.35 15 117  5.  0.91  1.9  3.4 CW   40099236 1
    h 19990521 2034 07.96 37 29.20 118 49.79   4.56 D 3.22 24 121  5.  0.06  0.2  0.6 BW   40099237 1

    Test4 Triggers:                   subnets                       contains
                          duration   tot/coinc                      hypocenters
    19990521 20:25:25.00   25 sec     1 /  1    EVENT ID: 027223
    19990521 20:25:59.00  151 sec     8 /  7    EVENT ID: 027224    a,b
    19990521 20:29:03.00  173 sec     8 /  7    EVENT ID: 027225    c,d,e
    19990521 20:32:01.00  404 sec    44 / 27    EVENT ID: 027226    f,g,h+

    I left carlsubtrig running on campbell on real-time data using the parameters
    from test4.  Check its results (after 19990521 2200 UTC) against ew03/ew04.
    5/24  Turns out campbell's system clock was not set to network time.  All
    of the carltrig events are 2 minutes later than the ew03,ew04 events for
    the weekend.  I started ntp on campbell at approx. 19990524 1700 UTC.

    ===========================================================================
    Rick gave me the following other examples of "too long multiples" from ew03:
             tot#/#coinc  time range to save:
    025164     49/29
    025195     88/47
    025268     44/28      19990515 23:33:48.00  590 sec
    025393     32/12      19990516 08:12:11.00  509 sec
    025470     30/19      19990516 12:16:48.00  498 sec
    025777
    026089
    026090
    026169

    ===========================================================================
    Next test case:

    src     id      original trigger    dur     time range saved for testing
    ----   -----  ---------------------------   ----------------------------
    ew03   25393  19990516 08:13:11.00  419 s   start:  19990516 08:12:11.00  509 sec
                                                  end:  19990516 08:20:40.00
                                                 file:  19990516.0812.tbuf

    Original Hypocenters:
    a 19990516 0813 06.24 36  3.84 120 28.10   2.16 D 1.43  6 165 12.  0.03  1.2  7.5 CW   40097977 1
    b 19990516 0815 20.36 37 29.57 118 49.78   5.51 D 3.00 26 125  5.  0.05  0.3  0.6 BW   40097978 1

    Original Trigger:              subnets                        contains
                        duration   tot/coinc                      hypocenters
    19990516 08:13:11.00  419 sec  32 / 12      EVENT ID: 025393  a,b

    Test1: carlsubtrig parameters same as ew03.
    -------------------------------------------
    NetTriggerDur    10   # Nominal duration (sec) for the base network trigger
    SubnetContrib    15   # Addition contribution (sec) to network trigger.
                          #   duration for each concurrently triggered subnet.
    MaxDuration     500   # Maximum duration (sec) allowed for network trigger.

    Let's first make sure we duplicate ew03's original results!
    Well, it wasn't quite as long, but it was indeed one big trigger
    containing both events.

    Test1 Hypocenters:
    a 19990524 1625 31.02 36  3.84 120 28.10   2.16 D 1.43  6 165 12.  0.03  1.2  7.5 CW   40099437 1
    b 19990524 1627 45.14 37 29.57 118 49.78   5.52 D 3.00 26 125  5.  0.05  0.3  0.6 BW   40099438 1

    Test1 Triggers:                   subnets                        contains
                          duration   tot/coinc                      hypocenters
    19990524 16:25:35.00  392 sec    31 / 12      EVENT ID: 027722  a,b

    Test4: new carlsubtrig params (increase NetTriggerDur,reduce SubnetContrib)
    ---------------------------------------------------------------------------
    NetTriggerDur    17   # Nominal duration (sec) for the base network trigger
    SubnetContrib     8   # Addition contribution (sec) to network trigger.
                          #   duration for each concurrently triggered subnet.
    MaxDuration     500   # Maximum duration (sec) allowed for network trigger.

    Maintain one-subnet duration!  Decrease multi-subnet durations.

    Test4 Hypocenters:
    a 19990524 1702 22.92 36  3.84 120 28.10   2.16 D 1.43  6 165 12.  0.03  1.2  7.5 CW   40099440 1
    b 19990524 1704 37.04 37 29.57 118 49.78   5.52 D 3.00 26 125  5.  0.05  0.3  0.6 BW   40099441 1

    Test4 Triggers:                   subnets                        contains
                          duration   tot/coinc                      hypocenters
    19990524 17:02:27.00   75 sec     6 /  6      EVENT ID: 027723  a
    19990524 17:03:58.00  293 sec    26 / 16      EVENT ID: 027724  ?,b

    Looks like we triggered 40 sec early for the bigger event.  The initial
    duration for the 2nd trigger was 25 sec (1 subnet); there was probabaly
    a little eq that binder didn't get.  At least it split the 1.4 and 3.0
    into two different triggers.

    ===========================================================================
    Next test case:

    src     id      original trigger    dur     time range saved for testing
    ----   -----  ---------------------------   ----------------------------
    ew03   25470  19990516 08:13:11.00  419 s   start:  19990516 12:16:48.00  498 sec
                                                  end:  19990516 12:25:06.00
                                                 file:  19990516.1216.tbuf

    Original Hypocenters:
    a 19990516 1217 44.19 37 31.58 118 48.85   4.09 D 0.81  6 137  6.  0.02  0.4  2.6 CW   40098078 1
    b 19990516 1218 09.54 37 31.11 118 49.49   4.29 D 2.99 25 116  5.  0.06  0.2  0.6 BW   40098079 1

    Original Trigger:              subnets                        contains
                        duration   tot/coinc                      hypocenters
    19990516 12:17:48.00  408 sec  30 / 19     EVENT ID: 025470   a,b

    I don't think there's any way I can separate these two because the larger
    event occurs within the minimum one-subnet event duration.  This will
    cause the TotalNetDuration to increase due to case 2 (in the introductory
    discussion).  Let's just try Test4 anyway (but I expect it to fail).

    Test4: new carlsubtrig params (increase NetTriggerDur,reduce SubnetContrib)
    ---------------------------------------------------------------------------
    NetTriggerDur    17   # Nominal duration (sec) for the base network trigger
    SubnetContrib     8   # Addition contribution (sec) to network trigger.
                          #   duration for each concurrently triggered subnet.
    MaxDuration     500   # Maximum duration (sec) allowed for network trigger.

    Maintain one-subnet duration!  Decrease multi-subnet durations.

    Test4 Hypocenters:
    a 19990524 1723 17.06 37 31.58 118 48.85   4.09 D 0.81  6 137  6.  0.02  0.4  2.6 CW   40099443 1
    b 19990524 1723 42.41 37 31.11 118 49.49   4.30 D 2.99 25 116  5.  0.06  0.2  0.6 BW   40099444 1

    Test4 Triggers:                   subnets                        contains
                          duration   tot/coinc                      hypocenters
    19990524 17:23:20.00  303 sec    29 / 20      EVENT ID: 027725  a,b
    19990524 17:28:35.00   66 sec     5 /  5      EVENT ID: 027726  ?

    Just as I thought, I couldn't seperate events a & b into individual triggers.
    However, we did generate a second trigger in Test4 that in the original
    ew03 run was probably buried at the end of the 408 second trigger.  Binder
    didn't stack this event, but it may have been fairly large, as it
    triggered 5 subnets!  So, there was indeed an improvement.

    ===========================================================================

### Download in other formats:

- Plain Text