ATLAS Offline Software
Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Static Private Member Functions | Private Attributes | Static Private Attributes | List of all members
ZdcRecChannelTool Class Reference

#include <ZdcRecChannelTool.h>

Inheritance diagram for ZdcRecChannelTool:
Collaboration diagram for ZdcRecChannelTool:

Public Member Functions

 ZdcRecChannelTool (const std::string &type, const std::string &name, const IInterface *parent)
 
virtual ~ZdcRecChannelTool ()
 
virtual StatusCode initialize ()
 
virtual StatusCode finalize ()
 
int makeRawFromDigits (const ZdcDigitsCollection &data_collection, ZdcRawChannelCollection &raw_collection)
 
int getCalibration (const ZdcDigitsCollection &data_collection, ZdcRawChannelCollection &raw_collection)
 
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 sysInitialize () override
 Perform system initialization for an algorithm. 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
 

Static Public Member Functions

static const InterfaceID & interfaceID ()
 AlgTool InterfaceID. More...
 

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

int getTimingCFD (const Identifier &id, const std::vector< std::vector< int > > &wfm)
 
int getTimingSinc (const Identifier &id, const std::vector< std::vector< int > > &wfm)
 
int getTimingSinc2 (const Identifier &id, const std::vector< std::vector< int > > &wfm)
 
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...
 

Static Private Member Functions

static double fx (double x0, void *params)
 
static double fxCallback (double d, void *v)
 

Private Attributes

gsl_interp_accel * m_interp_acc {}
 
gsl_spline * m_spline {}
 
std::vector< float > m_cfd_result
 
int m_bwl_time_resolution
 
std::vector< float > m_wfm_bwl
 
std::vector< float > m_bwl_vpeak
 
std::vector< float > m_bwl_tpeak
 
std::vector< float > m_bwl_vpeak2
 
std::vector< float > m_bwl_tpeak2
 
unsigned int m_nsamples
 
float m_sample_time
 
int m_delta
 
int m_zeroSupress
 
float m_cfd_fraction
 
float m_cfd_delay
 
const ZdcID * m_zdcId {}
 
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
 

Static Private Attributes

static const int s_FADC_SATURATION = 1022
 
static const int s_CFD_FRACTION = -3
 
static const int s_SAMPLING_TIME = 25000
 
static const int s_HALF_SAMPLING_TIME = 12500
 
static const int s_GAIN_RATIO = 1
 

Detailed Description

Definition at line 71 of file ZdcRecChannelTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ ZdcRecChannelTool()

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

Definition at line 52 of file ZdcRecChannelTool.cxx.

54  :
55  AthAlgTool(type, name, parent),
56  m_bwl_time_resolution(100),
57  m_nsamples(7),
58  m_sample_time (25.),
59  m_cfd_fraction (0.35),
60  m_cfd_delay (20.)
61 
62 {
63  //Declare properties here...
64 
65  declareInterface<ZdcRecChannelTool>(this);
66 
67  declareProperty("ZeroSuppress", m_zeroSupress = 0,
68  "Supress channels with only 0");
69 
70  declareProperty("DeltaPeak", m_delta = 5,
71  "Minimum difference between min and max to be considered a signal");
72 
73 
74 }

◆ ~ZdcRecChannelTool()

ZdcRecChannelTool::~ZdcRecChannelTool ( )
virtual

Definition at line 78 of file ZdcRecChannelTool.cxx.

79 {
80 }

Member Function Documentation

◆ declareGaudiProperty() [1/4]

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

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

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

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

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

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

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

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

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

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

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation 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.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

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

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation 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.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

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

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

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

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation 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.

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

◆ declareProperty() [5/6]

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

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation 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.

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

◆ declareProperty() [6/6]

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

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )
inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::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

◆ finalize()

StatusCode ZdcRecChannelTool::finalize ( )
virtual

Definition at line 138 of file ZdcRecChannelTool.cxx.

139 {
140  msg(MSG::INFO) << "Finalizing " << name() << endmsg;
141 
142  gsl_spline_free (m_spline);
143  gsl_interp_accel_free (m_interp_acc);
144 
145  return StatusCode::SUCCESS;
146 }

◆ fx()

double ZdcRecChannelTool::fx ( double  x0,
void *  params 
)
staticprivate

Definition at line 617 of file ZdcRecChannelTool.cxx.

618  {
619  ZdcRecChannelTool* p = static_cast<ZdcRecChannelTool*>(params);
620  //std::cout << " MMMMMMMMMMMMAAAAAAAAA " << p->m_cfd_delay << std::endl;
621  return gsl_spline_eval (p->m_spline, x0, p->m_interp_acc);
622  }

◆ fxCallback()

static double ZdcRecChannelTool::fxCallback ( double  d,
void *  v 
)
inlinestaticprivate

Definition at line 105 of file ZdcRecChannelTool.h.

106  {
107  CCallbackHolder* h = static_cast<CCallbackHolder*>(v);
108  return h->cls->fx(d, h->cls);
109  }

◆ getCalibration()

int ZdcRecChannelTool::getCalibration ( const ZdcDigitsCollection &  data_collection,
ZdcRawChannelCollection &  raw_collection 
)

Definition at line 457 of file ZdcRecChannelTool.cxx.

460 {
461  int ncha = 0;
462  for (const ZdcDigits* digits_p : mydata) {
463  digits_p->identify();
464  ncha++;
465  }
466  int ncal = ChannelCollection.size();
467  msg(MSG::DEBUG) << " Zdc ---> Calibration for " << ncal << " channels " << endmsg;
468  return ncha;
469 }

◆ getTimingCFD()

int ZdcRecChannelTool::getTimingCFD ( const Identifier &  id,
const std::vector< std::vector< int > > &  wfm 
)
private

Definition at line 473 of file ZdcRecChannelTool.cxx.

474 {
475 
476  unsigned int i = 0;
477  //unsigned int j = 0;
478  unsigned int nsamples_local = 14;
479  double v[14];
480  double t[14];
481  double td[14];
482  double vd[14];
483  double vf[14];
484 
485  std::vector<float> vaux;
486  unsigned int tmin = 0;
487  unsigned int tmax = 0;
488 
489  double x_lo = 0.0;
490  double x_hi = 200.0;
491  float r = 0;
492  int status = 0;
493 
494  int iterations = 0;
495  int max_iterations = 100;
496 
497  // Thread-safe according to https://www.gnu.org/software/gsl/doc/html/roots.html
498  const gsl_root_fsolver_type *T ATLAS_THREAD_SAFE = gsl_root_fsolver_brent;
499  gsl_root_fsolver *s = gsl_root_fsolver_alloc (T);
500  gsl_function F;
501 
502  CCallbackHolder cc{};
503 
504 
505  std::vector<int> y;
506 
507  std::vector<std::vector<int> >::const_iterator vit;
508  std::vector<int>::iterator it;
509 
510  std:: vector<float> result;
511 
512  m_cfd_result.clear();
513 
514  msg(MSG::DEBUG) << "--> ZDC : ZdcRecChannelTool::getTimingCFD: Id " << id.getString() ;
515 
516  for (i=0;i<14;i++) {
517  t[i] = 0.;
518  td[i] = 0.;
519  v[i] = 0.;
520  vd[i] = 0.;
521  vf[i] = 0.;
522  }
523 
524  for (vit = wfm.begin(); vit<wfm.end(); ++vit) {
525  y = *vit;
526  i = 4;
527  for (it=y.begin();it != y.end();++it) {
528  v[i] = *it;
529  i++;
530  }
531  for (i=0;i<nsamples_local;i++) {
532  t[i] = m_sample_time*i ; //TODO: Do it only once somewhere else
533  td[i] = t[i]-m_cfd_delay;
534  }
535 
536  gsl_spline_init (m_spline, t, v, nsamples_local);
537 
538 
539 
540  for (i=0;i<nsamples_local;i++) {
541  vd[i] = gsl_spline_eval (m_spline, td[i], m_interp_acc);
542  vf[i] = vd[i] - v[i]* m_cfd_fraction;
543  }
544 
545  //gsl_spline_free (m_spline); // Not sure if this is necessary
546 
547  gsl_spline_init (m_spline, t, vf, nsamples_local);
548 
549  cc.cls = this;
550  F.function = &ZdcRecChannelTool::fxCallback;
551  F.params = &cc;
552 
553 
554 
555 
556  //Get the range for root finding
557  vaux.assign(vf,vf+14);
558  tmin = std::min_element(vaux.begin(), vaux.end()) - vaux.begin();
559  tmax = std::max_element(vaux.begin(), vaux.end()) - vaux.begin();
560  if (tmin < tmax)
561  {
562  x_lo = t[tmin];
563  x_hi = t[tmax];
564  }
565  else
566  {
567  x_lo = t[tmax];
568  x_hi = t[tmin];
569  }
570  //std::cout << "::::::::::tmin,tmax " << x_lo << " " << x_hi << std::endl;
571 
572  /*
573  std::cout << "::::::::::V " << id.getString() << " "<< j << " " ;
574  for (i=0;i<nsamples_local;i++) {
575  std::cout << v[i] << " " ;
576  }
577  std::cout << std::endl;
578 
579  std::cout << "::::::::::VD " << id.getString() << " "<< j << " " ;
580  for (i=0;i<nsamples_local;i++) {
581  std::cout << vd[i] << " " ;
582  }
583  std::cout << std::endl;
584 
585  std::cout << "::::::::::VF " << id.getString() << " "<< j << " " ;
586  for (i=0;i<nsamples_local;i++) {
587  std::cout << vf[i] << " " ;
588  }
589  std::cout << std::endl;
590  */
591  //j++;
592 
593  gsl_root_fsolver_set (s, &F, x_lo, x_hi);
594 
595  iterations = 0;
596  do
597  {
598  iterations++;
599  status = gsl_root_fsolver_iterate(s);
600  r = gsl_root_fsolver_root(s);
601  x_lo = gsl_root_fsolver_x_lower(s);
602  x_hi = gsl_root_fsolver_x_upper(s);
603  status = gsl_root_test_interval (x_lo, x_hi, 0, 0.001);
604  //if (status == GSL_SUCCESS) printf ("Converged:\n");
605  }
606 
607  while (status == GSL_CONTINUE && iterations < max_iterations);
608  //std::cout << " COnvergence ----> " << iterations << " " << r << " " << x_lo << " " << x_hi << std::endl;
609  m_cfd_result.push_back(r-100); //subtraction accounts for zero padding
610  }
611  gsl_root_fsolver_free (s);
612  return 0;
613 }

◆ getTimingSinc()

int ZdcRecChannelTool::getTimingSinc ( const Identifier &  id,
const std::vector< std::vector< int > > &  wfm 
)
private

Definition at line 633 of file ZdcRecChannelTool.cxx.

