ALFA_PileUpTool Class Reference

#include <ALFA_PileUpTool.h>

Inheritance diagram for ALFA_PileUpTool:

Collaboration diagram for ALFA_PileUpTool:

Public Member Functions
	ALFA_PileUpTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual StatusCode	initialize () override final
virtual StatusCode	finalize () override final
virtual StatusCode	prepareEvent (const EventContext &ctx, const unsigned int nInputEvents) override final
	code taken from ZDC; author (highly) probably John Chapman called before the subevts loop. More...
virtual StatusCode	processBunchXing (int bunchXing, SubEventIterator bSubEvents, SubEventIterator eSubEvents) override final
	called for each active bunch-crossing to process current SubEvents bunchXing is in ns More...
virtual StatusCode	mergeEvent (const EventContext &ctx) override final
	return false if not interested in certain xing times (in ns) implemented by default in PileUpToolBase as FirstXing<=bunchXing<=LastXing More...
virtual StatusCode	processAllSubEvents (const EventContext &ctx) override final

Private Member Functions
StatusCode	recordCollection (ServiceHandle< StoreGateSvc > &evtStore, const std::string &key_digitCnt)
void	ALFA_MD_info (const ALFA_HitCollection *)
void	ALFA_MD_info (TimedHitCollection< ALFA_Hit > &)
StatusCode	fill_MD_DigitCollection (CLHEP::HepRandomEngine *)
StatusCode	recordODCollection (ServiceHandle< StoreGateSvc > &evtStore, const std::string &key_ODdigitCnt)
void	ALFA_OD_info (const ALFA_ODHitCollection *)
void	ALFA_OD_info (TimedHitCollection< ALFA_ODHit > &)
StatusCode	fill_OD_DigitCollection (CLHEP::HepRandomEngine *)
StatusCode	XTalk ()

Private Attributes
ServiceHandle< PileUpMergeSvc >	m_mergeSvc {this, "mergeSvc", "PileUpMergeSvc", ""}
ServiceHandle< IAthRNGSvc >	m_randomSvc {this, "RndmSvc", "AthRNGSvc", ""}
Gaudi::Property< std::string >	m_randomStreamName {this, "RandomStreamName", "ALFARndEng", ""}
double	m_E_fib [8][20][64]
double	m_E_ODfib [8][2][3][30]
std::string	m_SimHitCollectionName
std::string	m_SimODHitCollectionName
std::string	m_key_DigitCollection
std::string	m_key_ODDigitCollection
ALFA_DigitCollection *	m_digitCollection
ALFA_ODDigitCollection *	m_ODdigitCollection
double	m_sigma0
double	m_sigma1
double	m_meanE_dep
double	m_meanN_photo
double	m_mean
double	m_stdDev
double	m_AmplitudeCut
ALFA_HitCollection *	m_mergedALFA_HitList
ALFA_ODHitCollection *	m_mergedALFA_ODHitList
double	m_cross_talk [8][127]
int	m_fibres [8][20][64]
std::ifstream	m_fXTalk
std::stringstream	m_s
std::string	m_filename

structors and AlgTool implementation
virtual bool	toProcess (int bunchXing) const override
	the method this base class helps implementing More...
virtual StatusCode	processAllSubEvents (const EventContext &ctx)=0
	dummy implementation to allow compilation while all Digitization packages are migrated to use this new interface. More...
virtual bool	filterPassed () const override
	dummy implementation of passing filter More...
virtual void	resetFilter () override
	dummy implementation of filter reset More...
Gaudi::Property< int >	m_firstXing
Gaudi::Property< int >	m_lastXing
Gaudi::Property< int >	m_vetoPileUpTruthLinks
bool	m_filterPassed {true}

Detailed Description

Definition at line 36 of file ALFA_PileUpTool.h.

Constructor & Destructor Documentation

ALFA_PileUpTool()

ALFA_PileUpTool::ALFA_PileUpTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 49 of file ALFA_PileUpTool.cxx.

                                                          :
  PileUpToolBase(type, name, parent),  
  m_SimHitCollectionName ("ALFA_HitCollection"),
  m_SimODHitCollectionName ("ALFA_ODHitCollection"),
  m_key_DigitCollection ("ALFA_DigitCollection"),
  m_key_ODDigitCollection ("ALFA_ODDigitCollection"),
  m_digitCollection (nullptr),
  m_ODdigitCollection (nullptr),
  m_sigma0          (0.1),
  m_sigma1          (1.0),
  m_meanE_dep       (0.0863), //MeV 
  m_meanN_photo     (4.11), 
  m_mean            (0.0), 
  m_stdDev          (1.0),
  m_AmplitudeCut    (0.24),
  m_mergedALFA_HitList (nullptr),
  m_mergedALFA_ODHitList (nullptr)
{ 
  // Properties  
  declareProperty("SimHitCollection", m_SimHitCollectionName, "Name of the input ALFA Sim Hit Collection of simulated hits");
  declareProperty("SimODHitCollection", m_SimODHitCollectionName, "Name of the input ALFA Sim OD Hit Collection of simulated hits");  
 
  declareProperty("ALFA_DigitCollection", m_key_DigitCollection, "Name of the Collection to hold the output from the ALFA main detector digitization"); 
  declareProperty("ALFA_ODDigitCollection", m_key_ODDigitCollection, "Name of the Collection to hold the output from the ALFA OD digitization");
  
 
//  declareProperty("fillRootTree", m_fillRootTree);  
   
  
  declareProperty("sigma0", m_sigma0); 
  declareProperty("sigma1", m_sigma1);
  declareProperty("meanE_dep", m_meanE_dep);
  declareProperty("meanN_photo", m_meanN_photo);
  declareProperty("mean", m_mean);
  declareProperty("stdDev", m_stdDev);
  declareProperty("AmplitudeCut", m_AmplitudeCut);  
//  declareProperty("fillRootTree", m_fillRootTree);
}

