FixLArElecCalib.cxx File Reference

#include "FixLArElecCalib.h"
#include "Gaudi/Property.h"
#include "GaudiKernel/IToolSvc.h"
#include "CaloIdentifier/CaloCell_SuperCell_ID.h"
#include "CaloIdentifier/LArEM_ID.h"
#include "LArIdentifier/LArOnlineID.h"
#include "LArIdentifier/LArOnline_SuperCellID.h"
#include "StoreGate/StoreGateSvc.h"
#include "LArRawConditions/LArDAC2uAMC.h"
#include "LArRawConditions/LArShape32MC.h"
#include "LArRawConditions/LAruA2MeVMC.h"
#include "LArRawConditions/LArfSamplMC.h"
#include "LArRawConditions/LArMinBiasMC.h"
#include "LArRawConditions/LArMinBiasAverageMC.h"
#include "LArRawConditions/LArNoiseMC.h"
#include "LArRawConditions/LArRampMC.h"
#include "LArRawConditions/LArAutoCorrMC.h"
#include "CaloDetDescr/ICaloSuperCellIDTool.h"
#include "TTree.h"
#include "TFile.h"
#include "AthenaPoolUtilities/AthenaAttributeList.h"
#include "CoralBase/Blob.h"
#include <fstream>
#include <cstdint>

Go to the source code of this file.

Functions
StatusCode FixLArElecCalib::stop	ATLAS_NOT_THREAD_SAFE ()
	Install fatal handler with default options.
StatusCode FixLArElecCalib::fix3	ATLAS_NOT_THREAD_SAFE (const LArOnOffIdMapping *cabling)
StatusCode FixLArElecCalib::updateEMfSampl	ATLAS_NOT_THREAD_SAFE (const std::string &filename, const LArOnOffIdMapping *cabling)

Function Documentation

ATLAS_NOT_THREAD_SAFE() [1/3]

StatusCode FixLArElecCalib::stop ATLAS_NOT_THREAD_SAFE ( )

inline

Install fatal handler with default options.

This is meant to be easy to call from python via ctypes.

Definition at line 160 of file FixLArElecCalib.cxx.

                                                        {
 
   ATH_MSG_INFO ( " in stop " );
 
  const EventContext& ctx = Gaudi::Hive::currentContext();  
 
   SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey, ctx};
   const LArOnOffIdMapping* cabling{*cablingHdl};
   if(!cabling) {
     ATH_MSG_ERROR( "Do not have cabling mapping from key " << m_cablingKey.key() );
     return StatusCode::FAILURE;
   }
   SG::ReadCondHandle<LArCalibLineMapping> clHdl{m_CLKey, ctx};
   const LArCalibLineMapping *clCont {*clHdl};
   if(!clCont) {
     ATH_MSG_ERROR( "Do not have calib line mapping !!!" );
     return StatusCode::FAILURE;
   }
 
 
    if(m_fixFlag==1) 
    return fix1(cabling); 
    if(m_fixFlag==2) 
    return fix2(cabling); 
    if(m_fixFlag==3) 
    return fix3(cabling); 
    if(m_fixFlag==4) 
    return fix4(cabling); 
    if(m_fixFlag==5) 
    return fix5(cabling); 
    if(m_fixFlag==6) 
    return fix6(cabling); 
    if(m_fixFlag==7) 
    return fix7(cabling); 
    if(m_fixFlag==8) 
    return fix8(cabling); 
    if(m_fixFlag==9)
    return fix9(cabling); 
    if(m_fixFlag==10)
    return fix10(cabling); 
    if(m_fixFlag==11)
    return fix11(cabling); 
    if(m_fixFlag==12)
    return fix12(cabling); 
    if(m_fixFlag==13)
    return fix13(cabling); 
    if(m_fixFlag==14)
    return fix14(cabling); 
    if(m_fixFlag==15)
    return fix15(); 
    if(m_fixFlag==16)
    return fix16(); 
    if(m_fixFlag==17)
    return fix17(cabling); 
    if(m_fixFlag==18)
    return fix18(clCont); 
 
    return StatusCode::SUCCESS; 
 
}

