[u/mrichter/AliRoot.git] / RICH / menu.C

void Show()
{  
//How to get number of events:  
  Info("RICH/menu.C::Show","3 ways to get number of events 1=%i 2=%f 3=%f",
                          gAlice->GetEventsPerRun(),
                          al->TreeE()->GetEntries(),
                          gAlice->TreeE()->GetEntries());
  rl->LoadHits();
  for(Int_t iEventN=0;iEventN<gAlice->GetEventsPerRun();iEventN++){// loop on events
    Int_t iNparticles=gAlice->GetEvent(iEventN);
    Int_t iNtracks=rl->TreeH()->GetEntries();
    Info("RICH/menu.C::Show","Event %i contains %i particles in total while %i are primary",
                                     iEventN,    iNparticles,                iNtracks);
    for(Int_t iTrackN=0;iTrackN<iNtracks;iTrackN++){ // loop on tracks
      rl->TreeH()->GetEntry(iTrackN);
      al->Stack()->Particle(iTrackN)->Print();
      Info("RICH/menu.C::Show","track %i has %i hits",iTrackN,r->Hits()->GetEntries());
//      for(AliRICHhit *pRichHit=(AliRICHhit*)gpRich->FirstHit(-1);pRichHit;pRichHit=(AliRICHhit*)gpRich->NextHit()){
//        pRichHit->Dump();        
//              TVector3 armV3=pRICH->ToArm(mrsV3);
//              cout<<"After\n";armV3.Dump();
//      }//loop on hits of given track
    }// loop on tracks
  }//loop on events
}//void Show()

void menu(Int_t iNevents=5)// How many events to generate.
{ 
  Info("RICH/menu.C","%i event(s) are requested",iNevents);
   
  TString runString="gAlice->Run(";  runString=runString+iNevents+");";

  TControlBar *pMenu = new TControlBar("vertical","RICH main");
       
  pMenu->AddButton("Debug ON-OFF",     "Debug();",   "Switch debug on-off");   
  if(GetAlice()){//it's from file, reconstruct
    pMenu->AddButton("Show","Show()","Shows the structure of events in files");
    pMenu->AddButton("Hits2SDigits","r->Hits2SDigits()","Perform first phase converstion");
    pMenu->AddButton("RingViewer","RingViewer()","Show rings with reconstructed info");
  }else{//it's aliroot, simulate
    pMenu->AddButton("Run",           runString.Data(),       "Process!");
  }
  pMenu->AddButton("Geo",          "Geo()",                 "Geomentry submenu");
  pMenu->AddButton("Browser",      "new TBrowser;",         "Start ROOT TBrowser");
  pMenu->AddButton("Quit",         ".q",                    "Close session");
  pMenu->Show();
  a=gAlice;//for manual manipulation convinience
}//void menu(Int_t iNevents)


void Debug()
{
  if(gAlice->GetDebug()){
    Info("RICH/menu.C::Debug","OFF");
    gAlice->SetDebug(0);
    if(r)r->SetDebug(0);
    if(t)t->SetDebug(0);
    AliLoader::SetDebug(0);
  }else{
    Info("RICH/menu.C::Debug","ON");
    gAlice->SetDebug(1);
    if(r)r->SetDebug(1);
    if(t)t->SetDebug(1);
    AliLoader::SetDebug(1);
  }
}//void Debug()

AliRun *a;
AliRunLoader *al;
AliLoader *rl,*tl,*il;

AliRICH *r;
AliTPC *t;
AliITS *i;
AliPHOS *p;

