Commit cb9311d6 authored by Lennart Huth's avatar Lennart Huth
Browse files

added a DUT option for the millepede alignment

parent 87276e6c
......@@ -15,6 +15,7 @@ using namespace std;
AlignmentMillepede::AlignmentMillepede(Configuration config, std::vector<std::shared_ptr<Detector>> detectors)
: Module(std::move(config), std::move(detectors)) {
m_excludeDUT = m_config.get<bool>("excludeDUT", false);
m_numberOfTracksForAlignment = m_config.get<size_t>("number_of_tracks", 20000);
m_dofs = m_config.getArray<bool>("dofs", {});
m_nIterations = m_config.get<size_t>("iterations", 5);
......@@ -99,11 +100,11 @@ void AlignmentMillepede::finalise() {
size_t nPlanes = num_detectors();
for(const auto& det : get_detectors()) {
if(det->isDUT()) {
if(det->isDUT() && m_excludeDUT) {
nPlanes--;
}
}
LOG(INFO) << "Aligning " << nPlanes << "planes";
const size_t nParameters = 6 * nPlanes;
for(unsigned int iteration = 0; iteration < m_nIterations; ++iteration) {
// Define the constraint equations.
......@@ -160,6 +161,7 @@ void AlignmentMillepede::finalise() {
if(converg < m_convergence)
break;
}
// save the new geometry?
}
//=============================================================================
......
......@@ -153,6 +153,8 @@ namespace corryvreckan {
std::vector<unsigned int> m_fixedPlanes;
/// Flag to fix all degrees of freedom or only the translations.
bool m_fix_all;
/// It can be also reasonable to include the DUT in the alignemnt
bool m_excludeDUT;
};
} // namespace corryvreckan
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment