|
ATLAS Offline Software
|
#include <MuonCalibIntSagittaTool.h>
|
| MuonCalibIntSagittaTool (const std::string &name) |
|
virtual | ~MuonCalibIntSagittaTool ()=default |
|
virtual StatusCode | initialize () override |
| Dummy implementation of the initialisation function. More...
|
|
virtual CorrectionCode | applyCorrection (MCP::MuonObj &mu) const =0 |
| Declare the interface that the class provides. More...
|
|
virtual void | print () const =0 |
| Print the state of the tool. More...
|
|
virtual void | print () const |
| Print the state of the tool. More...
|
|
ServiceHandle< StoreGateSvc > & | evtStore () |
| The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
const ServiceHandle< StoreGateSvc > & | evtStore () const |
| The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
const ServiceHandle< StoreGateSvc > & | detStore () const |
| The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
virtual StatusCode | sysInitialize () override |
| Perform system initialization for an algorithm. More...
|
|
virtual StatusCode | sysStart () override |
| Handle START transition. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | inputHandles () const override |
| Return this algorithm's input handles. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | outputHandles () const override |
| Return this algorithm's output handles. More...
|
|
Gaudi::Details::PropertyBase & | declareProperty (Gaudi::Property< T > &t) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc="none") |
| Declare a new Gaudi property. More...
|
|
void | updateVHKA (Gaudi::Details::PropertyBase &) |
|
MsgStream & | msg () const |
|
MsgStream & | msg (const MSG::Level lvl) const |
|
bool | msgLvl (const MSG::Level lvl) const |
|
|
StatusCode | calcSystematicVariation (const SystematicSet &systConfig, ParameterSetSagitta ¶m) const |
|
void | applySagittaCorrection (double &pt, const std::vector< double > &correction, const int &charge, const bool &isData) const |
|
std::vector< double > | getCorrections (const MCP::TrackCalibObj &mu) const |
|
double | statCombCorrection (const MCP::MuonObj &mu, double corrIDpT, double corrMEpT, double CBpT) const |
|
CorrectionCode | applyStatCombination (AmgVector(5) parsID, const AmgSymMatrix(5)&covID, AmgVector(5) parsMS, const AmgSymMatrix(5)&covMS, int charge, AmgVector(5) &parsCB, AmgSymMatrix(5) &covCB, double &chi2) const |
|
void | renounceArray (SG::VarHandleKeyArray &handlesArray) |
| remove all handles from I/O resolution More...
|
|
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > | renounce (T &h) |
|
void | extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps) |
| Add StoreName to extra input/output deps as needed. More...
|
|
|
SG::ReadHandleKey< xAOD::EventInfo > | m_eventInfo {this, "EventInfoContName", "EventInfo", "event info key"} |
|
Gaudi::Property< std::string > | m_release {this, "release", "Recs2022_03_20", "Release"} |
|
Gaudi::Property< bool > | m_doEtaSagittaSys {this, "doEtaSagittaSys", false, "Do Eta dependant systematic system"} |
|
Gaudi::Property< bool > | m_applyCorrectionOnData {this, "applyCorrectionOnData", true, "If to apply sagitta corrections on data, or take the full effect as systematic"} |
|
Gaudi::Property< bool > | m_doDirectCBCalib {this, "doDirectCBCalib", true, "Apply sagitta corrections on CB tracks, otherwise, on ID+MS"} |
|
Gaudi::Property< bool > | m_useFixedRho {this, "useFixedRho", true, "Force Fixed Rho"} |
|
Gaudi::Property< float > | m_fixedRhoValue {this, "fixedRhoValue", 1, "Value for rho, if it is forced"} |
|
Gaudi::Property< std::string > | m_sysScheme {this, "systematicScheme", "Corr_Scale", "Systematic scheme to be configured"} |
|
Gaudi::Property< float > | m_extraRebiasSys {this, "extraRebias", 0.00002, "Additional systematic applied for the residual bias observed"} |
|
SystematicsCache< ParameterSetSagitta > | m_Parameters {this} |
|
const ParameterSetSagitta * | m_currentParameters {nullptr} |
|
std::map< MCP::DataYear, SagittaCorrConstMap > | m_sagittaCorrConst_ID |
|
std::map< MCP::DataYear, SagittaCorrConstMap > | m_sagittaCorrConst_ME |
|
std::map< MCP::DataYear, SagittaCorrConstMap > | m_sagittaCorrConst_CB |
|
std::map< MCP::DataYear, SagittaCorrConstMap > | m_sagittaCorrConst_mcID |
|
std::map< MCP::DataYear, SagittaCorrConstMap > | m_sagittaCorrConst_mcME |
|
std::map< MCP::DataYear, SagittaCorrConstMap > | m_sagittaCorrConst_mcCB |
|
Definition at line 27 of file MuonCalibIntSagittaTool.h.
◆ SagittaCorrConstMap
◆ StoreGateSvc_t
◆ MuonCalibIntSagittaTool()
CP::MuonCalibIntSagittaTool::MuonCalibIntSagittaTool |
( |
const std::string & |
name | ) |
|
◆ ~MuonCalibIntSagittaTool()
virtual CP::MuonCalibIntSagittaTool::~MuonCalibIntSagittaTool |
( |
| ) |
|
|
virtualdefault |
◆ affectingSystematics()
SystematicSet CP::MuonCalibIntSagittaTool::affectingSystematics |
( |
| ) |
const |
|
overrideprivatevirtual |
the list of all systematics this tool can be affected by
Sagitta Systematis
Implements CP::IReentrantSystematicsTool.
Definition at line 421 of file MuonCalibIntSagittaTool.cxx.
428 result.insert(SystematicVariation(
"MUON_SAGITTA_RHO", 1));
429 result.insert(SystematicVariation(
"MUON_SAGITTA_RHO", -1));
434 result.insert(SystematicVariation(
"MUON_SAGITTA_RESBIAS_NEGETA", 1));
435 result.insert(SystematicVariation(
"MUON_SAGITTA_RESBIAS_NEGETA", -1));
437 result.insert(SystematicVariation(
"MUON_SAGITTA_RESBIAS_POSETA", 1));
438 result.insert(SystematicVariation(
"MUON_SAGITTA_RESBIAS_POSETA", -1));
443 result.insert(SystematicVariation(
"MUON_SAGITTA_RESBIAS", 1));
444 result.insert(SystematicVariation(
"MUON_SAGITTA_RESBIAS", -1));
447 result.insert(SystematicVariation(
"MUON_SAGITTA_DATASTAT", 1));
448 result.insert(SystematicVariation(
"MUON_SAGITTA_DATASTAT", -1));
450 result.insert(SystematicVariation(
"MUON_SAGITTA_GLOBAL", 1));
451 result.insert(SystematicVariation(
"MUON_SAGITTA_GLOBAL", -1));
453 result.insert(SystematicVariation(
"MUON_SAGITTA_PTEXTRA", 1));
454 result.insert(SystematicVariation(
"MUON_SAGITTA_PTEXTRA", -1));
◆ applyCorrection()
Declare the interface that the class provides.
Apply the correction on a modifyable object
◆ applySagittaCorrection()
void CP::MuonCalibIntSagittaTool::applySagittaCorrection |
( |
double & |
pt, |
|
|
const std::vector< double > & |
correction, |
|
|
const int & |
charge, |
|
|
const bool & |
isData |
|
) |
| const |
|
protected |
◆ applyStatCombination()
Definition at line 359 of file MuonCalibIntSagittaTool.cxx.
363 parsID[4] = std::abs(parsID[4]);
364 parsMS[4] = std::abs(parsMS[4]);
368 if (!matID.isInvertible()) {
376 if (!matMS.isInvertible()) {
383 Eigen::FullPivLU<
AmgSymMatrix(5)> matCB(weightID + weightMS);
384 if (!matCB.isInvertible()) {
388 covCB = matCB.inverse();
391 Eigen::FullPivLU<
AmgSymMatrix(5)> matSum(covID + covMS);
392 if (!matSum.isInvertible()) {
400 chi2 = diffPars.transpose() * invCovSum * diffPars;
403 parsCB = covCB * (weightID * parsID + weightMS * parsMS);
407 parsCB[2] -= 2. *
M_PI;
408 else
if (parsCB[2] < -
M_PI)
409 parsCB[2] += 2. *
M_PI;
410 return CorrectionCode::Ok;
◆ applySystematicVariation()
StatusCode CP::MuonCalibIntSagittaTool::applySystematicVariation |
( |
const SystematicSet & |
systConfig | ) |
|
|
overrideprivatevirtual |
effects: configure this tool for the given list of systematic variations.
any requested systematics that are not affecting this tool will be silently ignored (unless they cause other errors). failures: systematic unknown failures: requesting multiple variations on the same systematic (e.g. up & down) failures: requesting an unsupported variation on an otherwise supported systematic (e.g. a 2 sigma variation and the tool only supports 1 sigma variations) failures: unsupported combination of supported systematic failures: other tool specific errors
Implements CP::ISystematicsTool.
Definition at line 546 of file MuonCalibIntSagittaTool.cxx.
◆ calcSystematicVariation()
Definition at line 461 of file MuonCalibIntSagittaTool.cxx.
471 SystematicVariation syst = systConfig.getSystematicByBaseName(
"MUON_SAGITTA_RHO");
475 else if (!syst.empty())
return StatusCode::FAILURE;
478 syst = systConfig.getSystematicByBaseName(
"MUON_SAGITTA_RESBIAS");
481 else if (syst == SystematicVariation(
"MUON_SAGITTA_RESBIAS", -1)) param.SagittaBias =
MCP::SystVariation::Up;
482 else if (!syst.empty())
return StatusCode::FAILURE;
485 syst = systConfig.getSystematicByBaseName(
"MUON_SAGITTA_GLOBAL");
488 else if (syst == SystematicVariation(
"MUON_SAGITTA_GLOBAL", -1)) param.SagittaGlobal =
MCP::SystVariation::Up;
489 else if (!syst.empty())
return StatusCode::FAILURE;
492 syst = systConfig.getSystematicByBaseName(
"MUON_SAGITTA_PTEXTRA");
495 else if (syst == SystematicVariation(
"MUON_SAGITTA_PTEXTRA", -1)) param.SagittaPtExtra =
MCP::SystVariation::Up;
496 else if (!syst.empty())
return StatusCode::FAILURE;
499 syst = systConfig.getSystematicByBaseName(
"MUON_SAGITTA_RESBIAS_POSETA");
501 if (syst == SystematicVariation(
"MUON_SAGITTA_RESBIAS_POSETA", 1))
506 else if (syst == SystematicVariation(
"MUON_SAGITTA_RESBIAS_POSETA", -1))
511 else if (!syst.empty())
return StatusCode::FAILURE;
514 syst = systConfig.getSystematicByBaseName(
"MUON_SAGITTA_RESBIAS_NEGETA");
516 if (syst == SystematicVariation(
"MUON_SAGITTA_RESBIAS_NEGETA", 1))
521 else if (syst == SystematicVariation(
"MUON_SAGITTA_RESBIAS_NEGETA", -1))
526 else if (!syst.empty())
return StatusCode::FAILURE;
529 syst = systConfig.getSystematicByBaseName(
"MUON_SAGITTA_DATASTAT");
531 if (syst == SystematicVariation(
"MUON_SAGITTA_DATASTAT", 1)) param.SagittaDataStat =
MCP::SystVariation::Up;
532 else if (syst == SystematicVariation(
"MUON_SAGITTA_DATASTAT", -1)) param.SagittaDataStat =
MCP::SystVariation::Down;
533 else if (!syst.empty())
return StatusCode::FAILURE;
536 ATH_MSG_DEBUG(
"Systematic variation's parameters, SagittaRho: " << param.SagittaRho);
537 ATH_MSG_DEBUG(
"Systematic variation's parameters, SagittaBias: " << param.SagittaBias);
538 ATH_MSG_DEBUG(
"Systematic variation's parameters, SagittaGlobal: " << param.SagittaGlobal);
539 ATH_MSG_DEBUG(
"Systematic variation's parameters, SagittaPtExtra: " << param.SagittaPtExtra);
540 ATH_MSG_DEBUG(
"Systematic variation's parameters, SagittaDataStat: " << param.SagittaDataStat);
541 ATH_MSG_DEBUG(
"Systematic variation's parameters, SagittaEtaSlice: " << param.SagittaEtaSlice);
543 return StatusCode::SUCCESS;
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleBase
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 245 of file AthCommonDataStore.h.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleKey
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 221 of file AthCommonDataStore.h.
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This is the generic version, for types that do not derive from SG::VarHandleKey
. It just forwards to the base class version of declareProperty
.
Definition at line 333 of file AthCommonDataStore.h.
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This dispatches to either the generic declareProperty
or the one for VarHandle/Key/KeyArray.
Definition at line 352 of file AthCommonDataStore.h.
◆ declareProperty() [6/6]
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ extraDeps_update_handler()
Add StoreName to extra input/output deps as needed.
use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given
◆ getCorrections()
Definition at line 148 of file MuonCalibIntSagittaTool.cxx.
150 double eta = trk.eta;
151 double pT = trk.calib_pt;
170 std::vector<double> corrections;
179 corrections.push_back(corr);
191 ATH_MSG_ERROR(
"Sagitta correction is not applied to data, yet Eta dependant systematics are requested. This configuration is not supported");
206 corrections.push_back(corr);
218 corrections.push_back(corr);
228 double corr = deltas *
scale;
231 corrections.push_back(corr);
245 double deltas = 0.00002;
246 if (
eta > 2 || (
eta > -2 &&
eta < -1.05)) {
248 corr += std::abs(450.0 - 45) / 100 * deltas;
250 corr += std::abs(
pT - 45) / 100 * deltas;
252 if (
eta < -2 || (eta < 2 && eta> 1.5)) {
254 corr += std::abs(450.0 - 45) / 200 * deltas;
256 corr += std::abs(
pT - 45) / 200 * deltas;
260 if (eta < 0 && eta> -0.5) corr += 2.1*deltas;
261 else if (
eta < -1.05) corr += 1.1*deltas;
262 else if (
eta > 0.5 ) corr += 0.8*deltas;
267 corrections.push_back(corr*
scale);
278 corr *= (0.5 *
scale);
280 corrections.push_back(corr);
◆ getKey()
Get the (hashed) key of an object that is in the event store.
This is a bit of a special one. StoreGateSvc
and xAOD::TEvent
both provide ways for getting the SG::sgkey_t
key for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.
In order to allow tools to efficiently perform this operation, they can use this helper function.
- See also
- asg::AsgTool::getName
- Parameters
-
ptr | The bare pointer to the object that the event store should know about |
- Returns
- The hashed key of the object in the store. If not found, an invalid (zero) key.
Definition at line 119 of file AsgTool.cxx.
121 #ifdef XAOD_STANDALONE
124 return evtStore()->event()->getKey( ptr );
127 return (
proxy ==
nullptr ? 0 :
proxy->sgkey() );
128 #endif // XAOD_STANDALONE
◆ getName()
const std::string & asg::AsgTool::getName |
( |
const void * |
ptr | ) |
const |
|
inherited |
Get the name of an object that is / should be in the event store.
This is a bit of a special one. StoreGateSvc
and xAOD::TEvent
both provide ways for getting the std::string
name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.
In order to allow tools to efficiently perform this operation, they can use this helper function.
- See also
- asg::AsgTool::getKey
- Parameters
-
ptr | The bare pointer to the object that the event store should know about |
- Returns
- The string name of the object in the store. If not found, an empty string.
Definition at line 106 of file AsgTool.cxx.
108 #ifdef XAOD_STANDALONE
111 return evtStore()->event()->getName( ptr );
114 static const std::string
dummy =
"";
116 #endif // XAOD_STANDALONE
◆ getProperty()
template<class T >
const T* asg::AsgTool::getProperty |
( |
const std::string & |
name | ) |
const |
|
inherited |
Get one of the tool's properties.
◆ initialize()
StatusCode CP::MuonCalibIntSagittaTool::initialize |
( |
| ) |
|
|
overridevirtual |
Dummy implementation of the initialisation function.
It's here to allow the dual-use tools to skip defining an initialisation function. Since many are doing so...
Reimplemented from asg::AsgTool.
Definition at line 23 of file MuonCalibIntSagittaTool.cxx.
32 for(
const auto&
year: MCP::dataYearList)
52 return StatusCode::FAILURE;
55 if (
registry.registerSystematics(*
this) != StatusCode::SUCCESS)
58 return StatusCode::FAILURE;
61 return StatusCode::SUCCESS;
◆ inputHandles()
Return this algorithm's input handles.
We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.
◆ isAffectedBySystematic()
◆ msg() [1/2]
◆ msg() [2/2]
◆ msg_level_name()
const std::string & asg::AsgTool::msg_level_name |
( |
| ) |
const |
|
inherited |
A deprecated function for getting the message level's name.
Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:
MSG::name( msg().level() )
This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.
- Returns
- The string name of the current minimum message level that's printed
Definition at line 101 of file AsgTool.cxx.
◆ msgLvl()
◆ outputHandles()
Return this algorithm's output handles.
We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.
◆ print() [1/2]
void asg::AsgTool::print |
( |
| ) |
const |
|
virtualinherited |
Print the state of the tool.
Implements asg::IAsgTool.
Reimplemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, top::TopObjectSelection, JetDumper, CP::BaseTrackVertexAssociationTool, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
Definition at line 131 of file AsgTool.cxx.
◆ print() [2/2]
virtual void asg::IAsgTool::print |
( |
| ) |
const |
|
pure virtualinherited |
Print the state of the tool.
Implemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, top::TopObjectSelection, asg::AsgTool, JetDumper, CP::BaseTrackVertexAssociationTool, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
◆ recommendedSystematics()
SystematicSet CP::MuonCalibIntSagittaTool::recommendedSystematics |
( |
| ) |
const |
|
overrideprivatevirtual |
◆ renounce()
◆ renounceArray()
◆ statCombCorrection()
double CP::MuonCalibIntSagittaTool::statCombCorrection |
( |
const MCP::MuonObj & |
mu, |
|
|
double |
corrIDpT, |
|
|
double |
corrMEpT, |
|
|
double |
CBpT |
|
) |
| const |
|
protected |
Definition at line 296 of file MuonCalibIntSagittaTool.cxx.
303 if(
mu.ID.calib_pt == 0)
return CBpT;
304 if(
mu.ME.calib_pt == 0)
return CBpT;
305 if(corrIDpT == 0)
return CBpT;
306 if(corrMEpT == 0)
return CBpT;
308 double chi2Nom = -999;
321 using TLV =
ROOT::Math::LorentzVector<
ROOT::Math::PtEtaPhiM4D<
double>>;
322 TLV tlv{
mu.ID.calib_pt,
mu.ID.eta,
mu.ID.phi,
mu.ID.mass};
324 if(tlv.P() == 0) parsID[4] = 1e12;
325 else parsID[4] = 1.0 / (tlv.P() * 1
e3);
327 tlv.SetCoordinates(
mu.ME.calib_pt,
mu.ME.eta,
mu.ME.phi,
mu.ME.mass);
329 if(tlv.P() == 0) parsMS[4] = 1e12;
330 else parsMS[4] = 1.0 / (tlv.P() * 1
e3);
332 CorrectionCode SysCorrCode =
applyStatCombination(parsID, covID, parsMS, covMS,
mu.CB.calib_charge, parsCBNom, covCBNom, chi2Nom);
338 tlv.SetCoordinates(corrIDpT,
mu.ID.eta,
mu.ID.phi,
mu.ID.mass);
340 if(tlv.P() == 0) parsID[4] = 1e12;
341 else parsID[4] = 1.0 / (tlv.P() * 1
e3);
343 tlv.SetCoordinates(corrMEpT,
mu.ME.eta,
mu.ME.phi,
mu.ME.mass);
345 if(tlv.P() == 0) parsMS[4] = 1e12;
346 else parsMS[4] = 1.0 / (tlv.P() * 1
e3);
348 SysCorrCode =
applyStatCombination(parsID, covID, parsMS, covMS,
mu.CB.calib_charge, parsCBCorr, covCBCorr, chi2Nom);
352 double statCombPtNom =
std::sin(parsCBNom[3]) / std::abs(parsCBNom[4]);
353 double statCombPtCorr =
std::sin(parsCBCorr[3]) / std::abs(parsCBCorr[4]);
354 double corrCBpT = CBpT * (statCombPtCorr / statCombPtNom);
◆ sysInitialize()
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ updateVHKA()
◆ m_applyCorrectionOnData
Gaudi::Property<bool> CP::MuonCalibIntSagittaTool::m_applyCorrectionOnData {this, "applyCorrectionOnData", true, "If to apply sagitta corrections on data, or take the full effect as systematic"} |
|
protected |
◆ m_currentParameters
◆ m_detStore
◆ m_doDirectCBCalib
Gaudi::Property<bool> CP::MuonCalibIntSagittaTool::m_doDirectCBCalib {this, "doDirectCBCalib", true, "Apply sagitta corrections on CB tracks, otherwise, on ID+MS"} |
|
protected |
◆ m_doEtaSagittaSys
Gaudi::Property<bool> CP::MuonCalibIntSagittaTool::m_doEtaSagittaSys {this, "doEtaSagittaSys", false, "Do Eta dependant systematic system"} |
|
protected |
◆ m_eventInfo
◆ m_evtStore
◆ m_extraRebiasSys
Gaudi::Property<float> CP::MuonCalibIntSagittaTool::m_extraRebiasSys {this, "extraRebias", 0.00002, "Additional systematic applied for the residual bias observed"} |
|
protected |
◆ m_fixedRhoValue
Gaudi::Property<float> CP::MuonCalibIntSagittaTool::m_fixedRhoValue {this, "fixedRhoValue", 1, "Value for rho, if it is forced"} |
|
protected |
◆ m_Parameters
◆ m_release
Gaudi::Property<std::string> CP::MuonCalibIntSagittaTool::m_release {this, "release", "Recs2022_03_20", "Release"} |
|
protected |
◆ m_sagittaCorrConst_CB
◆ m_sagittaCorrConst_ID
◆ m_sagittaCorrConst_mcCB
◆ m_sagittaCorrConst_mcID
◆ m_sagittaCorrConst_mcME
◆ m_sagittaCorrConst_ME
◆ m_sysScheme
Gaudi::Property<std::string> CP::MuonCalibIntSagittaTool::m_sysScheme {this, "systematicScheme", "Corr_Scale", "Systematic scheme to be configured"} |
|
protected |
◆ m_useFixedRho
Gaudi::Property<bool> CP::MuonCalibIntSagittaTool::m_useFixedRho {this, "useFixedRho", true, "Force Fixed Rho"} |
|
protected |
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
std::vector< Identifier > ID
Scalar eta() const
pseudorapidity method
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
#define ATH_MSG_VERBOSE(x)
#define AmgSymMatrix(dim)
virtual void setOwner(IDataHandleHolder *o)=0
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
@ Error
Some error happened during the object correction.
double chi2(TH1 *h0, TH1 *h1)
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
double charge(const T &p)
@ Ok
The correction was done successfully.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
std::map< SagittaCorrection, std::shared_ptr< CalibContainer > > createSagittaCorrMap(DataYear dataYear, TrackType type, const std::string &recommendationPath, const std::string &correctionType)
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
static SystematicRegistry & getInstance()
Get the singleton instance of the registry for the curren thread.