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

ByteStream converter for L1TriggerResult. More...

#include <L1TriggerResultByteStreamCnv.h>

Inheritance diagram for L1TriggerResultByteStreamCnv:
Collaboration diagram for L1TriggerResultByteStreamCnv:

Public Member Functions

 L1TriggerResultByteStreamCnv (ISvcLocator *svcLoc)
 Standard constructor.
virtual ~L1TriggerResultByteStreamCnv ()
 Standard destructor.
virtual StatusCode initialize () override
virtual StatusCode finalize () override
virtual StatusCode createObj (IOpaqueAddress *pAddr, DataObject *&pObj) override
 Create xAOD (L1TriggerResult) from ByteStream.
virtual StatusCode createRep (DataObject *pObj, IOpaqueAddress *&pAddr) override
 Create ByteStream from xAOD (L1TriggerResult).
long repSvcType () const override
 return repSvcType
void addTrack (TIDA::Track *t)
const std::vector< TIDA::Track * > & tracks () const
void clear ()
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 ()
 Storage type used by this converter.
static CLID classID ()
 CLID of the class of the L1TriggerResult converted by this converter (xAOD::TrigCompositeContainer).

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 printRob (const OFFLINE_FRAGMENTS_NAMESPACE_WRITE::ROBFragment &rob) const
 Helper method for debug printouts.
void initMessaging () const
 Initialize our message level and MessageSvc.

Private Attributes

ServiceHandle< IByteStreamEventAccessm_ByteStreamEventAccess {"ByteStreamCnvSvc", "L1TriggerResultByteStreamCnv"}
 Helper to obtain the RawEvent pointer.
ToolHandle< IL1TriggerByteStreamToolm_muonEncoderTool {"MuonRoIByteStreamTool/L1MuonBSEncoderTool"}
 Encoder tools for L1Muon RoIs (one writing RoIB ROB, one writing DAQ ROB).
ToolHandle< IL1TriggerByteStreamToolm_muonEncoderToolDaq {"MuonRoIByteStreamTool/L1MuonBSEncoderToolDAQ"}
ToolHandle< IL1TriggerByteStreamToolm_ctpResultEncoderTool {"CTPResultByteStreamTool/CTPResultBSEncoderTool"}
 Encoder tool for CTP result.
ToolHandle< IL1TriggerByteStreamToolm_efexEncoderTool {"eFexByteStreamTool/eFexBSEncoderTool"}
 eFEX tool
ToolHandle< IL1TriggerByteStreamToolm_jfexEncoderTool {"jFexRoiByteStreamTool/jFexBSEncoderTool"}
 jFEX tool
ToolHandle< IL1TriggerByteStreamToolm_gfexEncoderTool {"gFexByteStreamTool/gFexBSEncoderTool"}
 gFEX tool
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.
std::atomic_flag m_initialized ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
 Messaging initialized (initMessaging).

Detailed Description

ByteStream converter for L1TriggerResult.

Definition at line 27 of file L1TriggerResultByteStreamCnv.h.

Constructor & Destructor Documentation

◆ L1TriggerResultByteStreamCnv()

L1TriggerResultByteStreamCnv::L1TriggerResultByteStreamCnv ( ISvcLocator * svcLoc)

Standard constructor.

Definition at line 29 of file L1TriggerResultByteStreamCnv.cxx.

29 :
30 Converter(storageType(), classID(), svcLoc),
31 AthMessaging(msgSvc(), "L1TriggerResultByteStreamCnv") {}
AthMessaging()
Default constructor:
static CLID classID()
CLID of the class of the L1TriggerResult converted by this converter (xAOD::TrigCompositeContainer).
static long storageType()
Storage type used by this converter.

◆ ~L1TriggerResultByteStreamCnv()

L1TriggerResultByteStreamCnv::~L1TriggerResultByteStreamCnv ( )
virtual

Standard destructor.

Definition at line 36 of file L1TriggerResultByteStreamCnv.cxx.

36{}

Member Function Documentation

◆ addTrack()

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

Definition at line 43 of file Converter.h.

