ATLAS Offline Software
Loading...
Searching...
No Matches
Muon::MdtCalibJsonDumpAlg Class Reference

#include <MdtCalibJsonDumpAlg.h>

Inheritance diagram for Muon::MdtCalibJsonDumpAlg:

Public Member Functions

virtual ~MdtCalibJsonDumpAlg ()
virtual StatusCode initialize () override final
virtual StatusCode execute (const EventContext &ctx) override final
 Execute method.
virtual StatusCode finalize () override final
 AthHistogramAlgorithm (const std::string &name, ISvcLocator *pSvcLocator)
 Constructor with parameters:
virtual StatusCode sysInitialize ()
 Initialization method invoked by the framework.
const ServiceHandle< ITHistSvc > & histSvc () const
 The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
virtual bool filterPassed (const EventContext &ctx) const
 Get filter decision:
virtual void setFilterPassed (bool state, const EventContext &ctx) const
 Set filter decision:
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Protected Member Functions

virtual bool isReEntrant () const override final
 Legacy algorithms are not thread-safe.
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
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.
StatusCode configAthHistogramming (const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix)
 To be called to fill the internal configuration.
TH1 * bookGetPointer (const TH1 &hist, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 *hist, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 &histRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TTree * bookGetPointer (const TTree &treeRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
TGraph * bookGetPointer (const TGraph &graphRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
TEfficiency * bookGetPointer (const TEfficiency &eff, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency *eff, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency &effRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (const TH1 &hist, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 *hist, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 &histRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (const TTree &treeRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
StatusCode book (const TGraph &graphRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
StatusCode book (const TEfficiency &eff, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency *eff, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency &effRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TH1 * hist (std::string_view histName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered histograms of any type.
TH2 * hist2d (std::string_view histName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered 2-d histograms.
TH3 * hist3d (std::string_view histName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered 3-d histograms.
TTree * tree (std::string_view treeName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered TTrees.
TGraph * graph (std::string_view graphName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered TGraphs.
TEfficiency * efficiency (std::string_view effName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered TEfficiency.

Private Types

using CalibParamSorter = MuonCalib::CalibParamSorter
using SingleTubeCalib = CalibParamSorter::SingleTubeCalib
using RtGrouper = std::map<const MuonCalib::MdtRtRelation*, std::set<Identifier>, CalibParamSorter>
using T0PerChamb = std::map<const SingleTubeCalib*, std::set<Identifier>, CalibParamSorter>
using T0Grouper = std::map<Identifier, T0PerChamb>
typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
typedef uint32_t hash_t
 typedef for the internal hash
typedef std::map< const hash_t, TH1 * > HistMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TEfficiency * > EffMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TTree * > TreeMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TGraph * > GraphMap_t
 Typedef for convenience.

Private Member Functions

std::string dump (const std::vector< double > &values) const
std::string dump (const double v) const
std::string toString (const EventIDRange &range) const
StatusCode dumpRtRelations (const EventContext &ctx, const RtGrouper &rtRelMap, const EventIDRange &eventRange)
 writes the RT relations into the JSON & output root file
StatusCode dumpTubeT0s (const EventContext &ctx, const T0Grouper &t0Map, const EventIDRange &eventRange)
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>
std::string buildBookingString (std::string_view histName, std::string_view tDir, std::string_view stream, bool usePrefixPostfix=false) const
 Method to build individual booking string.
hash_t hash (std::string_view histName) const
 Method to calculate a 32-bit hash from a string.

Private Attributes

ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
SG::ReadCondHandleKey< MuonCalib::MdtCalibDataContainer > m_readKey
Gaudi::Property< std::string > m_rtJSON {this, "RtJSON", "RtRelations.json"}
Gaudi::Property< std::string > m_t0JSON {this, "TubeT0JSON", "TubeT0.json"}
Gaudi::Property< bool > m_savePerIOV {this, "JsonPerIOV", true}
Gaudi::Property< unsigned > m_whiteSpace {this, "WhiteSpaces", 2}
Gaudi::Property< unsigned > m_precision {this, "FloatPrecision", 4}
std::vector< EventIDRange > m_seenIDs {}
MuonVal::MuonTesterTree m_rtDumpTree {"RtCalibConstants", "MDTRTCALIBDUMP"}
MuonVal::ScalarBranch< std::string > & m_rt_type {m_rtDumpTree.newScalar<std::string>("rtType")}
 string defining the type of the rt function
MuonVal::VectorBranch< double > & m_rt_pars {m_rtDumpTree.newVector<double>("rtParams")}
 parameters of the rt-function
MuonVal::ScalarBranch< std::string > & m_tr_type {m_rtDumpTree.newScalar<std::string>("trType", "")}
 String defining the type of the tr function.
MuonVal::VectorBranch< double > & m_tr_pars {m_rtDumpTree.newVector<double>("trParams")}
 parameters of the tr-function
MuonVal::ScalarBranch< std::string > & m_rt_resoType {m_rtDumpTree.newScalar<std::string>("resoType")}
 string defining the type of the rt resolution function
MuonVal::VectorBranch< double > & m_rt_resoPars {m_rtDumpTree.newVector<double>("resoParams")}
 parameters of the rt resolution function
MuonVal::VectorBranch< std::string > & m_rt_stName {m_rtDumpTree.newVector<std::string>("stationName")}
 Station name string of the associated chambers.
MuonVal::VectorBranch< short > & m_rt_stEta {m_rtDumpTree.newVector<short>("stationEta")}
 Station eta of the associated chambers.
MuonVal::VectorBranch< unsigned short > & m_rt_stPhi {m_rtDumpTree.newVector<unsigned short>("stationPhi")}
 Station phi of the associated chamber.
MuonVal::VectorBranch< unsigned short > & m_rt_stMl {m_rtDumpTree.newVector<unsigned short>("multiLayer")}
 Multi layer of the associated chamber.
MuonVal::EventIDBranch m_rt_iov_start {m_rtDumpTree, "startIOV"}
 Data mebmers to save the IOV range (Start point).
MuonVal::EventIDBranch m_rt_iov_end {m_rtDumpTree, "endIOV"}
 Data mebmers to save the IOV range (End point).
MuonVal::MuonTesterTree m_t0DumpTree {"T0CalibConstants", "MDTT0CALIBDUMP"}
MuonVal::ScalarBranch< std::string > & m_t0_stName {m_t0DumpTree.newScalar<std::string>("stationName")}
MuonVal::ScalarBranch< short > & m_t0_stEta {m_t0DumpTree.newScalar<short>("stationEta")}
MuonVal::ScalarBranch< unsigned short > & m_t0_stPhi {m_t0DumpTree.newScalar<unsigned short>("stationPhi")}
MuonVal::ScalarBranch< float > & m_t0_t0 {m_t0DumpTree.newScalar<float>("t0")}
MuonVal::ScalarBranch< float > & m_t0_adc {m_t0DumpTree.newScalar<float>("adc")}
MuonVal::ScalarBranch< unsigned short > & m_t0_code {m_t0DumpTree.newScalar<unsigned short>("code")}
MuonVal::VectorBranch< unsigned short > & m_t0_multiLayer {m_t0DumpTree.newVector<unsigned short>("multiLayer")}
MuonVal::VectorBranch< unsigned short > & m_t0_tubeLayer {m_t0DumpTree.newVector<unsigned short>("tubeLayer")}
MuonVal::VectorBranch< unsigned short > & m_t0_tube {m_t0DumpTree.newVector<unsigned short>("tube")}
MuonVal::EventIDBranch m_t0_iov_start {m_t0DumpTree, "startIOV"}
 Data mebmers to save the IOV range (Start point).
MuonVal::EventIDBranch m_t0_iov_end {m_t0DumpTree, "endIOV"}
 Data mebmers to save the IOV range (End point).
ServiceHandle< ITHistSvc > m_histSvc
 Default constructor: AthHistogramAlgorithm();.
std::string m_prefix
 Name of the ROOT output stream (file).
std::string m_rootDir
 Name of the ROOT directory.
std::string m_histNamePrefix
 The prefix for the histogram THx name.
std::string m_histNamePostfix
 The postfix for the histogram THx name.
std::string m_histTitlePrefix
 The prefix for the histogram THx title.
std::string m_histTitlePostfix
 The postfix for the histogram THx title.
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default).
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared
HistMap_t m_histMap
 The map of histogram names to their pointers.
EffMap_t m_effMap
 The map of histogram names to their pointers.
TreeMap_t m_treeMap
 The map of TTree names to their pointers.
GraphMap_t m_graphMap
 The map of TGraph names to their pointers.
std::string m_streamName
 Name of the ROOT output stream (file).
std::string m_name
 Instance name.
MsgStream m_msg
 Cached Message Stream.

Detailed Description

Definition at line 19 of file MdtCalibJsonDumpAlg.h.

Member Typedef Documentation

◆ CalibParamSorter

◆ EffMap_t

typedef std::map< const hash_t, TEfficiency* > AthHistogramming::EffMap_t
privateinherited

Typedef for convenience.

Definition at line 199 of file AthHistogramming.h.

◆ GraphMap_t

typedef std::map< const hash_t, TGraph* > AthHistogramming::GraphMap_t
privateinherited

Typedef for convenience.

Definition at line 213 of file AthHistogramming.h.

◆ hash_t

typedef uint32_t AthHistogramming::hash_t
privateinherited

typedef for the internal hash

Definition at line 171 of file AthHistogramming.h.

◆ HistMap_t

typedef std::map< const hash_t, TH1* > AthHistogramming::HistMap_t
privateinherited

Typedef for convenience.

Definition at line 192 of file AthHistogramming.h.

◆ RtGrouper

◆ SingleTubeCalib

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

◆ T0Grouper

Definition at line 33 of file MdtCalibJsonDumpAlg.h.

◆ T0PerChamb

Definition at line 32 of file MdtCalibJsonDumpAlg.h.

◆ TreeMap_t

typedef std::map< const hash_t, TTree* > AthHistogramming::TreeMap_t
privateinherited

Typedef for convenience.

Definition at line 206 of file AthHistogramming.h.

Constructor & Destructor Documentation

◆ ~MdtCalibJsonDumpAlg()

Muon::MdtCalibJsonDumpAlg::~MdtCalibJsonDumpAlg ( )
virtualdefault

Member Function Documentation

◆ AthHistogramAlgorithm()

AthHistogramAlgorithm::AthHistogramAlgorithm ( const std::string & name,
ISvcLocator * pSvcLocator )

Constructor with parameters:

Definition at line 40 of file AthHistogramAlgorithm.cxx.

32 :
33 ::AthAlgorithm( name, pSvcLocator ),
34 ::AthHistogramming( name ),
35 m_histSvc ( "THistSvc/THistSvc", name )
36{
37 //
38 // Property declaration
39 //
40 // The histogram service
41 declareProperty("THistSvc",
42 m_histSvc = ServiceHandle<ITHistSvc>("THistSvc/THistSvc", name),
43 "Handle to a THistSvc instance: it will be used to write "
44 "ROOT objects to ROOT files" );
45 // declareProperty("THistService", m_histSvc, "The THistSvc" );
46
47 declareProperty("RootStreamName", m_prefix = "/ANALYSIS", "Name of the output ROOT stream (file) that the THistSvc uses");
48 declareProperty("RootDirName", m_rootDir = "",
49 "Name of the ROOT directory inside the ROOT file where the histograms will go");
50
51 declareProperty( "HistNamePrefix", m_histNamePrefix = "", "The prefix for the histogram THx name" );
52 declareProperty( "HistNamePostfix", m_histNamePostfix = "", "The postfix for the histogram THx name" );
53
54 declareProperty( "HistTitlePrefix", m_histTitlePrefix = "", "The prefix for the histogram THx title" );
55 declareProperty( "HistTitlePostfix", m_histTitlePostfix = "", "The postfix for the histogram THx title" );
56}
AthAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Constructor.
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
std::string m_histNamePrefix
The prefix for the histogram THx name.
ServiceHandle< ITHistSvc > m_histSvc
Default constructor: AthHistogramAlgorithm();.
std::string m_prefix
Name of the ROOT output stream (file).
std::string m_histNamePostfix
The postfix for the histogram THx name.
std::string m_histTitlePostfix
The postfix for the histogram THx title.
std::string m_rootDir
Name of the ROOT directory.
std::string m_histTitlePrefix
The prefix for the histogram THx title.
AthHistogramming(const std::string &name)
Constructor with parameters:

◆ book() [1/8]

StatusCode AthHistogramming::book ( const TEfficiency & eff,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 329 of file AthHistogramming.h.

330{
331 // We need to create a non-const clone
332 TEfficiency* effClone = dynamic_cast< TEfficiency* >( eff.Clone() );
333 if ( !effClone ) {
334 m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone" << endmsg;
335 return StatusCode::FAILURE;
336 }
337 return this->book( *effClone, tDir, stream );
338}
#define endmsg
StatusCode book(const TH1 &hist, std::string_view tDir="", std::string_view stream="")
Simplify the booking and registering (into THistSvc) of histograms.
MsgStream m_msg
Cached Message Stream.

◆ book() [2/8]

StatusCode AthHistogramming::book ( const TGraph & graphRef,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

◆ book() [3/8]

StatusCode AthHistogramming::book ( const TH1 & hist,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 297 of file AthHistogramming.h.

298{
299 // We need to create a non-const clone
300 TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
301 if ( !histClone ) {
302 m_msg << MSG::ERROR << "Couldn't create a TH1 clone" << endmsg;
303 return StatusCode::FAILURE;
304 }
305 return this->book( *histClone, tDir, stream );
306}
TH1 * hist(std::string_view histName, std::string_view tDir="", std::string_view stream="")
Simplify the retrieval of registered histograms of any type.

◆ book() [4/8]

StatusCode AthHistogramming::book ( const TTree & treeRef,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 397 of file AthHistogramming.h.

398{
399 // Call the other Book method and see if it returns a valid pointer
400 TTree* treePointer = this->bookGetPointer( treeRef, tDir, stream );
401 if ( treePointer )
402 {
403 return StatusCode::SUCCESS;
404 }
405 else
406 {
407 return StatusCode::FAILURE;
408 }
409}
TH1 * bookGetPointer(const TH1 &hist, std::string_view tDir="", std::string_view stream="")
Simplify the booking and registering (into THistSvc) of histograms.

◆ book() [5/8]

StatusCode AthHistogramming::book ( TEfficiency & effRef,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 350 of file AthHistogramming.h.

351{
352 // Call the other Book method and see if it returns a valid pointer
353 TEfficiency* effPointer = this->bookGetPointer( effRef, tDir, stream );
354 if ( !effPointer ) {
355 m_msg << MSG::ERROR << "Couldn't book a TEfficiency" << endmsg;
356 return StatusCode::FAILURE;
357 }
358 return StatusCode::SUCCESS;
359}

◆ book() [6/8]

StatusCode AthHistogramming::book ( TEfficiency * eff,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 340 of file AthHistogramming.h.

341{
342 if ( !eff ) {
343 m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency" << endmsg;
344 return StatusCode::FAILURE;
345 }
346 return this->book( *eff, tDir, stream );
347}

◆ book() [7/8]

StatusCode AthHistogramming::book ( TH1 & histRef,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 318 of file AthHistogramming.h.

319{
320 // Call the other Book method and see if it returns a valid pointer
321 TH1* histPointer = this->bookGetPointer( histRef, tDir, stream );
322 if ( !histPointer ) {
323 m_msg << MSG::ERROR << "Couldn't book a TH1" << endmsg;
324 return StatusCode::FAILURE;
325 }
326 return StatusCode::SUCCESS;
327}

◆ book() [8/8]

StatusCode AthHistogramming::book ( TH1 * hist,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 308 of file AthHistogramming.h.

309{
310 if ( !hist ) {
311 m_msg << MSG::ERROR << "Got a zero pointer to a TH1" << endmsg;
312 return StatusCode::FAILURE;
313 }
314 return this->book( *hist, tDir, stream );
315}

◆ bookGetPointer() [1/8]

TEfficiency * AthHistogramming::bookGetPointer ( const TEfficiency & eff,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 275 of file AthHistogramming.h.

276{
277 // We need to create a non-const clone
278 TEfficiency* histClone = dynamic_cast< TEfficiency* >( hist.Clone() );
279 if ( !histClone ) {
280 m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone in bookGetPointer" << endmsg;
281 return 0;
282 }
283 return this->bookGetPointer( *histClone, tDir, stream );
284
285}

◆ bookGetPointer() [2/8]

TGraph * AthHistogramming::bookGetPointer ( const TGraph & graphRef,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

Definition at line 374 of file AthHistogramming.cxx.

375{
376 // Get a pointer
377 const TGraph* graphPointer = &graphRef;
378
379 // Check that we got a valid pointer
380 if ( !graphPointer )
381 {
382 m_msg << MSG::WARNING
383 << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
384 << "!" << endmsg;
385 return NULL;
386 }
387
388 // Modify the name and title according to the prefixes of this classes instance
389 std::string graphName = graphPointer->GetName();
390 const std::string graphTitle = graphPointer->GetTitle();
391
392 // Check if the hash for this graphName already exists, i.e., if we have a hash collision
393 const hash_t graphHash = this->hash(graphName);
394 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
395 if ( it != m_graphMap.end() ) // It does exist!
396 {
397 m_msg << MSG::WARNING
398 << "Detected a hash collision. The hash for the TGraph with name=" << graphName
399 << " already exists and points to a TGraph with name=" << it->second->GetName()
400 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
401 return NULL;
402 }
403
404 // Create a clone that has the new name
405 TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
406 if( !graphClone )
407 {
408 m_msg << MSG::WARNING
409 << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
410 << "!" << endmsg;
411 return NULL;
412 }
413 graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
414
415 // Massage the final string to book things
416 std::string bookingString = buildBookingString( graphName, tDir, stream );
417
418 // Register the TGraph into the THistSvc
419 if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
420 {
421 m_msg << MSG::WARNING
422 << "Problem registering TGraph with name " << graphName
423 << ", title " << graphTitle
424 << " in " << m_name << "!" << endmsg;
425 return NULL;
426 }
427
428 // Also register it in the local map of string to pointer
429 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
430
431 return graphClone;
432}
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
hash_t hash(std::string_view histName) const
Method to calculate a 32-bit hash from a string.
std::string buildBookingString(std::string_view histName, std::string_view tDir, std::string_view stream, bool usePrefixPostfix=false) const
Method to build individual booking string.
std::string m_histNamePostfix
The postfix for the histogram THx name.
uint32_t hash_t
typedef for the internal hash
std::string m_histTitlePostfix
The postfix for the histogram THx title.
std::string m_histTitlePrefix
The prefix for the histogram THx title.
std::string m_name
Instance name.
std::string m_histNamePrefix
The prefix for the histogram THx name.
GraphMap_t m_graphMap
The map of TGraph names to their pointers.

◆ bookGetPointer() [3/8]

TH1 * AthHistogramming::bookGetPointer ( const TH1 & hist,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 254 of file AthHistogramming.h.

255{
256 // We need to create a non-const clone
257 TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
258 if ( !histClone ) {
259 m_msg << MSG::ERROR << "Couldn't create a TH1 clone in bookGetPointer" << endmsg;
260 return 0;
261 }
262 return this->bookGetPointer( *histClone, tDir, stream );
263
264}

◆ bookGetPointer() [4/8]

TTree * AthHistogramming::bookGetPointer ( const TTree & treeRef,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 268 of file AthHistogramming.cxx.

269{
270 // Get a pointer
271 const TTree* treePointer = &treeRef;
272
273 // Check that we got a valid pointer
274 if ( !treePointer )
275 {
276 m_msg << MSG::WARNING
277 << "We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" << m_name
278 << "!" << endmsg;
279 return NULL;
280 }
281
282 // Modify the name and title according to the prefixes of this classes instance
283 std::string treeName = treePointer->GetName();
284 const std::string treeTitle = treePointer->GetTitle();
285
286 // Check if the hash for this treeName already exists, i.e., if we have a hash collision
287 const hash_t treeHash = this->hash(treeName);
288 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
289 if ( it != m_treeMap.end() ) // It does exist!
290 {
291 m_msg << MSG::WARNING
292 << "Detected a hash collision. The hash for the TTree with name=" << treeName
293 << " already exists and points to a TTree with name=" << it->second->GetName()
294 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
295 return NULL;
296 }
297
298 // Create a clone that has the new name
299 TTree* treeClone = dynamic_cast< TTree* >( treePointer->Clone(treeName.c_str()) );
300 if( !treeClone )
301 {
302 m_msg << MSG::WARNING
303 << "We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" << m_name
304 << "!" << endmsg;
305 return NULL;
306 }
307 treeClone->SetTitle (treeTitle.c_str());
308
309 // Massage the final string to book things
310 std::string bookingString = buildBookingString( treeName, tDir, stream );
311
312 // Register the TTree into the THistSvc
313 if ( !((histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
314 {
315 m_msg << MSG::WARNING
316 << "Problem registering TTree with name " << treeName
317 << ", title " << treeTitle
318 << " in " << m_name << "!" << endmsg;
319 return NULL;
320 }
321
322 // Also register it in the local map of string to pointer
323 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
324
325 return treeClone;
326}
TreeMap_t m_treeMap
The map of TTree names to their pointers.

◆ bookGetPointer() [5/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency & effRef,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 139 of file AthHistogramming.cxx.

140{
141 std::string originalName{effRef.GetName()};
142 const auto effName = splitName(originalName).second;
143 const std::string effTitle(effRef.GetTitle());
144 std::string bookingString = buildBookingString(originalName, tDir, stream);
145 std::string finalEffName{m_histNamePrefix};
146 finalEffName += effName;
147 finalEffName += m_histNamePostfix;
148 effRef.SetTitle((m_histTitlePrefix + effTitle + m_histTitlePostfix).c_str());
149 effRef.SetName(finalEffName.c_str());
150
151 const hash_t effHash = this->hash(effName);
152 EffMap_t::const_iterator it = m_effMap.find(effHash);
153 if (it != m_effMap.end()) {
154 m_msg << MSG::WARNING
155 << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
156 << " already exists and points to a TEfficiency with name=" << it->second->GetName()
157 << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
158 return nullptr;
159 }
160
161 if (!histSvc()->regEfficiency(bookingString, &effRef).isSuccess()) {
162 m_msg << MSG::WARNING
163 << "Problem registering TEfficiency with name " << effName
164 << ", name prefix " << m_histNamePrefix
165 << ", title " << effTitle
166 << ", title prefix " << m_histTitlePrefix
167 << ", and title postfix " << m_histTitlePostfix
168 << " in " << m_name << "!" << endmsg;
169 return nullptr;
170 }
171
172 m_effMap.insert(m_effMap.end(), std::pair<const hash_t, TEfficiency*>(effHash, &effRef));
173
174 return &effRef;
175}
EffMap_t m_effMap
The map of histogram names to their pointers.

◆ bookGetPointer() [6/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency * eff,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 287 of file AthHistogramming.h.

288{
289 if ( !hist ) {
290 m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency in bookGetPointer" << endmsg;
291 return 0;
292 }
293 return this->bookGetPointer( *hist, tDir, stream );
294}

◆ bookGetPointer() [7/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 & histRef,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 97 of file AthHistogramming.cxx.

97 {
98
99 const std::string originalName{histRef.GetName()};
100 const auto histName = splitName(originalName).second;
101
102 const std::string histTitle(histRef.GetTitle());
103 std::string bookingString = buildBookingString(originalName, tDir, stream);
104
105 std::string finalHistName{m_histNamePrefix};
106 finalHistName += histName;
107 finalHistName += m_histNamePostfix;
108
109 histRef.SetTitle((m_histTitlePrefix + histTitle + m_histTitlePostfix).c_str());
110 histRef.SetName(finalHistName.c_str());
111
112 const hash_t histHash = this->hash(histName);
113 HistMap_t::const_iterator it = m_histMap.find(histHash);
114 if (it != m_histMap.end()) {
115 m_msg << MSG::WARNING
116 << "Detected a hash collision. The hash for the histogram with name=" << histName
117 << " already exists and points to a histogram with name=" << it->second->GetName()
118 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
119 return nullptr;
120 }
121
122 if (!histSvc()->regHist(bookingString, &histRef).isSuccess()) {
123 m_msg << MSG::WARNING
124 << "Problem registering histogram with name " << histName
125 << ", name prefix " << m_histNamePrefix
126 << ", title " << histTitle
127 << ", title prefix " << m_histTitlePrefix
128 << ", and title postfix " << m_histTitlePostfix
129 << " in " << m_name << "!" << endmsg;
130 return nullptr;
131 }
132
133 m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, &histRef));
134
135 return &histRef;
136}
HistMap_t m_histMap
The map of histogram names to their pointers.

◆ bookGetPointer() [8/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 * hist,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 266 of file AthHistogramming.h.

267{
268 if ( !hist ) {
269 m_msg << MSG::ERROR << "Got a zero pointer to a TH1 in bookGetPointer" << endmsg;
270 return 0;
271 }
272 return this->bookGetPointer( *hist, tDir, stream );
273}

◆ buildBookingString()

std::string AthHistogramming::buildBookingString ( std::string_view histName,
std::string_view tDir,
std::string_view stream,
bool usePrefixPostfix = false ) const
privateinherited

Method to build individual booking string.

Definition at line 500 of file AthHistogramming.cxx.

501 {
502 if (tDir.empty()) {
503 tDir = m_rootDir;
504 }
505 if (stream.empty()) {
507 }
508 const auto [histDir, baseName] = splitName(histName);
509 std::string bookingString{"/"};
510 bookingString += stream;
511 bookingString += '/';
512 bookingString += tDir;
513 if (!histDir.empty()) {
514 bookingString += '/';
515 bookingString += histDir;
516 }
517 bookingString += '/';
518 if (usePrefixPostfix) {
519 bookingString += m_histNamePrefix;
520 }
521 bookingString += baseName;
522 if (usePrefixPostfix) {
523 bookingString += m_histNamePostfix;
524 }
525 const auto tail = std::ranges::unique(
526 bookingString, [](char lhs, char rhs) { return lhs == '/' && rhs == '/'; });
527 bookingString.erase(tail.begin(), tail.end());
528 return bookingString;
529}
std::string m_rootDir
Name of the ROOT directory.
std::string m_streamName
Name of the ROOT output stream (file).
std::string tail(std::string s, const std::string &pattern)
tail of a string

◆ configAthHistogramming()

StatusCode AthHistogramming::configAthHistogramming ( const ServiceHandle< ITHistSvc > & histSvc,
const std::string & prefix,
const std::string & rootDir,
const std::string & histNamePrefix,
const std::string & histNamePostfix,
const std::string & histTitlePrefix,
const std::string & histTitlePostfix )
protectedinherited

To be called to fill the internal configuration.

Definition at line 64 of file AthHistogramming.cxx.

68{
71 m_rootDir = rootDir;
72 m_histNamePrefix = histNamePrefix;
73 m_histNamePostfix = histNamePostfix;
74 m_histTitlePrefix = histTitlePrefix;
75 m_histTitlePostfix = histTitlePostfix;
76
77 return StatusCode::SUCCESS;
78}
ServiceHandle< ITHistSvc > m_histSvc
Pointer to the THistSvc (event store by default).

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl,
const SG::VarHandleKeyType &  )
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

◆ dump() [1/2]

std::string Muon::MdtCalibJsonDumpAlg::dump ( const double v) const
private

Definition at line 214 of file MdtCalibJsonDumpAlg.cxx.

214 {
215 std::stringstream sstr{};
216 if (std::abs(v)< std::pow(0.1, m_precision + 1)) {
217 sstr<<0;
218 } else {
219 sstr<<std::setprecision(m_precision + std::max(std::floor(std::log10(std::abs(v))), 1.))<<v;
220 }
221 return sstr.str();
222 }
Gaudi::Property< unsigned > m_precision

◆ dump() [2/2]

std::string Muon::MdtCalibJsonDumpAlg::dump ( const std::vector< double > & values) const
private

Definition at line 204 of file MdtCalibJsonDumpAlg.cxx.

204 {
205 std::stringstream sstr{};
206 sstr<<"[";
207 for (const auto [count, v] : Acts::enumerate(values)) {
208 sstr<<dump(v);
209 if (count + 1 != values.size())sstr<<", ";
210 }
211 sstr<<"]";
212 return sstr.str();
213 }
std::string dump(const std::vector< double > &values) const
int count(std::string s, const std::string &regx)
count how many occurances of a regx are in a string
Definition hcg.cxx:148

◆ dumpRtRelations()

StatusCode Muon::MdtCalibJsonDumpAlg::dumpRtRelations ( const EventContext & ctx,
const RtGrouper & rtRelMap,
const EventIDRange & eventRange )
private

writes the RT relations into the JSON & output root file

write the RT relation

write the TR relation

Write the resolution

Definition at line 140 of file MdtCalibJsonDumpAlg.cxx.

140 {
141 std::ofstream outFileRT{toString(eventRange) + m_rtJSON};
142 if (!outFileRT.good()) {
143 ATH_MSG_FATAL("Failed to write "<<m_rtJSON);
144 return StatusCode::FAILURE;
145 }
146 outFileRT<<"[\n";
147 const MdtIdHelper& idHelper{m_idHelperSvc->mdtIdHelper()};
148 nlohmann::json rtDump{};
149 for (const auto [counter, payload] : Acts::enumerate(rtRelMap)) {
150 outFileRT<<whiteSpace(m_whiteSpace)<<"{\n";
151 outFileRT<<whiteSpace(m_whiteSpace)<<" \"chambers\": [\n";
152 const MdtRtRelation* rtRel = payload.first;
153 const std::set<Identifier>& chambers = payload.second;
154
155 for (const auto [chCounter, chambId] : Acts::enumerate(chambers)) {
156 outFileRT<<whiteSpace(2*m_whiteSpace)<<"{"
157 <<std::format("\"station\": \"{:}\", \"eta\": {:2d}, \"phi\": {:1d}, \"ml\": {:1d}",
158 idHelper.stationNameString(idHelper.stationName(chambId)), idHelper.stationEta(chambId),
159 idHelper.stationPhi(chambId), idHelper.multilayer(chambId))<<"}";
160
161 m_rt_stName.push_back(idHelper.stationNameString(idHelper.stationName(chambId)));
162 m_rt_stEta.push_back(idHelper.stationEta(chambId));
163 m_rt_stPhi.push_back(idHelper.stationPhi(chambId));
164 m_rt_stMl.push_back(idHelper.multilayer(chambId));
165
166 if (chCounter +1 != chambers.size()) {
167 outFileRT<<",\n";
168 } else {
169 outFileRT<<"],\n";
170 }
171 }
173 outFileRT<<whiteSpace(m_whiteSpace+1)<<"\"rtRelation\": {\n";
174 outFileRT<<whiteSpace(2*m_whiteSpace)<<"\"type\": \""<<rtRel->rt()->name()<<"\",\n";
175 m_rt_type = rtRel->rt()->name();
176 m_rt_pars = rtRel->rt()->parameters();
177 outFileRT<<whiteSpace(2*m_whiteSpace)<<"\"params\": "<<dump(rtRel->rt()->parameters())<<"},\n";
179 if (rtRel->tr()->name() != "TrRelationLookUp") {
180 outFileRT<<whiteSpace(m_whiteSpace+1)<<"\"trRelation\": {\n";
181 outFileRT<<whiteSpace(2*m_whiteSpace)<<"\"type\": \""<<rtRel->tr()->name()<<"\",\n";
182 outFileRT<<whiteSpace(2*m_whiteSpace)<<"\"params\": "<<dump(rtRel->tr()->parameters())<<"},\n";
183 m_tr_type = rtRel->tr()->name();
184 m_tr_pars = rtRel->tr()->parameters();
185 }
187 outFileRT<<whiteSpace(m_whiteSpace+1)<<"\"rtReso\": {\n";
188 outFileRT<<whiteSpace(2*m_whiteSpace)<<"\"type\": \""<<rtRel->rtRes()->name()<<"\",\n";
189 outFileRT<<whiteSpace(2*m_whiteSpace)<<"\"params\": "<<dump(rtRel->rtRes()->parameters())<<"}\n";
190 m_rt_resoType = rtRel->rtRes()->name();
191 m_rt_resoPars = rtRel->rtRes()->parameters();
192 if (counter +1 != rtRelMap.size()) {
193 outFileRT<<whiteSpace(m_whiteSpace)<<"},\n";
194 } else {
195 outFileRT<<whiteSpace(m_whiteSpace)<<"}\n";
196 }
197 m_rt_iov_start = eventRange.start();
198 m_rt_iov_end = eventRange.stop();
199 ATH_CHECK(m_rtDumpTree.fill(ctx));
200 }
201 outFileRT<<"]\n";
202 return StatusCode::SUCCESS;
203 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_FATAL(x,...)
int multilayer(const Identifier &id) const
Access to components of the ID.
const ParVec & parameters() const
Definition CalibFunc.h:40
virtual std::string name() const =0
const ITrRelation * tr() const
t(r) relationship
const IRtResolution * rtRes() const
resolution
const IRtRelation * rt() const
rt relation
int stationEta(const Identifier &id) const
const std::string & stationNameString(const Identifier &id) const
int stationPhi(const Identifier &id) const
int stationName(const Identifier &id) const
MuonVal::VectorBranch< double > & m_rt_resoPars
parameters of the rt resolution function
Gaudi::Property< std::string > m_rtJSON
MuonVal::ScalarBranch< std::string > & m_rt_type
string defining the type of the rt function
MuonVal::VectorBranch< unsigned short > & m_rt_stPhi
Station phi of the associated chamber.
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
MuonVal::EventIDBranch m_rt_iov_end
Data mebmers to save the IOV range (End point).
MuonVal::VectorBranch< double > & m_rt_pars
parameters of the rt-function
MuonVal::ScalarBranch< std::string > & m_tr_type
String defining the type of the tr function.
std::string toString(const EventIDRange &range) const
MuonVal::MuonTesterTree m_rtDumpTree
MuonVal::VectorBranch< unsigned short > & m_rt_stMl
Multi layer of the associated chamber.
MuonVal::EventIDBranch m_rt_iov_start
Data mebmers to save the IOV range (Start point).
Gaudi::Property< unsigned > m_whiteSpace
MuonVal::ScalarBranch< std::string > & m_rt_resoType
string defining the type of the rt resolution function
MuonVal::VectorBranch< short > & m_rt_stEta
Station eta of the associated chambers.
MuonVal::VectorBranch< double > & m_tr_pars
parameters of the tr-function
MuonVal::VectorBranch< std::string > & m_rt_stName
Station name string of the associated chambers.

◆ dumpTubeT0s()

StatusCode Muon::MdtCalibJsonDumpAlg::dumpTubeT0s ( const EventContext & ctx,
const T0Grouper & t0Map,
const EventIDRange & eventRange )
private

Definition at line 223 of file MdtCalibJsonDumpAlg.cxx.

223 {
224 std::ofstream outFileT0{toString(eventRange) + m_t0JSON};
225 if (!outFileT0.good()) {
226 ATH_MSG_FATAL("Failed to write "<<m_t0JSON);
227 return StatusCode::FAILURE;
228 }
229 outFileT0<<"[\n";
230 const MdtIdHelper& idHelper{m_idHelperSvc->mdtIdHelper()};
231 for (const auto [count, payload] : Acts::enumerate(t0Map)) {
232 outFileT0<<whiteSpace(m_whiteSpace)<<"{\n";
233 const Identifier& detId = payload.first;
234 const T0PerChamb& calibConstants = payload.second;
235
236 m_t0_stName = idHelper.stationNameString(idHelper.stationName(detId));
237 m_t0_stEta = idHelper.stationEta(detId);
238 m_t0_stPhi = idHelper.stationPhi(detId);
239
240 outFileT0<<whiteSpace(m_whiteSpace+1)
241 <<std::format("\"station\": \"{:}\", \"eta\": {:2d}, \"phi\": {:1d}, \"calibConstants\":[\n",
242 idHelper.stationNameString(idHelper.stationName(detId)),
243 idHelper.stationEta(detId), idHelper.stationPhi(detId));
244 for (const auto [calibCounter, calibPayload] : Acts::enumerate(calibConstants)) {
245 outFileT0<<whiteSpace(2*m_whiteSpace)<<"{\"t0\": "<<dump(calibPayload.first->t0)<<", ";
246 outFileT0<<"\"adc\": "<<dump(calibPayload.first->adcCal)<<", ";
247 outFileT0<<"\"code\": "<<dump(calibPayload.first->statusCode)<<", ";
248 outFileT0<<"\"tubes\": [";
249 Identifier refId = (*calibPayload.second.begin());
250
251 m_t0_t0= calibPayload.first->t0;
252 m_t0_adc = calibPayload.first->adcCal;
253 m_t0_code = calibPayload.first->statusCode;
254 std::vector<u_short> tubes{};
255 auto dumpTubeLayer = [&tubes, &outFileT0, &idHelper, &refId](const Identifier& newLayId) {
256 outFileT0<<std::format("{{\"ml\": {:1d}, \"tl\": {:1d}, \"no\": \"{:}\"}}",
257 idHelper.multilayer(refId), idHelper.tubeLayer(refId),
258 rangeString(tubes));
259 tubes.clear();
260 refId = newLayId;
261 };
262 for (const Identifier& tubeId : calibPayload.second){
263 const int ml = idHelper.multilayer(tubeId);
264 const int tl = idHelper.tubeLayer(tubeId);
265 m_t0_multiLayer.push_back(ml);
266 m_t0_tubeLayer.push_back(tl);
267 m_t0_tube.push_back(idHelper.tube(tubeId));
268 if (tl != idHelper.tubeLayer(refId) || ml != idHelper.multilayer(refId)) {
269 dumpTubeLayer(tubeId);
270 outFileT0<<", ";
271 }
272 tubes.push_back(idHelper.tube(tubeId));
273 }
274 dumpTubeLayer(refId);
275 outFileT0<<"]";
276 if (calibCounter +1 != calibConstants.size()){
277 outFileT0<<"},\n";
278 } else {
279 outFileT0<<"}]"<<std::endl;
280 }
281 }
282
283 if (count +1 != t0Map.size()){
284 outFileT0<<whiteSpace(m_whiteSpace)<<"},\n";
285 } else {
286 outFileT0<<whiteSpace(m_whiteSpace)<<"}"<<std::endl;
287 }
288 m_t0_iov_start = eventRange.start();
289 m_t0_iov_end = eventRange.stop();
290 ATH_CHECK(m_t0DumpTree.fill(ctx));
291 }
292 outFileT0<<"]\n";
293 return StatusCode::SUCCESS;
294 }
std::map< const SingleTubeCalib *, std::set< Identifier >, CalibParamSorter > T0PerChamb
MuonVal::ScalarBranch< short > & m_t0_stEta
Gaudi::Property< std::string > m_t0JSON
MuonVal::VectorBranch< unsigned short > & m_t0_tubeLayer
MuonVal::ScalarBranch< float > & m_t0_adc
MuonVal::ScalarBranch< float > & m_t0_t0
MuonVal::EventIDBranch m_t0_iov_start
Data mebmers to save the IOV range (Start point).
MuonVal::ScalarBranch< unsigned short > & m_t0_code
MuonVal::ScalarBranch< unsigned short > & m_t0_stPhi
MuonVal::VectorBranch< unsigned short > & m_t0_multiLayer
MuonVal::VectorBranch< unsigned short > & m_t0_tube
MuonVal::MuonTesterTree m_t0DumpTree
MuonVal::ScalarBranch< std::string > & m_t0_stName
MuonVal::EventIDBranch m_t0_iov_end
Data mebmers to save the IOV range (End point).

◆ efficiency()

TEfficiency * AthHistogramming::efficiency ( std::string_view effName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered TEfficiency.

Definition at line 214 of file AthHistogramming.cxx.

215{
216 // Build a 32 bit hash out of the name
217 const hash_t effHash = this->hash(effName);
218
219 // See if this entry exists in the map
220 EffMap_t::const_iterator it = m_effMap.find( effHash );
221 if ( it == m_effMap.end() ) // It doesn't exist!
222 { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
223 // Massage the final string to book things
224 std::string bookingString = buildBookingString(effName, tDir, stream ,false);
225
226 TEfficiency* effPointer{};
227 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
228 {
229 // Book things
230 std::string bookingString = buildBookingString( effName, tDir, stream, true );
231
232 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
233 {
234 m_msg << MSG::WARNING
235 << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
236 << bookingString
237 << " or with name " << effName
238 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
239 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
240 return NULL;
241 }
242 // If we get to here, we actually found the TEfficiency in the THistSvc.
243 // So let's add it to the local cache map and return its pointer
244 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
245 return effPointer;
246 }
247 // If we get to here, we actually found the TEfficiency in the THistSvc.
248 // So let's add it to the local cache map and return its pointer
249 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
250 return effPointer;
251 }
252
253 // Return the pointer to the TEfficiency that we got from the local cache map
254 return it->second;
255}
std::pair< StatusCode, TEfficiency * > getEfficiency(ITHistSvc &svc, const std::string &name)

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode Muon::MdtCalibJsonDumpAlg::execute ( const EventContext & ctx)
finaloverridevirtual

Execute method.

Implements AthAlgorithm.

Definition at line 80 of file MdtCalibJsonDumpAlg.cxx.

80 {
81 SG::ReadCondHandle calibHandle{m_readKey, ctx};
82 ATH_CHECK(calibHandle.isValid());
83
84 if (std::ranges::find(m_seenIDs, calibHandle.getRange()) != m_seenIDs.end()) {
85 ATH_MSG_INFO("IOV "<<calibHandle.getRange()<<" has already been dumped");
86 return StatusCode::SUCCESS;
87 }
88 m_seenIDs.emplace_back(calibHandle.getRange());
89 const MdtIdHelper &idHelper{m_idHelperSvc->mdtIdHelper()};
90
91 T0Grouper t0Groups{};
92 const CalibParamSorter sorter{std::pow(0.1, m_precision +2)};
93 RtGrouper rtGroups{sorter};
94 unsigned int rtCounter{0}, t0Counter{0};
95 for (auto itr = idHelper.detectorElement_begin(); itr != idHelper.detectorElement_end(); ++itr) {
96 const Identifier &detId{*itr};
97 const MdtFullCalibData *calibData = calibHandle->getCalibData(detId, msgStream());
98 if (!calibData || !calibData->rtRelation) {
99 continue;
100 }
101 ++rtCounter;
102 ATH_MSG_VERBOSE("Group rt constants of "<<m_idHelperSvc->toStringDetEl(detId));
103 rtGroups[calibData->rtRelation.get()].insert(detId);
104 const int ml = idHelper.multilayer(detId);
105
106 T0PerChamb& chambT0{t0Groups.insert(std::make_pair(idHelper.elementID(detId), T0PerChamb{sorter})).first->second};
107
108 for (int layer = 1; layer <= idHelper.tubeLayerMax(detId); ++layer) {
109 for (int tube = 1; tube <= idHelper.tubeMax(detId); ++tube) {
110 bool valid{false};
111 const Identifier chId{idHelper.channelID(detId, ml, layer, tube, valid)};
112 if (!valid) {
113 continue;
114 }
115 const MuonCalib::MdtTubeCalibContainer::SingleTubeCalib* singleCalib = calibData->tubeCalib->getCalib(chId);
116 if (!singleCalib) {
117 continue;
118 }
119 chambT0[singleCalib].insert(chId);
120 ++t0Counter;
121 }
122 }
123 }
124 ATH_CHECK(dumpRtRelations(ctx, rtGroups, calibHandle.getRange()));
125 ATH_CHECK(dumpTubeT0s(ctx, t0Groups, calibHandle.getRange()));
126
127 unsigned t0GrpCounter{0};
128 for (const auto&[chId, t0Channels] : t0Groups) {
129 t0GrpCounter+=t0Channels.size();
130 }
131 if ((rtCounter == 0) or (t0Counter == 0)){
132 ATH_MSG_ERROR("MdtCalibJsonDumpAlg::execute: Counter is zero in numerator");
133 return StatusCode::FAILURE;
134 }
135 ATH_MSG_INFO(std::format("Grouped {:d} / {:d} ({:.2f}%) rt relations & {:d}/ {:d} ({:.2f}%) t0 calibration constants",
136 rtGroups.size(), rtCounter, (100.* rtGroups.size() / rtCounter),
137 t0GrpCounter, t0Counter,(100.*t0GrpCounter / t0Counter)));
138 return StatusCode::SUCCESS;
139 }
#define ATH_MSG_ERROR(x,...)
#define ATH_MSG_VERBOSE(x,...)
#define ATH_MSG_INFO(x,...)
Obj * get()
Get (non-const) access to the underlying object.
Identifier elementID(int stationName, int stationEta, int stationPhi) const
static int tubeLayerMax()
int tubeMax() const
Identifier channelID(int stationName, int stationEta, int stationPhi, int multilayer, int tubeLayer, int tube) const
const SingleTubeCalib * getCalib(const Identifier &tubeId) const
return calibration constants of a single tube
const_id_iterator detectorElement_begin() const
Iterators over full set of ids.
const_id_iterator detectorElement_end() const
StatusCode dumpRtRelations(const EventContext &ctx, const RtGrouper &rtRelMap, const EventIDRange &eventRange)
writes the RT relations into the JSON & output root file
std::map< Identifier, T0PerChamb > T0Grouper
MuonCalib::CalibParamSorter CalibParamSorter
SG::ReadCondHandleKey< MuonCalib::MdtCalibDataContainer > m_readKey
std::vector< EventIDRange > m_seenIDs
std::map< const MuonCalib::MdtRtRelation *, std::set< Identifier >, CalibParamSorter > RtGrouper
StatusCode dumpTubeT0s(const EventContext &ctx, const T0Grouper &t0Map, const EventIDRange &eventRange)
const EventIDRange & getRange()
@ layer
Definition HitInfo.h:79
list valid
Definition calibdata.py:44

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps)
protectedinherited

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

◆ extraOutputDeps()

const DataObjIDColl & AthCommonAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 89 of file AthCommonAlgorithm.cxx.

54{
55 // If we didn't find any symlinks to add, just return the collection
56 // from the base class. Otherwise, return the extended collection.
57 if (!m_extendedExtraObjects.empty()) {
59 }
61}
Common base class for algorithms.

◆ filterPassed()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx) const
inlinevirtualinherited

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

93 {
94 return execState( ctx ).filterPassed();
95 }
virtual bool filterPassed(const EventContext &ctx) const
Get filter decision:

◆ finalize()

StatusCode Muon::MdtCalibJsonDumpAlg::finalize ( )
finaloverridevirtual

Definition at line 75 of file MdtCalibJsonDumpAlg.cxx.

75 {
76 ATH_CHECK(m_rtDumpTree.write());
77 ATH_CHECK(m_t0DumpTree.write());
78 return StatusCode::SUCCESS;
79 }

◆ graph()

TGraph * AthHistogramming::graph ( std::string_view graphName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered TGraphs.

Definition at line 438 of file AthHistogramming.cxx.

439{
440 // Build a 32 bit hash out of the name
441 const hash_t graphHash = this->hash(graphName);
442
443 // See if this entry exists in the map
444 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
445 if ( it == m_graphMap.end() ) // It doesn't exist!
446 { // Let's see into the THistSvc if somebody else has registered the TGraph...
447
448 // Massage the final string to book things
449 std::string bookingString = buildBookingString( graphName, tDir, stream, false);
450
451 TGraph* graphPointer(nullptr);
452 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
453 {
454 // Massage the final string to book things
455 std::string bookingString = buildBookingString( graphName, tDir, stream, true );
456
457 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
458 {
459 m_msg << MSG::WARNING
460 << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
461 << bookingString
462 << " or with name " << graphName
463 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
464 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
465 return nullptr;
466 }
467 // If we get to here, we actually found the TGraph in the THistSvc.
468 // So let's add it to the local cache map and return its pointer
469 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
470 return graphPointer;
471 }
472 // If we get to here, we actually found the TGraph in the THistSvc.
473 // So let's add it to the local cache map and return its pointer
474 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
475 return graphPointer;
476 }
477
478
479 // Return the pointer to the TGraph that we got from the local cache map
480 return it->second;
481}
std::pair< StatusCode, TGraph * > getGraph(ITHistSvc &svc, const std::string &name)

◆ hash()

AthHistogramming::hash_t AthHistogramming::hash ( std::string_view histName) const
inlineprivateinherited

Method to calculate a 32-bit hash from a string.

Definition at line 422 of file AthHistogramming.h.

423{
424 const uint64_t hash64 = CxxUtils::crc64( std::string{histName} );
425 return (hash_t)(hash64 & 0xFFFFFFFF);
426}
uint64_t crc64(const CRCTable &table, const char *data, size_t data_len)
Find the CRC-64 of a string,.
Definition crc64.cxx:695
std::uint64_t hash64(const void *data, std::size_t size)
Passthrough to XXH3_64bits.
Definition XXH.cxx:9

◆ hist()

TH1 * AthHistogramming::hist ( std::string_view histName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered histograms of any type.

Definition at line 182 of file AthHistogramming.cxx.

183 {
184
185 const auto histBaseName = splitName(histName).second;
186 std::string bookingString = buildBookingString(histName, tDir, stream, false);
187
188 const hash_t histHash = this->hash(histBaseName);
189
190 HistMap_t::const_iterator it = m_histMap.find(histHash);
191 if (it == m_histMap.end()) {
192 TH1* histPointer(nullptr);
193
194 if (!histSvc()->getHist(bookingString, histPointer).isSuccess()) {
195 std::string prefixedBookingString = buildBookingString(histName, tDir, stream, true);
196 if (!histSvc()->getHist(prefixedBookingString, histPointer).isSuccess()) {
197 m_msg << MSG::WARNING
198 << "Problem retrieving the histogram with name (including pre- and post-fixes) "
199 << prefixedBookingString
200 << " or with name " << histName
201 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
202 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
203 return nullptr;
204 }
205 }
206
207 m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, histPointer));
208 return histPointer;
209 }
210
211 return it->second;
212}
std::pair< StatusCode, TH1 * > getHist(ITHistSvc &svc, const std::string &name, size_t index=0)

◆ hist2d()

TH2 * AthHistogramming::hist2d ( std::string_view histName,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the retrieval of registered 2-d histograms.

Definition at line 363 of file AthHistogramming.h.

364{
365 // Get the TH1 pointer
366 TH1* th1Pointer = this->hist(histName, tDir, stream);
367 if ( !th1Pointer )
368 {
369 m_msg << MSG::ERROR
370 << "Cannot get a 2-d histogram with name " << histName
371 << "... will probably seg-fault!" << endmsg;
372 return NULL;
373 }
374 // If the TH1 pointer is valid, simply return the dynamic_cast
375 return dynamic_cast<TH2*>( th1Pointer );
376}

◆ hist3d()

TH3 * AthHistogramming::hist3d ( std::string_view histName,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the retrieval of registered 3-d histograms.

Definition at line 380 of file AthHistogramming.h.

381{
382 // Get the TH1 pointer
383 TH1* th1Pointer = this->hist(histName, tDir, stream);
384 if ( !th1Pointer )
385 {
386 m_msg << MSG::ERROR
387 << "Cannot get a 3-d histogram with name " << histName
388 << "... will probably seg-fault!" << endmsg;
389 return NULL;
390 }
391 // If the TH1 pointer is valid, simply return the dynamic_cast
392 return dynamic_cast<TH3*>( th1Pointer );
393}

◆ histSvc()

const ServiceHandle< ITHistSvc > & AthHistogramAlgorithm::histSvc ( ) const
inlineinherited

The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc.

Definition at line 113 of file AthHistogramAlgorithm.h.

114{
115 return m_histSvc;
116}

◆ initialize()

StatusCode Muon::MdtCalibJsonDumpAlg::initialize ( )
finaloverridevirtual

Definition at line 59 of file MdtCalibJsonDumpAlg.cxx.

59 {
60 ATH_CHECK(m_readKey.initialize());
61 ATH_CHECK(m_idHelperSvc.retrieve());
62 if (m_rtJSON.value().empty() || m_t0JSON.value().empty()) {
63 ATH_MSG_FATAL("No output file has been given " << m_rtJSON<<" "<<m_t0JSON);
64 return StatusCode::FAILURE;
65 }
66 if (m_t0JSON.value() == m_rtJSON.value()){
67 ATH_MSG_FATAL("You can't dump the two streams into the same JSON "<<m_t0JSON<<" & "<<m_rtJSON);
68 return StatusCode::FAILURE;
69 }
70 ATH_CHECK(m_rtDumpTree.init(this));
71 ATH_CHECK(m_t0DumpTree.init(this));
72 ATH_MSG_INFO("Dump calibration constants in " << m_rtJSON);
73 return StatusCode::SUCCESS;
74 }

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::inputHandles ( ) const
overridevirtualinherited

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.

◆ isClonable()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const
inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

68 {
69 return true;
70 }

◆ isReEntrant()

virtual bool AthAlgorithm::isReEntrant ( ) const
inlinefinaloverrideprotectedvirtualinherited

Legacy algorithms are not thread-safe.

Definition at line 47 of file AthAlgorithm.h.

47{ return false; }

◆ msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::msgLvl ( const MSG::Level lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles ( ) const
overridevirtualinherited

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.

◆ renounce()

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 > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
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)

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

virtual void AthCommonAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state,
const EventContext & ctx ) const
inlinevirtualinherited

Set filter decision:

Definition at line 99 of file AthCommonAlgorithm.h.

99 {
101 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const
Set filter decision:

◆ sysExecute()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx)
overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

41{
42 return BaseAlg::sysExecute (ctx);
43}

◆ sysInitialize()

StatusCode AthHistogramAlgorithm::sysInitialize ( )
virtualinherited

Initialization method invoked by the framework.

This method is responsible for any bookkeeping of initialization required by the framework itself. It will in turn invoke the initialize() method of the derived algorithm, and of any sub-algorithms which it creates.

Reimplemented from AthCommonAlgorithm< Gaudi::Algorithm >.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 75 of file AthHistogramAlgorithm.cxx.

76{
77 // ---- stolen from GaudiKernel/Algorithm::sysInitialize -------
78 // Bypass the initialization if the algorithm
79 // has already been initialized.
80 if ( Gaudi::StateMachine::INITIALIZED <= FSMState() ) return StatusCode::SUCCESS;
81
82 // Set the Algorithm's properties
83 bindPropertiesTo( serviceLocator()->getOptsSvc() );
84
85 // Bypass the initialization if the algorithm is disabled.
86 // Need to do this after bindPropertiesTo.
87 if ( !isEnabled( ) ) return StatusCode::SUCCESS;
88
89 // ---- stolen from GaudiKernel/Algorithm::sysInitialize ------- END ---
90
91
92 // Get the THistSvc
93 ATH_CHECK ( histSvc().retrieve() );
94
95 // Configure the underlying AthHistogramming helper
100
101 // Print some setup information into the log file
102 ATH_MSG_DEBUG ("Initializing " << name() << "...");
103 ATH_MSG_DEBUG (" using THistService = " << m_histSvc );
104 ATH_MSG_DEBUG (" using RootStreamName = " << m_prefix );
105 ATH_MSG_DEBUG (" using RootDirName = " << m_rootDir );
106 ATH_MSG_DEBUG (" using HistNamePrefix = " << m_histNamePrefix );
107 ATH_MSG_DEBUG (" using HistNamePostfix = " << m_histNamePostfix );
108 ATH_MSG_DEBUG (" using HistTitlePrefix = " << m_histTitlePrefix );
109 ATH_MSG_DEBUG (" using HistTitlePostfix = " << m_histTitlePostfix );
110
111
112 // re-direct to base class...
114}
#define ATH_MSG_DEBUG(x,...)
virtual StatusCode sysInitialize() override
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
StatusCode configAthHistogramming(const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix)
To be called to fill the internal configuration.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ toString()

std::string Muon::MdtCalibJsonDumpAlg::toString ( const EventIDRange & range) const
private

Definition at line 295 of file MdtCalibJsonDumpAlg.cxx.

295 {
296 if (m_savePerIOV){
297 std::stringstream sstr{};
298 sstr<<range;
299 std::string outStr{sstr.str()};
300 outStr.erase(std::remove_if(outStr.begin(),outStr.end(),
301 [](const unsigned char c){
302 return !std::isalnum(c) && c!='-' && c !='_';
303 }), outStr.end());
304 return outStr+"_";
305 }
306 return {};
307 }
Gaudi::Property< bool > m_savePerIOV
DataModel_detail::iterator< DVL > remove_if(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end, Predicate pred)
Specialization of remove_if for DataVector/List.

◆ tree()

TTree * AthHistogramming::tree ( std::string_view treeName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered TTrees.

Definition at line 333 of file AthHistogramming.cxx.

334{
335 // Build a 32 bit hash out of the name
336 const hash_t treeHash = this->hash(treeName);
337 // See if this entry exists in the map
338 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
339 if ( it == m_treeMap.end() ) // It doesn't exist!
340 { // Let's see into the THistSvc if somebody else has registered the TTree...
341 // Massage the final string to book things
342 std::string bookingString = buildBookingString( treeName, tDir, stream);
343
344 TTree* treePointer(NULL);
345 if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
346 {
347 m_msg << MSG::WARNING
348 << "Problem retrieving the TTree with name " << treeName
349 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
350 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
351 return NULL;
352 }
353 // If we get to here, we actually found the TTree in the THistSvc.
354 // So let's add it to the local cache map and return its pointer
355 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
356 return treePointer;
357 }
358
359 // Return the pointer to the TTree that we got from the local cache map
360 return it->second;
361}
std::pair< StatusCode, TTree * > getTree(ITHistSvc &svc, const std::string &name)

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

◆ m_effMap

EffMap_t AthHistogramming::m_effMap
privateinherited

The map of histogram names to their pointers.

Definition at line 202 of file AthHistogramming.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

◆ m_graphMap

GraphMap_t AthHistogramming::m_graphMap
privateinherited

The map of TGraph names to their pointers.

Definition at line 216 of file AthHistogramming.h.

◆ m_histMap

HistMap_t AthHistogramming::m_histMap
privateinherited

The map of histogram names to their pointers.

Definition at line 195 of file AthHistogramming.h.

◆ m_histNamePostfix

std::string AthHistogramAlgorithm::m_histNamePostfix
privateinherited

The postfix for the histogram THx name.

Definition at line 97 of file AthHistogramAlgorithm.h.

◆ m_histNamePrefix

std::string AthHistogramAlgorithm::m_histNamePrefix
privateinherited

The prefix for the histogram THx name.

Definition at line 94 of file AthHistogramAlgorithm.h.

◆ m_histSvc

ServiceHandle<ITHistSvc> AthHistogramAlgorithm::m_histSvc
privateinherited

Default constructor: AthHistogramAlgorithm();.

a handle on the Hist/TTree registration service

Definition at line 83 of file AthHistogramAlgorithm.h.

◆ m_histTitlePostfix

std::string AthHistogramAlgorithm::m_histTitlePostfix
privateinherited

The postfix for the histogram THx title.

Definition at line 103 of file AthHistogramAlgorithm.h.

◆ m_histTitlePrefix

std::string AthHistogramAlgorithm::m_histTitlePrefix
privateinherited

The prefix for the histogram THx title.

Definition at line 100 of file AthHistogramAlgorithm.h.

◆ m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> Muon::MdtCalibJsonDumpAlg::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
private

Definition at line 45 of file MdtCalibJsonDumpAlg.h.

45{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

◆ m_msg

MsgStream AthHistogramming::m_msg
privateinherited

Cached Message Stream.

Definition at line 242 of file AthHistogramming.h.

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 239 of file AthHistogramming.h.

◆ m_precision

Gaudi::Property<unsigned> Muon::MdtCalibJsonDumpAlg::m_precision {this, "FloatPrecision", 4}
private

Definition at line 54 of file MdtCalibJsonDumpAlg.h.

54{this, "FloatPrecision", 4};

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix
privateinherited

Name of the ROOT output stream (file).

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_readKey

SG::ReadCondHandleKey<MuonCalib::MdtCalibDataContainer> Muon::MdtCalibJsonDumpAlg::m_readKey
private
Initial value:
{this, "CalibKey", "MdtCalibConstants",
"Conditions object containing the calibrations"}

Definition at line 46 of file MdtCalibJsonDumpAlg.h.

46 {this, "CalibKey", "MdtCalibConstants",
47 "Conditions object containing the calibrations"};

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir
privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_rt_iov_end

MuonVal::EventIDBranch Muon::MdtCalibJsonDumpAlg::m_rt_iov_end {m_rtDumpTree, "endIOV"}
private

Data mebmers to save the IOV range (End point).

Definition at line 83 of file MdtCalibJsonDumpAlg.h.

83{m_rtDumpTree, "endIOV"};

◆ m_rt_iov_start

MuonVal::EventIDBranch Muon::MdtCalibJsonDumpAlg::m_rt_iov_start {m_rtDumpTree, "startIOV"}
private

Data mebmers to save the IOV range (Start point).

Definition at line 81 of file MdtCalibJsonDumpAlg.h.

81{m_rtDumpTree, "startIOV"};

◆ m_rt_pars

MuonVal::VectorBranch<double>& Muon::MdtCalibJsonDumpAlg::m_rt_pars {m_rtDumpTree.newVector<double>("rtParams")}
private

parameters of the rt-function

Definition at line 63 of file MdtCalibJsonDumpAlg.h.

63{m_rtDumpTree.newVector<double>("rtParams")};

◆ m_rt_resoPars

MuonVal::VectorBranch<double>& Muon::MdtCalibJsonDumpAlg::m_rt_resoPars {m_rtDumpTree.newVector<double>("resoParams")}
private

parameters of the rt resolution function

Definition at line 71 of file MdtCalibJsonDumpAlg.h.

71{m_rtDumpTree.newVector<double>("resoParams")};

◆ m_rt_resoType

MuonVal::ScalarBranch<std::string>& Muon::MdtCalibJsonDumpAlg::m_rt_resoType {m_rtDumpTree.newScalar<std::string>("resoType")}
private

string defining the type of the rt resolution function

Definition at line 69 of file MdtCalibJsonDumpAlg.h.

69{m_rtDumpTree.newScalar<std::string>("resoType")};

◆ m_rt_stEta

MuonVal::VectorBranch<short>& Muon::MdtCalibJsonDumpAlg::m_rt_stEta {m_rtDumpTree.newVector<short>("stationEta")}
private

Station eta of the associated chambers.

Definition at line 75 of file MdtCalibJsonDumpAlg.h.

75{m_rtDumpTree.newVector<short>("stationEta")};

◆ m_rt_stMl

MuonVal::VectorBranch<unsigned short>& Muon::MdtCalibJsonDumpAlg::m_rt_stMl {m_rtDumpTree.newVector<unsigned short>("multiLayer")}
private

Multi layer of the associated chamber.

Definition at line 79 of file MdtCalibJsonDumpAlg.h.

79{m_rtDumpTree.newVector<unsigned short>("multiLayer")};

◆ m_rt_stName

MuonVal::VectorBranch<std::string>& Muon::MdtCalibJsonDumpAlg::m_rt_stName {m_rtDumpTree.newVector<std::string>("stationName")}
private

Station name string of the associated chambers.

Definition at line 73 of file MdtCalibJsonDumpAlg.h.

73{m_rtDumpTree.newVector<std::string>("stationName")};

◆ m_rt_stPhi

MuonVal::VectorBranch<unsigned short>& Muon::MdtCalibJsonDumpAlg::m_rt_stPhi {m_rtDumpTree.newVector<unsigned short>("stationPhi")}
private

Station phi of the associated chamber.

Definition at line 77 of file MdtCalibJsonDumpAlg.h.

77{m_rtDumpTree.newVector<unsigned short>("stationPhi")};

◆ m_rt_type

MuonVal::ScalarBranch<std::string>& Muon::MdtCalibJsonDumpAlg::m_rt_type {m_rtDumpTree.newScalar<std::string>("rtType")}
private

string defining the type of the rt function

Definition at line 61 of file MdtCalibJsonDumpAlg.h.

61{m_rtDumpTree.newScalar<std::string>("rtType")};

◆ m_rtDumpTree

MuonVal::MuonTesterTree Muon::MdtCalibJsonDumpAlg::m_rtDumpTree {"RtCalibConstants", "MDTRTCALIBDUMP"}
private

Definition at line 59 of file MdtCalibJsonDumpAlg.h.

59{"RtCalibConstants", "MDTRTCALIBDUMP"};

◆ m_rtJSON

Gaudi::Property<std::string> Muon::MdtCalibJsonDumpAlg::m_rtJSON {this, "RtJSON", "RtRelations.json"}
private

Definition at line 49 of file MdtCalibJsonDumpAlg.h.

49{this, "RtJSON", "RtRelations.json"};

◆ m_savePerIOV

Gaudi::Property<bool> Muon::MdtCalibJsonDumpAlg::m_savePerIOV {this, "JsonPerIOV", true}
private

Definition at line 51 of file MdtCalibJsonDumpAlg.h.

51{this, "JsonPerIOV", true};

◆ m_seenIDs

std::vector<EventIDRange> Muon::MdtCalibJsonDumpAlg::m_seenIDs {}
private

Definition at line 57 of file MdtCalibJsonDumpAlg.h.

57{};

◆ m_streamName

std::string AthHistogramming::m_streamName
privateinherited

Name of the ROOT output stream (file).

Definition at line 220 of file AthHistogramming.h.

◆ m_t0_adc

MuonVal::ScalarBranch<float>& Muon::MdtCalibJsonDumpAlg::m_t0_adc {m_t0DumpTree.newScalar<float>("adc")}
private

Definition at line 91 of file MdtCalibJsonDumpAlg.h.

91{m_t0DumpTree.newScalar<float>("adc")};

◆ m_t0_code

MuonVal::ScalarBranch<unsigned short>& Muon::MdtCalibJsonDumpAlg::m_t0_code {m_t0DumpTree.newScalar<unsigned short>("code")}
private

Definition at line 92 of file MdtCalibJsonDumpAlg.h.

92{m_t0DumpTree.newScalar<unsigned short>("code")};

◆ m_t0_iov_end

MuonVal::EventIDBranch Muon::MdtCalibJsonDumpAlg::m_t0_iov_end {m_t0DumpTree, "endIOV"}
private

Data mebmers to save the IOV range (End point).

Definition at line 100 of file MdtCalibJsonDumpAlg.h.

100{m_t0DumpTree, "endIOV"};

◆ m_t0_iov_start

MuonVal::EventIDBranch Muon::MdtCalibJsonDumpAlg::m_t0_iov_start {m_t0DumpTree, "startIOV"}
private

Data mebmers to save the IOV range (Start point).

Definition at line 98 of file MdtCalibJsonDumpAlg.h.

98{m_t0DumpTree, "startIOV"};

◆ m_t0_multiLayer

MuonVal::VectorBranch<unsigned short>& Muon::MdtCalibJsonDumpAlg::m_t0_multiLayer {m_t0DumpTree.newVector<unsigned short>("multiLayer")}
private

Definition at line 94 of file MdtCalibJsonDumpAlg.h.

94{m_t0DumpTree.newVector<unsigned short>("multiLayer")};

◆ m_t0_stEta

MuonVal::ScalarBranch<short>& Muon::MdtCalibJsonDumpAlg::m_t0_stEta {m_t0DumpTree.newScalar<short>("stationEta")}
private

Definition at line 87 of file MdtCalibJsonDumpAlg.h.

87{m_t0DumpTree.newScalar<short>("stationEta")};

◆ m_t0_stName

MuonVal::ScalarBranch<std::string>& Muon::MdtCalibJsonDumpAlg::m_t0_stName {m_t0DumpTree.newScalar<std::string>("stationName")}
private

Definition at line 86 of file MdtCalibJsonDumpAlg.h.

86{m_t0DumpTree.newScalar<std::string>("stationName")};

◆ m_t0_stPhi

MuonVal::ScalarBranch<unsigned short>& Muon::MdtCalibJsonDumpAlg::m_t0_stPhi {m_t0DumpTree.newScalar<unsigned short>("stationPhi")}
private

Definition at line 88 of file MdtCalibJsonDumpAlg.h.

88{m_t0DumpTree.newScalar<unsigned short>("stationPhi")};

◆ m_t0_t0

MuonVal::ScalarBranch<float>& Muon::MdtCalibJsonDumpAlg::m_t0_t0 {m_t0DumpTree.newScalar<float>("t0")}
private

Definition at line 90 of file MdtCalibJsonDumpAlg.h.

90{m_t0DumpTree.newScalar<float>("t0")};

◆ m_t0_tube

MuonVal::VectorBranch<unsigned short>& Muon::MdtCalibJsonDumpAlg::m_t0_tube {m_t0DumpTree.newVector<unsigned short>("tube")}
private

Definition at line 96 of file MdtCalibJsonDumpAlg.h.

96{m_t0DumpTree.newVector<unsigned short>("tube")};

◆ m_t0_tubeLayer

MuonVal::VectorBranch<unsigned short>& Muon::MdtCalibJsonDumpAlg::m_t0_tubeLayer {m_t0DumpTree.newVector<unsigned short>("tubeLayer")}
private

Definition at line 95 of file MdtCalibJsonDumpAlg.h.

95{m_t0DumpTree.newVector<unsigned short>("tubeLayer")};

◆ m_t0DumpTree

MuonVal::MuonTesterTree Muon::MdtCalibJsonDumpAlg::m_t0DumpTree {"T0CalibConstants", "MDTT0CALIBDUMP"}
private

Definition at line 85 of file MdtCalibJsonDumpAlg.h.

85{"T0CalibConstants", "MDTT0CALIBDUMP"};

◆ m_t0JSON

Gaudi::Property<std::string> Muon::MdtCalibJsonDumpAlg::m_t0JSON {this, "TubeT0JSON", "TubeT0.json"}
private

Definition at line 50 of file MdtCalibJsonDumpAlg.h.

50{this, "TubeT0JSON", "TubeT0.json"};

◆ m_tr_pars

MuonVal::VectorBranch<double>& Muon::MdtCalibJsonDumpAlg::m_tr_pars {m_rtDumpTree.newVector<double>("trParams")}
private

parameters of the tr-function

Definition at line 67 of file MdtCalibJsonDumpAlg.h.

67{m_rtDumpTree.newVector<double>("trParams")};

◆ m_tr_type

MuonVal::ScalarBranch<std::string>& Muon::MdtCalibJsonDumpAlg::m_tr_type {m_rtDumpTree.newScalar<std::string>("trType", "")}
private

String defining the type of the tr function.

Empty if it's a lookup table

Definition at line 65 of file MdtCalibJsonDumpAlg.h.

65{m_rtDumpTree.newScalar<std::string>("trType", "")};

◆ m_treeMap

TreeMap_t AthHistogramming::m_treeMap
privateinherited

The map of TTree names to their pointers.

Definition at line 209 of file AthHistogramming.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_whiteSpace

Gaudi::Property<unsigned> Muon::MdtCalibJsonDumpAlg::m_whiteSpace {this, "WhiteSpaces", 2}
private

Definition at line 53 of file MdtCalibJsonDumpAlg.h.

53{this, "WhiteSpaces", 2};

The documentation for this class was generated from the following files: