ultralytics · pderrenger · Mar 28, 2025 · Mar 28, 2025
diff --git a/DETECTOR/fcnDefinePixels.m b/DETECTOR/fcnDefinePixels.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [input, flags] = fcnDefinePixels(input,flags,handles)
 if ~flags.update.detectorgeometry;  return;  end
 fprintf('Defining Detector Geometry... '); tic

diff --git a/DETECTOR/fcncylinder.m b/DETECTOR/fcncylinder.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [P, input] = fcncylinder(input)
 P = cell(1,6);
 d2r = pi/180;

diff --git a/ML/ANTINEUTRINO/fcnantineutrino.m b/ML/ANTINEUTRINO/fcnantineutrino.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function results = fcnantineutrino(input,output,handles,G1,plotflag)
 zv=zeros(1,10); results.xhat=zv; results.true=zv;  if plotflag; try closeallexcept(handles.GUI.figure1); catch; end; tic; end
 

diff --git a/ML/ANTINEUTRINO/fcngetcandidates.m b/ML/ANTINEUTRINO/fcngetcandidates.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [C, XC, TC, EC, EFC, str] = fcngetcandidates(input,N,X,T,F,particleName)
 %C=candidates, XC = candidate estimates, TC = candidate truths, EC = candidate estimate errors (xhat-xtrue)
 %N=number of photons, X=xhat estimate, T=xtrue, E=error (xhat-xtrue), F=estimator failure flag, str='antineutrino'

diff --git a/ML/FASTNEUTRON/E2neutronE.m b/ML/FASTNEUTRON/E2neutronE.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function Eneutron = E2neutronE(Evis)
 P = [-0.000169585232527955,0.00783271796930060,-0.167346627765554,2.19492596963189,-19.7793109964175,129.859482522855,-642.770385257028,2448.57093493363,-7265.02204696498,16885.4910975076,-30764.0514123043,43767.5233347055,-48213.5281542030,40568.9020684766,-25563.4422468919,11732.7579233951,-3777.11961099886,813.059509362973,-112.429581800018,13.4656322313572,-0.0243103851590198];
 

diff --git a/ML/FASTNEUTRON/doubleScatterVerification.m b/ML/FASTNEUTRON/doubleScatterVerification.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function results = doubleScatterVerification(results,input,G1,PE,handles,x,s,particle,plotflag)
 MTCflag = ischecked(handles.GUI.realdataflag);
 fiberCaptureFraction=.5;  E2PE=input.Material(5).yield*input.cube.QEmean*fiberCaptureFraction; %MeV to PE

diff --git a/ML/FASTNEUTRON/fcnfastneutron.m b/ML/FASTNEUTRON/fcnfastneutron.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function results = fcnfastneutron(input,output,handles,PE,G1,plotflag)
 MTCflag = ischecked(handles.GUI.realdataflag);
 results = [];  nz = sum(output(1).N);  if nz<19 || nz>1E4;  return;  end;  if plotflag; closeallexcept(handles.GUI.figure1); end

diff --git a/ML/FASTNEUTRON/fcnneutronke.m b/ML/FASTNEUTRON/fcnneutronke.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function ke = fcnneutronke(x)
 %x=[xyz t E; xyz t E];
 neutronmass = 939.565378; %MeV/c^2

diff --git a/ML/FASTNEUTRON/fcnneutrontheta.m b/ML/FASTNEUTRON/fcnneutrontheta.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function theta = fcnneutrontheta(ke0,e)
 ep = e(1); %energy of first proton hit
 en = ke0; %energy of neutron before first proton hit

diff --git a/ML/FASTNEUTRON/fcnplotcone.m b/ML/FASTNEUTRON/fcnplotcone.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [X,Y,Z,conevecs] = fcnplotcone(a,b,angle)
 %function plots a cone with origin 'a', pointing toward 'b', half angle 'angle'
 vec = b - a;

