|
|
- clear all
- close all
- %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% fastP %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
- cut1=2700;cut2=18400;
- xMAX=8154;
- xf=(cut1-xMAX:1:cut2-xMAX);
- xcmfp=0;
- xcmfm=0;
- ymeanP=xf*0;
- ymeanM=xf*0;
- %%%%%%%%% legge i file con eterodina positiva fast %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
- limBKGD=7000;
- filecommon='caf2_maggio11_1757';
- figure(50)
- fi=dir([filecommon '*']);
- index1=[];
- index2=[];
-
- for kk=[1:10];
- index1=[index1,kk];
- kk
- file=[filecommon '_' num2str(kk)];
- str=['load ' file '.dat'];eval(str);
- str=['x=' file '(:,1);'];eval(str);
- str=['y=' file '(:,2);'];eval(str);
- str=['N=' file '(:,3);'];eval(str);
-
- nn=find(x>cut1 & x<cut2);x=x(nn);y=y(nn);
- nn=find(x<limBKGD);y=y-mean(y(nn)); %y=y/max(y);
-
- [A,B]=butter(20,0.2,'low');
- yf=filter(A,B,y);
- nn=find(yf>0.85*max(yf));xcm=sum(x(nn).*yf(nn))/sum(yf(nn));x=x-xcm;
- y=interp1(x,yf,xf,'linear',NaN);
-
- figure(333);
- plot(xf,y);ylim([-0.02 0.02])
-
- ymeanP=ymeanP+y;
- xcmfp=xcmfp+xcm;
- str=['clear ' file ];eval(str);
- figure(50);
- hold on;plot(xf,y,'g');
- end
-
-
- hold on
- plot(xf,(ymeanM+ymeanP)/(length(index1)+length(index2)),'r');
- hold off
-
-
- YY=(ymeanM+ymeanP)/(length(index1)+length(index2));
- jj=find(isnan(YY)==0); xfnew=xf(jj); Ynew=YY(jj);
-
- jj=find(isnan(YY)==0);
- figure(101)
- xx=xf(jj)*0.52733333; yy=YY(jj)/max(YY(jj));
- plot(xx,yy)
-
-
- nn1=find(xx<-1000); nn2=find(xx>3000); nn=[nn1 nn2]; P=polyfit(xx(nn),yy(nn),1);
- Ybkg=polyval(P,xx);
-
- figure(11)
- subplot(2,1,1); semilogy(xx,yy-Ybkg,xx,yy)
- subplot(2,1,2); plot(xx,yy,xx,Ybkg)
-
- if 0,
- AA=[xx' (yy-Ybkg)']; str=['save ' filecommon '.dat AA -ascii;'];eval(str);
- end;
-
|
