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

#include <TBBeamQualityMC.h>

Inheritance diagram for TBBeamQualityMC:
Collaboration diagram for TBBeamQualityMC:

Public Member Functions

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

Protected Member Functions

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

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

StatusCode getXcryoYtable (float &x, float &y, float &eBeam)
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

bool m_readFileforXcryo
 Get Xcryo and Ytable from a text file. More...
 
int m_nRun
 
float m_bm_cut_x
 
float m_bm_cut_y
 
bool m_check_trackpar
 
bool m_check_primary
 
std::vector< int > m_scint_prim
 
bool m_check_veto
 
bool m_check_clus
 
bool m_check_trackreco
 
std::string m_clusterCollName
 
DataObjIDColl m_extendedExtraObjects
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 11 of file TBBeamQualityMC.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ TBBeamQualityMC()

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

Definition at line 17 of file TBBeamQualityMC.cxx.

17  :
18  AthAlgorithm(name, pSvcLocator),
19  m_readFileforXcryo(true),
20  m_nRun(0)
21 {
22  declareProperty("ReadFileforXcryo", m_readFileforXcryo);
23  declareProperty("CheckTrackParams",m_check_trackpar=false);
24  declareProperty("TrackCutX",m_bm_cut_x=35.);
25  declareProperty("TrackCutY",m_bm_cut_y=35.);
26  declareProperty("CheckPrimaryTrack",m_check_primary=true);
27  declareProperty("ScintForPrimaryTrack",m_scint_prim);
28  declareProperty("CheckVetoScint",m_check_veto=true);
29  declareProperty("CheckClusters",m_check_clus=false);
30  declareProperty("CheckTrackReco",m_check_trackreco=false);
31  declareProperty("ClusterCollectionName",m_clusterCollName="CaloTopoClusters");
32 }

◆ ~TBBeamQualityMC()

virtual TBBeamQualityMC::~TBBeamQualityMC ( )
inlinevirtual

Definition at line 16 of file TBBeamQualityMC.h.

16 {};

Member Function Documentation

◆ declareGaudiProperty() [1/4]

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

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

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

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [3/4]

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

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

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

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

Definition at line 199 of file AthCommonDataStore.h.

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

◆ declareProperty() [1/6]

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

Declare a new Gaudi property.

Parameters
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< Algorithm > >::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< Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType  
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::declareProperty ( Gaudi::Property< T > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
148  }

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

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

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

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

StatusCode TBBeamQualityMC::execute ( )
overridevirtual

Print an informatory message:

Definition at line 39 of file TBBeamQualityMC.cxx.

