Changeset 261 in Idarraga for allpix/src/AllPixRun.cc

Timestamp:

Nov 21, 2011, 7:07:58 PM (13 years ago)

Author:

idarraga

Message:

 

File:

• allpix/src/AllPixRun.cc (modified) (3 diffs)

allpix/src/AllPixRun.cc

@@ -9 +9 @@
 #include "G4SDManager.hh"
 #include "G4DigiManager.hh"
+#include "G4Trajectory.hh"
 
 // digits, frames
@@ -29 +30 @@
  * This constructor is called once per run
  */
-AllPixRun::AllPixRun(TString ofp, TString dataset, TString tempDir){
+AllPixRun::AllPixRun(AllPixDetectorConstruction * det, TString ofp, TString dataset, TString tempDir){
+
+        m_detectorPtr = det;
 
         // Call for an instance to write.  Need to know how
@@ -387 +390 @@
         *m_lciobridge_dut_f << "R " << GetRunID() << endl;
 
+        // check event for information about the track
+        //G4TrajectoryContainer * tC = evt->GetTrajectoryContainer();
+        //G4Trajectory * tr;
+        G4TrajectoryContainer * tC = evt->GetTrajectoryContainer();
+        G4cout << "Trajectories : " << tC->GetVector()->size() << G4endl;
+        TrajectoryVector * tV = tC->GetVector();
+        TrajectoryVector::iterator tI = tV->begin();
+        TrajectoryVector::iterator tIE = tV->end();
+
+        G4Trajectory * tr;
+        G4TrajectoryPoint* tP;
+        for ( ; tI != tIE ; tI++ ) {
+                tr = static_cast<G4Trajectory*>(*tI);
+                //G4cout << "    " << tr->GetPDGEncoding() << G4endl;
+                //G4cout << "       " << tr->GetPointEntries() << G4endl;
+                tP = static_cast<G4TrajectoryPoint*>((*tI)->GetPoint(0));
+
+                // Check which of the points is at the interior of the sensor
+                // Use the G4VSolid::Inside(const G4ThreeVector& p)
+                // Find first the Silicon G4VSolids (sensors).
+
+
+                //G4VSolid * sol = m_detectorPtr->GetVSolidDetector(100);
+                //G4cout << "       " << tP->GetPosition().z()/mm << G4endl;
+/*
+                if(sol->Inside(tP->GetPosition()) == kInside){
+                        G4cout << "       inside !!" << G4endl;
+                } else {
+                        G4cout << "       outside !!" << G4endl;
+                }
+*/
+        }
+
         // Ok, I have to match the right detector m_frames with the digit collection
         for (G4int i = 0 ; i < nDC ; i++) {