Changeset 400 for trunk/liveplot

Timestamp:

Dec 4, 2009 5:05:32 PM (14 years ago)

Author:

toby

Message:

# on 2001/06/29 18:10:45, toby did:
Major revision: add bkgedit functionality

File:

• trunk/liveplot (modified) (11 diffs, 3 props)

trunk/liveplot

@@ -12 +12 @@
     destroy .
 }
+set program [file tail $argv0]
+#set program bkgedit
+
 if {[lindex $argv 1] == ""} {
     set hst 1
@@ -50 +53 @@
 set expgui(autotick) 0
 set expgui(pixelregion) 5
+# location for web pages, if not found locally
+set expgui(website) www.ncnr.nist.gov/xtal/software/expgui
 set peakinfo(obssym) scross
 set peakinfo(obssize) 1.0
@@ -201 +206 @@
     }
     $box element show [lsort -decreasing [$box element show]]
+    global program
+    if {$program == "bkgedit"}  bkghstInit
 }
    
@@ -217 +224 @@
 #       set input [open "| $expgui(gsasexe)/hstdump $expnam  < histdump.inp" w+]
 ###########################################################################
-    # use histdump for right now
+    # use histdmp for histogram info
     set input [open histdump$hst.inp w]
     puts $input "$expnam"
@@ -843 +850 @@
     }
 }