Member Function Documentation

ALFA_MD_info() [1/2]

void ALFA_PileUpTool::ALFA_MD_info ( const ALFA_HitCollection *  ALFA_HitCollection )

private

Definition at line 399 of file ALFA_PileUpTool.cxx.

{
 
// cleaning
  
  for ( int l = 0; l < 8; l++ )
  {
    for ( int i = 0; i < 20; i++ )
    {
      for ( int j = 0; j < 64;  j++ )
    {
    
        m_E_fib[l][i][j] = 0.;
        m_fibres[l][i][j] = 0;
            
    
    } 
    }
  }  
               
// cleaning end 
 
  int fiber, plate, sign, station;
 
  ALFA_HitConstIter it    = ALFA_HitCollection->begin();
  ALFA_HitConstIter itend = ALFA_HitCollection->end();
  
  for (; it != itend; ++it) {
 
    station = it->GetStationNumber();
    plate   = it->GetPlateNumber();
    fiber   = it->GetFiberNumber();
    sign    = it->GetSignFiber();
    
    ATH_MSG_DEBUG ("station=" << station << ", plate= "<< plate << ", fiber=" << fiber << ", sign=" << sign << ", dep energy=" << it->GetEnergyDeposit());    
 
    m_E_fib[station-1][2*(plate-1)+(1-sign)/2][fiber-1] += it->GetEnergyDeposit();
  
 
  }
}

ALFA_MD_info() [2/2]

void ALFA_PileUpTool::ALFA_MD_info ( TimedHitCollection< ALFA_Hit > &  tHitALFA )

private

Definition at line 355 of file ALFA_PileUpTool.cxx.

{ 
 
// cleaning
  
  for ( int l = 0; l < 8; l++ )
  {
    for ( int i = 0; i < 20; i++ )
    {
      for ( int j = 0; j < 64;  j++ )
    {
    
        m_E_fib[l][i][j] = 0.;
        m_fibres[l][i][j] = 0;
    
    } 
    }
  }  
               
// cleaning end 
 
  int fiber, plate, sign, station;
 
  TimedHitCollection<ALFA_Hit>  thpc = tHitALFA;
  TimedHitCollection<ALFA_Hit>::const_iterator i, e, it;
 
  while (thpc.nextDetectorElement(i, e)) for (it = i; it != e; ++it) {
       
    station = (*it)->GetStationNumber();
    plate   = (*it)->GetPlateNumber();
    fiber   = (*it)->GetFiberNumber();
    sign    = (*it)->GetSignFiber();
    
    ATH_MSG_DEBUG ("station=" << station << ", plate= "<< plate << ", fiber=" << fiber << ", sign=" << sign << ", dep energy=" << ((*it)->GetEnergyDeposit()));    
 
    m_E_fib[station-1][2*(plate-1)+(1-sign)/2][fiber-1] += ((*it)->GetEnergyDeposit());
      
  }  
  
}

ALFA_OD_info() [1/2]

void ALFA_PileUpTool::ALFA_OD_info ( const ALFA_ODHitCollection *  ALFA_ODHitCollection )

private

Definition at line 490 of file ALFA_PileUpTool.cxx.

{
  
  // cleaning
  
  for ( int l = 0; l < 8; l++ ){
    for ( int i = 0; i < 2; i++ ){
      for ( int j = 0; j < 3; j++ ){
        for ( int k = 0; k < 30; k++ ){
      
            m_E_ODfib[l][i][j][k] = 0.; 
      
      }
      }
    }
  }               
  
// cleaning end
 
  int fiber, plate, sign, side, station; 
  
  
  ALFA_ODHitConstIter it    = ALFA_ODHitCollection->begin();
  ALFA_ODHitConstIter itend = ALFA_ODHitCollection->end();
 
  for (; it != itend; ++it) {
 
     
    station = it->GetStationNumber();
    side    = it->GetODSide();    
    plate   = it->GetPlateNumber();
    fiber   = it->GetFiberNumber();
    sign    = it->GetSignFiber();
 
       
    ATH_MSG_DEBUG ("station=" << station << ", side=" << side << ", plate= "<< plate << ", fiber=" << fiber << ", sign=" << sign << ", dep energy=" << it->GetEnergyDeposit());
    
    if (sign==0) m_E_ODfib[station-1][side-1][plate-1][fiber+15] += it->GetEnergyDeposit();
    else m_E_ODfib[station-1][side-1][plate-1][fiber] += it->GetEnergyDeposit();
    
  }
}

