Changeset 418 in PSPA for Interface_Web/trunk/pspaWT/sources/controler/src/softwareTransport.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/softwareTransport.cc (modified) (5 diffs)

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

@@ -14 +14 @@
 {
   nameOfSoftware_ = nomDeLogiciel("transport");
+  registerElement(nomdElements::beam,TBoolOk);
+  registerElement(nomdElements::drift,TBoolOk);
+  registerElement(nomdElements::bend,TBoolOk);
+  registerElement(nomdElements::soleno,TBoolOk);
+  registerElement(nomdElements::fit,TBoolOk);
+  registerElement(nomdElements::snapshot,TBoolIgnore);
 }
 
@@ -52 +58 @@
 
   // Dans le cas où le premier element n'est pas beam
-  if ( premierEl->getNomdElement().getElementType() != beam ) {
+  if ( premierEl->getNomdElement().getElementType() != nomdElements::beam ) {
     // s'il n'y a pas faisceau courant, erreur
     if ( beamBefore == NULL ) {
@@ -85 +91 @@
     {
       elPtr = dataManager_->getElementPointerFromNumero(k);
-      cout << "debug:: element [" << k << "] " << elPtr->getNomdElement().getElementName() << endl;
-      outfile << elPtr->transportOutputFlow();
+      cout << "debug:: element [" << k << "] " << elPtr->getNomdElement().getExpandedName() << endl;
+      //      outfile << elPtr->transportOutputFlow();
+      outfile << elementsData(elPtr->parametersToSoftware());
     }
 
@@ -150 +157 @@
 
 
-      if ( ptr->is_accepted_by_software(nameOfSoftware_) != TBoolOk ) {
+      // if ( ptr->is_accepted_by_software(nameOfSoftware_) != TBoolOk ) {
+      if ( doAcceptElement(ptr->getNomdElement().getElementType())  != TBoolOk ) {
+
         // si l'element doit etre ignore, on renvoie sur le diag precedent
         dataManager_->updateCurrentDiagnostic(false);
@@ -241 +250 @@
   return moments.readFromTransportOutput(fichierStream);
 }
+
+string softwareTransport::elementsData(const vector< pair<string, vector<string> > >& donnees) const {
+  unsigned k;
+  cout << " PASSAGE softwareTransport::elementsData " << endl;
+  if ( donnees.at(0).first != "labelsGenericSpecific" ) {
+    cout << " softwareTransport::elementsData ERROR : element badly defined " << endl;
+    return string();
+  }
+  string genericName = donnees.at(0).second.at(0);
+  if ( genericName == "beam" ) return beamData(donnees);
+  if ( genericName == "drift" ) return driftData(donnees);
+  if ( genericName == "solnd" ) return solenoData(donnees);
+  if ( genericName == "bend" ) return bendData(donnees);
+  if ( genericName == "fit" ) return fitData(donnees);
+  return string();
+}
+
+string softwareTransport::beamData(const vector< pair<string, vector<string> > >& donnees) const {
+
+  cout << " PASSAGE softwareTransport::beamData " << endl;
+  double x = 0.0;
+  double xp = 0.0;
+  double y = 0.0;
+  double yp = 0.0;
+  double dl = 0.0;
+  double del = 0.0;
+  double p0 = 0.0;
+  string specificName = "";
+
+  unsigned k;
+  for ( k=0; k < donnees.size(); k++) {
+    if ( donnees.at(k).first == "transverseCoordinates" ) { 
+      x = atof( donnees.at(k).second.at(0).c_str() );
+      xp = atof( donnees.at(k).second.at(1).c_str() );
+      y = atof( donnees.at(k).second.at(2).c_str() );
+      yp = atof( donnees.at(k).second.at(3).c_str() );
+    } else if ( donnees.at(k).first == "longitudinalCoordinates" ) {
+      dl = atof( donnees.at(k).second.at(0).c_str() );
+      del = atof( donnees.at(k).second.at(1).c_str() );
+    } else if ( donnees.at(k).first == "momentum" ) {
+      p0 = atof( donnees.at(k).second.at(0).c_str() );
+    } else if ( donnees.at(k).first == "labelsGenericSpecific" ) {
+      specificName = donnees.at(k).second.at(1);
+    }
+  }
+  ostringstream sortie;
+    beam2Moments moments(x, xp, y, yp, dl, del);
+    sortie << specificName << ":" << moments.writeToTransportInput(p0) << endl;
+    return sortie.str(); 
+  return sortie.str();
+}
+
+string softwareTransport::driftData(const vector< pair<string, vector<string> > >& donnees) const {
+  cout << " PASSAGE softwareTransport::driftData " << endl;
+
+  string lenght = "0.0";
+  string specificName = "";
+
+  unsigned k;
+  for ( k=0; k < donnees.size(); k++) {
+    if ( donnees.at(k).first == "lenghtAperture" ) { 
+      lenght = donnees.at(k).second.at(0);
+    } else if ( donnees.at(k).first == "labelsGenericSpecific" ) {
+      specificName = donnees.at(k).second.at(1);
+    }
+  }
+  ostringstream sortie;
+    // longueur en m
+  sortie << specificName << ":" << "  DRIFT, L=" << 0.01*atof(lenght.c_str()) << ";" << endl;
+  return sortie.str();
+}
+
+
+string softwareTransport::solenoData(const vector< pair<string, vector<string> > >& donnees) const {
+  cout << " softwareTransport::solenoData " << endl;
+  string lenght = "0.0";
+  string B0 = "0.0";
+  string varyB0 = "0";
+  string specificName = "";
+
+
+  unsigned k;
+  for ( k=0; k < donnees.size(); k++) {
+    if ( donnees.at(k).first == "lenghtAperture" ) { 
+      lenght = donnees.at(k).second.at(0);
+    } else if ( donnees.at(k).first == "field" ) {
+      B0 = donnees.at(k).second.at(0);
+    } else if ( donnees.at(k).first == "isBvariable" ) {
+      varyB0 = donnees.at(k).second.at(0);
+    } else if ( donnees.at(k).first == "labelsGenericSpecific" ) {
+      specificName = donnees.at(k).second.at(1);
+    }
+  }
+    ostringstream sortie;
+    // on passe la longueur en metres 
+    sortie << specificName << ":" << "  SOLENOID, L=" << 0.01*atof (lenght.c_str()) << ", B=" << B0 << ";" << endl;
+    if ( atoi(varyB0.c_str()) ) {
+      sortie <<  "  VARY, " << specificName+"[B]"  << ";" << endl;
+    }
+    return sortie.str();
+
+}
+
+
+string softwareTransport::bendData(const vector< pair<string, vector<string> > >& donnees) const {
+
+  cout << " softwareTransport::bendData " << endl;
+
+  string lenght = "0.0";
+  string radius = "0.0";
+  string angleDeg = "0.0";
+  string beta1 = "0.0";
+  string beta2 = "0.0";
+
+
+
+  string specificName = "";
+
+  unsigned k;
+  for ( k=0; k < donnees.size(); k++) {
+    if ( donnees.at(k).first == "lenghtAperture" ) { 
+      lenght = donnees.at(k).second.at(0);
+    } else if ( donnees.at(k).first == "angleDegre" ) {
+      angleDeg = donnees.at(k).second.at(0);
+    } else if ( donnees.at(k).first == "radius" ) {
+      radius = donnees.at(k).second.at(0);
+     } else if ( donnees.at(k).first == "rotatedFaces" ) {
+      beta1 = donnees.at(k).second.at(0);
+      beta2 = donnees.at(k).second.at(1);
+    } else if ( donnees.at(k).first == "labelsGenericSpecific" ) {
+      specificName = donnees.at(k).second.at(1);
+    }   
+  }
+
+    ostringstream sortie;
+    sortie << specificName << ":" << "  SBEND, RADIUS=" << radius << ", ANGLE=" << angleDeg << ", E1=" << beta1 << ", E2=" << beta2 << ";" << endl;
+    return sortie.str();
+
+}
+
+
+string softwareTransport::fitData(const vector< pair<string, vector<string> > >& donnees) const {
+  cout << " softwareTransport::fitData " << endl;
+  string variableName = "";
+  string variableValue = "0.0";
+  string tolerance = "0.0";
+  string specificName = "";
+
+  unsigned k;
+  for ( k=0; k < donnees.size(); k++) {
+    if ( donnees.at(k).first == "variable") { 
+      variableName = donnees.at(k).second.at(0);
+      variableValue = donnees.at(k).second.at(1);
+      tolerance = donnees.at(k).second.at(2);
+    } else if ( donnees.at(k).first == "labelsGenericSpecific" ) {
+      specificName = donnees.at(k).second.at(1);
+    }   
+  }
+    ostringstream sortie;
+    sortie << specificName << ":" << "  FIT," << variableName << " = " << variableValue << ", TOLER = " << tolerance << " ;" << endl;
+    if ( variableName == "R11" || variableName == "R22" || variableName == "R12" || variableName == "R21" ) {
+        sortie << " PRINT, TRANS;" << endl;
+    }
+    return sortie.str();
+}