Changeset 3501 in Sophya for trunk/Cosmo/SimLSS/cmvsinxsx.cc

Timestamp:

Jul 18, 2008, 6:50:43 PM (17 years ago)

Author:

cmv

Message:

change libMinuit2Base libMinuit2 for roo 5.20 cmv 18/07/2008

File:

• trunk/Cosmo/SimLSS/cmvsinxsx.cc (modified) (6 diffs)

trunk/Cosmo/SimLSS/cmvsinxsx.cc

@@ -39 +39 @@
  bool only_test = false;
  double pr_test = 1.e99;
-
+ int nvec=0;
+ double veccent=0., vecmax = -2.;
+
+ //--- Read arguments
  char c;
- while((c = getopt(narg,arg,"htp:n:")) != -1) {
+ while((c = getopt(narg,arg,"htp:n:v:")) != -1) {
   switch (c) {
   case 'n' :
@@ -52 +55 @@
     sscanf(optarg,"%lf",&pr_test);
     break;
+  case 'v' :
+    sscanf(optarg,"%d,%lf",&nvec,&vecmax);
+    break;
   case 'h' :
   default :
-    cout<<"cmvsinxsx [-t] [-p eps_print] [-n N] [dt,tmin,tmax (rad)]"<<endl;
+    cout<<"cmvsinxsx [-t] [-p eps_print] [-n N] [-v nv,(+-)vmax (rad)]"<<endl;
+    cout<<"          [dt,tmin,tmax (rad)]"<<endl;
     return -1;
   }
@@ -62 +69 @@
  cout<<"N="<<N<<" dt="<<dt<<" tmin="<<tmin<<" tmax="<<tmax<<" rad"<<endl;
  cout<<"only_test="<<only_test<<" pr_test="<<pr_test<<endl;
+ if(vecmax<0.) vecmax *= -M_PI;
+ cout<<"nv="<<nvec<<" vmax="<<vecmax<<" -> ["<<-vecmax/2<<","<<vecmax/2<<"]"<<endl;
  if(tmax<=tmin) return -1;
 
+ //--- Create NTuple
  const int nvar = 11;
  char *vname[nvar] = {"t","s","sn","sx","sx2","snx","snx2","asx","asx2","asnx","asnx2"};
@@ -69 +79 @@
  NTuple nt(nvar,vname);
 
+ //--- Fill NTuple and do precision test
  long n = 0;
  double diffsx=0., diffsx2=0., diffsnx=0., diffsnx2=0.;
@@ -115 +126 @@
  cout<<"     (sNx/x)^2 = "<<diffsnx2<<" = "<<diffsnx2/(N*N)<<" * "<<N<<" * "<<N<<endl;
 
- if(!only_test) {
-   cout<<"writing ppf"<<endl;
-   POutPersist pos("cmvsinxsx.ppf");
-   pos.PutObject(nt,"nt");
+ //--- Fill vector for FFT
+ TVector<r_8> vsx,vsx2,vsnx,vsnx2;
+ if(nvec>1) {
+   vsx.ReSize(nvec); vsx2.ReSize(nvec); vsnx.ReSize(nvec); vsnx2.ReSize(nvec);
+   for(int i=0;i<nvec;i++) {
+     double t = veccent -vecmax/2. + i*vecmax/(double)nvec;
+     vsx(i)   = SinXsX(t);
+     vsx2(i)  = SinXsX_Sqr(t);
+     vsnx(i)  = SinNXsX(t,N);
+     vsnx2(i) = SinNXsX_Sqr(t,N);
+   }
  }
+
+ //--- Write PPF
+ POutPersist pos("cmvsinxsx.ppf");
+ if(!only_test) pos.PutObject(nt,"nt");
+ if(vsx.Size()>0) pos.PutObject(vsx,"vsx");
+ if(vsx2.Size()>0) pos.PutObject(vsx2,"vsx2");
+ if(vsnx.Size()>0) pos.PutObject(vsnx,"vsnx");
+ if(vsnx2.Size()>0) pos.PutObject(vsnx2,"vsnx2");
 
  return 0;
@@ -161 +187 @@
 n/plot nt.snx2-asnx2%t ! ! "connectpoints"
 
+#-------
+openppf cmvsinxsx.ppf
+
+zone 2 2
+n/plot vsx.val%n ! ! 
+n/plot vsx2.val%n ! ! 
+n/plot vsnx.val%n ! ! 
+n/plot vsnx2.val%n ! !
+
+zone
+fftforw vsx fvsx
+fftforw vsx2 fvsx2
+fftforw vsnx fvsnx
+fftforw vsnx2 fvsnx2
+
+zone
+n/plot fvsx.val%n ! ! "connectpoints"
+n/plot fvsx2.val%n ! ! "connectpoints"
+n/plot fvsnx.val%n ! ! "connectpoints"
+n/plot fvsnx2.val%n ! ! "connectpoints"
+
 */