43 {
44 m_tracks.push_back(t);
45 }
std::vector< TIDA::Track * > m_tracks
Definition Converter.h:264

◆ classID()

CLID L1TriggerResultByteStreamCnv::classID ( )
static

CLID of the class of the L1TriggerResult converted by this converter (xAOD::TrigCompositeContainer).

Definition at line 244 of file L1TriggerResultByteStreamCnv.cxx.

◆ clear()

void Converter::clear ( )
inlineinherited

Definition at line 51 of file Converter.h.

51{m_tracks.clear();}

◆ createObj()

StatusCode L1TriggerResultByteStreamCnv::createObj ( IOpaqueAddress * pAddr,
DataObject *& pObj )
overridevirtual

Create xAOD (L1TriggerResult) from ByteStream.

Definition at line 100 of file L1TriggerResultByteStreamCnv.cxx.

100 {
101 ATH_MSG_ERROR("L1TriggerResult cannot be created directly from ByteStream!"
102 << " Use the L1TriggerResultMaker algorithm instead");
103 return StatusCode::FAILURE;
104}
#define ATH_MSG_ERROR(x)

◆ createRep()

StatusCode L1TriggerResultByteStreamCnv::createRep ( DataObject * pObj,
IOpaqueAddress *& pAddr )
overridevirtual

Create ByteStream from xAOD (L1TriggerResult).

Definition at line 109 of file L1TriggerResultByteStreamCnv.cxx.