diff --git a/ML/FASTNEUTRON/fcnplotneutronMC.m b/ML/FASTNEUTRON/fcnplotneutronMC.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = fcnplotneutronMC(input,MC,MeVVector)
 close all; clc
 input.neutron = load('dEfit.mat');

diff --git a/ML/FASTNEUTRON/fcnplotneutronTS.m b/ML/FASTNEUTRON/fcnplotneutronTS.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = fcnplotneutronTS(input,MC,tsv)
 if nargin==1
       [fname, pname] = uigetfile('*.mat','Pick a TS file'); 

diff --git a/ML/FASTNEUTRON/testOrderStatistic.m b/ML/FASTNEUTRON/testOrderStatistic.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = testOrderStatistic(input)
 A=input.table.smearExp;
 

diff --git a/ML/FERMAT/fcnfermatpointvectorized.m b/ML/FERMAT/fcnfermatpointvectorized.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [fx, gx, F, sourceprobability, F1] = fcnfermatpointvectorized(k,X)
 %X = [xyztw...xyztw]
 gx      = []; %gradient

diff --git a/ML/FERMAT/fcnsolidangle.m b/ML/FERMAT/fcnsolidangle.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function sa = fcnsolidangle(r,A)
 %A=area of surface (units^2)
 %r=radius (units)

diff --git a/ML/FERMAT/fcnsolidanglevectorized.m b/ML/FERMAT/fcnsolidanglevectorized.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [fall, r] = fcnsolidanglevectorized(k,X)
 %mode = 1; %1=simple, 2=reflect
 np = size(X,1);

diff --git a/ML/FIBER/fiberMuonFitter.m b/ML/FIBER/fiberMuonFitter.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function results = fiberMuonFitter(input,output,flags,PE,G1,handles,x,s,plotflag)
 results=[];
 [~, si] = sortrows(x,4);  x=x(si,:);  s=s(si,:);  %sort by time

diff --git a/ML/FIBER/fiberNeutronFitter.m b/ML/FIBER/fiberNeutronFitter.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function results = fiberNeutronFitter(input,output,flags,PE,G1,handles,plotflag)
 results.protonprotonflag = [0 0];
 results.minangleerror = 0;

diff --git a/ML/FIBER/fitLineND.m b/ML/FIBER/fitLineND.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [mu,V] = fitLineND(x,w)
 %x = [n x nd] vector, nd=number dimensions
 if numel(x(:,1))==1

diff --git a/ML/FIBER/weightedcov.m b/ML/FIBER/weightedcov.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function C = weightedcov(x,w)
 %x = NxM, weights = Nx1
 [N,M]=size(x);

diff --git a/ML/GAMMA/fcnKleinNishina.m b/ML/GAMMA/fcnKleinNishina.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = fcnKleinNishina()
 clc; clear all;  close all
 me = 0.510998910; %electron rest mass (MeV) or (MeV/c^2)

diff --git a/ML/GAMMA/plotSkyMap.m b/ML/GAMMA/plotSkyMap.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = plotSkyMap(F,v,angles,np,particleName)
 nr = 100;
 nc = 200;

diff --git a/ML/MUON/fcnmuon.m b/ML/MUON/fcnmuon.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function results = fcnmuon(input,output,handles,photons,G1,plotflag)
 results=[];  MTCflag=ischecked(handles.GUI.realdataflag);  if plotflag; closeallexcept([handles.GUI.figure1, findobj('Name','Event Viewer')]); end
 nz = output(1).Nsum;  if nz<30 || nz>30E4;  return;  end

diff --git a/ML/POISSON/fcnMLenergy.m b/ML/POISSON/fcnMLenergy.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [x, fx] = fcnMLenergy(input,p0,zN,x0) %vectorized by default
 if sum(zN>0)<3; x=0; fx=0; return; end
 k = fcnoptimizerk(input,zN);  k.reflections = 0;

