Skip to content
Snippets Groups Projects
Commit a7b45302 authored by fangwx@ihep.ac.cn's avatar fangwx@ihep.ac.cn
Browse files

change G4ThreeVector to DDRec Vector3D

parent d1eea3bf
No related branches found
No related tags found
No related merge requests found
gaudi_subdir(DCHDigi v0r0)
find_package(Geant4 REQUIRED ui_all vis_all)
include(${Geant4_USE_FILE})
find_package(DD4hep COMPONENTS DDG4 REQUIRED)
find_package(EDM4HEP REQUIRED )
message("EDM4HEP_INCLUDE_DIRS: ${EDM4HEP_INCLUDE_DIR}")
......
......@@ -8,7 +8,7 @@
#include "DD4hep/Detector.h"
#include <DD4hep/Objects.h>
#include "G4ThreeVector.hh"
#include "DDRec/Vector3D.h"
#include <array>
......@@ -90,12 +90,12 @@ StatusCode DCHDigiAlg::execute()
float min_distance = 999 ;
for(unsigned int i=0; i< simhit_size; i++)
{
G4ThreeVector west(0,0,0);
G4ThreeVector east(0,0,0);
G4ThreeVector pos(iter->second.at(i).getPosition()[0], iter->second.at(i).getPosition()[1], iter->second.at(i).getPosition()[2]);
G4ThreeVector numerator = (east-west).cross(west-pos) ;
G4ThreeVector denominator = (east-west) ;
float distance = numerator.mag()/denominator.mag() ;
dd4hep::rec::Vector3D west(0,0,0);
dd4hep::rec::Vector3D east(0,0,0);
dd4hep::rec::Vector3D pos(iter->second.at(i).getPosition()[0], iter->second.at(i).getPosition()[1], iter->second.at(i).getPosition()[2]);
dd4hep::rec::Vector3D numerator = (east-west).cross(west-pos) ;
dd4hep::rec::Vector3D denominator = (east-west) ;
float distance = numerator.r()/denominator.r() ;
std::cout<<"distance="<<distance<<std::endl;
if(distance < min_distance){
min_distance = distance;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment