Skip to content
Snippets Groups Projects
Commit 1db278ad authored by Stefano Rosati's avatar Stefano Rosati
Browse files

update MM clusters uncertainty

parent 2a842b55
No related branches found
No related tags found
6 merge requests!58791DataQualityConfigurations: Modify L1Calo config for web display,!46784MuonCondInterface: Enable thread-safety checking.,!46776Updated LArMonitoring config file for WD to match new files produced using MT,!45405updated ART test cron job,!42417Draft: DIRE and VINCIA Base Fragments for Pythia 8.3,!39435SWEEP of Fixes in NSW segment finding and in treatment of MM stereo clusters (Sweep of MR !38871, !39241 and !39423)
......@@ -11,7 +11,6 @@ Muon::SimpleMMClusterBuilderTool::SimpleMMClusterBuilderTool(const std::string&
AthAlgTool(t,n,p) {
declareInterface<IMMClusterBuilderTool>(this);
declareProperty("writeStripProperties", m_writeStripProperties = true ); // true for debugging; needs to become false for large productions
declareProperty("useErrorParametrization", m_useErrorParametrization = true);
declareProperty("maxHoleSize", m_maxHoleSize = 1);
}
......@@ -164,23 +163,15 @@ StatusCode Muon::SimpleMMClusterBuilderTool::getClusters(std::vector<Muon::MMPre
Amg::Vector2D clusterLocalPosition(weightedPosX,posY);
double covX = MMprds[j].localCovariance()(Trk::locX, Trk::locX);
///
/// memory allocated dynamically for the PrepRawData is managed by Event Store
///
Amg::MatrixX* covN = new Amg::MatrixX(1,1);
covN->setIdentity();
if(!m_useErrorParametrization) {
(*covN)(0,0) = 6.*(nmerge + 1.)*covX;
if(nmerge<=1) (*covN)(0,0) = covX;
} else {
double localUncertainty = 0.074+0.66*theta-0.15*theta*theta;
if ( m_idHelperSvc->mmIdHelper().isStereo(MMprds[i].identify()) ) {
localUncertainty = 10.;
}
(*covN)(0,0) = localUncertainty * localUncertainty;
}
/// keeping 50% looser uncertainties until the final calibrations will be in place
double localUncertainty = 1.5*(0.074+0.66*theta-0.15*theta*theta);
(*covN)(0,0) = localUncertainty * localUncertainty;
ATH_MSG_VERBOSE(" make merged prepData at strip " << m_idHelperSvc->mmIdHelper().channel(MMprds[j].identify()) << " nmerge " << nmerge << " sqrt covX " << sqrt((*covN)(0,0)));
///
......
......@@ -35,7 +35,6 @@ namespace Muon
ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};
bool m_writeStripProperties;
bool m_useErrorParametrization;
uint m_maxHoleSize;
};
......
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