Changeset 400 in svn for trunk/Utilities/Hector/src
- Timestamp:
- May 20, 2009, 10:33:21 AM (15 years ago)
- Location:
- trunk/Utilities/Hector/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Utilities/Hector/src/H_BeamParticle.cc
r388 r400 461 461 } 462 462 463 void H_BeamParticle::getPath(const int x_or_y, const string& filename) const { 463 void H_BeamParticle::getPath(const int x_or_y, const string& filename) const{ 464 /// @param x_or_y = 0(1) draws the x(y) component; 464 465 ofstream outfile(filename.c_str()); 465 int index; 466 467 int index; 466 468 if(x_or_y==0) {index = INDEX_X;} else {index = INDEX_Y;} 467 469 468 470 vector<TVectorD>::const_iterator position_i; 469 471 for(position_i = positions.begin(); position_i < positions.end(); position_i++) { 470 outfile << (*position_i)[INDEX_S] <<"\t\t" << (*position_i)[index] << endl;472 outfile << (*position_i)[index] << "\t" << (*position_i)[INDEX_S] << endl; 471 473 } 472 474 outfile.close(); 473 475 } 476 474 477 475 478 // should be removed later, to keep only computePath(const H_AbstractBeamLine & , const bool) -
trunk/Utilities/Hector/src/H_RecRPObject.cc
r281 r400 64 64 {if(ss1==ss2) cout<<"<H_RecRPObject> WARNING : detectors are on same position"<<endl; 65 65 } 66 67 68 H_RecRPObject::H_RecRPObject(const float ss1, const float ss2, const H_AbstractBeamLine& beam) : emin(0), emax(-1), x1(0), x2(0), y1(0), y2(0), s1(ss1), s2(ss2), 69 txip(NOT_YET_COMPUTED), tyip(NOT_YET_COMPUTED), energy(NOT_YET_COMPUTED), q2(NOT_YET_COMPUTED), pt(NOT_YET_COMPUTED), 70 thebeam(beam.clone()), 71 f_1(new TF1("f_1","[0] + [1]*x + [2]*x*x ",emin,emax)), 72 f_2(new TF1("f_2","[0] + [1]*x + [2]*x*x ",emin,emax)), 73 g_1(new TF1("g_1","[0] + [1]*x + [2]*x*x ",emin,emax)), 74 g_2(new TF1("g_2","[0] + [1]*x + [2]*x*x ",emin,emax)), 75 d_1(new TF1("d_1","[0] + [1]*x + [2]*x*x ",emin,emax)), 76 d_2(new TF1("d_2","[0] + [1]*x + [2]*x*x ",emin,emax)), 77 k_1(new TF1("k_1","[0] + [1]*x + [2]*x*x ",emin,emax)), 78 k_2(new TF1("k_2","[0] + [1]*x + [2]*x*x ",emin,emax)), 79 l_1(new TF1("l_1","[0] + [1]*x + [2]*x*x ",emin,emax)), 80 l_2(new TF1("l_2","[0] + [1]*x + [2]*x*x ",emin,emax)) 81 {if(ss1==ss2) cout<<"<H_RecRPObject> WARNING : detectors are on same position"<<endl; 82 } 83 66 84 67 85 H_RecRPObject::H_RecRPObject(const H_RecRPObject& r): … … 351 369 return os; 352 370 } 371
Note:
See TracChangeset
for help on using the changeset viewer.