/*
 * StarData.java
 *
 * Created on June 13, 2002, 11:00 PM
 */

//package readstardump;

import java.io.*;
import jas.hist.*;


/**
 *
 * @author  karlen
 * @version 
 */
public class StarDumpData implements Rebinnable1DHistogramData {

    static final int nChan = 1152;
    int[][] data = new int[nChan][];
    int nBin;
    boolean goodData = false;
    int channel = 0;
    
    /** Creates new StarData */
    public StarDumpData(File readFile, int nBin) {
        this.nBin = nBin;
        for (int i=0; i<nChan; i++) {
            data[i] = new int[nBin];
        }
        
        // read in the data: format depends on file extension //
        String name = readFile.getName();
        String ext = name.substring(name.lastIndexOf(".")+1);
        if (ext.equals("dat")) {
            goodData = readData(readFile);
        } else {
            goodData = readIrqData(readFile);
        }
    }
            
    boolean readData(File readFile) {
        try {

            BufferedReader in = new BufferedReader(new FileReader(readFile));            
            String line = in.readLine();
            // skip header information:
            while (!line.substring(0,9).equals("c0001000:")) {
                line = in.readLine();
            }
            
            int iCh = -1;
            int iBin = 0;
            int nDat = nBin*nChan;
            int start = 6;
            for (int i=0; i<nDat; i++) {
                iCh++;
                if (iCh >= nChan) {
                    iCh = 0;
                    iBin++;
                }
                start += 5;
                if (start > 47) {
                    line = in.readLine();
                    start = 11;
                }
                data[iCh][iBin] = Integer.parseInt(line.substring(start,start+4),16);
            }
        } catch (Exception e) {
            System.err.println("Exception " +e+ " while reading file.");
            return false;
        }                
        return true;
    }

    boolean readIrqData(File readFile) {
        try {
	    // Read STAR IRQ data
            BufferedReader in = new BufferedReader(new FileReader(readFile));            
            String line = in.readLine();

            // skip header information:
            while (line.length() < 40 || !line.substring(0,10).equals("0xc0001000")) {
		if (line.length() > 40 && line.substring(0,10).equals("0xc0000030")) {
		    //System.out.println("  Check for dead face in "+line);
		    if (line.substring(17,19).equals("de") &&
			line.substring(26,28).equals("ad") &&
			line.substring(35,37).equals("fa") &&
			line.substring(44,46).equals("ce")) {
			//System.out.println("  Found dead face");
		    } else {
			System.out.println("  Did NOT find dead face in "+line);
			return false;
		    }
		}
                line = in.readLine();
            }
            
            int iCh = -1;
            int iBin = 0;
            int nDat = nBin*nChan;
            int start = 2;
            for (int i=0; i<nDat; i++) {
                iCh++;
                if (iCh >= nChan) {
                    iCh = 0;
                    iBin++;
                }
                if (iCh%2 == 0) {
                    start += 13;
                } else {
                    start -= 4;
                }
                if (start > 62) {
                    line = in.readLine();
                    start = 15;
                }
                data[iCh][iBin] = Integer.parseInt(line.substring(start,start+4),16);
            }
        } catch (Exception e) {
            System.err.println("Exception " +e+ " while reading file.");
            return false;
        }                
        return true;
    }
    
    public boolean isGood() {
        return goodData;
    }
    
    public int getLoc(int iLoc, String half, int iPin, int iBin) {
	int iHalf = 0;
	if (half.equals("B")) iHalf = 1;
	return getLoc(iLoc,iHalf,iPin,iBin);
    }
    public int getLoc(int iLoc, int iHalf, int iPin, int iBin) {
	channel = mapSequence[iHalf][iLoc] + iPin*72;
	return data[channel][iBin];
    }
    void setChannel(int channel) {
        this.channel = channel;
    }
    
    void setChannel(int iLoc, int iHalf, int iPin) {
        channel = mapSequence[iHalf][iLoc] + iPin*72;
    }
    
    public double[][] rebin(int rBins,double rMin,double rMax,
                            boolean wantErrors,boolean hurry) {
        double[][] result = new double[3][nBin];
        for (int i=0; i<nBin; i++) {
            result[0][i] = (double) data[channel][i];
        }
        return result;
    }
    
    public int getBins() {
        return nBin;
    }
    
    public boolean isRebinnable() {
        return false;
    }
    
    public int getAxisType() {
        return DOUBLE;
    }
    
    public java.lang.String[] getAxisLabels() {
        return null;
    }
    
    public double getMax() {
        return (double) nBin;
    }
    
    public double getMin() {
        return 0.;
    }
    
    public java.lang.String getTitle() {
        return "Channel " + channel;
    }
    
    static final int[][] mapSequence = {{0,1,2,3,4,5,18,19,20,54,55,56,57,58,59,6,7,8,            // A
			           42,43,44,45,46,47,60,61,62,30,31,32,33,34,35,48,49,50},
				  {21,22,23,36,37,38,39,40,41,9,10,11,24,25,26,27,28,29,    // B
				   63,64,65,12,13,14,15,16,17,51,52,53,66,67,68,69,70,71}};
    
}