Changeset 3932 in Sophya for trunk/Cosmo/RadioBeam/mdish.cc

Timestamp:

Dec 23, 2010, 7:33:37 PM (15 years ago)

Author:

ansari

Message:

amelioration repicon.cc (decodage arguments) , Reza 23/12/2010

File:

• trunk/Cosmo/RadioBeam/mdish.cc (modified) (2 diffs)

trunk/Cosmo/RadioBeam/mdish.cc

@@ -183 +183 @@
   SetThetaPhiRange();
   SetRespHisNBins();
+  SetPrtLevel();
   mcnt_=0;
 }
@@ -207 +208 @@
 
   double dtet = thetamax_/(double)ntet_;
-  double dphi = phimax_/(double)ntet_;
+  double dphi = phimax_/(double)nphi_;
+  cout << " MultiDish::GetResponse() - ThetaMax=" << thetamax_ << " NTheta=" << ntet_ 
+       << " PhiMax=" <<  phimax_ << " NPhi=" << nphi_ << endl;
 
   double sumw = 0.;
-  for(int kt=0; kt<ntet_; kt++) 
+  for(int kt=0; kt<ntet_; kt++) {
+    double theta=(double)kt*dtet;
     for(int jp=0; jp<nphi_; jp++) {
-      sumw += CumulResp(rd, (double)kt*dtet, (double)jp*dphi);
-      sumw += CumulResp(rd, (double)kt*dtet, -(double)jp*dphi);
-      sumw += CumulResp(rd, (double)kt*dtet, (double)jp*dphi+M_PI);
-      sumw += CumulResp(rd, (double)kt*dtet, -(double)jp*dphi-M_PI);
+      double phi=(double)jp*dphi;
+      sumw += CumulResp(rd, theta, phi);
+      if (theta<1.e-9) continue;
+      sumw += CumulResp(rd, theta, -phi);
+      sumw += CumulResp(rd, theta, phi+M_PI);
+      sumw += CumulResp(rd, theta, -(phi+M_PI));
     }
+    if (prtlev_>0) cout << "  MultiDish::GetResponse() done ktheta=" << kt << " / MaxNTheta= " << ntet_ << endl; 
+  }
   double kx1 = DeuxPI*(dishes_[0].DiameterX())/lambda_;
   double ky1 = DeuxPI*(dishes_[0].DiameterY())/lambda_;