CaloTTMgrDetDescrCnv Class Reference

This class is a converter for the CaloTTDescrManager which is stored in the detector store. More...

#include <CaloTTMgrDetDescrCnv.h>

Inheritance diagram for CaloTTMgrDetDescrCnv:

Collaboration diagram for CaloTTMgrDetDescrCnv:

Public Member Functions
virtual long int	repSvcType () const override
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual StatusCode	createObj (IOpaqueAddress *pAddr, DataObject *&pObj) override
	CaloTTMgrDetDescrCnv (ISvcLocator *svcloc)
virtual StatusCode	fillObjRefs (IOpaqueAddress *pAddr, DataObject *pObj)
virtual StatusCode	createRep (DataObject *pObj, IOpaqueAddress *&pAddr)
virtual StatusCode	fillRepRefs (IOpaqueAddress *pAddr, DataObject *pObj)
const ServiceHandle< StoreGateSvc > &	detStore () const
	Handle to DetectorStore.
void	addTrack (TIDA::Track *t)
const std::vector< TIDA::Track * > &	tracks () const
void	clear ()
void	selectTracks (const TrigInDetTrackCollection *trigtracks)
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 ()
static CLID	classID ()

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	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
ServiceHandle< StoreGateSvc >	m_detStore
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

This class is a converter for the CaloTTDescrManager which is stored in the detector store.

This class derives from DetDescrConverter which is a converter of the DetDescrCnvSvc. This converter creates a manager object and adds descriptors and detector elements to the manager. This objects are either created or accessed from the detector store.

Warning

The hadronic part of the FCAL calorimeter calorimeter is constructed with 4 channels in eta. These 4 channels actually correspond to 2 channels in eta FOR EACH of the hadronic samplings, i.e. FCAL2 and FCAL3. Therefore the geometry of these channels is wrong.

Definition at line 33 of file CaloTTMgrDetDescrCnv.h.

Constructor & Destructor Documentation

CaloTTMgrDetDescrCnv()

CaloTTMgrDetDescrCnv::CaloTTMgrDetDescrCnv

(

ISvcLocator *

svcloc

)

Definition at line 296 of file CaloTTMgrDetDescrCnv.cxx.

    :
    DetDescrConverter(ClassID_traits<CaloTTDescrManager>::ID(), svcloc)
{}

Member Function Documentation

addTrack()

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

inlineinherited

Definition at line 45 of file Converter.h.

                              {
    m_tracks.push_back(t);
  }

classID()

CLID CaloTTMgrDetDescrCnv::classID ( )

static

Definition at line 291 of file CaloTTMgrDetDescrCnv.cxx.

                              { 
    return ClassID_traits<CaloTTDescrManager>::ID(); 
}

clear()

void Converter::clear ( )

inlineinherited

Definition at line 53 of file Converter.h.

{m_tracks.clear();} 

createObj()

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

overridevirtual

Implements DetDescrConverter.

Definition at line 75 of file CaloTTMgrDetDescrCnv.cxx.

