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 <fstream>
#include "TTree.h"
#include "TFile.h"
#include "AthenaPoolUtilities/AthenaAttributeList.h"
#include "CoralBase/Blob.h"

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 ( void )

inline

Install fatal handler with default options.

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

Install fatal handler with default options.

getLorentzAngle() Read LorentzAngle from HIST and write out into local DB

getBSErrors() Read BSErrors from Monitoring HIST and write out into local DB

getEfficiency() Read Efficiency from Monitoring HIST and write out into local DB

getRawOccupancy() Read RawOccupancy from Monitoring HIST and write out into local DB

getNoiseOccupancy() Read NoiseOccupancy from HIST and write out into local DB

getNoisyStrip() Find noisy strips from hitmaps and write out into xml/db formats

beginning of the loop of channels

bad bit newly found

known bad bit

for low noisy cells

for high noisy cells

0.01 is used to scale "PER" to the same order of magnitude to "SIG"

smaller deviation: distorted

checking TmaxAmp, Not mixed with MaxAmp and Width

channel information output

Only dead or distorted, or short known BCs are considered below.

index of bc

Definition at line 158 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 441 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 582 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.c_str() );
   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;
}