Changeset 905 in Sophya

Timestamp:

Apr 13, 2000, 12:12:47 PM (25 years ago)

Author:

ansari

Message:

Sophie: adding the dipole here also

File:

• trunk/SophyaProg/PMixer/tgsky.cc (modified) (4 diffs)

trunk/SophyaProg/PMixer/tgsky.cc

@@ -29 +29 @@
     cout << "   - TypSky = 1  Params= m,sig  , OOFNoise(sig) + m" << endl;
     cout << "   - TypSky = 2  Params= K,a,b,m  K*cos(a*teta)*sin(b*phi)+m" << endl;   
+    cout << "   - TypSky = 3  Params= theta,phi DIPOLE" << endl;   
     exit(0);
     }
@@ -40 +41 @@
   int typ = atoi(arg[2]);
   double m,s,a,b,K;
+  double Theta,Phi;
+  Theta = 0.;
+  Phi   = 0.;
   m = 0.;
   s = 1.;
   a = b = 1.;
   K = 1.;
-  if ((typ < 0) || (typ > 2)) typ = 0;
+  if ((typ < 0) || (typ > 3)) typ = 0;
   if (typ < 2)  { 
     sscanf(arg[3],"%lg,%lg",&m,&s);
@@ -51 +55 @@
     else cout << "  ---> NoiseGenerator() " << endl;
   }
-  else { 
+  else if(typ==2) { 
     sscanf(arg[3],"%lg,%lg,%lg,%lg",&K,&a,&b,&m);
     cout << "  TypSky= " << typ << " K= " << m << " a= " << a << " b=" 
          << b << " m= " << m << endl;
     cout << "    ---> K * cos(a*teta) * sin(b*teta) + m" << endl;
+    }
+  else if(typ==3) { 
+    sscanf(arg[3],"%lg,%lg",&Theta,&Phi);
+    cout << "  TypSky= " << typ << " Theta= " << Theta << " phi= " << Phi ;
+    cout << "    ---> DIPOLE " << endl;
     }
 
@@ -66 +75 @@
     }
   }
-  else {
+  else if(typ==3)
+    {
+      UnitVector vd(Theta,Phi);
+      UnitVector vc(Theta,Phi);
+      for(int i=0; i<sph.NbPixels(); i++)
+        {
+          double Thetar,Phir;
+          sph.PixThetaPhi(i,Thetar,Phir);      
+          vc.SetThetaPhi(Thetar,  Phir);
+          sph(i) += vd.Psc(vc);
+        }
+    }
+  else{
     NoiseGenerator * ng;
     if (typ == 0) ng = new NoiseGenerator(s);