ATLAS Offline Software
Loading...
Searching...
No Matches
TrigTRTHTHCounter Class Reference

TrigTRTHTHCounter gives back the number of TRT high-threshold hits based on the code of HIPsFilterAlg by Sergey Burdin @migrated from TrigTRTHTHCounter.h & cxx. More...

#include <TrigTRTHTHCounter.h>

Inheritance diagram for TrigTRTHTHCounter:
Collaboration diagram for TrigTRTHTHCounter:

Public Member Functions

 TrigTRTHTHCounter (const std::string &name, ISvcLocator *pSvcLocator)
virtual StatusCode initialize () override
virtual StatusCode execute (const EventContext &ctx) const override
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
virtual bool filterPassed (const EventContext &ctx) const
virtual void setFilterPassed (bool state, const EventContext &ctx) const
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
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.

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

const TRT_IDm_trtHelper {nullptr}
 TRT ID helper.
float m_maxCaloEta {1.7}
Gaudi::Property< float > m_etaHalfWidth {this, "EtaHalfWidth", 0.1 , "subsection of RoI to retrieve fewer TRT hits"}
Gaudi::Property< float > m_phiHalfWidth {this, "PhiHalfWidth", 0.1 , "subsection of RoI to retrieve fewer TRT hits"}
Gaudi::Property< bool > m_doFullScan {this, "doFullScan", false , "Whether to use all RoI"}
Gaudi::Property< float > m_roadWidth {this, "RoadWidth", 4. , "Width of road in mm"}
Gaudi::Property< int > m_nBinCoarse {this, "nBinCoarse", 14 , "Number of coarse bins used while phi centering"}
Gaudi::Property< int > m_nBinFine {this, "nBinFine", 14 , "Number of fine bins used while phi centering"}
Gaudi::Property< float > m_wedgeMinEta {this, "WedgeMinEta", 0 , "Min eta for wedge algorithm"}
Gaudi::Property< float > m_roadMaxEta {this, "RoadMaxEta", 1.06 , "Max eta for road algorithm (barrel only)"}
Gaudi::Property< int > m_wedgeNBin {this, "WedgeNBin", 5 , "Number of fine bins to consider in the wedge algorithm"}
SG::ReadHandleKey< TrigRoiDescriptorCollectionm_roiCollectionKey
SG::ReadHandleKey< InDet::TRT_DriftCircleContainer > m_trtDCContainerKey
SG::WriteHandleKey< xAOD::TrigRNNOutputContainerm_trigRNNOutputKey
ToolHandle< GenericMonitoringToolm_monTool { this, "MonTool", "", "Monitoring tool" }
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default)
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

TrigTRTHTHCounter gives back the number of TRT high-threshold hits based on the code of HIPsFilterAlg by Sergey Burdin @migrated from TrigTRTHTHCounter.h & cxx.

Author
Jiri Masik Jiri..nosp@m.Masi.nosp@m.k@hep.nosp@m..man.nosp@m.chest.nosp@m.er.a.nosp@m.c.uk

Definition at line 42 of file TrigTRTHTHCounter.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ TrigTRTHTHCounter()

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

Definition at line 38 of file TrigTRTHTHCounter.cxx.

39 : AthReentrantAlgorithm(name, pSvcLocator)
40{}

Member Function Documentation

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const
overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64{
65 return 0;
66}

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode TrigTRTHTHCounter::execute ( const EventContext & ctx) const
overridevirtual

Definition at line 63 of file TrigTRTHTHCounter.cxx.

