fact.features.Leakage Maven / Gradle / Ivy

Go to download

Show more of this group Show more artifacts with this name
Show all versions of fact-tools Show documentation

The FACT telescope is a Cherenkov telescope on the island of La Palma. This set of tools is dedicated to read, parse and process raw data produced by the FACT telescope.

The newest version!

package fact.features;

import fact.Utils;
import fact.hexmap.FactCameraPixel;
import fact.hexmap.FactPixelMapping;
import org.slf4j.Logger;
import org.slf4j.LoggerFactory;
import stream.Data;
import stream.Processor;
import stream.annotations.Parameter;

public class Leakage implements Processor {
	static Logger log = LoggerFactory.getLogger(Leakage.class);

    @Parameter(required = true)
	private String shower;
    @Parameter(required = true)
	private String weights;
    @Parameter(required = true)
	private String leakage1OutputKey;
    @Parameter(required = true)
	private String leakage2OutputKey;

    FactPixelMapping pixelMap = FactPixelMapping.getInstance();

	@Override
	public Data process(Data input) {
		Utils.mapContainsKeys( input, shower, weights);
	
		int[] 	showerPixel = (int[])input.get(shower);
		double[] photonCharge = (double[]) input.get(weights);
		
		
		double size = 0;
	
	    double leakageBorder          = 0;
	    double leakageSecondBorder    = 0;

	    for (int pix: showerPixel)
	    {
	    	size += photonCharge[pix];
	        if (isBorderPixel(pix) )
	        {
	            leakageBorder          += photonCharge[pix];
	            leakageSecondBorder    += photonCharge[pix];
	        }
	        else if (isSecondBorderPixel(pix))
	        {
	            leakageSecondBorder    += photonCharge[pix];
	        }
	    }
	    leakageBorder          = leakageBorder        / size;
	    leakageSecondBorder    = leakageSecondBorder  / size;

		
		input.put(leakage1OutputKey , leakageBorder);
		input.put(leakage2OutputKey , leakageSecondBorder);
		return input;
		
		
	}
	
	//this is of course not the most efficient solution
	boolean isSecondBorderPixel(int pix){
		for(FactCameraPixel nPix: pixelMap.getNeighboursFromID(pix))
		{
			if(isBorderPixel(nPix.id)){
				return true;
			}
		}
		return false;
	}
	boolean isBorderPixel(int pix){
        return pixelMap.getNeighboursFromID(pix).length < 6;
	}

	
	public String getShower() {
		return shower;
	}
	public void setShower(String shower) {
		this.shower = shower;
	}
	

	public String getWeights() {
		return weights;
	}
	public void setWeights(String weights) {
		this.weights = weights;
	}

	public String getLeakage1OutputKey() {
		return leakage1OutputKey;
	}

	public void setLeakage1OutputKey(String leakage1OutputKey) {
		this.leakage1OutputKey = leakage1OutputKey;
	}

	public String getLeakage2OutputKey() {
		return leakage2OutputKey;
	}

	public void setLeakage2OutputKey(String leakage2OutputKey) {
		this.leakage2OutputKey = leakage2OutputKey;
	}

	

}