#########################
Script for data treatment
#########################

::

   proj.u=[1,0,0]; proj.v=[0,1,0]; proj.uoffset=[0,0,0,0]; proj.type='rrr';


   %==Background subtraction (including cuts, replication, binary operation)==

   my_slice=cut_sqw(sqw_file,proj,[2,0.1,6],[-1.1,-0.9],[-0.1,0.1],[0,10,1000]);%make a slice from which you wish to subtract a background
   %this makes most sense for QE slices, but can be done for QQ slices too.
   plot(my_slice);

   my_bg=cut(my_slice,[5,6],[]);%take a cut from this slice, going along the energy axis
   plot(my_bg);

   %Now replicate this 1d background into a 2d slice. Note that we must operate on dnd data objects now, because by doing this subtraction we lose
   %contact with the true original data
   my_bg_rep=replicate(d1d(my_bg),d2d(my_slice));%we replicate the 1d cut to 2d using the original slice as the template (for Q range etc)

   plot(my_bg_rep)

   my_slice_subtracted=d2d(my_slice) - my_bg_rep;%Just use the binary operation '-' to subtract the background slice
   plot(my_slice_subtracted);


   %=========== Symmetrisation =============================

   my_slice2=cut_sqw(sqw_file,proj,[-5,0.1,3],[-3,0.1,3],[-0.1,0.1],[500,600]);
   plot(my_slice2);

   %Fold along vertical (look at the proj axes above) - we specify the fold using a plane. The plane is specified by providing two non-parallel vectors that lie in it
   %the cross-product of the vectors define the plane normal. The position of the plane w.r.t. the origin is given by the 3rd vector argument
   %in this case (and most cases) it is not offset from the origin

   my_sym=symmetrise_sqw(my_slice2,[0,1,0],[0,0,1],[0,0,0]);
   plot(my_sym);

   %Two folds along diagonals
   my_sym2=symmetrise_sqw(my_slice2,[1,1,0],[0,0,1],[0,0,0]);
   my_sym2=symmetrise_sqw(my_sym2,[-1,1,0],[0,0,1],[0,0,0]);
   plot(my_sym2);

   %========== Bose factor correction ==============

   my_slice_bose=bose(my_slice,300);%data taken at 300K

   plot(my_slice_bose);