109 {
110 ATH_MSG_VERBOSE("start of " << __FUNCTION__);
111
112 const EventContext& ctx = Gaudi::Hive::currentContext();
113
114 // Cast the DataObject to L1TriggerResult
115 xAOD::TrigCompositeContainer* l1TriggerResult = nullptr;
116 bool castSuccessful = SG::fromStorable(pObj, l1TriggerResult);
117 if (!castSuccessful || !l1TriggerResult) {
118 ATH_MSG_ERROR("Failed to convert DataObject to xAOD::TrigCompositeContainer for L1TriggerResult");
119 return StatusCode::FAILURE;
120 }
121
122 // Obtain the RawEventWrite (aka eformat::write::FullEventFragment) pointer
123 RawEventWrite* re = m_ByteStreamEventAccess->getRawEvent();
124 if (!re) {
125 ATH_MSG_ERROR("Failed to obtain a pointer to RawEventWrite");
126 return StatusCode::FAILURE;
127 }
128 ATH_MSG_VERBOSE("Obtained RawEventWrite pointer = " << re);
129
130 // ===== CTPResult encoding =================
131
132 // Only perform encoding if tool was enabled, hence only when xAOD::CTPResult is used instead of ROIB::CTPResult (as part of ROIB::RoIBResult)
133 if (m_ctpResultEncoderTool.isEnabled()) {
134
135 // Extract the CTPResult SG key stored as a detail of the TrigComposite
136 const xAOD::TrigComposite* l1tr = l1TriggerResult->at(0);
137 std::string ctpKey;
138 l1tr->getDetail<std::string>("CTPResultKey", ctpKey);
139
140 // Use the retrieved key to read the CTPResult from the event store
141 SG::ReadHandle<xAOD::CTPResult> result(ctpKey, ctx);
142 ATH_CHECK(result.isValid());
143
144 // Helpful lambda function to convert from vectors of 32-bit words to bitsets
145 auto wordsToBitset = [](const std::vector<uint32_t>& words) {
146 std::bitset<512> out;
147 for (size_t i = 0; i < words.size(); ++i) {
148 uint32_t w = words[i];
149 for (size_t b = 0; b < 32; ++b) {
150 if (w & (1u << b)) {
151 out.set(i * 32 + b);
152 }
153 }
154 }
155 return out;
156 };
157
158 std::bitset<512> tbp = wordsToBitset(result->getTBPWords());
159 std::bitset<512> tap = wordsToBitset(result->getTAPWords());
160 std::bitset<512> tav = wordsToBitset(result->getTAVWords());
161
162 // Number of L1 words computed from TBP/TAP/TAV
163 constexpr size_t wordSize = 32;
164 constexpr size_t wordsPerSet = 512 / wordSize;
165 constexpr size_t numWords = 3 * wordsPerSet;
166
167 // Allocate per-event storage
168 std::vector<uint32_t> l1BitsData(numWords, 0);
169
170 // Fill l1BitsData from TBP/TAP/TAV
171 size_t iset{0};
172 for (const std::bitset<512>& bset : {tbp, tap, tav}) {
173 const std::string sbits = bset.to_string();
174 const size_t iWordOutputStart = wordsPerSet * iset;
175 for (size_t iWordInSet = 0; iWordInSet < wordsPerSet; ++iWordInSet) {
176 const size_t bwordPos = 512 - (iWordInSet + 1) * wordSize;
177 std::bitset<wordSize> bword{sbits.substr(bwordPos, wordSize)};
178 l1BitsData[iWordOutputStart + iWordInSet] = static_cast<uint32_t>(bword.to_ulong());
179 }
180 ++iset;
181 }
182 const uint32_t triggerType = result->triggerType();
183
184 // Update RawEventWrite
185 re->lvl1_trigger_info(l1BitsData.size(), l1BitsData.data());
186 re->lvl1_trigger_type(static_cast<uint8_t>(triggerType & 0xFF));
187
188 // Check there is no CTPResult ReadHandleKey in the L1TriggerResult when CTP bytestream encoding tool is not enabled
189 } else {
190 const xAOD::TrigComposite* l1tr = l1TriggerResult->at(0);
191 std::string ctpKey;
192 if (l1tr->getDetail<std::string>("CTPResultKey", ctpKey)) {
193 ATH_MSG_WARNING("L1TriggerResult contains CTPResult ReadHandleKey but CTP bytestream encoding tool is not enabled");
194 }
195 }
196
197 // ===== MuonRoI + eFex encoding =================
198
199 for (ToolHandle<IL1TriggerByteStreamTool>& tool : {std::reference_wrapper(m_muonEncoderTool), std::reference_wrapper(m_muonEncoderToolDaq), std::reference_wrapper(m_efexEncoderTool), std::reference_wrapper(m_jfexEncoderTool), std::reference_wrapper(m_gfexEncoderTool), std::reference_wrapper(m_ctpResultEncoderTool)}) {
200 if (not tool.isEnabled()) {continue;}
201 std::vector<WROBF*> muon_robs;
202 ATH_CHECK(tool->convertToBS(muon_robs, l1TriggerResult, ctx)); // TODO: find a way to avoid ThreadLocalContext
203 ATH_MSG_DEBUG(tool.name() << " created " << muon_robs.size() << " L1Muon ROB Fragments");
204 for (WROBF* rob : muon_robs) {
205 printRob(*rob);
206 // Set LVL1 Trigger Type from the full event
207 rob->rod_lvl1_type(re->lvl1_trigger_type());
208 // Set LVL1 ID from the full event
209 rob->rod_lvl1_id(re->lvl1_id());
210 // Add the ROBFragment to the full event
211 re->append(rob);
212 ATH_MSG_DEBUG("Added ROB fragment 0x" << MSG::hex << rob->source_id() << MSG::dec << " to the output raw event");
213 }
214 }
215
216 // Placeholder for other systems: L1Topo, gFEX
217
218 // Create a ByteStreamAddress for L1TriggerResult
219 if ( pAddr != nullptr ) pAddr->release();
220 ByteStreamAddress* bsAddr = new ByteStreamAddress(classID(), pObj->registry()->name(), "");
221 pAddr = static_cast<IOpaqueAddress*>(bsAddr);
222 pAddr->addRef();
223
224 ATH_MSG_VERBOSE("end of " << __FUNCTION__);
225 return StatusCode::SUCCESS;
226}
const std::regex re(r_e)
#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)
OFFLINE_FRAGMENTS_NAMESPACE_WRITE::ROBFragment WROBF
OFFLINE_FRAGMENTS_NAMESPACE_WRITE::FullEventFragment RawEventWrite
data type for writing raw event
Definition RawEvent.h:39
#define sbits(u, n)
Intrinsic functions which do not work correctly due to differences in byte ordering.
const T * at(size_type n) const
Access an element, as an rvalue.
void printRob(const OFFLINE_FRAGMENTS_NAMESPACE_WRITE::ROBFragment &rob) const
Helper method for debug printouts.
ToolHandle< IL1TriggerByteStreamTool > m_ctpResultEncoderTool
Encoder tool for CTP result.
ToolHandle< IL1TriggerByteStreamTool > m_efexEncoderTool
eFEX tool
ToolHandle< IL1TriggerByteStreamTool > m_muonEncoderTool
Encoder tools for L1Muon RoIs (one writing RoIB ROB, one writing DAQ ROB).
ToolHandle< IL1TriggerByteStreamTool > m_jfexEncoderTool
jFEX tool
ServiceHandle< IByteStreamEventAccess > m_ByteStreamEventAccess
Helper to obtain the RawEvent pointer.
ToolHandle< IL1TriggerByteStreamTool > m_gfexEncoderTool
gFEX tool
ToolHandle< IL1TriggerByteStreamTool > m_muonEncoderToolDaq
bool getDetail(const std::string &name, TYPE &value) const
Get an TYPE detail from the object.
bool fromStorable(DataObject *pDObj, T *&pTrans, bool quiet=false, IRegisterTransient *irt=0, bool isConst=true)
TrigCompositeContainer_v1 TrigCompositeContainer
Declare the latest version of the container.
TrigComposite_v1 TrigComposite
Declare the latest version of the class.
setBGCode tap
setEventNumber uint32_t