63 {
64
65 ATH_MSG_DEBUG( "Executing " <<name());
66
67 auto trigRNNOutputColl = SG::makeHandle (m_trigRNNOutputKey, ctx);
68 ATH_CHECK(trigRNNOutputColl.record (std::make_unique<xAOD::TrigRNNOutputContainer>(),
69 std::make_unique<xAOD::TrigRNNOutputAuxContainer>()));
70
71 ATH_MSG_DEBUG( "Made WriteHandle " << m_trigRNNOutputKey );
72
73 auto roiCollection = SG::makeHandle(m_roiCollectionKey, ctx);
74 ATH_MSG_DEBUG( "Made handle " << m_roiCollectionKey );
75
76 if (roiCollection->size()==0) {
77 ATH_MSG_DEBUG(" RoI collection size = 0");
78 return StatusCode::SUCCESS;
79 }
80
81 const TrigRoiDescriptor* roiDescriptor = *(roiCollection->begin());
82
83 ATH_MSG_DEBUG(" RoI ID = " << (roiDescriptor)->roiId()
84 << ": Eta = " << (roiDescriptor)->eta()
85 << ", Phi = " << (roiDescriptor)->phi());
86
87
88 std::vector<float> trththits{0,-999,0,-999,-999,-999};
89
90 //Vectors to hold the count of total and HT TRT hits in the coarse bins
91 std::vector<int> count_httrt_c(m_nBinCoarse);
92 std::vector<int> count_tottrt_c(m_nBinCoarse);
93
94 //Vectors to hold the count of total and HT TRT hits in the fine bins
95 std::vector<int> count_httrt(3*m_nBinFine);
96 std::vector<int> count_tottrt(3*m_nBinFine);
97
98 //Vector to hold TRT hits that are within RoI
99 std::vector<TRT_hit> hit;
100
101 //Sanity check of the ROI size
102 double deltaEta= std::abs(roiDescriptor->etaPlus()-roiDescriptor->etaMinus());
103 double deltaPhi= std::abs(CxxUtils::deltaPhi(roiDescriptor->phiPlus(),roiDescriptor->phiMinus()));
104
105 ATH_MSG_DEBUG( "roiDescriptor->etaPlus() in TrigTRTHTHCounter:"<<roiDescriptor->etaPlus());
106 ATH_MSG_DEBUG( "roiDescriptor->etaMinus() in TrigTRTHTHCounter:"<<roiDescriptor->etaMinus());
107 ATH_MSG_DEBUG( "deltaEta in TrigTRTHTHCounter:"<<deltaEta);
108
109 float phiTolerance = 0.001;
110 float etaTolerance = 0.001;
111
112 if((m_etaHalfWidth - deltaEta/2.) > etaTolerance)
113 ATH_MSG_WARNING ( "ROI eta range too small : " << deltaEta);
114
115 if((m_phiHalfWidth - deltaPhi/2.) > phiTolerance)
116 ATH_MSG_WARNING ( "ROI phi range too small : " << deltaPhi);
117
118 float coarseWedgeHalfWidth = m_phiHalfWidth/m_nBinCoarse;
119 float fineWedgeHalfWidth = coarseWedgeHalfWidth/m_nBinFine;
120
121 //Code will only proceed if the RoI eta is not too large; used to limit rate from endcap
122 if ( std::abs(roiDescriptor->eta())<=m_maxCaloEta ){
123
124 SG::ReadHandle<InDet::TRT_DriftCircleContainer> trtDCC(m_trtDCContainerKey, ctx);
125 ATH_MSG_DEBUG( "Made handle " << m_trtDCContainerKey );
126
127 if (trtDCC->size() == 0){
128 return StatusCode::SUCCESS; // Exit early if there are no tracks
129 }
130
131 InDet::TRT_DriftCircleContainer::const_iterator trtdriftContainerItr = trtDCC->begin();
132 InDet::TRT_DriftCircleContainer::const_iterator trtdriftContainerItrE = trtDCC->end();
133
134 for (; trtdriftContainerItr != trtdriftContainerItrE; ++trtdriftContainerItr) {
135
136 InDet::TRT_DriftCircleCollection::const_iterator trtItr = (*trtdriftContainerItr)->begin();
137 InDet::TRT_DriftCircleCollection::const_iterator trtEnd = (*trtdriftContainerItr)->end();
138
139 for(; trtItr!=trtEnd; ++trtItr){
140
141 // find out which detector element the hit belongs to
142 const InDetDD::TRT_BaseElement *det = (*trtItr)->detectorElement();
143 Identifier ID = (*trtItr)->identify();
144 const Amg::Vector3D& strawcenter = det->strawCenter(m_trtHelper->straw(ID));
145
146 bool hth = false;
147 float hphi = strawcenter.phi();
148 float heta = strawcenter.eta();
149 float R = strawcenter.perp();
150
151 if ((*trtItr)->highLevel()) hth = true;
152
153 //hit needs to be stored
154 if(hth_eta_match((roiDescriptor)->eta(), heta, m_etaHalfWidth))
155 hit.push_back(make_hit(hphi,R,hth));
156
157
158 //First, define coarse wedges in phi, and count the TRT hits in these wedges
159 int countbin=0;
160 if(std::abs(CxxUtils::deltaPhi(hphi, static_cast<float>(roiDescriptor->phi()))) < 0.1){
161 float startValue = roiDescriptor->phi() - m_phiHalfWidth + coarseWedgeHalfWidth;
162 float endValue = roiDescriptor->phi() + m_phiHalfWidth;
163 float increment = 2*coarseWedgeHalfWidth;
164 for(float roibincenter = startValue; roibincenter < endValue; roibincenter += increment){
165 if (std::abs(CxxUtils::deltaPhi(hphi,roibincenter))<=coarseWedgeHalfWidth && hth_eta_match((roiDescriptor)->eta(), heta, m_etaHalfWidth)) {
166 if(hth) count_httrt_c.at(countbin) += 1.;
167 count_tottrt_c.at(countbin) += 1.;
168 break; //the hit has been assigned to one of the coarse wedges, so no need to continue the for loop
169 }
170 countbin++;
171 }
172 }
173 ATH_MSG_VERBOSE ( "timeOverThreshold=" << (*trtItr)->timeOverThreshold()
174 << " highLevel=" << (*trtItr)->highLevel()
175 << " rawDriftTime=" << (*trtItr)->rawDriftTime()
176 << " barrel_ec=" << m_trtHelper->barrel_ec(ID)
177 << " phi_module=" << m_trtHelper->phi_module(ID)
178 << " layer_or_wheel=" << m_trtHelper->layer_or_wheel(ID)
179 << " straw_layer=" << m_trtHelper->straw_layer(ID)
180 << " straw=" << m_trtHelper->straw(ID)
181 << " scR=" << det->strawCenter(m_trtHelper->straw(ID)).perp()
182 << " scPhi=" << hphi
183 << " scEta=" << heta);
184 } // end of loop over TRT drift circle collection
185 } //end of loop over TRT drift circle container
186 } //end of check to see if RoI eta is not too large
187
188 //Now figure out which of the coarse bins in phi has the max number of HT TRT hits
189 int maxHits = 0; //used to keep track of the max number of HT TRT hits in a coarse bin
190 int dist = 0; //used to keep track of which coarse bin has the max number of HT TRT hits
191
192 for (size_t iw = 0; iw < count_httrt_c.size(); iw++){
193 if(maxHits <= count_httrt_c[iw]){
194 maxHits = count_httrt_c[iw];
195 dist = iw;
196 }
197 }
198
199 //Value of dist can range between 0 and (m_nBinCoarse-1)
200 float center_pos_phi=roiDescriptor->phi()+(2*dist+1-m_nBinCoarse)*coarseWedgeHalfWidth;
201
202 //Now, define fine wedges in phi, centered around the best coarse wedge, and count the TRT hits in these fine wedges
203 for(size_t v=0;v<hit.size();v++){
204 int countbin=0;
205 if(std::abs(CxxUtils::deltaPhi(hit[v].phi, center_pos_phi)) < 0.01){
206 float startValue = center_pos_phi - 3*coarseWedgeHalfWidth + fineWedgeHalfWidth;
207 float endValue = center_pos_phi + 3*coarseWedgeHalfWidth;
208 float increment = 2*fineWedgeHalfWidth;
209 for(float roibincenter = startValue; roibincenter < endValue; roibincenter += increment){
210 if (std::abs(CxxUtils::deltaPhi(hit[v].phi,roibincenter))<=fineWedgeHalfWidth) {
211 if(hit[v].isHT) count_httrt.at(countbin) += 1.;
212 count_tottrt.at(countbin) += 1.;
213 break; //the hit has been assigned to one of the fine wedges, so no need to continue the for loop
214 }
215 countbin++;
216 }
217 }
218 }
219
220 //Now figure out which of the fine bins in phi has the max number of HT TRT hits
221 maxHits = 0; //used to keep track of the max number of HT TRT hits in a fine bin
222 dist = 0; //used to keep track of which fine bin has the max number of HT TRT hits
223
224 for (size_t iw = 0; iw < count_httrt.size(); iw++){
225 if(maxHits <= count_httrt[iw]){
226 maxHits = count_httrt[iw];
227 dist = iw;
228 }
229 }
230
231 //Value of dist can range between 0 and (3*m_nBinFine-1)
232 center_pos_phi+=(2*dist+1-3*m_nBinFine)*fineWedgeHalfWidth;
233
234 //Count the number of total and HT TRT hits for the road algorithm
235 int trthit=0, trthit_ht=0;
236 for(size_t v=0;v<hit.size();v++){
237 if (dist2COR(hit[v].R,hit[v].phi,center_pos_phi)<=m_roadWidth){
238 if(hit[v].isHT) trthit_ht+=1;
239 trthit+=1;
240 }
241 }
242
243 if (trthit!=0&&(std::abs(roiDescriptor->eta())<=m_roadMaxEta)){
244 trththits[0] = trthit_ht;
245 trththits[1] = (float)trthit_ht/trthit;
246 }
247
248 //Count the number of total and HT TRT hits for the wedge algorithm
249 trthit = 0;
250 trthit_ht = 0;
251
252 for (int k=0;k<(2*m_wedgeNBin+1);k++){
253 int binNumber = dist+k-m_wedgeNBin;
254 //Even though we are supposed to use 2*m_wedgeNBin+1 fine bins,
255 //need to make sure that binNumber is sensible
256 if(binNumber >= 0 && binNumber < (int)count_httrt.size()){
257 trthit += count_tottrt.at(binNumber);
258 trthit_ht += count_httrt.at(binNumber);
259 }
260 }
261
262 if (trthit!=0&&(std::abs(roiDescriptor->eta())>=m_wedgeMinEta)){
263 trththits[2] = trthit_ht;
264 trththits[3] = (float)trthit_ht/trthit;
265 }
266
267 trththits[4]=roiDescriptor->eta();
268 trththits[5]=roiDescriptor->phi();
269
270 ATH_MSG_VERBOSE ( "trthits with road algorithm : " << trththits[0]);
271 ATH_MSG_VERBOSE ( "fHT with road algorithm : " << trththits[1]);
272 ATH_MSG_VERBOSE ( "trthits with wedge algorithm : " << trththits[2]);
273 ATH_MSG_VERBOSE ( "fHT with wedge algorithm : " << trththits[3]);
274 ATH_MSG_VERBOSE ( "ROI eta : " << trththits[4]);
275
276 //Writing to xAOD
277 xAOD::TrigRNNOutput* rnnOutput = new xAOD::TrigRNNOutput();
278 trigRNNOutputColl->push_back(rnnOutput);
279 rnnOutput->setRnnDecision(trththits);
280
281 ATH_MSG_DEBUG("REGTEST: returning an xAOD::TrigRNNOutputContainer with size "<< trigRNNOutputColl->size() << ".");
282
283 return StatusCode::SUCCESS;
284}
Scalar eta() const
pseudorapidity method
Scalar deltaPhi(const MatrixBase< Derived > &vec) const
Scalar phi() const
phi method
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_WARNING(x)
#define ATH_MSG_DEBUG(x)
Athena::TPCnvVers::Current TrigRoiDescriptor
float dist2COR(float R, float phi1, float phi2)
bool hth_eta_match(float caleta, float trteta, float etaWindow)
TRT_hit make_hit(float phi, float R, bool isHT)
virtual double etaMinus() const override final
gets eta at zMinus
virtual double etaPlus() const override final
gets eta at zedPlus
virtual double phi() const override final
Methods to retrieve data members.
virtual double phiMinus() const override final
gets phiMinus
virtual double eta() const override final
virtual double phiPlus() const override final
gets phiPlus
Gaudi::Property< int > m_wedgeNBin
Gaudi::Property< float > m_etaHalfWidth
Gaudi::Property< int > m_nBinCoarse
Gaudi::Property< float > m_roadWidth
SG::ReadHandleKey< TrigRoiDescriptorCollection > m_roiCollectionKey
const TRT_ID * m_trtHelper
TRT ID helper.
Gaudi::Property< float > m_phiHalfWidth
SG::WriteHandleKey< xAOD::TrigRNNOutputContainer > m_trigRNNOutputKey
Gaudi::Property< float > m_roadMaxEta
Gaudi::Property< float > m_wedgeMinEta
SG::ReadHandleKey< InDet::TRT_DriftCircleContainer > m_trtDCContainerKey
Gaudi::Property< int > m_nBinFine
void setRnnDecision(const std::vector< float > &d)
Eigen::Matrix< double, 3, 1 > Vector3D
double R(const INavigable4Momentum *p1, const double v_eta, const double v_phi)
T deltaPhi(T phiA, T phiB)
Return difference phiA - phiB in range [-pi, pi].
Definition phihelper.h:42
double deltaEta(const I4Momentum &p1, const I4Momentum &p2)
Computes efficiently .
Definition P4Helpers.h:66
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
ID
//////////////////////////////////////// JetAlgorithmType::ID defines most common physics jet finding...
TrigRNNOutput_v2 TrigRNNOutput
Define the latest version of the RingerRings class.
setTeId setLumiBlock roiId

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::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 & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

