function trainCG(pres)

% get features and labels
[f,y] = sampleDetector(@detector,pres);
f=f'; y=y';
% normalize features to unit variance
fstd = std(f);
fstd = fstd + (fstd==0);
f = f ./ repmat(fstd,size(f,1),1);
% fit the model
fprintf(2,'Fitting model...\n');
beta = logist2(y,f);
% save the result
save(sprintf('beta_cg_%s.txt',pres),'fstd','beta','-ascii');

function [f] = detector(im)
cg = detCG(im);
a = max(squeeze(cg(:,:,2,:)),[],3);
b = max(squeeze(cg(:,:,3,:)),[],3);
a = a(:);
b = b(:);
f = [ ones(size(a)) a b ]';