// Class AliAstrolab
// Virtual lab to correlate measurements with astrophysical phenomena.
//
+// This class is derived from TTask, but the only reason for this
+// is to enable this class to serve as a base class for other TTask
+// derived classes (e.g. AliEventSelector) without the need for
+// multiple virtual inheritance.
+// So, this AliAstrolab class itself does not provide any TTask
+// related functionality.
+//
// The lab can be given a terrestrial location via the usual longitude
// and latitude specifications.
// Since this class is derived from AliTimestamp, a lab can also be
// so-called reference signals.
// This class provides facilities (e.g. MatchRefSignal) to check
// correlations of the stored measurement with these reference signals.
+// Also graphical facilities (e.g. DisplaySignals) are available to
+// provide skymaps in various projections.
//
// Coding example :
// ----------------
// Float_t pos[3]={1,90.-b,l};
// r.SetVector(pos,"sph","deg");
// lab.SetUT(1989,7,30,8,14,16,0,0);
-// lab->SetSignal(&r,"gal","M",0,-1,"GRB890730");
+// lab.SetSignal(&r,"gal","M",0,-1,"GRB890730");
//
// // List all stored objects
// lab.ListSignals("equ","M",5);
// }
// }
//
+// // Display all stored objects in a skymap (Hammer projection)
+// lab.DisplaySignals("equ","M",&ts,"ham",1);
//
//--- Author: Nick van Eijndhoven 15-mar-2007 Utrecht University
//- Modified: NvE $Date$ Utrecht University
///////////////////////////////////////////////////////////////////////////
+#include <cstdlib>
#include "AliAstrolab.h"
#include "Riostream.h"
ClassImp(AliAstrolab) // Class implementation to enable ROOT I/O
-AliAstrolab::AliAstrolab(const char* name,const char* title) : TNamed(name,title),AliTimestamp()
+AliAstrolab::AliAstrolab(const char* name,const char* title) : TTask(name,title),AliTimestamp()
{
// Default constructor
fBias=0;
fGal=0;
fIndices=0;
+ fMeridian=-999;
+ fProj="none";
+ fCanvas=0;
+ fHist=0;
+ fMarkers=0;
}
///////////////////////////////////////////////////////////////////////////
AliAstrolab::~AliAstrolab()
delete fIndices;
fIndices=0;
}
+ if (fHist)
+ {
+ delete fHist;
+ fHist=0;
+ }
+ if (fMarkers)
+ {
+ delete fMarkers;
+ fMarkers=0;
+ }
+ if (fCanvas)
+ {
+ delete fCanvas;
+ fCanvas=0;
+ }
}
///////////////////////////////////////////////////////////////////////////
-AliAstrolab::AliAstrolab(const AliAstrolab& t) : TNamed(t),AliTimestamp(t)
+AliAstrolab::AliAstrolab(const AliAstrolab& t) : TTask(t),AliTimestamp(t)
{
// Copy constructor
fBias=0;
fGal=0;
fIndices=0;
+ fMeridian=-999;
+ fProj="none";
+ fCanvas=0;
+ fHist=0;
+ fMarkers=0;
}
///////////////////////////////////////////////////////////////////////////
void AliAstrolab::Data(Int_t mode,TString u)
//
// frame = "equ" ==> Equatorial coordinates with right ascension (a) and declination (d),
// where the "sph" components of r correspond to theta=(pi/2)-d and phi=a.
-// "gal" ==> Galactic coordinates with longitude (l) and lattitude (b).
+// "gal" ==> Galactic coordinates with longitude (l) and latitude (b).
// where the "sph" components of r correspond to theta=(pi/2)-b and phi=l.
-// "ecl" ==> Ecliptic coordinates with longitude (l) and lattitude (b),
+// "ecl" ==> Ecliptic coordinates with longitude (l) and latitude (b),
// where the "sph" components of r correspond to theta=(pi/2)-b and phi=l.
// "hor" ==> Horizontal coordinates at the AliAstrolab location, where the "sph"
// components of r correspond to theta=zenith angle and phi=pi-azimuth.
-// "icr" ==> ICRS coordinates with longitude (l) and lattitude (b),
+// "icr" ==> ICRS coordinates with longitude (l) and latitude (b),
// where the "sph" components of r correspond to theta=(pi/2)-b and phi=l.
// "loc" ==> Local coordinates at the AliAstrolab location, where the "sph"
// components of r correspond to the usual theta and phi angles.
//
// frame = "equ" ==> Equatorial coordinates with right ascension (a) and declination (d),
// where the "sph" components of r correspond to theta=(pi/2)-d and phi=a.
-// "gal" ==> Galactic coordinates with longitude (l) and lattitude (b).
+// "gal" ==> Galactic coordinates with longitude (l) and latitude (b).
// where the "sph" components of r correspond to theta=(pi/2)-b and phi=l.
-// "ecl" ==> Ecliptic coordinates with longitude (l) and lattitude (b),
+// "ecl" ==> Ecliptic coordinates with longitude (l) and latitude (b),
// where the "sph" components of r correspond to theta=(pi/2)-b and phi=l.
// "hor" ==> Horizontal coordinates at the AliAstrolab location, where the "sph"
// components of r correspond to theta=zenith angle and phi=pi-azimuth.
-// "icr" ==> ICRS coordinates with longitude (l) and lattitude (b),
+// "icr" ==> ICRS coordinates with longitude (l) and latitude (b),
// where the "sph" components of r correspond to theta=(pi/2)-b and phi=l.
// "loc" ==> Local coordinates at the AliAstrolab location, where the "sph"
// components of r correspond to the usual theta and phi angles.
//
// frame = "equ" ==> Equatorial coordinates with right ascension (a) and declination (d),
// where the "sph" components of r correspond to theta=(pi/2)-d and phi=a.
-// "gal" ==> Galactic coordinates with longitude (l) and lattitude (b).
+// "gal" ==> Galactic coordinates with longitude (l) and latitude (b).
// where the "sph" components of r correspond to theta=(pi/2)-b and phi=l.
-// "ecl" ==> Ecliptic coordinates with longitude (l) and lattitude (b),
+// "ecl" ==> Ecliptic coordinates with longitude (l) and latitude (b),
// where the "sph" components of r correspond to theta=(pi/2)-b and phi=l.
// "hor" ==> Horizontal coordinates at the AliAstrolab location, where the "sph"
// components of r correspond to theta=zenith angle and phi=pi-azimuth.
-// "icr" ==> ICRS coordinates with longitude (l) and lattitude (b),
+// "icr" ==> ICRS coordinates with longitude (l) and latitude (b),
// where the "sph" components of r correspond to theta=(pi/2)-b and phi=l.
// "loc" ==> Local coordinates at the AliAstrolab location, where the "sph"
// components of r correspond to the usual theta and phi angles.
//
// frame = "equ" ==> Equatorial coordinates with right ascension (a) and declination (d).
//
-// "gal" ==> Galactic coordinates with longitude (l) and lattitude (b).
+// "gal" ==> Galactic coordinates with longitude (l) and latitude (b).
//
-// "ecl" ==> Ecliptic coordinates with longitude (l) and lattitude (b).
+// "ecl" ==> Ecliptic coordinates with longitude (l) and latitude (b).
//
// "hor" ==> Horizontal azimuth and altitude coordinates at the AliAstrolab location.
//
-// "icr" ==> ICRS coordinates with longitude (l) and lattitude (b).
+// "icr" ==> ICRS coordinates with longitude (l) and latitude (b).
//
// "loc" ==> Local spherical angles theta and phi at the AliAstrolab location.
//
// mode = "M" --> Input coordinates are the mean values
// "T" --> Input coordinates are the true values
//
+// The input parameter "mode" also determines which type of time and
+// local hour angle will appear in the printout.
+//
+// mode = "M" --> Mean Sidereal Time (MST) and Local Mean Hour Angle (LMHA)
+// "T" --> Apparent Sidereal Time (AST) and Local Apparent Hour Angle (LAHA)
+//
// The input parameter "jref" allows printing of a so-called "reference" signal.
// These reference signals may serve to check space-time event coincidences with the
// stored measurement (e.g. coincidence of the measurement with transient phenomena).
if (!sx) return;
+ // Local Hour Angle of the signal
+ Double_t lha=GetHourAngle("A",ts,jref);
+ TString slha="LAHA";
+ if (mode=="M" || mode=="m")
+ {
+ lha=GetHourAngle("M",ts,jref);
+ slha="LMHA";
+ }
+
TString name=sx->GetName();
if (name != "") cout << name.Data() << " ";
a=r.GetX(3,"sph","rad");
cout << "Equatorial (" << mode.Data() <<") a : "; PrintAngle(a,"rad","hms",ndig);
cout << " d : "; PrintAngle(d,"deg","dms",ndig);
+ cout << " " << slha.Data() << " : "; PrintAngle(lha,"deg","hms",ndig);
return;
}
l=r.GetX(3,"sph","deg");
cout << "Galactic l : "; PrintAngle(l,"deg","deg",ndig);
cout << " b : "; PrintAngle(b,"deg","deg",ndig);
+ cout << " " << slha.Data() << " : "; PrintAngle(lha,"deg","hms",ndig);
return;
}
a=r.GetX(3,"sph","rad");
cout << "ICRS l : "; PrintAngle(a,"rad","hms",ndig);
cout << " b : "; PrintAngle(d,"deg","dms",ndig);
+ cout << " " << slha.Data() << " : "; PrintAngle(lha,"deg","hms",ndig);
return;
}
a=r.GetX(3,"sph","deg");
cout << "Ecliptic l : "; PrintAngle(a,"deg","deg",ndig);
cout << " b : "; PrintAngle(d,"deg","deg",ndig);
+ cout << " " << slha.Data() << " : "; PrintAngle(lha,"deg","hms",ndig);
return;
}
}
cout << "Horizontal azi : "; PrintAngle(azi,"deg","deg",ndig);
cout << " alt : "; PrintAngle(alt,"deg","deg",ndig);
+ cout << " " << slha.Data() << " : "; PrintAngle(lha,"deg","hms",ndig);
return;
}
Double_t phi=r.GetX(3,"sph","deg");
cout << "Local-frame phi : "; PrintAngle(phi,"deg","deg",ndig);
cout << " theta : "; PrintAngle(theta,"deg","deg",ndig);
+ cout << " " << slha.Data() << " : "; PrintAngle(lha,"deg","hms",ndig);
return;
}
}
{
// Print data of the stored signal with the specified name in user specified coordinates
// at the specific timestamp ts.
-// In case such stored signal was available or one of the input arguments was
+// In case no such stored signal was available or one of the input arguments was
// invalid, no printout is produced.
//
// The argument "ndig" specifies the number of digits for the fractional
//
// frame = "equ" ==> Equatorial coordinates with right ascension (a) and declination (d).
//
-// "gal" ==> Galactic coordinates with longitude (l) and lattitude (b).
+// "gal" ==> Galactic coordinates with longitude (l) and latitude (b).
//
-// "ecl" ==> Ecliptic coordinates with longitude (l) and lattitude (b).
+// "ecl" ==> Ecliptic coordinates with longitude (l) and latitude (b).
//
// "hor" ==> Horizontal azimuth and altitude coordinates at the AliAstrolab location.
//
-// "icr" ==> ICRS coordinates with longitude (l) and lattitude (b).
+// "icr" ==> ICRS coordinates with longitude (l) and latitude (b).
//
// "loc" ==> Local spherical angles theta and phi at the AliAstrolab location.
//
// mode = "M" --> Input coordinates are the mean values
// "T" --> Input coordinates are the true values
//
+// The input parameter "mode" also determines which type of time and
+// local hour angle will appear in the printout.
+//
+// mode = "M" --> Mean Sidereal Time (MST) and Local Mean Hour Angle (LMHA)
+// "T" --> Apparent Sidereal Time (AST) and Local Apparent Hour Angle (LAHA)
+//
// Note : In case ts=0 the current timestamp of the lab will be taken.
Int_t j=GetSignalIndex(name);
// will appear as 03.4 on the output.
// Due to computer accuracy, precision on the pico-arcsecond level may get lost.
//
-// The default value is ndig=1;
+// The default value is ndig=1.
//
// Note : The angle info is printed without additional spaces or "endline".
// This allows the print to be included in various composite output formats.
//
// frame = "equ" ==> Equatorial coordinates with right ascension (a) and declination (d).
//
-// "gal" ==> Galactic coordinates with longitude (l) and lattitude (b).
+// "gal" ==> Galactic coordinates with longitude (l) and latitude (b).
//
-// "ecl" ==> Ecliptic coordinates with longitude (l) and lattitude (b).
+// "ecl" ==> Ecliptic coordinates with longitude (l) and latitude (b).
//
// "hor" ==> Horizontal azimuth and altitude coordinates at the AliAstrolab location.
//
-// "icr" ==> ICRS coordinates with longitude (l) and lattitude (b).
+// "icr" ==> ICRS coordinates with longitude (l) and latitude (b).
//
// "loc" ==> Local spherical angles theta and phi at the AliAstrolab location.
//
// mode = "M" --> Input coordinates are the mean values
// "T" --> Input coordinates are the true values
//
-// Note : In case ts=0 the current timestamp of the lab will be taken.
+// The input parameter "mode" also determines which type of time and
+// local hour angle will appear in the listing.
+//
+// mode = "M" --> Mean Sidereal Time (MST) and Local Mean Hour Angle (LMHA)
+// "T" --> Apparent Sidereal Time (AST) and Local Apparent Hour Angle (LAHA)
Int_t iprint=0;
AliTimestamp* tx=0;
+ Int_t dform=1;
+ if (mode=="T" || mode=="t") dform=-1;
+
if (fXsig)
{
- AliTimestamp* tx=fXsig->GetTimestamp();
+ tx=fXsig->GetTimestamp();
if (!tx) tx=(AliTimestamp*)this;
cout << " *AliAstrolab::ListSignals* List of all stored signals." << endl;
- cout << " The measurement under investigation." << endl;
- cout << " --- Timestamp of the actual observation ---" << endl;
- cout << " "; tx->Date(1);
- cout << " --- Location of the actual observation ---" << endl;
- cout << " "; PrintSignal(frame,mode,tx,ndig); cout << endl;
+ cout << " === The measurement under investigation ===" << endl;
+ cout << " Timestamp of the actual observation" << endl;
+ tx->Date(dform,fToffset);
+ cout << " Location of the actual observation" << endl;
+ cout << " "; PrintSignal(frame,mode,tx,ndig); cout << endl;
iprint=1;
}
{
cout << " *AliAstrolab::ListRefSignals* List of all stored signals." << endl;
tx=(AliTimestamp*)this;
- cout << " --- Current timestamp of the laboratory ---" << endl;
- cout << " "; tx->Date(1);
+ cout << " Current timestamp of the laboratory" << endl;
+ tx->Date(dform,fToffset);
iprint=1;
}
if (iprint==1)
{
- cout << " --- All stored reference signals according to the above timestamp ---" << endl;
+ cout << " === All stored reference signals according to the above timestamp ===" << endl;
iprint=2;
}
- cout << " Index : " << i << " "; PrintSignal(frame,mode,tx,ndig,i); cout << endl;
+ cout << " Index : " << i << " "; PrintSignal(frame,mode,tx,ndig,i); cout << endl;
}
}
///////////////////////////////////////////////////////////////////////////
// Specification of the orientations of the local-frame axes.
// The input arguments represent the angles (in degrees) of the local-frame axes
// w.r.t. a so called Master Reference Frame (MRF), with the same convention
-// as the input arguments of TRrotMatix::SetAngles.
+// as the input arguments of TRotMatrix::SetAngles.
//
// The right handed Master Reference Frame is defined as follows :
// Z-axis : Points to the local Zenith
// "deg" : input angle provided in degrees
// "dms" : input angle provided in dddmmss.sss
// "hms" : input angle provided in hhmmss.sss
+// "hrs" : input angle provided in fractional hours
//
// out = "rad" : output angle provided in radians
// "deg" : output angle provided in degrees
// "dms" : output angle provided in dddmmss.sss
// "hms" : output angle provided in hhmmss.sss
+// "hrs" : output angle provided in fractional hours
if (in==out) return a;
if (in=="rad") b*=180./pi;
+ if (in=="hrs") b*=15.;
+
if (in=="dms")
{
word=Int_t(b);
if (out=="rad") b*=pi/180.;
+ if (out=="hrs") b/=15.;
+
if (out=="dms")
{
ddd=Int_t(b);
//
// The input parameter "jref" allows the user to specify so-called "reference" signals.
// These reference signals may be used to check space-time event coincidences with the
-// stored regular signal (e.g. coincidence of the stored signal with transient phenomena).
+// stored measurement (e.g. coincidence of the measurement with transient phenomena).
//
-// jref = 0 --> Use the regular signal
+// jref = 0 --> Use the stored measurement
// j --> Use the reference signal at the j-th position (j=1 is first)
//
// Default value is jref=0.
// Get corrected right ascension and declination for the specified timestamp.
Double_t a,d;
- GetSignal(a,d,"T",ts,jref);
+ if (mode=="M" || mode=="m") GetSignal(a,d,"M",ts,jref);
+ if (mode=="A" || mode=="a") GetSignal(a,d,"T",ts,jref);
// Convert coordinates to fractional degrees.
a=ConvertAngle(a,"hms","deg");
return fIndices;
}
///////////////////////////////////////////////////////////////////////////
+void AliAstrolab::DisplaySignal(TString frame,TString mode,AliTimestamp* ts,Int_t jref,TString proj,Int_t clr)
+{
+// Display a stored signal in a user specified coordinate projection
+// at the specific timestamp ts.
+//
+// In case no stored signal was available or one of the input arguments was
+// invalid, no display is produced.
+//
+// Note : In case ts=0 the current timestamp of the lab will be taken.
+//
+// The input parameter "frame" allows the user to specify the frame to which
+// the coordinates refer. Available options are :
+//
+// frame = "equ" ==> Equatorial coordinates with right ascension (a) and declination (d).
+//
+// "gal" ==> Galactic coordinates with longitude (l) and latitude (b).
+//
+// "ecl" ==> Ecliptic coordinates with longitude (l) and latitude (b).
+//
+// "hor" ==> Horizontal azimuth and altitude coordinates at the AliAstrolab location.
+//
+// "icr" ==> ICRS coordinates with longitude (l) and latitude (b).
+//
+// "loc" ==> Local spherical angles theta and phi at the AliAstrolab location.
+//
+// In case the coordinates are the equatorial right ascension and declination (a,d),
+// they can represent so-called "mean" and "true" values.
+// The distinction between these two representations is the following :
+//
+// mean values : (a,d) are only corrected for precession and not for nutation
+// true values : (a,d) are corrected for both precession and nutation
+//
+// The input parameter "mode" allows the user to specifiy either "mean" or "true"
+// values for the input in case of equatorial (a,d) coordinates.
+//
+// mode = "M" --> Input coordinates are the mean values
+// "T" --> Input coordinates are the true values
+//
+// The input parameter "jref" allows display of a so-called "reference" signal.
+// These reference signals may serve to check space-time event coincidences with the
+// stored measurement (e.g. coincidence of the measurement with transient phenomena).
+//
+// jref = 0 --> Display of the measurement
+// j --> Display of the j-th reference signal
+//
+// Default value is jref=0.
+//
+// The input parameter "proj" allows the user to specify the desired projection.
+// The available projection modes are :
+//
+// cyl : Cylindrical projection plotted with colored markers
+// cylh : Cylindrical projection plotted in a 2-D histogram
+// ham : Hammer equal area projection plotted with colored markers
+// hamh : Hammer equal area projection plotted in a 2-D histogram
+// ait : Aitoff projection plotted with colored markers
+// aith : Aitoff projection plotted in a 2-D histogram
+// mer : Mercator projection plotted with colored markers
+// merh : Mercator projection plotted in a 2-D histogram
+// ang : Straight cos(b) vs. l plot with colored markers
+// angh : Straight cos(b) vs. l plot in a 2-D histogram
+//
+// Note : The ang(h) plot allows for easy identification of a homogeneous distribution.
+//
+// The input argument "clr" allows to clear (1) the display before drawing or not (0).
+//
+// The default values are : jref=0, proj="ham" and clr=0.
+//
+// This routine is based on the work by Garmt de Vries-Uiterweerd.
+
+ Ali3Vector r;
+ AliSignal* sx=GetSignal(r,frame,mode,ts,jref);
+
+ if (!sx) return;
+
+ // The generic input angles (in rad) for the projections
+ Double_t theta=0;
+ Double_t phi=0;
+
+ Double_t pi=acos(-1.);
+
+ if (frame=="equ" || frame=="gal" || frame=="icr" || frame=="ecl" || frame=="loc")
+ {
+ theta=(pi/2.)-r.GetX(2,"sph","rad");
+ phi=r.GetX(3,"sph","rad");
+ }
+
+ if (frame=="hor")
+ {
+ theta=(pi/2.)-r.GetX(2,"sph","rad");
+ phi=pi-r.GetX(3,"sph","rad");
+ }
+
+ // Automatic choice of central meridian if not selected by the user
+ if (fMeridian<-pi)
+ {
+ if (frame=="gal")
+ {
+ fMeridian=0;
+ }
+ else
+ {
+ fMeridian=pi;
+ }
+ }
+
+ // Obtain the projected (x,y) position
+ Double_t x=0;
+ Double_t y=0;
+ Project(phi,theta,proj,x,y);
+
+ Int_t hist=0;
+ if (proj=="hamh" || proj=="aith" || proj=="merh" || proj=="cylh" || proj=="angh") hist=1;
+
+ // Update the display for this signal position
+
+ // Create a new canvas if needed
+ if (!fCanvas) fCanvas=new TCanvas("AliAstrolab","Skymap");
+
+ // Construct the title string for this map
+ TString title="Projection : ";
+ title+=frame;
+ title+=" ";
+ title+=proj;
+ title+=" Center : ";
+ Int_t ang,h,m,s,d;
+ if (frame=="equ" || frame=="icr")
+ {
+ ang=int(ConvertAngle(fMeridian,"rad","hms"));
+ h=ang/10000;
+ ang=ang%10000;
+ m=ang/100;
+ s=ang%100;
+ title+=h;
+ title+="h ";
+ title+=m;
+ title+="m ";
+ title+=s;
+ title+="s";
+ }
+ else
+ {
+ ang=int(ConvertAngle(fMeridian,"rad","dms"));
+ d=ang/10000;
+ ang=ang%10000;
+ m=ang/100;
+ s=ang%100;
+ title+=d;
+ title+="d ";
+ title+=m;
+ title+="' ";
+ title+=s;
+ title+="\"";
+ }
+
+ if (!hist) // 2-D Marker display (i.e. not a histogram)
+ {
+ // Remove existing markers, grid and outline from display if needed
+ if (clr==1 || proj!=fProj)
+ {
+ if (fMarkers)
+ {
+ delete fMarkers;
+ fMarkers=0;
+ }
+ fCanvas->Clear();
+ fProj=proj;
+ }
+
+ if (!fMarkers)
+ {
+ fMarkers=new TObjArray();
+ fMarkers->SetOwner();
+ // Set canvas range, header and axes
+ if (proj=="mer")
+ {
+ fCanvas->Range(-2.2,-6.,2.2,6.);
+ TText* header=new TText;
+ fMarkers->Add(header);
+ header->DrawText(-1.4,5.5,title.Data());
+ TLine* line=new TLine(-2,0,2,0);
+ fMarkers->Add(line);
+ line->Draw();
+ line=new TLine(0,5.5,0,-5.5);
+ fMarkers->Add(line);
+ line->Draw();
+ }
+ else
+ {
+ fCanvas->Range(-2.2,-1.2,2.2,1.2);
+ TText* header=new TText;
+ fMarkers->Add(header);
+ header->DrawText(-1.4,1.1,title.Data());
+ TLine* line=new TLine(-2,0,2,0);
+ fMarkers->Add(line);
+ line->Draw();
+ line=new TLine(0,1,0,-1);
+ fMarkers->Add(line);
+ line->Draw();
+ }
+ // Draw the outline
+ if (proj=="ham" || proj=="ait")
+ {
+ // Draw ellips outline with axes
+ TEllipse* outline=new TEllipse(0,0,2,1);
+ fMarkers->Add(outline);
+ outline->Draw();
+ }
+ }
+ TMarker* marker=0;
+ if (!jref) // The measurement
+ {
+ marker=new TMarker(x,y,29);
+ marker->SetMarkerColor(kRed);
+ }
+ else // The reference signal(s)
+ {
+ marker=new TMarker(x,y,20);
+ marker->SetMarkerColor(kBlue);
+ }
+ marker->SetMarkerSize(1);
+ fMarkers->Add(marker);
+ marker->Draw();
+ }
+ else if (hist==1) // 2-D display via histogram
+ {
+ // Reset the histogram if needed
+ if (clr==1 || proj!=fProj)
+ {
+ fCanvas->Clear();
+ if (!fHist) fHist=new TH2F();
+ fHist->Reset();
+ fHist->SetMarkerStyle(20);
+ fHist->SetMarkerSize(1);
+ fHist->SetMarkerColor(kBlue);
+ fHist->SetNameTitle("SkyMap",title.Data());
+ if (proj=="merh")
+ {
+ fHist->SetBins(100,-2.2,2.2,100,-5.5,5.5);
+ }
+ else
+ {
+ fHist->SetBins(100,-2.2,2.2,100,-1.1,1.1);
+ }
+ fProj=proj;
+ }
+ fHist->Fill(x,y);
+ fHist->Draw();
+ }
+}
+///////////////////////////////////////////////////////////////////////////
+void AliAstrolab::DisplaySignal(TString frame,TString mode,AliTimestamp* ts,TString name,TString proj,Int_t clr)
+{
+// Display the stored signal with the specified name in a user specified
+// coordinate projection at the specific timestamp ts.
+//
+// Note : In case ts=0 the current timestamp of the lab will be taken.
+//
+// In case no such stored signal was available or one of the input arguments was
+// invalid, no display is produced.
+//
+// The input parameter "frame" allows the user to specify the frame to which
+// the coordinates refer. Available options are :
+//
+// frame = "equ" ==> Equatorial coordinates with right ascension (a) and declination (d).
+//
+// "gal" ==> Galactic coordinates with longitude (l) and latitude (b).
+//
+// "ecl" ==> Ecliptic coordinates with longitude (l) and latitude (b).
+//
+// "hor" ==> Horizontal azimuth and altitude coordinates at the AliAstrolab location.
+//
+// "icr" ==> ICRS coordinates with longitude (l) and latitude (b).
+//
+// "loc" ==> Local spherical angles theta and phi at the AliAstrolab location.
+//
+// In case the coordinates are the equatorial right ascension and declination (a,d),
+// they can represent so-called "mean" and "true" values.
+// The distinction between these two representations is the following :
+//
+// mean values : (a,d) are only corrected for precession and not for nutation
+// true values : (a,d) are corrected for both precession and nutation
+//
+// The input parameter "mode" allows the user to specifiy either "mean" or "true"
+// values for the input in case of equatorial (a,d) coordinates.
+//
+// mode = "M" --> Input coordinates are the mean values
+// "T" --> Input coordinates are the true values
+//
+// The input parameter "proj" allows the user to specify the desired projection.
+// The available projection modes are :
+//
+// cyl : Cylindrical projection plotted with colored markers
+// cylh : Cylindrical projection plotted in a 2-D histogram
+// ham : Hammer equal area projection plotted with colored markers
+// hamh : Hammer equal area projection plotted in a 2-D histogram
+// ait : Aitoff projection plotted with colored markers
+// aith : Aitoff projection plotted in a 2-D histogram
+// mer : Mercator projection plotted with colored markers
+// merh : Mercator projection plotted in a 2-D histogram
+// ang : Straight cos(b) vs. l plot with colored markers
+// angh : Straight cos(b) vs. l plot in a 2-D histogram
+//
+// Note : The ang(h) plot allows for easy identification of a homogeneous distribution.
+//
+// The input argument "clr" allows to clear (1) the display before drawing or not (0).
+//
+// The default values are : proj="ham" and clr=0.
+//
+// This routine is based on the work by Garmt de Vries-Uiterweerd.
+
+ Int_t j=GetSignalIndex(name);
+ if (j>=0) DisplaySignal(frame,mode,ts,j,proj,clr);
+}
+///////////////////////////////////////////////////////////////////////////
+void AliAstrolab::DisplaySignals(TString frame,TString mode,AliTimestamp* ts,TString proj,Int_t clr)
+{
+// Display of all stored signals in user specified coordinate projection
+// at the specific timestamp ts.
+// In case ts=0, the timestamp of the actual recording of the stored measurement
+// under investigation will be used.
+// In case no (timestamp of the) actual measurement is available,
+// the current timestamp of the lab will be taken.
+// In case no stored signal is available or one of the input arguments is
+// invalid, no display is produced.
+//
+// The input parameter "frame" allows the user to specify the frame to which
+// the coordinates refer. Available options are :
+//
+// frame = "equ" ==> Equatorial coordinates with right ascension (a) and declination (d).
+//
+// "gal" ==> Galactic coordinates with longitude (l) and latitude (b).
+//
+// "ecl" ==> Ecliptic coordinates with longitude (l) and latitude (b).
+//
+// "hor" ==> Horizontal azimuth and altitude coordinates at the AliAstrolab location.
+//
+// "icr" ==> ICRS coordinates with longitude (l) and latitude (b).
+//
+// "loc" ==> Local spherical angles theta and phi at the AliAstrolab location.
+//
+// In case the coordinates are the equatorial right ascension and declination (a,d),
+// they can represent so-called "mean" and "true" values.
+// The distinction between these two representations is the following :
+//
+// mean values : (a,d) are only corrected for precession and not for nutation
+// true values : (a,d) are corrected for both precession and nutation
+//
+// The input parameter "mode" allows the user to specifiy either "mean" or "true"
+// values for the input in case of equatorial (a,d) coordinates.
+//
+// mode = "M" --> Input coordinates are the mean values
+// "T" --> Input coordinates are the true values
+//
+// The input parameter "proj" allows the user to specify the desired projection.
+// The available projection modes are :
+//
+// cyl : Cylindrical projection plotted with colored markers
+// cylh : Cylindrical projection plotted in a 2-D histogram
+// ham : Hammer equal area projection plotted with colored markers
+// hamh : Hammer equal area projection plotted in a 2-D histogram
+// ait : Aitoff projection plotted with colored markers
+// aith : Aitoff projection plotted in a 2-D histogram
+// mer : Mercator projection plotted with colored markers
+// merh : Mercator projection plotted in a 2-D histogram
+// ang : Straight cos(b) vs. l plot with colored markers
+// angh : Straight cos(b) vs. l plot in a 2-D histogram
+//
+// Note : The ang(h) plot allows for easy identification of a homogeneous distribution.
+//
+// The input argument "clr" allows to clear (1) the display before drawing or not (0).
+//
+// The default values are : proj="ham" and clr=0.
+//
+// This routine is based on the work by Garmt de Vries-Uiterweerd.
+
+ AliTimestamp* tx=ts;
+ if (fXsig && !tx) tx=fXsig->GetTimestamp();
+ if (!tx) tx=(AliTimestamp*)this;
+
+ // Display all stored reference signals
+ if (fRefs)
+ {
+ for (Int_t i=1; i<=fRefs->GetSize(); i++)
+ {
+ DisplaySignal(frame,mode,tx,i,proj,clr);
+ clr=0; // No display clear for subsequent signals
+ }
+ }
+
+ // Display the measurement
+ DisplaySignal(frame,mode,tx,0,proj,clr);
+}
+///////////////////////////////////////////////////////////////////////////
+void AliAstrolab::SetCentralMeridian(Double_t phi,TString u)
+{
+// Set the central meridian for the sky display.
+// Setting a value smaller than -pi will induce automatic meridian setting
+// in the display.
+// By default the central meridian is set at -999 in the constructor.
+//
+// The string argument "u" allows to choose between different angular units
+// u = "rad" : angle provided in radians
+// "deg" : angle provided in degrees
+// "dms" : angle provided in dddmmss.sss
+// "hms" : angle provided in hhmmss.sss
+//
+// The default is u="deg".
+//
+// This routine is based on the work by Garmt de Vries-Uiterweerd.
+
+ fMeridian=ConvertAngle(phi,u,"rad");
+}
+///////////////////////////////////////////////////////////////////////////
+void AliAstrolab::Project(Double_t l,Double_t b,TString proj,Double_t& x,Double_t& y)
+{
+// Generic interface for projection of a (long,lat) pair onto a (x,y) pair.
+//
+// Meaning of the arguments :
+// --------------------------
+// l : Longitude (e.g. right ascension)
+// b : Latitude (e.g. declination)
+// proj : Projection mode (e.g. "ham")
+// x : The resulting x coordinate for the selected projection
+// y : The resulting x coordinate for the selected projection
+//
+// The available projection modes are :
+//
+// cyl : Cylindrical projection plotted with colored markers
+// cylh : Cylindrical projection plotted in a 2-D histogram
+// ham : Hammer equal area projection plotted with colored markers
+// hamh : Hammer equal area projection plotted in a 2-D histogram
+// ait : Aitoff projection plotted with colored markers
+// aith : Aitoff projection plotted in a 2-D histogram
+// mer : Mercator projection plotted with colored markers
+// merh : Mercator projection plotted in a 2-D histogram
+// ang : Straight cos(b) vs. l plot with colored markers
+// angh : Straight cos(b) vs. l plot in a 2-D histogram
+//
+// Note : The ang(h) plot allows for easy identification of a homogeneous distribution.
+//
+// This routine is based on the work by Garmt de Vries-Uiterweerd.
+
+ Double_t pi=acos(-1.);
+
+ // Subtract central meridian from longitude
+ l-=fMeridian;
+
+ // Take l between -180 and 180 degrees
+ while (l>pi)
+ {
+ l-=2*pi;
+ }
+ while (l<-pi)
+ {
+ l+=2*pi;
+ }
+
+ x=0;
+ y=0;
+
+ // Convert (l,b) to (x,y) with -2 < x <= 2
+ if (proj=="cyl" || proj=="cylh") ProjectCylindrical(l,b,x,y);
+ if (proj=="ham" || proj=="hamh") ProjectHammer(l,b,x,y);
+ if (proj=="ait" || proj=="aith") ProjectAitoff(l,b,x,y);
+ if (proj=="mer" || proj=="merh") ProjectMercator(l,b,x,y);
+ if (proj=="ang" || proj=="angh")
+ {
+ x=2*l/pi;
+ y=cos(b);
+ }
+}
+///////////////////////////////////////////////////////////////////////////
+void AliAstrolab::ProjectCylindrical(Double_t l,Double_t b,Double_t& x, Double_t& y)
+{
+// Cylindrical projection of a (long,lat) coordinate pair
+//
+// This routine is based on the work by Garmt de Vries-Uiterweerd.
+
+ Double_t pi=acos(-1.);
+ x=2*l/pi;
+ y=2*b/pi;
+}
+///////////////////////////////////////////////////////////////////////////
+void AliAstrolab::ProjectHammer(Double_t l,Double_t b,Double_t& x,Double_t& y)
+{
+// Hammer-Aitoff projection of a (long,lat) coordinate pair.
+// This is an equal-area projection.
+//
+// This routine is based on the work by Garmt de Vries-Uiterweerd.
+
+ Double_t k=1./sqrt(1.+cos(b)*cos(l/2.));
+ x=2.*k*cos(b)*sin(l/2.);
+ y=k*sin(b);
+}
+///////////////////////////////////////////////////////////////////////////
+void AliAstrolab::ProjectAitoff(Double_t l,Double_t b,Double_t& x,Double_t& y)
+{
+// Aitoff projection of a (long,lat) coordinate pair.
+//
+// This routine is based on the work by Garmt de Vries-Uiterweerd.
+
+ Double_t pi=acos(-1.);
+ x=0;
+ y=0;
+ Double_t k=acos(cos(b)*cos(l/2.));
+ if(sin(k)!=0)
+ {
+ x=4.*k*cos(b)*sin(l/2.)/(pi*sin(k));
+ y=2.*k*sin(b)/(pi*sin(k));
+ }
+}
+///////////////////////////////////////////////////////////////////////////
+void AliAstrolab::ProjectMercator(Double_t l,Double_t b,Double_t& x,Double_t& y)
+{
+// Mercator projection of a (long,lat) coordinate pair
+//
+// This routine is based on the work by Garmt de Vries-Uiterweerd.
+
+ Double_t pi=acos(-1.);
+ x=2.*l/pi;
+ y=0;
+ if (b >= 89.*pi/180.) b=89.*pi/180.;
+ if (b <= -89.*pi/180.) y=-89.*pi/180.;
+ if (fabs(y) < 1) y=log((1.+sin(b))/(1.-sin(b)))/2.;
+}
+///////////////////////////////////////////////////////////////////////////
TObject* AliAstrolab::Clone(const char* name) const
{
// Make a deep copy of the current object and provide the pointer to the copy.
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff