Fork me on GitHub

source: svn/trunk/Utilities/HepMC/src/PdfInfo.cc@ 606

Last change on this file since 606 was 571, checked in by cp3-support, 13 years ago

upgrade HepMC to version 2.06.05

File size: 3.6 KB
Line 
1//--------------------------------------------------------------------------
2//
3// PdfInfo.cc
4// Author: Lynn Garren
5//
6// Implement operator >> and operator <<
7//
8// ----------------------------------------------------------------------
9
10#include <iostream>
11#include <ostream>
12#include <istream>
13#include <sstream>
14
15#include "PdfInfo.h"
16#include "StreamHelpers.h"
17#include "IO_Exception.h"
18
19namespace HepMC {
20
21std::ostream & operator << ( std::ostream & os, PdfInfo const * pdf)
22{
23 if ( !os ) {
24 std::cerr << "operator << for PdfInfo: !os, "
25 << " setting badbit" << std::endl;
26 os.clear(std::ios::badbit);
27 return os;
28 }
29 os << 'F';
30 // PdfInfo* is set to 0 by default
31 if ( !pdf ) {
32 detail::output( os, 0 );
33 detail::output( os, 0 );
34 detail::output( os, 0. );
35 detail::output( os, 0. );
36 detail::output( os, 0. );
37 detail::output( os, 0. );
38 detail::output( os, 0. );
39 detail::output( os, 0 );
40 detail::output( os, 0 );
41 detail::output( os,'\n');
42 return os;
43 }
44 //
45 detail::output( os, pdf->id1() );
46 detail::output( os, pdf->id2() );
47 detail::output( os, pdf->x1() );
48 detail::output( os, pdf->x2() );
49 detail::output( os, pdf->scalePDF() );
50 detail::output( os, pdf->pdf1() );
51 detail::output( os, pdf->pdf2() );
52 detail::output( os, pdf->pdf_id1() );
53 detail::output( os, pdf->pdf_id2() );
54 detail::output( os,'\n');
55
56 return os;
57}
58
59std::istream & operator >> (std::istream & is, PdfInfo * pdf)
60{
61 // make sure the stream is valid
62 if ( !is ) {
63 std::cerr << "PdfInfo input stream setting badbit." << std::endl;
64 is.clear(std::ios::badbit);
65 return is;
66 }
67 //
68 // get the PdfInfo line
69 std::string line;
70 std::getline(is,line);
71 std::istringstream iline(line);
72 std::string firstc;
73 iline >> firstc;
74 // test to be sure the next entry is of type "F" then ignore it
75 if ( firstc != "F" ) {
76 std::cerr << "PdfInfo input stream invalid line type: "
77 << firstc << std::endl;
78 // this is non-recoverable, so throw here
79 throw IO_Exception("PdfInfo input stream encounterd invalid data");
80 }
81 // read values into temp variables, then create a new PdfInfo object
82 int id1 =0, id2 =0, pdf_id1=0, pdf_id2=0;
83 double x1 = 0., x2 = 0., scale = 0., pdf1 = 0., pdf2 = 0.;
84 iline >> id1 ;
85 if(!iline) throw IO_Exception("PdfInfo input stream encounterd invalid data");
86 // check now for empty PdfInfo line
87 if( id1 == 0 ) return is;
88 // continue reading
89 iline >> id2 ;
90 if(!iline) throw IO_Exception("PdfInfo input stream encounterd invalid data");
91 iline >> x1 ;
92 if(!iline) throw IO_Exception("PdfInfo input stream encounterd invalid data");
93 iline >> x2 ;
94 if(!iline) throw IO_Exception("PdfInfo input stream encounterd invalid data");
95 iline >> scale ;
96 if(!iline) throw IO_Exception("PdfInfo input stream encounterd invalid data");
97 iline >> pdf1 ;
98 if(!iline) throw IO_Exception("PdfInfo input stream encounterd invalid data");
99 iline >> pdf2;
100 if(!iline) throw IO_Exception("PdfInfo input stream encounterd invalid data");
101 // check to see if we are at the end of the line
102 if( !iline.eof() ) {
103 iline >> pdf_id1 ;
104 if(!iline) throw IO_Exception("PdfInfo input stream encounterd invalid data");
105 iline >> pdf_id2;
106 if(!iline) throw IO_Exception("PdfInfo input stream encounterd invalid data");
107 }
108 pdf->set_id1( id1 );
109 pdf->set_id2( id2 );
110 pdf->set_pdf_id1( pdf_id1 );
111 pdf->set_pdf_id2( pdf_id2 );
112 pdf->set_x1( x1 );
113 pdf->set_x2( x2 );
114 pdf->set_scalePDF( scale );
115 pdf->set_pdf1( pdf1 );
116 pdf->set_pdf2( pdf2 );
117
118 return is;
119}
120
121} // HepMC
Note: See TracBrowser for help on using the repository browser.