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

Contains methods to do the conversion from bytestream format to TileRawChannel stored in a TileRawChannelContainer and vice-versa making use of the TileRawChannelContByteStreamTool. More...

#include <TileRawChannelContByteStreamCnv.h>

Inheritance diagram for TileRawChannelContByteStreamCnv:
Collaboration diagram for TileRawChannelContByteStreamCnv:

Public Types

typedef TileRawChannelContByteStreamTool BYTESTREAMTOOL

Public Member Functions

 TileRawChannelContByteStreamCnv (ISvcLocator *svcloc)
virtual StatusCode initialize () override
virtual StatusCode createObjConst (IOpaqueAddress *pAddr, DataObject *&pObj) const override
 Create the transient representation of an object.
virtual StatusCode createRepConst (DataObject *pObj, IOpaqueAddress *&pAddr) const override
 Convert the transient object to the requested representation.
virtual StatusCode finalize () override
virtual long repSvcType () const override
 Storage type and class ID.
virtual StatusCode createObj (IOpaqueAddress *pAddress, DataObject *&refpObject) override final
 Create the transient representation of an object.
virtual StatusCode createRep (DataObject *pObject, IOpaqueAddress *&refpAddress) override final
 Convert the transient object to the requested representation.
const std::string & name () const
void addTrack (TIDA::Track *t)
const std::vector< TIDA::Track * > & tracks () const
void clear ()
void selectTracks (const TrigInDetTrackCollection *trigtracks)
void selectTracks (const Rec::TrackParticleContainer *trigtracks)
void selectTracks (const Analysis::MuonContainer *muontracks)
bool msgLvl (const MSG::Level lvl) const
 Test the output level.
MsgStream & msg () const
 The standard message stream.
MsgStream & msg (const MSG::Level lvl) const
 The standard message stream.
void setLevel (MSG::Level lvl)
 Change the current logging level.

Static Public Member Functions

static long storageType ()
static const CLIDclassID ()

Protected Member Functions

double phiCorr (double phi)
void ipCorr (double d0, double z0, double &d0c, double &z0c, double phi0, double eta, double pt)

Protected Attributes

double m_beamX
double m_beamY
double m_beamZ
std::vector< TIDA::Track * > m_tracks

Private Member Functions

void initMessaging () const
 Initialize our message level and MessageSvc.

Private Attributes

ToolHandle< TileRawChannelContByteStreamToolm_tool
ServiceHandle< IByteStreamCnvSvcm_byteStreamCnvSvc
ServiceHandle< IROBDataProviderSvcm_robSvc
 Pointer to IROBDataProviderSvc.
ToolHandle< TileROD_Decoderm_decoder
 Pointer to TileROD_Decoder.
const TileHid2RESrcIDm_hid2re
 Pointer to TileHid2RESrcID.
Athena::RecyclableDataQueue< TileMutableRawChannelContainer > m_queue ATLAS_THREAD_SAFE
 Queue of data objects to recycle.
std::string m_name
std::string m_nm
 Message source name.
boost::thread_specific_ptr< MsgStream > m_msg_tls
 MsgStream instance (a std::cout like with print-out levels)
std::atomic< IMessageSvc * > m_imsg { nullptr }
 MessageSvc pointer.
std::atomic< MSG::Level > m_lvl { MSG::NIL }
 Current logging level.

Detailed Description

Contains methods to do the conversion from bytestream format to TileRawChannel stored in a TileRawChannelContainer and vice-versa making use of the TileRawChannelContByteStreamTool.

Definition at line 39 of file TileRawChannelContByteStreamCnv.h.

Member Typedef Documentation

◆ BYTESTREAMTOOL

Constructor & Destructor Documentation

◆ TileRawChannelContByteStreamCnv()

TileRawChannelContByteStreamCnv::TileRawChannelContByteStreamCnv ( ISvcLocator * svcloc)

Definition at line 36 of file TileRawChannelContByteStreamCnv.cxx.