+# evaluate the Chebyshev polynomial with coefficients A at point x
+# coordinates are rescaled from $xmin=-1 to $xmax=1
+proc chebeval {A x xmin xmax} {
+    set xs [expr {-1 + 2 * (1.*$x - $xmin) / (1.*$xmax - 1.*$xmin)}]
+    set Tpp 0
+    set Tp 0
+    set total 0
+    foreach a $A {
+        if {$Tpp == $Tp && $Tp == 0} {
+            set T 1
+        } elseif {$Tpp == 0} {
+            set T $xs
+        } else {        
+            set T [expr {2. * $xs * $Tp - $Tpp}]
+        }
+        set total [expr {$total + $a * $T}]
+        set Tpp $Tp
+        set Tp $T
+    }
+    return $total
+}
+
+# determine a very approximate set of Chebyshev coefficients of order n 
+# to compute Y from X (fast but not very good)
+proc chebgen {X Y xmin xmax n} {
+    if {[llength $X] != [llength $Y]} return
+    set xnorm [expr {2. / ($xmax - $xmin)}]
+    set pi [expr {2*asin(1)}]
+    set a(0) 0.
+    for {set i 1} {$i < $n} {incr i} {set a($i) 0.} 
+    foreach x1 $X x2 [lrange $X 1 end] y1 $Y y2 [lrange $Y 1 end] {
+        if {$x2 == ""} break
+        set xs1 [expr {-1 + ($x1 - $xmin) * $xnorm}]
+        set th1 [expr {acos(-1 + ($x1 - $xmin) * $xnorm)}]
+        set xs2 [expr {-1 + ($x2 - $xmin) * $xnorm}]
+        set th2 [expr {acos(-1 + ($x2 - $xmin) * $xnorm)}]
+        set thbar [expr {($th1 + $th2)/2.}]
+        set dth [expr {$th1 - $th2}]
+        set xsbar [expr {cos($thbar)}]
+        set ybar [expr {($xsbar - $xs1) / ($xs2 - $xs1) * ($y2 - $y1) + $y1}]
+        # seems to work better starting with just 2 terms
+#       for {set i 0} {$i < $n} {incr i} 
+        for {set i 0} {$i < 2} {incr i} {
+            set a($i) [expr {$a($i) + $ybar * cos($i*$thbar) * $dth}]
+        }
+    }
+    set A [expr {$a(0) / $pi}]
+    for {set i 1} {$i < $n} {incr i} {
+        lappend A [expr {2 * $a($i) / $pi}]
+    } 
+    return $A
+}
+
+# disable the improve fit button
+proc bkgResetFit {} {
+    .bkg.f.fit2 config -state disabled
+}
+
+# enable the improve fit button, if appropriate
+proc bkgMoreFit {} {
+    global cheblist
+    if {[llength $cheblist] < 2} {bkgResetFit;return}
+    .bkg.f.fit2 config -state normal
+}
+
+# perform a Gauss-Newton fit to optimize Chebyshev coefficients A
+proc chebGN {X Y S A xmin xmax "damp 0.75"} {
+    # Gauss-Newton
+    if {[llength $X] != [llength $Y]} return
+    set xnorm [expr {2. / ($xmax - $xmin)}]
+    # denominator
+    set sum2 0.
+    foreach x $X s $S {
+        set xs [expr {-1 + (1.*$x - $xmin) * $xnorm}]
+        set Tpp 0
+        set Tp 0
+        foreach a1 $A {
+            if {$Tpp == $Tp && $Tp == 0} {
+                set T 1
+            } elseif {$Tpp == 0} {
+                set T $xs
+            } else {    
+                set T [expr {2. * $xs * $Tp - $Tpp}]
+            }
+            set sum2 [expr {$sum2 + $T /($s*$s)}]
+            set Tpp $Tp
+            set Tp $T
+        }
+    }
+    # Evaluate Ycalc & sum(delta2)
+    set sumd2 0.
+    foreach x $X y $Y {
+#       set xs [expr {-1 + (1.*$x - $xmin) * $xnorm}]
+        lappend Ycalc [set yc [chebeval $A $x $xmin $xmax]]
+        set sumd2 [expr {$sumd2 + ($y - $yc)*($y - $yc)}]
+    }
+    set k -1
+    foreach a $A {incr k; set sum($k) 0}
+    foreach x $X y $Y yc $Ycalc s $S {
+        set xs [expr {-1 + (1.*$x - $xmin) * $xnorm}]
+        set Tpp 0
+        set Tp 0
+        set k -1
+        foreach a $A {
+            incr k
+            if {$Tpp == $Tp && $Tp == 0} {
+                set T 1
+            } elseif {$Tpp == 0} {
+                set T $xs
+            } else {    
+                set T [expr {2. * $xs * $Tp - $Tpp}]
+            }
+            set sum($k) [expr {$sum($k) + ($T * ($yc - $y))/($s*$s)}]
+            set Tpp $Tp
+            set Tp $T
+        }
+    }
+    set sumd2r $sumd2
+    set d $damp
+    # compute new cheb terms
+    while {$d > $damp/32} {
+        set k -1
+        set Anew {}
+        foreach a $A {
+            incr k
+            lappend Anew [expr {$a - $d*($sum($k) / $sum2)}]
+        }
+        # Evaluate new Ycalc & sum(delta2)
+        set sumd2r 0.
+        foreach x $X y $Y {
+#           set xs [expr {-1 + (1.*$x - $xmin) * $xnorm}]
+            set yc [chebeval $Anew $x $xmin $xmax]
+            set sumd2r [expr {$sumd2r + ($y - $yc)*($y - $yc)}]
+        }
+        # are these shifts an improvement?
+        if {$sumd2r < $sumd2} {
+            # are we converged?
+            if {[expr {($sumd2-$sumd2r)/$sumd2}] < 0.0001} {return ""}
+            return $Anew
+        }
+        set d [expr {$d/2.}]
+    }
+    return ""
+}
+
+# change the binding of the mouse, based on the selected mode
+proc bkgEditMode {b} {
+    global zoomcommand box
+    # get binding
+    set bindtag $box
+    catch {
+        if {[bind bltZoomGraph] != ""} {
+            set bindtag bltZoomGraph
+        }
+    }
+    # save the zoom command
+    if [catch {set zoomcommand}] {
+        set zoomcommand [bind $bindtag <1>]
+        .bkg.f.fit1 config -state disabled
+        .bkg.f.fit2 config -state disabled
+        .bkg.f.terms config -state disabled
+    }
+    foreach c {1 2 3} {
+        if {$c == $b} {
+            .bkg.l.b$c config -relief sunken
+        } else {
+            .bkg.l.b$c config -relief raised
+        }
+    }
+    if {$b == 2} {
+        bind $bindtag <1> "bkgAddPoint %x %y"
+        .g config -cursor arrow
+    } elseif {$b == 3} {
+        bind $bindtag <1> "bkgDelPoint %x %y"
+        .g config -cursor circle
+    } else {
+        bind $bindtag <1> $zoomcommand
+        .g config -cursor crosshair
+    }
+}
+
+# plot the background points
+proc bkgPointPlot {} {
+    global bkglist termmenu chebterms expnam hst tmin tmax
+    set l {}
+    set fp [open $expnam.bkg$hst w]
+    puts $fp "y p h e $hst b ! fixed background points for use in BKGEDIT"
+    foreach p $bkglist {
+        puts $fp "i\t$p\t0.0"
+        append l " $p"
+    }
+    if {[llength $bkglist] > 0} {
+        puts $fp "i\t[expr $tmin*0.99] [lindex [lindex $bkglist 0] 1]\t0.0"
+        puts $fp "i\t[expr $tmax*1.01] [lindex [lindex $bkglist end] 1]\t0.0"
+    }
+    close $fp
+    .g element config 12 -data $l
+    if {[set l [llength $bkglist]] > 3} {
+        .bkg.f.fit1 config -state normal
+        .bkg.f.terms config -state normal
+        $termmenu delete 0 end
+        set imax {}
+        for {set i 2} {$i <= $l/1.5} {incr i 2} {
+            $termmenu insert end radiobutton -label $i \
+                    -variable chebterms  -command {bkgMoreFit}
+            set imax $i
+        }
+        if {$imax < $chebterms} {set chebterms $imax}
+    } else {
+        .bkg.f.fit1 config -state disabled
+        .bkg.f.fit2 config -state disabled
+        .bkg.f.terms config -state disabled
+        set chebterms 2
+    }
+}
+
+# add a bkg point at screen coordinates x,y
+proc bkgAddPoint {x y} {
+    global bkglist tmin tmax
+    set xy [.g invtransform $x $y]
+    set x [lindex $xy 0]
+    if {$x < $tmin} {set x $tmin}
+    if {$x > $tmax} {set x $tmax}
+    lappend bkglist [list $x [lindex $xy 1]]
+    set bkglist [lsort -real -index 0  $bkglist]
+    bkgMoreFit
+    bkgFillPoints
+    bkgPointPlot
+}
+
+# delete the bkg point closest to screen coordinates x,y
+proc bkgDelPoint {x y} {
+    global bkglist
+    set closest {}
+    set dist2 {}
+    set i -1
+    foreach p $bkglist {
+        incr i
+        set sxy [eval .g transform $p]
+        if {$closest == ""} {
+            set closest $i
+            set dist2 0
+            foreach v1 $sxy v2 "$x $y" {
+                set dist2 [expr {$dist2 + ($v1 - $v2)*($v1 - $v2)}]
+            }
+        } else {
+            set d2 0
+            foreach v1 $sxy v2 "$x $y" {
+                set d2 [expr {$d2 + ($v1 - $v2)*($v1 - $v2)}]
+            }
+            if {$d2 < $dist2} {
+                set closest $i
+                set dist2 $d2
+            }           
+        }
+    }
+    set bkglist [lreplace $bkglist $closest $closest]
+    bkgMoreFit
+    bkgPointPlot
+    bkgFillPoints
+}
+
+# initialize the background plot
+proc bkghstInit {} {
+    global bkglist tmin tmax hst expnam cheblist chebterms
+    set tmin [histinfo $hst tmin]
+    set tmax [histinfo $hst tmax]
+    if {[catch {expr $tmin}] || [catch {expr $tmax}]} {
+        tk_dialog .err "MIN/MAX Error" "Error -- Unable read tmin or tmax (has POWPREF been run?" \
+                error 0 Quit
+        destroy .
+    }
+
+    set bkglist {}
+    if [file exists $expnam.bkg$hst] {
+        catch {
+            set fp [open $expnam.bkg$hst r]
+            gets $fp line
+            while {[gets $fp line]>=0} {
+                set x [lindex $line 1]
+                set y [lindex $line 2]
+                if {$x >= $tmin && $x <= $tmax} {
+                    lappend bkglist [list $x $y]
+                }
+            }
+        }
+        close $fp
+    }
+
+    bkgEditMode 1
+    bkgPointPlot
+    bkgFillPoints
+    set cheblist ""
+    bkgResetFit
+    BkgFillCheb
+    set chebterms 2
+}
+
+# fit a Chebyshev polynomial to the selected background points
+proc bkgFit {termlist button} {
+    global bkglist chebterms cheblist
+    $button config -relief sunken
+    update
+    foreach p $bkglist {
+        lappend S 1.
+        foreach v $p var {X Y} {
+            lappend $var $v
+        }
+    }
+    global tmin tmax
+    if {[llength $termlist] < 2} {
+        # get a starting point
+        set termlist [chebgen $X $Y $tmin $tmax $chebterms]
+        # plot it
+        set calcb {}
+        foreach x [xvec range 0 end] {
+            lappend calcb [chebeval $termlist $x $tmin $tmax]
+        }
+        .g element configure 11 -xdata xvec -ydata $calcb
+        update
+    } elseif {[llength $termlist] < $chebterms} {
+        while {[llength $termlist] < $chebterms} {
+            lappend termlist 0.
+        }
+    } elseif {[llength $termlist] > $chebterms} {
+        set termlist [lrange $termlist 0 [expr $chebterms -1]]
+    }
+    # iterate 
+    for {set i 1} {$i < 20} {incr i} {
+        set termlist1 [chebGN $X $Y $S $termlist $tmin $tmax]
+        # have we converged?
+        if {$termlist1 == ""} {
+            bkgResetFit
+            set cheblist $termlist
+            BkgFillCheb
+            bkgFillPoints
+            $button config -relief raised
+            return 
+        }
+        set termlist $termlist1
+        set calcb {}
+        foreach x [xvec range 0 end] {
+            lappend calcb [chebeval $termlist $x $tmin $tmax]
+        }
+        .g element configure 11 -xdata xvec -ydata $calcb
+        update
+    }
+    set cheblist $termlist
+    BkgFillCheb
+    bkgFillPoints
+    bkgMoreFit
+    $button config -relief raised
+}
+
+# put the Chebyshev coefficients into edit widgets
+proc BkgFillCheb {} {
+    global cheblist
+    global chebedit
+    catch {destroy .bkg.canvas.fr}
+    set top [frame .bkg.canvas.fr]
+    .bkg.canvas create window 0 0 -anchor nw -window $top
+    # delete trace on chebedit
+    foreach v [ trace vinfo chebedit] {
+        eval trace vdelete chebedit $v
+    }
+    if {[llength $cheblist] == 0} {
+        grid [label $top.0 -text "(no terms defined)"] -col 1 -row 1
+        .bkg.cw config -state disabled
+    } else {
+        set i -1
+        .bkg.cw config -state normal
+        foreach c $cheblist {
+            incr i
+            grid [frame $top.$i -relief groove -bd 3] -col $i -row 1
+            grid [label $top.$i.l -text "[expr 1+$i]"] -col 1 -row 1
+            grid [entry $top.$i.e -textvariable chebedit($i) -width 13] \
+                    -col 2 -row 1
+            set chebedit($i) $c
+        }
+        trace variable chebedit w "BkgRecalcCheb $top" 
+    }
+    update idletasks
+    set sizes [grid bbox $top]
+    .bkg.canvas config -scrollregion $sizes -height [lindex $sizes 3]
+}
+
+# respond to edits made to Chebyshev terms
+proc BkgRecalcCheb {top var i mode} {
+    global chebedit cheblist
+    if [catch {expr $chebedit($i)}] {
+        $top.$i.e config -fg red
+    } else {
+        $top.$i.e config -fg black
+        set cheblist [lreplace $cheblist $i $i $chebedit($i)]
+        global tmin tmax
+        # plot it
+        set calcb {}
+        foreach x [xvec range 0 end] {
+            lappend calcb [chebeval $cheblist $x $tmin $tmax]
+        }
+        .g element configure 11 -xdata xvec -ydata $calcb
+        update
+        bkgMoreFit
+    }
+}
+
+# put the bkg points into edit widgets
+proc bkgFillPoints {} {
+    global bkglist tmin tmax bkgedit
+    # delete trace on bkgedit
+    foreach v [ trace vinfo bkgedit] {
+        eval trace vdelete bkgedit $v
+    }
+    catch {destroy .bkg.bc.fr}
+    set top [frame .bkg.bc.fr]
+    .bkg.bc create window 0 0 -anchor nw -window $top
+    if {[llength $bkglist] == 0} {
+        grid [label $top.0 -text "(no points defined)"] -col 1 -row 1
+    } else {
+        set i -1
+        foreach p $bkglist {
+            incr i
+            grid [frame $top.$i -relief groove -bd 3] -col $i -row 1
+            grid [label $top.$i.l -text "[expr 1+$i]"] -col 1 -rowspan 2 -row 1
+            grid [entry $top.$i.ex -textvariable bkgedit(x$i) -width 13] \
+                    -col 2 -row 1
+            grid [entry $top.$i.ey -textvariable bkgedit(y$i) -width 13] \
+                    -col 2 -row 2
+            foreach val $p var {x y} {
+                set bkgedit(${var}$i) $val
+            }
+        }
+        trace variable bkgedit w "BkgRecalcBkg $top" 
+    }
+    update idletasks
+    set sizes [grid bbox $top]
+    .bkg.bc config -scrollregion $sizes -height [lindex $sizes 3]
+}
+
+# respond to edits made to bkg points
+proc BkgRecalcBkg {top var i mode} {
+    global bkgedit bkglist tmin tmax
+    regexp {(.)([0-9]*)} $i junk var num
+    if [catch {expr $bkgedit($i)}] {
+        $top.$num.e$var config -fg red
+    } else {
+        $top.$num.e$var config -fg black
+        set p [lindex $bkglist $num]
+        if {$var == "x"} {
+            set x $bkgedit($i)
+            if {$x < $tmin} {set x $tmin}
+            if {$x > $tmax} {set x $tmax}
+            set bkglist [lreplace $bkglist $num $num \
+                    [list $x [lindex $p 1]]]
+        } else {
+            set bkglist [lreplace $bkglist $num $num \
+                    [list [lindex $p 0] $bkgedit($i)]]
+        }
+    }
+        bkgPointPlot
+}
+
+# save the Chebyshev terms in the .EXP file
+proc bkgChebSave {} {
+    global hst cheblist expgui Revision expmap expnam
+    histinfo $hst backtype set 1
+    histinfo $hst backterms set [llength $cheblist]
+    set num 0
+    foreach v $cheblist {
+        set var "bterm[incr num]"
+        histinfo $hst $var set $v
+    }
+    histinfo $hst bref set 0
+    # add a history record
+    exphistory add " BKGEDIT [lindex $Revision 1] [lindex $expmap(Revision) 1] -- [clock format [clock seconds]]"
+    # now save the file
+    expwrite $expnam.EXP
+}
+
 
 source [file join $expgui(scriptdir) gsascmds.tcl]