ATLAS_NOT_THREAD_SAFE() [2/3]

StatusCode FixLArElecCalib::fix17 ATLAS_NOT_THREAD_SAFE

(

const LArOnOffIdMapping *

cabling

)

Definition at line 443 of file FixLArElecCalib.cxx.

                                                                                        {
 
    ATH_MSG_INFO ( " in fix3() " );
    std::string filename ;
 
    // updated for 13.0.20
    if(m_g4Phys=="QGSP_EMV"){
      filename = "fsampl_1302_qgsp-emv.txt";
    } 
    else
      if( m_g4Phys=="QGSP_BERT"){
    filename = "fsampl_1303_qgsp-bert.txt";
      }
      else {
    ATH_MSG_ERROR ( " wrong option  "<<m_g4Phys );
    return StatusCode::FAILURE; 
      }
 
    // Fix3 is for updating the FCAL conditions data after IdFix7
    // Input should be MC Conditions data with DetDescrVersion=ATLAS-DC3-05
    
    ATH_CHECK( updateEMfSampl(filename, cabling ) );
    ATH_CHECK( updateHADfSampl(cabling) );
    return StatusCode::SUCCESS;
}

ATLAS_NOT_THREAD_SAFE() [3/3]

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

Definition at line 584 of file FixLArElecCalib.cxx.

                                                                                                                             {
 
   ATH_MSG_INFO ( " in updateEMfSampl() " );
 
   // this method updates the EM fSample, assuming symmetry. 
   // input is the text file.
 
   const LArfSamplMC * fsampl_c = nullptr;
   ATH_CHECK( detStore()->retrieve(fsampl_c) );
   LArfSamplMC* fsampl = const_cast<LArfSamplMC*>(fsampl_c);
 
   // read in the file
 
   ATH_MSG_INFO(" opening file "<< filename );
   std::ifstream infile( filename);
   if(! infile.good() ) 
     {
       ATH_MSG_ERROR(" fail to open file ");
       return StatusCode::FAILURE; 
     }
   int   det,samp,reg,eta; 
   float value ; 
   int n = 0; 
   while ( infile>>det>>samp>>reg>>eta>>value )
    {
    ++n; 
    ATH_MSG_INFO(" det,samp,reg,eta,value="
                     <<det<<" " 
                     <<samp<<" " 
                     <<reg<<" " 
                     <<eta<<" " 
                     <<value<<" " );
 
    Identifier id ; 
    if ( det==1 && samp==1 && reg==0 && eta==0 ){
       // eta=0 for strip, not connected, but keep it
           // Removed disconnected channels: 03/2016 RDS 
           // ATH_MSG_INFO(" disconnected strip ");
       // id = m_em_idhelper->disc_channel_id(det,samp,reg,eta,0); 
        }else
     id = m_em_idhelper->channel_id(det,samp,reg,eta,0); 
 
    HWIdentifier hid = cabling->createSignalChannelID(id);
    const LArfSamplComplete::LArCondObj & t = fsampl->get(hid,0); 
    std::string id_str = m_online_idhelper->print_to_string(hid); 
    ATH_MSG_INFO(" online id = "<<id_str);
 
    if( t.isEmpty() ) 
    {
      ATH_MSG_ERROR(" No existing conditions data " );
          ATH_MSG_ERROR( m_em_idhelper->print_to_string(id) );
      // return StatusCode::FAILURE ;
      continue ; 
    }
    ATH_MSG_INFO(" Old fsampl = "<< t.m_fSampl<< " " );
 
    LArfSamplComplete::LArCondObj& t2 = const_cast<LArfSamplComplete::LArCondObj&>(t); 
    t2.m_fSampl = value ; 
 
    const LArfSamplComplete::LArCondObj & t3 = fsampl->get(hid,0); 
    ATH_MSG_INFO(" New fSampl = "<< t3.m_fSampl<< " " );
    }
 
   ATH_MSG_INFO("  Number of entries changes =  " <<n );
   return StatusCode::SUCCESS;
}