37 : AthConstConverter(storageType(), classID(), svcloc, "TileRawChannelContByteStreamCnv")
38 , m_tool("TileRawChannelContByteStreamTool")
39 , m_byteStreamCnvSvc("ByteStreamCnvSvc", name())
40 , m_robSvc("ROBDataProviderSvc", name())
41 , m_decoder("TileROD_Decoder")
42 , m_hid2re(0)
43{
44}
AthConstConverter(long storage_type, const CLID &class_type, ISvcLocator *svc, const std::string &name)
const std::string & name() const
ToolHandle< TileROD_Decoder > m_decoder
Pointer to TileROD_Decoder.
ServiceHandle< IROBDataProviderSvc > m_robSvc
Pointer to IROBDataProviderSvc.
ServiceHandle< IByteStreamCnvSvc > m_byteStreamCnvSvc
ToolHandle< TileRawChannelContByteStreamTool > m_tool
const TileHid2RESrcID * m_hid2re
Pointer to TileHid2RESrcID.

Member Function Documentation

◆ addTrack()

void Converter::addTrack ( TIDA::Track * t)
inlineinherited

Definition at line 45 of file Converter.h.

45 {
46 m_tracks.push_back(t);
47 }
std::vector< TIDA::Track * > m_tracks
Definition Converter.h:321

◆ classID()

const CLID & TileRawChannelContByteStreamCnv::classID ( )
static

◆ clear()

void Converter::clear ( )
inlineinherited

Definition at line 53 of file Converter.h.

53{m_tracks.clear();}

◆ createObj()

StatusCode AthConstConverter::createObj ( IOpaqueAddress * pAddress,
DataObject *& refpObject )
finaloverridevirtualinherited

Create the transient representation of an object.

Non-const version; just calls the const version.

Definition at line 44 of file AthConstConverter.cxx.

46{
47 return createObjConst (pAddress, refpObject);
48}
virtual StatusCode createObjConst(IOpaqueAddress *pAddress, DataObject *&refpObject) const
Create the transient representation of an object.

◆ createObjConst()

StatusCode TileRawChannelContByteStreamCnv::createObjConst ( IOpaqueAddress * pAddress,
DataObject *& refpObject ) const
overridevirtual

Create the transient representation of an object.

The default implementation is a no-op.

Reimplemented from AthConstConverter.

Definition at line 70 of file TileRawChannelContByteStreamCnv.cxx.

