Changeset 1148 in Sophya for trunk/SigPredictor/sigcalctools.cc

Timestamp:

Aug 29, 2000, 12:03:24 PM (25 years ago)

Author:

ansari

Message:

mise a jour

File:

• trunk/SigPredictor/sigcalctools.cc (modified) (7 diffs)

trunk/SigPredictor/sigcalctools.cc

@@ -6 +6 @@
 #include <fstream>
 #ifdef __MWERKS__
-   #include "mwerksmath.h"
+    
    #include "unixmac.h"
 #endif
@@ -35 +35 @@
 
 SigCalcTool::SigCalcTool(AbsLightSource* pLightSrc, AbsLobeNoPolar* pLobeNoPolar, 
-        SpectralResponse* pFilt):pLSrc(pLightSrc),pLobe(pLobeNoPolar),pFilter(pFilt)
-{       
+        SpectralResponse* pFilt):pLSrc(pLightSrc)
+{       pLobe=pLobeNoPolar;
+        pFilter=pFilt;
+        
         SigCalcToolInit(); 
 }
@@ -94 +96 @@
 }
 
-double SigCalcTool::compPixel(double theta,double phi) {
-        UnitVector VP(theta,phi);
-        UnitVector VYbidon=VP.VperpPhi();         
-// Compute unit vector perpendicular to Vpoin at same theta
-        return compPixel(VP,VYbidon);
-}
+
 
 double SigCalcTool::compPixel(UnitVector& VP, UnitVector& VdirectY){
@@ -273 +270 @@
 }
 
-double SigCalcTool::NormKelvinRayleighJeans()
-{       
-        double tempeCNoir=10000.;       
-                                        // Kelvin
-        double CutFreq=1.380662e-23*tempeCNoir/6.626176e-34/5.;
-        if(FreqMax>1.380662e-23*tempeCNoir/6.626176e-34/5.) 
-        {       cerr<< "RaleighJeans approximation is not valid for this frequency"<<endl;
-                cerr<< "Frequency: "<< FreqMax<<" in SigCalcTool::NormRayleighJeans"<<endl;
-        }
-        
-        LightBlackBody CorpsNoir(tempeCNoir, RAngComp);
-        SigCalcTool ToolRJ(&CorpsNoir,pLobe,pFilter);
-        double puissNorm = ToolRJ.compPixel(M_PI/2.,M_PI);      // Un pixel au hasard
-        return tempeCNoir/puissNorm;                            // Kelvin RaleighJeans/(Watt/m2)
- 
-}
- 
-double SigCalcTool::NormKelvinCMB()
-{       
-        double deltatempeCNoir=1.;                                      // Kelvin
-        LightNormTCMB DeltaCorpsNoir(deltatempeCNoir, RAngComp);
-        SigCalcTool ToolDeltaCMB(&DeltaCorpsNoir,pLobe,pFilter);
-        double puissNorm = ToolDeltaCMB.compPixel(M_PI/2.,M_PI);        // Un pixel au hasard
-        return deltatempeCNoir/puissNorm;                       // KelvinCMB/(Watt/m2)
- 
-}
-
 /*
 double SigCalcTool::CalcInBandPower(double theta, double phi)
@@ -322 +292 @@
   else AngRes=RAngComp*pLobe->ResolutionCurve(angle);
   return AngRes;
-}
-
-
-
-double SigCalcTool::max(double a, double b) const{
-  if(a>b) return a;
-  else return b;
-}
-
- double SigCalcTool::min(double a, double b) const{
-  if(a<b) return a;
-  else return b;
 }
 
@@ -432 +390 @@
   
   return SizeInteg;
-}
- 
-double SigCalcTool::diffSolidAng(double ang1,double ang2) const
-{       double returnVal;               // Steradians
-        if(!pLSrc->IsPtSourceS()) returnVal= fabs(2*M_PI*(cos(ang1)-cos(ang2)));
-                // Cas d'une source Žtendue.
-        else  returnVal= 1.;
-                // Cas d'une source ponctuelle. Sont flux est exprime en W/m2/Hz
-                // Pas d'angle solide. 
-        return returnVal;
-                
 }
 
@@ -513 +460 @@
 }
 
+double SigCalcTool::diffSolidAng(double ang1,double ang2) const
+{       double returnVal;               // Steradians
+        if(!pLSrc->IsPtSourceS()) returnVal= fabs(2*M_PI*(cos(ang1)-cos(ang2)));
+                // Cas d'une source Žtendue.
+        else  returnVal= 1.;
+                // Cas d'une source ponctuelle. Sont flux est exprime en W/m2/Hz
+                // Pas d'angle solide. 
+        return returnVal;
+}
+
 template void addInInBandPowerMap(PixelMap<float>& Map, SigCalcTool& tool);
 template void addInInBandPowerMap(PixelMap<double>& Map, SigCalcTool& tool);