ALFA_Raw2Digit Class Reference

#include <ALFA_Raw2Digit.h>

Inheritance diagram for ALFA_Raw2Digit:

Collaboration diagram for ALFA_Raw2Digit:

Public Member Functions
	ALFA_Raw2Digit (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
virtual StatusCode	sysStart () override
	Handle START transition. More...
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
MsgStream &	msg (const MSG::Level lvl) const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	recordCollection ()
StatusCode	recordODCollection ()
StatusCode	mapping ()
StatusCode	ChanPlot (uint32_t pot_n, uint32_t iLayer, uint32_t MarChan, uint32_t &chan_i, uint32_t &chan_j)
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
std::string	m_ALFA_RawDataCollectionKey
int	m_nEvents
int	m_mbID
int	m_pmfID
int	m_fibID
int	m_MAPMTChan
int	m_MapChan
int	m_pot
int	m_side
int	m_plate
int	m_fiber
int	m_ODpot
int	m_ODside
int	m_ODplate
int	m_ODfiber
int	m_maroc2fiber [8][20][64]
int	m_maroc2mapmt [8][20][64]
int	m_mapmt2maroc [8][20][64]
int	m_pmf2layer [8][23]
int	m_OD_pmf2layer [8][4]
int	m_OD_pmf_maroc2fiber [8][4][64]
int	m_OD_pmf_maroc2side [8][4][64]
int	m_mb2det [8]
TH2I *	m_All_Trigger
TH2D *	m_All_Trigger_norm
TH1F *	m_act_lay_h [8]
TH1F *	m_hit_lay_h [8]
uint32_t	m_lumi_block
uint32_t	m_lumi_block_old
uint16_t	m_charge_1 [8]
uint16_t	m_charge_2 [8]
std::vector< bool >	m_trigger_pattern
bool	m_pot_b [8]
bool	m_layer_b [8][20]
bool	m_fiber_b [8][20][64]
bool	m_trigger_pattern_b [8][16]
Int_t	m_n_hits_lay [8][20]
Int_t	m_event_no
TTree *	m_tree
TFile *	m_rootOutput
TString	m_rootOutputFileName
Int_t	m_pmf_chan_hit_counter [8][23][64]
Int_t	m_hit_counter [8]
ALFA_DigitCollection *	m_digitCollection
ALFA_Digit *	m_digitObject
ALFA_ODDigitCollection *	m_ODdigitCollection
ALFA_ODDigit *	m_ODdigitObject
std::string	m_key_DigitCollection
std::string	m_key_ODDigitCollection
std::ifstream	m_inChan
std::ifstream	m_inPMF
std::ifstream	m_inOD
std::ifstream	m_inDet
std::string	m_mapname
std::stringstream	m_s
std::string	m_h_name1
std::string	m_h_name2
std::string	m_h_name3
std::string	m_h_loc1
std::string	m_h_loc2
std::string	m_h_loc3
std::string	m_CanName
int	m_WordId_count
uint32_t	m_chan_i
uint32_t	m_chan_j
std::string	m_strMeasuredDataType
std::ofstream	m_toFile
DataObjIDColl	m_extendedExtraObjects
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 51 of file ALFA_Raw2Digit.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

ALFA_Raw2Digit()

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

Definition at line 12 of file ALFA_Raw2Digit.cxx.

: AthAlgorithm(name, pSvcLocator)
{
 
    m_key_DigitCollection  = "ALFA_DigitCollection";
    m_key_ODDigitCollection  = "ALFA_ODDigitCollection";
 
    declareProperty("ALFA_RawDataCollectionKey", m_ALFA_RawDataCollectionKey = "ALFA_RawData");
 
    //allowed values (m_strMeasuredDataType="testbeam", "tunnel"), default value is "tunnel"
    //this property is needed due to different mapping in the testbeam and the tunnel data
    declareProperty("MeasuredDataType", m_strMeasuredDataType = "tunnel");
 
 
    m_All_Trigger_norm = nullptr;
    m_All_Trigger = nullptr;
    m_ODdigitCollection = nullptr;
    m_ODdigitObject = nullptr;
    m_digitCollection = nullptr;
    m_digitObject = nullptr;
//  m_histSvc = NULL;
    m_rootOutput = nullptr;
    m_tree = nullptr;
 
    m_nEvents = 0;
    m_mbID      = 0;
    m_pmfID     = 0;
    m_fibID     = 0;
    m_MAPMTChan = 0;
    m_MapChan   = 0;
    m_pot     = 0;
    m_side    = 0;
    m_plate   = 0;
    m_fiber   = 0;
    m_ODpot   = 0;
    m_ODside  = 0;
    m_ODplate = 0;
    m_ODfiber = 0;
 
    memset(&m_act_lay_h, 0, sizeof(m_act_lay_h));
    memset(&m_hit_lay_h, 0, sizeof(m_hit_lay_h));
    m_lumi_block     = 0;
    m_lumi_block_old = 0;
    m_event_no       = 0;
    m_WordId_count   = 0;
    m_chan_i         = 0;
    m_chan_j         = 0;
}

Member Function Documentation

ChanPlot()

StatusCode ALFA_Raw2Digit::ChanPlot ( uint32_t  pot_n, uint32_t  iLayer, uint32_t  MarChan, uint32_t &  chan_i, uint32_t &  chan_j  )

private

Definition at line 741 of file ALFA_Raw2Digit.cxx.

{
    // returs values in range 1-8; input is 0-63!
//  if (int(fmod(double(m_maroc2mapmt[pot_n][MarChan]+1),double(8)))>0)
    if (int(fmod(double(m_maroc2mapmt[pot_n][iLayer][MarChan]+1),double(8)))>0)
    {
//      chan_i=int((m_maroc2mapmt[pot_n][MarChan]+1)/8)+1;
//      chan_j=9-int(fmod(double(m_maroc2mapmt[pot_n][MarChan]+1),double(8)));
        chan_i=int((m_maroc2mapmt[pot_n][iLayer][MarChan]+1)/8)+1;
        chan_j=9-int(fmod(double(m_maroc2mapmt[pot_n][iLayer][MarChan]+1),double(8)));
    }
    else
    {
//      chan_i=int((m_maroc2mapmt[pot_n][MarChan]+1)/8);
        chan_i=int((m_maroc2mapmt[pot_n][iLayer][MarChan]+1)/8);
        chan_j=1;
    }
 
//  msg(MSG::FATAL)  << "in funct: m_fibID = " << MarChan << ", chan_i " << chan_i << ", chan_j = " << chan_j <<  endmsg;
    return StatusCode::SUCCESS;
}

declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleKeyArrayType &    )

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareGaudiProperty() [2/4]

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleType &    )

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
  }

declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  t, const SG::NotHandleType &    )

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(t);
  }

declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleBase &  hndl, const std::string &  doc, const SG::VarHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

  {
    this->declare(hndl.vhKey());
    hndl.vhKey().setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKey &  hndl, const std::string &  doc, const SG::VarHandleKeyType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

  {
    this->declare(hndl);
    hndl.setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKeyArray &  hndArr, const std::string &  doc, const SG::VarHandleKeyArrayType &    )

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

  {
 
    // std::ostringstream ost;
    // ost << Algorithm::name() << " VHKA declareProp: " << name 
    //     << " size: " << hndArr.keys().size() 
    //     << " mode: " << hndArr.mode() 
    //     << "  vhka size: " << m_vhka.size()
    //     << "\n";
    // debug() << ost.str() << endmsg;
 
    hndArr.setOwner(this);
    m_vhka.push_back(&hndArr);
 
    Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
    if (p != 0) {
      p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
    } else {
      ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                    << name);
    }
 
    return p;
 
  }

declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc, const SG::NotHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(name, property, doc);
  }

declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc = "none"  )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

  {
    typedef typename SG::HandleClassifier<T>::type htype;
    return declareProperty (name, property, doc, htype());
  }

declareProperty() [6/6]

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                 {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

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.

{ return m_detStore; }

evtStore() [1/2]

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.

{ return m_evtStore; }

evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode ALFA_Raw2Digit::execute

(

)

Definition at line 174 of file ALFA_Raw2Digit.cxx.

{
    ATH_MSG_DEBUG ("ALFA_Raw2Digit: BEGINING of execute");
 
    m_nEvents++;
    msg(MSG::DEBUG) << "event # "<< m_nEvents << endmsg;
//  m_toFile << "Event " << m_nEvents << std::endl;
 
    // retrieve the RawData RDO container from the TDS
    const ALFA_RawDataContainer* container = nullptr;
    StatusCode sc_read = evtStore()->retrieve(container, m_ALFA_RawDataCollectionKey);
 
    if ( sc_read != StatusCode::SUCCESS ) {
        msg(MSG::ERROR) << "ALFA_Raw2Digit: Could not find container" << endmsg;
        return StatusCode::SUCCESS;
    }
 
    else if ( sc_read == StatusCode::SUCCESS ) {
        msg(MSG::DEBUG)  << "ALFA_Raw2Digit: Raw data container retrieved" << endmsg;
    }
 
    msg(MSG::DEBUG) << "MD record_collection()" << endmsg;
 
    StatusCode sc = recordCollection();
    if  (sc.isFailure()) {
        msg(MSG::WARNING) << "ALFA_Raw2Digit: recordCollection failed" << endmsg;
        return StatusCode::SUCCESS;
    }
 
    msg(MSG::DEBUG) << "OD record_collection()" << endmsg;
 
    sc = recordODCollection();
    if  (sc.isFailure()) {
        msg(MSG::WARNING) << "ALFA_Raw2Digit, recordODCollection failed" << endmsg;
        return StatusCode::SUCCESS;
    }
 
 
    for (unsigned int i=0;i<8;i++) {
     m_pot_b[i] = false;
     for (unsigned int j=0;j<20;j++) {
         m_layer_b[i][j] = false;
         for (unsigned int k=0;k<64;k++) {
             m_fiber_b[i][j][k] = false;
         }
     }
     for (unsigned int k=0;k<16;k++) {
         m_trigger_pattern_b[i][k] = false;
     }
    }
 
 
    m_event_no = m_nEvents;
 
 
    m_lumi_block = container->GetLumiBlock();
//  std::cout << "m_lumi_block = " << m_lumi_block << std::endl;
 
    // loop over container start
    ALFA_RawDataContainer::const_iterator RawData_Collection_Beg = container->begin();
    ALFA_RawDataContainer::const_iterator RawData_Collection_End = container->end();
 
    //loop over collection (container) with hits (i.e. over 1 event)
    for(;RawData_Collection_Beg!=RawData_Collection_End; ++RawData_Collection_Beg)
    {
        m_mbID = ((*RawData_Collection_Beg)->GetMBId_POT());
        msg(MSG::DEBUG) << "ALFA_Raw2Digit: MBId POT= " << m_mbID <<  endmsg;
 
//      std::cout << "*******    Time stamp = " << (*RawData_Collection_Beg)->GetTimeStamp_POT() << std::endl;
//      time_t rawtime = (*RawData_Collection_Beg)->GetTimeStamp_POT();
//      struct tm * timeinfo;
//      timeinfo = gmtime ( &rawtime );
//      printf ( "The date/time is: %m_s", asctime (timeinfo) );
 
        m_pot_b[m_mbID-1] = true;
  //      m_event_no = ((*RawData_Collection_Beg)->GetEventCount_POT());
 
        m_charge_1[m_mbID-1] = ((*RawData_Collection_Beg)->Get_ADC1_POT());
        m_charge_2[m_mbID-1] = ((*RawData_Collection_Beg)->Get_ADC2_POT());
        for (unsigned int i=0;i<16;i++) {
      m_trigger_pattern_b[m_mbID-1][i] = ((*RawData_Collection_Beg)->Get_pattern_POT())[i];
//    msg(MSG::ERROR) << " pattern: bit " << i << ", value = "<< m_trigger_pattern_b[m_mbID-1][i] << " in MB = " << m_mbID-1 << endmsg;
        }
 
// loop over collection start
        ALFA_RawDataCollection::const_iterator p_RawData_Beg=(*RawData_Collection_Beg)->begin();//pointer to RawData data vector
        ALFA_RawDataCollection::const_iterator p_RawData_End=(*RawData_Collection_Beg)->end();
 
        for(;p_RawData_Beg!=p_RawData_End; ++p_RawData_Beg)
        {
//          msg(MSG::DEBUG) << " all, pmfId = " << (*p_RawData_Beg)->GetPMFId_PMF() <<  ", MBId = " << (*p_RawData_Beg)->GetMBId_PMF() << endmsg;
 
            if ((*p_RawData_Beg)->GetPMFId_PMF()>0 && ((*p_RawData_Beg)->GetWordId_PMF())<24)
            {
                m_pmfID= (*p_RawData_Beg)->GetPMFId_PMF();
 
                if (((*p_RawData_Beg)->GetWordId_PMF())==0) m_WordId_count=1;
                else  m_WordId_count=0;
 
                if (m_WordId_count==1){
                    msg(MSG::DEBUG)  << " pmfId = " << (*p_RawData_Beg)->GetPMFId_PMF() <<  ", MBId = " << (*p_RawData_Beg)->GetMBId_PMF() <<endmsg;
                }
 
//              msg(MSG::DEBUG) << " m_WordId_count " << m_WordId_count << endmsg;
 
                for (unsigned int i=0;i<16;i++)
                {
                    if ((*p_RawData_Beg)->Get_Chan(i)!=100)
                    {
                        msg(MSG::DEBUG)  <<"ALFA_Raw2Digit, maroc chanel # = " << (*p_RawData_Beg)->Get_Chan(i)<< endmsg;
 
                        m_fibID = (*p_RawData_Beg)->Get_Chan(i);
 
                        // in Sara Diglio Byte Stream Converter PMF are counted from 1, then OD detectors have values of PMF in the range 2-4; fibers start at 0
                        if ((m_pmfID == 1) || (m_pmfID > 4))
                        {
                            m_pot = m_mb2det[m_mbID-1];
//                          m_pot = m_mbID-1;
                            m_plate = m_pmf2layer[m_pot][m_pmfID-1];
                            // changed by Petr - 19.12.2012 -----------------------------
//                          m_fiber = 63-m_maroc2fiber[m_pot][m_fibID];
                            m_fiber = m_maroc2fiber[m_pot][m_plate][m_fibID];
 
                            m_layer_b[m_mbID-1][m_plate] = true;
                            m_fiber_b[m_mbID-1][m_plate][m_fiber] = true;
//                          msg(MSG::ERROR) << " bim ho, lumi = " << m_lumi_block << ", pot = "<< m_mbID-1 << ", plate " << m_plate << endmsg;
                            m_n_hits_lay[m_mbID-1][m_plate] = m_n_hits_lay[m_mbID-1][m_plate] + 1;
//                          msg(MSG::ERROR) << " number of hits " << m_n_hits_lay[m_mbID-1][m_plate] << endmsg;
 
//                          m_digitCollection->push_back(new ALFA_Digit(m_pot,m_plate,m_fiber));
                            m_digitCollection->push_back(new ALFA_Digit(m_mbID-1,m_plate,m_fiber));
 
                            // changed by Petr - 19.12.2012 -----------------------------
//                          m_MapChan = m_maroc2mapmt[m_pot][m_fibID];
                            m_MapChan = m_maroc2mapmt[m_pot][m_plate][m_fibID];
                            msg(MSG::DEBUG)  << "pot = " << m_pot << ", pmf = " << m_pmfID-1 << ", maroc = " << m_fibID << ", mapmt = " << m_MapChan <<  endmsg;
 
//                          PMFvsMAPMT[m_pot]->Fill(m_MapChan,m_pmfID-1);
//                          PMFvsMAROC[m_pot]->Fill(m_fibID,m_pmfID-1);
                            // changed by Petr - 19.12.2012 -----------------------------
//                          ChanPlot(m_pot,m_fibID,m_chan_i,m_chan_j);
//                          ChanPlot(m_pot,m_plate,m_fibID,m_chan_i,m_chan_j);
//                          MapmtHit[m_pot][m_pmfID-1]->Fill(m_chan_i,m_chan_j);
                            msg(MSG::DEBUG)  << "fibID = " << m_fibID << ", chan_i " << m_chan_i << ", chan_j = " << m_chan_j << endmsg;
 
                            m_pmf_chan_hit_counter[m_pot][m_pmfID-1][m_fibID]++;
                            m_hit_counter[m_pot]++;
 
                            msg(MSG::DEBUG)  << "ALFA_Raw2Digit: new digit (softsim numbering-starting from 0) in the collection, pot = " << m_pot << ", plate = " << m_plate << ", fiber = " << m_fiber << endmsg;
                        }
                        else
                        {
                            m_ODpot = m_mb2det[m_mbID-1];
//                          m_ODpot = m_mbID-1;
                            m_ODside  = m_OD_pmf_maroc2side[m_ODpot][m_pmfID-1][m_fibID];
                            m_ODplate = m_OD_pmf2layer[m_ODpot][m_pmfID-1];
                            m_ODfiber = m_OD_pmf_maroc2fiber[m_ODpot][m_pmfID-1][m_fibID];
 
                            // changed by Petr - 19.12.2012 ----------------------------- commented, is it needed by ODs? TODO: to correct m_MapChan for OD case
//                          m_MapChan = m_maroc2mapmt[m_ODpot][m_fibID];
//                          msg(MSG::DEBUG)  << " pot = " << m_ODpot << ", pmf = " << m_pmfID-1 << ", maroc = " << m_fibID << ", mapmt = " << m_MapChan << endmsg;
 
                            if (m_ODfiber!=98)
                            {
//                              m_ODfiber = (m_ODside==0)? m_ODfiber : 29-m_ODfiber;
 
//                              m_ODdigitCollection->push_back(new ALFA_ODDigit(m_ODpot,m_ODside,m_ODplate,m_ODfiber));
                                m_ODdigitCollection->push_back(new ALFA_ODDigit(m_mbID-1,m_ODside,m_ODplate,m_ODfiber));
 
//                              msg(MSG::DEBUG) << "pot, side, plate, fiber = " << m_ODpot << ", " << m_ODside << ", " << m_ODplate << ", " << m_ODfiber << endmsg;
 
                            }
 
 
//                          PMFvsMAPMT[m_ODpot]->Fill(m_MapChan,m_pmfID-1);
//                          PMFvsMAROC[m_ODpot]->Fill(m_fibID,m_pmfID-1);
 
                            // changed by Petr - 19.12.2012 ----------------------------- commented, is it needed by ODs? TODO: to correct ChanPlot for OD case
//                          ChanPlot(m_ODpot,m_fibID,m_chan_i,m_chan_j);
//                          MapmtHit[m_ODpot][m_pmfID-1]->Fill(m_chan_i,m_chan_j);
                            msg(MSG::DEBUG)  << "fibID = " << m_fibID << ", chan_i " << m_chan_i << ", chan_j = " << m_chan_j << endmsg;
 
                            m_pmf_chan_hit_counter[m_ODpot][m_pmfID-1][m_fibID]++;
                            m_hit_counter[m_ODpot]++;
 
                            msg(MSG::DEBUG)  << "ALFA_Raw2Digit:  new OD digit (softsim numbering-starting from 0) in the collection,  pot = " << m_ODpot << ", side = " << m_ODside << ", plate = " << m_ODplate << ", fiber = " << m_ODfiber << endmsg;
                        }
                    }
                        // */
                }
            } //if (PFM!=0)
            else
            {
                msg(MSG::DEBUG) << "ALFA_Raw2Digit: More than n Raw Data " <<endmsg;
                break;
            }
        } // loop raw data
    
        msg(MSG::DEBUG)  << "ALFA_Raw2Digit: finish loop over raw data " <<   endmsg;
    } // loop collection
 
    msg(MSG::DEBUG)  << "ALFA_Raw2Digit: finish loop over collection " <<   endmsg;
 
//  m_tree->Fill();
 
    return StatusCode::SUCCESS;
}

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.

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

finalize()

StatusCode ALFA_Raw2Digit::finalize

(

)

Definition at line 384 of file ALFA_Raw2Digit.cxx.

{
    ATH_MSG_DEBUG ("ALFA_Raw2Digit:finalize()");
 
//  m_toFile.close();
//  m_tree->Write();
 
/*      
    TCanvas *canvas = new TCanvas("MAPMT_activity", "MAPMT_activity");
//  TCanvas *tcDistance = new TCanvas("tcDistance", "tcDistance");
 
    gStyle->SetPalette(1);
 
    for (int j=0;j<8;j++){
        for (int k=0;k<23;k++){
            for (int m=1;m<9;m++){
                for (int n=1;n<9;n++){
                    MAPMT_activity_h[j][k]->SetBinContent(m,n,double(MapmtHit[j][k]->GetBinContent(m,n))/double(m_hit_counter[j]));
                }
            }
            for (int i=0;i<64;i++){
                MAPMT_activity[j]->SetPoint(k*64+i,k*64+i,double(m_pmf_chan_hit_counter[j][k][i])/double(m_hit_counter[j]));
            }
        }
 
        m_s.str("");
        m_s << j+1;
 
        m_h_name1 = "MAPMT_activity_p_";
        m_h_name1 += m_s.str();
 
        canvas->Clear();
        canvas->cd();
 
        MAPMT_activity[j]->SetTitle(m_h_name1.c_str());
        MAPMT_activity[j]->SetMarkerStyle(6);
        MAPMT_activity[j]->Draw("AP");
        MAPMT_activity[j]->SetName(m_h_name1.c_str());
 
        canvas->cd()->Update();
 
        m_h_name1 += ".eps";
        canvas->Print(m_h_name1.c_str());
    }
 
    delete canvas;
 
    for (int j=0;j<8;j++){
        TCanvas *POT_View = new TCanvas("POT_View","POT_View");
        POT_View->Divide(5,5);
 
        gStyle->SetPalette(1);
 
        for (unsigned int i=0;i<4;i++)
        {
            POT_View->cd(i+1);
            //  MapmtEff[i]->GetZaxis()->SetRangeUser(0,max_trig_eff_old[0]);
            MAPMT_activity_h[j][i]->Draw("COLZ,A");
            MAPMT_activity_h[j][i]->SetStats(0);
        }
    
        for (unsigned int i=4;i<11;i++)
        {
            POT_View->cd(i+2);
            //  MapmtEff[i]->GetZaxis()->SetRangeUser(0,max_trig_eff_old[0]);
            MAPMT_activity_h[j][i]->Draw("COLZ,A");
            MAPMT_activity_h[j][i]->SetStats(0);
        }
    
        for (unsigned int i=11;i<23;i++)
        {
            POT_View->cd(i+3);
            //  MapmtEff[i]->GetZaxis()->SetRangeUser(0,max_trig_eff_old[0]);
            MAPMT_activity_h[j][i]->Draw("COLZ,A");
            MAPMT_activity_h[j][i]->SetStats(0);
        }
      
        POT_View->cd();
 
        m_s.str("");
        m_s << j+1;
        m_CanName = "h_MAPMT_activity_p_";
        m_CanName += m_s.str();
        m_CanName+=".eps";
 
        POT_View->Print(m_CanName.c_str());
        POT_View->Close();
    }
  
*/
 
//  m_rootOutput->Write();
//  m_rootOutput->Close();
 
    return StatusCode::SUCCESS;
}

initialize()

StatusCode ALFA_Raw2Digit::initialize

(

)

Definition at line 61 of file ALFA_Raw2Digit.cxx.

{
    ATH_MSG_DEBUG ("ALFA_Raw2Digit:initialize()");
 
    if (m_strMeasuredDataType=="tunnel")
    {
        m_mb2det[0] = 1; m_mb2det[1] = 0; m_mb2det[2] = 7; m_mb2det[3] = 2;
        m_mb2det[4] = 4; m_mb2det[5] = 5; m_mb2det[6] = 6; m_mb2det[7] = 3;
    }
    else if (m_strMeasuredDataType=="testbeam")
    {
        m_mb2det[0] = 0; m_mb2det[1] = 1; m_mb2det[2] = 2; m_mb2det[3] = 3;
        m_mb2det[4] = 4; m_mb2det[5] = 5; m_mb2det[6] = 6; m_mb2det[7] = 7;
    }
 
    m_nEvents = 0;
    m_MAPMTChan = -1;
 
    for (int m=0;m<8;m++)
    {
        m_hit_counter[m] = 0;
        for (int k=0;k<23;k++)
        {
            for (int i=0;i<64;i++)
            {
                m_pmf_chan_hit_counter[m][k][i] = 0;
            }
        }
    }
 
    StatusCode sc = mapping();
 
    if (sc.isFailure()) { msg (MSG::FATAL) << " Could not open mapping " << endmsg; return sc; }
    else                  msg (MSG::DEBUG) << " mapping is fine " << endmsg;
 
//  sc = service("THistSvc", m_histSvc);
 
    // retrieve ThistSvc ...
//  if(sc.isFailure())
//  {
//      msg(MSG::FATAL) << "Unable to retrieve pointer to THistSvc" << endmsg;
//      return sc;
//  }
 
/*
    m_rootOutputFileName = "ALFA_toy_tree_tom.root";
    m_rootOutput = new TFile(m_rootOutputFileName,"RECREATE");
    m_tree = new TTree("MyTree", "ALFA_toy_tree");
 
    TBranch *branch;
    branch = m_tree->Branch("m_event_no", &m_event_no,"m_event_no/I"); branch->SetAutoDelete(kFALSE);
    branch = m_tree->Branch("m_lumi_block", &m_lumi_block, "m_lumi_block/I"); branch->SetAutoDelete(kFALSE);
    branch = m_tree->Branch("pot_b", m_pot_b, "pot_b[8]/O"); branch->SetAutoDelete(kFALSE);
    branch = m_tree->Branch("layer_b", m_layer_b, "layer_b[8][20]/O"); branch->SetAutoDelete(kFALSE);
    branch = m_tree->Branch("fiber_b", m_fiber_b, "fiber_b[8][20][64]/O"); branch->SetAutoDelete(kFALSE);
    branch = m_tree->Branch("trigger_pattern_b", m_trigger_pattern_b, "trigger_pattern[8][16]/O"); branch->SetAutoDelete(kFALSE);
    branch = m_tree->Branch("charge_1", &m_charge_1, "charge_1[8]/s"); branch->SetAutoDelete(kFALSE);
    branch = m_tree->Branch("charge_2", &m_charge_2, "charge_2[8]/s"); branch->SetAutoDelete(kFALSE);
*/
 
//  for (unsigned int j=0; j<8; j++)
//  {
//      m_h_name1 = "MAPMT_activity_p_";
//      m_h_name2 = "PMFvsMAPMT_p_";
//      m_h_name3 = "PMFvsMAROC_p_";
 
//      m_h_loc1 = "/AANT/MAPMT_activity_p_";
//      m_h_loc2 = "/AANT/PMFvsMAPMT_p_";
//      m_h_loc3 = "/AANT/PMFvsMAROC_p_";
 
//      m_s.str("");
//      m_s << j+1;
 
//      m_h_name1 += m_s.str();
//      m_h_name2 += m_s.str();
//      m_h_name3 += m_s.str();
 
//      m_h_loc1 += m_s.str();
//      m_h_loc2 += m_s.str();
//      m_h_loc3 += m_s.str();
     
/*  
        m_rootOutputFileName = "ALFA_toy_tree_tom.root";
        m_rootOutput = new TFile(m_rootOutputFileName,"RECREATE");
        m_tree = new TTree("MyTree", "ALFA_toy_tree");
 
        MAPMT_activity[j] = new TGraph(1472);
//      MAPMT_activity[j]->SetTitle("Trigger efficiency for Upper pot;MAROC Channel; Trigger efficiency");
        MAPMT_activity[j]->SetTitle(m_h_name1.c_str());
        MAPMT_activity[j]->SetMarkerStyle(6);
        MAPMT_activity[j]->Draw("AP");
        MAPMT_activity[j]->SetName(m_h_name1.c_str());
 
//      sc = m_histSvc->regHist(m_h_loc1.c_str(),MAPMT_activity[j]);
 
        PMFvsMAPMT[j]     = new TH2I(m_h_name2.c_str(),m_h_name2.c_str(),64,-0.5,63.5,23,-0.5,22.5);
        sc = m_histSvc->regHist(m_h_loc2.c_str(),PMFvsMAPMT[j]);
 
        PMFvsMAROC[j]     = new TH2I(m_h_name3.c_str(),m_h_name3.c_str(),64,-0.5,63.5,23,-0.5,22.5);
        sc = m_histSvc->regHist(m_h_loc3.c_str(),PMFvsMAROC[j]);
 
        for (unsigned int i=0; i<23; i++)
        {
            MapmtHit[j][i] = new TH2I(m_h_name1.c_str(),m_h_name1.c_str(),8,0.5,8.5,8,0.5,8.5);
            MAPMT_activity_h[j][i] = new TH2D("","",8,0.5,8.5,8,0.5,8.5);
            sc = m_histSvc->regHist(m_h_loc1.c_str(),MAPMT_activity_h[j][i]);
        }
*/
//  }
 
    return StatusCode::SUCCESS;
}

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.

mapping()

StatusCode ALFA_Raw2Digit::mapping ( )

private

Definition at line 508 of file ALFA_Raw2Digit.cxx.

{
    int MarChan      = 0;
    int MapChan      = 0;
    int FibChan      = 0;
    int iLayer       = 0;
    int PMFNum       = 0;
    int LAYNum       = 0;
    int OD_PMFNum    = 0;
    int OD_Dieter    = 0;
    int OD_LAYNum    = 0;
    int OD_Side      = 0;
    int OD_MarChan   = 0;
    int OD_FibChan   = 0;
    int OD_MaPmtChan = 0;
   
//  int MBnum, DETnum;
 
 
//  m_mapname = "./mapping/motheboard2detector.dat";
//  m_inDet.open(m_mapname.c_str());
//  msg(MSG::DEBUG) << "file name " << m_mapname.c_str() << endmsg;
//  for (unsigned int i=0;i<64;i++) {
//      m_inDet >> MBnum >> DETnum;
//      m_mb2det[MBnum] = DETnum;
//  }
 
    std::fill_n(&m_maroc2fiber[0][0][0], sizeof(m_maroc2fiber)/sizeof(int), -9999.0);
    std::fill_n(&m_maroc2mapmt[0][0][0], sizeof(m_maroc2mapmt)/sizeof(int), -9999.0);
    std::fill_n(&m_mapmt2maroc[0][0][0], sizeof(m_mapmt2maroc)/sizeof(int), -9999.0);
 
 
    for (unsigned int j=0; j<8; j++)
    {
        m_s.str("");
 
//      m_mapname = "mapping/MAROC_MAPMT_FIBER_MD";
 
        // changed by Petr - 19.12.2012 -------------------------------------------------
//      m_mapname = "MAROC_MAPMT_FIBER_MD";
        m_mapname = "MAROC_MAPMT_FIBER_NEW_MD_";
        m_s << j+1;
        m_mapname += m_s.str();
        m_mapname += ".dat";
        
        msg(MSG::DEBUG) << "file name " << m_mapname.c_str() << endmsg;
        
        // ************
        std::string filePath = PathResolver::find_file(m_mapname.c_str(),"DATAPATH", PathResolver::RecursiveSearch);
        if(filePath.length() == 0)
        {
            msg(MSG::FATAL) << " the mapping file MD maroc-mapmt \"" <<  m_mapname.c_str() << "\" not found in Datapath." << endmsg;
            throw std::runtime_error("FATAL: mapping MD maroc-mapmt not found in Datapath.");
        }
        else
        {
            msg(MSG::DEBUG) << "DEBUG: the mapping file MD maroc-mapmt \"" <<  m_mapname.c_str() << "\" found in Datapath." << endmsg;  
            msg(MSG::DEBUG) << "DEBUG: filePath =  " << filePath.c_str() << endmsg;
        }
                // **************   
        
//      m_inChan.open(m_mapname.c_str());
        m_inChan.open(filePath.c_str());        
        if (m_inChan.is_open())
        {
            // changed by Petr - 19.12.2012 ---------------------------------------------
//          for (unsigned int i=0;i<64;i++){
            for (unsigned int i=0;i<1280;i++){
//              m_inChan >> MarChan >> MapChan >> FibChan;
                m_inChan >> iLayer >> MarChan >> MapChan >> FibChan;
//              msg(MSG::DEBUG) << i << " maroc_mapmt_fiber, maroc = " << MarChan << ", map = "<< MapChan << ", fib = " << FibChan << endmsg;
                msg(MSG::DEBUG) << i << " maroc_mapmt_fiber, layer = " << iLayer << ", maroc = " << MarChan << ", map = "<< MapChan << ", fib = " << FibChan << endmsg;
    
                if ((MarChan < 0) || (MarChan > 63))
                {
                    std::string s1("MarChan out of range; file corrupted");
                    throw std::runtime_error(s1);   
                }
                
                if ((MapChan < 0) || (MapChan > 63))
                {
                    std::string s2("MapChan out of range; file corrupted");
                    throw std::runtime_error(s2);   
                }           
                
                if ((iLayer < 0) || (iLayer > 19))
                {
                    std::string s2("Layer out of range; file corrupted");
                    throw std::runtime_error(s2);
                }
 
//              m_maroc2fiber[j][MarChan] = FibChan;
//              m_maroc2mapmt[j][MarChan] = MapChan;
//              m_mapmt2maroc[j][MapChan] = MarChan;
                m_maroc2fiber[j][iLayer][MarChan] = FibChan;
                m_maroc2mapmt[j][iLayer][MarChan] = MapChan;
                m_mapmt2maroc[j][iLayer][MapChan] = MarChan;
            }
            m_inChan.close();
        }
        else
        {
            msg(MSG::ERROR) << "the file " << m_mapname.c_str() << " was not open" << endmsg;
            return StatusCode::FAILURE;
        }
        
//      m_mapname = "mapping/PMF_LAYER_MD";
        m_mapname = "PMF_LAYER_MD";     
        m_mapname += m_s.str();
        m_mapname += ".dat";
        msg(MSG::DEBUG) << "file name " << m_mapname.c_str() << endmsg;
        
        // ************
        filePath = PathResolver::find_file(m_mapname.c_str(),"DATAPATH", PathResolver::RecursiveSearch);
        if(filePath.length() == 0)
        {
            msg(MSG::FATAL) << " the mapping file PMF_LAYER_MD \"" <<  m_mapname.c_str() << "\" not found in Datapath." << endmsg;
            throw std::runtime_error("FATAL: mapping PMF_LAYER_MD not found in Datapath.");
        }
        else
        {
            msg(MSG::DEBUG) << "DEBUG: the mapping file PMF_LAYER_MD \"" <<  m_mapname.c_str() << "\" found in Datapath." << endmsg;    
            msg(MSG::DEBUG) << "DEBUG: filePath =  " << filePath.c_str() << endmsg;
        }
                // **************
 
        m_inPMF.open(filePath.c_str());
        if (m_inPMF.is_open())
        {
            for (unsigned int i=0;i<20;i++)
            {
                m_inPMF >> PMFNum >> LAYNum;
                msg(MSG::DEBUG) << i << " pmf2layer, pmf = " << PMFNum  << ", layer = "<< LAYNum << endmsg;
                
                if ((PMFNum < 0) || ( PMFNum > 22))
                {
                    std::string s3("PMFNum out of range; file corrupted");
                    throw std::runtime_error(s3);   
                }
                
                
                m_pmf2layer[j][PMFNum] = LAYNum;
    //          //   layer2pmf[j][LAYNum] = PMFNum;
            }
            m_inPMF.close();
        }
        else
        {
            msg(MSG::WARNING) << "the file " << m_mapname.c_str() << " was not open" << endmsg;
            return StatusCode::FAILURE; 
        }
 
//      m_mapname="mapping/OD_MAP";
        m_mapname="OD_MAP";     
        m_mapname += m_s.str();
        m_mapname+=".dat";
        msg(MSG::DEBUG) << "file name " << m_mapname.c_str() << endmsg;
 
        // ************
        filePath = PathResolver::find_file(m_mapname.c_str(),"DATAPATH", PathResolver::RecursiveSearch);
        if(filePath.length() == 0)
        {
            msg(MSG::FATAL) << " the mapping file OD_MAP \"" <<  m_mapname.c_str() << "\" not found in Datapath." << endmsg;
            throw std::runtime_error("FATAL: mapping OD_MAP not found in Datapath.");
        }
        else
        {
            msg(MSG::DEBUG) << "DEBUG: the mapping file OD_MAP \"" <<  m_mapname.c_str() << "\" found in Datapath." << endmsg;  
            msg(MSG::DEBUG) << "DEBUG: filePath =  " << filePath.c_str() << endmsg;
        }
                // **************
        
        m_inOD.open(filePath.c_str());
        if (m_inOD.is_open())
        {
            for(unsigned int jlay=0; jlay<3; jlay++)
            {
                for(unsigned int jside=0; jside<2; jside++)
                {
                    for(unsigned int jch=0; jch<32; jch++)
                    {
                        m_inOD >> OD_PMFNum >> OD_Dieter >> OD_LAYNum >> OD_Side >> OD_MaPmtChan >> OD_FibChan; // warning!!! MaPmt, not Maroc!
                        msg(MSG::DEBUG) << " pmf2layer, OD_PMFNum = " << OD_PMFNum  << ", OD_Dieter = "<< OD_Dieter << ", OD_LAYNum = " << OD_LAYNum << ", OD_Side = " << OD_Side << ", OD_MaPmtChan = " << OD_MaPmtChan << ", OD_FibChan = " << OD_FibChan << endmsg;
    
                        if ((OD_PMFNum < 2) || ( OD_PMFNum > 4))
                        {
                            std::string s4("OD_PMFNum out of range; file corrupted");
                            throw std::runtime_error(s4);   
                        }
                        
                        if ((OD_MaPmtChan < 1) || ( OD_MaPmtChan > 64))
                        {
                            std::string s5("OD_MaPmtChan out of range; file corrupted");
                            throw std::runtime_error(s5);   
                        }
                        
                        
//                      m_OD_pmf2layer[j][OD_PMFNum-1] = OD_LAYNum-1;
    //                  OD_layer2pmf[j][OD_LAYNum-1] = OD_PMFNum-1;
                        // changed by Petr - 19.12.2012 ---------------------------------
//                      OD_MarChan = m_mapmt2maroc[j][OD_MaPmtChan-1];
//                      OD_MarChan = m_mapmt2maroc[j][OD_LAYNum-1][OD_MaPmtChan-1];
                        //added to solve ATLAS coverity 13339 (only three OD layers are allowed)
                        if ((OD_LAYNum > 0) && ( OD_LAYNum < 4))
                        {
                            m_OD_pmf2layer[j][OD_PMFNum-1] = OD_LAYNum-1;
                            OD_MarChan = m_mapmt2maroc[j][OD_LAYNum-1][OD_MaPmtChan-1];
                        }
                        else
                        {
                            msg(MSG::DEBUG) << "OD_LAYNum out of bounds" << endmsg;
                        }
 
                        m_OD_pmf_maroc2fiber[j][OD_PMFNum-1][OD_MarChan] = OD_FibChan-1;
                        m_OD_pmf_maroc2side[j][OD_PMFNum-1][OD_MarChan] = OD_Side;
                    }
                }
            }
    
            m_inOD.close();
        }
        else
        {
            msg(MSG::WARNING) << "the file " << m_mapname.c_str() << " was not open" << endmsg;
            return StatusCode::FAILURE;
        }
    }
 
    return StatusCode::SUCCESS;
}

msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

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.

recordCollection()

StatusCode ALFA_Raw2Digit::recordCollection ( )

private

Definition at line 482 of file ALFA_Raw2Digit.cxx.

{
    msg(MSG::DEBUG) << " ALFA_Digitization::recordCollection " << endmsg;
 
    m_digitCollection = new ALFA_DigitCollection();
 
    StatusCode sc = evtStore()->record(m_digitCollection, m_key_DigitCollection);
    if (sc.isFailure()) { msg (MSG::FATAL) << " MD - Could not record the empty digit collection in StoreGate " << endmsg; return sc; }
    else                  msg (MSG::DEBUG) << " MD - Digit collection is recorded in StoreGate " << endmsg;
 
    return sc;
}

recordODCollection()

StatusCode ALFA_Raw2Digit::recordODCollection ( )

private

Definition at line 495 of file ALFA_Raw2Digit.cxx.

{
    msg(MSG::DEBUG) << " ALFA_Digitization::recordODCollection " << endmsg;
 
    m_ODdigitCollection = new ALFA_ODDigitCollection();
 
    StatusCode sc = evtStore()->record(m_ODdigitCollection, m_key_ODDigitCollection);
    if (sc.isFailure()) { msg (MSG::FATAL) << " OD - Could not record the empty digit collection in StoreGate " << endmsg; return sc; }
    else                  msg (MSG::DEBUG) << " OD - Digit collection is recorded in StoreGate " << endmsg;
 
    return sc;
}

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.

  {
    h.renounce();
    PBASE::renounce (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.

                                                          {
    handlesArray.renounce();
  }

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, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

                                       {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;
}

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.

updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_act_lay_h

TH1F* ALFA_Raw2Digit::m_act_lay_h[8]

private

Definition at line 131 of file ALFA_Raw2Digit.h.

m_ALFA_RawDataCollectionKey

std::string ALFA_Raw2Digit::m_ALFA_RawDataCollectionKey

private

Definition at line 76 of file ALFA_Raw2Digit.h.

m_All_Trigger

TH2I* ALFA_Raw2Digit::m_All_Trigger

private

Definition at line 128 of file ALFA_Raw2Digit.h.

m_All_Trigger_norm

TH2D* ALFA_Raw2Digit::m_All_Trigger_norm

private

Definition at line 129 of file ALFA_Raw2Digit.h.

m_CanName

std::string ALFA_Raw2Digit::m_CanName

private

Definition at line 167 of file ALFA_Raw2Digit.h.

m_chan_i

uint32_t ALFA_Raw2Digit::m_chan_i

private

Definition at line 170 of file ALFA_Raw2Digit.h.

m_chan_j

uint32_t ALFA_Raw2Digit::m_chan_j

private

Definition at line 170 of file ALFA_Raw2Digit.h.

m_charge_1

uint16_t ALFA_Raw2Digit::m_charge_1[8]

private

Definition at line 134 of file ALFA_Raw2Digit.h.

m_charge_2

uint16_t ALFA_Raw2Digit::m_charge_2[8]

private

Definition at line 134 of file ALFA_Raw2Digit.h.

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_digitCollection

ALFA_DigitCollection* ALFA_Raw2Digit::m_digitCollection

private

Definition at line 149 of file ALFA_Raw2Digit.h.

m_digitObject

ALFA_Digit* ALFA_Raw2Digit::m_digitObject

private

Definition at line 150 of file ALFA_Raw2Digit.h.

m_event_no

Int_t ALFA_Raw2Digit::m_event_no

private

Definition at line 138 of file ALFA_Raw2Digit.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_fiber

int ALFA_Raw2Digit::m_fiber

private

Definition at line 95 of file ALFA_Raw2Digit.h.

m_fiber_b

bool ALFA_Raw2Digit::m_fiber_b[8][20][64]

private

Definition at line 136 of file ALFA_Raw2Digit.h.

m_fibID

int ALFA_Raw2Digit::m_fibID

private

Definition at line 87 of file ALFA_Raw2Digit.h.

m_h_loc1

std::string ALFA_Raw2Digit::m_h_loc1

private

Definition at line 166 of file ALFA_Raw2Digit.h.

m_h_loc2

std::string ALFA_Raw2Digit::m_h_loc2

private

Definition at line 166 of file ALFA_Raw2Digit.h.

m_h_loc3

std::string ALFA_Raw2Digit::m_h_loc3

private

Definition at line 166 of file ALFA_Raw2Digit.h.

m_h_name1

std::string ALFA_Raw2Digit::m_h_name1

private

Definition at line 165 of file ALFA_Raw2Digit.h.

m_h_name2

std::string ALFA_Raw2Digit::m_h_name2

private

Definition at line 165 of file ALFA_Raw2Digit.h.

m_h_name3

std::string ALFA_Raw2Digit::m_h_name3

private

Definition at line 165 of file ALFA_Raw2Digit.h.

m_hit_counter

Int_t ALFA_Raw2Digit::m_hit_counter[8]

private

Definition at line 145 of file ALFA_Raw2Digit.h.

m_hit_lay_h

TH1F* ALFA_Raw2Digit::m_hit_lay_h[8]

private

Definition at line 132 of file ALFA_Raw2Digit.h.

m_inChan

std::ifstream ALFA_Raw2Digit::m_inChan

private

Definition at line 158 of file ALFA_Raw2Digit.h.

m_inDet

std::ifstream ALFA_Raw2Digit::m_inDet

private

Definition at line 161 of file ALFA_Raw2Digit.h.

m_inOD

std::ifstream ALFA_Raw2Digit::m_inOD

private

Definition at line 160 of file ALFA_Raw2Digit.h.

m_inPMF

std::ifstream ALFA_Raw2Digit::m_inPMF

private

Definition at line 159 of file ALFA_Raw2Digit.h.

m_key_DigitCollection

std::string ALFA_Raw2Digit::m_key_DigitCollection

private

Definition at line 155 of file ALFA_Raw2Digit.h.

m_key_ODDigitCollection

std::string ALFA_Raw2Digit::m_key_ODDigitCollection

private

Definition at line 156 of file ALFA_Raw2Digit.h.

m_layer_b

bool ALFA_Raw2Digit::m_layer_b[8][20]

private

Definition at line 136 of file ALFA_Raw2Digit.h.

m_lumi_block

uint32_t ALFA_Raw2Digit::m_lumi_block

private

Definition at line 133 of file ALFA_Raw2Digit.h.

m_lumi_block_old

uint32_t ALFA_Raw2Digit::m_lumi_block_old

private

Definition at line 133 of file ALFA_Raw2Digit.h.

m_MapChan

int ALFA_Raw2Digit::m_MapChan

private

Definition at line 90 of file ALFA_Raw2Digit.h.

m_mapmt2maroc

int ALFA_Raw2Digit::m_mapmt2maroc[8][20][64]

private

Definition at line 108 of file ALFA_Raw2Digit.h.

m_MAPMTChan

int ALFA_Raw2Digit::m_MAPMTChan

private

Definition at line 89 of file ALFA_Raw2Digit.h.

m_mapname

std::string ALFA_Raw2Digit::m_mapname

private

Definition at line 163 of file ALFA_Raw2Digit.h.

m_maroc2fiber

int ALFA_Raw2Digit::m_maroc2fiber[8][20][64]

private

Definition at line 106 of file ALFA_Raw2Digit.h.

m_maroc2mapmt

int ALFA_Raw2Digit::m_maroc2mapmt[8][20][64]

private

Definition at line 107 of file ALFA_Raw2Digit.h.

m_mb2det

int ALFA_Raw2Digit::m_mb2det[8]

private

Definition at line 118 of file ALFA_Raw2Digit.h.

m_mbID

int ALFA_Raw2Digit::m_mbID

private

Definition at line 85 of file ALFA_Raw2Digit.h.

m_n_hits_lay

Int_t ALFA_Raw2Digit::m_n_hits_lay[8][20]

private

Definition at line 137 of file ALFA_Raw2Digit.h.

m_nEvents

int ALFA_Raw2Digit::m_nEvents

private

Definition at line 80 of file ALFA_Raw2Digit.h.

m_OD_pmf2layer

int ALFA_Raw2Digit::m_OD_pmf2layer[8][4]

private

Definition at line 113 of file ALFA_Raw2Digit.h.

m_OD_pmf_maroc2fiber

int ALFA_Raw2Digit::m_OD_pmf_maroc2fiber[8][4][64]

private

Definition at line 115 of file ALFA_Raw2Digit.h.

m_OD_pmf_maroc2side

int ALFA_Raw2Digit::m_OD_pmf_maroc2side[8][4][64]

private

Definition at line 116 of file ALFA_Raw2Digit.h.

m_ODdigitCollection

ALFA_ODDigitCollection* ALFA_Raw2Digit::m_ODdigitCollection

private

Definition at line 152 of file ALFA_Raw2Digit.h.

m_ODdigitObject

ALFA_ODDigit* ALFA_Raw2Digit::m_ODdigitObject

private

Definition at line 153 of file ALFA_Raw2Digit.h.

m_ODfiber

int ALFA_Raw2Digit::m_ODfiber

private

Definition at line 100 of file ALFA_Raw2Digit.h.

m_ODplate

int ALFA_Raw2Digit::m_ODplate

private

Definition at line 99 of file ALFA_Raw2Digit.h.

m_ODpot

int ALFA_Raw2Digit::m_ODpot

private

Definition at line 97 of file ALFA_Raw2Digit.h.

m_ODside

int ALFA_Raw2Digit::m_ODside

private

Definition at line 98 of file ALFA_Raw2Digit.h.

m_plate

int ALFA_Raw2Digit::m_plate

private

Definition at line 94 of file ALFA_Raw2Digit.h.

m_pmf2layer

int ALFA_Raw2Digit::m_pmf2layer[8][23]

private

Definition at line 110 of file ALFA_Raw2Digit.h.

m_pmf_chan_hit_counter

Int_t ALFA_Raw2Digit::m_pmf_chan_hit_counter[8][23][64]

private

Definition at line 144 of file ALFA_Raw2Digit.h.

m_pmfID

int ALFA_Raw2Digit::m_pmfID

private

Definition at line 86 of file ALFA_Raw2Digit.h.

m_pot

int ALFA_Raw2Digit::m_pot

private

Definition at line 92 of file ALFA_Raw2Digit.h.

m_pot_b

bool ALFA_Raw2Digit::m_pot_b[8]

private

Definition at line 136 of file ALFA_Raw2Digit.h.

m_rootOutput

TFile* ALFA_Raw2Digit::m_rootOutput

private

Definition at line 140 of file ALFA_Raw2Digit.h.

m_rootOutputFileName

TString ALFA_Raw2Digit::m_rootOutputFileName

private

Definition at line 141 of file ALFA_Raw2Digit.h.

m_s

std::stringstream ALFA_Raw2Digit::m_s

private

Definition at line 164 of file ALFA_Raw2Digit.h.

m_side

int ALFA_Raw2Digit::m_side

private

Definition at line 93 of file ALFA_Raw2Digit.h.

m_strMeasuredDataType

std::string ALFA_Raw2Digit::m_strMeasuredDataType

private

Definition at line 179 of file ALFA_Raw2Digit.h.

m_toFile

std::ofstream ALFA_Raw2Digit::m_toFile

private

Definition at line 182 of file ALFA_Raw2Digit.h.

m_tree

TTree* ALFA_Raw2Digit::m_tree

private

Definition at line 139 of file ALFA_Raw2Digit.h.

m_trigger_pattern

std::vector<bool> ALFA_Raw2Digit::m_trigger_pattern

private

Definition at line 135 of file ALFA_Raw2Digit.h.

m_trigger_pattern_b

bool ALFA_Raw2Digit::m_trigger_pattern_b[8][16]

private

Definition at line 136 of file ALFA_Raw2Digit.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.

m_WordId_count

int ALFA_Raw2Digit::m_WordId_count

private

Definition at line 169 of file ALFA_Raw2Digit.h.

---

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

• ALFA_Raw2Digit.h
• ALFA_Raw2Digit.cxx