71{
72 ATH_MSG_DEBUG(" Executing createObj method ");
73
74 ByteStreamAddress *pRE_Addr;
75 pRE_Addr = dynamic_cast<ByteStreamAddress*>(pAddr);
76
77 if(!pRE_Addr) {
78 ATH_MSG_ERROR(" Can not cast to ByteStreamAddress ");
79 return StatusCode::FAILURE;
80 }
81
82 const std::string containerName(*(pRE_Addr->par()));
83 bool isTMDB(containerName == std::string("MuRcvRawChCnt"));
84
85 uint32_t newrob = 0x0;
86
87
88 std::vector<uint32_t> robid(1);
89 robid[0] = 0;
90 std::vector<const ROBDataProviderSvc::ROBF*> robf;
91
92 TileMutableRawChannelContainer* cont = m_queue.get (true);
93 ATH_CHECK( cont->status() );
94
95 if (isTMDB) {
97 } else {
99 }
100
101 std::unordered_map<uint32_t,int> bsflags;
103
104 // iterate over all collections in a container and fill them
105 for (IdentifierHash hash : cont->GetAllCurrentHashes()) {
106 TileRawChannelCollection* rawChannelCollection = cont->indexFindPtr (hash);
107 rawChannelCollection->clear();
108 TileRawChannelCollection::ID collID = rawChannelCollection->identify();
109
110 // find ROB
111 if (isTMDB) {
112 newrob = m_hid2re->getRobFromTileMuRcvFragID(collID);
113 } else {
114 newrob = m_hid2re->getRobFromFragID(collID);
115 }
116
117 if (newrob != robid[0]) {
118 robid[0] = newrob;
119 robf.clear();
120 m_robSvc->getROBData(Gaudi::Hive::currentContext(), robid, robf);
121 }
122
123 // unpack ROB data
124 if (robf.size() > 0 ) {
125 if (isTMDB) {// reid for TMDB 0x5x010x
126 m_decoder->fillCollection_TileMuRcv_RawChannel(robf[0], *rawChannelCollection);
127 } else {
128 m_decoder->fillCollection(robf[0], *rawChannelCollection, cont);
129 }
130 flag = cont->get_bsflags();
131 auto result = bsflags.insert(std::pair<uint32_t, int>(flag, 1));
132 if (result.second == false) result.first->second++;
133 } else {
134 ATH_MSG_DEBUG( "ROB for " << ((isTMDB)?"TMDB ":"") << "drawer 0x" << MSG::hex << collID << MSG::dec << " not found in BS" );
136 rawChannelCollection->setFragGlobalCRC(status);
137 ATH_MSG_DEBUG( "Status for " << ((isTMDB)?"TMDB ":"") << "drawer 0x" << MSG::hex << collID << " is 0x" << status << MSG::dec);
138 }
139 }
140
141 if (bsflags.size() > 1) {
142 int n=0;
143 for (const auto & elem : bsflags) {
144 if (elem.second > n) {
145 n = elem.second;
146 flag = elem.first;
147 }
148 }
149
150 if (flag != cont->get_bsflags()) {
151
152 uint32_t unit = ((flag & 0xC0000000) >> 30);
153 if ((flag & 0x30000000) < 0x30000000) unit += TileRawChannelUnit::OnlineOffset; // Online units in real data
154
155 ATH_MSG_DEBUG( "Changing units for " << ((isTMDB)?"TMDB ":"") << "RawChannelContainer from "
156 << cont->get_unit() << " to " << unit << MSG::hex
157 << " and BS flags from 0x" << cont->get_bsflags() << " to 0x" << flag << MSG::dec);
158
160 cont->set_bsflags(flag);
161 }
162 }
163
164 ATH_MSG_DEBUG( "Creating raw channel container: " << containerName );
165
166 TileRawChannelContainer* basecont = cont;
167 pObj = SG::asStorable( basecont );
168
169
170 return StatusCode::SUCCESS;
171}
const PlainObject unit() const
This is a plugin that makes Eigen look like CLHEP & defines some convenience methods.
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x)
#define ATH_MSG_DEBUG(x)
TileMutableDataContainer< TileRawChannelContainer > TileMutableRawChannelContainer
Athena::TPCnvVers::Old Athena::TPCnvVers::Old Athena::TPCnvVers::Current Athena::TPCnvVers::Old Athena::TPCnvVers::Current TileRawChannelContainer
virtual std::vector< IdentifierHash > GetAllCurrentHashes() const override final
Returns a collection of all hashes availiable in this IDC.
StatusCode status() const
Return the error status from the constructors.
Collection * indexFindPtr(IdentifierHash hash)
Look up a (non-const) collection via hash.
void setFragGlobalCRC(uint32_t globalcrc)
Various set methods.
void clear()
Clear everything for next event.
void set_bsflags(uint32_t bsflags)
uint32_t get_bsflags() const
DataObject * asStorable(SG::DataObjectSharedPtr< T > pObject)
bool flag
Definition master.py:29
status
Definition merge.py:16
setEventNumber uint32_t

◆ createRep()

StatusCode AthConstConverter::createRep ( DataObject * pObject,
IOpaqueAddress *& refpAddress )
finaloverridevirtualinherited

Convert the transient object to the requested representation.

Non-const version; just calls the const version.

Definition at line 56 of file AthConstConverter.cxx.

58{
59 return createRepConst (pObject, refpAddress);
60}
virtual StatusCode createRepConst(DataObject *pObject, IOpaqueAddress *&refpAddress) const
Convert the transient object to the requested representation.

