JetWebAIDAXMLPlotReader.java

Go to the documentation of this file.

package cedar.jetweb.xml;

import java.util.Vector;
import cedar.jetweb.model.plots.DataPoint;
import cedar.jetweb.model.plots.DataPlot;
import cedar.jetweb.model.plots.PredictedPlot;
import cedar.jetweb.JetWebException;
import hep.aida.*;
import hep.aida.ref.plotter.PlotterUtilities;

public class JetWebAIDAXMLPlotReader implements JetWebPlotReader {
    IDataPointSet m_data_ = null;
    ITree tree_ = null;
    String directory_ = "";
    IDataPointSetFactory dpsf_ = null;
    IHistogram1D hist_ = null;

    public JetWebAIDAXMLPlotReader(String filename) throws JetWebXMLException{
       
        //filename is an absolute path apparently
        try{
            IAnalysisFactory af = IAnalysisFactory.create();
            
            tree_ = af.createTreeFactory().create(filename,"xml",false,false);
            dpsf_ = af.createDataPointSetFactory(tree_);

        /*
     * This needs to reflect the directory structure of the tree created with the 
     * datapoints. Dont think the JAIDA xml stores this info so work needs to be done

        m_data = dpsf.create("/data","",2);
     */
        }catch (Exception e) {
            throw new JetWebXMLException(e);
        }
    }
    
    private Vector<DataPoint> getDataPoints() {
            
    Vector<DataPoint> dps = new Vector();

    /*
    int size = m_data_.size();

        
    for(int i=0;i<size;i++) {
        IDataPoint p = m_data_.point(i);
        double x = p.coordinate(0).value();
        double y = p.coordinate(1).value();
        double xerrup = p.coordinate(0).errorPlus();
        double xerrdown = p.coordinate(0).errorMinus();
        double yerrup = p.coordinate(1).errorPlus();
        double yerrdown = p.coordinate(1).errorMinus();

        System.out.println("x "+x+" xerr "+xerrup);
        System.out.println("y "+x+" yerr "+yerrup);
        
        DataPoint dp = new DataPoint();
        dp.setX(x);
        dp.setY(y);
        dp.setYUp(yerrup);
        dp.setYDown(yerrdown);
        dp.setBinWidth(1.0);
        dp.setChi2(0.0);
        dps.add(dp);

    }
    */

    IAxis axis = hist_.axis();
    for (int iBin=0 ; iBin<axis.bins(); iBin++){

        DataPoint dp = new DataPoint();
        dp.setX( (axis.binUpperEdge(iBin)+axis.binLowerEdge(iBin))/2.0 );
        dp.setY(hist_.binHeight(iBin + 1 ));//iBin));
        dp.setYUp(hist_.binError(iBin+1));
        dp.setYDown(hist_.binError(iBin+1));
        dp.setBinWidth(axis.binWidth(iBin));
        dp.setChi2(0.0);

        //System.out.println("adding data point: " + dp);

        dps.add(dp);
    }
        
    return dps;
    }
    
    public void setDirectory(String dir) throws JetWebException {
    directory_ = dir;
    }

    public void fillPlot(DataPlot plot) throws JetWebException {
    //set the title - already set by plotfactory
    //set xlabel - same
    //set ylabel -same
    
    try {

        String location = directory_ + "/h" + plot.getNumber();

        System.out.println("using histogram at " + location);

        hist_ = (IHistogram1D) tree_.find(directory_+"/h"+plot.getNumber());
               
        // Create IDataPointSets from the the above AIDA objects.
        m_data_   = dpsf_.create("histogram",hist_);
        
        plot.setDataPoints(getDataPoints());
        
        if(plot instanceof PredictedPlot) {
        PredictedPlot.adjustErrors((PredictedPlot)plot);
        }

        plot.setHasData(true);

    } catch (IllegalArgumentException e) {
        System.out.println("JetWebAIDAPlotReader: couldnt find "+directory_+"/h"+plot.getNumber());
        throw new JetWebException(e,
                      "JetWebAIDAPlotReader: couldnt find "
                      +directory_+"/h"+plot.getNumber());
    }


    }
}

Generated Wed Jan 17 09:14:27 GMT 2007