|
ATLAS Offline Software
|
#include <SiSpacePointsSeedMaker_Trigger.h>
|
|
| SiSpacePointsSeedMaker_Trigger (const std::string &, const std::string &, const IInterface *) |
|
virtual | ~SiSpacePointsSeedMaker_Trigger ()=default |
|
virtual StatusCode | initialize () override |
|
virtual StatusCode | finalize () override |
|
|
virtual void | newEvent (const EventContext &ctx, EventData &data, int iteration) const override |
|
virtual void | newRegion (const EventContext &ctx, EventData &data, const std::vector< IdentifierHash > &vPixel, const std::vector< IdentifierHash > &vSCT) const override |
|
virtual void | newRegion (const EventContext &ctx, EventData &data, const std::vector< IdentifierHash > &vPixel, const std::vector< IdentifierHash > &vSCT, const IRoiDescriptor &iRD) const override |
|
|
virtual void | find2Sp (EventData &data, const std::list< Trk::Vertex > &lv) const override |
| with two space points with or without vertex constraint More...
|
|
virtual void | find3Sp (const EventContext &ctx, EventData &data, const std::list< Trk::Vertex > &lv) const override |
| with three space points with or without vertex constraint More...
|
|
virtual void | find3Sp (const EventContext &ctx, EventData &data, const std::list< Trk::Vertex > &lv, const double *zVertex) const override |
| with three space points with or without vertex constraint with information about min and max Z of the vertex More...
|
|
virtual void | findVSp (const EventContext &ctx, EventData &data, const std::list< Trk::Vertex > &lv) const override |
| with variable number space points with or without vertex constraint Variable means (2,3,4,....) any number space points More...
|
|
|
produced accordingly methods find
|
virtual const SiSpacePointsSeed * | next (const EventContext &ctx, EventData &data) const override |
|
virtual void | writeNtuple (const SiSpacePointsSeed *seed, const Trk::Track *track, int seedType, long eventNumber) const override |
|
virtual bool | getWriteNtupleBoolProperty () const override |
|
|
|
SG::ReadHandleKey< SpacePointContainer > | m_spacepointsSCT {this, "SpacePointsSCTName", "SCT_SpacePoints", "SCT space points container"} |
|
SG::ReadHandleKey< SpacePointContainer > | m_spacepointsPixel {this, "SpacePointsPixelName", "PixelSpacePoints", "Pixel space points container"} |
|
SG::ReadHandleKey< SpacePointOverlapCollection > | m_spacepointsOverlap {this, "SpacePointsOverlapName", "OverlapSpacePoints"} |
|
SG::ReadCondHandleKey< InDet::BeamSpotData > | m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"} |
|
SG::ReadCondHandleKey< AtlasFieldCacheCondObj > | m_fieldCondObjInputKey |
|
|
BooleanProperty | m_pixel {this, "usePixel", true} |
|
BooleanProperty | m_sct {this, "useSCT", true} |
|
BooleanProperty | m_useOverlap {this, "useOverlapSpCollection", true} |
|
IntegerProperty | m_maxsize {this, "maxSize", 20000} |
|
IntegerProperty | m_maxsizeSP {this, "maxSizeSP", 1500} |
|
IntegerProperty | m_maxOneSize {this, "maxSeedsForSpacePoint", 3} |
|
UnsignedIntegerProperty | m_maxNumberVertices {this, "maxNumberVertices", 99} |
|
FloatProperty | m_r1min {this, "minRadius1", 0.} |
|
FloatProperty | m_r1minv {this, "minVRadius1", 0.} |
|
FloatProperty | m_r1max {this, "maxRadius1", 600.} |
|
FloatProperty | m_r1maxv {this, "maxVRadius1", 60.} |
|
FloatProperty | m_r2min {this, "minRadius2", 0.} |
|
FloatProperty | m_r2minv {this, "minVRadius2", 70.} |
|
FloatProperty | m_r2max {this, "maxRadius2", 600.} |
|
FloatProperty | m_r2maxv {this, "maxVRadius2", 200.} |
|
FloatProperty | m_r3min {this, "minRadius3", 0.} |
|
FloatProperty | m_drmin {this, "mindRadius", 10.} |
|
FloatProperty | m_drmax {this, "maxdRadius", 270.} |
|
FloatProperty | m_zmin {this, "minZ", -250.} |
|
FloatProperty | m_zmax {this, "maxZ", +250.} |
|
FloatProperty | m_r_rmax {this, "radMax", 600.} |
|
FloatProperty | m_r_rstep {this, "radStep", 2.} |
|
FloatProperty | m_dzver {this, "maxdZver", 5.} |
|
FloatProperty | m_dzdrver {this, "maxdZdRver", 0.02} |
|
FloatProperty | m_diver {this, "maxdImpact", 10.} |
|
FloatProperty | m_diverpps {this, "maxdImpactPPS", 1.7} |
|
FloatProperty | m_diversss {this, "maxdImpactSSS", 1000.} |
|
|
FloatProperty | m_etamax {this, "etaMax", 2.7} |
|
FloatProperty | m_ptmin {this, "pTmin", 500.} |
|
|
float | m_drminv {20.} |
|
|
bool | m_initialized {false} |
|
int | m_outputlevel {0} |
|
int | m_r_size {0} |
|
int | m_fNmax {0} |
|
int | m_fvNmax {0} |
|
int | m_rfz_b [SizeRFZ] {} |
|
int | m_rfz_t [SizeRFZ] {} |
|
int | m_rfz_ib [SizeRFZ][SizeI] {} |
|
int | m_rfz_it [SizeRFZ][SizeI] {} |
|
int | m_rfzv_n [SizeRFZV] {} |
|
int | m_rfzv_i [SizeRFZV][SizeIV] {} |
|
float | m_dzdrmin {0.} |
|
float | m_dzdrmax {0.} |
|
float | m_r3max {0.} |
|
float | m_ipt {0.} |
|
float | m_ipt2 {0.} |
|
float | m_COF {0.} |
|
float | m_sF {0.} |
|
float | m_sFv {0.} |
|
|
| SiSpacePointsSeedMaker_Trigger ()=delete |
|
| SiSpacePointsSeedMaker_Trigger (const SiSpacePointsSeedMaker_Trigger &)=delete |
|
SiSpacePointsSeedMaker_Trigger & | operator= (const SiSpacePointsSeedMaker_Trigger &)=delete |
|
MsgStream & | dumpConditions (EventData &data, MsgStream &out) const |
|
void | buildFrameWork () |
|
void | buildBeamFrameWork (EventData &data) const |
|
void | newOneSeed (EventData &data, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const float &, const float &) const |
|
void | fillLists (EventData &data) const |
|
void | production2Sp (EventData &data) const |
|
void | production3Sp (EventData &data) const |
|
void | production3Sp (EventData &data, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, int, int, int &) const |
|
void | production3SpTrigger (EventData &data, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, int, int, int &) const |
|
bool | newVertices (EventData &data, const std::list< Trk::Vertex > &) const |
|
void | findNext (EventData &data) const |
|
bool | isZCompatible (EventData &data, float &, float &, float &) const |
|
void | initializeEventData (EventData &data) const |
|
static MsgStream & | dumpEvent (EventData &data, MsgStream &out) |
|
static SiSpacePointForSeed * | newSpacePoint (EventData &data, const Trk::SpacePoint *const &) |
|
static void | newSeed (EventData &data, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const float &) |
|
static void | fillSeeds (EventData &data) |
|
static void | erase (EventData &data) |
|
static void | convertToBeamFrameWork (EventData &data, const Trk::SpacePoint *const &, float *) |
|
static float | dZVertexMin (EventData &data, float &) |
|
◆ Size
◆ SiSpacePointsSeedMaker_Trigger() [1/3]
InDet::SiSpacePointsSeedMaker_Trigger::SiSpacePointsSeedMaker_Trigger |
( |
const std::string & |
t, |
|
|
const std::string & |
n, |
|
|
const IInterface * |
p |
|
) |
| |
◆ ~SiSpacePointsSeedMaker_Trigger()
virtual InDet::SiSpacePointsSeedMaker_Trigger::~SiSpacePointsSeedMaker_Trigger |
( |
| ) |
|
|
virtualdefault |
◆ SiSpacePointsSeedMaker_Trigger() [2/3]
InDet::SiSpacePointsSeedMaker_Trigger::SiSpacePointsSeedMaker_Trigger |
( |
| ) |
|
|
privatedelete |
◆ SiSpacePointsSeedMaker_Trigger() [3/3]
◆ buildBeamFrameWork()
void InDet::SiSpacePointsSeedMaker_Trigger::buildBeamFrameWork |
( |
EventData & |
data | ) |
const |
|
private |
Definition at line 800 of file SiSpacePointsSeedMaker_Trigger.cxx.
805 double tx =
std::tan(beamSpotHandle->beamTilt(0));
806 double ty =
std::tan(beamSpotHandle->beamTilt(1));
808 double ph = std::atan2(ty,
tx);
809 double th = std::acos(1./std::sqrt(1.+
tx*
tx+ty*ty));
810 double sint =
sin(
th);
812 double sinp =
sin(ph);
813 double cosp =
cos(ph);
815 data.xbeam[0] =
static_cast<float>(cb.x());
816 data.xbeam[1] =
static_cast<float>(
cost*cosp*cosp+sinp*sinp);
817 data.xbeam[2] =
static_cast<float>(
cost*sinp*cosp-sinp*cosp);
818 data.xbeam[3] = -
static_cast<float>(sint*cosp );
820 data.ybeam[0] =
static_cast<float>(cb.y());
821 data.ybeam[1] =
static_cast<float>(
cost*cosp*sinp-sinp*cosp);
822 data.ybeam[2] =
static_cast<float>(
cost*sinp*sinp+cosp*cosp);
823 data.ybeam[3] = -
static_cast<float>(sint*sinp );
825 data.zbeam[0] =
static_cast<float>(cb.z());
826 data.zbeam[1] =
static_cast<float>(sint*cosp);
827 data.zbeam[2] =
static_cast<float>(sint*sinp);
828 data.zbeam[3] =
static_cast<float>(
cost);
◆ buildFrameWork()
void InDet::SiSpacePointsSeedMaker_Trigger::buildFrameWork |
( |
| ) |
|
|
private |
Definition at line 652 of file SiSpacePointsSeedMaker_Trigger.cxx.
672 const float sFmax =
static_cast<float>(NFmax)/
pi2;
673 const float sFmin = 100./60.;
677 else if (
m_sF < sFmin)
m_sF = sFmin;
683 const int NFtmax = 100;
684 const float sFvmax =
static_cast<float>(NFtmax)/
pi2;
◆ convertToBeamFrameWork()
◆ dump()
MsgStream & InDet::SiSpacePointsSeedMaker_Trigger::dump |
( |
EventData & |
data, |
|
|
MsgStream & |
out |
|
) |
| const |
|
overridevirtual |
◆ dumpConditions()
MsgStream & InDet::SiSpacePointsSeedMaker_Trigger::dumpConditions |
( |
EventData & |
data, |
|
|
MsgStream & |
out |
|
) |
| const |
|
private |
Definition at line 449 of file SiSpacePointsSeedMaker_Trigger.cxx.
452 std::string
s2;
for (
int i=0;
i<
n; ++
i)
s2.append(
" ");
s2.append(
"|");
454 std::string
s3;
for (
int i=0;
i<
n; ++
i)
s3.append(
" ");
s3.append(
"|");
456 std::string
s4;
for (
int i=0;
i<
n; ++
i)
s4.append(
" ");
s4.append(
"|");
458 std::string s5;
for (
int i=0;
i<
n; ++
i) s5.append(
" "); s5.append(
"|");
461 out<<
"|---------------------------------------------------------------------|"
475 <<std::setw(12)<<
m_sct
480 out<<
"| maxSizeSP | "
483 out<<
"| pTmin (mev) | "
484 <<std::setw(12)<<std::setprecision(5)<<
m_ptmin
487 <<std::setw(12)<<std::setprecision(5)<<
m_etamax
489 out<<
"| max radius SP | "
490 <<std::setw(12)<<std::setprecision(5)<<
m_r_rmax
492 out<<
"| radius step | "
493 <<std::setw(12)<<std::setprecision(5)<<
m_r_rstep
495 out<<
"| min Z-vertex position | "
496 <<std::setw(12)<<std::setprecision(5)<<
m_zmin
498 out<<
"| max Z-vertex position | "
499 <<std::setw(12)<<std::setprecision(5)<<
m_zmax
501 out<<
"| min radius first SP(3) | "
502 <<std::setw(12)<<std::setprecision(5)<<
m_r1min
504 out<<
"| min radius second SP(3) | "
505 <<std::setw(12)<<std::setprecision(5)<<
m_r2min
507 out<<
"| min radius last SP(3) | "
508 <<std::setw(12)<<std::setprecision(5)<<
m_r3min
510 out<<
"| max radius first SP(3) | "
511 <<std::setw(12)<<std::setprecision(4)<<
m_r1max
513 out<<
"| max radius second SP(3) | "
514 <<std::setw(12)<<std::setprecision(5)<<
m_r2max
516 out<<
"| max radius last SP(3) | "
517 <<std::setw(12)<<std::setprecision(5)<<
m_r3max
519 out<<
"| min radius first SP(2) | "
520 <<std::setw(12)<<std::setprecision(5)<<
m_r1minv
522 out<<
"| min radius second SP(2) | "
523 <<std::setw(12)<<std::setprecision(5)<<
m_r2minv
525 out<<
"| max radius first SP(2) | "
526 <<std::setw(12)<<std::setprecision(5)<<
m_r1maxv
528 out<<
"| max radius second SP(2) | "
529 <<std::setw(12)<<std::setprecision(5)<<
m_r2maxv
531 out<<
"| min space points dR | "
532 <<std::setw(12)<<std::setprecision(5)<<
m_drmin
534 out<<
"| max space points dR | "
535 <<std::setw(12)<<std::setprecision(5)<<
m_drmax
537 out<<
"| max dZ impact | "
538 <<std::setw(12)<<std::setprecision(5)<<
m_dzver
540 out<<
"| max dZ/dR impact | "
541 <<std::setw(12)<<std::setprecision(5)<<
m_dzdrver
543 out<<
"| max impact | "
544 <<std::setw(12)<<std::setprecision(5)<<
m_diver
546 out<<
"| max impact pps | "
547 <<std::setw(12)<<std::setprecision(5)<<
m_diverpps
549 out<<
"| max impact sss | "
550 <<std::setw(12)<<std::setprecision(5)<<
m_diversss
552 out<<
"|---------------------------------------------------------------------|"
554 out<<
"| Beam X center | "
555 <<std::setw(12)<<std::setprecision(5)<<
data.xbeam[0]
557 out<<
"| Beam Y center | "
558 <<std::setw(12)<<std::setprecision(5)<<
data.ybeam[0]
560 out<<
"| Beam Z center | "
561 <<std::setw(12)<<std::setprecision(5)<<
data.zbeam[0]
563 out<<
"| Beam X-axis direction | "
564 <<std::setw(12)<<std::setprecision(5)<<
data.xbeam[1]
565 <<std::setw(12)<<std::setprecision(5)<<
data.xbeam[2]
566 <<std::setw(12)<<std::setprecision(5)<<
data.xbeam[3]
568 out<<
"| Beam Y-axis direction | "
569 <<std::setw(12)<<std::setprecision(5)<<
data.ybeam[1]
570 <<std::setw(12)<<std::setprecision(5)<<
data.ybeam[2]
571 <<std::setw(12)<<std::setprecision(5)<<
data.ybeam[3]
573 out<<
"| Beam Z-axis direction | "
574 <<std::setw(12)<<std::setprecision(5)<<
data.zbeam[1]
575 <<std::setw(12)<<std::setprecision(5)<<
data.zbeam[2]
576 <<std::setw(12)<<std::setprecision(5)<<
data.zbeam[3]
578 out<<
"|---------------------------------------------------------------------|"
◆ dumpEvent()
MsgStream & InDet::SiSpacePointsSeedMaker_Trigger::dumpEvent |
( |
EventData & |
data, |
|
|
MsgStream & |
out |
|
) |
| |
|
staticprivate |
Definition at line 587 of file SiSpacePointsSeedMaker_Trigger.cxx.
589 out<<
"|---------------------------------------------------------------------|"
592 <<std::setw(12)<<
data.ns
594 out<<
"| data.nsaz | "
595 <<std::setw(12)<<
data.nsaz
597 out<<
"| data.nsazv | "
598 <<std::setw(12)<<
data.nsazv
601 <<std::setw(12)<<
data.l_seeds.size()
603 out<<
"|---------------------------------------------------------------------|"
◆ dZVertexMin()
float InDet::SiSpacePointsSeedMaker_Trigger::dZVertexMin |
( |
EventData & |
data, |
|
|
float & |
Z |
|
) |
| |
|
staticprivate |
◆ erase()
void InDet::SiSpacePointsSeedMaker_Trigger::erase |
( |
EventData & |
data | ) |
|
|
staticprivate |
◆ fillLists()
void InDet::SiSpacePointsSeedMaker_Trigger::fillLists |
( |
EventData & |
data | ) |
const |
|
private |
Definition at line 847 of file SiSpacePointsSeedMaker_Trigger.cxx.
854 if (!
data.r_map[
i])
continue;
855 r =
data.r_Sorted[
i].begin();
857 while (
r!=
data.r_Sorted[
i].end()) {
861 float F = (*r)->phi();
864 int f =
static_cast<int>(
F*
m_sF);
873 if (
Z>0.) {
Z < 250.?
z=5 :
z=6;}
874 else {
Z >-250.?
z=5 :
z=4;}
877 data.rfz_Sorted[
n].push_back(*
r);
880 if ((*r)->spacepoint->clusterList().second == 0 &&
z>=3 &&
z<=7) {
881 z<=4 ? z=0 : z>=6 ?
z=2 :
z=1;
891 data.rfzv_Sorted[
n].push_back(*
r);
894 data.r_Sorted[
i].erase(
r++);
◆ fillSeeds()
void InDet::SiSpacePointsSeedMaker_Trigger::fillSeeds |
( |
EventData & |
data | ) |
|
|
staticprivate |
Definition at line 1632 of file SiSpacePointsSeedMaker_Trigger.cxx.
1635 l =
data.mapOneSeeds.begin(),
1636 le =
data.mapOneSeeds.end ();
1638 for (;
l!=le; ++
l) {
1639 float q = (*l).first;
1642 if ((*
s0->spacePoints().rbegin())->clusterList().second) {
1643 (*
s0->spacePoints().begin())->clusterList().second ?
q+=1000. :
q+=10000.;
1646 if (
data.i_seede!=
data.l_seeds.end()) {
1649 data.mapSeeds.insert(std::make_pair(
q,
s));
1651 data.l_seeds.emplace_back(*s0);
1654 data.mapSeeds.insert(std::make_pair(
q,
s));
◆ finalize()
StatusCode InDet::SiSpacePointsSeedMaker_Trigger::finalize |
( |
| ) |
|
|
overridevirtual |
◆ find2Sp()
◆ find3Sp() [1/2]
◆ find3Sp() [2/2]
◆ findNext()
void InDet::SiSpacePointsSeedMaker_Trigger::findNext |
( |
EventData & |
data | ) |
const |
|
private |
◆ findVSp()
◆ getWriteNtupleBoolProperty()
bool InDet::SiSpacePointsSeedMaker_Trigger::getWriteNtupleBoolProperty |
( |
| ) |
const |
|
overridevirtual |
◆ initialize()
StatusCode InDet::SiSpacePointsSeedMaker_Trigger::initialize |
( |
| ) |
|
|
overridevirtual |
◆ initializeEventData()
void InDet::SiSpacePointsSeedMaker_Trigger::initializeEventData |
( |
EventData & |
data | ) |
const |
|
private |
◆ isZCompatible()
bool InDet::SiSpacePointsSeedMaker_Trigger::isZCompatible |
( |
EventData & |
data, |
|
|
float & |
Zv, |
|
|
float & |
R, |
|
|
float & |
T |
|
) |
| const |
|
private |
◆ newEvent()
void InDet::SiSpacePointsSeedMaker_Trigger::newEvent |
( |
const EventContext & |
ctx, |
|
|
EventData & |
data, |
|
|
int |
iteration |
|
) |
| const |
|
overridevirtual |
Definition at line 82 of file SiSpacePointsSeedMaker_Trigger.cxx.
91 double f[3], gP[3] ={10.,10.,0.};
98 if (fieldCondObj ==
nullptr) {
102 fieldCondObj->getInitializedCache (fieldCache);
107 data.K = 2./(300.*
f[2]);
109 data.K = 2./(300.* 5.);
112 data.i_spforseed =
data.l_spforseed.begin();
122 if (spacepointsPixel.isValid()) {
132 int ir =
static_cast<int>(sps->
radius()*irstep);
133 if (
ir>irmax)
ir = irmax;
134 data.r_Sorted[
ir].push_back(sps);
148 if (spacepointsSCT.isValid()) {
158 int ir =
static_cast<int>(sps->
radius()*irstep);
159 if (
ir>irmax)
ir = irmax;
160 data.r_Sorted[
ir].push_back(sps);
173 if (spacepointsOverlap.isValid()) {
182 int ir =
static_cast<int>(sps->
radius()*irstep);
183 if (
ir>irmax)
ir = irmax;
184 data.r_Sorted[
ir].push_back(sps);
◆ newOneSeed()
Definition at line 1505 of file SiSpacePointsSeedMaker_Trigger.cxx.
1512 data.OneSeeds[
data.nOneSeeds].erase();
1516 data.OneSeeds[
data.nOneSeeds].setZVertex(
static_cast<double>(
z));
1517 data.mapOneSeeds.insert(std::make_pair(
q, &(
data.OneSeeds[
data.nOneSeeds])));
1520 std::multimap<float,InDet::SiSpacePointsSeed*>::reverse_iterator
1521 l =
data.mapOneSeeds.rbegin();
1522 if ((*l).first <=
q)
return;
1529 s->setZVertex(
static_cast<double>(
z));
1531 i =
data.mapOneSeeds.insert(std::make_pair(
q,
s));
1533 for (++
i;
i!=
data.mapOneSeeds.end(); ++
i) {
1534 if ((*i).second==
s) {
1535 data.mapOneSeeds.erase(
i);
◆ newRegion() [1/2]
Definition at line 199 of file SiSpacePointsSeedMaker_Trigger.cxx.
204 data.trigger =
false;
210 double f[3], gP[3] ={10.,10.,0.};
217 if (fieldCondObj ==
nullptr) {
221 fieldCondObj->getInitializedCache (fieldCache);
226 data.K = 2./(300.*
f[2]);
228 data.K = 2./(300.* 5.);
231 data.i_spforseed =
data.l_spforseed.begin();
238 if (
m_pixel && !vPixel.empty()) {
241 if (spacepointsPixel.isValid()) {
246 const auto *
w = spacepointsPixel->indexFindPtr(
l);
247 if (
w==
nullptr)
continue;
252 int ir =
static_cast<int>(sps->
radius()*irstep);
253 if (
ir>irmax)
ir = irmax;
254 data.r_Sorted[
ir].push_back(sps);
265 if (
m_sct && !vSCT.empty()) {
268 if (spacepointsSCT.isValid()) {
273 const auto *
w = spacepointsSCT->indexFindPtr(
l);
274 if (
w==
nullptr)
continue;
279 int ir =
static_cast<int>(sps->
radius()*irstep);
280 if (
ir>irmax)
ir = irmax;
281 data.r_Sorted[
ir].push_back(sps);
◆ newRegion() [2/2]
Definition at line 296 of file SiSpacePointsSeedMaker_Trigger.cxx.
307 double dzdrmin = 1./
tan(2.*
atan(
exp(-IRD.etaMinus())));
308 double dzdrmax = 1./
tan(2.*
atan(
exp(-IRD.etaPlus ())));
310 data.zminB = IRD.zedMinus()-
data.zbeam[0];
311 data.zmaxB = IRD.zedPlus ()-
data.zbeam[0];
312 data.zminU =
data.zminB+550.f*dzdrmin;
313 data.zmaxU =
data.zmaxB+550.f*dzdrmax;
314 double fmax = IRD.phiPlus ();
315 double fmin = IRD.phiMinus();
316 if (fmin > fmax) fmin-=(2.*
M_PI);
317 data.ftrig = (fmin+fmax)*.5;
318 data.ftrigW = (fmax-fmin)*.5;
◆ newSeed()
◆ newSpacePoint()
◆ newVertices()
◆ next()
◆ operator=()
◆ production2Sp()
void InDet::SiSpacePointsSeedMaker_Trigger::production2Sp |
( |
EventData & |
data | ) |
const |
|
private |
Definition at line 938 of file SiSpacePointsSeedMaker_Trigger.cxx.
940 if (
data.nsazv<2)
return;
956 if (!
data.rfzv_map[
a])
continue;
957 r0 =
data.rfzv_Sorted[
a].begin();
958 r0e =
data.rfzv_Sorted[
a].end ();
967 for (;
r0!=r0e; ++
r0) {
969 float X = (*r0)->x();
970 float Y = (*r0)->y();
971 float R = (*r0)->radius();
974 float Z = (*r0)->z();
981 for (
int i=0;
i!=NB; ++
i) {
984 if (!
data.rfzv_map[an])
continue;
986 r =
data.rfzv_Sorted[an].begin();
987 re =
data.rfzv_Sorted[an].end ();
991 float Rb =(*r)->radius();
997 float dZ =
Z-(*r)->z();
999 if (Tz<m_dzdrmin || Tz>
m_dzdrmax)
continue;
1008 float dx =(*r)->x()-
X;
1009 float dy =(*r)->y()-
Y;
1012 float xy =
x*
x+
y*
y;
if (xy == 0.)
continue;
1016 float UR = Ut*
R+1.;
if (UR == 0.)
continue;
1019 if (std::abs(
B*
data.K) >
m_ipt*std::sqrt(1.f+
A*
A))
continue;
1021 newSeed(
data, (*r)->spacepoint, (*r0)->spacepoint, Zo);
1033 data.endlist =
true;
◆ production3Sp() [1/2]
void InDet::SiSpacePointsSeedMaker_Trigger::production3Sp |
( |
EventData & |
data | ) |
const |
|
private |
Definition at line 1040 of file SiSpacePointsSeedMaker_Trigger.cxx.
1042 data.mapSeeds.clear();
1043 if (
data.nsaz<3)
return;
1045 const int ZI[
SizeZ]= {5,6,7,8,9,10,4,3,2,1,0};
1061 if (!
data.rfz_map[
a])
continue;
1066 if (!
data.rfz_map[an])
continue;
1067 rb [NB] =
data.rfz_Sorted[an].begin();
1068 rbe[NB++] =
data.rfz_Sorted[an].end();
1073 if (!
data.rfz_map[an])
continue;
1074 rt [NT] =
data.rfz_Sorted[an].begin();
1075 rte[NT++] =
data.rfz_Sorted[an].end();
1079 if (!
data.endlist) {
1086 data.endlist =
true;
◆ production3Sp() [2/2]
Definition at line 1093 of file SiSpacePointsSeedMaker_Trigger.cxx.
1102 if (!
data.endlist) {
1104 data.endlist =
true;
1109 for (;
r0!=rbe[0]; ++
r0) {
1112 data.mapOneSeeds.erase(
data.mapOneSeeds.begin(),
data.mapOneSeeds.end());
1114 float R = (*r0)->radius();
1123 float X = (*r0)->x();
1124 float Y = (*r0)->y();
1125 float Z = (*r0)->z();
1130 for (
int i=0;
i!=NB; ++
i) {
1132 for (
r=
rb[
i];
r!=rbe[
i]; ++
r) {
1134 float Rb =(*r)->radius();
1144 if (dR >
m_drmax || (*r)->sur()==sur0)
continue;
1146 if ( !
pix && !(*r)->spacepoint->clusterList().second)
continue;
1148 float Tz = (
Z-(*r)->z())/dR;
1150 if (Tz < m_dzdrmin || Tz >
m_dzdrmax)
continue;
1166 for (
int i=0;
i!=NT; ++
i) {
1168 for (
r=rt[
i];
r!=rte[
i]; ++
r) {
1170 float Rt =(*r)->radius();
1178 if ( (*r)->sur()==sur0)
continue;
1180 float Tz = ((*r)->z()-
Z)/dR;
1182 if (Tz < m_dzdrmin || Tz >
m_dzdrmax)
continue;
1194 if (!(Nt-Nb))
continue;
1196 float covr0 = (*r0)->covr ();
1197 float covz0 = (*r0)->covz ();
1202 for (
int i=0;
i!=Nt; ++
i) {
1206 float dx = sp->
x()-
X;
1207 float dy = sp->
y()-
Y;
1208 float dz = sp->
z()-
Z;
1211 float r2 = 1.f/(
x*
x+
y*
y);
1212 float dr = std::sqrt(r2);
1214 if (
i < Nb) tz = -tz;
1230 float ipt2K = ipt2/K2;
1231 float ipt2C = ipt2*COF;
1232 float COFK = COF*K2;
1238 for (
int b=0;
b!=Nb; ++
b) {
1242 float Zob =
data.Zo[
b];
1243 float Tzb =
data.Tz[
b];
1244 float Rb2r =
data.R [
b]*covr0;
1245 float Rb2z =
data.R [
b]*covz0;
1246 float Erb =
data.Er[
b];
1247 float Vb =
data.V [
b];
1248 float Ub =
data.U [
b];
1249 float Tzb2 = (1.f+Tzb*Tzb);
1250 float CSA = Tzb2*COFK;
1251 float ICSA = Tzb2*ipt2C;
1253 float Iz = (dZ*dZ)/Tzb2;
1255 for (
int t=Nb;
t!=Nt; ++
t) {
1257 float Ts = .5*(Tzb+
data.Tz[
t]);
1260 if ( dT > ICSA)
continue;
1261 float dU =
data.U[
t]-Ub;
1262 if (dU == 0.)
continue;
1263 float A = (
data.V[
t]-Vb)/dU;
1267 if (B2 > ipt2K*
S2 || dT*
S2 > B2*CSA)
continue;
1268 float Im = std::abs((
A-
B*
R)*
R);
1270 if (Im > imc )
continue;
1273 if (
data.SP[
t]->spacepoint->clusterList().second && Im > imcs)
continue;
1280 nseed +=
data.mapOneSeeds.size();
◆ production3SpTrigger()
Definition at line 1295 of file SiSpacePointsSeedMaker_Trigger.cxx.
1303 constexpr
float pi2 = 2.*
M_PI;
1306 if (!
data.endlist) {
1308 data.endlist =
true;
1313 for (;
r0!=rbe[0]; ++
r0) {
1316 data.mapOneSeeds.erase(
data.mapOneSeeds.begin(),
data.mapOneSeeds.end());
1318 float R = (*r0)->radius();
1327 float X = (*r0)->x();
1328 float Y = (*r0)->y();
1329 float Z = (*r0)->z();
1334 for (
int i=0;
i!=NB; ++
i) {
1336 for (
r=
rb[
i];
r!=rbe[
i]; ++
r) {
1338 float Rb =(*r)->radius();
1348 if (dR >
m_drmax || (*r)->sur()==sur0)
continue;
1350 if ( !
pix && !(*r)->spacepoint->clusterList().second)
continue;
1354 float Tz = (
Z-(*r)->z())/dR;
1356 if (Zo < data.zminB || Zo >
data.zmaxB)
continue;
1357 float Zu =
Z+(550.f-
R)*Tz;
1358 if (Zu < data.zminU || Zu >
data.zmaxU)
continue;
1369 for (
int i=0;
i!=NT; ++
i) {
1371 for (
r=rt[
i];
r!=rte[
i]; ++
r) {
1373 float Rt =(*r)->radius();
1381 if ( (*r)->sur()==sur0)
continue;
1385 float Tz = ((*r)->z()-
Z)/dR;
1387 if (Zo < data.zminB || Zo >
data.zmaxB)
continue;
1388 float Zu =
Z+(550.f-
R)*Tz;
1389 if (Zu < data.zminU || Zu >
data.zmaxU)
continue;
1396 if (!(Nt-Nb))
continue;
1398 float covr0 = (*r0)->covr ();
1399 float covz0 = (*r0)->covz ();
1404 for (
int i=0;
i!=Nt; ++
i) {
1408 float dx = sp->
x()-
X;
1409 float dy = sp->
y()-
Y;
1410 float dz = sp->
z()-
Z;
1413 float r2 = 1.f/(
x*
x+
y*
y);
1414 float dr = std::sqrt(r2);
1416 if (
i < Nb) tz = -tz;
1432 float ipt2K = ipt2/K2;
1433 float ipt2C = ipt2*COF;
1434 float COFK = COF*K2;
1440 for (
int b=0;
b!=Nb; ++
b) {
1444 float Zob =
data.Zo[
b];
1445 float Tzb =
data.Tz[
b];
1446 float Rb2r =
data.R [
b]*covr0;
1447 float Rb2z =
data.R [
b]*covz0;
1448 float Erb =
data.Er[
b];
1449 float Vb =
data.V [
b];
1450 float Ub =
data.U [
b];
1451 float Tzb2 = (1.f+Tzb*Tzb);
1452 float CSA = Tzb2*COFK;
1453 float ICSA = Tzb2*ipt2C;
1455 float Iz = (dZ*dZ)/Tzb2;
1457 for (
int t=Nb;
t!=Nt; ++
t) {
1459 float Ts = .5f*(Tzb+
data.Tz[
t]);
1462 if ( dT > ICSA)
continue;
1463 float dU =
data.U[
t]-Ub;
1464 if (dU == 0.)
continue;
1465 float A = (
data.V[
t]-Vb)/dU;
1469 if (B2 > ipt2K*
S2 || dT*
S2 > B2*CSA)
continue;
1470 float Im = std::abs((
A-
B*
R)*
R);
1472 if (Im > imc )
continue;
1474 if (
data.SP[
t]->spacepoint->clusterList().second && Im > imcs)
continue;
1482 float x = 2.f*
B*
R-
A;
1485 if (
df >
data.ftrigW)
continue;
1490 nseed +=
data.mapOneSeeds.size();
◆ writeNtuple()
◆ m_beamSpotKey
◆ m_COF
float InDet::SiSpacePointsSeedMaker_Trigger::m_COF {0.} |
|
private |
◆ m_diver
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_diver {this, "maxdImpact", 10.} |
|
private |
◆ m_diverpps
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_diverpps {this, "maxdImpactPPS", 1.7} |
|
private |
◆ m_diversss
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_diversss {this, "maxdImpactSSS", 1000.} |
|
private |
◆ m_drmax
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_drmax {this, "maxdRadius", 270.} |
|
private |
◆ m_drmin
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_drmin {this, "mindRadius", 10.} |
|
private |
◆ m_drminv
float InDet::SiSpacePointsSeedMaker_Trigger::m_drminv {20.} |
|
private |
◆ m_dzdrmax
float InDet::SiSpacePointsSeedMaker_Trigger::m_dzdrmax {0.} |
|
private |
◆ m_dzdrmin
float InDet::SiSpacePointsSeedMaker_Trigger::m_dzdrmin {0.} |
|
private |
◆ m_dzdrver
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_dzdrver {this, "maxdZdRver", 0.02} |
|
private |
◆ m_dzver
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_dzver {this, "maxdZver", 5.} |
|
private |
◆ m_etamax
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_etamax {this, "etaMax", 2.7} |
|
private |
◆ m_fieldCondObjInputKey
Initial value:{this, "AtlasFieldCacheCondObj", "fieldCondObj",
"Name of the Magnetic Field conditions object key"}
Definition at line 143 of file SiSpacePointsSeedMaker_Trigger.h.
◆ m_fNmax
int InDet::SiSpacePointsSeedMaker_Trigger::m_fNmax {0} |
|
private |
◆ m_fvNmax
int InDet::SiSpacePointsSeedMaker_Trigger::m_fvNmax {0} |
|
private |
◆ m_initialized
bool InDet::SiSpacePointsSeedMaker_Trigger::m_initialized {false} |
|
private |
◆ m_ipt
float InDet::SiSpacePointsSeedMaker_Trigger::m_ipt {0.} |
|
private |
◆ m_ipt2
float InDet::SiSpacePointsSeedMaker_Trigger::m_ipt2 {0.} |
|
private |
◆ m_maxNumberVertices
UnsignedIntegerProperty InDet::SiSpacePointsSeedMaker_Trigger::m_maxNumberVertices {this, "maxNumberVertices", 99} |
|
private |
◆ m_maxOneSize
IntegerProperty InDet::SiSpacePointsSeedMaker_Trigger::m_maxOneSize {this, "maxSeedsForSpacePoint", 3} |
|
private |
◆ m_maxsize
IntegerProperty InDet::SiSpacePointsSeedMaker_Trigger::m_maxsize {this, "maxSize", 20000} |
|
private |
◆ m_maxsizeSP
IntegerProperty InDet::SiSpacePointsSeedMaker_Trigger::m_maxsizeSP {this, "maxSizeSP", 1500} |
|
private |
◆ m_outputlevel
int InDet::SiSpacePointsSeedMaker_Trigger::m_outputlevel {0} |
|
private |
◆ m_pixel
BooleanProperty InDet::SiSpacePointsSeedMaker_Trigger::m_pixel {this, "usePixel", true} |
|
private |
◆ m_ptmin
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_ptmin {this, "pTmin", 500.} |
|
private |
◆ m_r1max
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r1max {this, "maxRadius1", 600.} |
|
private |
◆ m_r1maxv
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r1maxv {this, "maxVRadius1", 60.} |
|
private |
◆ m_r1min
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r1min {this, "minRadius1", 0.} |
|
private |
◆ m_r1minv
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r1minv {this, "minVRadius1", 0.} |
|
private |
◆ m_r2max
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r2max {this, "maxRadius2", 600.} |
|
private |
◆ m_r2maxv
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r2maxv {this, "maxVRadius2", 200.} |
|
private |
◆ m_r2min
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r2min {this, "minRadius2", 0.} |
|
private |
◆ m_r2minv
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r2minv {this, "minVRadius2", 70.} |
|
private |
◆ m_r3max
float InDet::SiSpacePointsSeedMaker_Trigger::m_r3max {0.} |
|
private |
◆ m_r3min
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r3min {this, "minRadius3", 0.} |
|
private |
◆ m_r_rmax
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r_rmax {this, "radMax", 600.} |
|
private |
◆ m_r_rstep
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_r_rstep {this, "radStep", 2.} |
|
private |
◆ m_r_size
int InDet::SiSpacePointsSeedMaker_Trigger::m_r_size {0} |
|
private |
◆ m_rfz_b
int InDet::SiSpacePointsSeedMaker_Trigger::m_rfz_b[SizeRFZ] {} |
|
private |
◆ m_rfz_ib
int InDet::SiSpacePointsSeedMaker_Trigger::m_rfz_ib[SizeRFZ][SizeI] {} |
|
private |
◆ m_rfz_it
int InDet::SiSpacePointsSeedMaker_Trigger::m_rfz_it[SizeRFZ][SizeI] {} |
|
private |
◆ m_rfz_t
int InDet::SiSpacePointsSeedMaker_Trigger::m_rfz_t[SizeRFZ] {} |
|
private |
◆ m_rfzv_i
int InDet::SiSpacePointsSeedMaker_Trigger::m_rfzv_i[SizeRFZV][SizeIV] {} |
|
private |
◆ m_rfzv_n
int InDet::SiSpacePointsSeedMaker_Trigger::m_rfzv_n[SizeRFZV] {} |
|
private |
◆ m_sct
BooleanProperty InDet::SiSpacePointsSeedMaker_Trigger::m_sct {this, "useSCT", true} |
|
private |
◆ m_sF
float InDet::SiSpacePointsSeedMaker_Trigger::m_sF {0.} |
|
private |
◆ m_sFv
float InDet::SiSpacePointsSeedMaker_Trigger::m_sFv {0.} |
|
private |
◆ m_spacepointsOverlap
◆ m_spacepointsPixel
◆ m_spacepointsSCT
◆ m_useOverlap
BooleanProperty InDet::SiSpacePointsSeedMaker_Trigger::m_useOverlap {this, "useOverlapSpCollection", true} |
|
private |
◆ m_zmax
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_zmax {this, "maxZ", +250.} |
|
private |
◆ m_zmin
FloatProperty InDet::SiSpacePointsSeedMaker_Trigger::m_zmin {this, "minZ", -250.} |
|
private |
The documentation for this class was generated from the following files:
const std::pair< const PrepRawData *, const PrepRawData * > & clusterList() const
return the pair of cluster pointers by reference
JetConstituentVector::iterator iterator
int cost(std::vector< std::string > &files, node &n, const std::string &directory="", bool deleteref=false, bool relocate=false)
void buildBeamFrameWork(EventData &data) const
char data[hepevt_bytes_allocation_ATLAS]
def finalize(self)
_info( "content of StoreGate..." ) self.sg.dump()
static MsgStream & dumpEvent(EventData &data, MsgStream &out)
virtual const Amg::Vector3D & globalPosition() const override final
Interface method to get the global Position.
void getFieldZR(const double *ATH_RESTRICT xyz, double *ATH_RESTRICT bxyz, double *ATH_RESTRICT deriv=nullptr)
get B field valaue on the z-r plane at given position works only inside the solenoid.
static void newSeed(EventData &data, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const float &)
const float & covz() const
MsgStream & dumpConditions(EventData &data, MsgStream &out) const
void newOneSeed(EventData &data, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const float &, const float &) const
const std::string & key() const
Return the StoreGate ID for the referenced object.
bool isZCompatible(EventData &data, float &, float &, float &) const
SG::ReadCondHandleKey< AtlasFieldCacheCondObj > m_fieldCondObjInputKey
static void convertToBeamFrameWork(EventData &data, const Trk::SpacePoint *const &, float *)
void production3Sp(EventData &data) const
virtual void find3Sp(const EventContext &ctx, EventData &data, const std::list< Trk::Vertex > &lv) const override
with three space points with or without vertex constraint
const float & radius() const
bool newVertices(EventData &data, const std::list< Trk::Vertex > &) const
virtual void newRegion(const EventContext &ctx, EventData &data, const std::vector< IdentifierHash > &vPixel, const std::vector< IdentifierHash > &vSCT) const override
SG::ReadHandleKey< SpacePointContainer > m_spacepointsSCT
void set(const Trk::SpacePoint *const &, const float *)
SG::ReadHandleKey< SpacePointOverlapCollection > m_spacepointsOverlap
void initializeEventData(EventData &data) const
void production2Sp(EventData &data) const
static void erase(EventData &data)
::StatusCode StatusCode
StatusCode definition for legacy code.
const double r0
electron radius{cm}
virtual MsgStream & dump(EventData &data, MsgStream &out) const override
IntegerProperty m_maxsize
double R(const INavigable4Momentum *p1, const double v_eta, const double v_phi)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
void fillLists(EventData &data) const
bool solenoidOn() const
status of the magnets
StatusCode initialize(bool used=true)
Eigen::Matrix< double, 3, 1 > Vector3D
static float dZVertexMin(EventData &data, float &)
const float & covr() const
int ir
counter of the current depth
SG::ReadCondHandleKey< InDet::BeamSpotData > m_beamSpotKey
int m_rfz_it[SizeRFZ][SizeI]
UnsignedIntegerProperty m_maxNumberVertices
static SiSpacePointForSeed * newSpacePoint(EventData &data, const Trk::SpacePoint *const &)
const boost::regex re(r_e)
int m_rfz_ib[SizeRFZ][SizeI]
Local cache for magnetic field (based on MagFieldServices/AtlasFieldSvcTLS.h)
struct TBPatternUnitContext S2
BooleanProperty m_useOverlap
void findNext(EventData &data) const
IntegerProperty m_maxsizeSP
int m_rfzv_i[SizeRFZV][SizeIV]
IntegerProperty m_maxOneSize
This is a "hash" representation of an Identifier. This encodes a 32 bit index which can be used to lo...
void production3SpTrigger(EventData &data, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, int, int, int &) const
static void fillSeeds(EventData &data)
SG::ReadHandleKey< SpacePointContainer > m_spacepointsPixel