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] {}

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 32 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 50 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 400 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 356 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 491 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 445 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 537 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 629 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 326 of file ALFA_PileUpTool.cxx.

{ return StatusCode::SUCCESS; }

initialize()

StatusCode ALFA_PileUpTool::initialize ( )

finaloverridevirtual

Reimplemented from PileUpToolBase.

Definition at line 94 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 299 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 223 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 117 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 244 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 330 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 342 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 677 of file ALFA_PileUpTool.cxx.

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

Member Data Documentation

m_AmplitudeCut

double ALFA_PileUpTool::m_AmplitudeCut {}

private

Definition at line 105 of file ALFA_PileUpTool.h.

m_cross_talk

double ALFA_PileUpTool::m_cross_talk[8][127] {}

private

Definition at line 110 of file ALFA_PileUpTool.h.

m_digitCollection

ALFA_DigitCollection* ALFA_PileUpTool::m_digitCollection {}

private

Definition at line 95 of file ALFA_PileUpTool.h.

m_E_fib

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

private

Definition at line 86 of file ALFA_PileUpTool.h.

m_E_ODfib

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

private

Definition at line 87 of file ALFA_PileUpTool.h.

m_fibres

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

private

Definition at line 111 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_key_DigitCollection

std::string ALFA_PileUpTool::m_key_DigitCollection

private

Definition at line 92 of file ALFA_PileUpTool.h.

m_key_ODDigitCollection

std::string ALFA_PileUpTool::m_key_ODDigitCollection

private

Definition at line 93 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 103 of file ALFA_PileUpTool.h.

m_meanE_dep

double ALFA_PileUpTool::m_meanE_dep {}

private

Definition at line 101 of file ALFA_PileUpTool.h.

m_meanN_photo

double ALFA_PileUpTool::m_meanN_photo {}

private

Definition at line 102 of file ALFA_PileUpTool.h.

m_mergedALFA_HitList

ALFA_HitCollection* ALFA_PileUpTool::m_mergedALFA_HitList {}

private

Definition at line 107 of file ALFA_PileUpTool.h.

m_mergedALFA_ODHitList

ALFA_ODHitCollection* ALFA_PileUpTool::m_mergedALFA_ODHitList {}

private

Definition at line 108 of file ALFA_PileUpTool.h.

m_mergeSvc

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

private

Definition at line 82 of file ALFA_PileUpTool.h.

m_ODdigitCollection

ALFA_ODDigitCollection* ALFA_PileUpTool::m_ODdigitCollection {}

private

Definition at line 96 of file ALFA_PileUpTool.h.

m_randomStreamName

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

private

Definition at line 84 of file ALFA_PileUpTool.h.

m_randomSvc

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

private

Definition at line 83 of file ALFA_PileUpTool.h.

m_sigma0

double ALFA_PileUpTool::m_sigma0 {}

private

Definition at line 99 of file ALFA_PileUpTool.h.

m_sigma1

double ALFA_PileUpTool::m_sigma1 {}

private

Definition at line 100 of file ALFA_PileUpTool.h.

m_SimHitCollectionName

std::string ALFA_PileUpTool::m_SimHitCollectionName

private

Definition at line 90 of file ALFA_PileUpTool.h.

m_SimODHitCollectionName

std::string ALFA_PileUpTool::m_SimODHitCollectionName

private

Definition at line 91 of file ALFA_PileUpTool.h.

m_stdDev

double ALFA_PileUpTool::m_stdDev {}

private

Definition at line 104 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