DataPoint.java

Go to the documentation of this file.

package cedar.jetweb.model.plots;

import cedar.jetweb.JetWebException;

public class DataPoint implements Cloneable, Comparable {

    private double x_;// = null;
    private double y_;// = null;
    private double yUp_;// = null;
    private double yDown_;// = null;
    private double binWidth_;// = null;
    private double chi2_;// = null;
    private Integer pointNumber_;

    public DataPoint(){
    return;
    }

    public int compareTo(Object dpo)
    throws ClassCastException{
    
    if(! (dpo instanceof DataPoint)) {
        throw new ClassCastException
        ("DataPoint.compareTo: cannot compare object " + dpo + 
         " with DataPoint "+this.toString());
    }

    DataPoint dp = (DataPoint) dpo;

    if(pointNumber_ < dp.getNumber()) return -1;
    if(pointNumber_ > dp.getNumber()) return 1;

    if(x_<dp.getX()) return -1;
    if(x_>dp.getX()) return 1;
    return 0;
    } 

    public double[] toArray(){
    double[] array = {x_, y_, yUp_, yDown_, binWidth_, chi2_};
    return array;
    }
    public double getX(){
    return x_;
    }
    public double getY(){
    return y_;
    }
    public double getYUp(){
    return yUp_;
    }
    public double getYDown(){
    return yDown_;
    }
    public double getBinWidth(){
    return binWidth_;
    }
    public double getChi2(){
    return chi2_;
    }
    public DataPoint setX(double d){
        x_ = d;
    return this;
    }
    public DataPoint setY(double d){
    y_ = d;
    return this;
    }
    public DataPoint setYUp(double d){
    yUp_ = d;
    return this;
    }
    public DataPoint  setYDown(double d){
    yDown_ = d;
    return this;
    }
    public DataPoint  setBinWidth(double d){
    binWidth_ = d;
    return this;
    }

    public Integer getNumber(){
    return pointNumber_;
    }

    public DataPoint setNumber(Integer number){
    pointNumber_ = number;
    return this;
    }

    public DataPoint setChi2(double d){
    chi2_ = d;
    return this;
    }


    

    public Object clone(){

    DataPoint newPoint=new DataPoint();
    newPoint.setX(getX());
    newPoint.setY(getY());
    newPoint.setYUp(getYUp());
    newPoint.setYDown(getYDown());
    newPoint.setBinWidth(getBinWidth());
    newPoint.setChi2(getChi2());

    return newPoint;
    }
    public double getXMax(){
    return getX() + 0.5*getBinWidth();
    }
    public double getXMin(){
    return getX() - 0.5*getBinWidth();
    }
    public String toString(){
    StringBuffer b = new StringBuffer("DataPoint");
    b.append("\nX:");b.append(x_);
    b.append("\nY:");b.append(y_);
    b.append("\nY ERROR UP:");b.append(yUp_);
    b.append("\nY ERROR DOWN:");b.append(yDown_);
    b.append("\nBINWIDTH:");b.append(binWidth_);
    b.append("\nCHI2:");b.append(chi2_);
    return b.toString();
    }

    public void add(DataPoint p2){

    if (Math.abs(getX()-p2.getX())/(getX()+p2.getX())>0.001 ){
        System.out.println("x values not equal "+getX()+","+p2.getX());
    } else if (p2.getYUp()==0.0 && p2.getYDown()==0.0) {
    // No errors on new data - do nothing
    } else if (getYUp()==0.0 && getYDown()==0.0) {
        // No data in this point. Replace it by the new one
        setX(p2.getX());
        setY(p2.getY());
        setYUp(p2.getYUp());
        setYDown(p2.getYDown());

    } else {
        double sig1 = 0.5*(getYUp()+getYDown());
        //      double sig1sq = Math.pow( (getYUp()+getYDown())/2.0, 2);
        double sig1sq = sig1*sig1;
        double sig2 = 0.5 *(p2.getYUp()+p2.getYDown());
        double sig2sq = sig2 * sig2;;
        
        double w1 = 1.0/sig1sq;
        double w2 = 1.0/sig2sq;

        double invW1W2 = 1.0/(w1+w2);

        // Contents.
        double tmp1 = (getY()*w1+p2.getY()*w2) * invW1W2;
        // Errors
        //setYUp(1.0/Math.sqrt(w1+w2));
        setYUp(Math.sqrt(invW1W2));
        setYDown(getYUp());
        setY(tmp1);
    }                

    }
    public void fixZeroError(double lumi){

    if (lumi<=0.0) {
        System.out.println("DataPoint: attempt to fix errors with illegal lumi value="+lumi);
        return;
    }

    if (getY()==0.0 && getYUp()==0.0 && getYDown()==0.0) {
        setYUp(1.3/(getBinWidth()*lumi));
        setYDown(1.3/(getBinWidth()*lumi));
    }
    }
}

Generated Wed Jan 17 09:14:27 GMT 2007