90{
91 // If we didn't find any symlinks to add, just return the collection
92 // from the base class. Otherwise, return the extended collection.
93 if (!m_extendedExtraObjects.empty()) {
95 }
97}
An algorithm that can be simultaneously executed in multiple threads.

◆ filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx) const
inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

96 {
97 return execState( ctx ).filterPassed();
98 }
virtual bool filterPassed(const EventContext &ctx) const

◆ initialize()

StatusCode TrigTRTHTHCounter::initialize ( )
overridevirtual

Definition at line 42 of file TrigTRTHTHCounter.cxx.

43{
44 ATH_MSG_DEBUG ( "Initialising this TrigTRTHTHCounter: " << name());
45
46 // Get a TRT identifier helper
47 ATH_CHECK( detStore()->retrieve(m_trtHelper, "TRT_ID") );
48
49 if (!m_doFullScan){
50 ATH_MSG_INFO ( "PhiHalfWidth: " << m_phiHalfWidth << " EtaHalfWidth: "<< m_etaHalfWidth);
51 } else {
52 ATH_MSG_INFO ( "FullScan mode");
53 }
54
55 ATH_CHECK( m_roiCollectionKey.initialize() );
56 ATH_CHECK( m_trtDCContainerKey.initialize() );
57 ATH_CHECK( m_trigRNNOutputKey.initialize() );
58
59 return StatusCode::SUCCESS;
60}
#define ATH_MSG_INFO(x)
const ServiceHandle< StoreGateSvc > & detStore() const
Gaudi::Property< bool > m_doFullScan