Bool_t GetAlice()
{
  if(gAlice){//it's aliroot
    Info("RICH/menu.C::GetAlice","gAlice!=NULL, IT'S ALIROOT, EXECUTE Init()");
    gAlice->Init("~/my/AliceConfig.C");
    r=(AliRICH*)gAlice->GetDetector("RICH");
    t=(AliTPC*)gAlice->GetDetector("TPC");
    i=(AliITS*)gAlice->GetDetector("ITS");
    p=(AliPHOS*)gAlice->GetDetector("PHOS");    
    return kFALSE;
  }else{//it's root with ALICE libs loaded
    Info("RICH/menu.C::GetAlice","gAlice=0, GETTING IT FROM SIMULATED FILE.");
        
    if(!(al=AliRunLoader::Open("galice.root","AlicE","update"))) Fatal("RICH/menu.C::GetAlice","Can't get AliRunLoader");
    al->LoadgAlice();
    if(!gAlice) Fatal("RICH/menu.C::GetAlice","No gAlice in file");
//    a=al->GetAliRun();    
    al->LoadHeader();//loads events tree
    al->LoadKinematics();//loads the primaries info
//RICH      
    if(!(r=(AliRICH*)gAlice->GetDetector("RICH"))) Warning("RICH/menu.C::GetAlice","No RICH in file");
    if(!(rl=al->GetLoader("RICHLoader")))          Warning("RICH/menu.C::GetAlice","No RICH loader in file");        
//TPC            
    if(!(t=(AliTPC*)gAlice->GetDetector("TPC")))   Warning("RICH/menu.C::GetAlice","No TPC in file");
    if(!(tl=al->GetLoader("TPCLoader")))           Warning("RICH/menu.C::GetAlice","No TPC loader in file");    
        
    Info("RICH/menu.C::GetAlice","Run contains %i event(s)",gAlice->GetEventsPerRun());    
    return kTRUE;
  }       
}//void GetAlice()         

void RingViewer()
{
  gStyle->SetPalette(1);
  TCanvas *view=new TCanvas("Display","ALICE RICH Display",0,0,1200,750);
  
  TH2F *pH2=new TH2F("pH2F","RICH DISPLAY",160,0,160,144,0,144);
  pH2->SetStats(0);
  pH2->SetMaximum(100);

  Int_t Nevents = gAlice->GetEventsPerRun();
}

//______________________________________________________________________________
void Geo()
{
  TControlBar *pMenu = new TControlBar("vertical","RICH draw");
  pMenu->AddButton("Draw Front", "gMC->Gdraw(\"ALIC\", 0,0,0, 10,10, 0.01,0.01)","Draws ALIC volume in XY view");
  pMenu->AddButton("Draw Side",  "gMC->Gdraw(\"ALIC\",90,180, 0, 10,10, 0.01,0.01)","Draws ALIC volume in YZ view");
  pMenu->AddButton("Draw Top",   "gMC->Gdraw(\"ALIC\",90, 90, 0, 10,10, 0.03,0.03)","Draws ALIC volume in XZ view");
  pMenu->AddButton("ALICE Tree", "((TGeant3*)gMC)->Gdtree(\"ALIC\")","Draws ALICE tree");      
  pMenu->AddButton("RICH Tree",  "((TGeant3*)gMC)->Gdtree(\"RICH\")","Draws RICH tree");      
  pMenu->AddButton("Geo test",  "GeoTest()",   "Draw RICH geo as a macro");
  pMenu->AddButton("Print ref", "PrintGeo()",  "Print RICH chambers default position");
  pMenu->AddButton("Print act", "r->Print();", "Print RICH chambers default position");
  pMenu->Show();  
}//void Draw()