◆ createRepConst()

StatusCode TileRawChannelContByteStreamCnv::createRepConst ( DataObject * pObject,
IOpaqueAddress *& refpAddress ) const
overridevirtual

Convert the transient object to the requested representation.

The default implementation is a no-op.

Reimplemented from AthConstConverter.

Definition at line 173 of file TileRawChannelContByteStreamCnv.cxx.

174{
175 // convert TileRawChannels in the container into ByteStream
176
177 ATH_MSG_DEBUG(" Executing createRep method ");
178
179 // uint32_t runnum = re->header()->specific_part()->run_no();
180
181 std::string key("Tile");
182
183 // get Full Event Assembler
184 FullEventAssembler<TileHid2RESrcID> *fea = 0;
185 ATH_CHECK( m_byteStreamCnvSvc->getFullEventAssembler(fea, key) );
186
187 // create TileRawChannelContainer
188 TileRawChannelContainer* rccont(0) ;
189 SG::fromStorable(pObj, rccont );
190 if(!rccont){
191 ATH_MSG_ERROR(" Can not cast to TileRawChannelContainer ");
192 return StatusCode::FAILURE;
193 }
194
195 std::string name = pObj->registry()->name();
196
197 if ( pAddr != nullptr ) pAddr->release();
198 ByteStreamAddress* addr = new ByteStreamAddress(classID(), name, "");
199 pAddr = addr;
200 pAddr->addRef();
201
202 ATH_CHECK( m_tool->convert(rccont, fea) );
203
204 return StatusCode::SUCCESS;
205}
bool fromStorable(DataObject *pDObj, T *&pTrans, bool quiet=false, IRegisterTransient *irt=0, bool isConst=true)

◆ finalize()

StatusCode TileRawChannelContByteStreamCnv::finalize ( )
overridevirtual

Definition at line 207 of file TileRawChannelContByteStreamCnv.cxx.

208{
209 return Converter::finalize();
210}

◆ initialize()

StatusCode TileRawChannelContByteStreamCnv::initialize ( )
overridevirtual

Definition at line 50 of file TileRawChannelContByteStreamCnv.cxx.

50 {
51
52 ATH_MSG_DEBUG(" initialize ");
53
54 ATH_CHECK(Converter::initialize());
55
56 // Get ByteStreamCnvSvc
57 ATH_CHECK( m_byteStreamCnvSvc.retrieve() );
58
59 // retrieve Tool
60 ATH_CHECK( m_decoder.retrieve() );
61 m_hid2re = m_decoder->getHid2re();
62
63 ATH_CHECK( m_tool.retrieve() );
64
65 ATH_CHECK( m_robSvc.retrieve() );
66
67 return StatusCode::SUCCESS;
68}

◆ initMessaging()

void AthMessaging::initMessaging ( ) const
privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

40{
42 // If user did not set an explicit level, set a default
43 if (m_lvl == MSG::NIL) {
44 m_lvl = m_imsg ?
45 static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
46 MSG::INFO;
47 }
48}
std::string m_nm
Message source name.
std::atomic< IMessageSvc * > m_imsg
MessageSvc pointer.
std::atomic< MSG::Level > m_lvl
Current logging level.
IMessageSvc * getMessageSvc(bool quiet=false)

◆ ipCorr()

void Converter::ipCorr ( double d0,
double z0,
double & d0c,
double & z0c,
double phi0,
double eta,
double pt )
inlineprotectedinherited

Definition at line 257 of file Converter.h.

