Changeset 4521 for trunk/MagicSoft/Mars/mastro

Timestamp:

08/06/04 10:13:59 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mastro

Files:

• MAstroCatalog.cc (modified) (9 diffs)
• MAstroCatalog.h (modified) (1 diff)

trunk/MagicSoft/Mars/mastro/MAstroCatalog.cc

@@ -287 +287 @@
 Int_t MAstroCatalog::ReadXephem(TString catalog)
 {
-    SetBit(kHasChanged);
-
     gLog << inf << "Reading Xephem catalog: " << catalog << endl;
 
@@ -354 +352 @@
         dec0 *= TMath::Pi()/180;
 
-        MVector3 *star=new MVector3;
-        star->SetRaDec(ra0, dec0, mag);
-        star->SetName(name);
-
-        if (star->Angle(fRaDec)*TMath::RadToDeg()>fRadiusFOV)
-        {
-            delete star;
-            continue;
-        }
-
-        fList.Add(star);
-        add++;
+        if (AddObject(ra0, dec0, mag, name))
+            add++;
     }
     gLog << inf << "Read " << add << " out of " << cnt << " (Total max mag=" << maxmag << ")" << endl;
@@ -380 +368 @@
 Int_t MAstroCatalog::ReadNGC2000(TString catalog)
 {
-    SetBit(kHasChanged);
-
     gLog << inf << "Reading NGC2000 catalog: " << catalog << endl;
 
@@ -429 +415 @@
         const Double_t dec = MAstro::Dms2Rad(decd, decm, 0, decs);
 
-        MVector3 *star=new MVector3;
-        star->SetRaDec(ra, dec, mag);
-        star->SetName(row(0, 8));
-        if (star->Angle(fRaDec)*TMath::RadToDeg()>fRadiusFOV)
-        {
-            delete star;
-            continue;
-        }
-
-        fList.Add(star);
-        add++;
+        if (AddObject(ra, dec, mag, row(0,8)))
+            add++;
     }
 
@@ -455 +432 @@
 Int_t MAstroCatalog::ReadBSC(TString catalog)
 {
-    SetBit(kHasChanged);
-
     gLog << inf << "Reading Bright Star Catalog (BSC5) catalog: " << catalog << endl;
 
@@ -506 +481 @@
         const Double_t dec = MAstro::Dms2Rad(decd, decm, decs, decsgn);
 
-        MVector3 *star=new MVector3;
-        star->SetRaDec(ra, dec, mag);
-        star->SetName(row(4,9));
-        if (star->Angle(fRaDec)*TMath::RadToDeg()>fRadiusFOV)
-        {
-            delete star;
-            continue;
-        }
-
-        fList.Add(star);
-        add++;
+        if (AddObject(ra, dec, mag, row(4,9)))
+            add++;
     }
 
@@ -538 +504 @@
 Int_t MAstroCatalog::ReadHeasarcPPM(TString catalog, TString outname)
 {
-    SetBit(kHasChanged);
-
     gLog << inf << "Reading Heasarc PPM catalog: " << catalog << endl;
 
@@ -620 +584 @@
         const Double_t de0 = MAstro::Dms2Rad(d, m, s, sgn);
 
-        MVector3 *star=new MVector3;
-        star->SetRaDec(ra0, de0, mag);
-        star->SetName(name);
-        if (star->Angle(fRaDec)*TMath::RadToDeg()>fRadiusFOV)
-        {
-            delete star;
-            continue;
-        }
-
-        fList.Add(star);
+        if (!AddObject(ra0, de0, mag, name))
+            continue;
+
         add++;
 
         if (fout)
-            star->WriteBinary(*fout);
+            ((MVector3*)fList.Last())->WriteBinary(*fout);
     }
 
@@ -698 +655 @@
 // --------------------------------------------------------------------------
 //
-// Add an object to the star catalog manually.
-//
-Int_t MAstroCatalog::AddObject(Float_t ra, Float_t dec, Float_t mag, TString name)
-{
-    SetBit(kHasChanged);
-
-    MVector3 *star = new MVector3();
+// Add an object to the star catalog manually. Return true if the object
+// was added and false otherwise (criteria is the FOV)
+//
+Bool_t MAstroCatalog::AddObject(Float_t ra, Float_t dec, Float_t mag, TString name)
+{
+    MVector3 *star = new MVector3;
     star->SetRaDec(ra, dec, mag);
     star->SetName(name);
 
-    if (star->Angle(fRaDec)*TMath::RadToDeg()<=fRadiusFOV) {
-      fList.Add(star);
-    } else {
-      delete(star);
-    }
-    
+    if (star->Angle(fRaDec)*TMath::RadToDeg()>fRadiusFOV)
+    {
+        delete star;
+        return 0;
+    }
+
+    SetBit(kHasChanged);
+    fList.AddLast(star);
     return 1;
 }

trunk/MagicSoft/Mars/mastro/MAstroCatalog.h

@@ -175 +175 @@
     Int_t ReadHeasarcPPM(TString catalog = "heasarc_ppm.tdat", TString fout="");
     Int_t ReadCompressed(TString catalog);
-    Int_t AddObject(Float_t ra, Float_t dec, Float_t mag, TString name="");
+    Bool_t AddObject(Float_t ra, Float_t dec, Float_t mag, TString name="");
 
     void Print(Option_t *o="") const { fList.Print(); } // Print all stars