Ignore:
Timestamp:
Mar 1, 2013, 3:38:12 PM (11 years ago)
Author:
lemeur
Message:

complement dessin phase space + rationalisation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • Interface_Web/trunk/pspaWT/sources/controler/src/softwareParmela.cc

    r353 r354  
    163163bool softwareParmela::beamFromParmela(string workingDir,unsigned numeroParmel, double referencefrequency, vector<double>& centroid, bareParticle& refPart,vector<bareParticle>& particles )
    164164{
     165
    165166  string nomfilefais = workingDir + "parmdesz";
    166167  cout << " nom fichier desz : " << nomfilefais << endl;
     
    183184
    184185      if ( partic.ne == (int)numeroParmel ) {
    185         faisceau.push_back(partic);
    186186        if ( partic.np == 1 ) {
    187187          // en principe on est sur la particule de reference
     
    191191          }
    192192          phaseRef = partic.phi;
    193           TRIDVECTOR  posRef(partic.xx,partic.yy,0.0);
     193
     194          // le 'z' est 'absolu' (le long de la trajectoire)
     195          TRIDVECTOR  posRef(partic.xx,partic.yy,partic.z);
    194196          TRIDVECTOR betagammaRef(partic.xxp*partic.begamz, partic.yyp*partic.begamz, partic.begamz);
    195197          refPart = bareParticle(posRef, betagammaRef);
     
    200202            return false;
    201203          }
    202         }
     204        }
     205        faisceau.push_back(partic);
    203206      }
    204207    } //while
Note: See TracChangeset for help on using the changeset viewer.