◆ finalize()

StatusCode L1TriggerResultByteStreamCnv::finalize ( )
overridevirtual

Definition at line 77 of file L1TriggerResultByteStreamCnv.cxx.

77 {
78 ATH_MSG_VERBOSE("start of " << __FUNCTION__);
79 if (m_ByteStreamEventAccess.release().isFailure())
80 ATH_MSG_WARNING("Failed to release service " << m_ByteStreamEventAccess.typeAndName());
81 if (m_muonEncoderTool.isEnabled() && m_muonEncoderTool.release().isFailure())
82 ATH_MSG_WARNING("Failed to release tool " << m_muonEncoderTool.typeAndName());
83 if (m_muonEncoderToolDaq.isEnabled() && m_muonEncoderToolDaq.release().isFailure())
84 ATH_MSG_WARNING("Failed to release tool " << m_muonEncoderToolDaq.typeAndName());
85 if (m_ctpResultEncoderTool.isEnabled() && m_ctpResultEncoderTool.release().isFailure())
86 ATH_MSG_WARNING("Failed to release tool " << m_ctpResultEncoderTool.typeAndName());
87 if (m_efexEncoderTool.isEnabled() && m_efexEncoderTool.release().isFailure())
88 ATH_MSG_WARNING("Failed to release tool " << m_efexEncoderTool.typeAndName());
89 if (m_jfexEncoderTool.isEnabled() && m_jfexEncoderTool.release().isFailure())
90 ATH_MSG_WARNING("Failed to release tool " << m_jfexEncoderTool.typeAndName());
91 if (m_gfexEncoderTool.isEnabled() && m_gfexEncoderTool.release().isFailure())
92 ATH_MSG_WARNING("Failed to release tool " << m_gfexEncoderTool.typeAndName());
93 ATH_MSG_VERBOSE("end of " << __FUNCTION__);
94 return StatusCode::SUCCESS;
95}

◆ initialize()

StatusCode L1TriggerResultByteStreamCnv::initialize ( )
overridevirtual

Definition at line 41 of file L1TriggerResultByteStreamCnv.cxx.

