PredictedPlot.java

Go to the documentation of this file.

package cedar.jetweb.model.plots;

import java.io.File;
import java.util.Enumeration;
import java.util.Vector;

import cedar.jetweb.db.DBPlotManager;
import cedar.jetweb.model.*;
import cedar.jetweb.*;
import cedar.jetweb.generator.Generator;

public class PredictedPlot extends DataPlot {

    private int runSeriesId=-1;
    private DataPlot comparisonPlot;
    private RunSeries runSeries;
    private boolean tested=false;
    private boolean isValid;
    private boolean dirtyData;

    public PredictedPlot(DataPlot basePlot, 
             RunSeries runSeriesIn) 
    throws JetWebException {
    
    comparisonPlot = basePlot;
    csnId=basePlot.getId();
    runSeries=runSeriesIn;
    runSeriesId=runSeriesIn.getId();
    setFrom(basePlot);  
    dirtyData = true;
    cuts_ = basePlot.getCutCollection();
    }
    //accessors


    public void setDataPoints(Vector<DataPoint> newDataPoints){
    super.setDataPoints(newDataPoints);
    }

    public PredictedPlot setRunSeries(RunSeries rs){
    runSeries = rs;
    runSeriesId = rs.getId();
    return this;
    }


    public void setRunSeriesId(int id){ runSeriesId = id; }
    public int getRunSeriesId(){ return runSeriesId; }

    public RunSeries getRunSeries(){return runSeries; }

    public DataPlot getCompPlot(){ return comparisonPlot; }
    public void setCompPlot(DataPlot newplot){ comparisonPlot=newplot; }


    //public methods
    public boolean hasData() throws JetWebException {

    if (!hasdata){// && dirtyData) {

        dataPoints = DBPlotManager.getDataPoints(this);
        if (dataPoints.isEmpty()){
        hasdata = false;
        System.out.println("No data for MC plot " 
                   + number + ":" 
                   + title+" RSId="+runSeriesId);
        
        } else {
        hasdata = true;
        }
        dirtyData = false;
        }

    return hasdata;
    }
    public void toDB(int update) throws JetWebException {

        //check whether crossection with this plot id exists
    //if it does, update crossection record,otherwise create  
    //delete existing data plot (if exists)
    //insert new data plot records

    if (update==0){
        DBPlotManager.updateCrossSectionPoints(this);
    } else {
        DBPlotManager.insertData(this);
    }

    }
    //
    // Some plots give unrealistic errors for good reasons. This
    // method provides an option to replace/modify them to make them
    // more realistic.
    // This function is called directly after a plot is read from XML.
    //
    public static void adjustErrors(PredictedPlot plot) throws JetWebException {

    // Field and Stuart mean pT plots.

        
    // The pT away plot.
    if (plot.getId() == 115){

        for (DataPoint point : plot.getDataPoints()) {
        // Errors are dodgy for points above about 10 GeV.
        if (point.getX()>10 && point.getY()!=0){
            double y = point.getY();
            point.setYUp( y/2.0 + point.getX()/6.0 );
            point.setYDown(point.getYUp());
        }
        }

    } else if (plot.getId() == 115){

        for (DataPoint point : plot.getDataPoints()){

        // Errors are dodgy for points above about 10 GeV.
        if (point.getX()>10 && point.getY()!=0){
            double y = point.getY();
            point.setYUp( y/2.0 + 1.5 );
            point.setYDown(point.getYUp());
        }
        }
    }
    }
     public void add(PredictedPlot plot2) throws JetWebException {

         Enumeration<DataPoint> points  = this.getDataPoints().elements();
     Enumeration<DataPoint> points2 = plot2.getDataPoints().elements();
     while (points.hasMoreElements()){
         double sig1sq=0, sig2sq=0, w1=0, w2=0;

         DataPoint p1 = (DataPoint)points.nextElement();
         DataPoint p2 = (DataPoint)points2.nextElement();

         p1.add(p2);
     }
     }
    public void fixZeroErrors(double lumi) throws JetWebException {

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

    for (DataPoint point : getDataPoints()){
        point.fixZeroError(lumi);
    }
    }

    public boolean isPredictedBy(RunSeries rs){
    return (runSeriesId==rs.getId());
    }


    public MCProcessType getMCProcessType() throws JetWebException {

    if (proc==null) {
        proc = new MCProcessType(DBPlotManager.getProcId(runSeriesId));
    }
    return proc;

    }

    public boolean validMCProcessType() throws JetWebException {

    isValid = (getProcIdWanted()==getMCProcessType().getId());

    //System.out.println("Summary of validity test:"+procIdWanted+","
    //         +getMCProcessType().getId());
    //System.out.println("Plot ID:"+getId()+" "+runSeriesId);

    return isValid;
    }

    public int getProcIdWanted() throws JetWebException {

    if (procIdWanted!=0) return procIdWanted;

    // Get the runSeries if necessary.
    if ( (runSeries != null) && runSeries.getId()==runSeriesId ){
        // Do nothing. The run series exists & is up-to-date.
    } else {
        //runSeries = new RunSeries(runSeriesId);
        runSeries = RunSeries.Maker(runSeriesId);
    }

    return getProcIdWanted(runSeries.getGenerator());
    }


}

Generated Wed Jan 17 09:14:27 GMT 2007