ALFA_OD_info() [2/2]

void ALFA_PileUpTool::ALFA_OD_info ( TimedHitCollection< ALFA_ODHit > &  tODHitALFA )

private

Definition at line 444 of file ALFA_PileUpTool.cxx.

{
 
// cleaning
  
  for ( int l = 0; l < 8; l++ ){
    for ( int i = 0; i < 2; i++ ){
      for ( int j = 0; j < 3; j++ ){
        for ( int k = 0; k < 30; k++ ){
      
            m_E_ODfib[l][i][j][k] = 0.; 
      
      }
      }
    }
  }               
  
// cleaning end
 
  int fiber, plate, sign, side, station; 
 
 
  TimedHitCollection<ALFA_ODHit>  thpc = tODHitALFA;
  TimedHitCollection<ALFA_ODHit>::const_iterator i, e, it;
  
  while (thpc.nextDetectorElement(i, e)) for (it = i; it != e; ++it) {
     
    station = (*it)->GetStationNumber();
    side    = (*it)->GetODSide();    
    plate   = (*it)->GetPlateNumber();
    fiber   = (*it)->GetFiberNumber();
    sign    = (*it)->GetSignFiber();
 
       
    ATH_MSG_DEBUG ("station=" << station << ", side=" << side << ", plate= "<< plate << ", fiber=" << fiber << ", sign=" << sign << ", dep energy=" << ((*it)->GetEnergyDeposit()));
 
    if (sign==0) m_E_ODfib[station-1][side-1][plate-1][fiber+15] += ((*it)->GetEnergyDeposit());
    else m_E_ODfib[station-1][side-1][plate-1][fiber] += ((*it)->GetEnergyDeposit());
    
  }  
 
}

fill_MD_DigitCollection()

StatusCode ALFA_PileUpTool::fill_MD_DigitCollection ( CLHEP::HepRandomEngine *  rndEngine )

private

Definition at line 536 of file ALFA_PileUpTool.cxx.

                                                                                   {
    
  ATH_MSG_DEBUG(" ALFA_PileUpTool::fill_MD_DigitCollection()");
  
  double N_photo = 0.;
  //double N_CTphoto = 0.;
  //double sigma   = 0.;
  double amplitude = 0.; 
  //double amplitude_CT = 0.; 
 
  
  for ( int l = 0; l < 8; l++ ) 
  {
    for ( int i = 0; i < 20; i++ ) 
    {
      for ( int j = 0; j < 64;  j++ ) 
      { 
      
        N_photo   = CLHEP::RandPoisson::shoot(rndEngine,m_E_fib[l][i][j]*m_meanN_photo/m_meanE_dep);
        
        //sigma     = m_sigma1 * sqrt(N_photo);
        //double noise_1 = sigma * CLHEP::RandGaussZiggurat::shoot (rndEngine, m_mean, m_stdDev);
        //double noise_2 = m_sigma0 * CLHEP::RandGaussZiggurat::shoot (rndEngine, m_mean, m_stdDev);
       
        //amplitude =  N_photo 
        //             + noise_1 
        //             + noise_2;
            
        amplitude = CLHEP::RandGaussZiggurat::shoot(rndEngine, N_photo, sqrt(pow(m_sigma0,2)+N_photo*pow(m_sigma1,2)));
 
        //N_photo_CT   = CLHEP::RandPoisson::shoot(rndEngine,0.08*m_E_fib[l][i][j]*m_meanN_photo/m_meanE_dep);      
        //amplitude_CT = CLHEP::RandGaussQ::shoot (rndEngine, N_photo_CT, sqrt(pow(m_sigma0,2)+N_photo_CT*pow(m_sigma1,2)));
                 
                  
        if (amplitude >= m_AmplitudeCut )
        {
          ATH_MSG_DEBUG(" ALFA_Digitization::fillDigitCollection, amplitude " << amplitude);
          ATH_MSG_DEBUG(" station = " << l << ", plate= " << i << ", fiber=" << j );
          m_digitCollection->push_back(new ALFA_Digit(l,i,j));
        m_fibres[l][i][j] = m_fibres[l][i][j] + 1;
        }
      
      }
    }
  }
  
  double rand_fib;
  
    for ( int l = 0; l < 8; l++ ) 
    {
        for ( int i = 0; i < 20; i++ ) 
        {
            for ( int j = 0; j < 64;  j++ ) 
            {   
                if ( m_fibres[l][i][j] > 0 )
                {
                    for (int f = j+1; f < 64; f++)
                    {   
                        rand_fib  = CLHEP::RandFlat::shoot(rndEngine,0.,1.);
                        
                        if ( m_fibres[l][i][f] == 0)
                        {
                            if (rand_fib <= m_cross_talk[l][63+f-j]) 
                            {   
                                m_digitCollection->push_back(new ALFA_Digit(l,i,f));
                            }
                        }
                    }   
                                    
                    for (int f = j-1; f > -1; f--)
                    {   
                        rand_fib = CLHEP::RandFlat::shoot(rndEngine,0.,1.); 
                            
                        if ( m_fibres[l][i][f] == 0)
                        {
                            if (rand_fib <= m_cross_talk[l][63-(j-f)]) 
                            {   
                                m_digitCollection->push_back(new ALFA_Digit(l,i,f));
                            }
                        }
                    }
                }
            }
        }
    }
  return StatusCode::SUCCESS;
   
}  