+source [file join $expgui(scriptdir) readexp.tcl]
 
 # override options with locally defined values
@@ -887 +1374 @@
             error 0 "Limp ahead"
 }
+# modify zoom so that y2axis is not zoomed in for blt2.4u+
+catch {
+    regsub -all y2axis [info body blt::PushZoom] " " b1
+    proc blt::PushZoom {graph} $b1
+}
+
 $box element create 0 -xdata xvec -ydata wifdvec -color magenta \
         -line 3 -symbol none -label "Chi2" -mapy y2
+$box element create 1 -label bckgr -color green  -symbol none  
+$box element config 1 -xdata xvec -ydata bckvec
 $box element create 3 -color black -linewidth 0 -label Obs \
         -symbol $peakinfo(obssym) \
@@ -894 +1389 @@
 $box element create 2 -label Calc -color red  -symbol none  
 $box element create 4 -label diff -color blue  -symbol none  
+
+if {$program == "liveplot"} {
+    $box y2axis config -min 0 -title {Cumulative Chi Squared}
+} elseif {$program == "bkgedit"}  {
+    eval $box element config 0 $graph(ElementHideOption) 
+    eval $box y2axis config $graph(ElementHideOption) 
+    $box element config 0 -label ""
+    eval $box element config 1 $graph(ElementHideOption) 
+    $box element config 1 -label ""
+    eval $box element config 4 $graph(ElementHideOption) 
+    $box element config 4 -label ""
+    $box element create 11
+    $box element create 12
+    $box element configure 12  -color magenta  -pixels 12 \
+            -line 0 -symbol triangle -label "bkg pts"
+    $box element configure 11 -color blue \
+            -symbol none -label "Cheb fit" -dashes 5 -line 2
+    $box element show "3 2 11 12"
+}
 $box element config 3 -xdata xvec -ydata obsvec
 $box element config 2 -xdata xvec -ydata calcvec
 $box element config 4 -xdata xvec -ydata diffvec