634 {
635  unsigned int i = 0;
636  unsigned int k = 0;
637  int t = 0;
638  int tpeak = 0;
639  float tfpeak = 0.;
640  float vpeak = 0.;
641  float x = 0.;
642  float z = 0.;
643  float t_cor = 0. ;
644  int mSide = 0;
645  int mModule = 0;
646  int mType = 0;
647  int mChannel = 0;
648  int do_tcor = false;
649 
650 
651  std::vector<float> vt(7);
652  std::vector<int> y;
653 
654  std::vector<float>::iterator vf_it;
655  std::vector<std::vector<int> >::const_iterator vvi_it;
656 
657  //initialize the vectors
658  m_bwl_vpeak.resize(4,0);
659  m_bwl_tpeak.resize(4,0);
660 
661  //zId = id.get_identifier32().get_compact();
662  msg(MSG::DEBUG) << "--> ZDC : START OF MODIFICATION 2 " << endmsg ;
663  mSide = m_zdcId->side(id);
664  mModule = m_zdcId->module(id);
665  mType = m_zdcId->type(id);
666  mChannel = m_zdcId->channel(id);
667  msg(MSG::DEBUG) << "--> ZDC : ZdcRecChannelTool::getTimingSinc: "
668  << " id;Side;Module;Type;Channel: "
669  << id.getString() << ";"
670  << mSide << ";"
671  << mModule << ";"
672  << mType << ";"
673  << mChannel << endmsg;
674 
675  for (vvi_it = wfm.begin(); vvi_it<wfm.end(); ++vvi_it) {
676  //FIXME: Change to the method of ID identification
677  //if ((zId == 0xec000000) || (zId == 0xed000000) ||
678  // (zId == 0xec200000) || (zId == 0xed200000) ||
679  // (zId == 0xec400000) || (zId == 0xed400000) ||
680  // (zId == 0xec600000) || (zId == 0xed600000) )
681  // {
682  if (mType == 0) {
683 
684  y = *vvi_it;
685 
686  t = 0;
687  //Fill the interpolated vector
688  for (vf_it = m_wfm_bwl.begin(); vf_it != m_wfm_bwl.end(); ++vf_it ) {
689  z = 0.;
690  for (i=0;i<m_nsamples;i++) {
691  x = (TMath::Pi() *(t*0.1 - i*25.))/25. ;
692  if (fabs(x) < 1e-8) {
693  z = z + y[i]*1.0;
694  }
695  else {
696  z = z + y[i]*(sin(x)/(x));
697  }
698  }
699  *vf_it = z;
700  //if (id.get_identifier32().get_compact() == 0xec400000) {
701  //std::cout << "========== " << t << " " << *vf_it << std::endl;
702  //}
703  t++;
704  }
705 
706  //Get The maximum and the position of the peak
707  tpeak = std::max_element(m_wfm_bwl.begin(),m_wfm_bwl.end()) - m_wfm_bwl.begin();
708  vpeak = m_wfm_bwl[tpeak];
709 
710  //Account for 100 ps step
711  tfpeak = tpeak/10.;
712 
713 
714  //The correction for non-linear behavior
715  //Right now only two channels implemented:
716  //First hadronic modules for A and C sides
717  //zId = id.get_identifier32().get_compact();
718 
719  //Side C, Mod 0, Low Gain
720  //if (id.get_compact() == 0xec400000) {
721  if ( (mSide == -1) && (mType == 0) && (mModule == 0) && (mChannel == 0) ) {
722  t_cor = -777.277 +
723  33.685 * tfpeak +
724  -0.483265 * pow(tfpeak,2) +
725  0.00234842 * pow(tfpeak,3);
726  }
727 
728  //Side A, Mod 0, Low Gain
729  //if (id.get_compact() == 0xed400000) {
730  if ( (mSide == 1) && (mType == 0) && (mModule == 0) && (mChannel == 0) ) {
731  //piecewise fit
732  if (tfpeak <= 68) {
733  t_cor = -1749.27 +
734  80.5426 * tfpeak +
735  -1.229 * pow(tfpeak,2) +
736  0.006276 * pow(tfpeak,3);
737  }
738  else {
739  t_cor = -1965.52 +
740  84.0379 * tfpeak +
741  -1.19169 * pow(tfpeak,2) +
742  0.005659 * pow(tfpeak,3);
743  }
744  }
745 
746  if (do_tcor) m_bwl_tpeak[k] = t_cor;
747  else m_bwl_tpeak[k] = tpeak;
748  m_bwl_vpeak[k] = vpeak;
749  }
750  else {
751  m_bwl_tpeak[k] = -999;
752  m_bwl_vpeak[k] = -999;
753  }
754  k++;
755  }
756  return k;
757 }

◆ getTimingSinc2()

int ZdcRecChannelTool::getTimingSinc2 ( const Identifier &  id,
const std::vector< std::vector< int > > &  wfm 
)
private

Definition at line 764 of file ZdcRecChannelTool.cxx.