void GeoTest()
{

   TBRIK *pAliceBRIK=new TBRIK("aliceBRIK","ALICE mother volume","void",500,500,500);
   TBRIK *pArmBRIK=new TBRIK("armBRIK","RICH arm1","void",pRICH->GetSizeX(),pRICH->GetSizeY(),pRICH->GetSizeZ());
   
   TNode *pAliceNode=new TNode("aliceNode","Mother volume","aliceBRIK");
   pAliceNode->cd();

// ARM 1 LEFT      
   TRotation rot1;
   TVector3  v1(0,pRICH->GetOffset(),0);
   
         
   rot1.RotateX(pRICH->GetYZAngle()*kDegrad);        v1.RotateX(pRICH->GetYZAngle()*kDegrad);
   rot1.RotateZ(pRICH->GetXYAngle()*kDegrad);        v1.RotateZ(pRICH->GetXYAngle()*kDegrad);
   rot1.RotateZ(pRICH->GetRotAngle()*kDegrad);       v1.RotateZ(pRICH->GetRotAngle()*kDegrad);
         
   TRotMatrix *pArm1RotMatrix=new TRotMatrix("rotArm1","rotArm1",rot1.ThetaX()*kRaddeg, rot1.PhiX()*kRaddeg,
                                                                 rot1.ThetaY()*kRaddeg, rot1.PhiY()*kRaddeg,
                                                                 rot1.ThetaZ()*kRaddeg, rot1.PhiZ()*kRaddeg);

   TNode *pArm1Node=new TNode("arm1Node","Left arm","armBRIK",v1.X(),v1.Y(),v1.Z(),"rotArm1");
   arm1Node->SetLineColor(kRed);
   
// ARM 2 LEFT      
   TRotation rot2;
   TVector3  v2(0,pRICH->GetOffset(),0);
   
         
   rot2.RotateX(pRICH->GetYZAngle()*kDegrad);        v2.RotateX(pRICH->GetYZAngle()*kDegrad);
   rot2.RotateZ(-pRICH->GetXYAngle()*kDegrad);        v2.RotateZ(-pRICH->GetXYAngle()*kDegrad);
   rot2.RotateZ(pRICH->GetRotAngle()*kDegrad);        v2.RotateZ(pRICH->GetRotAngle()*kDegrad);
         
   TRotMatrix *pArm2RotMatrix=new TRotMatrix("rotArm2","rotArm2",rot2.ThetaX()*kRaddeg, rot2.PhiX()*kRaddeg,
                                                                 rot2.ThetaY()*kRaddeg, rot2.PhiY()*kRaddeg,
                                                                 rot2.ThetaZ()*kRaddeg, rot2.PhiZ()*kRaddeg);

   TNode *pArm2Node=new TNode("arm2Node","Left arm","armBRIK",v2.X(),v2.Y(),v2.Z(),"rotArm2");
   arm2Node->SetLineColor(kBlue);
   
   aliceNode->Draw();
}//void GeoTest()
//______________________________________________________________________________
void PrintGeo(Float_t rotDeg=0)
{
  AliRICHParam *p=new AliRICHParam;  
  Double_t r=p->Offset();
  Double_t kP=p->AngleXY();
  Double_t kT=p->AngleYZ();
  Double_t kRot;
  
  if(rotDeg==0)
    kRot=p->AngleRot();
  else
    kRot=rotDeg*deg;
        
  cout<<endl;
  Double_t phi=90*deg+kRot,theta=90*deg-kT;    
  Info("   menu for          0","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
                                 r,      theta*r2d,  phi*r2d,  
                                                               r*sin(theta)*cos(phi),
                                                                       r*sin(theta)*sin(phi),
                                                                               r*cos(theta));
  
  phi=90*deg+kRot+kP,theta=90*deg;
  Info("   menu for          1","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
                                 r,      theta*r2d,  phi*r2d,  
                                                               r*sin(theta)*cos(phi),
                                                                       r*sin(theta)*sin(phi),
                                                                               r*cos(theta));
  
  phi=90*deg+kRot,theta=90*deg;
  Info("   menu for          2","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
                                 r,      theta*r2d,  phi*r2d,  
                                                               r*sin(theta)*cos(phi),
                                                                       r*sin(theta)*sin(phi),
                                                                               r*cos(theta));

  
  phi=90*deg+kRot-kP,theta=90*deg;
  Info("   menu for          3","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
                                 r,      theta*r2d,  phi*r2d,  
                                                               r*sin(theta)*cos(phi),
                                                                       r*sin(theta)*sin(phi),
                                                                               r*cos(theta));


  phi=90*deg+kRot+kP,theta=90*deg+kT;
  Info("   menu for          4","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
                                 r,      theta*r2d,  phi*r2d,  
                                                               r*sin(theta)*cos(phi),
                                                                       r*sin(theta)*sin(phi),
                                                                               r*cos(theta));
  
  
  phi=90*deg+kRot,theta=90*deg+kT;
  Info("   menu for          5","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
                                 r,      theta*r2d,  phi*r2d,  
                                                                r*sin(theta)*cos(phi),
                                                                       r*sin(theta)*sin(phi),
                                                                               r*cos(theta));
  
  phi=90*deg+kRot-kP,theta=90*deg+kT;
  Info("   menu for          6","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
                                 r,      theta*r2d,  phi*r2d,  
                                                               r*sin(theta)*cos(phi),
                                                                       r*sin(theta)*sin(phi),
                                                                               r*cos(theta));

  delete p;
}//void PrintGeo()
Commit	Line	Data
d48cca74	1	void Show()
	2	{
	3	//How to get number of events:
853634d3	4	Info("RICH/menu.C::Show","3 ways to get number of events 1=%i 2=%f 3=%f",
d48cca74	5	gAlice->GetEventsPerRun(),
853634d3	6	al->TreeE()->GetEntries(),
d48cca74	7	gAlice->TreeE()->GetEntries());
853634d3	8	rl->LoadHits();
	9	for(Int_t iEventN=0;iEventN<gAlice->GetEventsPerRun();iEventN++){// loop on events
	10	Int_t iNparticles=gAlice->GetEvent(iEventN);
	11	Int_t iNtracks=rl->TreeH()->GetEntries();
	12	Info("RICH/menu.C::Show","Event %i contains %i particles in total while %i are primary",
	13	iEventN, iNparticles, iNtracks);
	14	for(Int_t iTrackN=0;iTrackN<iNtracks;iTrackN++){ // loop on tracks
	15	rl->TreeH()->GetEntry(iTrackN);
	16	al->Stack()->Particle(iTrackN)->Print();
	17	Info("RICH/menu.C::Show","track %i has %i hits",iTrackN,r->Hits()->GetEntries());
	18	// for(AliRICHhit pRichHit=(AliRICHhit)gpRich->FirstHit(-1);pRichHit;pRichHit=(AliRICHhit*)gpRich->NextHit()){
	19	// pRichHit->Dump();
	20	// TVector3 armV3=pRICH->ToArm(mrsV3);
	21	// cout<<"After\n";armV3.Dump();
	22	// }//loop on hits of given track
	23	}// loop on tracks
	24	}//loop on events
d48cca74	25	}//void Show()
	26
	27	void menu(Int_t iNevents=5)// How many events to generate.
	28	{
	29	Info("RICH/menu.C","%i event(s) are requested",iNevents);
	30
	31	TString runString="gAlice->Run("; runString=runString+iNevents+");";
	32
	33	TControlBar *pMenu = new TControlBar("vertical","RICH main");
	34
853634d3	35	pMenu->AddButton("Debug ON-OFF", "Debug();", "Switch debug on-off");
d48cca74	36	if(GetAlice()){//it's from file, reconstruct
d48cca74	37	pMenu->AddButton("Show","Show()","Shows the structure of events in files");
853634d3	38	pMenu->AddButton("Hits2SDigits","r->Hits2SDigits()","Perform first phase converstion");
d48cca74	39	pMenu->AddButton("RingViewer","RingViewer()","Show rings with reconstructed info");
d48cca74	40	}else{//it's aliroot, simulate
d48cca74	41	pMenu->AddButton("Run", runString.Data(), "Process!");
d48cca74	42	}
d48cca74	43	pMenu->AddButton("Geo", "Geo()", "Geomentry submenu");
	44	pMenu->AddButton("Browser", "new TBrowser;", "Start ROOT TBrowser");
	45	pMenu->AddButton("Quit", ".q", "Close session");
	46	pMenu->Show();
853634d3	47	a=gAlice;//for manual manipulation convinience
d48cca74	48	}//void menu(Int_t iNevents)
d48cca74	49
853634d3	50
853634d3	51	void Debug()
d48cca74	52	{
853634d3	53	if(gAlice->GetDebug()){
	54	Info("RICH/menu.C::Debug","OFF");
	55	gAlice->SetDebug(0);
	56	if(r)r->SetDebug(0);
	57	if(t)t->SetDebug(0);
	58	AliLoader::SetDebug(0);
	59	}else{
	60	Info("RICH/menu.C::Debug","ON");
	61	gAlice->SetDebug(1);
	62	if(r)r->SetDebug(1);
	63	if(t)t->SetDebug(1);
	64	AliLoader::SetDebug(1);
	65	}
	66	}//void Debug()