41 {
42 ATH_MSG_VERBOSE("start of " << __FUNCTION__);
44
45 // Check one property of the encoder tools to determine if the tools are configured in the job
46 const bool doMuon = not serviceLocator()->getOptsSvc().get("ToolSvc.L1MuonBSEncoderTool.ROBIDs").empty();
47 ATH_MSG_DEBUG("MUCTPI BS encoding is " << (doMuon ? "enabled" : "disabled"));
48 ATH_CHECK(m_muonEncoderTool.retrieve(EnableTool(doMuon)));
49
50 const bool doMuonDaq = not serviceLocator()->getOptsSvc().get("ToolSvc.L1MuonBSEncoderToolDAQ.ROBIDs").empty();
51 ATH_MSG_DEBUG("MUCTPI DAQ ROB encoding is " << (doMuonDaq ? "enabled" : "disabled"));
52 ATH_CHECK(m_muonEncoderToolDaq.retrieve(EnableTool(doMuonDaq)));
53
54 const bool doCTP = not serviceLocator()->getOptsSvc().get("ToolSvc.CTPResultBSEncoderTool.ROBIDs").empty();
55 ATH_MSG_DEBUG("CTP BS encoding is " << (doCTP ? "enabled" : "disabled"));
56 ATH_CHECK(m_ctpResultEncoderTool.retrieve(EnableTool(doCTP)));
57
58 const bool doEfex = not serviceLocator()->getOptsSvc().get("ToolSvc.eFexBSEncoderTool.ROBIDs").empty();
59 ATH_MSG_DEBUG("eFex BS encoding is " << (doEfex ? "enabled" : "disabled"));
60 ATH_CHECK(m_efexEncoderTool.retrieve(EnableTool(doEfex)));
61
62 const bool doJfex = not serviceLocator()->getOptsSvc().get("ToolSvc.jFexBSEncoderTool.ROBIDs").empty();
63 ATH_MSG_DEBUG("jFex BS encoding is " << (doJfex ? "enabled" : "disabled"));
64 ATH_CHECK(m_jfexEncoderTool.retrieve(EnableTool(doJfex)));
65
66 const bool doGfex = not serviceLocator()->getOptsSvc().get("ToolSvc.gFexBSEncoderTool.ROBIDs").empty();
67 ATH_MSG_DEBUG("gFex BS encoding is " << (doGfex ? "enabled" : "disabled"));
68 ATH_CHECK(m_gfexEncoderTool.retrieve(EnableTool(doGfex)));
69
70 ATH_MSG_VERBOSE("end of " << __FUNCTION__);
71 return StatusCode::SUCCESS;
72}

◆ 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)
int outputLevel(const IMessageSvc *ims, const std::string &source)

◆ ipCorr()

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

Definition at line 200 of file Converter.h.

200 {
201
202 double sn = sin(phi0);
203 double cs = cos(phi0);
204 double sd0 = (d0 != 0 ? d0/fabs(d0) : 1);
205 double spt = (pt != 0 ? pt/fabs(pt) : 1);
206
207 if (fabs(pt) >= 1*CLHEP::TeV) {
208
209 d0c = d0 + m_beamX*sn - m_beamY*cs;
210
211 } else {
212
213 double rc = fabs(pt)*15.0/(9.0*1.042);
214
215 double xc = (fabs(d0)-spt*sd0*rc)*cos(phi0+M_PI/2*sd0) - m_beamX;
216 double yc = (fabs(d0)-spt*sd0*rc)*sin(phi0+M_PI/2*sd0) - m_beamY;
217
218 double newphi;
219 double xd01,yd01,xd02,yd02;
220
221 if (xc == 0) {
222 xd01 = 0; yd01 = rc+yc;
223 xd02 = 0; yd02 = -rc+yc;
224 } else {
225 xd01 = xc+yc/xc*yc+sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd01 = yc/xc*xd01;
226 xd02 = xc+yc/xc*yc-sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd02 = yc/xc*xd02;
227 }
228
229 double r1 = sqrt(xd01*xd01+yd01*yd01);
230 double r2 = sqrt(xd02*xd02+yd02*yd02);
231
232 double phiV;
233
234 if (r1 < r2)
235 phiV = atan2(yd01,xd01);
236 else
237 phiV = atan2(yd02,xd02);
238
239 double phi1 = phiCorr(phiV+M_PI/2);
240 double phi2 = phiCorr(phiV-M_PI/2);
241
242 if (fabs(phiCorr(phi1-phi0))<=fabs(phiCorr(phi2-phi0)))
243 newphi = phi1;
244 else
245 newphi = phi2;
246
247 d0c = fabs(sqrt(xc*xc+yc*yc)-rc)*sin(phiV-newphi);
248
249 double theta=2*atan2(exp(-eta),1);
250 double theta_save=theta;
251 theta = theta - (1+spt)/2*M_PI;
252 if (theta>0) theta = theta_save;
253
254 double deltaz0= -spt*rc/tan(theta)*phiCorr(phiCorr(newphi)-phiCorr(phi0));
255 z0c = z0 + deltaz0;
256
257 }
258 }
#define M_PI
Scalar eta() const
pseudorapidity method
Scalar theta() const
theta method
static Double_t rc
double m_beamX
Definition Converter.h:263
double phiCorr(double phi)
Definition Converter.h:191
double m_beamY
Definition Converter.h:263
constexpr int pow(int x)
Definition conifer.h:27

