|
ATLAS Offline Software
|
Go to the documentation of this file.
19 , m_inputKey(
"TRT_DriftCircles")
26 , m_timeOverThreshold(0)
76 , m_h_timeOverThreshold(0)
78 , m_h_driftTimeValid(0)
109 , m_h_allStrCntr_x(0)
110 , m_h_allStrCntr_y(0)
111 , m_h_allStrCntr_z(0)
112 , m_h_allStrCntr_r(0)
113 , m_h_allStrCntr_phi(0)
114 , m_h_allStrAxis_x(0)
115 , m_h_allStrAxis_y(0)
116 , m_h_allStrAxis_z(0)
117 , m_h_allStrAxis_r(0)
118 , m_h_allStrAxis_phi(0)
121 , m_ntupleFileName(
"/ntuples/file1")
122 , m_ntupleDirName(
"/TRT_FastRDOAnalysis/")
123 , m_ntupleTreeName(
"/TRT_FastRDOAna")
124 ,
m_path(
"/TRT_FastRDOAnalysis/")
125 , m_thistSvc(
"THistSvc",
name)
402 return StatusCode::SUCCESS;
463 InDet::TRT_DriftCircleContainer::const_iterator trtDCcont_itr(p_trtDC_cont->begin());
464 const InDet::TRT_DriftCircleContainer::const_iterator trtDCcont_end(p_trtDC_cont->end());
465 for ( ; trtDCcont_itr != trtDCcont_end; ++trtDCcont_itr ) {
467 InDet::TRT_DriftCircleCollection::const_iterator trtDC_itr(p_trtDC_coll->begin());
468 const InDet::TRT_DriftCircleCollection::const_iterator trtDC_end(p_trtDC_coll->end());
469 for ( ; trtDC_itr != trtDC_end; ++trtDC_itr ) {
471 const unsigned int word((*trtDC_itr)->getWord());
473 const int trailEdge((*trtDC_itr)->trailingEdge());
474 const bool HL((*trtDC_itr)->highLevel());
475 const bool firstBin((*trtDC_itr)->firstBinHigh());
476 const bool lastBin((*trtDC_itr)->lastBinHigh());
477 const double ToT((*trtDC_itr)->timeOverThreshold());
478 const double rawDriftTime((*trtDC_itr)->rawDriftTime());
479 const bool driftTimeValid((*trtDC_itr)->driftTimeValid());
480 const bool noise((*trtDC_itr)->isNoise());
509 const int brl_ec((*detElement).type());
510 const Identifier trtID((*detElement).identify());
511 const unsigned long long trtID_int = trtID.
get_compact();
516 const int surfType(
static_cast<int>(surf.
type()));
519 const int bndsType(bnds.
type());
522 const float cntr_x(cntr.x());
523 const float cntr_y(cntr.y());
524 const float cntr_z(cntr.z());
525 const float cntr_r(cntr.perp());
526 const float cntr_phi(cntr.phi());
529 const float norm_x(
norm.x());
530 const float norm_y(
norm.y());
531 const float norm_z(
norm.z());
532 const float norm_r(
norm.perp());
536 const Trk::Surface& strSurf((*detElement).surface(trtID));
537 const int strSurfType(
static_cast<int>(strSurf.
type()));
540 const int strBndsType(strBnds.
type());
543 const float strCntr_x(strCntr.x());
544 const float strCntr_y(strCntr.y());
545 const float strCntr_z(strCntr.z());
546 const float strCntr_r(strCntr.perp());
547 const float strCntr_phi(strCntr.phi());
550 const float strNorm_x(strNorm.x());
551 const float strNorm_y(strNorm.y());
552 const float strNorm_z(strNorm.z());
553 const float strNorm_r(strNorm.perp());
554 const float strNorm_phi(strNorm.phi());
556 const unsigned int strawN((*detElement).nStraws());
557 const double& strawL((*detElement).strawLength());
558 const int strawD((*detElement).strawDirection());
561 for (
unsigned int i = 0;
i != strawN; ++
i) {
563 const float strawCenter_x(strawCenter.x());
564 const float strawCenter_y(strawCenter.y());
565 const float strawCenter_z(strawCenter.z());
566 const float strawCenter_r(strawCenter.perp());
567 const float strawCenter_phi(strawCenter.phi());
570 const float strawAxis_x(strawAxis.x());
571 const float strawAxis_y(strawAxis.y());
572 const float strawAxis_z(strawAxis.z());
573 const float strawAxis_r(strawAxis.perp());
574 const float strawAxis_phi(strawAxis.phi());
669 return StatusCode::SUCCESS;
673 return StatusCode::SUCCESS;
std::vector< unsigned long long > * m_trtID
std::vector< float > * m_strCntr_x
std::vector< bool > * m_driftTimeValid
std::vector< float > * m_allStrAxis_phi
std::vector< float > * m_cntr_r
std::vector< double > * m_rawDriftTime
std::vector< float > * m_allStrCntr_r
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
std::vector< float > * m_allStrAxis_r
ServiceHandle< ITHistSvc > m_thistSvc
std::vector< bool > * m_highLevel
std::vector< float > * m_strNorm_phi
value_type get_compact() const
Get the compact id.
virtual StatusCode execute() override final
double norm_phi(const double input)
std::vector< int > * m_strSurfType
std::vector< float > * m_strNorm_x
std::vector< double > * m_strawD
TH1 * m_h_timeOverThreshold
std::vector< int > * m_surfType
std::vector< float > * m_cntr_phi
std::vector< float > * m_allStrCntr_phi
virtual BoundsType type() const =0
Return the bounds type - for persistency optimization.
std::vector< float > * m_norm_x
std::vector< float > * m_cntr_x
std::vector< float > * m_cntr_y
::StatusCode StatusCode
StatusCode definition for legacy code.
std::vector< float > * m_allStrCntr_z
std::vector< bool > * m_lastBin
std::vector< int > * m_bndsType
std::string m_path
the path being used
std::vector< float > * m_strCntr_r
virtual StatusCode finalize() override final
std::vector< float > * m_cntr_z
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
std::vector< float > * m_strCntr_phi
virtual bool isValid() override final
Can the handle be successfully dereferenced?
std::vector< float > * m_allStrCntr_x
double rawDriftTime(unsigned int m_word)
std::vector< bool > * m_noise
std::vector< int > * m_strawN
virtual StatusCode initialize() override final
std::vector< float > * m_norm_y
std::vector< float > * m_allStrAxis_x
Trk::PrepRawDataCollection< TRT_DriftCircle > TRT_DriftCircleCollection
std::vector< double > * m_strawL
std::vector< float > * m_norm_r
Eigen::Matrix< double, 3, 1 > Vector3D
std::string m_ntupleDirName
std::vector< float > * m_norm_phi
std::vector< float > * m_strCntr_z
std::vector< int > * m_brl_ec
std::vector< float > * m_strCntr_y
std::vector< int > * m_strBndsType
std::string m_ntupleTreeName
SG::ReadHandleKey< InDet::TRT_DriftCircleContainer > m_inputKey
def TH1F(name, title, nxbins, bins_par2, bins_par3=None, path='', **kwargs)
std::vector< float > * m_strNorm_r
Handle class for reading from StoreGate.
std::vector< int > * m_driftTimeBin
std::vector< float > * m_norm_z
std::vector< double > * m_timeOverThreshold
std::vector< float > * m_allStrAxis_y
TRT_FastRDOAnalysis(const std::string &name, ISvcLocator *pSvcLocator)
constexpr virtual SurfaceType type() const =0
Returns the Surface type to avoid dynamic casts.
std::vector< bool > * m_firstBin
std::vector< int > * m_trailEdge
std::vector< float > * m_allStrAxis_z
std::string m_ntupleFileName
std::vector< unsigned int > * m_word
std::vector< float > * m_strNorm_y
std::vector< float > * m_strNorm_z
int driftTimeBin(unsigned int m_word)
std::vector< float > * m_allStrCntr_y