matlabProj

rawdata = importdata('input.txt');
%photon = rawdata(isphoton~=0,:);
x = rawdata(:,2); y = rawdata(:,3); xcos = rawdata(:,5); ycos = rawdata(:,6); E = rawdata(:,7);

xround = round(x);
yround = round(y);
xcosround = round(xcos.*100)./100;
ycosround = round(ycos.*100)./100;

xidx = unique(xround);
yidx = unique(yround);
xcosidx = unique(xcosround);
ycosidx = unique(ycosround);

xarray = -1:1;
yarray = -1:1;
xcosarray = -1:0.01:1;
xcosarray = round(xcosarray.*100)./100;
ycosarray = -1:0.01:1;
ycosarray = round(ycosarray.*100)./100;

Egrid = nan(length(xarray),length(yarray),length(xcosarray),length(ycosarray));
for xid = 1:length(xarray)
    for yid = 1:length(yarray)
        for xcosid = 1:length(xcosarray)
            for ycosid = 1:length(ycosarray)
                oneid = xround==xarray(xid) & yround==yarray(yid) & xcosround==xcosarray(xcosid) & ycosround==ycosarray(ycosid);
                oneid = oneid(oneid~=0);
                if(~isempty(oneid))
                    Egrid(xid,yid,xcosid,ycosid) = E(oneid(1));
                end
            end
        end
    end
end

nF = max(Egrid(~isnan(Egrid)));
Egrid = Egrid./nF;
Egrid = permute(Egrid,[4,3,2,1]);

Egrid(isnan(Egrid)) = -1;
tid = fopen('Egrid.bin','w');
fwrite(tid,Egrid,'float');
fclose(tid);