257 {
258
259 double sn = sin(phi0);
260 double cs = cos(phi0);
261 double sd0 = (d0 != 0 ? d0/fabs(d0) : 1);
262 double spt = (pt != 0 ? pt/fabs(pt) : 1);
263
264 if (fabs(pt) >= 1*CLHEP::TeV) {
265
266 d0c = d0 + m_beamX*sn - m_beamY*cs;
267
268 } else {
269
270 double rc = fabs(pt)*15.0/(9.0*1.042);
271
272 double xc = (fabs(d0)-spt*sd0*rc)*cos(phi0+M_PI/2*sd0) - m_beamX;
273 double yc = (fabs(d0)-spt*sd0*rc)*sin(phi0+M_PI/2*sd0) - m_beamY;
274
275 double newphi;
276 double xd01,yd01,xd02,yd02;
277
278 if (xc == 0) {
279 xd01 = 0; yd01 = rc+yc;
280 xd02 = 0; yd02 = -rc+yc;
281 } else {
282 xd01 = xc+yc/xc*yc+sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd01 = yc/xc*xd01;
283 xd02 = xc+yc/xc*yc-sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd02 = yc/xc*xd02;
284 }
285
286 double r1 = sqrt(xd01*xd01+yd01*yd01);
287 double r2 = sqrt(xd02*xd02+yd02*yd02);
288
289 double phiV;
290
291 if (r1 < r2)
292 phiV = atan2(yd01,xd01);
293 else
294 phiV = atan2(yd02,xd02);
295
296 double phi1 = phiCorr(phiV+M_PI/2);
297 double phi2 = phiCorr(phiV-M_PI/2);
298
299 if (fabs(phiCorr(phi1-phi0))<=fabs(phiCorr(phi2-phi0)))
300 newphi = phi1;
301 else
302 newphi = phi2;
303
304 d0c = fabs(sqrt(xc*xc+yc*yc)-rc)*sin(phiV-newphi);
305
306 double theta=2*atan2(exp(-eta),1);
307 double theta_save=theta;
308 theta = theta - (1+spt)/2*M_PI;
309 if (theta>0) theta = theta_save;
310
311 double deltaz0= -spt*rc/tan(theta)*phiCorr(phiCorr(newphi)-phiCorr(phi0));
312 z0c = z0 + deltaz0;
313
314 }
315 }
#define M_PI
Scalar eta() const
pseudorapidity method
Scalar theta() const
theta method
static Double_t rc
constexpr int pow(int base, int exp) noexcept
double m_beamX
Definition Converter.h:320
double phiCorr(double phi)
Definition Converter.h:248
double m_beamY
Definition Converter.h:320

◆ msg() [1/2]

MsgStream & AthMessaging::msg ( ) const
inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 163 of file AthMessaging.h.

164{
165 MsgStream* ms = m_msg_tls.get();
166 if (!ms) {
167 if (!m_initialized.test_and_set()) initMessaging();
168 ms = new MsgStream(m_imsg,m_nm);
169 m_msg_tls.reset( ms );
170 }
171
172 ms->setLevel (m_lvl);
173 return *ms;
174}
boost::thread_specific_ptr< MsgStream > m_msg_tls
MsgStream instance (a std::cout like with print-out levels)
void initMessaging() const
Initialize our message level and MessageSvc.

◆ msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level lvl) const
inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 178 of file AthMessaging.h.

179{ return msg() << lvl; }
MsgStream & msg() const
The standard message stream.

◆ msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl) const
inlineinherited

Test the output level.

Parameters
lvlThe message level to test against
Returns
boolean Indicating if messages at given level will be printed
Return values
trueMessages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

152{
153 if (m_lvl <= lvl) {
154 msg() << lvl;
155 return true;
156 } else {
157 return false;
158 }
159}

◆ name()

const std::string & AthConstConverter::name ( ) const
inlineinherited

Definition at line 67 of file AthConstConverter.h.

67{ return m_name; }

◆ phiCorr()

double Converter::phiCorr ( double phi)
inlineprotectedinherited

Definition at line 248 of file Converter.h.

248 {
249
250 if (phi < -TMath::Pi()) phi += 2*TMath::Pi();
251 if (phi > TMath::Pi()) phi -= 2*TMath::Pi();
252
253 return phi;
254 }
Scalar phi() const
phi method

◆ repSvcType()

