source: ZHANGProjects/ICOSIM/CPP/trunk/source/testLattice.cc @ 2

#include <iostream>
#include <vector>
#include <string>
#include "lattice.h"
using namespace std;


/*
####################################################################

Test file for the class Lattice.

####################################################################
*/

int main()
{

    vector <int> ips;

    Element elt1(-2.316794521, 0.537534425, 0.0712, 0.01725, 0, 0, "RECTANGLE", 103.4788975, 20.8994661, 0.027602708, 0, 1.231007397, 0, "MARKER", 0, 0, 0, "BEGI.10010", 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, "x");

    Element elt2(2.316066957, -0.534433046, 0.0712, 0.01725, 0, 0, "RECTANGLE", 103.4812508, 20.8903137, -0.027848221, 0, 1.230624298, 0, "QUADRUPOLE", 3.085, 0.004637558, 0.024132382, "QF.10010", 0, 0, 0, 3.085, 0, 0, 0, 0, 0, 0.044535174, 0, 0, 0, "x");

    StandardCollimator stdcolli(0.673651882, -2.183766647, 0.05, 0.05, 0, 0, "RECTANGLE", 23.65235333, 93.32238148, 0.021383636, 0, -0.338555786, 0, "RCOLLIMATOR", 0.6, 19.69294239, 19.76251791, "BRCV.51899", 0, 0, 0, 5213.4097, 0, 0, 0, 0, 0, 0, 0, 0, 0, "x", "standard", 0.002, 0.001, 0.1, 0.1, 0.2, 0, 0);

    FlukaCollimator flukacolli(1.288591765, -1.577339513, 0.078, 0.078, 0, 0, "ELLIPSE", 43.26312194, 56.37282712, -0.005548392, 0, -0.242857529, 0, "RCOLLIMATOR", 2, 1.980859479, 2.044620435, "TBSM.11672", 0, 0, 0, 531.7442, 0, 0, 0, 0, 0, 0, 0, 0, 0, "x", "fluka", 0.002, 0.001, 0.1, 0, 0.2, 0, 0);

    MagneticCollimator magnetcolli(0.735310344, -2.131043701, 0.0425, 0.0204, 0, 0, "RECTANGLE", 24.89598243, 90.02483262, 0.021833793, 0, -0.367904763, 0, "RCOLLIMATOR", 4.3, 2.266614904, 2.306252605, "TIDV.11892", 0, 0, 0, 604.9077, 0, 0, 0, 0, 0, 0, 0, 0, 0, "x", "magnetic", 0.002, 0.001, 0.1, 0.1, 0.2, 0, 0, 2, 0.0001, 574000, 193.6872);

    ips.push_back(2);
    ips.push_back(3);
    ips.push_back(4);

    Lattice lat(elt1, elt2, stdcolli, flukacolli, magnetcolli, ips);

    cout << "On affiche s: " << lat.reseau[0]->S << ", " << lat.reseau[1]->S << endl;

    cout << "On a " << lat.npart << " particules." << endl;

    Particle p1(1, 1, 1, 1, 1, 1, 1, 1, 1);
    Particle p2;
    Particle pref;

    //lat.trackensemblechrom(0);

    return 0;
}