source: Sophya/trunk/SophyaPI/PIext/pawexecut.cc@ 1059

Last change on this file since 1059 was 1059, checked in by ercodmgr, 25 years ago

re-activation methode HRebin pour HProf apres modifs HiStats cmv 7/7/2000
File size: 20.4 KB
Line 
1#include <stdio.h>
2#include <stdlib.h>
3#include <ctype.h>
4#include <iostream.h>
5#include <typeinfo>
6
7#include "strutil.h"
8#include "histos.h"
9#include "histos2.h"
10#include "hisprof.h"
11#include "ntuple.h"
12
13#include "pawexecut.h"
14#include "nobjmgr.h"
15#include "servnobjm.h"
16#include "pistdimgapp.h"
17
18#ifdef SANS_EVOLPLANCK
19#include "cvector.h"
20#include "matrix.h"
21#else
22#include "tmatrix.h"
23#include "tvector.h"
24#endif
25
26/* Reza + cmv 13/10/99 */
27
28uint_4 PAWExecutor::autoc_counter_ = 0;
29
30/* methode */
31PAWExecutor::PAWExecutor(PIACmd *piac, PIStdImgApp* app)
32: mApp(app)
33{
34string kw, usage;
35string hgrp = "pawCmd";
36
37kw = "reset";
38usage = "Reset histograms vectors or matrix";
39usage += "\n reset nameobj";
40piac->RegisterCommand(kw,usage,this,hgrp);
41
42kw = "n/plot";
43usage = "Plot NTuple variables a la paw";
44usage += "\n n/plot nameobj.x_exp [cut] [w_exp] [loop] [gratt]";
45usage += "\n n/plot nameobj.y_exp%x_exp [cut] [loop] [gratt]";
46usage += "\n n/plot nameobj.z_exp%y_exp%x_exp [cut] [loop] [gratt]";
47usage += "\n for default use ! , loop=i1[:i2[:di]]";
48usage += "\n Related commands: plot2dw plot3d";
49piac->RegisterCommand(kw,usage,this,hgrp);
50
51kw = "n/proj";
52usage = "Project NTuple in histogram (1D or 2D) a la paw";
53usage += "\n n/proj nameproj nameobj.x_exp [cut] [w_exp] [loop] [gratt]";
54usage += "\n n/proj nameproj nameobj.y_exp%x_exp [cut] [w_exp] [loop] [gratt]";
55usage += "\n for default use ! , loop=i1[:i2[:di]]";
56usage += "\n Related commands: projh1d projh2d projprof";
57piac->RegisterCommand(kw,usage,this,hgrp);
58
59kw = "h/integ";
60usage = "Integrate a 1D histogram";
61usage += "\n h/integ nameh1d [norm]";
62usage += "\n Related commands: h/deriv";
63piac->RegisterCommand(kw,usage,this,hgrp);
64
65kw = "h/deriv";
66usage = "Derivate a 1D histogram";
67usage += "\n h/deriv nameh1d";
68usage += "\n Related commands: h/integ";
69piac->RegisterCommand(kw,usage,this,hgrp);
70
71kw = "h/rebin";
72usage = "Rebin a 1D histogram or profile";
73usage += "\n h/rebin nameh1d nbin";
74piac->RegisterCommand(kw,usage,this,hgrp);
75
76kw = "h/cadd";
77usage = "Add a constant to an histogram";
78usage += "\n h/cadd namehisto val";
79usage += "\n Related commands: h/cmult h/oper";
80piac->RegisterCommand(kw,usage,this,hgrp);
81
82kw = "h/cmult";
83usage = "Multiply an histogram by a constant";
84usage += "\n h/cmult namehisto val";
85usage += "\n Related commands: h/cadd h/oper";
86piac->RegisterCommand(kw,usage,this,hgrp);
87
88kw = "h/oper";
89usage = "Operation on histograms";
90usage += "\n h/oper @ h1 h2 hres";
91usage += "\n hres = h1 @ h2 with @ = (+,-,*,/)";
92usage += "\n Related commands: h/cadd h/cmult";
93piac->RegisterCommand(kw,usage,this,hgrp);
94
95kw = "h/plot/2d";
96usage = "Specific plot for 2D histogrammes";
97usage += "\n h/plot/2d nameh2d show : infos on 2D histogramme";
98usage += "\n h/plot/2d nameh2d h [dopt] : plot 2D histogramme";
99usage += "\n h/plot/2d nameh2d px [dopt] : plot X projection";
100usage += "\n h/plot/2d nameh2d py [dopt] : plot Y projection";
101usage += "\n h/plot/2d nameh2d bx n [dopt] : plot X band number n";
102usage += "\n h/plot/2d nameh2d by n [dopt] : plot Y band number n";
103usage += "\n h/plot/2d nameh2d sx n [dopt] : plot X slice number n";
104usage += "\n h/plot/2d nameh2d sy n [dopt] : plot Y slice number n";
105usage += "\n n < 0 means Show Info";
106piac->RegisterCommand(kw,usage,this,hgrp);
107}
108
109/* methode */
110PAWExecutor::~PAWExecutor()
111{
112}
113
114/* methode */
115int PAWExecutor::Execute(string& kw, vector<string>& tokens)
116{
117if(kw == "reset") {
118 reset(tokens); return(0);
119} else if(kw == "n/plot") {
120 n_plot(tokens); return(0);
121} else if(kw == "n/proj") {
122 n_proj(tokens); return(0);
123} else if(kw == "h/integ") {
124 h_integ(tokens); return(0);
125} else if(kw == "h/deriv") {
126 h_deriv(tokens); return(0);
127} else if(kw == "h/rebin") {
128 h_rebin(tokens); return(0);
129} else if(kw == "h/cadd") {
130 h_cadd(tokens); return(0);
131} else if(kw == "h/cmult") {
132 h_cmult(tokens); return(0);
133} else if(kw == "h/oper") {
134 h_oper(tokens); return(0);
135} else if(kw == "h/plot/2d") {
136 h_plot_2d(tokens); return(0);
137} else return(1);
138}
139
140/* methode */
141void PAWExecutor::reset(vector<string>& tokens)
142// Reset d'histogrammes, vecteurs et matrices
143{
144if(tokens.size() < 1)
145 {cout<<"Usage: reset nameobj"<<endl; return;}
146NamedObjMgr omg;
147AnyDataObj* mobj = omg.GetObj(tokens[0]);
148if(mobj == NULL)
149 {cout<<"PAWExecutor::reset Error , Pas d'objet de nom "<<tokens[0]<<endl;
150 return;}
151string ctyp = typeid(*mobj).name();
152
153#ifdef SANS_EVOLPLANCK
154if(typeid(*mobj)==typeid(Vector)) {Vector* ob=(Vector*) mobj; ob->Zero();}
155else if(typeid(*mobj)==typeid(Matrix)) {Matrix* ob=(Matrix*) mobj; ob->Zero();}
156#else
157 if(typeid(*mobj)==typeid(Vector)) {Vector* ob=(Vector*) mobj; (*ob) = 0.; }
158// ob->DataBlock().Reset(0.);}
159 else if(typeid(*mobj)==typeid(Matrix)) {Matrix* ob=(Matrix*) mobj; (*ob) = 0.; }
160//ob->DataBlock().Reset(0.);}
161#endif
162else if(typeid(*mobj)==typeid(Histo)) {Histo* ob=(Histo*) mobj; ob->Zero();}
163else if(typeid(*mobj)==typeid(HProf)) {HProf* ob=(HProf*) mobj; ob->Zero();}
164else if(typeid(*mobj)==typeid(Histo2D)) {Histo2D* ob=(Histo2D*)mobj; ob->Zero();}
165else {
166 cout<<"PAWExecutor::reset Error , No reset possible on "<<ctyp<<endl;
167 return;
168}
169
170return;
171}
172
173/* methode */
174void PAWExecutor::n_plot(vector<string>& tokens)
175// Equivalent n/plot de paw
176// Plot 1D
177// n/plot nameobj.x_exp [cut] [w_exp] [gratt]
178// Plot 2D (plot2dw)
179// n/plot nameobj.y_exp%x_exp [cut] [w_exp] [gratt]
180// Plot 3D (plot3d)
181// n/plot nameobj.z_exp%y_exp%x_exp [cut] [gratt]
182{
183if(tokens.size() < 1) {
184 cout
185 <<"Usage: n/plot nameobj.x_exp [cut] [w_exp] [loop] [gratt] [nomh1]"<<endl
186 <<" n/plot nameobj.y_exp%x_exp [cut] [loop] [gratt]"<<endl
187 <<" n/plot nameobj.z_exp%y_exp%x_exp [cut] [loop] [gratt]"<<endl
188 <<" for default use ! , loop=i1[:i2[:di]]"<<endl;
189 return;
190}
191string nameobj,expx,expy,expz;
192int nvar = decodepawstring(tokens[0],nameobj,expx,expy,expz);
193string expcut = "1"; string expwt = "1."; string loop = "";
194string dopt = ""; string nameproj="";
195if(tokens.size()>=2) expcut = tokens[1]; if(expcut=="!") expcut="1";
196
197NamedObjMgr omg;
198Services2NObjMgr* srvo = omg.GetServiceObj();
199
200if(nvar<=0) {
201 cout<<"PAWExecutor::n_plot Error: bad coding "<<tokens[0]<<endl;
202} else if(nvar==1) { // c'est un plot 1D
203 if(tokens.size()>=3) expwt = tokens[2]; if(expwt=="!") expwt="1.";
204 if(tokens.size()>=4) loop = tokens[3]; if(loop=="!") loop="";
205 if(tokens.size()>=5) dopt = tokens[4]; if(dopt=="!") dopt="";
206 if(tokens.size()>=6) nameproj = tokens[5];
207 if(nameproj.length()<=0 || nameproj=="!") {
208 nameproj = "/autoc/paw_n_plot1D_";
209 AnyDataObj* mobj = omg.GetObj(nameproj);
210 if(mobj!=NULL) {
211 char buff[16]; autoc_counter_++; sprintf(buff,"%d",autoc_counter_);
212 nameproj += buff;
213 }
214 }
215 srvo->ProjectH1(nameobj,expx,expwt,expcut,nameproj,dopt,loop);
216} else if(nvar==2) { // c'est un plot 2D
217 if(tokens.size()>=3) loop = tokens[2]; if(loop=="!") loop="";
218 if(tokens.size()>=4) dopt = tokens[3];
219 string err = "";
220 srvo->DisplayPoints2D(nameobj,expx,expy,err,err,expcut,dopt,loop);
221} else { // c'est un plot 3D
222 if(tokens.size()>=3) loop = tokens[2]; if(loop=="!") loop="";
223 if(tokens.size()>=4) dopt = tokens[3];
224 srvo->DisplayPoints3D(nameobj,expx,expy,expz,expcut,dopt,loop);
225}
226
227return;
228}
229
230/* methode */
231void PAWExecutor::n_proj(vector<string>& tokens)
232// Equivalent n/proj de paw
233// Project NTuple in histogram a la paw
234// Dans un Histo 1D
235// n/proj nameproj nameobj.x_exp [cut] [w_exp] [gratt]
236// Dans un Histo 2D ou un HProf (dans ce cas nameproj doit etre cree).
237// n/proj nameproj nameobj.y_exp%x_exp [cut] [w_exp] [gratt]
238{
239if(tokens.size()<2)
240 {cout<<"Usage: n/proj nameproj nameobj.[y_exp%]x_exp [cut] [w_exp] [loop] [gratt]"<<endl
241 <<" for default use ! , loop=i1[:i2[:di]]"<<endl; return;}
242string nameproj = tokens[0];
243string nameobj,expx,expy,expz;
244int nvar = decodepawstring(tokens[1],nameobj,expx,expy,expz);
245string expcut = "1"; string expwt = "1."; string loop = ""; string dopt = "";
246if(tokens.size()>=3) expcut = tokens[2]; if(expcut=="!") expcut="1";
247if(tokens.size()>=4) expwt = tokens[3]; if(expwt=="!") expwt="1.";
248if(tokens.size()>=5) loop = tokens[4]; if(loop=="!") loop="";
249if(tokens.size()>=6) dopt = tokens[5];
250
251NamedObjMgr omg;
252Services2NObjMgr* srvo = omg.GetServiceObj();
253
254if(nvar<=0) {
255 cout<<"PAWExecutor::n_proj Error: bad coding "<<tokens[1]<<endl;
256} else if(nvar==1) {
257 // c'est une projection dans un histo 1D
258 srvo->ProjectH1(nameobj,expx,expwt,expcut,nameproj,dopt,loop);
259} else {
260 // c'est une projection dans un histo2D
261 // OU un HProf si nameproj est HProf un deja defini
262 AnyDataObj* mobj = omg.GetObj(nameproj);
263 if(mobj==NULL)
264 srvo->ProjectH2(nameobj,expx,expy,expwt,expcut,nameproj,dopt,loop);
265 else if(dynamic_cast<HProf*>(mobj))
266 srvo->ProjectHProf(nameobj,expx,expy,expwt,expcut,nameproj,dopt,loop);
267 else
268 srvo->ProjectH2(nameobj,expx,expy,expwt,expcut,nameproj,dopt,loop);
269}
270
271return;
272}
273
274/* methode */
275void PAWExecutor::h_integ(vector<string>& tokens)
276// Pour remplacer le contenu d'un histo 1D par son integrale
277{
278if(tokens.size()<1)
279 {cout<<"Usage: h/integ nameh1d [norm]"<<endl; return;}
280NamedObjMgr omg;
281AnyDataObj* mobj = omg.GetObj(tokens[0]);
282if(mobj==NULL)
283 {cout<<"PAWExecutor::h_integ Error: unknow object"<<tokens[0]<<endl;
284 return;}
285double norm = 1.;
286if(tokens.size()>=2) norm = atof(tokens[1].c_str());
287// attention: dynamic_cast<Histo*>(HProf)=Vrai!
288Histo* h1 = dynamic_cast<Histo*>(mobj);
289HProf* hp = dynamic_cast<HProf*>(mobj);
290if(hp || !h1)
291 {cout<<"PAWExecutor::h_integ Error: "<<tokens[0]<<" not an Histo"<<endl;
292 return;}
293h1->HInteg(norm);
294}
295
296/* methode */
297void PAWExecutor::h_deriv(vector<string>& tokens)
298// Pour remplacer le contenu d'un histo 1D par sa derivee
299{
300if(tokens.size()<1)
301 {cout<<"Usage: h/deriv nameh1d"<<endl; return;}
302NamedObjMgr omg;
303AnyDataObj* mobj = omg.GetObj(tokens[0]);
304if(mobj==NULL)
305 {cout<<"PAWExecutor::h_deriv Error: unknow object"<<tokens[0]<<endl;
306 return;}
307// attention: dynamic_cast<Histo*>(HProf)=Vrai!
308Histo* h1 = dynamic_cast<Histo*>(mobj);
309HProf* hp = dynamic_cast<HProf*>(mobj);
310if(hp || !h1)
311 {cout<<"PAWExecutor::h_deriv Error: "<<tokens[0]<<" not an Histo"<<endl;
312 return;}
313h1->HDeriv();
314}
315
316/* methode */
317void PAWExecutor::h_rebin(vector<string>& tokens)
318// Pour re-binner un histogramme 1D
319{
320if(tokens.size()<2)
321 {cout<<"Usage: h/rebin nameh1d nbin"<<endl; return;}
322NamedObjMgr omg;
323AnyDataObj* mobj = omg.GetObj(tokens[0]);
324if(mobj==NULL)
325 {cout<<"PAWExecutor::h_rebin Error: unknow object"<<tokens[0]<<endl;
326 return;}
327int nbin = atoi(tokens[1].c_str());
328Histo* h1 = dynamic_cast<Histo*>(mobj);
329// Ca marche aussi pour les HProf, HRebin a ete passe virtuel
330//HProf* hp = dynamic_cast<HProf*>(mobj); if(hp || !h1)
331if(!h1)
332 {cout<<"PAWExecutor::h_rebin Error: "<<tokens[0]<<" not an Histo"<<endl;
333 return;}
334h1->HRebin(nbin);
335}
336
337/* methode */
338void PAWExecutor::h_cadd(vector<string>& tokens)
339// Additionne une constante a un histogramme
340{
341if(tokens.size()<2)
342 {cout<<"Usage: h/cadd nameh1d val"<<endl; return;}
343NamedObjMgr omg;
344AnyDataObj* mobj = omg.GetObj(tokens[0]);
345if(mobj==NULL)
346 {cout<<"PAWExecutor::h_cadd Error: unknow object"<<tokens[0]<<endl;
347 return;}
348double val = atof(tokens[1].c_str());
349Histo* h1 = dynamic_cast<Histo*>(mobj);
350HProf* hp = dynamic_cast<HProf*>(mobj);
351Histo2D* h2 = dynamic_cast<Histo2D*>(mobj);
352if(h1 && !hp) *h1 += val;
353else if(h2) *h2 += val;
354else cout<<"PAWExecutor::h_cadd Error: "<<tokens[0]<<" not an Histo/Histo2D"<<endl;
355}
356
357/* methode */
358void PAWExecutor::h_cmult(vector<string>& tokens)
359// Multiplie un histogramme par une constante
360{
361if(tokens.size()<2)
362 {cout<<"Usage: h/cmult nameh1d val"<<endl; return;}
363NamedObjMgr omg;
364AnyDataObj* mobj = omg.GetObj(tokens[0]);
365if(mobj==NULL)
366 {cout<<"PAWExecutor::h_cmult Error: unknow object"<<tokens[0]<<endl;
367 return;}
368double val = atof(tokens[1].c_str());
369Histo* h1 = dynamic_cast<Histo*>(mobj);
370HProf* hp = dynamic_cast<HProf*>(mobj);
371Histo2D* h2 = dynamic_cast<Histo2D*>(mobj);
372if(h1 && !hp) *h1 *= val;
373else if(h2) *h2 *= val;
374else cout<<"PAWExecutor::h_cmult Error: "<<tokens[0]
375 <<" not an Histo/Histo2D"<<endl;
376}
377
378/* methode */
379void PAWExecutor::h_oper(vector<string>& tokens)
380// Equivalent h/oper/add sub,mul,div de paw
381// Operation entre 2 histogrammes
382// h/oper @ h1 h2 hres
383// hres = h1 @ h2 with @ = (+,-,*,/)
384{
385if(tokens.size()<4)
386 {cout<<"Usage: n/oper @ h1 h2 hres with @=(+,-,*,/)"<<endl;
387 return;}
388
389// Decode arguments
390const char * oper = tokens[0].c_str();
391if( oper[0]!='+' && oper[0]!='-' && oper[0]!='*' && oper[0]!='/' )
392 {cout<<"PAWExecutor::h_oper Error: unknow operation "<<oper<<endl;
393 return;}
394string h1name = tokens[1];
395string h2name = tokens[2];
396string h3name = tokens[3];
397
398// Get objects
399NamedObjMgr omg;
400AnyDataObj* mobjh1 = omg.GetObj(h1name);
401AnyDataObj* mobjh2 = omg.GetObj(h2name);
402if( mobjh1==NULL || mobjh2==NULL )
403 {cout<<"PAWExecutor::h_oper Error: unknow object(s) "<<h1name<<" or "<<h2name<<endl;
404 return;}
405AnyDataObj* mobjh3 = omg.GetObj(h3name);
406
407// Operations on HProf, only + is working
408if( dynamic_cast<HProf*>(mobjh1) != NULL ) {
409 if( oper[0]!='+' )
410 { cout<<"PAWExecutor::h_oper Error: operation "<<oper
411 <<" not implemented for HProf"<<endl; return;}
412 if( dynamic_cast<HProf*>(mobjh2) == NULL )
413 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h2\n"
414 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh2).name()<<endl;
415 return;}
416 HProf* h1 =(HProf*) mobjh1;
417 HProf* h2 =(HProf*) mobjh2;
418 if( h1->NBins() != h2->NBins() )
419 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h2 "
420 <<h1->NBins()<<" "<<h2->NBins()<<endl; return;}
421 HProf* h3 = NULL;
422 if( mobjh3 == NULL ) { // l'objet n'existe pas, on le cree
423 h3 = new HProf(*h1); h3->Zero(); omg.AddObj(h3,h3name);
424 } else { // l'objet existe
425 h3 = dynamic_cast<HProf*>(mobjh3);
426 if(h3 == NULL) // ce n'est pas un HProf
427 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h3\n"
428 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh3).name()<<endl;
429 return;}
430 if(h1->NBins() != h3->NBins())
431 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h3 "
432 <<h1->NBins()<<" "<<h3->NBins()<<endl; return;}
433 *h3 = *h1 + *h2;
434 h3->UpdateHisto();
435 }
436
437// Operations on Histo
438} else if( dynamic_cast<Histo*>(mobjh1) != NULL ) {
439 if( dynamic_cast<Histo*>(mobjh2) == NULL )
440 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h2\n"
441 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh2).name()<<endl;
442 return;}
443 Histo* h1 =(Histo*) mobjh1;
444 Histo* h2 =(Histo*) mobjh2;
445 if( h1->NBins() != h2->NBins() )
446 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h2 "
447 <<h1->NBins()<<" "<<h2->NBins()<<endl; return;}
448 Histo* h3 = NULL;
449 if( mobjh3 == NULL ) { // l'objet n'existe pas, on le cree
450 h3 = new Histo(*h1); h3->Zero(); omg.AddObj(h3,h3name);
451 } else { // l'objet existe
452 h3 = dynamic_cast<Histo*>(mobjh3);
453 if(h3 == NULL) // ce n'est pas un Histo
454 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h3\n"
455 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh3).name()<<endl;
456 return;}
457 if(h1->NBins() != h3->NBins())
458 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h3 "
459 <<h1->NBins()<<" "<<h3->NBins()<<endl; return;}
460 if( oper[0]=='+') *h3 = *h1 + *h2;
461 else if( oper[0]=='-') *h3 = *h1 - *h2;
462 else if( oper[0]=='*') *h3 = *h1 * *h2;
463 else if( oper[0]=='/') *h3 = *h1 / *h2;
464 }
465
466// Operations on Histo2D
467} else if( dynamic_cast<Histo2D*>(mobjh1) != NULL ) {
468 if( dynamic_cast<Histo2D*>(mobjh2) == NULL )
469 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h2\n"
470 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh2).name()<<endl;
471 return;}
472 Histo2D* h1 =(Histo2D*) mobjh1;
473 Histo2D* h2 =(Histo2D*) mobjh2;
474 if( h1->NBinX() != h2->NBinX() || h1->NBinY() != h2->NBinY() )
475 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h2 "
476 <<h1->NBinX()<<","<<h1->NBinY()<<" "
477 <<h2->NBinX()<<","<<h2->NBinY()<<endl; return;}
478 Histo2D* h3 = NULL;
479 if( mobjh3 == NULL ) { // l'objet n'existe pas, on le cree
480 h3 = new Histo2D(*h1); h3->Zero(); omg.AddObj(h3,h3name);
481 } else { // l'objet existe
482 h3 = dynamic_cast<Histo2D*>(mobjh3);
483 if(h3 == NULL) // ce n'est pas un Histo2D
484 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h3\n"
485 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh3).name()<<endl;
486 return;}
487 if( h1->NBinX() != h3->NBinX() || h1->NBinY() != h3->NBinY() )
488 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h3 "
489 <<h1->NBinX()<<","<<h1->NBinY()<<" "
490 <<h3->NBinX()<<","<<h3->NBinY()<<endl; return;}
491 if( oper[0]=='+') *h3 = *h1 + *h2;
492 else if( oper[0]=='-') *h3 = *h1 - *h2;
493 else if( oper[0]=='*') *h3 = *h1 * *h2;
494 else if( oper[0]=='/') *h3 = *h1 / *h2;
495 }
496
497// Doesn't work for other objects
498} else {
499 cout<<"PAWExecutor::h_oper Error: not implemented for "
500 <<typeid(*mobjh1).name()<<endl;
501 return;
502}
503
504return;
505}
506
507/* methode */
508void PAWExecutor::h_plot_2d(vector<string>& tokens)
509// plot for 2D histogramme: plot histo, bandx/y, slicex/y or projx/y
510{
511if(tokens.size()<2)
512 {cout<<"Usage: h/plot/2d nameh2d to_plot [n/s] [dopt]"<<endl; return;}
513NamedObjMgr omg;
514AnyDataObj* mobj = omg.GetObj(tokens[0]);
515if(mobj==NULL)
516 {cout<<"PAWExecutor::h_plot_2d Error: unknow object"<<tokens[0]<<endl;
517 return;}
518Histo2D* h2 = dynamic_cast<Histo2D*>(mobj);
519if(!h2)
520 {cout<<"PAWExecutor::h_plot_2d Error: "<<tokens[0]<<" not an Histo2D"<<endl;
521 return;}
522
523Histo* h1p = NULL; string nametoplot = "/autoc/h_plot_2d_h1";
524AnyDataObj* mobjh1 = omg.GetObj(nametoplot);
525// if(mobjh1!=NULL) omg.DelObj(nametoplot); $CHECK$ Reza 12/12/99 - Ne pas faire - autoclean
526
527string dopt = ""; if(tokens.size()>=3) dopt = tokens[2];
528if(tokens[1] == "show") {
529 h2->ShowProj();
530 h2->ShowBand(2);
531 h2->ShowSli(2);
532 return;
533} else if(tokens[1] == "h") {
534 nametoplot = tokens[0];
535} else if(tokens[1] == "px") {
536 if((h1p=h2->HProjX())) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
537 else {h2->ShowProj(); return;}
538} else if(tokens[1] == "py") {
539 if((h1p=h2->HProjY())) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
540 else {h2->ShowProj(); return;}
541} else {
542 if(tokens.size()<3)
543 {cout<<"Usage: h/plot/2d nameh2d bx/by/sx/sy n [dopt]"<<endl; return;}
544 int n = atoi(tokens[2].c_str());
545 dopt = ""; if(tokens.size()>=4) dopt = tokens[3];
546 if(tokens[1] == "bx") {
547 if((h1p=h2->HBandX(n))) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
548 else {h2->ShowBand(); return;}
549 } else if(tokens[1] == "by") {
550 if((h1p=h2->HBandY(n))) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
551 else {h2->ShowBand(); return;}
552 } else if(tokens[1] == "sx") {
553 if((h1p=h2->HSliX(n))) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
554 else {h2->ShowSli(); return;}
555 } else if(tokens[1] == "sy") {
556 if((h1p=h2->HSliY(n))) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
557 else {h2->ShowSli(); return;}
558 }
559}
560
561omg.DisplayObj(nametoplot,dopt);
562}
563
564/* methode */
565int PAWExecutor::decodepawstring(string tokens,string& nameobj
566 ,string& xexp,string& yexp,string& zexp)
567// Decodage general de "nameobj.xexp"
568// "nameobj.yexp%xexp"
569// "nameobj.zexp%yexp%xexp"
570// Return: nombre de variables trouvees, -1 si probleme
571{
572nameobj = ""; xexp= ""; yexp= ""; zexp= "";
573
574int lt = (int) tokens.length();
575if(lt<=0) return -1;
576
577// decodage de la chaine de type PAW.
578char *str = new char[lt+2];
579strcpy(str,tokens.c_str()); strip(str,'B',' '); lt = strlen(str);
580//cout<<"chaine1["<<lt<<"] :"<<str<<":"<<endl;
581char *c[3] = {NULL,NULL,NULL};
582int i, np=0; bool namefound = false;
583for(i=0;i<lt;i++) {
584 if(!namefound && str[i]=='.') {
585 str[i]='\0';
586 namefound=true;
587 c[np] = str+i+1; np++;
588 }
589 if( namefound && str[i]=='%') {
590 str[i]='\0';
591 if(np<3) {c[np] = str+i+1; np++;}
592 }
593}
594//cout<<"chaine2 :"; for(i=0;i<lt;i++) cout<<str[i]; cout<<":"<<endl;
595
596// Remplissage du nom et des variables
597nameobj = str;
598if(np==1) xexp=c[0];
599if(np==2) {yexp=c[0]; xexp=c[1];}
600if(np==3) {zexp=c[0]; yexp=c[1]; xexp=c[2];}
601//cout<<"pawstring str,c[0-2] "<<str<<" "<<c[0]<<" "<<c[1]<<" "<<c[2]<<endl;
602delete [] str;
603
604// Comptage des variables
605np = -1;
606if(nameobj.length()>0)
607 {np = 0; if(xexp.length()>0)
608 {np++; if(yexp.length()>0)
609 {np++; if(zexp.length()>0) np++;}}}
610cout<<"pawstring["<<np<<"] name="<<nameobj
611 <<" xexp="<<xexp<<" yexp="<<yexp<<" zexp="<<zexp<<endl;
612return np;
613}
Note: See TracBrowser for help on using the repository browser.