◆ inputHandles()

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

◆ isClonable()

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

bool AthCommonMsg< Gaudi::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< Gaudi::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< Gaudi::Algorithm > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
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)

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state,
const EventContext & ctx ) const
inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100 {
102 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const

◆ sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx)
overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

77{
78 return BaseAlg::sysExecute (ctx);
79}

◆ sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::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< Gaudi::Algorithm > >.

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

107 {
109
110 if (sc.isFailure()) {
111 return sc;
112 }
113
114 ServiceHandle<ICondSvc> cs("CondSvc",name());
115 for (auto h : outputHandles()) {
116 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
117 // do this inside the loop so we don't create the CondSvc until needed
118 if ( cs.retrieve().isFailure() ) {
119 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
120 return StatusCode::SUCCESS;
121 }
122 if (cs->regHandle(this,*h).isFailure()) {
124 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
125 << " with CondSvc");
126 }
127 }
128 }
129 return sc;
130}
#define ATH_MSG_ERROR(x)
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::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< Gaudi::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) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }

Member Data Documentation

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_doFullScan

Gaudi::Property<bool> TrigTRTHTHCounter::m_doFullScan {this, "doFullScan", false , "Whether to use all RoI"}
private

Definition at line 59 of file TrigTRTHTHCounter.h.

