Changeset 294 for ELYSE/HEAD/source/SteppingAction.cxx

Timestamp:

Feb 23, 2007, 2:22:42 PM (17 years ago)

Author:

campagne

Message:

Unuseful Sensitive Detector for the moment (JEC)

File:

• ELYSE/HEAD/source/SteppingAction.cxx (modified) (4 diffs)

ELYSE/HEAD/source/SteppingAction.cxx

@@ -20 +20 @@
 #include "ELYSE/TrackingAction.hh"
 #include "ELYSE/Trajectory.hh"
-#include "ELYSE/TrappingVolume.hh"
+//#include "ELYSE/TrappingVolume.hh"
 #include "ELYSE/TrackInformation.hh"
 
@@ -40 +40 @@
   G4VPhysicalVolume* thePostPV = thePostPoint->GetPhysicalVolume();
 
+  
+  //JEC START debug Optical surface properties
+//   G4StepPoint* pPreStepPoint  = theStep->GetPreStepPoint();
+//   G4cout << "(JEC): UserSteppingAction PrePoint Volume:" <<  pPreStepPoint->GetPhysicalVolume()->GetName()
+//       << " PostPoint Volume: " << thePostPV->GetName()
+//       << G4endl;
+//   G4LogicalSurface* Surface = 
+//     G4LogicalSkinSurface::GetSurface(pPreStepPoint->GetPhysicalVolume()->GetLogicalVolume());
+//   G4OpticalSurface* OpticalSurface = 0;
+//   if (Surface) OpticalSurface = (G4OpticalSurface*) Surface->GetSurfaceProperty();
+//   if (OpticalSurface) { 
+//     G4cout << "(JEC) UserSteppingAction: Optical Surface Found" << G4endl;
+//     OpticalSurface->DumpInfo();    
+//     G4MaterialPropertiesTable* aMaterialPropertiesTable;
+//     G4MaterialPropertyVector* Rindex;
+//     aMaterialPropertiesTable = OpticalSurface->GetMaterialPropertiesTable();
+    
+//     if (aMaterialPropertiesTable) {
+//      Rindex = aMaterialPropertiesTable->GetProperty("RINDEX");
+//      if (Rindex) {
+//        G4cout << "(JEC) Found RINDEX properties: " << G4endl;
+//        Rindex->DumpVector();
+//      } else {
+//        G4cout << "(JEC) NO  RINDEX properties " << G4endl;
+//      }
+//     } else {
+//       G4cout << "(JEC) NO  MATERIAL properties " << G4endl;
+//     }
+//   } else {
+//     G4cout << "(JEC) NO OpticalSurface " << G4endl;
+//   }
+  
   G4OpBoundaryProcessStatus boundaryStatus=Undefined;
   static G4OpBoundaryProcess* boundary=NULL;
@@ -74 +106 @@
     
     //------------Kill photons entering expHall
+
+//     G4cout << "(JEC) UserSteppingAction: Post Volume " << thePostPV->GetName()
+//         << G4endl;
+
     if(thePostPV->GetName()=="expHall") {
       theTrack->SetTrackStatus(fStopAndKill);}
@@ -94 +130 @@
 //                                   Absorption, Detection, NotAtBoundary,
 //                                   SameMaterial, StepTooSmall, NoRINDEX}