{
    ATH_MSG_INFO("in createObj: creating a CaloTTDescrManager object in the detector store");
 
    // Create a new CaloTTDescrManager
    CaloTTDescrManager* caloTTMgr = new CaloTTDescrManager(); 
 
    // Pass a pointer to the container to the Persistency service by reference.
    pObj = SG::asStorable(caloTTMgr);
 
    // Get idhelper from detector store and add to mgr
    const CaloLVL1_ID* lvl1_id = nullptr;
    ATH_CHECK(detStore()->retrieve(lvl1_id, "CaloLVL1_ID"));
    caloTTMgr->set_helper(lvl1_id);
    ATH_MSG_INFO("Set CaloLVL1_ID helper in CaloTTMgr ");
 
    // Get CaloDetDescrManager from condition store
    // to build geometry of trigger towers
    //
    // NB! The information retrieved from the CaloDetDescrManager is NOT sensitive
    //     to Calo alignment changes. Hence, it is OK to write CaloTTDescrManager
    //     into DetStore
    SG::ReadCondHandleKey<CaloDetDescrManager> caloMgrKey {"CaloDetDescrManager"};
    ATH_CHECK(caloMgrKey.initialize());
    SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{caloMgrKey};
    const CaloDetDescrManager* caloMgr = *caloMgrHandle;
 
    SmartIF<IToolSvc> toolSvc{service("ToolSvc")};
    ATH_CHECK(toolSvc.isValid());
 
    CaloTriggerTowerService*    ttSvc{nullptr};
    ATH_CHECK(toolSvc->retrieveTool("CaloTriggerTowerService",ttSvc));
 
    // Initialize the caloTT mgr
    //  We protect here in case this has been initialized elsewhere
 
    if (!caloTTMgr->is_initialized()) {
 
      ATH_MSG_DEBUG("Initializing CaloTTMgr from values in CaloTTMgrDetDescrCnv");
 
      int numberOfIdRegions=0;
      int numberOfDescrRegions=0;
      int nEmb=0;
      int nEmec=0;
      int nHec=0;
      int nFcal=0;
 
      // Initialize the manager ...
      const CaloLVL1_ID* lvl1_helper = lvl1_id;
      
      std::vector<Identifier>::const_iterator itId = lvl1_id->region_begin();
      std::vector<Identifier>::const_iterator itIdEnd = lvl1_id->region_end();
      
      for(; itId!=itIdEnd;++itId){
    Identifier regId = *itId;
    ++numberOfIdRegions;
    int posNeg   =  lvl1_helper->pos_neg_z(regId);
    //  int sampling = lvl1_helper->sampling(regId);
    int region   =  lvl1_helper->region(regId);
    int halfNphi = (lvl1_helper->phi_max(regId)+1)/2;
    int nLay     =  lvl1_helper->layer_max(regId)+1;
    
    // create Descriptors for both LAr and Tile
    CaloTTDescriptor* calo_descriptor = new CaloTTDescriptor() ;
    calo_descriptor->set(regId);
 
    if( 0 == region ) {
      //      calo_descriptor->set(0. ,2.5,0.1,-M_PI,M_PI,M_PI/halfNphi,posNeg,nLay);
      calo_descriptor->set(0. ,2.5,0.1,0.,2.*M_PI,M_PI/halfNphi,posNeg,nLay);
    } else if( 1 == region ) {
      //      calo_descriptor->set(2.5,3.1,0.2,-M_PI,M_PI,M_PI/halfNphi,posNeg,nLay);
      calo_descriptor->set(2.5,3.1,0.2,0.,2.*M_PI,M_PI/halfNphi,posNeg,nLay);
    } else if( 2 == region ) {
      //      calo_descriptor->set(3.1,3.2,0.1,-M_PI,M_PI,M_PI/halfNphi,posNeg,nLay);
      calo_descriptor->set(3.1,3.2,0.1,0.,2.*M_PI,M_PI/halfNphi,posNeg,nLay);
    } else if( 3 == region ) {
      //      calo_descriptor->set(3.2,4.9,0.4,-M_PI,M_PI,M_PI/halfNphi,posNeg,nLay);
      calo_descriptor->set(3.2,4.9,0.4,0.,2.*M_PI,M_PI/halfNphi,posNeg,nLay);
    }
 
    // Initialize 
    caloTTMgr->add(calo_descriptor);
    
    // now, create CaloTTDescrRegions for LAr only
    // Loop on TTs of this region and create CaloTTDescrRegion's (equiv. to cell's CaloDDE for TT)
    int maxEta=calo_descriptor->nEta();
    for(int iEta=0;iEta<maxEta;++iEta) {
        
      int maxPhi=calo_descriptor->nPhi();
      for(int iPhi=0;iPhi<maxPhi;++iPhi) {
          
        Identifier TTid = lvl1_helper->tower_id(regId,iEta,iPhi);
        if(!lvl1_helper->is_tile(TTid)) {
 
          int maxLay=calo_descriptor->nLay();
          for(int iLay=0;iLay<maxLay;++iLay) {
 
        Identifier layId = lvl1_helper->layer_id(TTid,iLay);
        bool isFCAL(0), isEC(0);
        if(lvl1_helper->is_fcal(layId)) {
          isFCAL=true;
        }
        else if(lvl1_helper->is_emec(layId) || lvl1_helper->is_hec(layId)) {
          isEC=true;
        }
        if (this->msgLvl(MSG::DEBUG)) {
          if(lvl1_helper->is_emb(layId) || lvl1_helper->is_barrel_end(layId) ) {
            ++nEmb;
             ATH_MSG_DEBUG(" Found EMB TT " << lvl1_helper->show_to_string(layId));
          }
          else if(lvl1_helper->is_emec(layId)) {
            ++nEmec;
             ATH_MSG_DEBUG(" Found EMEC TT " << lvl1_helper->show_to_string(layId));
          }
          else if(lvl1_helper->is_hec(layId)) {
            ++nHec;
             ATH_MSG_DEBUG(" Found HEC TT " << lvl1_helper->show_to_string(layId));
          }
          else {  // FCAL
            ++nFcal;
             ATH_MSG_DEBUG(" Found FCAL TT " << lvl1_helper->show_to_string(layId));
          }
        }
 
        double dEta=calo_descriptor->deta();
        double eta=calo_descriptor->eta_min()+(iEta+0.5)*dEta;
        eta*=calo_descriptor->sign_eta();
        
        double dPhi=calo_descriptor->dphi();
        // beware of the -pi,pi convention
        // fixed thanks to Denis Damazio when integrating LArRegionSelector
        //      float orig=-M_PI; 
        float orig=0.;
        double phi=calo_descriptor->phiMin() + (iPhi+0.5)*dPhi + orig;
        if(phi>=M_PI) {phi-=2.*M_PI;}
        if(phi<-M_PI) {phi+=2.*M_PI;}
        
        
        float rhoMin=99999.;
        float rhoMax=0.;
        float zMin=99999.;
        float zMax=-99999.;
        
        //int atlas_tb=0;
        std::vector<Identifier> vec = ttSvc->createCellIDvecLayer(layId);
        
        if(vec.size() > 0) {
          numberOfDescrRegions++;
          std::vector<Identifier>::const_iterator it     = vec.begin();
          std::vector<Identifier>::const_iterator it_end = vec.end();
          for(;it!=it_end;++it) {
            Identifier offId = (*it);
            const CaloDetDescrElement* caloDDE = 
                      caloMgr->get_element(offId);
            double rho  =caloDDE->r_raw();
            double cDrho=caloDDE->dr();
            double z    =caloDDE->z_raw();
            double cDz  =caloDDE->dz();
            
            if( (rho-cDrho/2.)<rhoMin ) {rhoMin=std::max(rho-cDrho/2.,0.);}
            if( (rho+cDrho/2.)>rhoMax ) {rhoMax=rho+cDrho/2.;}
            if( (z-cDz/2.)<zMin ) {zMin=z-cDz/2.;}
            if( (z+cDz/2.)>zMax ) {zMax=z+cDz/2.;}
          } // end loop on vector elements
        
        
          CaloTTDescrRegion* tt_region = new CaloTTDescrRegion(layId,calo_descriptor);
          if(isFCAL) {
            // FIX ME !?
            tt_region->set_cylindric(eta,phi,(zMax+zMin)/2.);
            tt_region->set_cylindric_size(dEta,dPhi,fabs(zMax-zMin));
          } else if(isEC) {
            tt_region->set_cylindric(eta,phi,(zMax+zMin)/2.);
            tt_region->set_cylindric_size(dEta,dPhi,fabs(zMax-zMin));
          } else {
            tt_region->set_spherical(eta,phi,(rhoMax+rhoMin)/2.);
            tt_region->set_spherical_size(dEta,dPhi,fabs(rhoMax-rhoMin));
          }
        
          caloTTMgr->add(tt_region);
 
        } // end condition of vec size
        else {
          ATH_MSG_DEBUG(" Found no cell for TT " << lvl1_helper->show_to_string(layId));
        }
          } // end loop on layers
        } // end condition against tile
      } // end loop on phi
    } // end loop on eta
      } // end loop on descr_regions
      
      // Set to initialized state only if descriptors have been found
      if (caloTTMgr->calo_descriptors_size () > 0) caloTTMgr->initialize();
    
      ATH_MSG_INFO(" Initialized CaloTTMgr, number of descr regions is " << numberOfDescrRegions);
      ATH_MSG_DEBUG(" including "
            << nEmb << " Em Barrel "
            << nEmec << " Em EC "
            << nHec << " HEC "
            << nFcal << " FCAL ");
      ATH_MSG_DEBUG(" number of helper regions= " << numberOfIdRegions);
    } // end of condition !is_initialized()
 
    return StatusCode::SUCCESS; 
}