◆ 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 167 of file AthMessaging.h.

168{
169 MsgStream* ms = m_msg_tls.get();
170 if (!ms) {
171 if (!m_initialized.test_and_set()) initMessaging();
172 ms = new MsgStream(m_imsg,m_nm);
173 m_msg_tls.reset( ms );
174 }
175
176 ms->setLevel (m_lvl);
177 return *ms;
178}
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 182 of file AthMessaging.h.

183{ 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 user did not set explicit message level we have to initialize
154 // the messaging and retrieve the default via the MessageSvc.
155 if (m_lvl==MSG::NIL && !m_initialized.test_and_set()) initMessaging();
156
157 if (m_lvl <= lvl) {
158 msg() << lvl;
159 return true;
160 } else {
161 return false;
162 }
163}

◆ phiCorr()

double Converter::phiCorr ( double phi)
inlineprotectedinherited

Definition at line 191 of file Converter.h.

191 {
192
193 if (phi < -TMath::Pi()) phi += 2*TMath::Pi();
194 if (phi > TMath::Pi()) phi -= 2*TMath::Pi();
195
196 return phi;
197 }
Scalar phi() const
phi method

◆ printRob()

void L1TriggerResultByteStreamCnv::printRob ( const OFFLINE_FRAGMENTS_NAMESPACE_WRITE::ROBFragment & rob) const
private

Helper method for debug printouts.

Definition at line 231 of file L1TriggerResultByteStreamCnv.cxx.

231 {
232 if (not msgLvl(MSG::DEBUG)) {return;}
233 const uint32_t ndata = rob.rod_ndata();
234 const uint32_t* data = rob.rod_data();
235 ATH_MSG_DEBUG("This ROB has " << ndata << " data words");
236 for (uint32_t i=0; i<ndata; ++i, ++data) {
237 ATH_MSG_DEBUG("--- 0x" << MSG::hex << *data << MSG::dec);
238 }
239}
bool msgLvl(const MSG::Level lvl) const
Test the output level.

◆ repSvcType()

long L1TriggerResultByteStreamCnv::repSvcType ( ) const
inlineoverride

return repSvcType

Definition at line 49 of file L1TriggerResultByteStreamCnv.h.

◆ selectTracks() [1/2]

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

Definition at line 120 of file Converter.h.

