Changeset 3435 in Sophya for trunk/SophyaPI

Timestamp:

Dec 11, 2007, 5:09:37 PM (18 years ago)

Author:

ansari

Message:

Ajout fonctions tetphi2mollX/Y() et longlat2mollX/Y() ds les fichiers d'expression ntuple generes par piapp - Reza 11/12/2007

File:

• trunk/SophyaPI/PIext/servnobjm.cc (modified) (2 diffs)

trunk/SophyaPI/PIext/servnobjm.cc

@@ -1094 +1094 @@
 fputs("#define gaurand() NorRand() \n\n", fip);
 
+fputs("/* Angle conversions   */\n", fip);
+fputs("double deg2rad(double x) { return (x*M_PI/180.); } \n", fip);
+fputs("double rad2deg(double x) { return (x*180./M_PI); } \n", fip);
+fputs("/* Theta(0..Pi, phi(0..2Pi) conversion to Molleweide X,Y   */\n", fip);
+fputs("double tetphi2mollX(double teta, double phi) { \n", fip);
+fputs(" if (phi<=M_PI) return(180.-sin(acos(2*teta/M_PI-1.))*((M_PI-phi)*180./M_PI)); \n", fip);
+fputs(" else return(180.+sin(acos(2.*teta/M_PI-1.))*((phi-M_PI)*180./M_PI)); } \n", fip);
+fputs("double tetphi2mollY(double teta) { \n", fip);
+fputs(" return (90.-(teta*180./M_PI)); } \n\n", fip);
+fputs("/* Longitude(0..360), Latitude(-90..90) conversion to Molleweide X,Y   */\n", fip);
+fputs("double longlat2mollX(double lng, double lat) { \n", fip);
+fputs(" if (lng<=180.) return(180.-sin(acos(lat/90.))*(180.-lng)); \n", fip);
+fputs(" else return(180.+sin(acos(lat/90.))*(lng-180.)); } \n", fip);
+fputs("double longlat2mollY(double lat) { return lat; } \n\n\n", fip);
+
 fputs("/* NTupleInterface Variable declaration - Generated by piapp    */\n", fip);
 fputs("/*    -- Services2NObjMgr::PrepareNTExpressionCFile()  --       */ \n", fip);
@@ -1356 +1371 @@
 fputs("#define gaurand() NorRand() \n\n", fip);
 
+fputs("/* Angle conversions   */\n", fip);
+fputs("double deg2rad(double x) { return (x*M_PI/180.); } \n", fip);
+fputs("double rad2deg(double x) { return (x*180./M_PI); } \n", fip);
+fputs("/* Theta(0..Pi, phi(0..2Pi) conversion to Molleweide X,Y   */\n", fip);
+fputs("double tetphi2mollX(double teta, double phi) { \n", fip);
+fputs(" if (phi<=M_PI) return(180.-sin(acos(2*teta/M_PI-1.))*((M_PI-phi)*180./M_PI)); \n", fip);
+fputs(" else return(180.+sin(acos(2.*teta/M_PI-1.))*((phi-M_PI)*180./M_PI)); } \n", fip);
+fputs("double tetphi2mollY(double teta) { \n", fip);
+fputs(" return (90.-(teta*180./M_PI)); } \n\n", fip);
+fputs("/* Longitude(0..360), Latitude(-90..90) conversion to Molleweide X,Y   */\n", fip);
+fputs("double longlat2mollX(double lng, double lat) { \n", fip);
+fputs(" if (lng<=180.) return(180.-sin(acos(lat/90.))*(180.-lng)); \n", fip);
+fputs(" else return(180.+sin(acos(lat/90.))*(lng-180.)); } \n", fip);
+fputs("double longlat2mollY(double lat) { return lat; } \n\n\n", fip);
+
 fputs("int expf_pia_dl_func(int_8_exprf _nl, double* _zz6qi_, double* _rx_6q_, double* _ry_6q_, double* _rz_6q_, double* _rt_6q_) \n{\n", fip);
 fprintf(fip,"%s \n", vardec.c_str());