createRep()

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

virtualinherited

Definition at line 20 of file DetDescrConverter.cxx.

                                                            {
    return StatusCode::SUCCESS;
}

detStore()

const ServiceHandle< StoreGateSvc > & DetDescrConverter::detStore ( ) const

inlineinherited

Handle to DetectorStore.

Definition at line 33 of file DetDescrConverter.h.

{ return m_detStore; }

fillObjRefs()

StatusCode DetDescrConverter::fillObjRefs ( IOpaqueAddress * pAddr, DataObject * pObj )

virtualinherited

Definition at line 15 of file DetDescrConverter.cxx.

                                                         {
    return StatusCode::SUCCESS;
}

fillRepRefs()

StatusCode DetDescrConverter::fillRepRefs ( IOpaqueAddress * pAddr, DataObject * pObj )

virtualinherited

Definition at line 25 of file DetDescrConverter.cxx.

                                                         {
    return StatusCode::SUCCESS;
}

finalize()

StatusCode CaloTTMgrDetDescrCnv::finalize ( )

overridevirtual

Definition at line 66 of file CaloTTMgrDetDescrCnv.cxx.

{
  ATH_MSG_DEBUG("in finalize");
  return StatusCode::SUCCESS;
}

initialize()

StatusCode CaloTTMgrDetDescrCnv::initialize ( )

