ATLAS Offline Software
Loading...
Searching...
No Matches
FixLArElecCalib Class Reference

fix electronic calibration constants for MC. More...

#include <FixLArElecCalib.h>

Inheritance diagram for FixLArElecCalib:

Public Member Functions

 FixLArElecCalib (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~FixLArElecCalib ()=default
virtual StatusCode initialize () override
virtual StatusCode execute () override
virtual StatusCode finalize () override
virtual StatusCode stop ATLAS_NOT_THREAD_SAFE () override
StatusCode fix1 (const LArOnOffIdMapping *cabling)
StatusCode fix2 (const LArOnOffIdMapping *cabling)
StatusCode fix3 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix4 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix5 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix6 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix7 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix8 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix9 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix10 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix11 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix12 (const LArOnOffIdMapping *cabling)
StatusCode fix13 (const LArOnOffIdMapping *cabling)
StatusCode fix14 (const LArOnOffIdMapping *cabling)
StatusCode fix15 ATLAS_NOT_THREAD_SAFE ()
StatusCode fix16 ATLAS_NOT_THREAD_SAFE ()
StatusCode fix17 ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode fix18 (const LArCalibLineMapping *clmap)
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
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

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.

Private Types

typedef std::pair< HWIdentifier, std::vector< float > > ROW_t
typedef std::vector< ROW_tVROW_t
typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

StatusCode ReadFile (const std::string &filename, const LArOnOffIdMapping *cabling, bool EM, bool withGain, int nvar)
StatusCode ReadFileAll (const std::string &filename, const LArOnOffIdMapping *cabling, bool withGain, int nvar)
StatusCode update_DAC2uA (const std::string &em_filename, const std::string &hec_filename)
template<class T>
StatusCode update_EM_HEC ATLAS_NOT_THREAD_SAFE (const std::string &em_filename, const std::string &hec_filename, const LArOnOffIdMapping *cabling, bool withGain, int nvar)
template<class T>
StatusCode update_All ATLAS_NOT_THREAD_SAFE (const std::string &filename, const LArOnOffIdMapping *cabling, bool withGain, int nvar)
void print_object (const std::string &msg, const LArNoiseMC::LArCondObj &obj)
void set_object (LArNoiseMC::LArCondObj &obj, const std::vector< float > &v)
void print_object (const std::string &msg, const LArAutoCorrMC::LArCondObj &obj)
void set_object (LArAutoCorrMC::LArCondObj &obj, const std::vector< float > &v)
void print_object (const std::string &msg, const LArRampMC::LArCondObj &obj)
void set_object (LArRampMC::LArCondObj &obj, const std::vector< float > &v)
void print_object (const std::string &msg, const LArShape32MC::LArCondObj &obj)
void set_object (LArShape32MC::LArCondObj &obj, const std::vector< float > &v)
void print_object (const std::string &msg, const LArSingleFloatP &obj)
void set_object (LArSingleFloatP &obj, const std::vector< float > &v)
StatusCode addMphysOverMcal (const LArOnOffIdMapping *cabling)
StatusCode fixDACuAMeV ATLAS_NOT_THREAD_SAFE ()
StatusCode updateEMfSampl ATLAS_NOT_THREAD_SAFE (const std::string &filename, const LArOnOffIdMapping *cabling)
StatusCode updateMinBias ATLAS_NOT_THREAD_SAFE (const std::string &filename, const LArOnOffIdMapping *cabling)
StatusCode updateEM_DACuAMeV ATLAS_NOT_THREAD_SAFE (const std::string &filename, const LArOnOffIdMapping *cabling)
StatusCode updateHADfSampl ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

SG::ReadCondHandleKey< LArOnOffIdMappingm_cablingKey {this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"}
SG::ReadCondHandleKey< LArCalibLineMappingm_CLKey {this, "CalibLineKey", "LArCalibLineMap", "calib line key"}
IntegerProperty m_fixFlag {this, "FixFlag", 1, "which fix to apply"}
StringProperty m_g4Phys {this, "G4Phys", "", "which physics list to use"}
StringProperty m_infile {this, "InputFile", "", "name of inpput file"}
const LArEM_IDm_em_idhelper
const LArHEC_IDm_hec_idhelper
const LArFCAL_IDm_fcal_idhelper
const LArOnlineIDm_online_idhelper
const LArEM_SuperCell_IDm_sem_idhelper
const LArHEC_SuperCell_IDm_shec_idhelper
const LArOnline_SuperCellIDm_sonline_idhelper
const CaloCell_SuperCell_IDm_scell_idhelper
VROW_t m_cache [3]
DataObjIDColl m_extendedExtraObjects
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

Detailed Description

fix electronic calibration constants for MC.

Author
Hong Ma hma@b.nosp@m.nl.g.nosp@m.ov

Purpose: read in the LAr electronic calib data, and fix them ready to be written out again as a new version.

This should be used with FixLArElecCalib_jobOptions.py to write a new POOL file.

Implemented fixes:

  1. Flag = 1 Replacing values for barrel in LArDAC2uAMC Input barrel_dac2ua_oct06.txt
  2. Flag=2 fix FCAL conditions data for new mapping in release 13.
  3. Flag=3 replace fsampl for 4.8.2 June 13,2007: FCAL numbers updated by P. Krieger. Junly 10, 2007: updated the number for 13.0.20
  4. Flag=4 replace uA2MeV DAC2uA with that in dac2ua_ua2mev_rel13.txt
  5. Flag=5 replace FCAL noise and minbiasRMS from Sven Menke's new number for release 13.

6 Flag=6 replace All FCAL noise and ADC2DAC, uA2MeV and DAC2uA. <HG/MG> = 9.96 +/- 0.01 <MG/LG> = 9.67 +/- 0.01

7 Flag=7 replace sampling fraction for QGSP_BERT from 13.0.30.

8 Flag=8 replace EM/HEC autoCorr, noise, dac2ua, mphys, ramp, shape, ua2mev

9 Flag=9 fix FCAL AutoCorr

10 Flag=10 fix FCAL ramp. The scaling factor were inverted. Ramp_M = Ramp_H*9.96 Ramp_L = Ramp_H*9.96*9.67

11 Flag=11 Update HEC Ramp

12 Flag=12 Fill the LArMinBiasMC and LArMinBiasAverageMC from ntuple created by LArMinBiasAlg (in CaloCondPhysAlgs)
13 Flag=13 Filling the OnlOflMap for SC

14 Flag=14 Fill the LArMinBiasAverageMC type in LArPileupAverage folder for BCID correction from ntuple created by LArMinBiasAlg (in CaloCondPhysAlgs)

15 Flag=15 Scale the MinBiasAverage by constant

16 Flag=16 Fix the EMEC-A FT 7 slot 11 OnOfIdMap cable swap

17 Flag=17 replace All shape

18 Flag=18 Fix the calib. line mapping for the EMF, where all channels in EMB-A FT 0 has CL=channel

Definition at line 35 of file FixLArElecCalib.h.

Member Typedef Documentation

◆ ROW_t

typedef std::pair<HWIdentifier,std::vector<float> > FixLArElecCalib::ROW_t
private

Definition at line 135 of file FixLArElecCalib.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

◆ VROW_t

typedef std::vector<ROW_t> FixLArElecCalib::VROW_t
private

Definition at line 136 of file FixLArElecCalib.h.

Constructor & Destructor Documentation

◆ FixLArElecCalib()

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

Definition at line 111 of file FixLArElecCalib.cxx.

111 :
112 AthAlgorithm(name,pSvcLocator),
113 m_em_idhelper(nullptr),
114 m_hec_idhelper(nullptr),
115 m_fcal_idhelper(nullptr),
116 m_online_idhelper(nullptr),
117 m_sem_idhelper(nullptr),
118 m_shec_idhelper(nullptr),
119 m_sonline_idhelper(nullptr),
120 m_scell_idhelper(nullptr)
121{
122
123}
AthAlgorithm()
Default constructor:
const LArOnline_SuperCellID * m_sonline_idhelper
const LArEM_SuperCell_ID * m_sem_idhelper
const LArHEC_ID * m_hec_idhelper
const LArOnlineID * m_online_idhelper
const LArFCAL_ID * m_fcal_idhelper
const CaloCell_SuperCell_ID * m_scell_idhelper
const LArEM_ID * m_em_idhelper
const LArHEC_SuperCell_ID * m_shec_idhelper

◆ ~FixLArElecCalib()

virtual FixLArElecCalib::~FixLArElecCalib ( )
virtualdefault

Member Function Documentation

◆ addMphysOverMcal()

StatusCode FixLArElecCalib::addMphysOverMcal ( const LArOnOffIdMapping * cabling)
private

Definition at line 1424 of file FixLArElecCalib.cxx.

1424 {
1425 // add the contain to DetStore.
1426 // fill 1 for all FCAL/HEC channels
1427
1428 ATH_MSG_INFO ( " addMphysOverMcal() " );
1429 // create LArMphysOverMcal with FCAL numbers = 1.
1430
1431 std::unique_ptr<LArMphysOverMcalMC> mphys = std::make_unique<LArMphysOverMcalMC>();
1432 mphys->setGroupingType(LArConditionsContainerBase::SingleGroup);
1433 ATH_CHECK( mphys->initialize() );
1434
1435 int n=0;
1436
1437 std::string filename("FCal_noise_minbias_adc2mev.txt");
1438
1439 std::ifstream infile(filename.c_str() ) ;
1440
1441 if(! infile.good() )
1442 {
1443 ATH_MSG_ERROR(" fail to open file "<<filename);
1444 return StatusCode::FAILURE;
1445 }
1446
1447
1448 char s[200];
1449
1450 ATH_MSG_INFO("Opened FCAL file");
1451 for(int i = 0;i<27;++i){
1452 infile.getline(s,200);
1453 ATH_MSG_INFO(s);
1454 }
1455
1456 std::string str_id;
1457
1458 float vol,noise_p,noise_h, noise_m,noise_l;
1459 float ramp_val,ua2MeV_val,dac2uA_val;
1460 std::string channel_type;
1461 while ( infile>>str_id >>vol>>noise_p >> noise_h>>noise_m>>noise_l
1462 >>channel_type>> ramp_val>>ua2MeV_val>>dac2uA_val )
1463 {
1464 const char* ch_id = str_id.c_str();
1465
1466 if(str_id.substr(0,1)!=std::string("A")){
1467 ATH_MSG_INFO(" skipping string"<<str_id);
1468 continue;
1469 }
1470
1471 int mod,phi,eta;
1472 sscanf(ch_id, "A%d.%d.%d", &mod,&phi,&eta);
1473 if(phi>7){
1474 ATH_MSG_INFO(" skipping phi"<<str_id<<" phi="<<phi);
1475 continue;
1476 }
1477
1478 ++n;
1479
1480 ATH_MSG_INFO(" Setting channel "<<str_id);
1481 Identifier id = m_fcal_idhelper->channel_id(2,mod,eta,phi);
1482 HWIdentifier hid = cabling->createSignalChannelID(id);
1483
1485 t.m_data =1. ;
1486 mphys->setPdata(hid,t,0);
1487 mphys->setPdata(hid,t,1);
1488 mphys->setPdata(hid,t,2);
1489 }
1490
1491 ATH_MSG_INFO(" MphysOverMcal added "<<n<<" FCAL channels");
1492
1493 bool EM=false ;
1494 bool withGain=false;
1495 int nvar = 1;
1496 std::string hec_filename("mphys_hec.txt");
1497 ATH_CHECK( ReadFile(hec_filename,cabling,EM,withGain,nvar) );
1498
1499 int igain=0;
1500 n=0;
1501 for( ROW_t &it : m_cache[igain])
1502 {
1503 HWIdentifier hid = it.first;
1505 t.m_data =1. ;
1506 mphys->setPdata(hid,t,0);
1507 mphys->setPdata(hid,t,1);
1508 mphys->setPdata(hid,t,2);
1509 ++n;
1510
1511 }
1512
1513 ATH_CHECK( detStore()->record(std::move(mphys),"LArMphysOverMcal") );
1514 ILArMphysOverMcal* imphys=nullptr;
1515 ATH_CHECK (detStore()->symLink(mphys.get(),imphys) );
1516
1517 ATH_MSG_ERROR(" Number of HEC channel added "<<n <<" per gain ");
1518 return StatusCode::SUCCESS ;
1519}
Scalar eta() const
pseudorapidity method
Scalar phi() const
phi method
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x)
#define ATH_MSG_INFO(x)
const ServiceHandle< StoreGateSvc > & detStore() const
StatusCode ReadFile(const std::string &filename, const LArOnOffIdMapping *cabling, bool EM, bool withGain, int nvar)
std::pair< HWIdentifier, std::vector< float > > ROW_t
str infile
Definition run.py:13

◆ ATLAS_NOT_THREAD_SAFE() [1/20]

StatusCode fix15 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( )

◆ ATLAS_NOT_THREAD_SAFE() [2/20]

StatusCode fix16 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( )

◆ ATLAS_NOT_THREAD_SAFE() [3/20]

StatusCode fixDACuAMeV FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( )
private

◆ ATLAS_NOT_THREAD_SAFE() [4/20]

virtual StatusCode stop FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( )
overridevirtual

◆ ATLAS_NOT_THREAD_SAFE() [5/20]

StatusCode fix3 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [6/20]

StatusCode fix4 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [7/20]

StatusCode fix5 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [8/20]

StatusCode fix6 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [9/20]

StatusCode fix7 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [10/20]

StatusCode fix8 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [11/20]

StatusCode fix9 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [12/20]

StatusCode fix10 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [13/20]

StatusCode fix11 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [14/20]

StatusCode fix17 FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)