virtual long TileRawChannelContByteStreamCnv::repSvcType ( ) const
inlineoverridevirtual

Storage type and class ID.

Definition at line 53 of file TileRawChannelContByteStreamCnv.h.

53{ return i_repSvcType(); }

◆ selectTracks() [1/3]

void Converter::selectTracks ( const Analysis::MuonContainer * muontracks)
inlineinherited

Definition at line 177 of file Converter.h.

177 {
178
179 static int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
180
181 Analysis::MuonContainer::const_iterator muonitr = muontracks->begin();
182 Analysis::MuonContainer::const_iterator muonend = muontracks->end();
183
184 while(muonitr!=muonend) {
185
186 bool hasInDet = (*muonitr)->hasInDetTrackParticle();
187 const Rec::TrackParticle *trackitr=NULL;
188 if(hasInDet) trackitr=(*muonitr)->inDetTrackParticle();
189 muonitr++;
190 if(!hasInDet) continue;
191
192 const Trk::MeasuredPerigee* measPer = trackitr->measuredPerigee();
193
194 CLHEP::HepVector perigeeParams = measPer->parameters();
195 double pT = measPer->pT();
196 double eta = measPer->eta();
197 double phi = perigeeParams[Trk::phi0];
198 double z0 = perigeeParams[Trk::z0];
199 double d0 = perigeeParams[Trk::d0];
200 //double theta = perigeeParams[Trk::theta];
201
202 double deta = 1;
203 double dphi = 1;
204 double dz0 = 1;
205 double dd0 = 1;
206 double dpT = 1;
207
208 // Check number of hits
209 // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really
210 // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
211 // stereo clusters making a spacepoint are two "hits"
212 const Trk::TrackSummary *summary = trackitr->trackSummary();
213 int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits);
215 int nSctHits = summary->get(Trk::numberOfSCTHits);
216 int nStrawHits = summary->get(Trk::numberOfTRTHits);
217 int nTrHits = summary->get(Trk::numberOfTRTHighThresholdHits);
218
219 int nSiHits = nPixelHits + nSctHits;
220
221 const Trk::FitQuality *quality = trackitr->fitQuality();
222 double chi2 = quality->chiSquared();
223 double dof = quality->quality->numberDoF();
224
225 unsigned bitmap = 0;
226
227 for ( int ih=0 ; ih<20 ; ih++ ) {
228 if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );
229 }
230
231 // Shift d0 and z0 according to beam position
232 ipCorr(d0, z0, d0, z0, phi, eta, pT);
233
234 // Create and save Track
235 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
236 deta, dphi, dz0, dd0, dpT,
237 nBlayerHits, nPixelHits, nSctHits, nSiHits,
238 nStrawHits, nTrHits, bitmap, 0,
239 -1) ;
240
241 addTrack(t);
242 }
243
244 }
void addTrack(TIDA::Track *t)
Definition Converter.h:45
void ipCorr(double d0, double z0, double &d0c, double &z0c, double phi0, double eta, double pt)
Definition Converter.h:257
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:838
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
const Trk::Perigee * measuredPerigee() const
Accessor method for Perigee.
int numberDoF() const
returns the number of degrees of freedom of the overall track or vertex fit as integer
Definition FitQuality.h:60
double chiSquared() const
returns the of the overall track fit
Definition FitQuality.h:56
double eta() const
Access method for pseudorapidity - from momentum.
double pT() const
Access method for transverse momentum.
const FitQuality * fitQuality() const
accessor function for FitQuality.
const TrackSummary * trackSummary() const
accessor function for TrackSummary.
double chi2(TH1 *h0, TH1 *h1)
float nSiHits(const U &p)
float nPixelHits(const U &p)
@ phi0
Definition ParamDefs.h:65
@ d0
Definition ParamDefs.h:63
@ z0
Definition ParamDefs.h:64
DetectorType
enumerates the various detector types currently accessible from the isHit() method.
@ numberOfPixelHits
number of pixel layers on track with absence of hits
@ numberOfBLayerHits
these are the hits in the 0th pixel layer?
@ numberOfTRTHighThresholdHits
total number of TRT hits which pass the high threshold
Perigee MeasuredPerigee

