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

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

cmv 7/7/2000

nomhistadapter.cc :

GetCopyObj() on fait un UpdateHisto() si c'est un HProf

servnobjm.cc :

ProjectHProf() on fait un UpdateHisto()

pihisto.cc :

Draw() on fait un UpdateHisto() si c'est un HProf

pawexecut.cc :

h/integ, h/deriv, h/rebin, h/cadd, h/cmult NON valable pour un HProf
h/oper teste et h3->Zero() ENLEVE pour permettre

h/oper @ h1 h3 h3
h/oper @ h3 h3 h3 ... etc...

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);
329HProf* hp = dynamic_cast<HProf*>(mobj);
330if(hp || !h1)
331 {cout<<"PAWExecutor::h_rebin Error: "<<tokens[0]<<" not an Histo"<<endl;
332 return;}
333h1->HRebin(nbin);
334}
335
336/* methode */
337void PAWExecutor::h_cadd(vector<string>& tokens)
338// Additionne une constante a un histogramme
339{
340if(tokens.size()<2)
341 {cout<<"Usage: h/cadd nameh1d val"<<endl; return;}
342NamedObjMgr omg;
343AnyDataObj* mobj = omg.GetObj(tokens[0]);
344if(mobj==NULL)
345 {cout<<"PAWExecutor::h_cadd Error: unknow object"<<tokens[0]<<endl;
346 return;}
347double val = atof(tokens[1].c_str());
348Histo* h1 = dynamic_cast<Histo*>(mobj);
349HProf* hp = dynamic_cast<HProf*>(mobj);
350Histo2D* h2 = dynamic_cast<Histo2D*>(mobj);
351if(h1 && !hp) *h1 += val;
352else if(h2) *h2 += val;
353else cout<<"PAWExecutor::h_cadd Error: "<<tokens[0]<<" not an Histo/Histo2D"<<endl;
354}
355
356/* methode */
357void PAWExecutor::h_cmult(vector<string>& tokens)
358// Multiplie un histogramme par une constante
359{
360if(tokens.size()<2)
361 {cout<<"Usage: h/cmult nameh1d val"<<endl; return;}
362NamedObjMgr omg;
363AnyDataObj* mobj = omg.GetObj(tokens[0]);
364if(mobj==NULL)
365 {cout<<"PAWExecutor::h_cmult Error: unknow object"<<tokens[0]<<endl;
366 return;}
367double val = atof(tokens[1].c_str());
368Histo* h1 = dynamic_cast<Histo*>(mobj);
369HProf* hp = dynamic_cast<HProf*>(mobj);
370Histo2D* h2 = dynamic_cast<Histo2D*>(mobj);
371if(h1 && !hp) *h1 *= val;
372else if(h2) *h2 *= val;
373else cout<<"PAWExecutor::h_cmult Error: "<<tokens[0]
374 <<" not an Histo/Histo2D"<<endl;
375}
376
377/* methode */
378void PAWExecutor::h_oper(vector<string>& tokens)
379// Equivalent h/oper/add sub,mul,div de paw
380// Operation entre 2 histogrammes
381// h/oper @ h1 h2 hres
382// hres = h1 @ h2 with @ = (+,-,*,/)
383{
384if(tokens.size()<4)
385 {cout<<"Usage: n/oper @ h1 h2 hres with @=(+,-,*,/)"<<endl;
386 return;}
387
388// Decode arguments
389const char * oper = tokens[0].c_str();
390if( oper[0]!='+' && oper[0]!='-' && oper[0]!='*' && oper[0]!='/' )
391 {cout<<"PAWExecutor::h_oper Error: unknow operation "<<oper<<endl;
392 return;}
393string h1name = tokens[1];
394string h2name = tokens[2];
395string h3name = tokens[3];
396
397// Get objects
398NamedObjMgr omg;
399AnyDataObj* mobjh1 = omg.GetObj(h1name);
400AnyDataObj* mobjh2 = omg.GetObj(h2name);
401if( mobjh1==NULL || mobjh2==NULL )
402 {cout<<"PAWExecutor::h_oper Error: unknow object(s) "<<h1name<<" or "<<h2name<<endl;
403 return;}
404AnyDataObj* mobjh3 = omg.GetObj(h3name);
405
406// Operations on HProf, only + is working
407if( dynamic_cast<HProf*>(mobjh1) != NULL ) {
408 if( oper[0]!='+' )
409 { cout<<"PAWExecutor::h_oper Error: operation "<<oper
410 <<" not implemented for HProf"<<endl; return;}
411 if( dynamic_cast<HProf*>(mobjh2) == NULL )
412 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h2\n"
413 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh2).name()<<endl;
414 return;}
415 HProf* h1 =(HProf*) mobjh1;
416 HProf* h2 =(HProf*) mobjh2;
417 if( h1->NBins() != h2->NBins() )
418 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h2 "
419 <<h1->NBins()<<" "<<h2->NBins()<<endl; return;}
420 HProf* h3 = NULL;
421 if( mobjh3 == NULL ) { // l'objet n'existe pas, on le cree
422 h3 = new HProf(*h1); h3->Zero(); omg.AddObj(h3,h3name);
423 } else { // l'objet existe
424 h3 = dynamic_cast<HProf*>(mobjh3);
425 if(h3 == NULL) // ce n'est pas un HProf
426 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h3\n"
427 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh3).name()<<endl;
428 return;}
429 if(h1->NBins() != h3->NBins())
430 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h3 "
431 <<h1->NBins()<<" "<<h3->NBins()<<endl; return;}
432 *h3 = *h1 + *h2;
433 h3->UpdateHisto();
434 }
435
436// Operations on Histo
437} else if( dynamic_cast<Histo*>(mobjh1) != NULL ) {
438 if( dynamic_cast<Histo*>(mobjh2) == NULL )
439 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h2\n"
440 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh2).name()<<endl;
441 return;}
442 Histo* h1 =(Histo*) mobjh1;
443 Histo* h2 =(Histo*) mobjh2;
444 if( h1->NBins() != h2->NBins() )
445 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h2 "
446 <<h1->NBins()<<" "<<h2->NBins()<<endl; return;}
447 Histo* h3 = NULL;
448 if( mobjh3 == NULL ) { // l'objet n'existe pas, on le cree
449 h3 = new Histo(*h1); h3->Zero(); omg.AddObj(h3,h3name);
450 } else { // l'objet existe
451 h3 = dynamic_cast<Histo*>(mobjh3);
452 if(h3 == NULL) // ce n'est pas un Histo
453 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h3\n"
454 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh3).name()<<endl;
455 return;}
456 if(h1->NBins() != h3->NBins())
457 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h3 "
458 <<h1->NBins()<<" "<<h3->NBins()<<endl; return;}
459 if( oper[0]=='+') *h3 = *h1 + *h2;
460 else if( oper[0]=='-') *h3 = *h1 - *h2;
461 else if( oper[0]=='*') *h3 = *h1 * *h2;
462 else if( oper[0]=='/') *h3 = *h1 / *h2;
463 }
464
465// Operations on Histo2D
466} else if( dynamic_cast<Histo2D*>(mobjh1) != NULL ) {
467 if( dynamic_cast<Histo2D*>(mobjh2) == NULL )
468 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h2\n"
469 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh2).name()<<endl;
470 return;}
471 Histo2D* h1 =(Histo2D*) mobjh1;
472 Histo2D* h2 =(Histo2D*) mobjh2;
473 if( h1->NBinX() != h2->NBinX() || h1->NBinY() != h2->NBinY() )
474 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h2 "
475 <<h1->NBinX()<<","<<h1->NBinY()<<" "
476 <<h2->NBinX()<<","<<h2->NBinY()<<endl; return;}
477 Histo2D* h3 = NULL;
478 if( mobjh3 == NULL ) { // l'objet n'existe pas, on le cree
479 h3 = new Histo2D(*h1); h3->Zero(); omg.AddObj(h3,h3name);
480 } else { // l'objet existe
481 h3 = dynamic_cast<Histo2D*>(mobjh3);
482 if(h3 == NULL) // ce n'est pas un Histo2D
483 {cout<<"PAWExecutor::h_oper Error: type mismatch between h1 and h3\n"
484 <<typeid(*mobjh1).name()<<" , "<<typeid(*mobjh3).name()<<endl;
485 return;}
486 if( h1->NBinX() != h3->NBinX() || h1->NBinY() != h3->NBinY() )
487 {cout<<"PAWExecutor::h_oper Error: size mismatch between h1, h3 "
488 <<h1->NBinX()<<","<<h1->NBinY()<<" "
489 <<h3->NBinX()<<","<<h3->NBinY()<<endl; return;}
490 if( oper[0]=='+') *h3 = *h1 + *h2;
491 else if( oper[0]=='-') *h3 = *h1 - *h2;
492 else if( oper[0]=='*') *h3 = *h1 * *h2;
493 else if( oper[0]=='/') *h3 = *h1 / *h2;
494 }
495
496// Doesn't work for other objects
497} else {
498 cout<<"PAWExecutor::h_oper Error: not implemented for "
499 <<typeid(*mobjh1).name()<<endl;
500 return;
501}
502
503return;
504}
505
506/* methode */
507void PAWExecutor::h_plot_2d(vector<string>& tokens)
508// plot for 2D histogramme: plot histo, bandx/y, slicex/y or projx/y
509{
510if(tokens.size()<2)
511 {cout<<"Usage: h/plot/2d nameh2d to_plot [n/s] [dopt]"<<endl; return;}
512NamedObjMgr omg;
513AnyDataObj* mobj = omg.GetObj(tokens[0]);
514if(mobj==NULL)
515 {cout<<"PAWExecutor::h_plot_2d Error: unknow object"<<tokens[0]<<endl;
516 return;}
517Histo2D* h2 = dynamic_cast<Histo2D*>(mobj);
518if(!h2)
519 {cout<<"PAWExecutor::h_plot_2d Error: "<<tokens[0]<<" not an Histo2D"<<endl;
520 return;}
521
522Histo* h1p = NULL; string nametoplot = "/autoc/h_plot_2d_h1";
523AnyDataObj* mobjh1 = omg.GetObj(nametoplot);
524// if(mobjh1!=NULL) omg.DelObj(nametoplot); $CHECK$ Reza 12/12/99 - Ne pas faire - autoclean
525
526string dopt = ""; if(tokens.size()>=3) dopt = tokens[2];
527if(tokens[1] == "show") {
528 h2->ShowProj();
529 h2->ShowBand(2);
530 h2->ShowSli(2);
531 return;
532} else if(tokens[1] == "h") {
533 nametoplot = tokens[0];
534} else if(tokens[1] == "px") {
535 if((h1p=h2->HProjX())) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
536 else {h2->ShowProj(); return;}
537} else if(tokens[1] == "py") {
538 if((h1p=h2->HProjY())) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
539 else {h2->ShowProj(); return;}
540} else {
541 if(tokens.size()<3)
542 {cout<<"Usage: h/plot/2d nameh2d bx/by/sx/sy n [dopt]"<<endl; return;}
543 int n = atoi(tokens[2].c_str());
544 dopt = ""; if(tokens.size()>=4) dopt = tokens[3];
545 if(tokens[1] == "bx") {
546 if((h1p=h2->HBandX(n))) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
547 else {h2->ShowBand(); return;}
548 } else if(tokens[1] == "by") {
549 if((h1p=h2->HBandY(n))) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
550 else {h2->ShowBand(); return;}
551 } else if(tokens[1] == "sx") {
552 if((h1p=h2->HSliX(n))) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
553 else {h2->ShowSli(); return;}
554 } else if(tokens[1] == "sy") {
555 if((h1p=h2->HSliY(n))) {Histo* h1=new Histo(*h1p); omg.AddObj(h1,nametoplot);}
556 else {h2->ShowSli(); return;}
557 }
558}
559
560omg.DisplayObj(nametoplot,dopt);
561}
562
563/* methode */
564int PAWExecutor::decodepawstring(string tokens,string& nameobj
565 ,string& xexp,string& yexp,string& zexp)
566// Decodage general de "nameobj.xexp"
567// "nameobj.yexp%xexp"
568// "nameobj.zexp%yexp%xexp"
569// Return: nombre de variables trouvees, -1 si probleme
570{
571nameobj = ""; xexp= ""; yexp= ""; zexp= "";
572
573int lt = (int) tokens.length();
574if(lt<=0) return -1;
575
576// decodage de la chaine de type PAW.
577char *str = new char[lt+2];
578strcpy(str,tokens.c_str()); strip(str,'B',' '); lt = strlen(str);
579//cout<<"chaine1["<<lt<<"] :"<<str<<":"<<endl;
580char *c[3] = {NULL,NULL,NULL};
581int i, np=0; bool namefound = false;
582for(i=0;i<lt;i++) {
583 if(!namefound && str[i]=='.') {
584 str[i]='\0';
585 namefound=true;
586 c[np] = str+i+1; np++;
587 }
588 if( namefound && str[i]=='%') {
589 str[i]='\0';
590 if(np<3) {c[np] = str+i+1; np++;}
591 }
592}
593//cout<<"chaine2 :"; for(i=0;i<lt;i++) cout<<str[i]; cout<<":"<<endl;
594
595// Remplissage du nom et des variables
596nameobj = str;
597if(np==1) xexp=c[0];
598if(np==2) {yexp=c[0]; xexp=c[1];}
599if(np==3) {zexp=c[0]; yexp=c[1]; xexp=c[2];}
600//cout<<"pawstring str,c[0-2] "<<str<<" "<<c[0]<<" "<<c[1]<<" "<<c[2]<<endl;
601delete [] str;
602
603// Comptage des variables
604np = -1;
605if(nameobj.length()>0)
606 {np = 0; if(xexp.length()>0)
607 {np++; if(yexp.length()>0)
608 {np++; if(zexp.length()>0) np++;}}}
609cout<<"pawstring["<<np<<"] name="<<nameobj
610 <<" xexp="<<xexp<<" yexp="<<yexp<<" zexp="<<zexp<<endl;
611return np;
612}
Note: See TracBrowser for help on using the repository browser.