/*- * Copyright 2015 Diamond Light Source Ltd. * * All rights reserved. This program and the accompanying materials * are made available under the terms of the Eclipse Public License v1.0 * which accompanies this distribution, and is available at * http://www.eclipse.org/legal/epl-v10.html */ package uk.ac.diamond.scisoft.analysis.processing.operations; import org.eclipse.dawnsci.analysis.api.processing.Atomic; import org.eclipse.dawnsci.analysis.api.processing.OperationData; import org.eclipse.dawnsci.analysis.api.processing.OperationException; import org.eclipse.dawnsci.analysis.api.processing.OperationRank; import org.eclipse.dawnsci.analysis.dataset.operations.AbstractOperation; import org.eclipse.january.IMonitor; import org.eclipse.january.dataset.Dataset; import org.eclipse.january.dataset.DatasetFactory; import org.eclipse.january.dataset.DatasetUtils; import org.eclipse.january.dataset.DoubleDataset; import org.eclipse.january.dataset.IDataset; import org.eclipse.january.dataset.IndexIterator; @Atomic public class GeneralDetectorErrorOperation extends AbstractOperation<GeneralDetectorErrorModel, OperationData> { @Override public String getId() { return "uk.ac.diamond.scisoft.analysis.processing.operations.GeneralDetectorErrorOperation"; } @Override protected OperationData process(IDataset input, IMonitor monitor) throws OperationException { DoubleDataset er = DatasetFactory.zeros(DoubleDataset.class, input.getShape()); Dataset in = DatasetUtils.convertToDataset(input.getSliceView()); IndexIterator i = in.getIterator(); double gain = model.getGain(); double gainF = model.getGainFluctuation(); double qE = model.getQuantumEfficiency(); double eNoise = model.getElectronicNoise(); double val = 0; while (i.hasNext()) { val = in.getElementDoubleAbs(i.index); val = (gain*gain+gainF*gainF)*qE*val+eNoise; if (val < 0) continue; er.setAbs(i.index, Math.sqrt(val)); } in.setErrors(er); return new OperationData(in); } @Override public OperationRank getInputRank() { return OperationRank.TWO; } @Override public OperationRank getOutputRank() { return OperationRank.TWO; } }