source: PSPA/Interface_Web/trunk/pspaWT/sources/controler/src/elementBend.cc @ 418

#include "elementBend.h"
#include "mixedTools.h"
#include "mathematicalConstants.h"
#include "PhysicalConstants.h"

elementBend::elementBend() :  abstractElement()
{
  setDefaultValues();
  setDefaults();
  elementName_ = nomdElements("bend");
  nbParam_ = elementName_.getElementNbParameters();
  parametersString_ = new string[nbParam_+1];
    // registerAcceptableSoftware(nomDeLogiciel::transport, TBoolOk);
    // registerAcceptableSoftware(nomDeLogiciel::parmela, TBoolOk);
}

  void elementBend::setDefaultValues()
  {
    defaultSpecificName_ = "bend";
    angleDegDef_ = 0.0;
    radiusDef_ = 1.0;
    beta1Def_ = 0.0;
    beta2Def_ = 0.0;
    momentumDef_ = 0.0;
  }

  void elementBend::setDefaults()
  {
    specificName_ = defaultSpecificName_;
    angleDeg_ = angleDegDef_;
    radius_ = radiusDef_;
    beta1_ = beta1Def_;
    beta2_ = beta2Def_;
    momentum_ = momentumDef_;
  }

 string* elementBend::getParametersString() const 
 {
   int compteur = -1;
   parametersString_[++compteur] = mixedTools::intToString(nbParam_);
   parametersString_[++compteur] = specificName_;
   parametersString_[++compteur] = mixedTools::doubleToString( angleDeg_);
   parametersString_[++compteur] = mixedTools::doubleToString(radius_);
   parametersString_[++compteur] = mixedTools::doubleToString(momentum_);
   parametersString_[++compteur] = mixedTools::doubleToString(aperture_);
   parametersString_[++compteur] = mixedTools::doubleToString(beta1_);
   parametersString_[++compteur] = mixedTools::doubleToString(beta2_);
   if ( compteur != nbParam_ ) {
     cerr << " elementBend::getParametersString() : ERROR nr pf parameters dosnt match " << endl;
     return NULL;
   }
   return parametersString_;
 }

 double elementBend::getLenghtOfElement() const 
{
 return 1.0e2 * radius_ * angleDeg_ * DEG_TO_RADIANS; 
}


 void  elementBend::setParametersString(string* param)
 {
        if ( param == NULL )
        {
          cerr << "  elementBend::setParametersString parameters empty parameter set";
          return;
        }

   int compteur = -1;
   int nbparam = atoi(param[++compteur].c_str());
    if ( nbparam != nbParam_ ) 
        {
          cerr << "  elementBend::setParametersString parameters do not match for a BEND" << endl;
          return;
        }
    specificName_ = param[++compteur];
    angleDeg_ =  atof(param[++compteur].c_str());
    radius_ = atof(param[++compteur].c_str());
    momentum_ = atof(param[++compteur].c_str());
    aperture_  = atof(param[++compteur].c_str());
    beta1_ = atof(param[++compteur].c_str());
    beta2_ = atof(param[++compteur].c_str());
    // calcul de la longueur : rayon en m ; longueur en cm
    lenghtElem_ = getLenghtOfElement();
    cout << " elementBend::setParametersString calcule longueur = " << lenghtElem_ << endl;
 }

// string elementBend::parmelaOutputFlow() const
//   {
//     ostringstream sortie;
//     // il faut entrer l'energie cinetique
//     double ecin = momentum_/EREST_MeV;
//     ecin = ecin*ecin + 1.;
//     ecin = EREST_MeV*(sqrt(ecin) - 1.0);
//     sortie << "BEND /l=" << lenghtElem_ << "  / aper=" << aperture_ << "  / iout=1 / wr=" << ecin << " /alpha=" << angleDeg_ << " / beta1=" << beta1_;
//     sortie << " / beta2=" << beta2_  << endl;

//     return sortie.str();
//   }


// string elementBend::transportOutputFlow() const
//   {
//     ostringstream sortie;
//     sortie << specificName_ << ":" << "  SBEND, RADIUS=" << radius_ << ", ANGLE=" << angleDeg_ << ", E1=" << beta1_ << ", E2=" << beta2_ << ";" << endl;
//     return sortie.str();
//   }


string elementBend::FileOutputFlow() const
  {
    ostringstream sortie;
    //    sortie << elementName_.getElementType() << endl;
    sortie << elementName_.getGenericLabel() << endl;
    sortie  << specificName_ << endl;
    sortie << angleDeg_ << "  " << radius_ << " " <<  aperture_  << endl;
    sortie << beta1_ << "  " << beta2_ << " " << momentum_ << endl;
    return sortie.str();
  }

vector< pair<string, vector<string> > > elementBend::parametersToSoftware () const{
  vector< pair<string, vector<string> > >  sortie;
  sortie.push_back( pair<string, vector<string> >("labelsGenericSpecific",  vector<string>() ) );
  sortie.back().second.push_back(elementName_.getGenericLabel());
  sortie.back().second.push_back(specificName_);
  sortie.push_back( pair<string, vector<string> >("lenghtAperture",  vector<string>() ) );
  sortie.back().second.push_back(mixedTools::doubleToString( lenghtElem_));
  sortie.back().second.push_back(mixedTools::doubleToString( aperture_));
  sortie.push_back( pair<string, vector<string> >("angleDegre",  vector<string>() ) );
  sortie.back().second.push_back(mixedTools::doubleToString( angleDeg_));
  sortie.push_back( pair<string, vector<string> >("radius",  vector<string>() ) );
  sortie.back().second.push_back(mixedTools::doubleToString( radius_));
  sortie.push_back( pair<string, vector<string> >("rotatedFaces",  vector<string>() ) );
  sortie.back().second.push_back(mixedTools::doubleToString( beta1_));
  sortie.back().second.push_back(mixedTools::doubleToString( beta2_));
  sortie.push_back( pair<string, vector<string> >("momentum",  vector<string>() ) );
  sortie.back().second.push_back(mixedTools::doubleToString( momentum_));
  return sortie;
}

void elementBend::FileInput(ifstream& ifs)
{
  ifs >> specificName_;
  ifs >> angleDeg_ >> radius_ >>  aperture_;
  ifs >> beta1_ >> beta2_ >> momentum_;
    lenghtElem_ = getLenghtOfElement();
    cout << " elementBend::FileInput calcule longueur = " << lenghtElem_ << endl;
}


string elementBend::print() {
    string txt = "";

    txt += specificName_;
    txt += "\nangle (deg.) : ";
    txt += mixedTools::doubleToString(angleDeg_);
    txt += "\nradius (m) : ";
    txt += mixedTools::doubleToString(radius_);
    txt += "\nmomentum (MeV/c) : ";
    txt += mixedTools::doubleToString(momentum_);
    txt += "\nlength (cm) : ";
    txt += mixedTools::doubleToString(lenghtElem_);
    txt += "\naperture (cm) : ";
    txt += mixedTools::doubleToString(aperture_);
    txt += "\npole face rotations (deg.) entrance : ";
    txt += mixedTools::doubleToString(beta1_);
    txt += "\nexit : ";
    txt += mixedTools::doubleToString(beta2_);
    
    return txt;
}