fill_OD_DigitCollection()

StatusCode ALFA_PileUpTool::fill_OD_DigitCollection ( CLHEP::HepRandomEngine *  rndEngine )

private

Definition at line 628 of file ALFA_PileUpTool.cxx.

                                                                                   {
    
  ATH_MSG_DEBUG(" ALFA_PileUpTool::fill_OD_DigitCollection ");
  
  double N_photo = 0.;
  double sigma   = 0.;
  double amplitude = 0.;
  
  for ( int l = 0; l < 8; l++ ) 
  {
    for ( int i = 0; i < 2; i++ ) 
    {
      for ( int j = 0; j < 3; j++ ) 
    {
        for (int k = 0; k < 30; k++) 
      { 
      
        N_photo   = CLHEP::RandPoisson::shoot(rndEngine,m_E_ODfib[l][i][j][k]*m_meanN_photo/m_meanE_dep);
        sigma     = m_sigma1 * sqrt(N_photo);
         
          double noise_1 = sigma * CLHEP::RandGaussZiggurat::shoot (rndEngine, m_mean, m_stdDev);
          double noise_2 = m_sigma0 * CLHEP::RandGaussZiggurat::shoot (rndEngine, m_mean, m_stdDev);
         
         
          amplitude =  N_photo 
                   + noise_1 
                 + noise_2;
      
            
        if (amplitude >= m_AmplitudeCut )
        {
          ATH_MSG_DEBUG(" ALFA_Digitization::fill_OD_DigitCollection, amplitude " << amplitude);
          ATH_MSG_DEBUG(" station = " << l << ", side = " << i << ", plate = " << j << ", fiber = " << k);
 
          m_ODdigitCollection->push_back(new ALFA_ODDigit(l,i,j,k));
        }
 
      
        }
      }
    }
  }
    
  return StatusCode::SUCCESS;
}

filterPassed()

virtual bool PileUpToolBase::filterPassed ( ) const

inlineoverridevirtualinherited

dummy implementation of passing filter

Definition at line 49 of file PileUpToolBase.h.

{ return m_filterPassed; }

finalize()

StatusCode ALFA_PileUpTool::finalize ( )

finaloverridevirtual

Definition at line 325 of file ALFA_PileUpTool.cxx.

{ return StatusCode::SUCCESS; }

initialize()

StatusCode ALFA_PileUpTool::initialize ( )

finaloverridevirtual

Reimplemented from PileUpToolBase.

Definition at line 93 of file ALFA_PileUpTool.cxx.

                                      {
  
  ATH_MSG_DEBUG ("ALFA_digiti initialize() called"); 
  
  ATH_MSG_DEBUG (" Pedestal: " << m_AmplitudeCut);
  
  ATH_CHECK (m_randomSvc.retrieve());
  ATH_MSG_DEBUG ( "Retrieved RandomNumber Service" );
 
  ATH_CHECK (m_mergeSvc.retrieve());
  ATH_MSG_DEBUG("Retrieved PileUpMergeSvc");
 
  m_mergedALFA_HitList = new ALFA_HitCollection();
  m_mergedALFA_ODHitList = new ALFA_ODHitCollection();
  
  return StatusCode::SUCCESS;
}

mergeEvent()

StatusCode ALFA_PileUpTool::mergeEvent ( const EventContext &  ctx )

finaloverridevirtual

return false if not interested in certain xing times (in ns) implemented by default in PileUpToolBase as FirstXing<=bunchXing<=LastXing

called at the end of the subevts loop. Not (necessarily) able to access SubEvents