39  {
41  ATH_MSG_DEBUG ( "in execute()" );
42 
43  // Retrieve Event Info
44  const TBEventInfo* theEventInfo = nullptr;
45  ATH_CHECK( evtStore()->retrieve(theEventInfo,"TBEventInfo") );
46 
47  // Fill run header
48  m_nRun = theEventInfo->getRunNum();
49  float beamMom = theEventInfo->getBeamMomentum();
50  float xCryo = -9999;
51  float yTable = -9999;
52  if(m_readFileforXcryo) {
53  if (!this->getXcryoYtable(xCryo, yTable, beamMom)) {
54  ATH_MSG_ERROR ( "xCryo and yTable are not found" );
55  return StatusCode::FAILURE;
56  }
57  } else {
58  xCryo = theEventInfo->getCryoX();
59  yTable = theEventInfo->getTableY();
60  }
61  ATH_MSG_VERBOSE ( "xCryo = " << xCryo );
62  ATH_MSG_VERBOSE ( "yTable = " << yTable );
63  ATH_MSG_VERBOSE ( "nRun = " << m_nRun );
64  ATH_MSG_VERBOSE ( "beamMom = " << beamMom );
65 
66 
67  if(m_check_trackpar) { // check for track parameters
68 
69  TBTrack *track = nullptr;
70  ATH_CHECK( evtStore()->retrieve(track,"Track") );
71 
72  const float zCalo = 30000.;
73  float beam_coor_x = track->getUslope()*zCalo + track->getUintercept() + xCryo;
74  float beam_coor_y = track->getVslope()*zCalo + track->getVintercept();
75 
76  bool sx1 = (beam_coor_x > xCryo - m_bm_cut_x) && (beam_coor_x < xCryo + m_bm_cut_x);
77  bool sx2 = (std::log(track->getChi2_u())>-10 && track->getChi2_u()!=1000);
78  bool sx3 = (beam_coor_x!=0 || track->getUintercept()!=0 || track->getUslope()!=0 );
79 
80  bool sy1 = (beam_coor_y > yTable - m_bm_cut_y) && (beam_coor_y < yTable + m_bm_cut_y);
81  bool sy2 = (std::log(track->getChi2_v())>-10 && track->getChi2_v()!=1000);
82  bool sy3 = (beam_coor_y!=0 || track->getVintercept()!=0 || track->getVslope()!=0 );
83 
84  if(! (sx1 && sx2 && sx3 && sy1 && sy2 && sy3)) {
85  ATH_MSG_DEBUG ( "CheckTrackParams failed: "<<sx1<<" "<<sx2<<" "<<sx3<<" / "<<sy1<<" "<<sy2<<" "<<sy3);
86  //return StatusCode::FAILURE;
87  //return StatusCode::RECOVERABLE;
88  setFilterPassed(false);
89  return StatusCode::SUCCESS;
90  }
91  }
92 
93  if(m_check_primary) { // check for primary energy deposition in vector of scintilators
94 
95  std::vector<bool> has_energy (m_scint_prim.size());
96  for(int ll=0; ll<(int)m_scint_prim.size(); ++ll) has_energy[ll]=false;
97 
98  const LArG4H6FrontHitCollection *frontcoll = nullptr;
99  ATH_CHECK( evtStore()->retrieve(frontcoll,"Front::Hits") );
100 
101  int scnum;
102  unsigned i;
103  for (const LArG4H6FrontHit* hit : *frontcoll) {
104  scnum = hit->GetSC();
105  if(scnum <= 0) continue; // not a scintilator hit
106  for(i=0; i<m_scint_prim.size(); ++i) {
107  ATH_MSG_DEBUG ( "Scint: "<<scnum<<", track: "<<hit->GetTrackID()<<", energy: "<<hit->GetEdep());
108  if( (scnum==m_scint_prim[i]) && (hit->GetTrackID() == 1) && (hit->GetEdep() > 0.)) {
109  has_energy[i] = true;
110  }
111  }
112  }
113 
114  const LArG4H6FrontHitCollection *movecoll = nullptr;
115  ATH_CHECK( evtStore()->retrieve(movecoll,"Movable::Hits") );
116  for (const LArG4H6FrontHit* hit : *movecoll) {
117  scnum = hit->GetSC();
118  if(scnum <= 0) continue; // not a scintilator hit
119  for(i=0; i<m_scint_prim.size(); ++i) {
120  ATH_MSG_DEBUG ( "Scint: "<<scnum<<", track: "<<hit->GetTrackID()<<", energy: "<<hit->GetEdep());
121  if( (scnum==m_scint_prim[i]) && (hit->GetTrackID() == 1) && (hit->GetEdep() > 0.)) {
122  has_energy[i] = true;
123  }
124  }
125  }
126  for(int ll=0; ll<(int)m_scint_prim.size(); ++ll) has_energy[0] = has_energy[0] && has_energy[ll];
127  if(!has_energy[0]) {
128  ATH_MSG_DEBUG ( "CheckPrimaryTrack failed, no energy deposit in all asked scint." );
129  //return StatusCode::FAILURE;
130  //return StatusCode::RECOVERABLE;
131  setFilterPassed(false);
132  return StatusCode::SUCCESS;
133  }
134  }
135 
136  if(m_check_veto) { // check if there is a hit in veto scint.
137  int scnum;
138  const LArG4H6FrontHitCollection *movecoll = nullptr;
139  ATH_CHECK( evtStore()->retrieve(movecoll,"Movable::Hits") );
140  for (const LArG4H6FrontHit* hit : *movecoll) {
141  scnum = hit->GetSC();
142  if(scnum <= 0) continue;
143  if(hit->GetSC() == 5 && (hit->GetTrackID() == 1)) {
144  ATH_MSG_DEBUG ( "CheckVeto failed ");
145  //return StatusCode::FAILURE;
146  //return StatusCode::RECOVERABLE;
147  setFilterPassed(false);
148  return StatusCode::SUCCESS;
149  }
150  }
151  }
152 
153  if(m_check_clus) { // Rejected if no good cluster
154  const CaloClusterContainer* cc = nullptr;
156  bool haveit=false;
157  for (const CaloCluster * theCluster : *cc) {
158  double cL = theCluster->getMomentValue(CaloClusterMoment::CENTER_LAMBDA);
159  double eD = theCluster->getMomentValue(CaloClusterMoment::FIRST_ENG_DENS);
160  unsigned int cSize = theCluster->getClusterSize();
161  if(cL != 0 && eD > 0. && cSize > 1) {
162  haveit = true;
163  break;
164  }
165  }
166  if (!haveit) {
167  ATH_MSG_DEBUG ( "CheckClusters failed ");
168  setFilterPassed(false);
169  return StatusCode::SUCCESS;
170  }
171  }
172 
173 
174  if(m_check_trackreco) { // Rejected if bad track reco
175  const TBTrack* mytrack = nullptr;
176  const TBEventInfo* mytbeinfo = nullptr;
177  ATH_CHECK( evtStore()->retrieve(mytrack,"Track" ) );
178  ATH_CHECK( evtStore()->retrieve(mytbeinfo,"TBEventInfo" ) );
179  double chi2 = mytrack->getChi2_u();
180  if(chi2 <= 0 || std::log(chi2)<=-10 || chi2 == 1000 || mytrack->getUslope() == 0 || mytrack->getUintercept() == 0 || mytbeinfo->getCryoX() == mytrack->getCryoHitu()) {
181  ATH_MSG_DEBUG ( "CheckTrackReco in X failed ");
182  ATH_MSG_DEBUG ( chi2 << " / " << std::log(chi2) << " :: " << mytrack->getUslope() << " / " << mytrack->getUintercept() << " / " << mytbeinfo->getCryoX() << " / " << mytrack->getCryoHitu() );
183  setFilterPassed(false);
184  return StatusCode::SUCCESS;
185  }
186  chi2 = mytrack->getChi2_v();
187  if(chi2 <= 0 || std::log(chi2)<=-10 || chi2 == 1000 || mytrack->getVslope() == 0 || mytrack->getVintercept() == 0 || mytbeinfo->getTableY() == mytrack->getCryoHitv()) {
188  ATH_MSG_DEBUG ( "CheckTrackReco in Y failed ");
189  ATH_MSG_DEBUG ( chi2 << " / " << std::log(chi2) << " :: " << mytrack->getVslope() << " / " << mytrack->getVintercept() << " / " << mytbeinfo->getTableY() << " / " << mytrack->getCryoHitv() );
190  setFilterPassed(false);
191  return StatusCode::SUCCESS;
192  }
193 
194  }
195  setFilterPassed(true);
196  return StatusCode::SUCCESS;
197 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

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

◆ finalize()

StatusCode TBBeamQualityMC::finalize ( )
overridevirtual

Definition at line 199 of file TBBeamQualityMC.cxx.

199  {
200  ATH_MSG_DEBUG ( "in finalize()" );
201  return StatusCode::SUCCESS;
202 }

◆ getXcryoYtable()

StatusCode TBBeamQualityMC::getXcryoYtable ( float &  x,
float &  y,
float &  eBeam 
)
private

Definition at line 205 of file TBBeamQualityMC.cxx.

205  {
206 
207  // not good to put here (just a workaround) - joh
208  std::string m_txtFileWithXY = "xcryo_ytable.txt";
209 
210  ATH_MSG_DEBUG ( "in getXcryoYtable(float x, float y)" );
211  std::ifstream xyFile;
212  std::string line;
213  std::string filename = PathResolver::find_file(m_txtFileWithXY, "DATAPATH");
214  xyFile.open(filename.c_str());
215  if (!xyFile.is_open()) {
216  ATH_MSG_ERROR ( "File " << m_txtFileWithXY << " fail to open in $DATAPATH");
217  return StatusCode::FAILURE;
218  }
219  while ( getline(xyFile, line, '\n') ) {
220  int run;
221  std::istringstream buf(line);
222  e = 0;
223  buf >> run >> x >> y >> e;
224  ATH_MSG_DEBUG ( "m_nRun,run,x,y,e= "<<m_nRun<<" "<<run<<" "<<x<<" "<<y<<" "<<e);
225  if (run==m_nRun && xyFile.good()) return StatusCode::SUCCESS;
226  }
227  return StatusCode::FAILURE;
228 }

◆ initialize()

StatusCode TBBeamQualityMC::initialize ( )
overridevirtual

Definition at line 34 of file TBBeamQualityMC.cxx.

34  {
35  ATH_MSG_DEBUG ( "in initialize()" );
36  return StatusCode::SUCCESS;
37 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::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< Algorithm >::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< Algorithm > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< Algorithm > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

66  {
68 
69  if (sc.isFailure()) {
70  return sc;
71  }
72  ServiceHandle<ICondSvc> cs("CondSvc",name());
73  for (auto h : outputHandles()) {
74  if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
75  // do this inside the loop so we don't create the CondSvc until needed
76  if ( cs.retrieve().isFailure() ) {
77  ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
78  return StatusCode::SUCCESS;
79  }
80  if (cs->regHandle(this,*h).isFailure()) {
81  sc = StatusCode::FAILURE;
82  ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
83  << " with CondSvc");
84  }
85  }
86  }
87  return sc;
88 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

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_bm_cut_x

float TBBeamQualityMC::m_bm_cut_x
private

Definition at line 28 of file TBBeamQualityMC.h.

◆ m_bm_cut_y

float TBBeamQualityMC::m_bm_cut_y
private

Definition at line 29 of file TBBeamQualityMC.h.

◆ m_check_clus

bool TBBeamQualityMC::m_check_clus
private

Definition at line 37 of file TBBeamQualityMC.h.

◆ m_check_primary

bool TBBeamQualityMC::m_check_primary
private

Definition at line 32 of file TBBeamQualityMC.h.

◆ m_check_trackpar

bool TBBeamQualityMC::m_check_trackpar
private

Definition at line 30 of file TBBeamQualityMC.h.

◆ m_check_trackreco

bool TBBeamQualityMC::m_check_trackreco
private

Definition at line 38 of file TBBeamQualityMC.h.

◆ m_check_veto

bool TBBeamQualityMC::m_check_veto
private

Definition at line 35 of file TBBeamQualityMC.h.

◆ m_clusterCollName

std::string TBBeamQualityMC::m_clusterCollName
private

Definition at line 40 of file TBBeamQualityMC.h.

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_nRun

int TBBeamQualityMC::m_nRun
private

Definition at line 26 of file TBBeamQualityMC.h.

◆ m_readFileforXcryo

bool TBBeamQualityMC::m_readFileforXcryo
private

Get Xcryo and Ytable from a text file.

Definition at line 24 of file TBBeamQualityMC.h.

◆ m_scint_prim

std::vector<int> TBBeamQualityMC::m_scint_prim
private

Definition at line 33 of file TBBeamQualityMC.h.

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
TBBeamQualityMC::m_scint_prim
std::vector< int > m_scint_prim
Definition: TBBeamQualityMC.h:33
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
TBTrack::getUslope
double getUslope() const
Definition: TBTrack.h:52
TBBeamQualityMC::m_check_clus
bool m_check_clus
Definition: TBBeamQualityMC.h:37
checkFileSG.line
line
Definition: checkFileSG.py:75
TBBeamQualityMC::m_check_trackpar
bool m_check_trackpar
Definition: TBBeamQualityMC.h:30
TBTrack
Definition: TBTrack.h:20
detail::ll
long long ll
Definition: PrimitiveHelpers.h:47
PathResolver::find_file
static std::string find_file(const std::string &logical_file_name, const std::string &search_path, SearchType search_type=LocalSearch)
Definition: PathResolver.cxx:251
CaloCellPos2Ntuple.int
int
Definition: CaloCellPos2Ntuple.py:24
AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
TBTrack::getUintercept
double getUintercept() const
Definition: TBTrack.h:54
CaloClusterContainer
Storable container for CaloCluster.
Definition: Calorimeter/CaloEvent/CaloEvent/CaloClusterContainer.h:37
TBEventInfo::getRunNum
unsigned int getRunNum() const
Definition: TBEventInfo.h:67
TBBeamQualityMC::m_bm_cut_x
float m_bm_cut_x
Definition: TBBeamQualityMC.h:28
TBTrack::getCryoHitv
double getCryoHitv() const
Definition: TBTrack.h:67
run
int run(int argc, char *argv[])
Definition: ttree2hdf5.cxx:28
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
CaloClusterMoment::FIRST_ENG_DENS
@ FIRST_ENG_DENS
First Moment in E/V.
Definition: CaloClusterMoment.h:56
x
#define x
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
TBEventInfo::getTableY
float getTableY() const
Definition: TBEventInfo.h:72
TBBeamQualityMC::m_check_trackreco
bool m_check_trackreco
Definition: TBBeamQualityMC.h:38
AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
AthCommonDataStore
Definition: AthCommonDataStore.h:52
TBTrack::getChi2_u
double getChi2_u() const
Definition: TBTrack.h:49
AthAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthAlgorithm.cxx:66
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
lumiFormat.i
int i
Definition: lumiFormat.py:85
python.DecayParser.buf
buf
print ("=> [%s]"cmd)
Definition: DecayParser.py:27
TBTrack::getVslope
double getVslope() const
Definition: TBTrack.h:53
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
TBTrack::getCryoHitu
double getCryoHitu() const
Definition: TBTrack.h:66
chi2
double chi2(TH1 *h0, TH1 *h1)
Definition: comparitor.cxx:523
CaloCluster
Principal data class for CaloCell clusters.
Definition: Calorimeter/CaloEvent/CaloEvent/CaloCluster.h:79
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
TBTrack::getChi2_v
double getChi2_v() const
Definition: TBTrack.h:50
run
Definition: run.py:1
TBBeamQualityMC::m_clusterCollName
std::string m_clusterCollName
Definition: TBBeamQualityMC.h:40
fitman.sy2
sy2
Definition: fitman.py:540
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
TBBeamQualityMC::m_check_veto
bool m_check_veto
Definition: TBBeamQualityMC.h:35
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
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
TBBeamQualityMC::getXcryoYtable
StatusCode getXcryoYtable(float &x, float &y, float &eBeam)
Definition: TBBeamQualityMC.cxx:205
TBTrack::getVintercept
double getVintercept() const
Definition: TBTrack.h:55
TBBeamQualityMC::m_readFileforXcryo
bool m_readFileforXcryo
Get Xcryo and Ytable from a text file.
Definition: TBBeamQualityMC.h:24
AthAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Definition: AthAlgorithm.h:79
fitman.sx2
sx2
Definition: fitman.py:537
a
TList * a
Definition: liststreamerinfos.cxx:10
y
#define y
h
TBEventInfo
Definition: TBEventInfo.h:27
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
CaloClusterMoment::CENTER_LAMBDA
@ CENTER_LAMBDA
Shower depth at Cluster Centroid.
Definition: CaloClusterMoment.h:50
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
CaloCellTimeCorrFiller.filename
filename
Definition: CaloCellTimeCorrFiller.py:24
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
TBEventInfo::getBeamMomentum
float getBeamMomentum() const
Definition: TBEventInfo.h:68
TBBeamQualityMC::m_bm_cut_y
float m_bm_cut_y
Definition: TBBeamQualityMC.h:29
AthAlgorithm::AthAlgorithm
AthAlgorithm()
Default constructor:
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
LArG4H6FrontHit
Definition: LArG4H6FrontHit.h:23
TBBeamQualityMC::m_check_primary
bool m_check_primary
Definition: TBBeamQualityMC.h:32
xAOD::track
@ track
Definition: TrackingPrimitives.h:512
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
TBEventInfo::getCryoX
float getCryoX() const
Definition: TBEventInfo.h:70
AthenaHitsVector
Definition: AthenaHitsVector.h:39
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
fitman.k
k
Definition: fitman.py:528
python.handimod.cc
int cc
Definition: handimod.py:523
TBBeamQualityMC::m_nRun
int m_nRun
Definition: TBBeamQualityMC.h:26
ServiceHandle< ICondSvc >