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

Timestamp:

Apr 20, 2005, 7:03:24 PM (20 years ago)

Author:

cmv

Message:

1-/ pour n/proj pas de display si histo existe ET pas d'option graphiques donnees
2-/ gestion correcte de la superposition d'histo 1D dans n/plot nt.x

L'histo pour le 2ieme display est dimensionne comme le premier.

(ajout de l'option graphique keepbin)

cmv 20/04/2005

File:

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

trunk/SophyaPI/PIext/pawexecut.cc

@@ -19 +19 @@
 #include "nomgadapter.h"
 #include "pistdimgapp.h"
+#include "pihisto.h"
 
 #ifdef SANS_EVOLPLANCK
@@ -45 +46 @@
 
 kw = "n/plot";
-usage = "Plot NTuple variables a la paw (alias n/pl)";
+usage = "Plot NTuple variables \"a la paw\" (alias n/pl)";
 usage += "\n n/plot nameobj.x_exp [cut] [w_exp] [loop] [gratt]";
 usage += "\n n/plot nameobj.y_exp%x_exp [cut] [loop] [gratt]";
 usage += "\n n/plot nameobj.z_exp%y_exp%x_exp [cut] [loop] [gratt]";
 usage += "\n        for default use ! , loop=i1[:i2[:di]]";
+usage += "\n for 1 dimensional (1D) projection:";
+usage += "\n   use graphic option \"keepbin\" to superimpose";
+usage += "\n   1D distribution with the same binning";
 usage += "\n  Related commands: plot2dw plot3d";  
 piac->RegisterCommand(kw,usage,this,hgrp);
@@ -58 +62 @@
 usage += "\n n/proj nameproj nameobj.y_exp%x_exp [cut] [w_exp] [loop] [gratt]";
 usage += "\n        for default use ! , loop=i1[:i2[:di]]";
+usage += "\n for 1 dimensional (1D) projection:";
+usage += "\n   no display is performed if \"nameproj\" is an existing histogram";
+usage += "\n   unless a graphic option \"gratt\" is given";
 usage += "\n  Related commands: projh1d projh2d projprof exptovec";  
 piac->RegisterCommand(kw,usage,this,hgrp);
@@ -375 +382 @@
       nameproj += buff;
     }
+    // --- Si option "keepbin" on projete avec le meme binning que l'histo 1D precedent
+    if(dopt.find("keepbin")<dopt.size()) {
+      PIStdImgApp* piimapp = omg.GetImgApp();
+      if(piimapp) {
+        PIBaseWdg* pibwdg = piimapp->CurrentBaseWdg();
+        if(pibwdg) {
+          int nbdrw = pibwdg->NbDrawers();
+          if(nbdrw>0) {
+            for(int i=nbdrw-1;i>=0;i--) {
+              PIDrawer* pidwr = pibwdg->GetDrawer(i);
+              PIHisto* pih = NULL;
+              if( (pih = dynamic_cast<PIHisto *>(pidwr))==NULL ) continue;
+              Histo* h = pih->Histogram();
+              if(h==NULL) continue;
+              if(h->NBins()<1) continue;
+              Histo* hsame = new Histo(h->XMin(),h->XMax(),h->NBins());
+              omg.AddObj(hsame,nameproj);
+              // on force donc le display sur le meme plot
+              dopt += " same";
+              break;
+            }
+          }
+        }
+      }
+    }
+    // --- 
   }
   srvo->ProjectH1(nameobj,expx,expwt,expcut,nameproj,dopt,loop);