Definition at line 298 of file ALFA_PileUpTool.cxx.

                                                             {
 
 
  ALFA_MD_info(m_mergedALFA_HitList);
 
  // Prepare RNG Service
  ATHRNG::RNGWrapper* rngWrapper = m_randomSvc->getEngine(this, m_randomStreamName);
  rngWrapper->setSeed( m_randomStreamName, ctx );
  CLHEP::HepRandomEngine* rngEngine = rngWrapper->getEngine(ctx);
  StatusCode sc = fill_MD_DigitCollection(rngEngine);
  if  (sc.isFailure()) {
    ATH_MSG_WARNING ("ALFA_PileUpTool::fill_MD_DigitCollection() failed");  
    return StatusCode::SUCCESS;     
  }
    
  
  ALFA_OD_info(m_mergedALFA_ODHitList);
  
  sc = fill_OD_DigitCollection(rngEngine);
  if  (sc.isFailure()) {
    ATH_MSG_WARNING ("ALFA_PileUpTool::fill_OD_DigitCollection() failed");  
    return StatusCode::SUCCESS;     
  }   
  
  return StatusCode::SUCCESS;
}

prepareEvent()

StatusCode ALFA_PileUpTool::prepareEvent ( const EventContext &  ctx, const unsigned int  nInputEvents  )

finaloverridevirtual

code taken from ZDC; author (highly) probably John Chapman called before the subevts loop.

Not (necessarily) able to access SubEvents

Definition at line 222 of file ALFA_PileUpTool.cxx.

                                                                                              {
 
  ATH_MSG_DEBUG ( "ALFA_PileUpTool::prepareEvent() called for " << nInputEvents << " input events" );
 
  StatusCode sc = recordCollection(this->evtStore(), m_key_DigitCollection);
 
  if (sc.isFailure()) { ATH_MSG_FATAL ( " ALFA_PileUpTool::prepareEvent(): Could not record the empty digit container in StoreGate " ); return sc; }
  else                { ATH_MSG_DEBUG ( " ALFA_PileUpTool::prepareEvent(): Digit container is recorded in StoreGate " ); }
  
  sc = recordODCollection(this->evtStore(), m_key_ODDigitCollection);
 
  if (sc.isFailure()) { ATH_MSG_FATAL ( " ALFA_PileUpTool::prepareEvent(): Could not record the empty digit OD container in StoreGate " ); return sc; }
  else                { ATH_MSG_DEBUG ( " ALFA_PileUpTool::prepareEvent(): Digit OD container is recorded in StoreGate " ); }
  
  m_mergedALFA_HitList->clear();
  m_mergedALFA_ODHitList->clear();
  
  return StatusCode::SUCCESS;
}

processAllSubEvents() [1/2]

StatusCode ALFA_PileUpTool::processAllSubEvents ( const EventContext &  ctx )

finaloverridevirtual

Reimplemented from PileUpToolBase.