◆ selectTracks() [2/3]

void Converter::selectTracks ( const Rec::TrackParticleContainer * trigtracks)
inlineinherited

Definition at line 111 of file Converter.h.

111 {
112
113 Rec::TrackParticleContainer::const_iterator trackitr = trigtracks->begin();
114 Rec::TrackParticleContainer::const_iterator trackend = trigtracks->end();
115
116 static int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
117
118 while ( trackitr!=trackend ) {
119
120 const Trk::MeasuredPerigee* measPer = (*trackitr)->measuredPerigee();
121
122 CLHEP::HepVector perigeeParams = measPer->parameters();
123 double pT = measPer->pT();
124 double eta = measPer->eta();
125 double phi = perigeeParams[Trk::phi0];
126 double z0 = perigeeParams[Trk::z0];
127 double d0 = perigeeParams[Trk::d0];
128 //double theta = perigeeParams[Trk::theta];
129
130 double deta = 1;
131 double dphi = 1;
132 double dz0 = 1;
133 double dd0 = 1;
134 double dpT = 1;
135
136 // Check number of hits
137 // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really
138 // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
139 // stereo clusters making a spacepoint are two "hits"
140 const Trk::TrackSummary *summary = (*trackitr)->trackSummary();
141 int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits);
143 int nSctHits = summary->get(Trk::numberOfSCTHits);
144 int nStrawHits = summary->get(Trk::numberOfTRTHits);
145 int nTrHits = summary->get(Trk::numberOfTRTHighThresholdHits);
146 int nSiHits = nPixelHits + nSctHits;
147 bool expectBL = false; // Not stored for Rec::TrackParticle
148
149 const Trk::FitQuality *quality = (*trackitr)->fitQuality();
150 double chi2 = quality->chiSquared();
151 double dof = quality->numberDoF();
152
153 unsigned bitmap = 0;
154
155 for ( int ih=0 ; ih<20 ; ih++ ) {
156 if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );
157 }
158
159 // Shift d0 and z0 according to beam position
160 ipCorr(d0, z0, d0, z0, phi, eta, pT);
161
162 // Create and save Track
163 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
164 deta, dphi, dz0, dd0, dpT,
165 nBlayerHits, nPixelHits, nSctHits, nSiHits,
166 nStrawHits, nTrHits, bitmap, 0,
167 -1,
168 expectBL) ;
169
170 addTrack(t);
171 trackitr++;
172 }
173
174 }

◆ selectTracks() [3/3]

void Converter::selectTracks ( const TrigInDetTrackCollection * trigtracks)
inlineinherited

Definition at line 56 of file Converter.h.

56 {
57
58 TrigInDetTrackCollection::const_iterator trackitr = trigtracks->begin();
59 TrigInDetTrackCollection::const_iterator trackend = trigtracks->end();
60
61 while ( trackitr!=trackend ) {
62
63 double eta = (*trackitr)->param()->eta();
64 double phi = (*trackitr)->param()->phi0();
65 double z0 = (*trackitr)->param()->z0();
66 double pT = (*trackitr)->param()->pT();
67 double d0 = (*trackitr)->param()->a0();
68 //double theta = 2*atan2(exp(-(*trackitr)->param()->eta()),1);
69
70 double deta = (*trackitr)->param()->eeta();
71 double dphi = (*trackitr)->param()->ephi0();
72 double dz0 = (*trackitr)->param()->ez0();
73 double dpT = (*trackitr)->param()->epT();
74 double dd0 = (*trackitr)->param()->ea0();
75
76 int algoid = (*trackitr)->algorithmId();
77
78 int nBlayerHits = ((*trackitr)->HitPattern() & 0x1);
79 int nPixelHits = 2*(*trackitr)->NPixelSpacePoints(); // NB: for comparison with offline
80 int nSctHits = 2*(*trackitr)->NSCT_SpacePoints(); // a spacepoint is 2 "hits"
81 int nStrawHits = (*trackitr)->NStrawHits();
82 int nTrHits = (*trackitr)->NTRHits();
83
84 int nSiHits = nPixelHits + nSctHits;
85 bool expectBL = false; //not available with TrigInDetTrack
86
87 unsigned hitPattern = (*trackitr)->HitPattern();
88 unsigned multiPattern = 0;
89
90 double chi2 = (*trackitr)->chi2();
91 double dof = 0;
92
93 // Shift d0 and z0 according to beam position
94 ipCorr(d0, z0, d0, z0, phi, eta, pT);
95
96 // Create and save Track
97 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
98 deta, dphi, dz0, dd0, dpT,
99 nBlayerHits, nPixelHits, nSctHits, nSiHits,
100 nStrawHits, nTrHits, hitPattern, multiPattern,
101 algoid,
102 expectBL ) ;
103
104 addTrack(t);
105 trackitr++;
106 }
107 }