overridevirtual

Definition at line 39 of file CaloTTMgrDetDescrCnv.cxx.

{
  ATH_MSG_DEBUG("in initialize");
  // First call parent init
  ATH_CHECK(DetDescrConverter::initialize());
    
    // The following is an attempt to "bootstrap" the loading of a
    // proxy for CaloTTDescrManager into the detector store. However,
    // CaloTTMgrDetDescrCnv::initialize is NOT called by the conversion
    // service.  So for the moment, this cannot be use. Instead the
    // DetDescrCnvSvc must do the bootstrap from a parameter list.
 
 
//      // Add CaloTTDescrManager proxy as entry point to the detector store
//      // - this is ONLY needed for the manager of each system
//      sc = addToDetStore(classID(), "CaloTTDescrManager");
//      if (sc.isFailure()) {
//      log << MSG::FATAL << "Unable to add proxy for CaloTTDescrManager to the Detector Store!" << endmsg;
//      return StatusCode::FAILURE;
//      } else {}
 
  return StatusCode::SUCCESS;
}

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.

{
  m_imsg = Athena::getMessageSvc();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

ipCorr()

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

inlineprotectedinherited

Definition at line 257 of file Converter.h.

                                                                                                  {
    
    double sn  = sin(phi0);
    double cs  = cos(phi0);
    double sd0 = (d0 != 0 ? d0/fabs(d0) : 1);
    double spt = (pt != 0 ? pt/fabs(pt) : 1);
    
    if (fabs(pt) >= 1*CLHEP::TeV) {
 
      d0c = d0 + m_beamX*sn - m_beamY*cs;
 
    } else {
      
      double rc = fabs(pt)*15.0/(9.0*1.042);
      
      double xc = (fabs(d0)-spt*sd0*rc)*cos(phi0+M_PI/2*sd0) - m_beamX;
      double yc = (fabs(d0)-spt*sd0*rc)*sin(phi0+M_PI/2*sd0) - m_beamY;
      
      double newphi;
      double xd01,yd01,xd02,yd02;
      
      if (xc == 0) {
    xd01 = 0; yd01 = rc+yc;
    xd02 = 0; yd02 = -rc+yc; 
      } else {
    xd01 = xc+yc/xc*yc+sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd01 = yc/xc*xd01;
    xd02 = xc+yc/xc*yc-sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd02 = yc/xc*xd02;
      }
      
      double r1 = sqrt(xd01*xd01+yd01*yd01);
      double r2 = sqrt(xd02*xd02+yd02*yd02);
      
      double phiV;
      
      if (r1 < r2)
    phiV = atan2(yd01,xd01); 
      else
    phiV = atan2(yd02,xd02); 
      
      double phi1 = phiCorr(phiV+M_PI/2);
      double phi2 = phiCorr(phiV-M_PI/2);
      
      if (fabs(phiCorr(phi1-phi0))<=fabs(phiCorr(phi2-phi0)))
    newphi = phi1;
      else
    newphi = phi2;
      
      d0c = fabs(sqrt(xc*xc+yc*yc)-rc)*sin(phiV-newphi);
      
      double theta=2*atan2(exp(-eta),1);
      double theta_save=theta;
      theta = theta - (1+spt)/2*M_PI;
      if (theta>0) theta = theta_save;
      
      double deltaz0= -spt*rc/tan(theta)*phiCorr(phiCorr(newphi)-phiCorr(phi0));
      z0c = z0 + deltaz0;
      
    }
  }

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.

{
  MsgStream* ms = m_msg_tls.get();
  if (!ms) {
    if (!m_initialized.test_and_set()) initMessaging();
    ms = new MsgStream(m_imsg,m_nm);
    m_msg_tls.reset( ms );
  }
 
  ms->setLevel (m_lvl);
  return *ms;
}

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.

{ return msg() << lvl; }

msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

{
  // If user did not set explicit message level we have to initialize
  // the messaging and retrieve the default via the MessageSvc.
  if (m_lvl==MSG::NIL && !m_initialized.test_and_set()) initMessaging();
 
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

phiCorr()

double Converter::phiCorr ( double phi )

inlineprotectedinherited

Definition at line 248 of file Converter.h.

                            {
    
    if (phi < -TMath::Pi()) phi += 2*TMath::Pi();
    if (phi >  TMath::Pi()) phi -= 2*TMath::Pi();
    
    return phi;
  }

repSvcType()

long int CaloTTMgrDetDescrCnv::repSvcType ( ) const

overridevirtual

Definition at line 33 of file CaloTTMgrDetDescrCnv.cxx.

{
  return (storageType());
}

selectTracks() [1/3]

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

inlineinherited

Definition at line 177 of file Converter.h.

                                                               {
    
    static int hpmap[20] = { 0, 1, 2,  7, 8, 9,  3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
    
    Analysis::MuonContainer::const_iterator muonitr = muontracks->begin();
    Analysis::MuonContainer::const_iterator muonend = muontracks->end();
    
    while(muonitr!=muonend) { 
      
      bool hasInDet = (*muonitr)->hasInDetTrackParticle();
      const Rec::TrackParticle *trackitr=NULL; 
      if(hasInDet) trackitr=(*muonitr)->inDetTrackParticle();
      muonitr++;
      if(!hasInDet) continue;
      
      const Trk::MeasuredPerigee* measPer = trackitr->measuredPerigee();
      
      CLHEP::HepVector perigeeParams = measPer->parameters();    
      double pT    = measPer->pT(); 
      double eta   = measPer->eta();
      double phi   = perigeeParams[Trk::phi0];
      double z0    = perigeeParams[Trk::z0];
      double d0    = perigeeParams[Trk::d0];
      //double theta = perigeeParams[Trk::theta];
      
      double deta = 1;
      double dphi = 1;
      double dz0  = 1;
      double dd0  = 1;
      double dpT  = 1;    
      
      // Check number of hits
      // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really 
      // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
      // stereo clusters making a spacepoint are two "hits"
      const Trk::TrackSummary *summary = trackitr->trackSummary();
      int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits); 
      int nPixelHits  = 2*summary->get(Trk::numberOfPixelHits);  
      int nSctHits    = summary->get(Trk::numberOfSCTHits); 
      int nStrawHits  = summary->get(Trk::numberOfTRTHits);
      int nTrHits     = summary->get(Trk::numberOfTRTHighThresholdHits);
      
      int nSiHits     = nPixelHits + nSctHits;
      
      const Trk::FitQuality *quality   = trackitr->fitQuality();
      double chi2 = quality->chiSquared();
      double dof  = quality->quality->numberDoF();
 
      unsigned bitmap = 0;
      
      for ( int ih=0 ; ih<20 ; ih++ ) {
    if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );    
      }
      
      // Shift d0 and z0 according to beam position
      ipCorr(d0, z0, d0, z0, phi, eta, pT);
      
      // Create and save Track  
      TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
                                 deta, dphi, dz0, dd0, dpT,
                                 nBlayerHits, nPixelHits, nSctHits, nSiHits,
                                 nStrawHits, nTrHits, bitmap, 0,
                                 -1) ;
      
      addTrack(t);
    }
    
  }

