LArRawScContByteStreamCnv Class Reference

#include <LArRawScContByteStreamCnv.h>

Inheritance diagram for LArRawScContByteStreamCnv:

Collaboration diagram for LArRawScContByteStreamCnv:

Public Member Functions
	LArRawScContByteStreamCnv (ISvcLocator *svcloc)
virtual	~LArRawScContByteStreamCnv ()
virtual StatusCode	initialize ()
virtual StatusCode	createObj (IOpaqueAddress *pAddr, DataObject *&pObj)
virtual StatusCode	createRep (DataObject *pObj, IOpaqueAddress *&pAddr)
virtual long	repSvcType () const
	Storage type and class ID. More...
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)

Static Public Member Functions
static long	storageType ()
static const 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 Attributes
MsgStream *	m_log
LArLATOMEDecoder *	m_scTool = nullptr
ByteStreamCnvSvc *	m_ByteStreamEventAccess
IROBDataProviderSvc *	m_rdpSvc
StoreGateSvc *	m_storeGate

Detailed Description

Definition at line 32 of file LArRawScContByteStreamCnv.h.

Constructor & Destructor Documentation

LArRawScContByteStreamCnv()

LArRawScContByteStreamCnv::LArRawScContByteStreamCnv

(

ISvcLocator *

svcloc

)

Definition at line 37 of file LArRawScContByteStreamCnv.cxx.

                                                                        :
  Converter(storageType(), classID(),svcloc), m_log(NULL),m_ByteStreamEventAccess(NULL),m_rdpSvc(NULL),m_storeGate(NULL){}

~LArRawScContByteStreamCnv()

LArRawScContByteStreamCnv::~LArRawScContByteStreamCnv ( )

virtual

Definition at line 40 of file LArRawScContByteStreamCnv.cxx.

                                                      {
 if (m_log) 
    delete m_log;
}

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()

const CLID & LArRawScContByteStreamCnv::classID ( )

static

Definition at line 45 of file LArRawScContByteStreamCnv.cxx.

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

clear()

void Converter::clear ( )

inlineinherited

Definition at line 53 of file Converter.h.

{m_tracks.clear();} 

createObj()

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

virtual

Definition at line 110 of file LArRawScContByteStreamCnv.cxx.

{//Convert Digits from ByteStream to StoreGate
  (*m_log) << MSG::VERBOSE << "Executing CreateObj method for  LArLATOMEHeaderContainer" << endmsg;
 
  if (!m_rdpSvc)
    {(*m_log) << MSG::ERROR << " ROBDataProviderSvc not loaded. Can't read ByteStream." << endmsg;
     return StatusCode::FAILURE;
    }
  ByteStreamAddress *pRE_Addr;
  pRE_Addr = dynamic_cast<ByteStreamAddress*>(pAddr); //Cast from OpaqueAddress to ByteStreamAddress
  if (!pRE_Addr)
    {(*m_log) << MSG::ERROR << "dynamic_cast of OpaqueAdress to ByteStreamAddress failed!" << endmsg;
     return StatusCode::FAILURE;
    }
 
  // patch for the HLT usage
  std::vector<const OFFLINE_FRAGMENTS_NAMESPACE::ROBFragment*> robFrags;
  m_rdpSvc->getROBData( LArByteStream::s_allROBIDs_LATOME, robFrags );
  const RawEvent* re = m_rdpSvc->getEvent();
  if (!re)
    {(*m_log) << MSG::ERROR << "Could not get raw event from ByteStreamInputSvc" << endmsg;
     return StatusCode::FAILURE;
    }
  const std::string& key = *(pAddr->par()); // Get key used in the StoreGateSvc::retrieve function
  // get gain and pass to convert function.
  LArDigitContainer* adc_coll=0;
  LArDigitContainer* adc_bas_coll=0;
  LArRawSCContainer* et_coll=0;
  LArRawSCContainer* et_id_coll=0;
  LArLATOMEHeaderContainer* header_coll=0;
 
  if(key=="SC_ET"){
    et_coll=new LArRawSCContainer();
  }
  else if(key=="SC_ET_ID"){
    et_id_coll=new LArRawSCContainer();
  }
    
  (*m_log) << MSG::DEBUG << "Converting LArRawSCContainer (from ByteStream). key=" << key << endmsg; 
  StatusCode sc;
 
  //Supercell readout
  sc=m_scTool->convert(robFrags,nullptr,adc_coll, adc_bas_coll, et_coll, et_id_coll, header_coll);
  if (sc!=StatusCode::SUCCESS)
    (*m_log) << MSG::WARNING << "Conversion tool returned an error. LAr SC containers might be empty." << endmsg;
  
   if(key=="SC_ET"){
    pObj = SG::asStorable(et_coll);
  }
  else if(key=="SC_ET_ID"){
    pObj = SG::asStorable(et_id_coll);
  }
 
  
  return StatusCode::SUCCESS;
 
}