d48cca74	67
d48cca74	68	AliRun *a;
853634d3	69	AliRunLoader *al;
	70	AliLoader rl,tl,*il;
	71
	72	AliRICH *r;
	73	AliTPC *t;
	74	AliITS *i;
	75	AliPHOS *p;
	76
d48cca74	77	Bool_t GetAlice()
	78	{
	79	if(gAlice){//it's aliroot
853634d3	80	Info("RICH/menu.C::GetAlice","gAlice!=NULL, IT'S ALIROOT, EXECUTE Init()");
	81	gAlice->Init("~/my/AliceConfig.C");
	82	r=(AliRICH*)gAlice->GetDetector("RICH");
	83	t=(AliTPC*)gAlice->GetDetector("TPC");
	84	i=(AliITS*)gAlice->GetDetector("ITS");
	85	p=(AliPHOS*)gAlice->GetDetector("PHOS");
d48cca74	86	return kFALSE;
d48cca74	87	}else{//it's root with ALICE libs loaded
853634d3	88	Info("RICH/menu.C::GetAlice","gAlice=0, GETTING IT FROM SIMULATED FILE.");
	89
	90	if(!(al=AliRunLoader::Open("galice.root","AlicE","update"))) Fatal("RICH/menu.C::GetAlice","Can't get AliRunLoader");
	91	al->LoadgAlice();
d48cca74	92	if(!gAlice) Fatal("RICH/menu.C::GetAlice","No gAlice in file");
853634d3	93	// a=al->GetAliRun();
	94	al->LoadHeader();//loads events tree
	95	al->LoadKinematics();//loads the primaries info
	96	//RICH
	97	if(!(r=(AliRICH*)gAlice->GetDetector("RICH"))) Warning("RICH/menu.C::GetAlice","No RICH in file");
	98	if(!(rl=al->GetLoader("RICHLoader"))) Warning("RICH/menu.C::GetAlice","No RICH loader in file");
	99	//TPC
	100	if(!(t=(AliTPC*)gAlice->GetDetector("TPC"))) Warning("RICH/menu.C::GetAlice","No TPC in file");
	101	if(!(tl=al->GetLoader("TPCLoader"))) Warning("RICH/menu.C::GetAlice","No TPC loader in file");
	102
	103	Info("RICH/menu.C::GetAlice","Run contains %i event(s)",gAlice->GetEventsPerRun());
d48cca74	104	return kTRUE;
	105	}
	106	}//void GetAlice()
	107
853634d3	108	void RingViewer()
	109	{
	110	gStyle->SetPalette(1);
	111	TCanvas *view=new TCanvas("Display","ALICE RICH Display",0,0,1200,750);
	112
	113	TH2F *pH2=new TH2F("pH2F","RICH DISPLAY",160,0,160,144,0,144);
	114	pH2->SetStats(0);
	115	pH2->SetMaximum(100);
d48cca74	116
853634d3	117	Int_t Nevents = gAlice->GetEventsPerRun();
853634d3	118	}
d48cca74	119
853634d3	120	//______________________________________________________________________________
d48cca74	121	void Geo()
	122	{
	123	TControlBar *pMenu = new TControlBar("vertical","RICH draw");
	124	pMenu->AddButton("Draw Front", "gMC->Gdraw(\"ALIC\", 0,0,0, 10,10, 0.01,0.01)","Draws ALIC volume in XY view");
	125	pMenu->AddButton("Draw Side", "gMC->Gdraw(\"ALIC\",90,180, 0, 10,10, 0.01,0.01)","Draws ALIC volume in YZ view");
	126	pMenu->AddButton("Draw Top", "gMC->Gdraw(\"ALIC\",90, 90, 0, 10,10, 0.03,0.03)","Draws ALIC volume in XZ view");
	127	pMenu->AddButton("ALICE Tree", "((TGeant3*)gMC)->Gdtree(\"ALIC\")","Draws ALICE tree");
	128	pMenu->AddButton("RICH Tree", "((TGeant3*)gMC)->Gdtree(\"RICH\")","Draws RICH tree");
	129	pMenu->AddButton("Geo test", "GeoTest()", "Draw RICH geo as a macro");
	130	pMenu->AddButton("Print ref", "PrintGeo()", "Print RICH chambers default position");
	131	pMenu->AddButton("Print act", "r->Print();", "Print RICH chambers default position");
	132	pMenu->Show();
	133	}//void Draw()
	134
	135	void GeoTest()
	136	{
	137
	138	TBRIK *pAliceBRIK=new TBRIK("aliceBRIK","ALICE mother volume","void",500,500,500);
d48cca74	139	TBRIK *pArmBRIK=new TBRIK("armBRIK","RICH arm1","void",pRICH->GetSizeX(),pRICH->GetSizeY(),pRICH->GetSizeZ());
	140
	141	TNode *pAliceNode=new TNode("aliceNode","Mother volume","aliceBRIK");
	142	pAliceNode->cd();
	143
	144	// ARM 1 LEFT
	145	TRotation rot1;
	146	TVector3 v1(0,pRICH->GetOffset(),0);
	147
	148
	149	rot1.RotateX(pRICH->GetYZAngle()kDegrad); v1.RotateX(pRICH->GetYZAngle()kDegrad);
	150	rot1.RotateZ(pRICH->GetXYAngle()kDegrad); v1.RotateZ(pRICH->GetXYAngle()kDegrad);
	151	rot1.RotateZ(pRICH->GetRotAngle()kDegrad); v1.RotateZ(pRICH->GetRotAngle()kDegrad);
	152
	153	TRotMatrix pArm1RotMatrix=new TRotMatrix("rotArm1","rotArm1",rot1.ThetaX()kRaddeg, rot1.PhiX()*kRaddeg,
	154	rot1.ThetaY()kRaddeg, rot1.PhiY()kRaddeg,
	155	rot1.ThetaZ()kRaddeg, rot1.PhiZ()kRaddeg);
	156
	157	TNode *pArm1Node=new TNode("arm1Node","Left arm","armBRIK",v1.X(),v1.Y(),v1.Z(),"rotArm1");
	158	arm1Node->SetLineColor(kRed);
	159
	160	// ARM 2 LEFT
	161	TRotation rot2;
	162	TVector3 v2(0,pRICH->GetOffset(),0);
	163
	164
	165	rot2.RotateX(pRICH->GetYZAngle()kDegrad); v2.RotateX(pRICH->GetYZAngle()kDegrad);
	166	rot2.RotateZ(-pRICH->GetXYAngle()kDegrad); v2.RotateZ(-pRICH->GetXYAngle()kDegrad);
	167	rot2.RotateZ(pRICH->GetRotAngle()kDegrad); v2.RotateZ(pRICH->GetRotAngle()kDegrad);
	168
	169	TRotMatrix pArm2RotMatrix=new TRotMatrix("rotArm2","rotArm2",rot2.ThetaX()kRaddeg, rot2.PhiX()*kRaddeg,
	170	rot2.ThetaY()kRaddeg, rot2.PhiY()kRaddeg,
	171	rot2.ThetaZ()kRaddeg, rot2.PhiZ()kRaddeg);
	172
	173	TNode *pArm2Node=new TNode("arm2Node","Left arm","armBRIK",v2.X(),v2.Y(),v2.Z(),"rotArm2");
	174	arm2Node->SetLineColor(kBlue);
	175
	176	aliceNode->Draw();
	177	}//void GeoTest()