◆ ATLAS_NOT_THREAD_SAFE() [15/20]

StatusCode updateHADfSampl FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const LArOnOffIdMapping * cabling)
private

◆ ATLAS_NOT_THREAD_SAFE() [16/20]

template<class T>
StatusCode update_EM_HEC FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const std::string & em_filename,
const std::string & hec_filename,
const LArOnOffIdMapping * cabling,
bool withGain,
int nvar )
private

◆ ATLAS_NOT_THREAD_SAFE() [17/20]

StatusCode updateEMfSampl FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const std::string & filename,
const LArOnOffIdMapping * cabling )
private

◆ ATLAS_NOT_THREAD_SAFE() [18/20]

StatusCode updateMinBias FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const std::string & filename,
const LArOnOffIdMapping * cabling )
private

◆ ATLAS_NOT_THREAD_SAFE() [19/20]

StatusCode updateEM_DACuAMeV FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const std::string & filename,
const LArOnOffIdMapping * cabling )
private

◆ ATLAS_NOT_THREAD_SAFE() [20/20]

template<class T>
StatusCode update_All FixLArElecCalib::ATLAS_NOT_THREAD_SAFE ( const std::string & filename,
const LArOnOffIdMapping * cabling,
bool withGain,
int nvar )
private

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< 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 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< 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< Algorithm > >::detStore ( ) const
inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< 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 FixLArElecCalib::execute ( )
overridevirtual

Definition at line 143 of file FixLArElecCalib.cxx.

143 {
144
145 const EventContext& ctx = Gaudi::Hive::currentContext();
146
147 SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey, ctx};
148 const LArOnOffIdMapping* cabling{*cablingHdl};
149 if(!cabling) {
150 ATH_MSG_WARNING( "Do not have cabling mapping from key " << m_cablingKey.key() );
151 }
152 if(m_fixFlag==13)
153 return fix13(cabling);
154
155 return StatusCode::SUCCESS;
156}
#define ATH_MSG_WARNING(x)
StatusCode fix13(const LArOnOffIdMapping *cabling)
SG::ReadCondHandleKey< LArOnOffIdMapping > m_cablingKey
IntegerProperty m_fixFlag

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< 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 & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 50 of file AthAlgorithm.cxx.

51{
52 // If we didn't find any symlinks to add, just return the collection
53 // from the base class. Otherwise, return the extended collection.
54 if (!m_extendedExtraObjects.empty()) {
56 }
57 return Algorithm::extraOutputDeps();
58}
DataObjIDColl m_extendedExtraObjects

◆ finalize()

virtual StatusCode FixLArElecCalib::finalize ( )
inlineoverridevirtual

Definition at line 43 of file FixLArElecCalib.h.

43{return StatusCode::SUCCESS;}

◆ fix1()

StatusCode FixLArElecCalib::fix1 ( const LArOnOffIdMapping * cabling)

Definition at line 219 of file FixLArElecCalib.cxx.

219 {
220
221 ATH_MSG_INFO ( " in fix1() " );
222
223 // Fix1 is for updating the EM DAC2uA, assuming symmetry.
224 // Input should be MC Conditions data with DetDescrVersion=ATLAS-DC3-05
225
226 // Pointer to StoreGate
227 const LArEM_ID* em_idhelper = nullptr;
228 ATH_CHECK( detStore()->retrieve(em_idhelper) );
229
230 const LArOnlineID* online_idhelper = nullptr;
231 ATH_CHECK( detStore()->retrieve(online_idhelper) );
232
233 // retrieve DetDescrManager
234 SmartIF<IToolSvc> toolSvc{Gaudi::svcLocator()->service("ToolSvc")};
235 ATH_CHECK(toolSvc.isValid());
236
237 LArDAC2uAMC * dac2ua = nullptr;
238 ATH_CHECK( detStore()->retrieve(dac2ua) );
239
240 // read in the file
241
242 std::ifstream infile("barrel_dac2ua_oct06.txt") ;
243 if(! infile.good() )
244 {
245 ATH_MSG_ERROR(" fail to open file ");
246 return StatusCode::FAILURE;
247 }
248
249
250 int det,samp,reg,eta;
251 float value ;
252 int n = 0;
253 while ( infile>>det>>samp>>reg>>eta>>value )
254 {
255 ++n;
256 ATH_MSG_DEBUG(" det,samp,reg,eta,value="
257 <<det<<" "
258 <<samp<<" "
259 <<reg<<" "
260 <<eta<<" "
261 <<value<<" " );
262
263 Identifier id ;
264 if ( det==1 && samp==1 && reg==0 && eta==0 ){
265 // eta=0 for strip, not connected, but keep it
266 // Removed disconnected channels: 03/2016 RDS
267 // ATH_MSG_DEBUG(" disconnected strip ");
268 // id = em_idhelper->disc_channel_id(det,samp,reg,eta,0);
269 }else
270 id = em_idhelper->channel_id(det,samp,reg,eta,0);
271
272 HWIdentifier hid = cabling->createSignalChannelID(id);
273 LArDAC2uAComplete::LArCondObj & t = dac2ua->get(hid,0);
274 ATH_MSG_DEBUG(" online id = "<<online_idhelper->print_to_string(hid));
275
276 if( t.isEmpty() )
277 {
278 ATH_MSG_WARNING(" No existing conditions data " );
279 continue ;
280 }
281 ATH_MSG_DEBUG(" Old DAC2uA = "<< t.m_data<< " " );
282
283 t.m_data = value ;
284
285 const LArDAC2uAComplete::LArCondObj & t3 = dac2ua->get(hid,0);
286 ATH_MSG_DEBUG(" New DAC2uA = "<< t3.m_data << " " );
287
288 }
289
290 ATH_MSG_DEBUG(" Number of entries changes = " <<n );
291 return StatusCode::SUCCESS;
292}
#define ATH_MSG_DEBUG(x)
std::string print_to_string(Identifier id, const IdContext *context=0) const
or provide the printout in string form
ConstReference get(const HWIdentifier id, unsigned int gain=0) const
get data with online identifier
Identifier channel_id(const ExpandedIdentifier &exp_id) const
Build a cell identifier from an expanded identifier.
std::vector< LUCID_RawData_p1 > t3
retrieve(aClass, aKey=None)
Definition PyKernel.py:110