120 {
121
122 static int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
123
124 Analysis::MuonContainer::const_iterator muonitr = muontracks->begin();
125 Analysis::MuonContainer::const_iterator muonend = muontracks->end();
126
127 while(muonitr!=muonend) {
128
129 bool hasInDet = (*muonitr)->hasInDetTrackParticle();
130 const Rec::TrackParticle *trackitr=NULL;
131 if(hasInDet) trackitr=(*muonitr)->inDetTrackParticle();
132 muonitr++;
133 if(!hasInDet) continue;
134
135 const Trk::MeasuredPerigee* measPer = trackitr->measuredPerigee();
136
137 CLHEP::HepVector perigeeParams = measPer->parameters();
138 double pT = measPer->pT();
139 double eta = measPer->eta();
140 double phi = perigeeParams[Trk::phi0];
141 double z0 = perigeeParams[Trk::z0];
142 double d0 = perigeeParams[Trk::d0];
143 //double theta = perigeeParams[Trk::theta];
144
145 double deta = 1;
146 double dphi = 1;
147 double dz0 = 1;
148 double dd0 = 1;
149 double dpT = 1;
150
151 // Check number of hits
152 // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really
153 // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
154 // stereo clusters making a spacepoint are two "hits"
155 const Trk::TrackSummary *summary = trackitr->trackSummary();
156 int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits);
158 int nSctHits = summary->get(Trk::numberOfSCTHits);
159 int nStrawHits = summary->get(Trk::numberOfTRTHits);
160 int nTrHits = summary->get(Trk::numberOfTRTHighThresholdHits);
161
162 int nSiHits = nPixelHits + nSctHits;
163
164 const Trk::FitQuality *quality = trackitr->fitQuality();
165 double chi2 = quality->chiSquared();
166 double dof = quality->quality->numberDoF();
167
168 unsigned bitmap = 0;
169
170 for ( int ih=0 ; ih<20 ; ih++ ) {
171 if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );
172 }
173
174 // Shift d0 and z0 according to beam position
175 ipCorr(d0, z0, d0, z0, phi, eta, pT);
176
177 // Create and save Track
178 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
179 deta, dphi, dz0, dd0, dpT,
180 nBlayerHits, nPixelHits, nSctHits, nSiHits,
181 nStrawHits, nTrHits, bitmap, 0,
182 -1) ;
183
184 addTrack(t);
185 }
186
187 }
void addTrack(TIDA::Track *t)
Definition Converter.h:43
void ipCorr(double d0, double z0, double &d0c, double &z0c, double phi0, double eta, double pt)
Definition Converter.h:200
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)
Perigee MeasuredPerigee
@ 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

◆ selectTracks() [2/2]

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

Definition at line 54 of file Converter.h.

54 {
55
57 Rec::TrackParticleContainer::const_iterator trackend = trigtracks->end();
58
59 static int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
60
61 while ( trackitr!=trackend ) {
62
63 const Trk::MeasuredPerigee* measPer = (*trackitr)->measuredPerigee();
64
65 CLHEP::HepVector perigeeParams = measPer->parameters();
66 double pT = measPer->pT();
67 double eta = measPer->eta();
68 double phi = perigeeParams[Trk::phi0];
69 double z0 = perigeeParams[Trk::z0];
70 double d0 = perigeeParams[Trk::d0];
71 //double theta = perigeeParams[Trk::theta];
72
73 double deta = 1;
74 double dphi = 1;
75 double dz0 = 1;
76 double dd0 = 1;
77 double dpT = 1;
78
79 // Check number of hits
80 // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really
81 // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
82 // stereo clusters making a spacepoint are two "hits"
83 const Trk::TrackSummary *summary = (*trackitr)->trackSummary();
84 int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits);
86 int nSctHits = summary->get(Trk::numberOfSCTHits);
87 int nStrawHits = summary->get(Trk::numberOfTRTHits);
89 int nSiHits = nPixelHits + nSctHits;
90 bool expectBL = false; // Not stored for Rec::TrackParticle
91
92 const Trk::FitQuality *quality = (*trackitr)->fitQuality();
93 double chi2 = quality->chiSquared();
94 double dof = quality->numberDoF();
95
96 unsigned bitmap = 0;
97
98 for ( int ih=0 ; ih<20 ; ih++ ) {
99 if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );
100 }
101
102 // Shift d0 and z0 according to beam position
103 ipCorr(d0, z0, d0, z0, phi, eta, pT);
104
105 // Create and save Track
106 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
107 deta, dphi, dz0, dd0, dpT,
108 nBlayerHits, nPixelHits, nSctHits, nSiHits,
109 nStrawHits, nTrHits, bitmap, 0,
110 -1,
111 expectBL) ;
112
113 addTrack(t);
114 trackitr++;
115 }
116
117 }

