JPhotonPathIO.cc File Reference

#include <iostream>
#include <fstream>
#include "TRandom3.h"
#include "JMarkov/JPhotonPath.hh"
#include "JMarkov/JPhotonPathReader.hh"
#include "JMarkov/JPhotonPathWriter.hh"

Go to the source code of this file.

Functions
void	printPolyline (JPolyline3D &pl)
int	main (int argc, char **argv)

Function Documentation

printPolyline()

void printPolyline

(

JPolyline3D &

pl

)

Definition at line 13 of file JPhotonPathIO.cc.

                                      {
  for( int i=0 ; i<(int)pl.size()-1 ; ++i ) 
    cout << "(" << pl[i] << ")-" ;
  if( pl.size()>0 ) cout << "(" << pl.back() << ")" ;
}

main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 25 of file JPhotonPathIO.cc.

                                  {
  // initialization
  TRandom* ran = new TRandom3(0) ;
  const int npaths = 5 ;
  char fname[100] = "out.paths" ;
 
  // write some example paths to a file
  cout << "Writing some example paths to '" << fname << "'." << endl ;
 
  JMARKOV::JPhotonPathWriter writer ;
  writer.open(fname) ;
  for( int n=0 ; n<npaths ; ++n ) {
    // create path
    int nscat = ran->Integer(6) ; // how often the photon scatters
    JPhotonPath path(nscat) ;
    // loop over the vertices
    for( JPhotonPath::iterator it=path.begin() ; it!=path.end() ; ++it ) {
      double x = ran->Integer(50) ;
      double y = ran->Integer(50) ;
      double z = ran->Integer(50) ;
      it->setPosition( JGEOMETRY3D::JVector3D(x,y,z) ) ;
    }
    // print the path
    printPolyline( path ) ;
    cout << endl ;
    // write the path
    writer.put(path) ;
  }
  writer.close() ;
  cout << endl ;
 
  // read paths from a file
  JMARKOV::JPhotonPathReader reader ;
  reader.open("out.paths") ;
 
  int nread = 0 ;
  JPhotonPath* p = NULL ;
  cout << "Reading file" << endl ;
  while( reader.hasNext() ) {
    p = reader.next() ;
    printPolyline( *p ) ;
    cout << endl ;
    ++nread ;
  }
 
  cout << "Done reading file. Read " << nread << " paths from it." << endl ;
  cout << endl ;
  cout << "Done!" << endl ;
 
  return 0 ;
}