Definition at line 116 of file ALFA_PileUpTool.cxx.

                                                                       {
 
  ATH_MSG_DEBUG ("ALFA_PileUpTool::processAllSubEvents()");
 
  using TimedALFAHitCollList = PileUpMergeSvc::TimedList<ALFA_HitCollection>::type;
  using TimedALFAODHitCollList = PileUpMergeSvc::TimedList<ALFA_ODHitCollection>::type;
  
  TimedALFAHitCollList tHitCollList;
  TimedALFAODHitCollList tODHitCollList;
  
  
  if (m_mergeSvc->retrieveSubEvtsData(m_SimHitCollectionName, tHitCollList).isFailure()) { 
 
    ATH_MSG_FATAL ( "Could not fill TimedALFAHitCollList" ); return StatusCode::FAILURE; 
  }
  else { ATH_MSG_DEBUG ( "Retrieved TimedALFAHitCollList" ); }
  
  ATH_MSG_DEBUG ( "PileUp: Merge " << tHitCollList.size() << " ALFA_HitCollection with key" << m_SimHitCollectionName );
  
  
  if (m_mergeSvc->retrieveSubEvtsData(m_SimODHitCollectionName, tODHitCollList).isFailure()) { 
 
    ATH_MSG_FATAL ( "Could not fill TimedALFAODHitCollList" ); return StatusCode::FAILURE; 
  }
  else { ATH_MSG_DEBUG ( "Retrieved TimedALFAODHitCollList" ); }
  
  ATH_MSG_DEBUG ( "PileUp: Merge " << tODHitCollList.size() << " ALFA_ODHitCollection with key " << m_SimODHitCollectionName );
  
  
  TimedHitCollection<ALFA_Hit> tALFAhit;
 
  TimedALFAHitCollList::iterator iHitColl (tHitCollList.begin());
  TimedALFAHitCollList::iterator eHitColl (tHitCollList.end());
  
  while (iHitColl != eHitColl) {
    
    const ALFA_HitCollection* tmpColl(iHitColl->second);
    
    tALFAhit.insert(iHitColl->first, tmpColl);
    
    ATH_MSG_DEBUG ( " ALFA_HitCollection found with " << tmpColl->size() << " hits " << iHitColl->first );
    
    ++iHitColl;
  }
  
  TimedHitCollection<ALFA_ODHit> tALFAODhit;
 
  TimedALFAODHitCollList::iterator iODHitColl  (tODHitCollList.begin());
  TimedALFAODHitCollList::iterator eODHitColl  (tODHitCollList.end());
  
  while (iODHitColl != eODHitColl) {
    
    const ALFA_ODHitCollection* tmpColl(iODHitColl->second);
    
    tALFAODhit.insert(iODHitColl->first, tmpColl);
    
    ATH_MSG_DEBUG ( " ALFA_ODHitCollection found with " << tmpColl->size() << " hits " << iODHitColl->first );
    
    ++iODHitColl;
  }
 
  // preparing containers 
  
  if (recordCollection(this->evtStore(), m_key_DigitCollection).isFailure()) { 
    
    ATH_MSG_FATAL ( " ALFA_PileUpTool::processAllSubEvents(): Could not record the empty ALFA digit container in StoreGate " ); return StatusCode::FAILURE; 
  }
  else { ATH_MSG_DEBUG ( " ALFA_PileUpTool::processAllSubEvents(): ALFA Digit container is recorded in StoreGate " ); }
  
  
  if (recordODCollection(this->evtStore(), m_key_ODDigitCollection).isFailure()) { 
    
    ATH_MSG_FATAL ( " ALFA_PileUpTool::processAllSubEvents(): Could not record the empty ALFA OD digit container in StoreGate " ); return StatusCode::FAILURE; 
  }
  else { ATH_MSG_DEBUG ( " ALFA_PileUpTool::processAllSubEvents(): ALFA OD Digit container is recorded in StoreGate " ); }  
 
 
  // filling containers 
  
  ALFA_MD_info(tALFAhit);
 
  // Prepare RNG Service
  ATHRNG::RNGWrapper* rngWrapper = m_randomSvc->getEngine(this, m_randomStreamName);
  rngWrapper->setSeed( m_randomStreamName, ctx );
  CLHEP::HepRandomEngine* rngEngine = rngWrapper->getEngine(ctx);
  StatusCode sc = fill_MD_DigitCollection(rngEngine);
  if  (sc.isFailure()) {
    ATH_MSG_WARNING ("ALFA_PileUpTool::fill_MD_DigitCollection() failed");  
    return StatusCode::SUCCESS;     
  }
  
    
  ALFA_OD_info(tALFAODhit);
  
  sc = fill_OD_DigitCollection(rngEngine);
  if  (sc.isFailure()) {
    ATH_MSG_WARNING ("ALFA_PileUpTool::fill_OD_DigitCollection() failed");  
    return StatusCode::SUCCESS;     
  }  
  
  return StatusCode::SUCCESS;
}

processAllSubEvents() [2/2]

virtual StatusCode IPileUpTool::processAllSubEvents

inherited

dummy implementation to allow compilation while all Digitization packages are migrated to use this new interface.

processBunchXing()

StatusCode ALFA_PileUpTool::processBunchXing ( int  bunchXing, SubEventIterator  bSubEvents, SubEventIterator  eSubEvents  )

finaloverridevirtual

called for each active bunch-crossing to process current SubEvents bunchXing is in ns

Reimplemented from PileUpToolBase.

Definition at line 243 of file ALFA_PileUpTool.cxx.

                                                                              {
  ATH_MSG_DEBUG ( "ALFA_PileUpTool::processBunchXing() " << bunchXing );
  SubEventIterator iEvt = bSubEvents;
  for (; iEvt!=eSubEvents; ++iEvt) {
    StoreGateSvc& seStore = *iEvt->ptr()->evtStore();
    //PileUpTimeEventIndex thisEventIndex = PileUpTimeEventIndex(static_cast<int>(iEvt->time()),iEvt->index()); not used 
    ATH_MSG_VERBOSE("SubEvt StoreGate " << seStore.name() << " :"
                    << " bunch crossing : " << bunchXing
                    << " time offset : " << iEvt->time()
                    << " event number : " << iEvt->ptr()->eventNumber()
                    << " run number : " << iEvt->ptr()->runNumber()
                    );
  
    
    const ALFA_HitCollection* tmpHitColl = nullptr;
   
    if (!seStore.retrieve(tmpHitColl, m_SimHitCollectionName).isSuccess()) {
      
      ATH_MSG_ERROR ("SubEvent ALFA_HitCollection not found in StoreGate " << seStore.name() );
      
      return StatusCode::FAILURE;
    }
    
    ATH_MSG_DEBUG ("SubEvent, ALFA_HitCollection found with " << tmpHitColl->size() << " hits");
    
    ALFA_HitCollection::const_iterator iHitColl = tmpHitColl->begin();
    ALFA_HitCollection::const_iterator eHitColl = tmpHitColl->end();
   
    for (; iHitColl!=eHitColl; ++iHitColl) m_mergedALFA_HitList->push_back((*iHitColl));
    
    const ALFA_ODHitCollection* tmpODHitColl = nullptr;
   
    if (!seStore.retrieve(tmpODHitColl, m_SimODHitCollectionName).isSuccess()) {
      
      ATH_MSG_ERROR ( "SubEvent ALFA_ODHitCollection not found in StoreGate " << seStore.name() );
      
      return StatusCode::FAILURE;
    }
    
    ATH_MSG_DEBUG ( "ALFA_ODHitCollection found with " << tmpODHitColl->size() << " hits" );
    
    ALFA_ODHitCollection::const_iterator iODHitColl = tmpODHitColl->begin();
    ALFA_ODHitCollection::const_iterator eODHitColl = tmpODHitColl->end();
   
    for (; iODHitColl!=eODHitColl; ++iODHitColl) m_mergedALFA_ODHitList->push_back((*iODHitColl));
  }
  
  return StatusCode::SUCCESS;
}

