TPC/SCDCorrections/genRho.cc

0001 #include <iostream>
0002 #include <fstream>
0003 #include <string>
0004 #include "QPileUpToy.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;
0012   int kNAzimuthalSteps;
0013   int kNLongitudinalSteps;
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.close();
0022   kNRadialSteps = (kOutterRadius-kInnerRadius)*10;
0023   kNAzimuthalSteps = 18;
0024   kNLongitudinalSteps = kHalfLength*20;
0025
0026
0027   QPileUp *initialdensity;
0028   initialdensity = new QPileUpToy(1.0/76628.0/*gasfactor [Vs]*/, 425.0 /*multiplicity*/, 5e+4/*rate [Hz]*/, 6); // PHENIX
0029   //initialdensity = new QPileUpToy(1.0/76628.0/*gasfactor [Vs]*/, 950.0 /*multiplicity*/, 5e+4/*rate [Hz]*/, 20, 1.5); // ALICE
0030   //initialdensity = new QPileUpToy(1.0/76628.0/*gasfactor [Vs]*/, 425.0 /*multiplicity*/, 15e+3/*rate [Hz]*/, 0); // STAR
0031   initialdensity->SetDebugLevel(2);
0032   initialdensity->OutputFileName(kFileNameRoot);
0033   initialdensity->TPCDimensions( kInnerRadius, kOutterRadius, kHalfLength );
0034   initialdensity->TPCGridSize( kNRadialSteps, kNAzimuthalSteps, kNLongitudinalSteps );
0035   initialdensity->Make();
0036
0037   delete initialdensity;
0038   return 0;
0039 }