diff --git a/ML/POISSON/fcnMLpoint.m b/ML/POISSON/fcnMLpoint.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [x, fx] = fcnMLpoint(input,output,x0) %vectorized by default
 k = fcnoptimizerk(input,output);  k.reflections = 1;
 zN = k.nhpp;

diff --git a/ML/POISSON/runMLpoint.m b/ML/POISSON/runMLpoint.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function results = runMLpoint(input,output,handles,G1,plotflag)
 results=[];  MTCflag = ischecked(handles.GUI.realdataflag); if MTCflag; A=input.MTC.A; ei=input.eventNumber+1; end; if plotflag; try closeallexcept(handles.GUI.figure1); catch; end; end; tic
 

diff --git a/ML/fcnBuffer2Memory.m b/ML/fcnBuffer2Memory.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [input, output] = fcnBuffer2Memory(input,flags,handles,PE)
 MTCflag = ischecked(handles.GUI.realdataflag);
 

diff --git a/ML/fcnanalogvoltage.m b/ML/fcnanalogvoltage.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [pdf, strip, pid] = fcnanalogvoltage(input,flags,PE,t,ystrip,i,addnoiseflag)
 pmt = input.cube.pmt;
 

diff --git a/ML/fcnguessinitenergy.m b/ML/fcnguessinitenergy.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
  function e0hat = fcnguessinitenergy(k,zN,p)
 %p = sqrt(zN)'*input.cube.all.xyz/sum(sqrt(zN)); %position guess
 k.reflections=0;

diff --git a/ML/fcnoptimizerk.m b/ML/fcnoptimizerk.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function k = fcnoptimizerk(input, output)
 %MTCflag = ischecked(evalin('base','handles.GUI.realdataflag'));
 MTC766flag = false;

diff --git a/ML/fcnoptimizerkr.m b/ML/fcnoptimizerkr.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function k = fcnoptimizerkr(input,output)
 %reduced version of optimizer constants using only active pixels
 k = fcnoptimizerk(input,output(1).N);

diff --git a/ML/oneBookend.m b/ML/oneBookend.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function A = oneBookend(F,I,D,PEs,t0,gain)
 MTC766flag=false;
 if MTC766flag

diff --git a/MOMENTUM/momentum.m b/MOMENTUM/momentum.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 clc; close; format long g
 
 MeV = 11;

diff --git a/MOMENTUM/momentum1event.m b/MOMENTUM/momentum1event.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 clc; clear; close; format long g
 MeV = 4;
 e = 4; %event number

diff --git a/MOMENTUM/plotPointStats.m b/MOMENTUM/plotPointStats.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = plotPointStats(G)
 
 p1 = zeros(1,3);

diff --git a/MTC/DataAnalysis/example_mtcb_commands.sh b/MTC/DataAnalysis/example_mtcb_commands.sh
@@ -1,3 +1,5 @@
+# Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 # LOGIN TO SERVER:
 ssh -X -p 25260 m033372@mtc-b.phys.hawaii.edu
 

diff --git a/MTC/DataAnalysis/fcnloadMTCQT.m b/MTC/DataAnalysis/fcnloadMTCQT.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [input, D, es] = fcnloadMTCQT(pathname,input,flags,handles,plotflag,ei)
 if exist('ei','var') && ~isempty(ei) || flags.status.MC;  MCflag=true; else MCflag=false; end
 if MCflag; plotflag = 0; end;  fprintf('Loading MTC data... '); tic

diff --git a/MTC/DataAnalysis/fcnprocess1run.m b/MTC/DataAnalysis/fcnprocess1run.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [A, Cout] = fcnprocess1run(A,handles,plotflag)
 try closeallexcept(handles.GUI.figure1); end %#ok<TRYNC>
 if ~exist('plotflag','var') || isempty(plotflag); plotflag=false; end

