Changeset 3058 in Sophya for trunk/SophyaPI/PIext/pawexecut.cc

Timestamp:

Aug 13, 2006, 3:15:57 PM (19 years ago)

Author:

cmv

Message:

pawexecutor: copie dans Histo des projections assiciees a un Histo2D cmv 13/8/2006

File:

• trunk/SophyaPI/PIext/pawexecut.cc (modified) (4 diffs)

trunk/SophyaPI/PIext/pawexecut.cc

@@ -218 +218 @@
 kw = "h/plot/2d";
 usage = "Specific plot for 2D histogrammes";
-usage += "\n h/plot/2d nameh2d show        : infos on 2D histogramme";
+usage += "\n h/plot/2d nameh2d show        : infos on 2D histogramme projections";
 usage += "\n h/plot/2d nameh2d h    [dopt] : plot 2D histogramme";
 usage += "\n h/plot/2d nameh2d px   [dopt] : plot X projection";
@@ -227 +227 @@
 usage += "\n h/plot/2d nameh2d sy n [dopt] : plot Y slice number n";
 usage += "\n                      n < 0 means Show Info";
+piac->RegisterCommand(kw,usage,this,hgrp);
+
+kw = "h/2d/geth";
+usage = "Get histo projections, bandes and slices";
+usage += "\n h/2d/geth nameh2d         : infos on 2D histogramme projections";
+usage += "\n h/2d/geth h nameh2d px    : copy X projection into histo h";
+usage += "\n h/2d/geth h nameh2d py    : copy Y projection into histo h";
+usage += "\n h/2d/geth h nameh2d bx n  : copy X band number n into histo h";
+usage += "\n h/2d/geth h nameh2d by n  : copy Y band number n into histo h";
+usage += "\n h/2d/geth h nameh2d sx n  : copy X slice number n into histo h";
+usage += "\n h/2d/geth h nameh2d sy n  : copy Y slice number n into histo h";
+usage += "\n  Related commands: h/plot/2d";  
 piac->RegisterCommand(kw,usage,this,hgrp);
 
@@ -346 +358 @@
 } else if(kw == "h/plot/2d") {
   h_plot_2d(tokens); return(0);
+} else if(kw == "h/2d/geth") {
+  h_2d_geth(tokens); return(0);
 } else if(kw == "h/put_vec") {
   h_put_vec(tokens); return(0);
@@ -1837 +1851 @@
 
 /* methode */
+void PAWExecutor::h_2d_geth(vector<string>& tokens)
+// copy bandx/y, slicex/y or projx/y from an 2D histo into a 1D histo
+{
+ NamedObjMgr omg;
+
+ if(tokens.size()==1) {
+   AnyDataObj* mobj = omg.GetObj(tokens[0]);
+   Histo2D* h2 = dynamic_cast<Histo2D*>(mobj);
+   if(!h2) return;
+   h2->ShowProj();
+   h2->ShowBand(2);
+   h2->ShowSli(2);
+   return;
+ }
+
+ if(tokens.size()<3)
+   {cout<<"Usage: h/2d/geth h nameh2d [px,py,bx n,by n,sx n,sy n]"<<endl; return;}
+
+ string nameh1 = tokens[0];
+ string nameh2 = tokens[1];
+ string proj   = tokens[2];
+ int_4 nump = (tokens.size()>3)? atoi(tokens[3].c_str()): -1;
+
+ // Decodage Histo2D
+ AnyDataObj* mobj = omg.GetObj(tokens[1]);
+ if(mobj==NULL)
+   {cout<<"PAWExecutor::h_2d_geth Error: unknow object"<<tokens[1]<<endl;
+   return;}
+ Histo2D* h2 = dynamic_cast<Histo2D*>(mobj);
+ if(!h2)
+   {cout<<"PAWExecutor::h_plot_2d Error: "<<tokens[1]<<" not an Histo2D"<<endl;
+   return;}
+
+ // Remplissage histo 1D avec la projection demandee
+ Histo *h1p;
+ if(proj == "px") {  
+   if((h1p=h2->HProjX())) {Histo* h=new Histo(*h1p); omg.AddObj(h,nameh1);}
+     else h2->ShowProj();
+ } else if(proj == "py") {
+   if((h1p=h2->HProjY())) {Histo* h=new Histo(*h1p); omg.AddObj(h,nameh1);}
+     else h2->ShowProj();
+ } else if(proj == "bx"  && nump>=0) {
+   if((h1p=h2->HBandX(nump))) {Histo* h=new Histo(*h1p); omg.AddObj(h,nameh1);}
+     else h2->ShowBand();
+ } else if(proj == "by"  && nump>=0) {
+   if((h1p=h2->HBandY(nump))) {Histo* h=new Histo(*h1p); omg.AddObj(h,nameh1);}
+     else h2->ShowBand();
+ } else if(proj == "sx"  && nump>=0) {
+   if((h1p=h2->HSliX(nump))) {Histo* h=new Histo(*h1p); omg.AddObj(h,nameh1);}
+     else h2->ShowSli();
+ } else if(proj == "sy"  && nump>=0) {
+   if((h1p=h2->HSliY(nump))) {Histo* h=new Histo(*h1p); omg.AddObj(h,nameh1);}
+     else h2->ShowSli();
+ } else {
+   cout<<"PAWExecutor::h_2d_geth Error: Unknown proj name "<<proj
+       <<" or bad projection number "<<nump<<endl;
+ }
+
+ return;
+}
+
+/* methode */
 int_4 PAWExecutor::decodepawstring(string tokens,string& nameobj
                                 ,string& xexp,string& yexp,string& zexp)