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 39 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 739 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 172 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 382 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 59 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 506 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,"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,"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,"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 480 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 493 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 94 of file ALFA_Raw2Digit.h.

m_ALFA_RawDataCollectionKey

std::string ALFA_Raw2Digit::m_ALFA_RawDataCollectionKey

private

Definition at line 55 of file ALFA_Raw2Digit.h.

m_All_Trigger

TH2I* ALFA_Raw2Digit::m_All_Trigger {}

private

Definition at line 91 of file ALFA_Raw2Digit.h.

m_All_Trigger_norm

TH2D* ALFA_Raw2Digit::m_All_Trigger_norm {}

private

Definition at line 92 of file ALFA_Raw2Digit.h.

m_CanName

std::string ALFA_Raw2Digit::m_CanName

private

Definition at line 129 of file ALFA_Raw2Digit.h.

m_chan_i

uint32_t ALFA_Raw2Digit::m_chan_i {}

private

Definition at line 132 of file ALFA_Raw2Digit.h.

m_chan_j

uint32_t ALFA_Raw2Digit::m_chan_j {}

private

Definition at line 132 of file ALFA_Raw2Digit.h.

m_charge_1

uint16_t ALFA_Raw2Digit::m_charge_1[8] {}

private

Definition at line 97 of file ALFA_Raw2Digit.h.

m_charge_2

uint16_t ALFA_Raw2Digit::m_charge_2[8] {}

private

Definition at line 97 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 111 of file ALFA_Raw2Digit.h.

m_digitObject

ALFA_Digit* ALFA_Raw2Digit::m_digitObject {}

private

Definition at line 112 of file ALFA_Raw2Digit.h.

m_event_no

Int_t ALFA_Raw2Digit::m_event_no {}

private

Definition at line 101 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 69 of file ALFA_Raw2Digit.h.

m_fiber_b

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

private

Definition at line 99 of file ALFA_Raw2Digit.h.

m_fibID

int ALFA_Raw2Digit::m_fibID {}

private

Definition at line 61 of file ALFA_Raw2Digit.h.

m_h_loc1

std::string ALFA_Raw2Digit::m_h_loc1

private

Definition at line 128 of file ALFA_Raw2Digit.h.

m_h_loc2

std::string ALFA_Raw2Digit::m_h_loc2

private

Definition at line 128 of file ALFA_Raw2Digit.h.

m_h_loc3

std::string ALFA_Raw2Digit::m_h_loc3

private

Definition at line 128 of file ALFA_Raw2Digit.h.

m_h_name1

std::string ALFA_Raw2Digit::m_h_name1

private

Definition at line 127 of file ALFA_Raw2Digit.h.

m_h_name2

std::string ALFA_Raw2Digit::m_h_name2

private

Definition at line 127 of file ALFA_Raw2Digit.h.

m_h_name3

std::string ALFA_Raw2Digit::m_h_name3

private

Definition at line 127 of file ALFA_Raw2Digit.h.

m_hit_counter

Int_t ALFA_Raw2Digit::m_hit_counter[8] {}

private

Definition at line 108 of file ALFA_Raw2Digit.h.

m_hit_lay_h

TH1F* ALFA_Raw2Digit::m_hit_lay_h[8] {}

private

Definition at line 95 of file ALFA_Raw2Digit.h.

m_inChan

std::ifstream ALFA_Raw2Digit::m_inChan

private

Definition at line 120 of file ALFA_Raw2Digit.h.

m_inDet

std::ifstream ALFA_Raw2Digit::m_inDet

private

Definition at line 123 of file ALFA_Raw2Digit.h.

m_inOD

std::ifstream ALFA_Raw2Digit::m_inOD

private

Definition at line 122 of file ALFA_Raw2Digit.h.

m_inPMF

std::ifstream ALFA_Raw2Digit::m_inPMF

private

Definition at line 121 of file ALFA_Raw2Digit.h.

m_key_DigitCollection

std::string ALFA_Raw2Digit::m_key_DigitCollection

private

Definition at line 117 of file ALFA_Raw2Digit.h.

m_key_ODDigitCollection

std::string ALFA_Raw2Digit::m_key_ODDigitCollection

private

Definition at line 118 of file ALFA_Raw2Digit.h.

m_layer_b

bool ALFA_Raw2Digit::m_layer_b[8][20] {}

private

Definition at line 99 of file ALFA_Raw2Digit.h.

m_lumi_block

uint32_t ALFA_Raw2Digit::m_lumi_block {}

private

Definition at line 96 of file ALFA_Raw2Digit.h.

m_lumi_block_old

uint32_t ALFA_Raw2Digit::m_lumi_block_old {}

private

Definition at line 96 of file ALFA_Raw2Digit.h.

m_MapChan

int ALFA_Raw2Digit::m_MapChan {}

private

Definition at line 64 of file ALFA_Raw2Digit.h.