◆ fix12()

StatusCode FixLArElecCalib::fix12 ( const LArOnOffIdMapping * cabling)

Definition at line 1997 of file FixLArElecCalib.cxx.

1997 {
1998
1999 ATH_MSG_INFO ( " in fix12() " );
2000
2001 // Fix12 is for filling the MinBiasRMS and MinBiasAverage from ntuple
2002
2003 // Create new container
2004 std::unique_ptr<LArMinBiasMC> minbias = std::make_unique<LArMinBiasMC>();
2005 ATH_CHECK( minbias->setGroupingType("Single", msg()) );
2006 ATH_CHECK( minbias->initialize() );
2007 std::unique_ptr<LArMinBiasAverageMC> minbias_av = std::make_unique<LArMinBiasAverageMC>();
2008 ATH_CHECK( minbias_av->setGroupingType("Single", msg()) );
2009 ATH_CHECK( minbias_av->initialize() );
2010 //
2011 std::unique_ptr<TFile> fin= std::make_unique<TFile>("ntuple_av.root");
2012 TTree *tin=dynamic_cast<TTree*>(fin->Get("m_tree"));
2013 if (!tin) {
2014 ATH_MSG_ERROR("Can't read tree m_tree");
2015 return StatusCode::FAILURE;
2016 }
2017 const int nindex=1833;
2018 int ncell;
2019 int identifier[nindex];
2020 int layer[nindex];
2021 int region[nindex];
2022 int ieta[nindex];
2023 float eta[nindex];
2024 double average[nindex];
2025 double rms[nindex];
2026 TBranch *b_ncell;
2027 TBranch *b_identifier;
2028 TBranch *b_layer;
2029 TBranch *b_region;
2030 TBranch *b_ieta;
2031 TBranch *b_eta;
2032 TBranch *b_average;
2033 TBranch *b_rms;
2034 tin->SetMakeClass(1);
2035 tin->SetBranchAddress("ncell", &ncell, &b_ncell);
2036 tin->SetBranchAddress("identifier", identifier, &b_identifier);
2037 tin->SetBranchAddress("layer", layer, &b_layer);
2038 tin->SetBranchAddress("region", region, &b_region);
2039 tin->SetBranchAddress("ieta", ieta, &b_ieta);
2040 tin->SetBranchAddress("eta", eta, &b_eta);
2041 tin->SetBranchAddress("average", average, &b_average);
2042 tin->SetBranchAddress("rms", rms, &b_rms);
2043 tin->GetEntry(0);
2044
2045 if(ncell>nindex) ncell=nindex;
2046 for(int icell=0; icell<ncell; ++icell) {
2047
2048 Identifier32 id32(identifier[icell]);
2049 Identifier id(id32);
2050 HWIdentifier hid = cabling->createSignalChannelID(id);
2051
2052 minbias->set(hid, rms[icell]);
2053 minbias_av->set(hid, average[icell]);
2054
2055 }
2056
2057 ATH_MSG_INFO(" number of channels in intuple ="<<ncell);
2058
2059 ATH_MSG_INFO ( "Stored container " << minbias->totalNumberOfConditions() << " conditions, key LArMinBias " );
2060 ATH_CHECK( detStore()->record(std::move(minbias),"LArMinBias") );
2061 //ATH_CHECK( detStore()->symLink(minbias.get(), dynamic_cast<ILArMinBias*>(minbias.get())) );
2063 ATH_MSG_INFO ( "Stored container " << minbias_av->totalNumberOfConditions() << " conditions, key LArMinBiasAverage " );
2064 ATH_CHECK( detStore()->record(std::move(minbias_av),"LArMinBiasAverage") );
2065 //ATH_CHECK( detStore()->symLink(minbias_av.get(), dynamic_cast<ILArMinBiasAverage*>(minbias_av.get())) );
2067
2068 return StatusCode::SUCCESS;
2069}
MsgStream & msg() const
@ layer
Definition HitInfo.h:79

◆ fix13()

StatusCode FixLArElecCalib::fix13 ( const LArOnOffIdMapping * cabling)

Definition at line 1726 of file FixLArElecCalib.cxx.