recordCollection()

StatusCode ALFA_PileUpTool::recordCollection ( ServiceHandle< StoreGateSvc > &  evtStore, const std::string &  key_digitCnt  )

private

Definition at line 329 of file ALFA_PileUpTool.cxx.

{
  ATH_MSG_DEBUG ("ALFA_Digitization::recordCollection"); 
    
  m_digitCollection = new ALFA_DigitCollection();
  
  StatusCode sc = evtStore->record(m_digitCollection, key_digitCnt);
 
  return sc;
}

recordODCollection()

StatusCode ALFA_PileUpTool::recordODCollection ( ServiceHandle< StoreGateSvc > &  evtStore, const std::string &  key_ODdigitCnt  )

private

Definition at line 341 of file ALFA_PileUpTool.cxx.

                                                                                                                     {
 
  ATH_MSG_DEBUG ("ALFA_Digitization::recordODCollection");  
    
  m_ODdigitCollection = new ALFA_ODDigitCollection();
  
  StatusCode sc = evtStore->record(m_ODdigitCollection, key_ODdigitCnt);
  
  return sc;
}

resetFilter()

virtual void PileUpToolBase::resetFilter ( )

inlineoverridevirtualinherited

dummy implementation of filter reset

Reimplemented in MergeTruthJetsTool.

Definition at line 51 of file PileUpToolBase.h.

{ m_filterPassed=true; }

toProcess()

virtual bool PileUpToolBase::toProcess ( int  bunchXing ) const

inlineoverridevirtualinherited

the method this base class helps implementing

Reimplemented in MergeHijingParsTool, and MergeTrackRecordCollTool.

Definition at line 32 of file PileUpToolBase.h.

                                                       {
    //closed interval [m_firstXing,m_lastXing]
    return !((m_firstXing > bunchXing) || (bunchXing > m_lastXing));
  }

XTalk()

StatusCode ALFA_PileUpTool::XTalk ( )

private

Definition at line 675 of file ALFA_PileUpTool.cxx.

  { 
    for (unsigned int j=0; j<8; j++){
        
        m_s.str("");
            
        m_filename = "Xtalk_station";
        m_s << j+1;
        m_filename += m_s.str();
        m_filename += ".txt";
        
        ATH_MSG_DEBUG("file name " << m_filename.c_str() );
        
        std::string filePath = PathResolver::find_file(m_filename,"DATAPATH", PathResolver::RecursiveSearch);
        
        if(filePath.length() == 0)
        {
            ATH_MSG_FATAL(" XTalk file " <<  m_filename.c_str() << " not found in Datapath");
            throw std::runtime_error("FATAL: mapping MD maroc-mapmt not found in Datapath.");
        }
        
        else
        {
            ATH_MSG_DEBUG("the XTALK file \"" <<  m_filename.c_str() << "\" found in Datapath");
            ATH_MSG_DEBUG("filePath =  " << filePath.c_str() );
        }
 
        m_fXTalk.open(filePath.c_str());
        
        if (m_fXTalk.is_open())
        {
            for (unsigned int i=0;i<127;i++)
            {               
                m_fXTalk >> m_cross_talk[j][i];
 
                ATH_MSG_DEBUG("cross_talk[" << j << "][" << i << "] = " << m_cross_talk[j][i]);
 
            }
            m_fXTalk.close();
        }
        else
        {
            ATH_MSG_ERROR("the file " << m_filename.c_str() << " was not open");
              return StatusCode::FAILURE;
        }   
    }       
    return StatusCode::SUCCESS;
  }

Member Data Documentation

m_AmplitudeCut

double ALFA_PileUpTool::m_AmplitudeCut

private

Definition at line 149 of file ALFA_PileUpTool.h.

m_cross_talk

double ALFA_PileUpTool::m_cross_talk[8][127]

private

Definition at line 162 of file ALFA_PileUpTool.h.

m_digitCollection