59{this, "doFullScan", false , "Whether to use all RoI"};

◆ m_etaHalfWidth

Gaudi::Property<float> TrigTRTHTHCounter::m_etaHalfWidth {this, "EtaHalfWidth", 0.1 , "subsection of RoI to retrieve fewer TRT hits"}
private

Definition at line 57 of file TrigTRTHTHCounter.h.

57{this, "EtaHalfWidth", 0.1 , "subsection of RoI to retrieve fewer TRT hits"};

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_maxCaloEta

float TrigTRTHTHCounter::m_maxCaloEta {1.7}
private

Definition at line 54 of file TrigTRTHTHCounter.h.

54{1.7};

◆ m_monTool

ToolHandle< GenericMonitoringTool > TrigTRTHTHCounter::m_monTool { this, "MonTool", "", "Monitoring tool" }
private

Definition at line 83 of file TrigTRTHTHCounter.h.

83{ this, "MonTool", "", "Monitoring tool" };

◆ m_nBinCoarse

Gaudi::Property<int> TrigTRTHTHCounter::m_nBinCoarse {this, "nBinCoarse", 14 , "Number of coarse bins used while phi centering"}
private

Definition at line 61 of file TrigTRTHTHCounter.h.

61{this, "nBinCoarse", 14 , "Number of coarse bins used while phi centering"};