1726 {
1727
1728 ATH_MSG_INFO ( " in fix13() " );
1729 // update EM/HEC 7 data objects.
1730
1731 ToolHandle<ICaloSuperCellIDTool> scidtool("CaloSuperCellIDTool");
1732 if ( scidtool.retrieve().isFailure() ) {
1733 ATH_MSG_ERROR ( " Could not retrieve scitool " );
1734 }
1735
1736 std::set<Identifier> scidset;
1737 std::map<HWIdentifier,int> sslot_schannel_idx;
1738 std::vector<unsigned int> OrderOnlHash_OffID;
1739 std::vector<unsigned int> OrderOffHash_OnlID;
1740 OrderOnlHash_OffID.resize(40000,0);
1741 OrderOffHash_OnlID.resize(40000,0);
1742
1743 Identifier::size_type channel_size = m_em_idhelper->channel_hash_max();
1744 unsigned int min_onl_hash = 999999;
1745 unsigned int max_onl_hash = 0;
1746 unsigned int min_off_hash = 999999;
1747 unsigned int max_off_hash = 0;
1748 std::ofstream hashes("AllHashChannels_All.txt");
1749 std::ofstream of("AllChannels_EM.txt");
1750
1751 for(size_t i = 0; i < channel_size ; i++) {
1752 Identifier chid = m_em_idhelper->channel_id(i);
1753 bool isEM = ( m_em_idhelper->is_em_barrel( chid ) || m_em_idhelper->is_em_endcap( chid ) ) ;
1754 if ( ! isEM ) continue;
1755 Identifier SCID = scidtool->offlineToSuperCellID(chid);
1756 if ( scidset.find(SCID) == scidset.end() && (m_scell_idhelper->calo_cell_hash( SCID ).value() < 99999) )
1757 scidset.insert(SCID);
1758 else continue;
1759 HWIdentifier hwid = cabling->createSignalChannelID(chid);
1760 if ( i == 0 ){
1761 of << "Off ID\t\tSCID\t\tOnl ID\t\tFT\tslot\tB-E pos_neg\tSamp\teta\tphi\tFEB_ID\t\tSHWID\t" << std::endl;
1762 }
1763 of << chid.get_identifier32().get_compact() << "\t0x" <<
1764 SCID.get_identifier32().get_compact() << "\t0x" <<
1765 hwid.get_identifier32().get_compact() ;
1766 Identifier regId = m_em_idhelper->region_id( chid );
1767 int reg = m_em_idhelper->region( chid );
1768 IdentifierHash regHash = m_em_idhelper->region_hash( regId );
1769 int samp = m_em_idhelper->sampling( chid );
1770 int feedthrough = m_online_idhelper->feedthrough( hwid );
1771 int bar_ec = m_online_idhelper->barrel_ec( hwid );
1772 int pos_neg = m_online_idhelper->pos_neg( hwid );
1773 int slot = m_online_idhelper->slot( hwid );
1774 float etaG = m_em_idhelper->etaGranularity(regHash);
1775 float phiG = m_em_idhelper->phiGranularity(regHash);
1776 float eta0 = m_em_idhelper->eta0(regHash);
1777 float phi0 = m_em_idhelper->phi0(regHash);
1778 float eta = etaG*m_em_idhelper->eta( chid ) + eta0;
1779 float phi = phiG*m_em_idhelper->phi( chid ) + phi0;
1780 char etaChar[10];
1781 char phiChar[10];
1782 sprintf(etaChar,"%5.4f",eta);
1783 sprintf(phiChar,"%5.4f",phi);
1784 of << std::dec << "\t" << feedthrough << "\t" << slot;
1785 of << "\t" << bar_ec << " " << pos_neg << "\t\t" << samp << "\t";
1786 of << etaChar << " " << phiChar << "\t0x";
1787 HWIdentifier feb_id = m_online_idhelper->feb_Id(hwid);
1788 int sslot=1;
1789 bool InSlot= m_em_idhelper->is_em_endcap_inner ( chid );
1790 bool OutSlot= m_em_idhelper->is_em_endcap_outer ( chid ) ;
1791 bool secondOutSlot=false;
1792 if ( OutSlot ) {
1793 if ( (samp==1) && (reg==5) ) secondOutSlot=true;
1794 if ( (samp==2) && (reg==1) && ( m_em_idhelper->eta( chid )>=39 ) )
1795 secondOutSlot=true;
1796 if ( (samp==3) && (reg==0) && ( m_em_idhelper->eta( chid )>=18 ) )
1797 secondOutSlot=true;
1798 }
1799 HWIdentifier sfeb_id = m_sonline_idhelper->feb_Id(bar_ec,pos_neg,feedthrough,sslot);
1800
1801 if ( InSlot ) {
1802 sfeb_id = m_sonline_idhelper->feb_Id(bar_ec,pos_neg,feedthrough-1,sslot+1);
1803 }
1804 if ( secondOutSlot ) {
1805 sfeb_id = m_sonline_idhelper->feb_Id(bar_ec,pos_neg,feedthrough,sslot+1);
1806 }
1807 if ( sslot_schannel_idx.find(sfeb_id) == sslot_schannel_idx.end() )
1808 sslot_schannel_idx[sfeb_id]=0;
1809 else sslot_schannel_idx[sfeb_id]++;
1810 HWIdentifier shwid = m_sonline_idhelper->channel_Id(sfeb_id,sslot_schannel_idx[sfeb_id]);
1811 unsigned int onl_hash = m_sonline_idhelper->channel_Hash( shwid ).value();
1812 if ( onl_hash > max_onl_hash ) max_onl_hash = onl_hash;
1813 if ( onl_hash < min_onl_hash ) min_onl_hash = onl_hash;
1814 unsigned int off_hash = m_scell_idhelper->calo_cell_hash( SCID ).value();
1815 if ( off_hash > max_off_hash ) max_off_hash = off_hash;
1816 if ( off_hash < min_off_hash ) min_off_hash = off_hash;
1817 of << std::hex;
1818 of << feb_id.get_identifier32().get_compact();
1819 of << "\t0x" << shwid.get_identifier32().get_compact();
1820 of << " " << m_online_idhelper->channel_name(hwid) << " " << m_sonline_idhelper->channel_name(shwid);
1821 of << std::dec;
1822 of << " " << onl_hash << " " << off_hash;
1823 of << std::endl;
1824
1825 if ( m_scell_idhelper->calo_cell_hash( SCID ).value() < 40000 )
1826 OrderOffHash_OnlID[m_scell_idhelper->calo_cell_hash( SCID ).value()] = shwid.get_identifier32().get_compact();
1827 if ( m_sonline_idhelper->channel_Hash( shwid ).value() < 40000 )
1828 OrderOnlHash_OffID[m_sonline_idhelper->channel_Hash( shwid ).value()] = SCID.get_identifier32().get_compact();
1829 }
1830 of << "Hash check. Onl : " << min_onl_hash << " " << max_onl_hash;
1831 of << "; Off : " << min_off_hash << " " << max_off_hash << std::endl;;
1832 of.close();
1833
1834 channel_size = m_hec_idhelper->channel_hash_max();
1835 min_onl_hash = 999999;
1836 max_onl_hash = 0;
1837 min_off_hash = 999999;
1838 max_off_hash = 0;
1839 std::ofstream of1("AllChannels_HEC.txt");
1840
1841 of1 << "Off ID\t\tSCID\t\tOnl ID\t\tFT\tslot\tB-E pos_neg\tSamp\teta\tphi\tFEB_ID\t\tSHWID\t" << std::endl;
1842 for(size_t i = 0; i < channel_size ; i++) {
1843 Identifier chid = m_hec_idhelper->channel_id(i);
1844 if ( ! ( m_hec_idhelper->is_lar_hec( chid ) ) ) continue;
1845 Identifier SCID = scidtool->offlineToSuperCellID(chid);
1846 if ( scidset.find(SCID) == scidset.end() )
1847 scidset.insert(SCID);
1848 else continue;
1849 HWIdentifier hwid = cabling->createSignalChannelID(chid);
1850 of1 << chid.get_identifier32().get_compact() << "\t0x" <<
1851 SCID.get_identifier32().get_compact() << "\t0x" <<
1852 hwid.get_identifier32().get_compact() ;
1853 Identifier regId = m_hec_idhelper->region_id( chid );
1854 IdentifierHash regHash = m_hec_idhelper->region_hash( regId );
1855 int samp = m_hec_idhelper->sampling( chid );
1856 int feedthrough = m_online_idhelper->feedthrough( hwid );
1857 int bar_ec = m_online_idhelper->barrel_ec( hwid );
1858 int pos_neg = m_online_idhelper->pos_neg( hwid );
1859 int slot = m_online_idhelper->slot( hwid );
1860 float etaG = m_hec_idhelper->etaGranularity(regHash);
1861 float phiG = m_hec_idhelper->phiGranularity(regHash);
1862 float eta0 = m_hec_idhelper->eta0(regHash);
1863 float phi0 = m_hec_idhelper->phi0(regHash);
1864 float eta = etaG*m_hec_idhelper->eta( chid ) + eta0;
1865 float phi = phiG*m_hec_idhelper->phi( chid ) + phi0;
1866 char etaChar[10];
1867 char phiChar[10];
1868 sprintf(etaChar,"%5.4f",eta);
1869 sprintf(phiChar,"%5.4f",phi);
1870 of1 << std::dec << "\t" << feedthrough << "\t" << slot;
1871 of1 << "\t" << bar_ec << " " << pos_neg << "\t\t" << samp << "\t";
1872 of1 << etaChar << " " << phiChar << "\t0x";
1873 HWIdentifier feb_id = m_online_idhelper->feb_Id(hwid);
1874 int sslot=1;
1875 int subtract=0;
1876 HWIdentifier sfeb_id = m_sonline_idhelper->feb_Id(bar_ec,pos_neg,feedthrough,sslot);
1877 if ( sslot_schannel_idx.find(sfeb_id) == sslot_schannel_idx.end() )
1878 sslot_schannel_idx[sfeb_id]=0;
1879 else sslot_schannel_idx[sfeb_id]++;
1880 if ( (feedthrough==6) && (bar_ec==1) && (slot>9) ) {
1881 sfeb_id = m_sonline_idhelper->feb_Id(bar_ec,pos_neg,feedthrough,sslot+1);
1882 subtract=192;
1883 }
1884 HWIdentifier shwid = m_sonline_idhelper->channel_Id(sfeb_id,sslot_schannel_idx[sfeb_id]-subtract);
1885 unsigned int onl_hash = m_sonline_idhelper->channel_Hash( shwid ).value();
1886 if ( onl_hash > max_onl_hash ) max_onl_hash = onl_hash;
1887 if ( onl_hash < min_onl_hash ) min_onl_hash = onl_hash;
1888 unsigned int off_hash = m_scell_idhelper->calo_cell_hash( SCID ).value();
1889 if ( off_hash > max_off_hash ) max_off_hash = off_hash;
1890 if ( off_hash < min_off_hash ) min_off_hash = off_hash;
1891 of1 << std::hex;
1892 of1 << feb_id.get_identifier32().get_compact();
1893 of1 << "\t0x" << shwid.get_identifier32().get_compact();
1894 of1 << " " << m_online_idhelper->channel_name(hwid) << " " << m_sonline_idhelper->channel_name(shwid);
1895 of1 << std::dec;
1896 of1 << " " << onl_hash << " " << off_hash;
1897 of1 << std::endl;
1898
1899 OrderOffHash_OnlID[m_scell_idhelper->calo_cell_hash( SCID ).value()] = shwid.get_identifier32().get_compact();
1900 OrderOnlHash_OffID[m_sonline_idhelper->channel_Hash( shwid ).value()] = SCID.get_identifier32().get_compact();
1901 }
1902
1903 channel_size = m_fcal_idhelper->channel_hash_max();
1904 min_onl_hash = 999999;
1905 max_onl_hash = 0;
1906 min_off_hash = 999999;
1907 max_off_hash = 0;
1908 std::ofstream of3("AllChannels_FCAL.txt");
1909
1910 of3 << "Off ID\t\tSCID\t\tOnl ID\t\tFT\tslot\tB-E pos_neg\tSamp\teta\tphi\tFEB_ID\t\tSHWID\t" << std::endl;
1911 for(size_t i = 0; i < channel_size ; i++) {
1912 Identifier chid = m_fcal_idhelper->channel_id(i);
1913 if ( ! ( m_fcal_idhelper->is_lar_fcal( chid ) ) ) continue;
1914 Identifier SCID = scidtool->offlineToSuperCellID(chid);
1915 if ( scidset.find(SCID) == scidset.end() )
1916 scidset.insert(SCID);
1917 else continue;
1918 HWIdentifier hwid = cabling->createSignalChannelID(chid);
1919 of3 << chid.get_identifier32().get_compact() << "\t0x" <<
1920 SCID.get_identifier32().get_compact() << "\t0x" <<
1921 hwid.get_identifier32().get_compact() ;
1922 int samp = 0;
1923 int feedthrough = m_online_idhelper->feedthrough( hwid );
1924 int bar_ec = m_online_idhelper->barrel_ec( hwid );
1925 int pos_neg = m_online_idhelper->pos_neg( hwid );
1926 int slot = m_online_idhelper->slot( hwid );
1927 of3 << std::dec << "\t" << feedthrough << "\t" << slot;
1928 of3 << "\t" << bar_ec << " " << pos_neg << "\t\t" << samp << "\t";
1929 HWIdentifier feb_id = m_online_idhelper->feb_Id(hwid);
1930 int sslot=1;
1931 HWIdentifier sfeb_id = m_sonline_idhelper->feb_Id(bar_ec,pos_neg,feedthrough,sslot);
1932 if ( (feedthrough==6) && (bar_ec==1) && (slot>9) ) {
1933 sfeb_id = m_sonline_idhelper->feb_Id(bar_ec,pos_neg,feedthrough,sslot+1);
1934 }
1935 if ( sslot_schannel_idx.find(sfeb_id) == sslot_schannel_idx.end() )
1936 sslot_schannel_idx[sfeb_id]=0;
1937 else sslot_schannel_idx[sfeb_id]++;
1938 HWIdentifier shwid = m_sonline_idhelper->channel_Id(sfeb_id,sslot_schannel_idx[sfeb_id]);
1939 unsigned int onl_hash = m_sonline_idhelper->channel_Hash( shwid ).value();
1940 if ( onl_hash > max_onl_hash ) max_onl_hash = onl_hash;
1941 if ( onl_hash < min_onl_hash ) min_onl_hash = onl_hash;
1942 unsigned int off_hash = m_scell_idhelper->calo_cell_hash( SCID ).value();
1943 if ( off_hash > max_off_hash ) max_off_hash = off_hash;
1944 if ( off_hash < min_off_hash ) min_off_hash = off_hash;
1945 of3 << std::hex;
1946 of3 << feb_id.get_identifier32().get_compact();
1947 of3 << "\t0x" << shwid.get_identifier32().get_compact();
1948 of3 << " " << m_online_idhelper->channel_name(hwid) << " " << m_sonline_idhelper->channel_name(shwid);
1949 of3 << std::dec;
1950 of3 << " " << onl_hash << " " << off_hash;
1951 of3 << std::endl;
1952
1953 OrderOffHash_OnlID[m_scell_idhelper->calo_cell_hash( SCID ).value()] = shwid.get_identifier32().get_compact();
1954 OrderOnlHash_OffID[m_sonline_idhelper->channel_Hash( shwid ).value()] = SCID.get_identifier32().get_compact();
1955 }
1956
1957 of3 << "Hash check. Onl : " << min_onl_hash << " " << max_onl_hash;
1958 of3 << "; Off : " << min_off_hash << " " << max_off_hash << std::endl;;
1959 of3.close();
1960
1961
1962
1963 const uint32_t onlHashMax=m_sonline_idhelper->channelHashMax();
1964 coral::AttributeListSpecification *spec_onOff = new coral::AttributeListSpecification();
1965 spec_onOff->extend("OnlineHashToOfflineId", "blob");
1966 spec_onOff->extend<unsigned>("version");
1967 std::unique_ptr<AthenaAttributeList> al_onOff = std::make_unique<AthenaAttributeList>(*spec_onOff);
1968 coral::Blob& blobOnOff=(*al_onOff)["OnlineHashToOfflineId"].data<coral::Blob>();
1969 (*al_onOff)["version"].setValue(0U);
1970 blobOnOff.resize(onlHashMax*sizeof(uint32_t));
1971 uint32_t* pBlobOnOff=static_cast<uint32_t*>(blobOnOff.startingAddress());
1972 unsigned nConn=0;
1973
1974 const uint32_t emptyId=Identifier().get_identifier32().get_compact();
1975 for(size_t ii=0;ii<onlHashMax;ii++) {
1976 if (OrderOnlHash_OffID[ii] != 0) {
1977 pBlobOnOff[ii]=OrderOnlHash_OffID[ii];
1978 nConn++;
1979 }
1980 else
1981 pBlobOnOff[ii]=emptyId;
1982 }
1983
1984 ATH_MSG_INFO ( "HashMax=" << onlHashMax << ", connected=" << nConn );
1985
1986 ATH_CHECK( detStore()->record(std::move(al_onOff),"/LAR/IdentifierOfl/OnOffIdMap_SC") );
1987
1988 hashes << "idx \t Off2OnlId \t Onl2OffId" << std::endl;
1989 for(size_t ii=0;ii<40000;ii++)
1990 if ( (OrderOffHash_OnlID[ii] != 0) || ( OrderOnlHash_OffID[ii] != 0 ) )
1991 hashes << std::dec << ii << std::hex << " \t " << OrderOffHash_OnlID[ii] << " \t\t " << OrderOnlHash_OffID[ii] << std::endl;
1992 hashes.close();
1993
1994 return StatusCode::SUCCESS;
1995}
value_type get_compact() const
Get the compact id.
Identifier32 get_identifier32() const
Get the 32-bit version Identifier, will be invalid if >32 bits needed.
setEventNumber uint32_t