853634d3	178	//______________________________________________________________________________
d48cca74	179	void PrintGeo(Float_t rotDeg=0)
	180	{
	181	AliRICHParam *p=new AliRICHParam;
	182	Double_t r=p->Offset();
	183	Double_t kP=p->AngleXY();
	184	Double_t kT=p->AngleYZ();
	185	Double_t kRot;
	186
	187	if(rotDeg==0)
	188	kRot=p->AngleRot();
	189	else
	190	kRot=rotDeg*deg;
	191
	192	cout<<endl;
	193	Double_t phi=90deg+kRot,theta=90deg-kT;
	194	Info(" menu for 0","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
853634d3	195	r, thetar2d, phir2d,
d48cca74	196	rsin(theta)cos(phi),
	197	rsin(theta)sin(phi),
	198	r*cos(theta));
	199
	200	phi=90deg+kRot+kP,theta=90deg;
	201	Info(" menu for 1","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
853634d3	202	r, thetar2d, phir2d,
d48cca74	203	rsin(theta)cos(phi),
	204	rsin(theta)sin(phi),
	205	r*cos(theta));
	206
	207	phi=90deg+kRot,theta=90deg;
	208	Info(" menu for 2","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
853634d3	209	r, thetar2d, phir2d,
d48cca74	210	rsin(theta)cos(phi),
	211	rsin(theta)sin(phi),
	212	r*cos(theta));
	213
	214
	215	phi=90deg+kRot-kP,theta=90deg;
	216	Info(" menu for 3","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
853634d3	217	r, thetar2d, phir2d,
d48cca74	218	rsin(theta)cos(phi),
	219	rsin(theta)sin(phi),
	220	r*cos(theta));
	221
	222
	223	phi=90deg+kRot+kP,theta=90deg+kT;
	224	Info(" menu for 4","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
853634d3	225	r, thetar2d, phir2d,
d48cca74	226	rsin(theta)cos(phi),
	227	rsin(theta)sin(phi),
	228	r*cos(theta));
	229
	230
	231	phi=90deg+kRot,theta=90deg+kT;
	232	Info(" menu for 5","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
853634d3	233	r, thetar2d, phir2d,
d48cca74	234	rsin(theta)cos(phi),
	235	rsin(theta)sin(phi),
	236	r*cos(theta));
	237
	238	phi=90deg+kRot-kP,theta=90deg+kT;
	239	Info(" menu for 6","r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f",
853634d3	240	r, thetar2d, phir2d,
d48cca74	241	rsin(theta)cos(phi),
	242	rsin(theta)sin(phi),
	243	r*cos(theta));
	244
	245	delete p;
	246	}//void PrintGeo()