Changeset 312 in PSPA for Interface_Web/trunk/pspaWT/sources/controler/src/softwareGenerator.cc

Timestamp:

Feb 5, 2013, 2:23:23 PM (12 years ago)

Author:

lemeur

Message:

suppression designation elements par index + fin generator

File:

• Interface_Web/trunk/pspaWT/sources/controler/src/softwareGenerator.cc (modified) (5 diffs)

Interface_Web/trunk/pspaWT/sources/controler/src/softwareGenerator.cc

@@ -1 +1 @@
 #include "softwareGenerator.h"
 #include "abstractElement.h"
-#include "parmelaParticle.h"
+#include "generatorParticle.h"
+#include "mathematicalConstants.h"
+#include "PhysicalConstants.h"
 #include "dataManager.h"
 
@@ -21 +23 @@
 bool softwareGenerator::createInputFile(string inputFileName, particleBeam* beamBefore, unsigned int numeroDeb, unsigned int numeroFin, string workingDir)
 {
-
   if ( numeroDeb != numeroFin ) return false;
 
   abstractElement* elPtr;
-  elPtr = dataManager_->getElementPointerFromIndex(numeroDeb);
+  elPtr = dataManager_->getElementPointerFromNumero(numeroDeb);
 
   if ( elPtr->getNomdElement().getElementType() != RFgun ) {
@@ -31 +32 @@
     return false;
   }
-
   ofstream outfile;
   string name = workingDir + inputFileName;
@@ -39 +39 @@
     return false;
   }
-
   outfile << "&INPUT" << endl;
   string fichier = "'" + workingDir + "faisceau.ini" + "'";
@@ -50 +49 @@
   return true;
 }