◆ fix14()

StatusCode FixLArElecCalib::fix14 ( const LArOnOffIdMapping * cabling)

Definition at line 2071 of file FixLArElecCalib.cxx.

2071 {
2072
2073 ATH_MSG_INFO ( " in fix14() " );
2074
2075 // Fix14 is for filling the LArPileupAverage from ntuple
2076
2077 const ILArfSampl* larfSampl = nullptr;
2078 ATH_CHECK( detStore()->retrieve(larfSampl,"") );
2079 // Will try to regenerate from scratch
2080 // Create new container
2081 std::unique_ptr<LArMinBiasAverageMC> minbias_av = std::make_unique<LArMinBiasAverageMC>();
2082 ATH_CHECK( minbias_av->setGroupingType("Single", msg()) );
2083 ATH_CHECK( minbias_av->initialize() );
2084 //
2085 std::unique_ptr<TFile> fin= std::make_unique<TFile>(m_infile.value().c_str());
2086 TTree *tin=dynamic_cast<TTree*>(fin->Get("m_tree"));
2087 if (!tin) {
2088 ATH_MSG_ERROR("Can't read tree m_tree");
2089 return StatusCode::FAILURE;
2090 }
2091 int ncell;
2092 std::vector<int> identifier(2862);
2093 std::vector<int> layer(2862);
2094 std::vector<int> region(2862);
2095 std::vector<int> ieta(2862);
2096 std::vector<float> eta(2862);
2097 std::vector<double> average(2862);
2098 std::vector<double> rms(2862);
2099 TBranch *b_ncell;
2100 TBranch *b_identifier;
2101 TBranch *b_layer;
2102 TBranch *b_region;
2103 TBranch *b_ieta;
2104 TBranch *b_eta;
2105 TBranch *b_average;
2106 TBranch *b_rms;
2107 tin->SetMakeClass(1);
2108 tin->SetBranchAddress("ncell", &ncell, &b_ncell);
2109 tin->SetBranchAddress("identifier", identifier.data(), &b_identifier);
2110 tin->SetBranchAddress("layer", layer.data(), &b_layer);
2111 tin->SetBranchAddress("region", region.data(), &b_region);
2112 tin->SetBranchAddress("ieta", ieta.data(), &b_ieta);
2113 tin->SetBranchAddress("eta", eta.data(), &b_eta);
2114 tin->SetBranchAddress("average", average.data(), &b_average);
2115 tin->SetBranchAddress("rms", rms.data(), &b_rms);
2116 tin->GetEntry(0);
2117
2118
2119 for(int icell=0; icell<ncell; ++icell) {
2120
2121 Identifier32 id32(identifier[icell]);
2122 Identifier id(id32);
2123 HWIdentifier hid = cabling->createSignalChannelID(id);
2124 float fsampl = larfSampl->FSAMPL(hid);
2125 minbias_av->set(hid, 6.31*average[icell]/fsampl);
2126
2127 }
2128
2129 ATH_MSG_INFO(" number of channels in intuple ="<<ncell);
2130
2131 ATH_MSG_INFO ( "Stored container " << minbias_av->totalNumberOfConditions() << " conditions, key LArPileupAverage " );
2132 ATH_CHECK( detStore()->record(std::move(minbias_av),"LArPileupAverage") );
2133 ATH_CHECK( detStore()->symLink(minbias_av.get(), dynamic_cast<ILArMinBiasAverage*>(minbias_av.get())) );
2134
2135 return StatusCode::SUCCESS;
2136}
StringProperty m_infile
virtual const float & FSAMPL(const HWIdentifier &id) const =0

◆ fix18()

