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

Timestamp:

Jun 21, 2013, 10:16:56 AM (11 years ago)

Author:

lemeur

Message:

definition des compatibilites des elements dans les software

File:

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

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

@@ -15 +15 @@
 {
   nameOfSoftware_ = nomDeLogiciel("generator");
+  registerElement(nomdElements::beam,TBoolOk);
 }
 
@@ -26 +27 @@
   elPtr = dataManager_->getElementPointerFromNumero(numeroDeb_);
 
-  if ( elPtr->getNomdElement().getElementType() != RFgun ) {
+  if ( elPtr->getNomdElement().getElementType() != nomdElements::RFgun ) {
     dataManager_->consoleMessage(" softwareGenerator::createInputFile : the element must be rfgun " );
     cerr << " softwareGenerator::createInputFile : the element must be rfgun" << endl;
@@ -45 +46 @@
   outfile << " Add=FALSE, N_add=0" << endl;
   outfile << " Species='electrons'" << endl;
-  outfile << elPtr->generatorOutputFlow() << endl;
+  //  outfile << elPtr->generatorOutputFlow() << endl;
+
+  outfile << elementsData(elPtr->parametersToSoftware()) << endl;
+
   outfile << "/" << endl;
   outfile.close();
@@ -231 +235 @@
   return true;
 }
+
+
+string softwareGenerator::elementsData(const vector< pair<string, vector<string> > >& donnees) const {
+  unsigned k;
+  cout << " PASSAGE softwareGenerator::elementsData " << endl;
+  if ( donnees.at(0).first != "labelsGenericSpecific" ) {
+    cout << " softwareGenerator::elementsData ERROR : element badly defined " << endl;
+    return string();
+  }
+  string genericName = donnees.at(0).second.at(0);
+  if ( genericName == "rfgun" ) return rfgunData(donnees);
+  return string();
+}
+
+string softwareGenerator::rfgunData(const vector< pair<string, vector<string> > >& donnees) const {
+
+  cout << " PASSAGE softwareGenerator::rfgunData " << endl;
+  string nmacrop = "0";
+  string sigma_t = "0.0";
+  string sigma_r = "0.0";
+  string emit_x = "0.0";
+  string emit_y = "0.0";
+  string E_cin = "0.0";
+  string sigma_E = "0.0";
+  string phaseStep = "0.0";
+  string totalCharge = "0.0";
+
+  unsigned k;
+  for ( k=1; k < donnees.size(); k++) {
+    if ( donnees.at(k).first == "nbMacroparticles" ) { 
+      nmacrop = donnees.at(k).second.at(0);
+    } else if ( donnees.at(k).first == "sigmasTR" ) {
+      sigma_t = donnees.at(k).second.at(0);
+      sigma_r = donnees.at(k).second.at(1);
+    } else if ( donnees.at(k).first == "emittancesXY" ) {
+      emit_x = donnees.at(k).second.at(0);
+      emit_y = donnees.at(k).second.at(1);
+    } else if ( donnees.at(k).first == "kineticE" ) {
+      E_cin = donnees.at(k).second.at(0);
+      sigma_E = donnees.at(k).second.at(1);
+    } else if ( donnees.at(k).first == "phaseStep" ) {
+      phaseStep = donnees.at(k).second.at(0);
+    } else if ( donnees.at(k).first == "totalCharge" ) {
+      totalCharge = donnees.at(k).second.at(0);
+    }
+  }
+  ostringstream sortie;
+
+
+    sortie << "Ipart=" << nmacrop << endl;
+
+    sortie << "Probe=.True." << endl;
+    sortie << "Noise_reduc=.T" << endl;
+    sortie << "Cathode=.T." << endl;
+    sortie << "Q_total=" << totalCharge << endl;
+    sortie << "Ref_zpos=0.0" << endl;
+    sortie << "Ref_clock=0.0" << endl;
+    sortie << "Ref_Ekin=" << E_cin << endl; // tjs en MeV
+    sortie << "Dist_z='g'" << endl;
+    sortie << "sig_clock=" << 1.0e-3*atof(sigma_t.c_str()) << endl; // passage de ps en ns
+    sortie <<  "Dist_pz='g', sig_Ekin=" << 1000.*atof(sigma_E.c_str()) << ",  emit_z=0.0 ,    cor_Ekin=0.0 " << endl;  // passage en keV
+    sortie << "Dist_x='gauss',  sig_x=" <<  10.*atof(sigma_r.c_str()) << endl;    // passage en mm
+    sortie << "Dist_px='g', Nemit_x=" << emit_x << ",   cor_px=0.0" << endl;
+    sortie << "Dist_y='gauss',  sig_y=" <<  10.*atof(sigma_r.c_str()) << endl;    // passage en mm
+    sortie << "Dist_py='g', Nemit_y=" << emit_y << ",   cor_py=0.0" << endl;
+ 
+  return sortie.str();
+}