13 #include "CLHEP/Geometry/Point3D.h"
14 #include "CLHEP/Units/SystemOfUnits.h"
17 #include "GaudiKernel/MsgStream.h"
47 static const double dRcut = 1.0e-7*
CLHEP::mm;
48 static const double dTcut = 1.0*
CLHEP::ns;
56 unsigned int endBC = 0;
57 unsigned int endId = 0;
58 unsigned int endHit = 0;
59 HepGeom::Point3D<double> lastEnd(0.0, 0.0, 0.0);
66 if ( trtHit->particleLink().barcode() != lastBarcode ||
idx - endBC > 65500) {
68 lastBarcode = trtHit->particleLink().barcode();
69 persCont->
m_barcode.push_back(lastBarcode);
71 persCont->
m_nBC.push_back(
idx - endBC);
76 if ( (
int)trtHit->GetParticleEncoding() != lastId ||
idx - endId > 65500) {
78 lastId = trtHit->GetParticleEncoding();
79 persCont->
m_id.push_back(lastId);
81 persCont->
m_nId.push_back(
idx - endId);
86 const HepGeom::Point3D<double> hitStart(trtHit->GetPreStepX(), trtHit->GetPreStepY(), trtHit->GetPreStepZ());
88 const double meanTime = trtHit->GetGlobalTime();
89 const double dTLast = fabs(meanTime - lastT);
90 const double dRLast = lastEnd.distance(hitStart);
95 if ( dRLast >= dRcut || dTLast >= dTcut ) {
107 const unsigned int strawId = trtHit->GetHitID();
108 persCont->
m_strawId1b.push_back( (
unsigned char)(strawId % 256) );
109 persCont->
m_strawId2b.push_back( (
unsigned short)(strawId / 256) );
110 if ( strawId>16777215 )
111 log << MSG::WARNING <<
"TRT_HitCollectionCnv: strawId > 2^24-1 cannot be persistified correctly! " <<
endmsg;
120 const double startR = sqrt( hitStart.x()*hitStart.x() + hitStart.y()*hitStart.y() );
121 unsigned short istartRflag;
134 const double startPhi = atan2( hitStart.y(), hitStart.x() );
135 persCont->
m_startPhi.push_back( (
unsigned char)( (startPhi+
M_PI)/(2.0*
M_PI)*256.0 ) );
151 istartZ = (istartZ << 1) | istartRflag;
152 persCont->
m_startZ.push_back( istartZ );
174 const HepGeom::Point3D<double> hitEnd(trtHit->GetPostStepX(), trtHit->GetPostStepY(), trtHit->GetPostStepZ());
175 const HepGeom::Point3D<double> hitLength = (hitEnd - hitStart);
194 double kinEne = trtHit->GetKineticEnergy() * 1.0e9;
195 double steplength = hitLength.distance() * 1.0e9;
196 if ( kinEne < 1.0 ) kinEne=1.0;
197 if ( steplength < 1.0 ) steplength=1.0;
198 if ( kinEne > 9.0e18 ) kinEne=9.0e18;
199 if ( steplength > 9.0e18 ) steplength=9.0e18;
200 const unsigned int kexponent = (
unsigned int)ceil(log10(kinEne)/0.30102999566398);
201 const unsigned int sexponent = (
unsigned int)ceil(log10(steplength)/0.30102999566398);
202 const unsigned int kmantissa = (
unsigned int)(kinEne/
pow(2.0,kexponent)*1024) - 512;
203 const unsigned int smantissa = (
unsigned int)(steplength/
pow(2.0,sexponent)*1024) - 512;
204 persCont->
m_kinEne.push_back( (kmantissa << 6) | kexponent );
205 persCont->
m_steplength.push_back( (smantissa << 6) | sexponent );
215 const double endR = sqrt( hitEnd.x()*hitEnd.x() + hitEnd.y()*hitEnd.y() );
216 unsigned short iendRflag;
222 persCont->
m_endR.push_back( (
unsigned char)(endR/(2.0*
CLHEP::mm)*256.0) );
230 const double endPhi = atan2( hitEnd.y(), hitEnd.x() );
231 persCont->
m_endPhi.push_back( (
unsigned char)( (endPhi+
M_PI)/(2.0*
M_PI)*256.0 ) );
239 unsigned short idZsign = (hitLength.z()>0.0) ? 1 : 0;
240 unsigned short imeanTime = ( meanTime < 75.0*
CLHEP::ns ) ? (
unsigned short)(meanTime/(75.0*
CLHEP::ns)*1024.0) : 1023;
241 if ( imeanTime == 1023 ) persCont->
m_meanTimeof.push_back( (
float)meanTime );
242 imeanTime = (imeanTime << 2) | (idZsign << 1) | iendRflag;
250 (
int)(abs(lastId)/100000) == 41 ||
251 (
int)(abs(lastId)/10000000) == 1
252 ) persCont->
m_hitEne.push_back( (
float)(trtHit->GetEnergyDeposit()) );
272 persCont->
m_nBC.push_back(
idx - endBC);
273 persCont->
m_nId.push_back(
idx - endId);
281 std::unique_ptr<TRTUncompressedHitCollection> trans(std::make_unique<TRTUncompressedHitCollection>(
"DefaultCollectionName",persObj->
m_nHits.size()));
283 return(trans.release());
294 unsigned int meanTimeofCount=0, startRCount=0, endRCount=0, hitEneCount=0;
295 unsigned int idxBC=0, idxId=0, endHit=0, endBC=0, endId=0;
305 for (
unsigned int i = 0;
i < persCont->
m_nHits.size();
i++ ) {
309 const unsigned int startHit = endHit;
317 const unsigned int strawId = i2*256+i1;
322 const unsigned int istartPhi = persCont->
m_startPhi[
i];
323 const double startPhi = -
M_PI + (istartPhi+0.5)*2.0*
M_PI/256.0;
328 const unsigned int istartZ = persCont->
m_startZ[
i] >> 1;
334 const unsigned int istartRflag = persCont->
m_startZ[
i] & 1;
340 if ( istartRflag == 1 ) {
344 const unsigned int istartR = persCont->
m_startR[startRCount++];
345 startR = (istartR+0.5)*2.0*
CLHEP::mm/256.0;
352 double startX = startR*
cos(startPhi);
353 double startY = startR*
sin(startPhi);
369 for (
unsigned int j = startHit; j < endHit; j++ ) {
371 if ( j >= endBC + persCont->
m_nBC[idxBC] ) endBC += persCont->
m_nBC[idxBC++];
372 if ( j >= endId + persCont->
m_nId[idxId] ) endId += persCont->
m_nId[idxId++];
377 const unsigned int imeanTime = persCont->
m_meanTime[j] >> 2;
378 double meanTime = (imeanTime+0.5)*75.0*
CLHEP::ns/1024.0;
379 if ( imeanTime == 1023 ) meanTime = (
double)persCont->
m_meanTimeof[meanTimeofCount++];
384 const unsigned int idZsign = (persCont->
m_meanTime[j] >> 1 ) & 1;
389 const unsigned int iendRflag = persCont->
m_meanTime[j] & 1;
394 const double hitEne = ( persCont->
m_id[idxId] == 22 ||
395 (
int)(abs(persCont->
m_id[idxId])/100000) == 41 ||
396 (
int)(abs(persCont->
m_id[idxId])/10000000) == 1
402 const unsigned int iendPhi = persCont->
m_endPhi[j];
403 double endPhi = -
M_PI + (iendPhi+0.5)*2.0*
M_PI/256.0;
409 if ( iendRflag==1 ) {
413 const unsigned int iendR = persCont->
m_endR[endRCount++];
421 double endX = endR*
cos(endPhi);
422 double endY = endR*
sin(endPhi);
434 const int kmantissa = persCont->
m_kinEne[j] >> 6;
436 const int kexponent = persCont->
m_kinEne[j] & 0x3F;
438 const double kinEne = (kmantissa+512.5)/1024 *
pow(2.0,kexponent) / 1.0e9;
439 double g4steplength = (smantissa+512.5)/1024 *
pow(2.0,sexponent) / 1.0e9;
440 if ( idZsign==0 ) g4steplength = -g4steplength;
445 double dX = endX-startX;
446 double dY = endY-startY;
448 double dXY2 = dX*dX+dY*dY;
449 double dL2 = g4steplength*g4steplength;
452 if (g4steplength<0.0) dZ=-dZ;
455 dX = dX * sqrt(dL2/dXY2);
456 dY = dY * sqrt(dL2/dXY2);
463 double endZ = startZ + dZ;
486 transCont->
Emplace( strawId, partLink, persCont->
m_id[idxId],
487 kinEne, hitEne, startX, startY, startZ,
488 endX, endY, endZ, meanTime );
494 startX = endXo; startY = endYo; startZ = endZ;