ALFA_DigitCollection* ALFA_PileUpTool::m_digitCollection

private

Definition at line 130 of file ALFA_PileUpTool.h.

m_E_fib

double ALFA_PileUpTool::m_E_fib[8][20][64]

private

Definition at line 115 of file ALFA_PileUpTool.h.

m_E_ODfib

double ALFA_PileUpTool::m_E_ODfib[8][2][3][30]

private

Definition at line 116 of file ALFA_PileUpTool.h.

m_fibres

int ALFA_PileUpTool::m_fibres[8][20][64]

private

Definition at line 163 of file ALFA_PileUpTool.h.

m_filename

std::string ALFA_PileUpTool::m_filename

private

Definition at line 167 of file ALFA_PileUpTool.h.

m_filterPassed

bool PileUpToolBase::m_filterPassed {true}

protectedinherited

Definition at line 60 of file PileUpToolBase.h.

m_firstXing

Gaudi::Property<int> PileUpToolBase::m_firstXing

protectedinherited

Initial value:

{this, "FirstXing", -999,
      "First bunch-crossing in which det is live"}

Definition at line 54 of file PileUpToolBase.h.

m_fXTalk

std::ifstream ALFA_PileUpTool::m_fXTalk

private

Definition at line 165 of file ALFA_PileUpTool.h.

m_key_DigitCollection

std::string ALFA_PileUpTool::m_key_DigitCollection

private

Definition at line 126 of file ALFA_PileUpTool.h.

m_key_ODDigitCollection

std::string ALFA_PileUpTool::m_key_ODDigitCollection

private

Definition at line 127 of file ALFA_PileUpTool.h.

m_lastXing

Gaudi::Property<int> PileUpToolBase::m_lastXing

protectedinherited

Initial value:

{this, "LastXing", 999,
      "Last bunch-crossing in which det is live"}

Definition at line 56 of file PileUpToolBase.h.

m_mean

double ALFA_PileUpTool::m_mean

private

Definition at line 147 of file ALFA_PileUpTool.h.

m_meanE_dep

double ALFA_PileUpTool::m_meanE_dep

private

Definition at line 145 of file ALFA_PileUpTool.h.

m_meanN_photo

double ALFA_PileUpTool::m_meanN_photo

private

Definition at line 146 of file ALFA_PileUpTool.h.

m_mergedALFA_HitList

ALFA_HitCollection* ALFA_PileUpTool::m_mergedALFA_HitList

private

Definition at line 159 of file ALFA_PileUpTool.h.

m_mergedALFA_ODHitList

ALFA_ODHitCollection* ALFA_PileUpTool::m_mergedALFA_ODHitList

private

Definition at line 160 of file ALFA_PileUpTool.h.

m_mergeSvc

ServiceHandle<PileUpMergeSvc> ALFA_PileUpTool::m_mergeSvc {this, "mergeSvc", "PileUpMergeSvc", ""}

private

Definition at line 111 of file ALFA_PileUpTool.h.

m_ODdigitCollection

ALFA_ODDigitCollection* ALFA_PileUpTool::m_ODdigitCollection

private

Definition at line 132 of file ALFA_PileUpTool.h.

m_randomStreamName

Gaudi::Property<std::string> ALFA_PileUpTool::m_randomStreamName {this, "RandomStreamName", "ALFARndEng", ""}

private

Definition at line 113 of file ALFA_PileUpTool.h.

m_randomSvc

ServiceHandle<IAthRNGSvc> ALFA_PileUpTool::m_randomSvc {this, "RndmSvc", "AthRNGSvc", ""}

private

Definition at line 112 of file ALFA_PileUpTool.h.

m_s

std::stringstream ALFA_PileUpTool::m_s

private

Definition at line 166 of file ALFA_PileUpTool.h.

m_sigma0

double ALFA_PileUpTool::m_sigma0

private

Definition at line 143 of file ALFA_PileUpTool.h.

m_sigma1

double ALFA_PileUpTool::m_sigma1

private

Definition at line 144 of file ALFA_PileUpTool.h.

m_SimHitCollectionName

std::string ALFA_PileUpTool::m_SimHitCollectionName

private

Definition at line 124 of file ALFA_PileUpTool.h.

m_SimODHitCollectionName

std::string ALFA_PileUpTool::m_SimODHitCollectionName

private

Definition at line 125 of file ALFA_PileUpTool.h.

m_stdDev

double ALFA_PileUpTool::m_stdDev

private

Definition at line 148 of file ALFA_PileUpTool.h.

m_vetoPileUpTruthLinks

Gaudi::Property<int> PileUpToolBase::m_vetoPileUpTruthLinks

protectedinherited

Initial value:

{this, "VetoPileUpTruthLinks", true,
      "Ignore links to suppressed pile-up truth"}

Definition at line 58 of file PileUpToolBase.h.

---

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

• ALFA_PileUpTool.h
• ALFA_PileUpTool.cxx