19 , m_inputKey(
"PixelClusters")
50 , m_siDetNtInPixLay(0)
52 , m_siDetHitDepthDir(0)
101 , m_h_siDetInPixLay(0)
102 , m_h_siDetNtInPixLay(0)
104 , m_h_siDetHitDepthDir(0)
105 , m_h_siDetHitPhiDir(0)
106 , m_h_siDetHitEtaDir(0)
114 , m_h_siDetMinWidth(0)
115 , m_h_siDetMaxWidth(0)
118 , m_h_siDetEtaPitch(0)
119 , m_h_siDetPhiPitch(0)
126 , m_ntupleFileName(
"/ntuples/file1")
127 , m_ntupleDirName(
"/PixelFastRDOAnalysis/")
128 , m_ntupleTreeName(
"/PixelFastRDOAna")
129 ,
m_path(
"/PixelFastRDOAnalysis/")
130 , m_thistSvc(
"THistSvc",
name)
420 return StatusCode::SUCCESS;
484 if(p_pixelClus_cont.
isValid()) {
486 InDet::PixelClusterContainer::const_iterator clusCont_itr(p_pixelClus_cont->begin());
487 const InDet::PixelClusterContainer::const_iterator clusCont_end(p_pixelClus_cont->end());
488 for ( ; clusCont_itr != clusCont_end; ++clusCont_itr ) {
490 InDet::PixelClusterCollection::const_iterator clus_itr(p_pixelClus_coll->begin());
491 const InDet::PixelClusterCollection::const_iterator clus_end(p_pixelClus_coll->end());
492 for ( ; clus_itr != clus_end; ++clus_itr ) {
494 const float omegax((*clus_itr)->omegax());
495 const float omegay((*clus_itr)->omegay());
496 const std::vector<int>& tot_v((*clus_itr)->totList());
497 const int totalTot((*clus_itr)->totalToT());
498 const std::vector<float>& charge_v((*clus_itr)->chargeList());
499 const float totalCharge((*clus_itr)->totalCharge());
500 const float eLoss((*clus_itr)->energyLoss());
501 const bool fake((*clus_itr)->isFake());
502 const bool ambig((*clus_itr)->isAmbiguous());
503 const bool split((*clus_itr)->isSplit());
504 const double splitProb1((*clus_itr)->splitProbability1());
505 const double splitProb2((*clus_itr)->splitProbability2());
507 const int LVL1A((*clus_itr)->LVL1A());
511 for (
size_t i = 0;
i != tot_v.size(); ++
i) {
512 m_tot->push_back(tot_v.at(
i));
516 for (
size_t j = 0; j != charge_v.size(); ++j) {
517 m_charge->push_back(charge_v.at(j));
546 const float siCol(siColRow.x());
547 const float siRow(siColRow.y());
549 const float siPhiR(siPhiRZ.x());
550 const float siZ(siPhiRZ.y());
552 const Amg::Vector3D& siPos((*clus_itr)->SiCluster::globalPosition());
553 const float siPos_x(siPos.x());
554 const float siPos_y(siPos.y());
555 const float siPos_z(siPos.z());
557 const bool siGangPix((*clus_itr)->SiCluster::gangedPixel());
560 const Identifier siDetID((*siDetEl).identify());
562 const bool siDetPix((*siDetEl).isPixel());
563 const bool siDetSCT((*siDetEl).isSCT());
564 const bool siDetBrl((*siDetEl).isBarrel());
565 const bool siDetEc((*siDetEl).isEndcap());
566 const bool siDetBlay((*siDetEl).isBlayer());
567 const bool siDetInPixLay((*siDetEl).isInnermostPixelLayer());
568 const bool siDetNtInPixLay((*siDetEl).isNextToInnermostPixelLayer());
569 const bool siDetDBM((*siDetEl).isDBM());
570 const double siDetHitDepthDir((*siDetEl).hitDepthDirection());
571 const double siDetHitPhiDir((*siDetEl).hitPhiDirection());
572 const double siDetHitEtaDir((*siDetEl).hitEtaDirection());
573 const double siDetMinR((*siDetEl).rMin());
574 const double siDetMaxR((*siDetEl).rMax());
575 const double siDetMinZ((*siDetEl).zMin());
576 const double siDetMaxZ((*siDetEl).zMax());
577 const double siDetMinPhi((*siDetEl).phiMin());
578 const double siDetMaxPhi((*siDetEl).phiMax());
579 const double siDetWidth((*siDetEl).width());
580 const double siDetMinWidth((*siDetEl).minWidth());
581 const double siDetMaxWidth((*siDetEl).maxWidth());
582 const double siDetLength((*siDetEl).length());
583 const double siDetThick((*siDetEl).thickness());
584 const double siDetEtaPitch((*siDetEl).etaPitch());
585 const double siDetPhiPitch((*siDetEl).phiPitch());
587 const unsigned long long siDetID_int = siDetID.
get_compact();
592 m_siZ->push_back(siZ);
658 const Identifier clusID((*clus_itr)->PrepRawData::identify());
659 const unsigned long long clusID_int = clusID.
get_compact();
661 const Amg::Vector2D& locPos((*clus_itr)->PrepRawData::localPosition());
662 float locpos_x(locPos.x());
663 float locpos_y(locPos.y());
665 const std::vector<Identifier>& rdoID_v((*clus_itr)->PrepRawData::rdoList());
670 for (
size_t k = 0;
k != rdoID_v.size(); ++
k) {
671 const unsigned long long rdoID_prd_int = rdoID_v.at(
k).get_compact();
687 return StatusCode::SUCCESS;
691 return StatusCode::SUCCESS;