+
+
+bool  softwareGenerator::execute(string inputFileName, unsigned int numeroDeb,unsigned int numeroFin,string workingDir,string& resul) {
+  ostringstream sortie;
+  bool ExecuteStatus = true;
+  resul.clear();
+  sortie << " EXECUTION DE GENERATOR DE l'ELEMENT " << numeroDeb << " A L'ELEMENT " << numeroFin << endl;
+
+  char buf[132];
+  string generatorJob = workingDir + "generator";
+  generatorJob += string("   ");
+  generatorJob += workingDir + inputFileName;
+
+
+  ofstream outfile;
+  string nameOut = workingDir + "generator.output";
+  string resultOfRun;
+  bool success = launchJob(generatorJob, resultOfRun);
+  sortie << resultOfRun << endl;
+  if ( !success) {
+    sortie << " launching of generator failed " << endl;
+    ExecuteStatus = false;
+  }
+  else {
+    cout << " execution generator MARCHE " << endl;
+    sortie << resultOfRun;
+    outfile.open(nameOut.c_str(), ios::out);
+    if (!outfile) 
+      {
+        sortie << " error first opening transport output stream " << nameOut << endl;
+        ExecuteStatus = false;
+      }
+    else {
+      // on copie la sortie dans un fichier 'generator.out'
+      outfile << resultOfRun << endl;
+      outfile.close();
+    }
+  }
+  resul =  sortie.str();
+  return ExecuteStatus;
+}
+
+bool  softwareGenerator::buildBeamAfterElements(unsigned int numeroDeb,unsigned int numeroFin, vector<particleBeam>& beams, string workingDir) {
+
+  bool result = true;
+  if ( numeroDeb != numeroFin ) return false;
+
+  beams.push_back(particleBeam());
+  vector<double> centroid;
+  bareParticle refPart;
+  vector<bareParticle> particles;
+  vector<bareParticle> particlesPassives; // on ne fait rien de ces particules pour l'instant
+  if (beamFromGenerator(string("faisceau.ini"),workingDir, centroid, refPart,particles, particlesPassives )) {
+    beams.back().setWithParticles(centroid, refPart,particles); 
+  } else {
+    result = false;
+  }
+  return result;
+}
+
+bool softwareGenerator::beamFromGenerator(string beamFileName, string workingDir, vector<double>& centroid, bareParticle& refPart,vector<bareParticle>& particles, vector<bareParticle>& passiveParticles ) {
+  unsigned  k;
+  FILE* filefais;
+  string nomfilefais = workingDir + beamFileName;
+  cout << " nom fichier sortie generator : " << nomfilefais << endl;
+  filefais = fopen(nomfilefais.c_str(), "r");
+  if ( filefais == (FILE*)0 ) {
+    cerr << " beamFromGenerator() erreur a l'ouverture du fichier" << nomfilefais  << endl;; 
+    return false;
+  }
+  else cout << " beamFromGenerator() : ouverture du fichier " << nomfilefais << endl; 
+
+  generatorParticle partic;
+  std::vector<generatorParticle> faisceau;
+  int nbProbPart =0;
+  double timeRef = 0.0;
+  double betagammaZRef = 0.0;
+  // lecture part. de reference 
+  if ( partic.readFromGeneratorFile(filefais) > 0 )
+    {
+         if ( partic.index != 1 ) {
+        cout << " ERROR softwareGenerator::beamFromGenerator : particles are not electrons, we have to reconsider this method " << endl;
+        return false;
+      }
+   
+      if ( partic.flag != -1 ) {
+        cout << " ATTENTION softwareGenerator::beamFromGenerator : flag different de -1 " << endl;
+      }
+      if (fabs(partic.xx) > EPSILON || fabs(partic.yy) > EPSILON || fabs(partic.px) > EPSILON  || fabs(partic.py) > EPSILON) {
+        printf(" ATTENTION softwareGenerator::beamFromGenerator : part. reference douteuse  \n");
+        partic.imprim();
+      }
+      timeRef = partic.clock;
+      TRIDVECTOR  posRef(100.*partic.xx,100.*partic.yy,100.*partic.zz); // en cm
+      betagammaZRef = partic.pz/ERESTeV;
+      // l'impulsion donnee par generator est en eV/c
+      TRIDVECTOR betagammaRef(partic.px/ERESTeV , partic.py/ERESTeV, betagammaZRef);
+      refPart = bareParticle(posRef, betagammaRef);
+
+      // seule la part. de reference a un pz absolu (les autres pat. on un pz relatif a la ref)
+      // je mets ici a zero, pour homogeneiser
+      partic.pz = 0.0;
+      faisceau.push_back(partic);
+    }
+
+  while( partic.readFromGeneratorFile(filefais) > 0 ) {
+    faisceau.push_back(partic);
+    if ( partic.flag != -1 ) nbProbPart++;
+  }
+
+  if ( faisceau.size() == 0) 
+    {
+      cerr << " softwareGenerator::beamFromGenerator echec lecture " << endl;
+      return false;
+    }
+  // pour l'instant on choisit un centroid nul;
+  centroid.clear();
+  centroid = vector<double>(6,0.0);
+
+  particles.clear();
+  passiveParticles.clear();
+  //  particles.resize(faisceau.size() - nbProbPart, bareParticle());
+  //  passiveParticles.resize(nbProbPart, bareParticle());
+  double x,y;
+  double deltaz;
+  double deltat;
+  TRIDVECTOR  pos;
+  TRIDVECTOR betagamma;
+  double pxPart;
+  double pyPart;
+  double pzPartRel;
+  double deltaPzPart;
+
+  for ( k=0; k < faisceau.size(); k++) {
+
+    pxPart = faisceau.at(k).px;
+    pyPart = faisceau.at(k).py;
+    pzPartRel = faisceau.at(k).pz;
+    deltaPzPart = faisceau.at(k).pz;
+    x=faisceau.at(k).xx;
+    y=faisceau.at(k).yy;
+
+    // tout ce qui suit sera a clarifier
+    double betaGammax = pxPart/ERESTeV;
+    double betaGammay = pyPart/ERESTeV;
+    double betaGammaz = betagammaZRef + pzPartRel/ERESTeV;
+    betagamma.setComponents(betaGammax, betaGammay, betaGammaz);
+    double gamma = sqrt(1.0 + betagamma.norm2());
+
+    // decalage temporel par rapport a la reference  
+    deltat = faisceau.at(k).clock - timeRef; // nanoseondes
+    double ds = CLIGHT_m_per_ns * deltat /gamma;
+    x += betaGammax * ds;  // en metres
+    y += betaGammay * ds;
+
+    // ici on neglige la difference entre gamma de la part. et gamma de la ref.
+    deltaz = (pzPartRel/ERESTeV) * CLIGHT_m_per_ns * deltat; // en metres
+
+    pos.setComponents(100.*x,100.*y,100.*deltaz);  // en cm
+    if ( faisceau.at(k).flag == -1 ) {
+      particles.push_back(bareParticle(pos,betagamma));
+    } else {
+      passiveParticles.push_back(bareParticle(pos,betagamma));
+    }
+  }
+
+
+  return true;
+}