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

space

parent f8da4a07
No related branches found
No related tags found
No related merge requests found
...@@ -84,10 +84,10 @@ static dd4hep::Ref_t create_detector(dd4hep::Detector& theDetector, ...@@ -84,10 +84,10 @@ static dd4hep::Ref_t create_detector(dd4hep::Detector& theDetector,
dd4hep::Tube det_chamber_solid(chamber_radius_min, chamber_radius_max, chamber_half_length); dd4hep::Tube det_chamber_solid(chamber_radius_min, chamber_radius_max, chamber_half_length);
dd4hep::Volume det_chamber_vol(det_name+"_chamber_vol", det_chamber_solid, chamber_mat); dd4hep::Volume det_chamber_vol(det_name+"_chamber_vol", det_chamber_solid, chamber_mat);
if ( x_det.isSensitive() ) { if ( x_det.isSensitive() ) {
det_chamber_vol.setRegion(theDetector,x_det.regionStr()); det_chamber_vol.setRegion(theDetector,x_det.regionStr());
det_chamber_vol.setLimitSet(theDetector,x_det.limitsStr()); det_chamber_vol.setLimitSet(theDetector,x_det.limitsStr());
det_chamber_vol.setSensitiveDetector(sens); det_chamber_vol.setSensitiveDetector(sens);
sd.setType("tracker"); sd.setType("tracker");
} }
// - wall // - wall
...@@ -185,7 +185,7 @@ static dd4hep::Ref_t create_detector(dd4hep::Detector& theDetector, ...@@ -185,7 +185,7 @@ static dd4hep::Ref_t create_detector(dd4hep::Detector& theDetector,
// | | // | |
// | F0 F1 F2 F3| // | F0 F1 F2 F3|
// ----------------------- // -----------------------
// if(layer_id == 0 || layer_id == 1 || layer_id == 2 || layer_id == 99) { if(layer_id == 0 || layer_id == 1 || layer_id == 2 || layer_id == 99) {
for(int icell=0; icell< numWire; icell++) { for(int icell=0; icell< numWire; icell++) {
double wire_phi = (icell+0.5)*layer_Phi + offset; double wire_phi = (icell+0.5)*layer_Phi + offset;
// - signal wire // - signal wire
...@@ -208,7 +208,7 @@ static dd4hep::Ref_t create_detector(dd4hep::Detector& theDetector, ...@@ -208,7 +208,7 @@ static dd4hep::Ref_t create_detector(dd4hep::Detector& theDetector,
} }
} }
// } }
// - place in det // - place in det
// - chamber // - chamber
......
...@@ -149,6 +149,7 @@ StatusCode DCHDigiAlg::execute() ...@@ -149,6 +149,7 @@ StatusCode DCHDigiAlg::execute()
TVector3 numerator = denominator.Cross(Wstart-pos) ; TVector3 numerator = denominator.Cross(Wstart-pos) ;
float tmp_distance = numerator.Mag()/denominator.Mag() ; float tmp_distance = numerator.Mag()/denominator.Mag() ;
//std::cout<<"tmp_distance="<<tmp_distance<<",x="<<pos.x()<<",y="<<pos.y()<<",z="<<pos.z()<<",mom="<<sim_hit_mom<<",pt="<<sim_hit_pt<<std::endl; //std::cout<<"tmp_distance="<<tmp_distance<<",x="<<pos.x()<<",y="<<pos.y()<<",z="<<pos.z()<<",mom="<<sim_hit_mom<<",pt="<<sim_hit_pt<<std::endl;
if(tmp_distance < min_distance){ if(tmp_distance < min_distance){
min_distance = tmp_distance; min_distance = tmp_distance;
pos_x = pos.x(); pos_x = pos.x();
...@@ -169,7 +170,7 @@ StatusCode DCHDigiAlg::execute() ...@@ -169,7 +170,7 @@ StatusCode DCHDigiAlg::execute()
m_n_sim ++ ; m_n_sim ++ ;
} }
} }
trkHit.setTime(min_distance*1e3/m_velocity);//m_velocity is um/ns, drift time in ns trkHit.setTime(min_distance*1e3/m_velocity);//m_velocity is um/ns, drift time in ns
trkHit.setEDep(tot_edep);// GeV trkHit.setEDep(tot_edep);// GeV
trkHit.setEdx (tot_edep/tot_length); // GeV/mm trkHit.setEdx (tot_edep/tot_length); // GeV/mm
......
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