diff --git a/MTC/DataAnalysis/firstfour.m b/MTC/DataAnalysis/firstfour.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function i=firstfour(w,i)
 i = i( w<(min(w)+4) );
 end
diff --git a/MTC/DataAnalysis/getRunCategories.m b/MTC/DataAnalysis/getRunCategories.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [run,category] = getRunCategories()
 X = GetGoogleSpreadsheet('1LEZYzBG8Ae_jQcYmxaTPxrfccSzXTMRNtyozElfS98k');
 X=X(3:end,[1 13]);  

diff --git a/MTC/DataAnalysis/lastfour.m b/MTC/DataAnalysis/lastfour.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function i=lastfour(w,i)
 i = i( w>(max(w)-4) );
 end
diff --git a/MTC/DataAnalysis/plotAnalysis.m b/MTC/DataAnalysis/plotAnalysis.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = plotAnalysis()
 clc; clear; close all 
 %MTC RUNS: https://docs.google.com/spreadsheets/d/1LEZYzBG8Ae_jQcYmxaTPxrfccSzXTMRNtyozElfS98k/edit?usp=sharing

diff --git a/MTC/DataAnalysis/removeSpikes.m b/MTC/DataAnalysis/removeSpikes.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function V=removeSpikes(V,threshold)
 if nargin==1; threshold=50; end
 

diff --git a/MTC/DataAnalysis/runAnalysis.m b/MTC/DataAnalysis/runAnalysis.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = runAnalysis()
 close all; clc; clear all %#ok<CLALL>
 load M2.mat;  input.MTC.A.x=[]; input.MTC.C=[]; %#ok<*STRNU>

diff --git a/MTC/Mapping/MTCSRCCH2pixelID.m b/MTC/Mapping/MTCSRCCH2pixelID.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function ID = MTCSRCCH2pixelID(SRCCH)
 %load('MTCmapping.mat');  %X = [SRCCH PMT PMTR PMTC ASICs ASICd Si CHs FACE];
 X=getmap;

diff --git a/MTC/Mapping/MTCpixelID2SRCCH.m b/MTC/Mapping/MTCpixelID2SRCCH.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [S, R, C, CH, PMT, PMTR, PMTC, As, Ad, Si, CHs] = MTCpixelID2SRCCH(i)
 load('MTCmapping.mat');  %X = [SRCCH PMT PMTR PMTC ASICs ASICd Si CHs FACE];
 %SRCCH(SRCCH(:,1)==0,   1)=102;

diff --git a/MTC/Mapping/fcnUHShigepixel2jocherpixel.m b/MTC/Mapping/fcnUHShigepixel2jocherpixel.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = fcnUHShigepixel2jocherpixel(input)
 C_uh2j = fcnENU2NED;  %x_j = C_uh2j*x_uh; 
 %UH - JOCHER

diff --git a/MTC/Mapping/fcnchannels2windows.m b/MTC/Mapping/fcnchannels2windows.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function a = fcnchannels2windows(a)
 %pa = 1536x32768x5 window format (for 5 pedestals)
 %pb = 786432x64x5 channel format

diff --git a/MTC/Mapping/fcnwindows2channels.m b/MTC/Mapping/fcnwindows2channels.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function a = fcnwindows2channels(a)
 %pa = 1536x32768x5 window format (for 5 pedestals)
 %pb = 786432x64x5 channel format

diff --git a/MTC/fcngmdistribution.m b/MTC/fcngmdistribution.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [mu,s,k,X] = fcngmdistribution(X)
 %fits data and tells you how many normals make up the underlying mixture distribution
 mu=[];

diff --git a/MTC/pedestals/convertData2Pedestals.m b/MTC/pedestals/convertData2Pedestals.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = convertData2Pedestals(input)
 A=input.MTC.A; input.MTC.A=[]; fname=A.filename;
 E = A.E-min(A.E)+1;

diff --git a/MTC/pedestals/fcnpruninglist.m b/MTC/pedestals/fcnpruninglist.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [b, i] = fcnpruninglist(a)
 if nargin==0; a=(1:1536)'; end
 

diff --git a/MTC/pedestals/subtractpedestals.m b/MTC/pedestals/subtractpedestals.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function V = subtractpedestals(pedestals,pedestalOutliers,pid,window,V)
 % load('MTCpedestal') %loads X 98304x71
 % pedestals = muvb(:,8:end);  %#ok<NODEF> %idxm=load('MTCmapping.mat');  b=b(idxm.X(:,1),:);

diff --git a/MTC/timebase offsets 2015/fcngettimebases.m b/MTC/timebase offsets 2015/fcngettimebases.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [] = fcngettimebases(input,handles,flags)
 closeallexcept(handles.GUI.figure1); clc;
 [pathname]=uigetdir([pwd '/*.*','MultiSelect','On'],'Select folder:'); 

diff --git a/NN/functions/conventionalFiberFit.m b/NN/functions/conventionalFiberFit.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function xhat = conventionalFiberFit(I, T, threshold)
 nw=size(I,2)/2;  vi{1}=1:nw; vi{2}=vi{1}+nw;
 if nargin==2

diff --git a/NN/functions/gatherTrainingData.m b/NN/functions/gatherTrainingData.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function results = gatherTrainingData(input,output,PE,G1)
 % mode='trainLAPPDNTC2';
 % switch mode

diff --git a/NN/functions/trainFiber64.m b/NN/functions/trainFiber64.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function trainFiber64(files)
 if nargin==0
     applyFlag=false;

diff --git a/NN/functions/trainNewDetector.m b/NN/functions/trainNewDetector.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function []=trainNewDetector()
 evalin('base','clc; clear all; close all;'); startclock=clock;
 

diff --git a/PHOTONS/cherenkov_magnitude.m b/PHOTONS/cherenkov_magnitude.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function wlpdf = cherenkov_magnitude(wl,ir,speed)
 %wl (nm)
 %speed (c)

diff --git a/PHOTONS/fcnBQF.m b/PHOTONS/fcnBQF.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function devis = fcnBQF(de,dx,BQF)
 %converts energy to visible energy using BQF 
 %dx = distance travelled in mm

diff --git a/PHOTONS/fcndetectorintercept.m b/PHOTONS/fcndetectorintercept.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [p2,dt] = fcndetectorintercept(input,p1,vel)
 %p1 is the vector origin, vel is the velocity vector
 %endPos is the vector intercept on the detector wall, dt is the time it took to get there

diff --git a/PHOTONS/fcne2wl.m b/PHOTONS/fcne2wl.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function wl = fcne2wl(e,ir) %energy (MeV)
 c = 299792458; %m/s, speed of light
 planck = 4.13566751691E-15; %eV*s

diff --git a/PHOTONS/fcngain2ratio.m b/PHOTONS/fcngain2ratio.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function ratio = fcngain2ratio(gain,type)
 %gain (dB)
 %ratio (amplitudeOut/amplitudeIn)

diff --git a/PHOTONS/fcninsidedetector.m b/PHOTONS/fcninsidedetector.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function i = fcninsidedetector(x, input, ec)
 %x is nx3 particle position
 %ec is edge cut (optional)

diff --git a/PHOTONS/fcnke2c.m b/PHOTONS/fcnke2c.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function c = fcnke2c(ke,mass)
 % kE in MeV
 % mass in MeV

diff --git a/PHOTONS/fcnphotondeconvolution.m b/PHOTONS/fcnphotondeconvolution.m
@@ -1,3 +1,5 @@
+% Ultralytics 🚀 AGPL-3.0 License - https://ultralytics.com/license
+
 function [V, QT] = fcnphotondeconvolution(input,flags,X)
 %signal = 1536x256
 if isempty(X); V=X; QT=[]; return; end