+
 if {$expgui(tcldump) != ""} {
-    $box element create 1 -label bckgr -color green  -symbol none  
-    $box element config 1 -xdata xvec -ydata bckvec
-    bind . <Shift-Button-1> "lblhkl $box %x"
     bind . <Key-h> "lblhkl $box %x"
     bind . <Key-H> "lblhkl $box %x"
 #    bind $box <Shift-Double-Button-1> "lblallhkl %W"
-    bind $box <Shift-Button-3> "delallhkllbl %W"
+    if {[bind bltZoomGraph] != ""} {
+        bind bltZoomGraph <Shift-Button-1> "lblhkl $box %x"
+        bind bltZoomGraph <Shift-Button-3> "delallhkllbl %W"
+    } else {
+        bind $box <Shift-Button-1> "lblhkl $box %x"
+        bind $box <Shift-Button-3> "delallhkllbl %W"
+    }
+} else {
+    $box element config 1 -label ""
+    eval $box element config 4 $graph(ElementHideOption) 
 }
 $box yaxis config -title {} 
 setlegend $box $graph(legend)
-ShowCumulativeChi2
-
-updateifnew
+
 frame .a -bd 3 -relief groove
 pack [menubutton .a.file -text File -underline 0 -menu .a.file.menu] -side left
@@ -953 +1472 @@
 }
 .a.options.menu add command -label "Obs symbol" -command setsymopts