StatusCode FixLArElecCalib::fix18 ( const LArCalibLineMapping * clmap)

Definition at line 2235 of file FixLArElecCalib.cxx.

2235 {
2236
2237 ATH_MSG_INFO ( " in fix18() " );
2238
2239 const uint32_t onlHashMax=m_online_idhelper->channelHashMax();
2240
2241 coral::AttributeListSpecification* spec_calib = new coral::AttributeListSpecification();
2242 spec_calib->extend("OnlineHashToCalibIds", "blob");
2243 spec_calib->extend<unsigned>("version");
2244 std::unique_ptr<AthenaAttributeList> al_calib = std::make_unique<AthenaAttributeList>(*spec_calib);
2245 coral::Blob& blobCalib=(*al_calib)["OnlineHashToCalibIds"].data<coral::Blob>();
2246 (*al_calib)["version"].setValue(0U);
2247 blobCalib.resize(onlHashMax*sizeof(uint32_t)*4); //Bigger than necessary
2248
2249 spec_calib->release();
2250 // cppcheck-suppress memleak
2251 spec_calib = nullptr;
2252
2253 uint32_t* pBlobCalib=static_cast<uint32_t*>(blobCalib.startingAddress());
2254 size_t calibIndex=0;
2255
2256
2257 // loop over all online hashes
2258 for (unsigned i=0;i<onlHashMax;++i) {
2259 const HWIdentifier hwid = m_online_idhelper->channel_Id(IdentifierHash(i));
2260
2261 const std::vector<HWIdentifier>& calibIDs=cabling->calibSlotLine(hwid);
2262 const size_t nCalibLines=calibIDs.size();
2263
2264 if (m_online_idhelper->isEMBchannel(hwid) && m_online_idhelper->pos_neg(hwid) == 1 && m_online_idhelper->feedthrough(hwid) == 0 ) { // EMB-A FT 0
2265 if(nCalibLines > 1) { //should not be in our channels
2266 ATH_MSG_ERROR("Wrong number of CL in our channels");
2267 return StatusCode::FAILURE;
2268 }
2269 if(nCalibLines==0) {
2270 pBlobCalib[calibIndex++]=0;
2271 } else {
2272 const HWIdentifier cmodule = m_online_idhelper->calib_module_Id(calibIDs[0]);
2273 const HWIdentifier newcl = m_online_idhelper->calib_channel_Id(cmodule, m_online_idhelper->channel(hwid) );
2274 if(!m_online_idhelper->isCalibration(newcl)) {
2275 ATH_MSG_ERROR("Something wrong, new CL not a calibration");
2276 ATH_MSG_ERROR("Channel: "<<hwid.get_identifier32().get_compact()<<" old CL: "<< calibIDs[0].get_identifier32().get_compact()<<" new CL: "<<newcl.get_identifier32().get_compact());
2277 return StatusCode::FAILURE;
2278 }
2279 ATH_MSG_INFO("Channel: "<<hwid.get_identifier32().get_compact()<<" old CL: "<< calibIDs[0].get_identifier32().get_compact()<<" new CL: "<<newcl.get_identifier32().get_compact());
2280 pBlobCalib[calibIndex++]=1;
2281 pBlobCalib[calibIndex++]=newcl.get_identifier32().get_compact();
2282 }
2283 } else {
2284 pBlobCalib[calibIndex++]=nCalibLines;
2285 for(uint32_t iCalib=0;iCalib<nCalibLines;++iCalib)
2286 pBlobCalib[calibIndex++]=calibIDs[iCalib].get_identifier32().get_compact();
2287 }
2288 }
2289 blobCalib.resize(calibIndex*sizeof(uint32_t)); //Size down to actual size
2290 ATH_MSG_INFO( "BlobSize CalibId:" << calibIndex);
2291
2292 // and now record the blob to SG
2293 ATH_CHECK(detStore()->record(std::move(al_calib),"/LAR/Identifier/CalibIdMap_EMF"));
2294
2295 return StatusCode::SUCCESS;
2296}

◆ fix2()

StatusCode FixLArElecCalib::fix2 ( const LArOnOffIdMapping * cabling)

Definition at line 295 of file FixLArElecCalib.cxx.

295 {
296
297 ATH_MSG_INFO ( " in fix2() " );
298
299 // Fix2 is for updating the FCAL conditions data after IdFix7
300 // Input should be MC Conditions data with DetDescrVersion=ATLAS-DC3-05
301
302 const LArFCAL_ID* fcal_idhelper = nullptr;
303 ATH_CHECK( detStore()->retrieve(fcal_idhelper) );
304
305 const LArOnlineID* online_idhelper = nullptr;
306 ATH_CHECK( detStore()->retrieve(online_idhelper) );
307
308 // retrieve DetDescrManager
309 SmartIF<IToolSvc> toolSvc{Gaudi::svcLocator()->service("ToolSvc")};
310 ATH_CHECK(toolSvc.isValid());
311
312 LArDAC2uAMC * dac2ua = nullptr;
313 ATH_CHECK( detStore()->retrieve(dac2ua) );
314
315 LAruA2MeVMC * ua2mev = nullptr;
316 ATH_CHECK( detStore()->retrieve(ua2mev) );
317
318 LArRampMC * ramp = nullptr;
319 ATH_CHECK( detStore()->retrieve(ramp) );
320
321 LArShape32MC * shape = nullptr;
322 ATH_CHECK( detStore()->retrieve(shape) );
323
324 LArNoiseMC * noise = nullptr;
325 ATH_CHECK( detStore()->retrieve(noise) );
326
327 LArfSamplMC * fsampl = nullptr;
328 ATH_CHECK( detStore()->retrieve(fsampl) );
329
330 LArMinBiasMC * minbias = nullptr;
331 ATH_CHECK( detStore()->retrieve(minbias) );
332
333 LArAutoCorrMC * ac = nullptr;
334 ATH_CHECK( detStore()->retrieve(ac) );
335
336
337 for( unsigned int gain=0;gain<3;++gain)
338 {
339 // loop over all the old data
341 LArRampMC::ConditionsMapIterator it_e = ramp->end(gain);
342
343 int n=0;
344 for (; it!=it_e;++it)
345 {
346 HWIdentifier hid = it.channelId();
347 if( (*it).isEmpty() ) continue;
348 if(cabling->isOnlineConnected(hid)){
349 continue;
350 }else
351 {
352 ATH_MSG_DEBUG(" unconnected channel"
353 << online_idhelper->print_to_string(hid) );
354 ac->get(hid,gain) = LArAutoCorrComplete::LArCondObj();
355 ramp->get(hid,gain) = LArRampComplete::LArCondObj();
356 shape->get(hid,gain) = LArShape32MC::LArCondObj();
357 noise->get(hid,gain) = LArNoiseComplete::LArCondObj();
358 if(gain==0){
359 dac2ua->get(hid,gain) = LArDAC2uAComplete::LArCondObj();
360 ua2mev->get(hid,gain) = LAruA2MeVComplete::LArCondObj();
361 fsampl->get(hid,gain) = LArfSamplComplete::LArCondObj();
362 minbias->get(hid,gain) = LArMinBiasComplete::LArCondObj();
363 }
364 ++n ;
365 }
366
367 }
368 ATH_MSG_DEBUG(" Gain="<<gain<<
369 " Number of entries removed = " <<n );
370 }
371
372
373 HWIdentifier hid_mod[3];
374 hid_mod[0] = online_idhelper->channel_Id(1,1,6,1,33);
375 hid_mod[1] = online_idhelper->channel_Id(1,1,6,10,12);
376 hid_mod[2] = online_idhelper->channel_Id(1,1,6,14,12);
377
378 // loop over all FCAL channels
379 std::vector<Identifier>::const_iterator it= fcal_idhelper->fcal_begin();
380 std::vector<Identifier>::const_iterator it_e= fcal_idhelper->fcal_end();
381 int n = 0;
382 for(;it!=it_e; ++it)
383 {
384 Identifier id = *it;
385 int pn = fcal_idhelper->pos_neg(id);
386 if(pn!=2) continue;
387 int phi = fcal_idhelper->phi(id);
388 if (phi>7) continue;
389
390 int module = fcal_idhelper->module(id)-1;
391
392 HWIdentifier hid = cabling->createSignalChannelID(id);
393
394 for(unsigned int gain = 0;gain<3;++gain){
395
396 LArRampComplete::LArCondObj& rampP = ramp->get(hid,gain) ;
397 if (!rampP.isEmpty()) continue ;
398
399 ATH_MSG_DEBUG(" channel needs repair "
400 << online_idhelper->print_to_string(hid) );
401 ATH_MSG_DEBUG ( fcal_idhelper->print_to_string(id) );
402 ATH_MSG_DEBUG ( "module = "<<module );
403
404 ++n;
405
406 LArRampComplete::LArCondObj ramp_std = ramp->get(hid_mod[module],gain) ;
407 ramp->setPdata(hid,ramp_std,gain);
408
409 LArAutoCorrComplete::LArCondObj ac_std = ac->get(hid_mod[module],gain) ;
410 ac->setPdata(hid,ac_std,gain);
411
412 LArShape32MC::LArCondObj shape_std = shape->get(hid_mod[module],gain) ;
413 shape->setPdata(hid,shape_std,gain);
414
415 LArNoiseComplete::LArCondObj noise_std = noise->get(hid_mod[module],gain) ;
416 noise->setPdata(hid,noise_std,gain);
417
418 if(gain==0){
419
420 LAruA2MeVComplete::LArCondObj ua2mev_std = ua2mev->get(hid_mod[module],gain) ;
421 ua2mev->setPdata(hid,ua2mev_std,gain);
422
423 LArDAC2uAComplete::LArCondObj dac2ua_std = dac2ua->get(hid_mod[module],gain) ;
424 dac2ua->setPdata(hid,dac2ua_std,gain);
425
426 LArfSamplComplete::LArCondObj fsampl_std = fsampl->get(hid_mod[module],gain) ;
427 fsampl->setPdata(hid,fsampl_std,gain);
428
429 LArMinBiasComplete::LArCondObj minbias_std = minbias->get(hid_mod[module],gain) ;
430 minbias->setPdata(hid,minbias_std,gain);
431
432 }
433
434 }
435 }
436
437 ATH_MSG_DEBUG( " Number of entries fixed " <<n );
438 return StatusCode::SUCCESS;
439}
void setPdata(const HWIdentifier id, const T &payload, unsigned int gain=0)
put payload in persistent data
ConstConditionsMapIterator begin(unsigned int gain) const
get iterator for all channels for a gain
ConstConditionsMapIterator end(unsigned int gain) const
end of all channels for this gain
id_iterator fcal_end() const
end iterator over full set of Fcal Identifiers for channels
id_iterator fcal_begin() const
begin iterator over full set of Fcal Identifiers for channels
int phi(const Identifier id) const
phi [0,15]
int pos_neg(const Identifier id) const
pos_neg : +/- 2 (A/C side)
HWIdentifier channel_Id(int barrel_ec, int pos_neg, int feedthrough, int slot, int channel) const
create channel identifier from fields
bool isEmpty() const
Definition LArRampP1.h:29