m_mapmt2maroc

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

private

Definition at line 79 of file ALFA_Raw2Digit.h.

m_MAPMTChan

int ALFA_Raw2Digit::m_MAPMTChan {}

private

Definition at line 63 of file ALFA_Raw2Digit.h.

m_mapname

std::string ALFA_Raw2Digit::m_mapname

private

Definition at line 125 of file ALFA_Raw2Digit.h.

m_maroc2fiber

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

private

Definition at line 77 of file ALFA_Raw2Digit.h.

m_maroc2mapmt

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

private

Definition at line 78 of file ALFA_Raw2Digit.h.

m_mb2det

int ALFA_Raw2Digit::m_mb2det[8] {}

private

Definition at line 87 of file ALFA_Raw2Digit.h.

m_mbID

int ALFA_Raw2Digit::m_mbID {}

private

Definition at line 59 of file ALFA_Raw2Digit.h.

m_n_hits_lay

Int_t ALFA_Raw2Digit::m_n_hits_lay[8][20] {}

private

Definition at line 100 of file ALFA_Raw2Digit.h.

m_nEvents

int ALFA_Raw2Digit::m_nEvents {}

private

Definition at line 57 of file ALFA_Raw2Digit.h.

m_OD_pmf2layer

int ALFA_Raw2Digit::m_OD_pmf2layer[8][4] {}

private

Definition at line 83 of file ALFA_Raw2Digit.h.

m_OD_pmf_maroc2fiber

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

private

Definition at line 84 of file ALFA_Raw2Digit.h.

m_OD_pmf_maroc2side

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

private

Definition at line 85 of file ALFA_Raw2Digit.h.

m_ODdigitCollection

ALFA_ODDigitCollection* ALFA_Raw2Digit::m_ODdigitCollection {}

private

Definition at line 114 of file ALFA_Raw2Digit.h.

m_ODdigitObject

ALFA_ODDigit* ALFA_Raw2Digit::m_ODdigitObject {}

private

Definition at line 115 of file ALFA_Raw2Digit.h.

m_ODfiber

int ALFA_Raw2Digit::m_ODfiber {}

private

Definition at line 74 of file ALFA_Raw2Digit.h.

m_ODplate

int ALFA_Raw2Digit::m_ODplate {}

private

Definition at line 73 of file ALFA_Raw2Digit.h.

m_ODpot

int ALFA_Raw2Digit::m_ODpot {}

private

Definition at line 71 of file ALFA_Raw2Digit.h.

m_ODside

int ALFA_Raw2Digit::m_ODside {}

private

Definition at line 72 of file ALFA_Raw2Digit.h.

m_plate

int ALFA_Raw2Digit::m_plate {}

private

Definition at line 68 of file ALFA_Raw2Digit.h.

m_pmf2layer

int ALFA_Raw2Digit::m_pmf2layer[8][23] {}

private

Definition at line 81 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 107 of file ALFA_Raw2Digit.h.

m_pmfID

int ALFA_Raw2Digit::m_pmfID {}

private

Definition at line 60 of file ALFA_Raw2Digit.h.

m_pot

int ALFA_Raw2Digit::m_pot {}

private

Definition at line 66 of file ALFA_Raw2Digit.h.

m_pot_b

bool ALFA_Raw2Digit::m_pot_b[8] {}

private

Definition at line 99 of file ALFA_Raw2Digit.h.

m_rootOutput

TFile* ALFA_Raw2Digit::m_rootOutput {}

private

Definition at line 103 of file ALFA_Raw2Digit.h.

m_rootOutputFileName

TString ALFA_Raw2Digit::m_rootOutputFileName {}

private

Definition at line 104 of file ALFA_Raw2Digit.h.

m_s

std::stringstream ALFA_Raw2Digit::m_s

private

Definition at line 126 of file ALFA_Raw2Digit.h.

m_side

int ALFA_Raw2Digit::m_side {}

private

Definition at line 67 of file ALFA_Raw2Digit.h.

m_strMeasuredDataType

std::string ALFA_Raw2Digit::m_strMeasuredDataType

private

Definition at line 135 of file ALFA_Raw2Digit.h.

m_toFile

std::ofstream ALFA_Raw2Digit::m_toFile

private

Definition at line 138 of file ALFA_Raw2Digit.h.

m_tree

TTree* ALFA_Raw2Digit::m_tree {}

private

Definition at line 102 of file ALFA_Raw2Digit.h.

m_trigger_pattern

std::vector<bool> ALFA_Raw2Digit::m_trigger_pattern

private

Definition at line 98 of file ALFA_Raw2Digit.h.

m_trigger_pattern_b

bool ALFA_Raw2Digit::m_trigger_pattern_b[8][16] {}

private

Definition at line 99 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 131 of file ALFA_Raw2Digit.h.

---

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

• ALFA_Raw2Digit.h
• ALFA_Raw2Digit.cxx