Light Ghosting

uniform sampler2D bgl_RenderedTexture;
int samples = 8;
float gain = 0.5;
float spread = 0.03;
float  threshold = 1.0;
int volumes = 4;

void main(void) {
    vec3 screen = texture2D(bgl_RenderedTexture, gl_TexCoord[0].st).rgb;
    vec2 deltacoord = (gl_TexCoord[0].st - vec2(0.5))*0.5;
 
    vec2 texcoord = -gl_TexCoord[0].st+1.0;
    float lum = 0.0;
 
    for (int i = 0; i < volumes; ++i) {      
        texcoord += deltacoord;
        for (int j = -samples; j < samples; ++j) {      
            vec2 newcoord = vec2(texcoord.x+(j*0.005), texcoord.y);
            float col = texture2D(bgl_RenderedTexture, newcoord).r;
            if (col >= threshold) {
                lum += spread/samples;
                } else {
                    lum += 0.0;
                    };
                };                                
        };
    float white = dot(screen, vec3(0.299, 0.587, 0.114)).r;
    float cval = ((0.5 - gl_TexCoord[0].x)*(0.5 - gl_TexCoord[0].x)) + ((0.5 - gl_TexCoord[0].y)*(0.5 - gl_TexCoord[0].y));
    float root = sqrt(cval);
 
    float stencil = max(lum , root/6);
    stencil = stencil-(root/6);
 
    gl_FragColor = vec4(screen+vec3(stencil*gain, stencil*gain, stencil*(gain*2)), 1.0);  
 
    }