selectTracks() [2/3]

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

inlineinherited

Definition at line 111 of file Converter.h.

                                                                   { 
    
    Rec::TrackParticleContainer::const_iterator  trackitr = trigtracks->begin();
    Rec::TrackParticleContainer::const_iterator  trackend = trigtracks->end();
    
    static int hpmap[20] = { 0, 1, 2,  7, 8, 9,  3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
 
    while ( trackitr!=trackend ) { 
        
      const Trk::MeasuredPerigee* measPer = (*trackitr)->measuredPerigee();
 
      CLHEP::HepVector perigeeParams = measPer->parameters();    
      double pT    = measPer->pT(); 
      double eta   = measPer->eta();
      double phi   = perigeeParams[Trk::phi0];
      double z0    = perigeeParams[Trk::z0];
      double d0    = perigeeParams[Trk::d0];
      //double theta = perigeeParams[Trk::theta];
    
      double deta = 1;
      double dphi = 1;
      double dz0  = 1;
      double dd0  = 1;
      double dpT  = 1;    
 
      // Check number of hits
      // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really 
      // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
      // stereo clusters making a spacepoint are two "hits"
      const Trk::TrackSummary *summary = (*trackitr)->trackSummary();
      int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits); 
      int nPixelHits  = 2*summary->get(Trk::numberOfPixelHits);  
      int nSctHits    = summary->get(Trk::numberOfSCTHits); 
      int nStrawHits  = summary->get(Trk::numberOfTRTHits);
      int nTrHits     = summary->get(Trk::numberOfTRTHighThresholdHits);
      int nSiHits     = nPixelHits + nSctHits;
      bool expectBL   = false; // Not stored for Rec::TrackParticle
 
      const Trk::FitQuality *quality   = (*trackitr)->fitQuality();
      double chi2 = quality->chiSquared();
      double dof  = quality->numberDoF();
 
      unsigned bitmap = 0;
      
      for ( int ih=0 ; ih<20 ; ih++ ) {
    if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );    
      }
 
      // Shift d0 and z0 according to beam position
      ipCorr(d0, z0, d0, z0, phi, eta, pT);
 
      // Create and save Track  
      TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof, 
                                 deta, dphi, dz0, dd0, dpT,
                                 nBlayerHits, nPixelHits, nSctHits, nSiHits,
                                 nStrawHits, nTrHits, bitmap, 0,
                                 -1,
                                 expectBL) ;
  
      addTrack(t);
      trackitr++;
    }
    
  }

