Changeset 162 for MEMPHYS


Ignore:
Timestamp:
Jun 16, 2006, 10:16:34 AM (18 years ago)
Author:
barrand
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed

  • MEMPHYS/HEAD/applications/MEMPHYS_analysis.cxx

    r161 r162  
    88
    99//////////////////////////////////////////////////////////////////////////////
    10 bool direction(
     10bool getXYZ(
    1111 AIDA::ITuple& aParent
     12,int aColumn
    1213,double& aX
    1314,double& aY
    … …  
    1718//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!//
    1819{
    19   AIDA::ITupleEntry* entry = (AIDA::ITupleEntry*)aParent.getObject(3);
     20  AIDA::ITupleEntry* entry = (AIDA::ITupleEntry*)aParent.getObject(aColumn);
    2021  if(!entry) return false;
    2122
    … …  
    3233  return true;
    3334}
     35//////////////////////////////////////////////////////////////////////////////
     36//////////////////////////////////////////////////////////////////////////////
    3437//////////////////////////////////////////////////////////////////////////////
    3538bool dump_tracks(
    … …  
    5255
    5356  tuple->start();
     57  int irow = 0;
    5458  while(tuple->next()) {
    5559
     60    int pId = tuple->getInt(0);
     61    int parent = tuple->getInt(1);
     62    int timeStart = tuple->getInt(2);
     63
    5664    double dx,dy,dz;
    57     if(!direction(*tuple,dx,dy,dz)) return false;
     65    if(!getXYZ(*tuple,3,dx,dy,dz)) return false;
    5866
     67    double mass = tuple->getDouble(4);
     68    double pTot = tuple->getDouble(5);
    5969    double ETot = tuple->getDouble(6);
    6070
    61     std::cout << " ETot " << ETot
    62               << " dx " << dx
    63               << " dy " << dy
    64               << " dz " << dz
    65               << std::endl;
     71    double px,py,pz;
     72    if(!getXYZ(*tuple,7,dx,dy,dz)) return false;
     73    double start_x,start_y,start_z;
     74    if(!getXYZ(*tuple,8,dx,dy,dz)) return false;
     75    double stop_x,stop_y,stop_z;
     76    if(!getXYZ(*tuple,9,dx,dy,dz)) return false;
     77
     78    int startVol = tuple->getDouble(10);
     79    int stopVol = tuple->getDouble(11);
     80
     81    std::cout << "----> Tk{"<<irow<<"}: "
     82              << " pId " << pId
     83              << " parent " << parent
     84                << " creation time " << timeStart
     85                << " Volumes " << startVol << " " << stopVol << "\n"
     86                << " Start Pos (" << start_x << "," << start_y << "," << start_z << ")\n"
     87                << " Stop Pos (" << stop_x << "," << stop_y << "," << stop_z << ")\n"
     88                << " dx,dy,dz " << dx << " " << dy << " " << dz << "\n"
     89                << " m " << mass
     90                << " ETot " << ETot
     91                << " pTot " << pTot
     92                << " px,py,pz " << px << " " << py << " " << pz << "\n"
     93                << std::endl;
     94
     95    irow++;
    6696  }
    6797  return true;
    … …  
    123153
    124154  //int nentries = tuple->rows();
    125   int nentries = 100000;
     155  //int nentries = 100000;
     156  int nentries = 40;
    126157  std::cout << header
    127158            << "traitements de " << nentries << " entrees"
Note: See TracChangeset for help on using the changeset viewer.