Changeset 161 for MEMPHYS

Timestamp:

Jun 16, 2006, 10:03:40 AM (18 years ago)

Author:

barrand

Message:

 

File:

• MEMPHYS/HEAD/applications/MEMPHYS_analysis.cxx (modified) (4 diffs)

MEMPHYS/HEAD/applications/MEMPHYS_analysis.cxx

@@ -1 @@
-/*
- * AIDA analysis program over MEMPHYS.root file.
- */
+// AIDA analysis program over MEMPHYS.root file.
 
 // AIDA :
@@ -9 +7 @@
 #include <iostream>
 
+//////////////////////////////////////////////////////////////////////////////
+bool direction(
+ AIDA::ITuple& aParent
+,double& aX
+,double& aY
+,double& aZ
+) 
+//////////////////////////////////////////////////////////////////////////////
+//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!//
+{
+  AIDA::ITupleEntry* entry = (AIDA::ITupleEntry*)aParent.getObject(3);
+  if(!entry) return false;
+
+  AIDA::ITuple* tuple = dynamic_cast<AIDA::ITuple*>(entry);
+  if(!tuple) return false;
+
+  tuple->start();
+  if(!tuple->next()) return false;
+
+  aX = tuple->getDouble(0);
+  aY = tuple->getDouble(1);
+  aZ = tuple->getDouble(2);
+
+  return true;
+}
+//////////////////////////////////////////////////////////////////////////////
+bool dump_tracks(
+ AIDA::ITuple& aParent
+) 
+//////////////////////////////////////////////////////////////////////////////
+//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!//
+{
+  AIDA::ITupleEntry* entry = (AIDA::ITupleEntry*)aParent.getObject(8);
+  if(!entry) return false;
+
+  AIDA::ITuple* tuple = dynamic_cast<AIDA::ITuple*>(entry);
+  if(!tuple) return false;
+
+  //if(nPart!=tracks->rows()) {
+  //  std::cout << "read: nPart / tracks mismatch " 
+  //              << nPart << " " << tracks->rows() << std::endl;
+  //  return false;
+  //}
+
+  tuple->start();
+  while(tuple->next()) {
+
+    double dx,dy,dz;
+    if(!direction(*tuple,dx,dy,dz)) return false;
+
+    double ETot = tuple->getDouble(6);
+
+    std::cout << " ETot " << ETot 
+              << " dx " << dx
+              << " dy " << dy
+              << " dz " << dz
+              << std::endl;
+  }
+  return true;
+}
 //////////////////////////////////////////////////////////////////////////////
 int main(
@@ -75 +133 @@
 
     int eventId = tuple->getInt(0);
-    int inputEvtId = tuple->getInt(1);
-    int interMode = tuple->getInt(2);
-    int vtxVol = tuple->getInt(3);
+    //int inputEvtId = tuple->getInt(1);
+    //int interMode = tuple->getInt(2);
+    //int vtxVol = tuple->getInt(3);
 
     int nPart = tuple->getInt(5);
-    int leptonIndex = tuple->getInt(6);
-    int protonIndex = tuple->getInt(7);
+    //int leptonIndex = tuple->getInt(6);
+    //int protonIndex = tuple->getInt(7);
 
     int nHits = tuple->getInt(9);
@@ -98 +156 @@
               << std::endl;
 
+    if(!dump_tracks(*tuple)) return false;
+
     irow++;
   }