765 {
766  int i = 0;
767  int wfmIndex = 0;
768  //float energy = 0.;
769  //float w = 0.;
770  //float timing = 0.;
771  int error = 0;
772  int warning = 0;
773 
774  double Slices[10];
775  int NSlice = 5;
776  double gain=1.;
777  double pedestal=0;
778  double CFD_frac = 1.0;
779  bool corr=0;
780 
781  std::vector<std::vector<int> >::const_iterator vit;
782  std::vector<int> y;
783 
784  int mType = 0;
785 
786 
787 
788  //Identifier::value_type zId;
789 
790  /*
791  Id: Side Module Number Module Type
792  0xec000000 0 (C) 0 0
793  0xec400000 0 (C) 1 0
794  0xec800000 0 (C) 2 0
795  0xecc00000 0 (C) 3 0
796  0xed000000 1 (A) 0 0
797  0xed400000 1 (A) 1 0
798  0xed800000 1 (A) 2 0
799  0xedc00000 1 (A) 3 0
800  */
801  //Initialize the vectors
802  m_bwl_vpeak2.resize(4,0);
803  m_bwl_tpeak2.resize(4,0);
804 
805  mType = m_zdcId->type(id);
806 
807 
808 
809  //This controls if we get timing from all channels or
810  //from only a subset
811  wfmIndex = 2;
812  i = 0;
813  for (vit = wfm.begin(); vit<wfm.end(); ++vit) {
814  if ( (i < wfmIndex) && (mType == 0) )
815  {
816  y = *vit;
817 
818  // This is for the sinx/x interpolation
819 
820  for(int I=0;I<NSlice;I++) {Slices[I]=y[I];}
821  ZdcSignalSinc zdcSignalSinc(NSlice);
822 
823  zdcSignalSinc.process(Slices,gain,pedestal,CFD_frac,corr);
824 
825 
826  m_bwl_vpeak2[i] = zdcSignalSinc.getAmp() ;
827  m_bwl_tpeak2[i] = zdcSignalSinc.getTime() ;
828 
829  error = zdcSignalSinc.getError();
830  warning = zdcSignalSinc.getWarning();
831 
832  msg(MSG::DEBUG) << "--> ZDC : ZdcRecChannelTool::getTimingSinc2: "
833  << error << ";" << warning << endmsg;
834 
835  }
836  else {
837  m_bwl_vpeak2[i] = -999 ;
838  m_bwl_tpeak2[i] = -999 ;
839  }
840  i++;
841  }
842  return 0;
843 
844 }

◆ initialize()

StatusCode ZdcRecChannelTool::initialize ( )
virtual

Definition at line 85 of file ZdcRecChannelTool.cxx.