◆ initialize()

StatusCode FixLArElecCalib::initialize ( )
overridevirtual

Definition at line 125 of file FixLArElecCalib.cxx.

125 {
126 ATH_MSG_INFO ( " in initialize " );
127
136
137 ATH_CHECK( m_cablingKey.initialize() );
138 ATH_CHECK( m_CLKey.initialize() );
139
140 return StatusCode::SUCCESS;
141}
SG::ReadCondHandleKey< LArCalibLineMapping > m_CLKey

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< 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.

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

bool AthCommonMsg< 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< 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.

◆ print_object() [1/5]

void FixLArElecCalib::print_object ( const std::string & msg,
const LArAutoCorrMC::LArCondObj & obj )
private

Definition at line 1301 of file FixLArElecCalib.cxx.

1302{
1303 if( obj.isEmpty()){
1304 ATH_MSG_INFO(" LArAutoCorrMC " << s << " is empty" );
1305 }else
1306 {
1307 msg()<<MSG::DEBUG<<" LArAutoCorrMC" << s ;
1308 for (unsigned int i =0 ; i<obj.m_vAutoCorr.size();++i)
1309 msg()<< " " << obj.m_vAutoCorr[i];
1310 msg()<<endmsg;
1311 }
1312 return ;
1313}
#define endmsg

◆ print_object() [2/5]

void FixLArElecCalib::print_object ( const std::string & msg,
const LArNoiseMC::LArCondObj & obj )
private

Definition at line 1281 of file FixLArElecCalib.cxx.

1282{
1283 if( obj.isEmpty()){
1284 ATH_MSG_INFO(" LArNoiseMC " << msg << " is empty" );
1285 }else
1286 ATH_MSG_DEBUG(" LArNoiseMC" << msg << obj.m_Noise );
1287 return ;
1288}

◆ print_object() [3/5]

void FixLArElecCalib::print_object ( const std::string & msg,
const LArRampMC::LArCondObj & obj )
private

Definition at line 1328 of file FixLArElecCalib.cxx.

1329{
1330 if( obj.isEmpty()){
1331 ATH_MSG_INFO(" LArRampMC " << s << " is empty" );
1332 }else
1333 {
1334 msg()<<MSG::DEBUG<<" LArRampMC" << s ;
1335 for (unsigned int i =0 ; i<obj.m_vRamp.size();++i)
1336 msg()<< " " << obj.m_vRamp[i];
1337 msg()<<endmsg;
1338 }
1339 return ;
1340}

◆ print_object() [4/5]

void FixLArElecCalib::print_object ( const std::string & msg,
const LArShape32MC::LArCondObj & obj )
private

Definition at line 1357 of file FixLArElecCalib.cxx.

1358{
1359 if( obj.isEmpty()){
1360 ATH_MSG_INFO(" LArShape32MC " << s << " is empty" );
1361 }else
1362 {
1363 msg()<<MSG::DEBUG<<" LArShape" << s ;
1364 for (unsigned int i =0 ; i<obj.m_vShape.size();++i)
1365 msg()<< " " << obj.m_vShape[i];
1366 msg()<<endmsg;
1367
1368 msg()<<MSG::DEBUG<<" LArShapeDer" << s ;
1369 for (unsigned int i =0 ; i<obj.m_vShapeDer.size();++i)
1370 msg()<< " " << obj.m_vShapeDer[i];
1371 msg()<<endmsg;
1372 }
1373 return ;
1374}

◆ print_object() [5/5]

void FixLArElecCalib::print_object ( const std::string & msg,
const LArSingleFloatP & obj )
private

Definition at line 1401 of file FixLArElecCalib.cxx.

1402{
1403 if( obj.isEmpty()){
1404 ATH_MSG_INFO(" LArSingleFloatP " << msg << " is empty" );
1405 }else
1406 ATH_MSG_DEBUG(" LArSingleFloatP " << msg << obj.m_data );
1407 return ;
1408}

◆ ReadFile()

StatusCode FixLArElecCalib::ReadFile ( const std::string & filename,
const LArOnOffIdMapping * cabling,
bool EM,
bool withGain,
int nvar )
private

Definition at line 1142 of file FixLArElecCalib.cxx.

1142 {
1143
1144 std::ifstream infile(filename.c_str() ) ;
1145
1146 if(! infile.good() )
1147 {
1148 ATH_MSG_ERROR(" fail to open file "<<filename );
1149 return StatusCode::FAILURE;
1150 }
1151
1152 ATH_MSG_INFO(" Opened file "<<filename );
1153
1154 m_cache[0].clear();
1155 m_cache[1].clear();
1156 m_cache[2].clear();
1157
1158 char s[200];
1159 infile.getline(s,200);
1160 ATH_MSG_INFO(" first line of the file "<<s);
1161
1162 int det=2;
1163 int samp,reg,eta;
1164 int first;
1165 while ( infile>>first )
1166 {
1167 if (EM) {
1168 det = first;
1169 infile>>samp>>reg>>eta;
1170 }else
1171 {
1172 samp=first;
1173 infile>>reg>>eta;
1174 }
1175
1176 int gain=0;
1177 if(withGain) infile>>gain ;
1178 std::vector<float> vfl;
1179 float x;
1180
1181 for (int i=0;i<nvar;++i)
1182 {
1183 infile>>x;
1184 vfl.push_back(x);
1185 }
1186 ATH_MSG_DEBUG(" det,samp,reg,eta,values="
1187 <<det<<" "
1188 <<samp<<" "
1189 <<reg<<" "
1190 <<eta<<" " );
1191
1192 Identifier id ;
1193 if (EM)
1194 {
1195 if ( det==1 && samp==1 && reg==0 && eta==0 ){
1196 // eta=0 for strip, not connected, but keep it
1197 // Removed disconnected channels: 03/2016 RDS
1198 ATH_MSG_DEBUG(" disconnected strip ");
1199 }else
1200 id = m_em_idhelper->channel_id(det,samp,reg,eta,0);
1201 }else //HEC
1202 {
1203 id = m_hec_idhelper->channel_id(2,samp,reg,eta,0);
1204 }
1205 HWIdentifier hid = cabling->createSignalChannelID(id);
1206
1207 if (gain >= 0 && gain < 3)
1208 m_cache[gain].push_back(ROW_t(hid,vfl));
1209 }
1210
1211 return StatusCode::SUCCESS;
1212}
#define x
bool first
Definition DeMoScan.py:534

◆ ReadFileAll()

StatusCode FixLArElecCalib::ReadFileAll ( const std::string & filename,
const LArOnOffIdMapping * cabling,
bool withGain,
int nvar )
private

Definition at line 1214 of file FixLArElecCalib.cxx.