-      switch(boundaryStatus){
-      case FresnelRefraction:
-        G4cout << "(JEC) Status of PostPoint: FresnelRefraction" << G4endl; break;
-      case FresnelReflection:
-        G4cout << "(JEC) Status of PostPoint: FresnelReflection" << G4endl; break;
-      case TotalInternalReflection:
-        G4cout << "(JEC) Status of PostPoint: TotalInternalReflection" << G4endl; break;
-       case LambertianReflection:
-        G4cout << "(JEC) Status of PostPoint: LambertianReflection" << G4endl; break;
-      case LobeReflection:
-        G4cout << "(JEC) Status of PostPoint: LobeReflection" << G4endl; break;
-      case SpikeReflection:
-        G4cout << "(JEC) Status of PostPoint: SpikeReflection" << G4endl; break;
-      case Absorption:
-        G4cout << "(JEC) Status of PostPoint: Absorption" << G4endl; break;
-     
-      case Detection: // Note, this assumes that the volume causing detection
-                      // is the interface between Tyvek and Air because it is the only one with
-                      // 100% efficiency
-        {
-          // Triger sensitive detector manually since photon is
-          // absorbed but status was Detection
-          G4SDManager* SDman = G4SDManager::GetSDMpointer();
-          //see DetectorConstruction::Construct for the PATH of the Sensitive Volume
-          TrappingVolume* trapVol = 0;
-          if(SDman) trapVol = (TrappingVolume*)SDman->FindSensitiveDetector("/ELYSE/TrappingVolume");
+//       switch(boundaryStatus){
+//       case FresnelRefraction:
+//              G4cout << "(JEC) Status of PostPoint: FresnelRefraction" << G4endl; break;
+//       case FresnelReflection:
+//      G4cout << "(JEC) Status of PostPoint: FresnelReflection" << G4endl; break;
+//       case TotalInternalReflection:
+//      G4cout << "(JEC) Status of PostPoint: TotalInternalReflection" << G4endl; break;
+//        case LambertianReflection:
+//      G4cout << "(JEC) Status of PostPoint: LambertianReflection" << G4endl; break;
+//       case LobeReflection:
+//      G4cout << "(JEC) Status of PostPoint: LobeReflection" << G4endl; break;
+//       case SpikeReflection:
+//      G4cout << "(JEC) Status of PostPoint: SpikeReflection" << G4endl; break;
+//       case Absorption:
+//      G4cout << "(JEC) Status of PostPoint: Absorption" << G4endl; break;
+//       case Detection: 
+//      G4cout << "(JEC) Status of PostPoint: Detection" << G4endl; break;
+//       case NotAtBoundary:
+//      G4cout << "(JEC) Status of PostPoint: NotAtBoundary" << G4endl; break;
+//       case SameMaterial:
+//      G4cout << "(JEC) Status of PostPoint: SameMaterial" << G4endl; break;
+//       case StepTooSmall:
+//      G4cout << "(JEC) Status of PostPoint: StepTooSmall" << G4endl; break;
+//       case NoRINDEX:
+//      G4cout << "(JEC) Status of PostPoint: NoRINDEX" << G4endl; break;
 
-          if(trapVol) { 
-            trapVol->ProcessHits_constStep(theStep,0);
+//       default:
+//      break;
 
-            //------------Save the hit into the track info
-
-            TrackInformation* aInfo = (TrackInformation*)theTrack->GetUserInformation();
-            G4cout << "(JEC) SteppingAction::opticalphoton detection = " 
-                   << " " << theTrack << " " << aInfo << G4endl; 
-            if(aInfo)aInfo->AddTrackStatusFlag(hitProduced);
-
-          } else {
-            G4cout << "(JEC) UserSteppingAction: FATAL: NO Sensitive Det. FOUND" << G4endl;
-            exit(0);
-          }
-          break;
-        }//case Detection
-
-      case NotAtBoundary:
-        G4cout << "(JEC) Status of PostPoint: NotAtBoundary" << G4endl; break;
-      case SameMaterial:
-        G4cout << "(JEC) Status of PostPoint: SameMaterial" << G4endl; break;
-      case StepTooSmall:
-        G4cout << "(JEC) Status of PostPoint: StepTooSmall" << G4endl; break;
-      case NoRINDEX:
-        G4cout << "(JEC) Status of PostPoint: NoRINDEX" << G4endl; break;
-
-      default:
-        break;
-
-      }//eo switch on boundary values
+//       }//eo switch on boundary values
 
     }//test on GeomBoundary