molecule.C File Reference

#include <iostream>
#include "TROOT.h"
#include "../dnadamage1/include/DNAVolumeType.hh"

Go to the source code of this file.

Classes
class	ThreeVector< T >
class	Molecule

Functions
	ClassImp (ThreeVector)
	ClassImp (Molecule)
std::vector< Molecule >	molecule ()

Function Documentation

ClassImp

(

ThreeVector

)

Referenced by ThreeVector< double >::operator==().

ClassImp

(

Molecule

)

std::vector<Molecule> molecule

(

)

Definition at line 87 of file molecule.C.

References file, fMolecules, util::size(), x, y, and z.

{
    std::vector<Molecule> fMolecules;
    double size;
    string name;
    ifstream file("VoxelStraight.fab2g4dna");
    if(!file.is_open())
    {
        string msg ="VoxelStraight.fab2g4dna could not be opened";
        throw std::invalid_argument(msg);
    }

   string line;
    while(getline(file, line) )
    {
        if(line.empty()) 
        {
            continue;
        }
         
        istringstream issLine(line);
        string firstItem;
        issLine >> firstItem;
        if("_Size" == firstItem)
        {
            issLine >> size;
        }
        else if("_pl" == firstItem)
        {
            string name;
            issLine >> name;

            string material;
            issLine >> material;

            int strand;
            issLine >> strand;

            int copyNumber;
            issLine >> copyNumber;

            double x;
            issLine >> x;

            double y;
            issLine >> y;

            double z;
            issLine >> z;

            Molecule molecule(name, 
                              copyNumber, 
                              ThreeVector<double>(x, y, z), 
                              strand);
            fMolecules.push_back(molecule);
        }
    }
    file.close();
    
    return fMolecules;
}