◆ setLevel()

void AthMessaging::setLevel ( MSG::Level lvl)
inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

29{
30 m_lvl = lvl;
31}

◆ storageType()

long TileRawChannelContByteStreamCnv::storageType ( )
static

Definition at line 48 of file TileRawChannelContByteStreamCnv.cxx.

static constexpr long storageType()

◆ tracks()

const std::vector< TIDA::Track * > & Converter::tracks ( ) const
inlineinherited

Definition at line 50 of file Converter.h.

50{return m_tracks;}

Member Data Documentation

◆ ATLAS_THREAD_SAFE

Athena::RecyclableDataQueue<TileMutableRawChannelContainer> m_queue TileRawChannelContByteStreamCnv::ATLAS_THREAD_SAFE
mutableprivate

Queue of data objects to recycle.

Definition at line 74 of file TileRawChannelContByteStreamCnv.h.

◆ m_beamX

double Converter::m_beamX
protectedinherited

Definition at line 320 of file Converter.h.

◆ m_beamY

double Converter::m_beamY
protectedinherited

Definition at line 320 of file Converter.h.

◆ m_beamZ

double Converter::m_beamZ
protectedinherited

Definition at line 320 of file Converter.h.

◆ m_byteStreamCnvSvc

ServiceHandle<IByteStreamCnvSvc> TileRawChannelContByteStreamCnv::m_byteStreamCnvSvc
private

Definition at line 61 of file TileRawChannelContByteStreamCnv.h.

◆ m_decoder

ToolHandle<TileROD_Decoder> TileRawChannelContByteStreamCnv::m_decoder
private

Pointer to TileROD_Decoder.

Definition at line 67 of file TileRawChannelContByteStreamCnv.h.

◆ m_hid2re

const TileHid2RESrcID* TileRawChannelContByteStreamCnv::m_hid2re
private

Pointer to TileHid2RESrcID.

Definition at line 70 of file TileRawChannelContByteStreamCnv.h.

◆ m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }
mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

135{ nullptr };

◆ m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }
mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

138{ MSG::NIL };

◆ m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls
mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

◆ m_name

std::string AthConstConverter::m_name
privateinherited

Definition at line 71 of file AthConstConverter.h.

◆ m_nm

std::string AthMessaging::m_nm
privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

◆ m_robSvc

ServiceHandle<IROBDataProviderSvc> TileRawChannelContByteStreamCnv::m_robSvc
private

Pointer to IROBDataProviderSvc.

Definition at line 64 of file TileRawChannelContByteStreamCnv.h.

◆ m_tool

ToolHandle<TileRawChannelContByteStreamTool> TileRawChannelContByteStreamCnv::m_tool
private

Definition at line 59 of file TileRawChannelContByteStreamCnv.h.

◆ m_tracks

std::vector<TIDA::Track*> Converter::m_tracks
protectedinherited

Definition at line 321 of file Converter.h.


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