Changeset 388 in svn

Timestamp:

May 16, 2009, 10:31:14 AM (16 years ago)

Author:

Xavier Rouby

Message:

modif Xav

Location:

trunk/Utilities/Hector

Files:

• include/H_BeamParticle.h (modified) (1 diff)
• src/H_BeamParticle.cc (modified) (2 diffs)

trunk/Utilities/Hector/include/H_BeamParticle.h

@@ -132 +132 @@
                 /// Returns the particle path in the beamline
                 TGraph * getPath(const int , const int ) const;
+                void getPath(const int , const string& ) const;
                 /// Computes the position of the particle at the end of each element of the beam.
                 void computePath(const H_AbstractBeamLine * beam, const bool NonLinear=true);

trunk/Utilities/Hector/src/H_BeamParticle.cc

@@ -26 +26 @@
 // c++ #includes
 #include <iostream>
+#include <fstream>
 #include <iomanip>
 #include <vector>
@@ -458 +459 @@
         delete [] graph;
         return ppath;
+}
+
+void H_BeamParticle::getPath(const int x_or_y, const string& filename) const {
+        ofstream outfile(filename.c_str());
+        int index;
+        if(x_or_y==0) {index = INDEX_X;} else {index = INDEX_Y;}
+
+        vector<TVectorD>::const_iterator position_i;
+        for(position_i = positions.begin(); position_i < positions.end(); position_i++) {
+                outfile << (*position_i)[INDEX_S] <<"\t\t" << (*position_i)[index] << endl;
+        }
+        outfile.close();
 }