◆ m_nBinFine

Gaudi::Property<int> TrigTRTHTHCounter::m_nBinFine {this, "nBinFine", 14 , "Number of fine bins used while phi centering"}
private

Definition at line 62 of file TrigTRTHTHCounter.h.

62{this, "nBinFine", 14 , "Number of fine bins used while phi centering"};

◆ m_phiHalfWidth

Gaudi::Property<float> TrigTRTHTHCounter::m_phiHalfWidth {this, "PhiHalfWidth", 0.1 , "subsection of RoI to retrieve fewer TRT hits"}
private

Definition at line 58 of file TrigTRTHTHCounter.h.

58{this, "PhiHalfWidth", 0.1 , "subsection of RoI to retrieve fewer TRT hits"};

◆ m_roadMaxEta

Gaudi::Property<float> TrigTRTHTHCounter::m_roadMaxEta {this, "RoadMaxEta", 1.06 , "Max eta for road algorithm (barrel only)"}
private

Definition at line 64 of file TrigTRTHTHCounter.h.

64{this, "RoadMaxEta", 1.06 , "Max eta for road algorithm (barrel only)"};

◆ m_roadWidth

Gaudi::Property<float> TrigTRTHTHCounter::m_roadWidth {this, "RoadWidth", 4. , "Width of road in mm"}
private

Definition at line 60 of file TrigTRTHTHCounter.h.

60{this, "RoadWidth", 4. , "Width of road in mm"};

◆ m_roiCollectionKey

SG::ReadHandleKey<TrigRoiDescriptorCollection> TrigTRTHTHCounter::m_roiCollectionKey
private
Initial value:
{ this,
"RoIs",
"rois",
"input RoI Collection name"}

Definition at line 67 of file TrigTRTHTHCounter.h.

67 { this,
68 "RoIs", // property name
69 "rois", // default value of StoreGate key
70 "input RoI Collection name"};

◆ m_trigRNNOutputKey

SG::WriteHandleKey<xAOD::TrigRNNOutputContainer> TrigTRTHTHCounter::m_trigRNNOutputKey
private
Initial value:
{ this,
"RNNOutputName",
"RNNOutput",
"output RNN container name "}

Definition at line 78 of file TrigTRTHTHCounter.h.

78 { this,
79 "RNNOutputName", // property name
80 "RNNOutput", // default value of StoreGate key
81 "output RNN container name "};

◆ m_trtDCContainerKey

SG::ReadHandleKey<InDet::TRT_DriftCircleContainer> TrigTRTHTHCounter::m_trtDCContainerKey
private
Initial value:
{this,
"TRT_DC_ContainerName",
"TRT_DriftCircleContainer",
"TRT DriftCircle Container"}

Definition at line 72 of file TrigTRTHTHCounter.h.

72 {this,
73 "TRT_DC_ContainerName",
74 "TRT_DriftCircleContainer",
75 "TRT DriftCircle Container"};

◆ m_trtHelper

const TRT_ID* TrigTRTHTHCounter::m_trtHelper {nullptr}
private

TRT ID helper.

Definition at line 52 of file TrigTRTHTHCounter.h.

52{nullptr};

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_wedgeMinEta

Gaudi::Property<float> TrigTRTHTHCounter::m_wedgeMinEta {this, "WedgeMinEta", 0 , "Min eta for wedge algorithm"}
private

Definition at line 63 of file TrigTRTHTHCounter.h.

63{this, "WedgeMinEta", 0 , "Min eta for wedge algorithm"};

◆ m_wedgeNBin

Gaudi::Property<int> TrigTRTHTHCounter::m_wedgeNBin {this, "WedgeNBin", 5 , "Number of fine bins to consider in the wedge algorithm"}
private

Definition at line 65 of file TrigTRTHTHCounter.h.

65{this, "WedgeNBin", 5 , "Number of fine bins to consider in the wedge algorithm"};

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