-if {$expgui(tcldump) != ""} {
+if {$expgui(tcldump) != "" && $program == "liveplot"} {
     .a.options.menu add cascade -label "X units" -menu .a.options.menu.xunits
     menu .a.options.menu.xunits
@@ -977 +1496 @@
             -variable graph(backsub) \
             -command {set cycle [getcycle];readdata .g}
+} else {
+    set graph(xunits) 0
 }
     
@@ -983 +1504 @@
         -command {setlegend $box $graph(legend)}
 .a.options.menu add command -label "Set PS output" -command setpostscriptout
-.a.options.menu add checkbutton -label "Raise on update" \
-        -variable graph(autoraise)
-.a.options.menu add checkbutton -label "Cumulative Chi2" \
-        -variable graph(chi2) -command ShowCumulativeChi2
-.a.options.menu add command -label "Save Options" -underline 1 \
-        -command "SaveOptions"
-$box y2axis config -min 0 -title {Cumulative Chi Squared}
-ShowCumulativeChi2
+if {$program == "liveplot"} {
+    .a.options.menu add checkbutton -label "Raise on update" \
+            -variable graph(autoraise)
+    .a.options.menu add checkbutton -label "Cumulative Chi2" \
+            -variable graph(chi2) -command ShowCumulativeChi2
+    .a.options.menu add command -label "Save Options" -underline 1 \
+            -command "SaveOptions"
+    ShowCumulativeChi2
+} elseif {$program == "bkgedit"}  {
+    catch {pack [frame .bkg -bd 3 -relief sunken] -side bottom -fill both}
+    grid [label .bkg.top -text "Background Point Editing"] \
+            -col 0 -row 0 -columnspan 4
+    grid [button .bkg.help -text Help -bg yellow \
+            -command "MakeWWWHelp liveplot.html bkgedit"] \
+            -column 5 -row 0 -rowspan 2 -sticky n
+    
+    grid [frame .bkg.l -bd 3 -relief groove] \
+            -col 0 -row 1 -columnspan 2 -sticky nse
+    grid [label .bkg.l.1 -text "Mouse click\naction"] -col 0 -row 0
+    foreach c {1 2 3} l {zoom add delete} {
+        grid [button .bkg.l.b$c -text $l -command "bkgEditMode $c"] \
+                -col $c -row 0
+    }
+    grid [frame .bkg.f -bd 3 -relief groove] \
+            -col 3 -row 1 -columnspan 2 -sticky nsw
+    grid [button .bkg.f.fit1 -text "Start\nFit" -command {bkgFit "" .bkg.f.fit1}] \
+            -col 1 -row 1
+    grid [button .bkg.f.fit2 -text "Improve\nFit" \
+            -command {bkgFit $cheblist .bkg.f.fit2}] -col 2 -row 1
+    grid [label .bkg.f.tl -text "with"] -col 3 -row 1
+    set termmenu [tk_optionMenu .bkg.f.terms chebterms 0]
+    grid .bkg.f.terms -col 4 -row 1
+    grid [label .bkg.f.tl1 -text "terms"] -col 5 -row 1
+
+    grid [frame .bkg.c1 -bd 3 -relief groove] \
+            -col 0 -row 5 -rowspan 2 -sticky nsew
+    grid [label .bkg.c1.1 -text "Chebyshev\nterms"] -col 0 -row 0
+    grid [canvas .bkg.canvas \
+            -scrollregion {0 0 5000 500} -width 0 -height 0 \
+            -xscrollcommand ".bkg.scroll set"] \
+            -column 1 -row 5 -columnspan 3 -sticky nsew
+    grid [scrollbar .bkg.scroll -command ".bkg.canvas xview" \
+            -orient horizontal] -column 1 -row 6 -columnspan 3 -sticky nsew
+    grid [button .bkg.cw -text "Save in EXP\nfile & Exit" \
+            -command "bkgChebSave;exit"] \
+            -col 4 -columnspan 2 -row 5 -rowspan 2 -sticky ns
+
+    grid [frame .bkg.bl -bd 3 -relief groove] \
+            -col 0 -row 3 -rowspan 2 -sticky nsew
+    grid [label .bkg.bl.1 -text "Background\npoints"] -col 0 -row 0
+    grid [canvas .bkg.bc \
+            -scrollregion {0 0 5000 500} -width 0 -height 0 \
+            -xscrollcommand ".bkg.bs set"] \
+            -column 1 -row 3 -columnspan 5 -sticky nsew
+    grid [scrollbar .bkg.bs -command ".bkg.bc xview" -orient horizontal] \
+            -column 1 -row 4 -columnspan 5 -sticky nsew
+
+    grid columnconfigure .bkg 1 -weight 1
+    grid columnconfigure .bkg 2 -weight 1
+    grid columnconfigure .bkg 3 -weight 1
+    grid rowconfigure .bkg 3 -weight 1
+    grid rowconfigure .bkg 5 -weight 1
+    .g config -title ""
+}
 
 pack [menubutton .a.help -text Help -underline 0 -menu .a.help.menu] -side right
@@ -1006 +1583 @@
 if [file exists $fl] {source $fl}
 
+if {$program == "bkgedit"}  {
+    expload $expnam.EXP
+    mapexp
+
+#    bkghstInit
+}
+updateifnew
 donewaitmsg
 trace variable peakinfo w plotdataupdate