TPC/SCDCorrections/genD.cc

0001 #include <iostream>
0002 #include <fstream>
0003 #include <string>
0004 #include "Langevin.h"
0005
0006 int main() {
0007   std::string kFileNameRoot="rho";
0008   float kInnerRadius = 30; //cm
0009   float kOutterRadius = 80; //cm
0010   float kHalfLength = 80; //cm
0011   int kNRadialSteps = 50;
0012   int kNAzimuthalSteps = 1;
0013   int kNLongitudinalSteps = 160;
0014
0015   std::cout << "fetching geo.dat..." << std::endl;
0016   std::ifstream ifile("geo.dat");
0017   ifile >> kFileNameRoot;
0018   ifile >> kInnerRadius;  //cm
0019   ifile >> kOutterRadius; //cm
0020   ifile >> kHalfLength;   //cm
0021   ifile >> kNRadialSteps;
0022   ifile >> kNAzimuthalSteps;
0023   ifile >> kNLongitudinalSteps;
0024   ifile.close();
0025
0026   Langevin *langevin = new Langevin();
0027   langevin->SetDebugLevel(2);
0028   langevin->OutputFileName(kFileNameRoot);
0029   langevin->TPCDimensions( kInnerRadius, kOutterRadius, kHalfLength );
0030   langevin->TPCGridSize( kNRadialSteps, kNAzimuthalSteps, kNLongitudinalSteps );
0031   langevin->Make();
0032
0033   delete langevin;
0034   return 0;
0035 }