86 {
87  msg(MSG::INFO) << "Initializing " << name() << endmsg;
88 
89  //Get the pedestal information for the channels.
90  //For now, this is a file; later on it will be moved to a database
91 
92 
93  // Bandwidth interpolation method initializations
94  // making sure we use only 50 ps steps for the resolution
95 
96  if (m_bwl_time_resolution < 50) m_bwl_time_resolution = 50;
97  m_bwl_time_resolution = floor(m_bwl_time_resolution/50.) * 50;
98 
99  int s = m_nsamples*(s_SAMPLING_TIME/m_bwl_time_resolution);
100 
101  m_wfm_bwl.resize(s);
102 
103  msg(MSG::INFO) << "Using a time step of "
104  << m_bwl_time_resolution
105  << "ps for the Bandwidth Limited Sin(x)/x Interpolation: "
106  << " Vector Size = "
107  << s
108  << endmsg;
109 
110 
111 
112 
113  // GSL Interpolation functions initializations
114  m_interp_acc = gsl_interp_accel_alloc ();
115  // Thread-safe according to https://www.gnu.org/software/gsl/doc/html/roots.html
116  const gsl_interp_type *interp_type ATLAS_THREAD_SAFE = gsl_interp_akima;
117  m_spline = gsl_spline_alloc (interp_type, 14); //FIXME: TWICE THE SAMPLES
118 
119  //Load Mapping
120  msg(MSG::DEBUG) << "--> ZDC : START OF MODIFICATION 0" << endmsg ;
121 
122  const ZdcID* zdcId = nullptr;
123  if (detStore()->retrieve( zdcId ).isFailure() ) {
124  msg(MSG::ERROR) << "execute: Could not retrieve ZdcID object from the detector store" << endmsg;
125  return StatusCode::FAILURE;
126  }
127  else {
128  msg(MSG::DEBUG) << "execute: retrieved ZdcID" << endmsg;
129  }
130  m_zdcId = zdcId;
131 
132  msg(MSG::DEBUG) << "--> ZDC : END OF MODIFICATION 0" << endmsg ;
133  return StatusCode::SUCCESS;
134 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

◆ interfaceID()

const InterfaceID & ZdcRecChannelTool::interfaceID ( )
static

AlgTool InterfaceID.

Definition at line 44 of file ZdcRecChannelTool.cxx.

45 {
46  return IID_IZdcRecChannelTool;
47 
48 }

◆ makeRawFromDigits()

int ZdcRecChannelTool::makeRawFromDigits ( const ZdcDigitsCollection &  data_collection,
ZdcRawChannelCollection &  raw_collection 
)

Traditional code bellow (Energy by peak or sum, time by ratio of samples

Definition at line 151 of file ZdcRecChannelTool.cxx.

155 {
156  //ZdcRawChannel *z;
157  Identifier id;
158  float soma = 0.;
159  float pico = 0.;
160  float tzc = 0.;
161  float tsr = 0.;
162  float t0 = 0.;
163  float t1 = 0.;
164  float m = 0.;
165 
166  //bool pp = false;
167  int i = 0;
168  //int j = 0;
169  int k = 0;
170  //int s = 0;
171  int ped = 0;
172  int imax = 0;
173 
174  int mSide = 0;
175  int mModule = 0;
176  int mType = 0;
177  int mChannel = 0;
178 
179  int ncha = 0;
180 
181 
182  //TODO use BOOST containers here, more flexible and inteligent, and it will
183  // nevertless stay inside this method, so no fussing with complex collections
184 
185  std::vector<std::vector<int> > wfm;
186  std::vector<std::vector<int> >::iterator vit;
187  std::vector<int>::iterator it;
188 
189  std::vector<float> energy_sum;
190  std::vector<float> energy_peak;
191  std::vector<float> time_cfd;
192  std::vector<float> time_sratio;
193  std::vector<float> chi;
194 
195  std::vector<int> v;
196  std::vector<int> vd;
197  std::vector<int> vi;
198 
199  std::vector<float> vcfd1;
200  std::vector<float> vcfd2;
201  std::vector<float> vcfd3;
202 
203  //Identifier::value_type zId;
204 
205 
206  // Main procedure: Reduce the data to energy and time, and store it into
207  // ZdcRawDataCollection
208 
209  //digits_p = *mydata.begin();
210  //nsamples = digits_p->m_nSamples;
211  //vi.resize(2*nsamples);
212  vi.resize(14);
213 
214  for (const ZdcDigits* digits_p : mydata) {
215 
216  msg(MSG::DEBUG) << "--> ZDC : START OF MODIFICATION " << endmsg ;
217  id = digits_p->identify();
218  mSide = m_zdcId->side(id);
219  mModule = m_zdcId->module(id);
220  mType = m_zdcId->type(id);
221  mChannel = m_zdcId->channel(id);
222 
223  msg(MSG::DEBUG) << "--> ZDC : ZdcRecChannelTool::makeRawFromDigits: "
224  << " id;Side;Module;Type;Channel: "
225  << id.getString() << ";"
226  << mSide << ";"
227  << mModule << ";"
228  << mType << ";"
229  << mChannel << endmsg;
230 
231 
232  //1) Check the high gain. If it has saturation (max = FADC_SATURATION)
233  // drop and use the low gain
234  //
235  wfm.clear();
236  wfm.resize(2);
237  energy_sum.clear();
238  energy_peak.clear();
239  time_cfd.clear();
240  time_sratio.clear();
241  chi.clear();
242 
243  wfm[0] = digits_p->get_digits_gain0_delay0();
244  if (wfm[0].empty()) {
245  msg(MSG::DEBUG) << "ZERO SIZE g0d0 at ID " << id.getString() << endmsg;
246  wfm[0].resize(7);
247  }
248 
249  wfm[1] = digits_p->get_digits_gain1_delay0();
250  if (wfm[1].empty()) {
251  msg(MSG::DEBUG) << "ZERO SIZE g1d0 at ID " << id.getString() << endmsg;
252  wfm[1].resize(7);
253  }
254 
255 
256  msg(MSG::DEBUG) << " ------- 1 " << endmsg;
257 
258  //2) Check to see if there are delayed information
259  // module type 0 -> (full Energy) delayed info
260  // module type 1 -> (segmented) do not have delayed info
261  // - use the identifiers here -
262  //
263  //if (((id.get_identifier32().get_compact() >> 21) & 1) == 0) {
264 
265  //zId = id.get_identifier32().get_compact();
266 
267  //Only type 0 (Total Module Energy) has delay information
268  //if ( ( (zId == 0xec000000) || (zId == 0xed000000) ||
269  // (zId == 0xec200000) || (zId == 0xed200000) ||
270  // (zId == 0xec400000) || (zId == 0xed400000) ||
271  // (zId == 0xec600000) || (zId == 0xed600000) ) ) {
272  if (mType == 0) {
273  //std::cout << "*** Resize 4 at Id " << zId << std::endl ;
274  wfm.resize(4);
275  wfm[2] = digits_p->get_digits_gain0_delay1();
276  wfm[3] = digits_p->get_digits_gain1_delay1();
277  }
278 
279 
280 
281 
282  msg(MSG::DEBUG) << " ------- 2 " << endmsg;
283  //3) Subtratct baseline pedestal
284  //We need to be carefull. Sometimes the vector size here is zero (PPM flaw) and
285  //the code crashs if we do not treat this.
286  i = 0;
287  for (vit = wfm.begin(); vit<wfm.end(); ++vit) {
288  if (vit->empty()) vit->resize(7);
289  ped = *vit->begin();
290  for (it=vit->begin(); it<vit->end();++it) {
291  (*it) -= ped;
292  //if ((i==1) || (i==3) ) (*it) = (*it) * s_GAIN_RATIO;
293  }
294  i++;
295  }
296  /*
297  if (((id.get_identifier32().get_compact() >> 21) & 1) == 0) {
298  if ( ( *(std::max_element(wfm[0].begin(),wfm[0].end()) ) > 300 ) ||
299  ( *(std::max_element(wfm[1].begin(),wfm[1].end()) ) > 300 ) ||
300  ( *(std::max_element(wfm[2].begin(),wfm[2].end()) ) > 300 ) ||
301  ( *(std::max_element(wfm[3].begin(),wfm[3].end()) ) > 300 ) ) {
302  //std::cout << "&&&&& OVFLW &&&&&" << std::endl;
303  msg(MSG::DEBUG) << "%%%%%%%% ID " << id.getString() << endmsg;
304  msg(MSG::DEBUG) << "%%%%%%%% g0d0 " << wfm[0] << endmsg;
305  msg(MSG::DEBUG) << "%%%%%%%% g1d0 " << wfm[1] << endmsg;
306  msg(MSG::DEBUG) << "%%%%%%%% g0d1 " << wfm[2] << endmsg;
307  msg(MSG::DEBUG) << "%%%%%%%% g1d1 " << wfm[3] << endmsg;
308 
309  }
310  }
311  else {
312  if ( ( *(std::max_element(wfm[0].begin(),wfm[0].end()) ) > 20 ) ||
313  (*(std::max_element(wfm[1].begin(),wfm[1].end()) ) > 20 ) ) {
314  //std::cout << "**** g1d0 OVFLW ****" << std::endl;
315  msg(MSG::DEBUG) << "******** ID " << id.getString() << endmsg;
316  msg(MSG::DEBUG) << "******** g0d0 " << wfm[0] << endmsg;
317  msg(MSG::DEBUG) << "******** g1d0 " << wfm[1] << endmsg;
318 
319  }
320  }
321  */
322  msg(MSG::DEBUG) << " ------- 3 " << endmsg;
323  //4) Calculate the quantities of interest. For now it will be Energy (by SUM and Peak)
324  // and timing (using undelayed, delayed and combined)
325  // Include a quality factor chi which will follow design definitions for
326  // chi < 0 -> something very bad (for example saturation)
327  // 0 >= chi >= 1 -> depends on the algorithm to be used; 1 is very good
328 
329  k = 0;
330 
332  // GSL for CFD Timing
334 
335  getTimingCFD(id, wfm);
336 
338  // Sinx/x Timing (GSL)
340  // Not yet needed
341  //getTimingSinc(id, wfm);
342 
343 
345  // Sinx/x Timing (Andrei's coding)
347 
348  getTimingSinc2(id, wfm);
349 
354  for (vit = wfm.begin(); vit<wfm.end(); ++vit) {
355  v = *vit;
356  soma = std::accumulate(v.begin(), v.end(), 0) ;
357  pico = *(std::max_element(v.begin(),v.end()) );
358  energy_sum.push_back ((float) soma);
359  energy_peak.push_back((float) pico);
360  if (pico >= (s_FADC_SATURATION -42)) {
361  chi.push_back(-1);
362  }
363  else {
364  chi.push_back(0); //nothing to say here
365  }
366 
367 
368  // Now, lets also calculate the timing by the ratio of the
369  // second to the third sample. Because the correction is module
370  // dependent, we postpone the calibration to the next stage
371  // Time will be
372  // t = (r-R0)/(R1-R0) ns
373  // where r = v[imax -1]/v[imax], (imax = peak)
374  // R0, R1 are constants that depends on the module
375 
376  // Get the position of the maximum
377  // Take care of div by 0
378  imax = std::max_element(v.begin(),v.end()) - v.begin();
379  if ((v[imax] != 0) && (imax>0)) {
380  tsr = ((float) v[imax-1])/v[imax];
381  }
382  else {
383  tsr = v[imax];
384  }
385 
386  time_sratio.push_back(tsr);
387 
388 
389  msg(MSG::DEBUG) << "--> ZDC : ZdcRecChannelTool " << id.getString() <<
390  " Type=" << k <<
391  " Energy Sum=" << soma <<
392  " Energy Peak=" << pico <<
393  " Chi=" << *(chi.end() -1) <<
394  " t0=" << t0 << " t1=" << t1 <<
395  " m=" << m <<
396  " Time by CFD=" << tzc <<
397  " Peak at imax = " << imax <<
398  " v[imax-1] = " << (imax > 0 ? v[imax-1] : 0) <<
399  " V[imax] = " << v[imax] <<
400  " Sample Ratio A1/A2=" << tsr << endmsg;
401  k++;
402  }
403 
404  msg(MSG::DEBUG) << " ------- 5 " << endmsg;
405  ZdcRawChannel *z = new ZdcRawChannel(id);
406 
407  msg(MSG::DEBUG) << "CFD EXACT***** ID "<< id.getString() << m_cfd_result << endmsg;
408  msg(MSG::DEBUG) << "CFD APPRO***** ID "<< id.getString() << m_bwl_tpeak << endmsg;
409 
410  z->setSize(3*k); //FIXME very important, but not smart move ...
411  /*
412  * What can go here:
413  *
414  * energy_peak -> max sample peak
415  * energy_sum -> sum of samples
416  * time_sratio -> from ratio of samples
417  *
418  * m_bwl_vpeak -> Energy peak from sinc interpolation
419  * m_bwl_tpeak -> Time at peak from sinc interpolation
420  *
421  * m_bwl_vpeak2 -> Energy peak from sinc interpolation (AP code)
422  * m_bwl_tpeak2 -> Time at peak from sinc interpolation (AP code)
423  *
424  * m_cfd_result -> timing by software CFD
425  */
426  for (i=0;i<k;i++) {
427  z->setEnergy (i, energy_peak[i]);
428  z->setTime (i, time_sratio[i]);
429  z->setChi (i, chi[i]);
430 
431  z->setEnergy (i+k, m_bwl_vpeak2[i]);
432  z->setTime (i+k, m_bwl_tpeak2[i]);
433  z->setChi (i+k, chi[i]);
434 
435  z->setEnergy (i+2*k, energy_sum[i]);
436  z->setTime (i+2*k, m_cfd_result[i]);
437  z->setChi (i+2*k, chi[i]);
438 
439  //z->setEnergy(i+k, energy_peak[i]);
440  //z->setTime(i,time_cfd[i]);
441 
442 
443  }
444  ncha++;
445  ChannelCollection.push_back(z);
446  msg(MSG::DEBUG) << " ------- 6 " << endmsg;
447  }
448  msg(MSG::DEBUG) << "--> ZDC : ZdcRecChannelTool ChannelCollection size " << ChannelCollection.size() << endmsg ;
449  return ncha;
450 }

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24  {
25  return this->msgStream();
26  }

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

◆ 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< AlgTool > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray &  handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::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< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase &  )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308  {
309  // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310  // << " size: " << m_vhka.size() << endmsg;
311  for (auto &a : m_vhka) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ m_bwl_time_resolution

int ZdcRecChannelTool::m_bwl_time_resolution
private

Definition at line 115 of file ZdcRecChannelTool.h.

◆ m_bwl_tpeak

std::vector<float> ZdcRecChannelTool::m_bwl_tpeak
private

Definition at line 118 of file ZdcRecChannelTool.h.

◆ m_bwl_tpeak2

std::vector<float> ZdcRecChannelTool::m_bwl_tpeak2
private

Definition at line 123 of file ZdcRecChannelTool.h.

◆ m_bwl_vpeak

std::vector<float> ZdcRecChannelTool::m_bwl_vpeak
private

Definition at line 117 of file ZdcRecChannelTool.h.

◆ m_bwl_vpeak2

std::vector<float> ZdcRecChannelTool::m_bwl_vpeak2
private

Definition at line 122 of file ZdcRecChannelTool.h.

◆ m_cfd_delay

float ZdcRecChannelTool::m_cfd_delay
private

Definition at line 142 of file ZdcRecChannelTool.h.

◆ m_cfd_fraction

float ZdcRecChannelTool::m_cfd_fraction
private

Definition at line 141 of file ZdcRecChannelTool.h.

◆ m_cfd_result

std::vector<float> ZdcRecChannelTool::m_cfd_result
private

Definition at line 111 of file ZdcRecChannelTool.h.

◆ m_delta

int ZdcRecChannelTool::m_delta
private

Definition at line 133 of file ZdcRecChannelTool.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_interp_acc

gsl_interp_accel* ZdcRecChannelTool::m_interp_acc {}
private

Definition at line 100 of file ZdcRecChannelTool.h.

◆ m_nsamples

unsigned int ZdcRecChannelTool::m_nsamples
private

Definition at line 130 of file ZdcRecChannelTool.h.

◆ m_sample_time

float ZdcRecChannelTool::m_sample_time
private

Definition at line 131 of file ZdcRecChannelTool.h.

◆ m_spline

gsl_spline* ZdcRecChannelTool::m_spline {}
private

Definition at line 101 of file ZdcRecChannelTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_wfm_bwl

std::vector<float> ZdcRecChannelTool::m_wfm_bwl
private

Definition at line 116 of file ZdcRecChannelTool.h.

◆ m_zdcId

const ZdcID* ZdcRecChannelTool::m_zdcId {}
private

Definition at line 156 of file ZdcRecChannelTool.h.

◆ m_zeroSupress

int ZdcRecChannelTool::m_zeroSupress
private

Definition at line 134 of file ZdcRecChannelTool.h.

◆ s_CFD_FRACTION

const int ZdcRecChannelTool::s_CFD_FRACTION = -3
staticprivate

Definition at line 144 of file ZdcRecChannelTool.h.

◆ s_FADC_SATURATION

const int ZdcRecChannelTool::s_FADC_SATURATION = 1022
staticprivate

Definition at line 137 of file ZdcRecChannelTool.h.

◆ s_GAIN_RATIO

const int ZdcRecChannelTool::s_GAIN_RATIO = 1
staticprivate

Definition at line 152 of file ZdcRecChannelTool.h.

◆ s_HALF_SAMPLING_TIME

const int ZdcRecChannelTool::s_HALF_SAMPLING_TIME = 12500
staticprivate

Definition at line 146 of file ZdcRecChannelTool.h.

◆ s_SAMPLING_TIME

const int ZdcRecChannelTool::s_SAMPLING_TIME = 25000
staticprivate

Definition at line 145 of file ZdcRecChannelTool.h.

The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
xAOD::iterator
JetConstituentVector::iterator iterator
Definition: JetConstituentVector.cxx:68
beamspotman.r
def r
Definition: beamspotman.py:676
python.SystemOfUnits.s
int s
Definition: SystemOfUnits.py:131
get_generator_info.result
result
Definition: get_generator_info.py:21
ZdcRecChannelTool::getTimingSinc2
int getTimingSinc2(const Identifier &id, const std::vector< std::vector< int > > &wfm)
Definition: ZdcRecChannelTool.cxx:764
python.SystemOfUnits.m
int m
Definition: SystemOfUnits.py:91
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
ZdcSignalSinc
Definition: ZdcSignalSinc.h:14
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
hist_file_dump.d
d
Definition: hist_file_dump.py:137
ZdcRecChannelTool::s_FADC_SATURATION
static const int s_FADC_SATURATION
Definition: ZdcRecChannelTool.h:137
ZdcID::module
int module(const Identifier &id) const
Definition: ZdcID.h:163
ZdcRecChannelTool::m_zdcId
const ZdcID * m_zdcId
Definition: ZdcRecChannelTool.h:156
accumulate
bool accumulate(AccumulateMap &map, std::vector< module_t > const &modules, FPGATrackSimMatrixAccumulator const &acc)
Accumulates an accumulator (e.g.
Definition: FPGATrackSimMatrixAccumulator.cxx:22
ZdcRecChannelTool::m_cfd_delay
float m_cfd_delay
Definition: ZdcRecChannelTool.h:142
conifer::pow
constexpr int pow(int x)
Definition: conifer.h:20
CaloCondBlobAlgs_fillNoiseFromASCII.gain
gain
Definition: CaloCondBlobAlgs_fillNoiseFromASCII.py:110
ALFA_EventTPCnv_Dict::t0
std::vector< ALFA_RawData_p1 > t0
Definition: ALFA_EventTPCnvDict.h:42
ZdcRecChannelTool::m_cfd_result
std::vector< float > m_cfd_result
Definition: ZdcRecChannelTool.h:111
ALFA_EventTPCnv_Dict::t1
std::vector< ALFA_RawDataCollection_p1 > t1
Definition: ALFA_EventTPCnvDict.h:43
skel.it
it
Definition: skel.GENtoEVGEN.py:423
ZdcDigits
Definition: ZdcDigits.h:28
ZdcRecChannelTool::s_SAMPLING_TIME
static const int s_SAMPLING_TIME
Definition: ZdcRecChannelTool.h:145
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
ZdcRecChannelTool::m_bwl_tpeak2
std::vector< float > m_bwl_tpeak2
Definition: ZdcRecChannelTool.h:123
ZdcRecChannelTool::m_cfd_fraction
float m_cfd_fraction
Definition: ZdcRecChannelTool.h:141
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ZdcRecChannelTool::m_bwl_tpeak
std::vector< float > m_bwl_tpeak
Definition: ZdcRecChannelTool.h:118
x
#define x
empty
bool empty(TH1 *h)
Definition: computils.cxx:294
mergePhysValFiles.end
end
Definition: DataQuality/DataQualityUtils/scripts/mergePhysValFiles.py:93
AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:95
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
CCallbackHolder
Definition: ZdcRecChannelTool.h:67
ZdcRecChannelTool::m_wfm_bwl
std::vector< float > m_wfm_bwl
Definition: ZdcRecChannelTool.h:116
ZdcRawChannel
Definition: ZdcRawChannel.h:24
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
ZdcRecChannelTool::m_nsamples
unsigned int m_nsamples
Definition: ZdcRecChannelTool.h:130
ZdcRecChannelTool
Definition: ZdcRecChannelTool.h:72
lumiFormat.i
int i
Definition: lumiFormat.py:92
z
#define z
Identifier
Definition: DetectorDescription/Identifier/Identifier/Identifier.h:32
endmsg
#define endmsg
Definition: AnalysisConfig_Ntuple.cxx:63
ZdcRecChannelTool::m_interp_acc
gsl_interp_accel * m_interp_acc
Definition: ZdcRecChannelTool.h:100
ZdcRecChannelTool::getTimingCFD
int getTimingCFD(const Identifier &id, const std::vector< std::vector< int > > &wfm)
Definition: ZdcRecChannelTool.cxx:473
test_pyathena.parent
parent
Definition: test_pyathena.py:15
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
ZdcRecChannelTool::m_sample_time
float m_sample_time
Definition: ZdcRecChannelTool.h:131
imax
int imax(int i, int j)
Definition: TileLaserTimingTool.cxx:33
AthAlgTool::AthAlgTool
AthAlgTool()
Default constructor:
ZdcRecChannelTool::fxCallback
static double fxCallback(double d, void *v)
Definition: ZdcRecChannelTool.h:105
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
id
SG::auxid_t id
Definition: Control/AthContainers/Root/debug.cxx:191
ZdcRecChannelTool::m_bwl_vpeak2
std::vector< float > m_bwl_vpeak2
Definition: ZdcRecChannelTool.h:122
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:192
validateBDTTau.vt
vt
Definition: validateBDTTau.py:43
ZdcRecChannelTool::m_delta
int m_delta
Definition: ZdcRecChannelTool.h:133
Example_ReadSampleNoise.ped
ped
Definition: Example_ReadSampleNoise.py:45
python.PyAthena.v
v
Definition: PyAthena.py:157
DiTauMassTools::MaxHistStrategyV2::e
e
Definition: PhysicsAnalysis/TauID/DiTauMassTools/DiTauMassTools/HelperFunctions.h:26
a
TList * a
Definition: liststreamerinfos.cxx:10
ZdcRecChannelTool::m_bwl_vpeak
std::vector< float > m_bwl_vpeak
Definition: ZdcRecChannelTool.h:117
y
#define y
h
ZdcRecChannelTool::m_bwl_time_resolution
int m_bwl_time_resolution
Definition: ZdcRecChannelTool.h:115
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
DEBUG
#define DEBUG
Definition: page_access.h:11
F
#define F(x, y, z)
Definition: MD5.cxx:112
ZdcID::side
int side(const Identifier &id) const
Values of different levels (failure returns 0)
Definition: ZdcID.h:157
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
ZdcID::type
int type(const Identifier &id) const
Definition: ZdcID.h:169
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:616
declareProperty
#define declareProperty(n, p, h)
Definition: BaseFakeBkgTool.cxx:15
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:790
merge.status
status
Definition: merge.py:17
ATLAS_THREAD_SAFE
#define ATLAS_THREAD_SAFE
Definition: checker_macros.h:211
I
#define I(x, y, z)
Definition: MD5.cxx:116
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
PowhegControl_ttFCNC_NLO.params
params
Definition: PowhegControl_ttFCNC_NLO.py:226
error
Definition: IImpactPoint3dEstimator.h:70
ZdcID
Definition: ZdcID.h:25
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
readCCLHist.float
float
Definition: readCCLHist.py:83
ZdcRecChannelTool::m_spline
gsl_spline * m_spline
Definition: ZdcRecChannelTool.h:101
TSU::T
unsigned long long T
Definition: L1TopoDataTypes.h:35
ZdcRecChannelTool::m_zeroSupress
int m_zeroSupress
Definition: ZdcRecChannelTool.h:134
fitman.k
k
Definition: fitman.py:528
python.handimod.cc
int cc
Definition: handimod.py:523
ZdcID::channel
int channel(const Identifier &id) const
Definition: ZdcID.h:175
Generated on Fri Jun 28 2024 22:11:49 for ATLAS Offline Software by   doxygen 1.8.18