createRep()

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

virtual

Definition at line 169 of file LArRawScContByteStreamCnv.cxx.

{
   return StatusCode::SUCCESS;
}

initialize()

StatusCode LArRawScContByteStreamCnv::initialize ( )

virtual

Definition at line 51 of file LArRawScContByteStreamCnv.cxx.

{
   StatusCode sc = Converter::initialize(); 
   if(StatusCode::SUCCESS!=sc) 
   { 
    return sc; 
   } 
   m_log = new MsgStream(msgSvc(),"LArLATOMEHeaderBSCnv");
   (*m_log) << MSG::DEBUG<< " initialize " <<endmsg; 
   IService* svc;
  //Get ByteStreamCnvSvc
  if(StatusCode::SUCCESS != serviceLocator()->getService("ByteStreamCnvSvc",svc)){
    (*m_log) << MSG::ERROR << " Can't get ByteStreamEventAccess interface " << endmsg;
    return StatusCode::FAILURE;
  }
  m_ByteStreamEventAccess=dynamic_cast<ByteStreamCnvSvc*>(svc);
  if (m_ByteStreamEventAccess==NULL)
    {
      (*m_log) <<MSG::ERROR<< "  LArRawScContByteStreamCnv: Can't cast to  ByteStreamCnvSvc " <<endmsg; 
      return StatusCode::FAILURE ;
    }
 
  //Get ByteStreamInputSvc (only necessary for reading of digits, not for writing and for channels)
  
  if(StatusCode::SUCCESS != serviceLocator()->getService("ROBDataProviderSvc",svc)){
    (*m_log) << MSG::WARNING << " Can't get ByteStreamInputSvc interface Reading of ByteStream Data not possible. " << endmsg;
    m_rdpSvc=0;
  }
  else {
    m_rdpSvc=dynamic_cast<IROBDataProviderSvc*>(svc);
    if(m_rdpSvc == 0 ) {
      (*m_log) <<MSG::ERROR<< "  LArRawScContByteStreamCnv: Can't cast to  ByteStreamInputSvc " <<endmsg; 
      return StatusCode::FAILURE;
    }
  }
 
  // retrieve ToolSvc
  IToolSvc* toolSvc;
 
  if(StatusCode::SUCCESS != service("ToolSvc",toolSvc)){
    (*m_log) << MSG::ERROR << " Can't get ToolSvc " << endmsg;
    return StatusCode::FAILURE;
  }
  
  //Get SuperCellTool
  sc=toolSvc->retrieveTool("LArLATOMEDecoder",m_scTool);
  if (sc.isFailure()) {
    (*m_log) << MSG::ERROR << "Can't get LArLATOMEDecoder Tool" << endmsg;
    return sc;
  } else {
     (*m_log) << MSG::INFO << "Retrieved LArLATOMEDecoder Tool" << endmsg;
  }
 
  //Get StoreGateSvc
  return service("StoreGateSvc", m_storeGate); 
 
}

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;
      
    }
  }

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()

virtual long LArRawScContByteStreamCnv::repSvcType ( ) const

inlinevirtual

Storage type and class ID.

Definition at line 43 of file LArRawScContByteStreamCnv.h.

{ return i_repSvcType(); }

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++; 
    }
  }

storageType()

static long LArRawScContByteStreamCnv::storageType ( )

inlinestatic

Definition at line 44 of file LArRawScContByteStreamCnv.h.

{ return ByteStreamAddress::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

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_ByteStreamEventAccess

ByteStreamCnvSvc* LArRawScContByteStreamCnv::m_ByteStreamEventAccess

private

Definition at line 51 of file LArRawScContByteStreamCnv.h.

m_log

MsgStream* LArRawScContByteStreamCnv::m_log

private

Definition at line 49 of file LArRawScContByteStreamCnv.h.

m_rdpSvc

IROBDataProviderSvc* LArRawScContByteStreamCnv::m_rdpSvc

private

Definition at line 52 of file LArRawScContByteStreamCnv.h.

m_scTool

LArLATOMEDecoder* LArRawScContByteStreamCnv::m_scTool = nullptr

private

Definition at line 50 of file LArRawScContByteStreamCnv.h.

m_storeGate

StoreGateSvc* LArRawScContByteStreamCnv::m_storeGate

private

Definition at line 53 of file LArRawScContByteStreamCnv.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:

• LArRawScContByteStreamCnv.h
• LArRawScContByteStreamCnv.cxx