selectTracks() [3/3]

void Converter::selectTracks ( const TrigInDetTrackCollection * trigtracks )

inlineinherited

Definition at line 56 of file Converter.h.

                                                                {
    
    TrigInDetTrackCollection::const_iterator  trackitr = trigtracks->begin();
    TrigInDetTrackCollection::const_iterator  trackend = trigtracks->end();
 
    while ( trackitr!=trackend ) { 
 
    double eta    = (*trackitr)->param()->eta();
    double phi    = (*trackitr)->param()->phi0();
    double z0     = (*trackitr)->param()->z0(); 
    double pT     = (*trackitr)->param()->pT(); 
    double d0     = (*trackitr)->param()->a0();
    //double theta  = 2*atan2(exp(-(*trackitr)->param()->eta()),1);
 
    double deta    = (*trackitr)->param()->eeta();
    double dphi    = (*trackitr)->param()->ephi0();
    double dz0     = (*trackitr)->param()->ez0(); 
    double dpT     = (*trackitr)->param()->epT(); 
    double dd0     = (*trackitr)->param()->ea0();
 
    int   algoid  = (*trackitr)->algorithmId();           
 
    int nBlayerHits = ((*trackitr)->HitPattern() & 0x1);
    int nPixelHits  = 2*(*trackitr)->NPixelSpacePoints();  // NB: for comparison with offline 
    int nSctHits    = 2*(*trackitr)->NSCT_SpacePoints();   //     a spacepoint is 2 "hits"
    int nStrawHits  = (*trackitr)->NStrawHits();
    int nTrHits     = (*trackitr)->NTRHits();
               
    int nSiHits     = nPixelHits + nSctHits;
    bool expectBL   = false;                               //not available with TrigInDetTrack
 
    unsigned hitPattern = (*trackitr)->HitPattern();
    unsigned multiPattern = 0;
 
    double chi2    = (*trackitr)->chi2();
    double dof     = 0;
 
    // Shift d0 and z0 according to beam position
    ipCorr(d0, z0, d0, z0, phi, eta, pT);
    
    // Create and save Track  
    TIDA::Track* t = new TIDA::Track(eta,  phi,  z0,  d0,  pT, chi2, dof, 
                                   deta, dphi, dz0, dd0, dpT, 
                                   nBlayerHits, nPixelHits, nSctHits, nSiHits, 
                                   nStrawHits, nTrHits, hitPattern, multiPattern, 
                                   algoid,
                                   expectBL ) ; 
 
    addTrack(t);
    trackitr++; 
    }
  }

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.

{
  m_lvl = lvl;
}

storageType()

long CaloTTMgrDetDescrCnv::storageType ( )

static

Definition at line 284 of file CaloTTMgrDetDescrCnv.cxx.

{
    return DetDescr_StorageType;
}

tracks()

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

inlineinherited

Definition at line 50 of file Converter.h.

{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 320 of file Converter.h.

m_beamY

double Converter::m_beamY

protectedinherited

Definition at line 320 of file Converter.h.

m_beamZ

double Converter::m_beamZ

protectedinherited

Definition at line 320 of file Converter.h.

m_detStore

ServiceHandle<StoreGateSvc> DetDescrConverter::m_detStore

privateinherited

Definition at line 42 of file DetDescrConverter.h.

m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ 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_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 321 of file Converter.h.

---

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

• CaloTTMgrDetDescrCnv.h
• CaloTTMgrDetDescrCnv.cxx