◆ 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 L1TriggerResultByteStreamCnv::storageType ( )
static

Storage type used by this converter.

Definition at line 248 of file L1TriggerResultByteStreamCnv.cxx.

248 {
250}
static constexpr long storageType()

◆ tracks()

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

Definition at line 48 of file Converter.h.

48{return m_tracks;}

Member Data Documentation

◆ ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
mutableprivateinherited

Messaging initialized (initMessaging).

Definition at line 141 of file AthMessaging.h.

◆ m_beamX

double Converter::m_beamX
protectedinherited

Definition at line 263 of file Converter.h.

◆ m_beamY

double Converter::m_beamY
protectedinherited

Definition at line 263 of file Converter.h.

◆ m_beamZ

double Converter::m_beamZ
protectedinherited

Definition at line 263 of file Converter.h.

◆ m_ByteStreamEventAccess

ServiceHandle<IByteStreamEventAccess> L1TriggerResultByteStreamCnv::m_ByteStreamEventAccess {"ByteStreamCnvSvc", "L1TriggerResultByteStreamCnv"}
private

Helper to obtain the RawEvent pointer.

Definition at line 56 of file L1TriggerResultByteStreamCnv.h.

56{"ByteStreamCnvSvc", "L1TriggerResultByteStreamCnv"};

◆ m_ctpResultEncoderTool

ToolHandle<IL1TriggerByteStreamTool> L1TriggerResultByteStreamCnv::m_ctpResultEncoderTool {"CTPResultByteStreamTool/CTPResultBSEncoderTool"}
private

Encoder tool for CTP result.

Definition at line 63 of file L1TriggerResultByteStreamCnv.h.

63{"CTPResultByteStreamTool/CTPResultBSEncoderTool"};

◆ m_efexEncoderTool

ToolHandle<IL1TriggerByteStreamTool> L1TriggerResultByteStreamCnv::m_efexEncoderTool {"eFexByteStreamTool/eFexBSEncoderTool"}
private

eFEX tool

Definition at line 65 of file L1TriggerResultByteStreamCnv.h.

65{"eFexByteStreamTool/eFexBSEncoderTool"};

◆ m_gfexEncoderTool

ToolHandle<IL1TriggerByteStreamTool> L1TriggerResultByteStreamCnv::m_gfexEncoderTool {"gFexByteStreamTool/gFexBSEncoderTool"}
private

gFEX tool

Definition at line 69 of file L1TriggerResultByteStreamCnv.h.

69{"gFexByteStreamTool/gFexBSEncoderTool"};

◆ m_imsg

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

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

135{ nullptr };

◆ m_jfexEncoderTool

ToolHandle<IL1TriggerByteStreamTool> L1TriggerResultByteStreamCnv::m_jfexEncoderTool {"jFexRoiByteStreamTool/jFexBSEncoderTool"}
private

jFEX tool

Definition at line 67 of file L1TriggerResultByteStreamCnv.h.

67{"jFexRoiByteStreamTool/jFexBSEncoderTool"};

◆ 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_muonEncoderTool

ToolHandle<IL1TriggerByteStreamTool> L1TriggerResultByteStreamCnv::m_muonEncoderTool {"MuonRoIByteStreamTool/L1MuonBSEncoderTool"}
private

Encoder tools for L1Muon RoIs (one writing RoIB ROB, one writing DAQ ROB).

Definition at line 60 of file L1TriggerResultByteStreamCnv.h.

60{"MuonRoIByteStreamTool/L1MuonBSEncoderTool"};

◆ m_muonEncoderToolDaq

ToolHandle<IL1TriggerByteStreamTool> L1TriggerResultByteStreamCnv::m_muonEncoderToolDaq {"MuonRoIByteStreamTool/L1MuonBSEncoderToolDAQ"}
private

Definition at line 61 of file L1TriggerResultByteStreamCnv.h.

61{"MuonRoIByteStreamTool/L1MuonBSEncoderToolDAQ"};

◆ m_nm

std::string AthMessaging::m_nm
privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

◆ m_tracks

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

Definition at line 264 of file Converter.h.


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