1214 {
1215
1216 std::ifstream infile(filename.c_str() ) ;
1217
1218 if(! infile.good() )
1219 {
1220 ATH_MSG_ERROR(" fail to open file "<<filename );
1221 return StatusCode::FAILURE;
1222 }
1223
1224 ATH_MSG_INFO(" Opened file "<<filename );
1225
1226 m_cache[0].clear();
1227 m_cache[1].clear();
1228 m_cache[2].clear();
1229
1230 char s[200];
1231 infile.getline(s,200);
1232 ATH_MSG_INFO(" first line of the file "<<s);
1233
1234 int det=0,samp=0,reg=0,eta=0;
1235 int first=0;
1236 while ( infile>>first>>samp>>reg>>eta ) {
1237 Identifier id ;
1238 int gain=0;
1239 float x;
1240 if(first <= 2) { // Calo 0,1,2 is EM
1241 det = first+1;
1242 if ( det==1 && samp==1 && reg==0 && eta==0 ){
1243 // eta=0 for strip, not connected, but keep it
1244 // Removed disconnected channels: 03/2016 RDS
1245 ATH_MSG_DEBUG(" disconnected strip ");
1246 // id = m_em_idhelper->disc_channel_id(det,samp,reg,eta,0);
1247 }else
1248 id = m_em_idhelper->channel_id(det,samp,reg,eta,0);
1249 } else if (first == 3) { // HEC
1250 id = m_hec_idhelper->channel_id(2,samp,reg,eta,0);
1251 } else if (first == 4) { // FCAL
1252 id = m_fcal_idhelper->channel_id(2,samp+1,eta,0);
1253 } else {
1254 ATH_MSG_WARNING(" Wrong calo number "<<first<<" skipping !!!");
1255 if(withGain) infile>>gain ;
1256 for (int i=0;i<nvar;++i) infile>>x;
1257 continue;
1258 }
1259 if(withGain) infile>>gain ;
1260 std::vector<float> vfl;
1261
1262 for (int i=0;i<nvar;++i)
1263 {
1264 infile>>x;
1265 vfl.push_back(x);
1266 }
1267 ATH_MSG_DEBUG(" first, det,samp,reg,eta,values="
1268 <<first<<" " <<det<<" " <<samp<<" " <<reg<<" " <<eta<<" " );
1269
1270 HWIdentifier hid = cabling->createSignalChannelID(id);
1271
1272 if (gain >= 0 && gain < 3) m_cache[gain].push_back(ROW_t(hid,vfl));
1273 else {
1274 ATH_MSG_WARNING(" Wrong gain " <<gain<<" skipping !");
1275 }
1276 }// while
1277
1278 return StatusCode::SUCCESS;
1279}

◆ 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< 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< Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ set_object() [1/5]

void FixLArElecCalib::set_object ( LArAutoCorrMC::LArCondObj & obj,
const std::vector< float > & v )
private

Definition at line 1315 of file FixLArElecCalib.cxx.

1316{
1317 if (v.size()<4 )
1318 {
1319 ATH_MSG_ERROR( "for LArAutoCorrMC, size of vector = "<< v.size()<<" is less then 4" );
1320 return ;
1321 }
1322
1323 obj.m_vAutoCorr.resize(4);
1324 std::copy_n(v.begin(), 4, obj.m_vAutoCorr.begin());
1325 return;
1326}

◆ set_object() [2/5]

void FixLArElecCalib::set_object ( LArNoiseMC::LArCondObj & obj,
const std::vector< float > & v )
private

Definition at line 1290 of file FixLArElecCalib.cxx.

1291{
1292 if (v.size()!=1)
1293 {
1294 ATH_MSG_ERROR ("for LArNoiseMC, size of vector = "<< v.size() <<" differs from one");
1295 return ;
1296 }
1297 obj.m_Noise = v[0];
1298 return;
1299}

◆ set_object() [3/5]

void FixLArElecCalib::set_object ( LArRampMC::LArCondObj & obj,
const std::vector< float > & v )
private

Definition at line 1342 of file FixLArElecCalib.cxx.

1343{
1344 if (v.size()!=1 )
1345 {
1346 ATH_MSG_ERROR ("for LArRampMC, size of vector = "<< v.size());
1347 return ;
1348 }
1349 obj.m_vRamp.resize(3);
1350 obj.m_vRamp[0]=0;
1351 obj.m_vRamp[1]=v[0];
1352 obj.m_vRamp[2]=0;
1353
1354 return;
1355}

◆ set_object() [4/5]

void FixLArElecCalib::set_object ( LArShape32MC::LArCondObj & obj,
const std::vector< float > & v )
private

Definition at line 1376 of file FixLArElecCalib.cxx.

1377{
1378
1379 if (v.size()!=64 && v.size()!=32 )
1380 {
1381 ATH_MSG_ERROR ("for LArShape32MC, size of vector = "<< v.size());
1382 return ;
1383 }
1384
1385 unsigned int size = 32;
1386 obj.m_vShape.resize(size);
1387 obj.m_vShapeDer.resize(size);
1388
1389 for(unsigned int i=0;i<size;++i){
1390 obj.m_vShape[i]=v[i];
1391 if(v.size() == 64 ) obj.m_vShapeDer[i]=v[i+size]; else obj.m_vShapeDer[i]=0.;
1392 }
1393
1394 return;
1395}

◆ set_object() [5/5]

void FixLArElecCalib::set_object ( LArSingleFloatP & obj,
const std::vector< float > & v )
private

Definition at line 1410 of file FixLArElecCalib.cxx.

1411{
1412 if (v.size()!=1)
1413 {
1414 ATH_MSG_ERROR( "for LArSingleFloatP, size of vector = "<< v.size() );
1415 return ;
1416 }
1417 obj.m_data = v[0];
1418 return;
1419}

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, and PyAthena::Alg.

Definition at line 66 of file AthAlgorithm.cxx.

66 {
68
69 if (sc.isFailure()) {
70 return sc;
71 }
72 ServiceHandle<ICondSvc> cs("CondSvc",name());
73 for (auto h : outputHandles()) {
74 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
75 // do this inside the loop so we don't create the CondSvc until needed
76 if ( cs.retrieve().isFailure() ) {
77 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
78 return StatusCode::SUCCESS;
79 }
80 if (cs->regHandle(this,*h).isFailure()) {
81 sc = StatusCode::FAILURE;
82 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
83 << " with CondSvc");
84 }
85 }
86 }
87 return sc;
88}
static Double_t sc
virtual StatusCode sysInitialize() override
Override sysInitialize.
AthCommonDataStore(const std::string &name, T... args)
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
::StatusCode StatusCode
StatusCode definition for legacy code.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< 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.

◆ update_DAC2uA()

StatusCode FixLArElecCalib::update_DAC2uA ( const std::string & em_filename,
const std::string & hec_filename )
private

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< 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 }
std::vector< SG::VarHandleKeyArray * > m_vhka

Member Data Documentation

◆ m_cablingKey

SG::ReadCondHandleKey<LArOnOffIdMapping> FixLArElecCalib::m_cablingKey {this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"}
private

Definition at line 67 of file FixLArElecCalib.h.

67{this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"};

◆ m_cache

VROW_t FixLArElecCalib::m_cache[3]
private

Definition at line 138 of file FixLArElecCalib.h.

◆ m_CLKey

SG::ReadCondHandleKey<LArCalibLineMapping> FixLArElecCalib::m_CLKey {this, "CalibLineKey", "LArCalibLineMap", "calib line key"}
private

Definition at line 68 of file FixLArElecCalib.h.

68{this, "CalibLineKey", "LArCalibLineMap", "calib line key"};

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_em_idhelper

const LArEM_ID* FixLArElecCalib::m_em_idhelper
private

Definition at line 123 of file FixLArElecCalib.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_fcal_idhelper

const LArFCAL_ID* FixLArElecCalib::m_fcal_idhelper
private

Definition at line 125 of file FixLArElecCalib.h.

◆ m_fixFlag

IntegerProperty FixLArElecCalib::m_fixFlag {this, "FixFlag", 1, "which fix to apply"}
private

Definition at line 70 of file FixLArElecCalib.h.

70{this, "FixFlag", 1, "which fix to apply"};

◆ m_g4Phys

StringProperty FixLArElecCalib::m_g4Phys {this, "G4Phys", "", "which physics list to use"}
private

Definition at line 71 of file FixLArElecCalib.h.

71{this, "G4Phys", "", "which physics list to use"};

◆ m_hec_idhelper

const LArHEC_ID* FixLArElecCalib::m_hec_idhelper
private

Definition at line 124 of file FixLArElecCalib.h.

◆ m_infile

StringProperty FixLArElecCalib::m_infile {this, "InputFile", "", "name of inpput file"}
private

Definition at line 72 of file FixLArElecCalib.h.

72{this, "InputFile", "", "name of inpput file"};

◆ m_online_idhelper

const LArOnlineID* FixLArElecCalib::m_online_idhelper
private

Definition at line 126 of file FixLArElecCalib.h.

◆ m_scell_idhelper

const CaloCell_SuperCell_ID* FixLArElecCalib::m_scell_idhelper
private

Definition at line 130 of file FixLArElecCalib.h.

◆ m_sem_idhelper

const LArEM_SuperCell_ID* FixLArElecCalib::m_sem_idhelper
private

Definition at line 127 of file FixLArElecCalib.h.

◆ m_shec_idhelper

const LArHEC_SuperCell_ID* FixLArElecCalib::m_shec_idhelper
private

Definition at line 128 of file FixLArElecCalib.h.

◆ m_sonline_idhelper

const LArOnline_SuperCellID* FixLArElecCalib::m_sonline_idhelper
private

Definition at line 129 of file FixLArElecCalib.h.

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


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