Measurement* getMeasurement(bhep::hit& hit) { EVector hit_pos(2,0); EVector meas_pos(3,0); meas_pos[0] = hit.x()[0]; meas_pos[1] = hit.x()[1]; meas_pos[2] = hit.x()[2]; hit_pos[0] = meas_pos[0]; hit_pos[1] = meas_pos[1]; //covariance and meastype hardwired for now. EMatrix cov(2,2,0); cov[0][0] = 1.; cov[1][1] = 1.; string meastype = "xy"; Measurement* me = new Measurement(); me->set_name(meastype); me->set_hv(HyperVector(hit_pos,cov)); me->set_name("volume", "Detector"); me->set_position( meas_pos ); //Add the hit energy deposit as a key to the Measurement. const dict::Key Edep = "E_dep"; const dict::Key EdepVal = hit.data("E_dep"); me->set_name(Edep, EdepVal); return me; }
//-------------------------------------------------------------------- RP::Trajectory* KFSvcManager::CreateTrajectory(std::vector<const Hit* > hits, std::vector<double> hitErrors) //-------------------------------------------------------------------- { log4cpp::Category& klog = log4cpp::Category::getRoot(); klog << log4cpp::Priority::DEBUG << " In CreateTrajectory --> " ; klog << log4cpp::Priority::DEBUG << "Cov Matrix --> " ; EVector m(3,0); fCov = EMatrix(3,3,0); for (size_t i=0; i<3; ++i) { fCov[i][i] = hitErrors[i]; } klog << log4cpp::Priority::DEBUG << "Cov Matrix --> " << fCov; // destroy measurements in fMeas container //stc_tools::destroy(fMeas); RP::measurement_vector fMeas; auto size = hits.size(); klog << log4cpp::Priority::DEBUG << " size of hits --> " << size; klog << log4cpp::Priority::DEBUG << " fill measurement vector --> " << size; for(auto hit : hits) { m[0]=hit->XYZ().X(); m[1]=hit->XYZ().Y(); m[2]=hit->XYZ().Z(); klog << log4cpp::Priority::DEBUG << "+++hit x = " << m[0] << " y = " << m[1] << " z = " << m[2]; klog << log4cpp::Priority::DEBUG << "Create a measurement " ; Measurement* meas = new Measurement(); string type=RP::xyz; meas->set_name(type); // the measurement type meas->set_hv(HyperVector(m,fCov,type)); // the HyperVector meas->set_deposited_energy(hit->Edep()); // the deposited energy meas->set_name(RP::setup_volume,"mother"); // the volume // the position of the measurement klog << log4cpp::Priority::DEBUG << "Create measurement plane " ; meas->set_position_hv(HyperVector(m,EMatrix(3,3,0),RP::xyz)); // Add the measurement to the vector fMeas.push_back(meas); } // add the measurements to the trajectory klog << log4cpp::Priority::INFO << " Measurement vector created " ; Trajectory* trj = new Trajectory(); trj->add_constituents(fMeas); klog << log4cpp::Priority::INFO << " Trajectory --> " << *trj ; return trj; }