McParticleCnv Class Reference

Concrete converter for the McParticle collection. More...

#include <McParticleCnv.h>

Inheritance diagram for McParticleCnv:

Public Member Functions
virtual	~McParticleCnv ()
	McParticleCnv (ISvcLocator *svc)
virtual StatusCode	DataObjectToTObject (DataObject *obj, RootAddress *addr)
	transformation to root
virtual StatusCode	TObjectToDataObject (DataObject *&obj)
	transformation from root
virtual	~McParticleCnv ()
	McParticleCnv (ISvcLocator *svc)
virtual StatusCode	DataObjectToTObject (DataObject *obj, RootAddress *addr)
	transformation to root
virtual StatusCode	TObjectToDataObject (DataObject *&obj)
	transformation from root
virtual	~McParticleCnv ()
	McParticleCnv (ISvcLocator *svc)
virtual StatusCode	DataObjectToTObject (DataObject *obj, RootAddress *addr)
	transformation to root
virtual StatusCode	TObjectToDataObject (DataObject *&obj)
	transformation from root
Public Member Functions inherited from RootEventBaseCnv
	RootEventBaseCnv (const CLID &clid, ISvcLocator *svc)
virtual	~RootEventBaseCnv ()
virtual long	repSvcType () const
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	declareObject (const std::string &fullPath, const CLID &clid, const std::string &treename, const std::string &branchname)
	Store TDS path to link a particular converter to an object on the TDS.
virtual StatusCode	createObj (IOpaqueAddress *addr, DataObject *&dat)
	Convert the persistent object to transient.
virtual StatusCode	createRep (DataObject *pObject, IOpaqueAddress *&refpAddress)
	Convert the transient object to the requested representation.
virtual StatusCode	fillRepRefs (IOpaqueAddress *pAddress, DataObject *pObject)
	Resolve the references of the converted object.
virtual StatusCode	fillObjRefs (IOpaqueAddress *pAddress, DataObject *pObject)
	Resolve the references of the converted object.
TObject *	getReadObject () const
	get the object to be read
	RootEventBaseCnv (const CLID &clid, ISvcLocator *svc)
virtual	~RootEventBaseCnv ()
virtual long	repSvcType () const
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	declareObject (const std::string &fullPath, const CLID &clid, const std::string &treename, const std::string &branchname)
	Store TDS path to link a particular converter to an object on the TDS.
virtual StatusCode	createObj (IOpaqueAddress *addr, DataObject *&dat)
	Convert the persistent object to transient.
virtual StatusCode	createRep (DataObject *pObject, IOpaqueAddress *&refpAddress)
	Convert the transient object to the requested representation.
virtual StatusCode	fillRepRefs (IOpaqueAddress *pAddress, DataObject *pObject)
	Resolve the references of the converted object.
virtual StatusCode	fillObjRefs (IOpaqueAddress *pAddress, DataObject *pObject)
	Resolve the references of the converted object.
TObject *	getReadObject () const
	get the object to be read
	RootEventBaseCnv (const CLID &clid, ISvcLocator *svc)
virtual	~RootEventBaseCnv ()
virtual long	repSvcType () const
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	declareObject (const std::string &fullPath, const CLID &clid, const std::string &treename, const std::string &branchname)
	Store TDS path to link a particular converter to an object on the TDS.
virtual StatusCode	createObj (IOpaqueAddress *addr, DataObject *&dat)
	Convert the persistent object to transient.
virtual StatusCode	createRep (DataObject *pObject, IOpaqueAddress *&refpAddress)
	Convert the transient object to the requested representation.
virtual StatusCode	fillRepRefs (IOpaqueAddress *pAddress, DataObject *pObject)
	Resolve the references of the converted object.
virtual StatusCode	fillObjRefs (IOpaqueAddress *pAddress, DataObject *pObject)
	Resolve the references of the converted object.
TObject *	getReadObject () const
	get the object to be read
Public Member Functions inherited from Converter< Ty1, Ty2 >
destination *	operator (const source &) const
destination *	operator (const source &) const
destination *	operator (const source &) const

Static Public Member Functions
static const CLID &	classID ()
static const CLID &	classID ()
static const CLID &	classID ()
Static Public Member Functions inherited from RootEventBaseCnv
static const unsigned char	storageType ()
	Storage type and class ID.
static const unsigned char	storageType ()
	Storage type and class ID.
static const unsigned char	storageType ()
	Storage type and class ID.

Friends
class	CnvFactory< McParticleCnv >

Additional Inherited Members
Public Types inherited from Converter< Ty1, Ty2 >
typedef Ty1	source
typedef Ty2	destination
typedef Ty1	source
typedef Ty2	destination
typedef Ty1	source
typedef Ty2	destination
Protected Member Functions inherited from Converter< Ty1, Ty2 >
virtual destination *	convert (const source &) const =0
virtual destination *	convert (const source &) const =0
virtual destination *	convert (const source &) const =0
Protected Attributes inherited from RootEventBaseCnv
IRootCnvSvc *	m_cnvSvc
std::vector< RootCnvSvc::Leaf >	m_leaves
RootInterface *	m_rootInterface
	pointer to the RootInterface
IDataProviderSvc *	m_eds
	pointer to eventdataservice
int	m_branchNr
	the branchNr of this converter for writing
int	m_branchNrDst
int	m_branchNrMc
int	m_branchNrRecon
int	m_branchNrEvtRec
int	m_branchNrEvtHeader
int	m_branchNrEvtNavigator
TObject *	m_objRead
	the object that was read
CLID	CLID_top
	the CLID of the upper converter if any
TArrayS *	m_branchNumbers
	array with number of branches for reading
std::string	m_rootBranchname
	root branchname (may be concatenated of severals)
std::string	m_rootTreename
	each converter knows it's treename
std::string	m_currentFileName
std::vector< void * >	m_adresses
	each converter knows the corresponding adresses
RootEvtSelector *	m_evtsel

Detailed Description

Concrete converter for the McParticle collection.

Definition at line 19 of file Event/RootCnvSvc/include/RootCnvSvc/Mc/McParticleCnv.h.

Constructor & Destructor Documentation

~McParticleCnv() [1/3]

virtual McParticleCnv::~McParticleCnv ( )

inlinevirtual

Definition at line 26 of file Event/RootCnvSvc/include/RootCnvSvc/Mc/McParticleCnv.h.

{};

McParticleCnv() [1/3]

McParticleCnv::McParticleCnv

(

ISvcLocator *

svc

)

Definition at line 24 of file McParticleCnv.cxx.

                                               : RootEventBaseCnv( classID(), svc ) {
  // Here we associate this converter with the /Event path on the TDS.
  MsgStream log( msgSvc(), "McParticleCnv" );
  // log << MSG::DEBUG << "Constructor called for " << objType() << endmsg;
  m_rootBranchname = "m_mcParticleCol";
  // declareObject(EventModel::MC::McParticleCol, objType(), m_rootTreename, m_rootBranchname);
  m_adresses.push_back( &m_mcParticleCol );
  m_mcParticleCol = 0;
}

~McParticleCnv() [2/3]

virtual McParticleCnv::~McParticleCnv ( )

inlinevirtual

Definition at line 26 of file InstallArea/x86_64-el9-gcc13-dbg/include/RootCnvSvc/Mc/McParticleCnv.h.

{};

McParticleCnv() [2/3]

McParticleCnv::McParticleCnv

(

ISvcLocator *

svc

)

~McParticleCnv() [3/3]

virtual McParticleCnv::~McParticleCnv ( )

inlinevirtual

Definition at line 26 of file InstallArea/x86_64-el9-gcc13-opt/include/RootCnvSvc/Mc/McParticleCnv.h.

{};

McParticleCnv() [3/3]

McParticleCnv::McParticleCnv

(

ISvcLocator *

svc

)

Member Function Documentation

classID() [1/3]

const CLID & McParticleCnv::classID ( )

inlinestatic

Definition at line 24 of file Event/RootCnvSvc/include/RootCnvSvc/Mc/McParticleCnv.h.

{ return CLID_McParticleCol; }

Referenced by RootCnvSvc::createConverter(), and McParticleCnv().

classID() [2/3]

const CLID & McParticleCnv::classID ( )

inlinestatic

Definition at line 24 of file InstallArea/x86_64-el9-gcc13-dbg/include/RootCnvSvc/Mc/McParticleCnv.h.

{ return CLID_McParticleCol; }

classID() [3/3]

const CLID & McParticleCnv::classID ( )

inlinestatic

Definition at line 24 of file InstallArea/x86_64-el9-gcc13-opt/include/RootCnvSvc/Mc/McParticleCnv.h.

{ return CLID_McParticleCol; }

DataObjectToTObject() [1/3]

StatusCode McParticleCnv::DataObjectToTObject ( DataObject * obj, RootAddress * addr )

virtual

transformation to root

Implements RootEventBaseCnv.

Definition at line 141 of file McParticleCnv.cxx.

                                                                                      {
 
  MsgStream log( msgSvc(), "McParticleCnv" );
  log << MSG::DEBUG << "McParticleCnv::DataObjectToTObject" << endmsg;
  StatusCode sc = StatusCode::SUCCESS;
 
  McParticleCol* mcParticleCnvTds = dynamic_cast<McParticleCol*>( obj );
  if ( !mcParticleCnvTds )
  {
    log << MSG::ERROR << "Could not downcast to McParticleCol" << endmsg;
    return StatusCode::FAILURE;
  }
 
  DataObject* evt;
  m_eds->findObject( EventModel::MC::Event, evt );
  if ( evt == NULL )
  {
    log << MSG::ERROR << "Could not get McEvent in TDS " << endmsg;
    return StatusCode::FAILURE;
  }
  McEvent* devtTds = dynamic_cast<McEvent*>( evt );
  if ( !devtTds )
  { log << MSG::ERROR << "McParticleCnv:Could not downcast to TDS McEvent" << endmsg; }
  IOpaqueAddress* addr;
 
  m_cnvSvc->getMcCnv()->createRep( evt, addr );
  TMcEvent* McEvt = m_cnvSvc->getMcCnv()->getWriteObject();
 
  const TObjArray* m_mcParticleCol = McEvt->getMcParticleCol();
  if ( !m_mcParticleCol ) { return sc; }
  McEvt->clearMcParticleCol(); // necessary in case there is I/O at the same time since array
                               // is static
 
  McParticleCol::const_iterator mcParticleTds;
 
  for ( mcParticleTds = mcParticleCnvTds->begin(); mcParticleTds != mcParticleCnvTds->end();
        mcParticleTds++ )
  {
    Int_t particleID = ( *mcParticleTds )->particleProperty();
    Int_t trackIndex = ( *mcParticleTds )->trackIndex();
 
    Int_t vertexIndex0 = ( *mcParticleTds )->vertexIndex0();
    Int_t vertexIndex1 = ( *mcParticleTds )->vertexIndex1();
    Int_t statusFlags  = ( *mcParticleTds )->statusFlags();
 
    HepLorentzVector initialPosition     = ( *mcParticleTds )->initialPosition();
    HepLorentzVector finalPosition       = ( *mcParticleTds )->finalPosition();
    HepLorentzVector initialFourMomentum = ( *mcParticleTds )->initialFourMomentum();
    // HepLorentzVector finalFourMomentum   = (*mcParticleTds)->finalFourMomentum();
 
    // cout << " HepID " << (*mcParticleTds)->particleProperty() << endl;
    // cout << " init pos " << initialPosition.x() << "  " << initialPosition.y() << "  " <<
    // initialPosition.z() << endl;
    Int_t mother = -99;
    if ( !( *mcParticleTds )->primaryParticle() )
      mother = ( ( *mcParticleTds )->mother() ).trackIndex();
 
    vector<Int_t>              daughters;
    SmartRefVector<McParticle> daughterList = ( *mcParticleTds )->daughterList();
    // SmartRefVector<McParticle>::iterator iter;
    for ( int iPar = 0; iPar < daughterList.size(); iPar++ )
    {
      // cout <<"daughter Index   " <<daughterList[iPar]->getTrackIndex()<<endl;
      daughters.push_back( daughterList[iPar]->trackIndex() );
    }
    /*
    for (int idau = 0; idau<daughters.size();idau++){
      cout <<"daughter Index      " <<daughters[idau]<<endl;
    }
 
    cout<<"###############################"<<endl;
    */
    TMcParticle* mcParticleRoot = new TMcParticle();
    // m_common.m_mcParticleMap[(*mcParticleTds)] = mcParticleRoot;
 
    mcParticleRoot->setParticleID( particleID );
    mcParticleRoot->setTrackIndex( trackIndex );
 
    mcParticleRoot->setVertexIndex0( vertexIndex0 );
    mcParticleRoot->setVertexIndex1( vertexIndex1 );
 
    mcParticleRoot->setStatusFlags( statusFlags );
 
    mcParticleRoot->setInitialPositionX( initialPosition.x() );
    mcParticleRoot->setInitialPositionY( initialPosition.y() );
    mcParticleRoot->setInitialPositionZ( initialPosition.z() );
    mcParticleRoot->setInitialPositionT( initialPosition.t() );
 
    mcParticleRoot->setFinalPositionX( finalPosition.x() );
    mcParticleRoot->setFinalPositionY( finalPosition.y() );
    mcParticleRoot->setFinalPositionZ( finalPosition.z() );
    mcParticleRoot->setFinalPositionT( finalPosition.t() );
 
    mcParticleRoot->setInitialMomentumX( initialFourMomentum.x() );
    mcParticleRoot->setInitialMomentumY( initialFourMomentum.y() );
    mcParticleRoot->setInitialMomentumZ( initialFourMomentum.z() );
    mcParticleRoot->setInitialMomentumE( initialFourMomentum.e() );
 
    // mcParticleRoot->setFinalMomentumX(finalFourMomentum.x());
    // mcParticleRoot->setFinalMomentumY(finalFourMomentum.y());
    // mcParticleRoot->setFinalMomentumZ(finalFourMomentum.z());
    // mcParticleRoot->setFinalMomentumE(finalFourMomentum.e());
 
    mcParticleRoot->setMother( mother );
    mcParticleRoot->setDaughters( daughters );
    /*
    vector<int> dau = mcParticleRoot->getDaughters();
 
    for (int idau = 0; idau<daughters.size();idau++){
      cout <<"daughter Index root      " <<dau[idau]<<endl;
    }
    */
 
    McEvt->addMcParticle( mcParticleRoot );
  }
  return StatusCode::SUCCESS;
}

DataObjectToTObject() [2/3]

virtual StatusCode McParticleCnv::DataObjectToTObject ( DataObject * obj, RootAddress * addr )

virtual

transformation to root

Implements RootEventBaseCnv.

DataObjectToTObject() [3/3]

virtual StatusCode McParticleCnv::DataObjectToTObject ( DataObject * obj, RootAddress * addr )

virtual

transformation to root

Implements RootEventBaseCnv.

TObjectToDataObject() [1/3]

StatusCode McParticleCnv::TObjectToDataObject ( DataObject *& obj )

virtual

transformation from root

Implements RootEventBaseCnv.

Definition at line 34 of file McParticleCnv.cxx.

                                                                       {
  // creation of TDS object from root object
 
  MsgStream log( msgSvc(), "McParticleCnv" );
  log << MSG::DEBUG << "McParticleCnv::TObjectToDataObject" << endmsg;
  StatusCode sc = StatusCode::SUCCESS;
 
  // create the TDS location for the MdcMc Collection
  McParticleCol* mcParticleTdsCol = new McParticleCol;
  refpObject                      = mcParticleTdsCol;
 
  // now convert
  if ( !m_mcParticleCol ) return sc;
  TIter        mcParticleIter( m_mcParticleCol );
  TMcParticle* mcParticleRoot = 0;
 
  vector<int> mothers;
  // vector<vector<int> > daughtList;
 
  while ( ( mcParticleRoot = (TMcParticle*)mcParticleIter.Next() ) )
  {
    unsigned int particleID = mcParticleRoot->getParticleID();
    unsigned int trackIndex = mcParticleRoot->getTrackIndex();
 
    int          vertexIndex0 = mcParticleRoot->getVertexIndex0();
    int          vertexIndex1 = mcParticleRoot->getVertexIndex1();
    unsigned int statusFlags  = mcParticleRoot->getStatusFlags();
 
    double xInitialPosition = mcParticleRoot->getInitialPositionX();
    double yInitialPosition = mcParticleRoot->getInitialPositionY();
    double zInitialPosition = mcParticleRoot->getInitialPositionZ();
    double tInitialPosition = mcParticleRoot->getInitialPositionT();
 
    double xFinalPosition = mcParticleRoot->getFinalPositionX();
    double yFinalPosition = mcParticleRoot->getFinalPositionY();
    double zFinalPosition = mcParticleRoot->getFinalPositionZ();
    double tFinalPosition = mcParticleRoot->getFinalPositionT();
 
    double xInitialMomentum = mcParticleRoot->getInitialMomentumX();
    double yInitialMomentum = mcParticleRoot->getInitialMomentumY();
    double zInitialMomentum = mcParticleRoot->getInitialMomentumZ();
    double eInitialMomentum = mcParticleRoot->getInitialMomentumE();
 
    // double xFinalMomentum =   mcParticleRoot->getFinalMomentumX();
    // double yFinalMomentum =   mcParticleRoot->getFinalMomentumY();
    // double zFinalMomentum =   mcParticleRoot->getFinalMomentumZ();
    // double eFinalMomentum =   mcParticleRoot->getFinalMomentumE();
 
    HepLorentzVector initialMomentum( xInitialMomentum, yInitialMomentum, zInitialMomentum,
                                      eInitialMomentum );
    HepLorentzVector initialPosition( xInitialPosition, yInitialPosition, zInitialPosition,
                                      tInitialPosition );
    // HepLorentzVector finalMomentum(xFinalMomentum, yFinalMomentum, zFinalMomentum,
    // eFinalMomentum);
    HepLorentzVector finalPosition( xFinalPosition, yFinalPosition, zFinalPosition,
                                    tFinalPosition );
 
    int         mother    = mcParticleRoot->getMother();
    vector<int> daughters = mcParticleRoot->getDaughters();
 
    mothers.push_back( mother );
    // daughtList.push_back(daughters);
 
    McParticle* mcParticleTds                    = new McParticle;
    m_common.m_rootMcParticleMap[mcParticleRoot] = mcParticleTds;
 
    mcParticleTds->initialize( particleID, statusFlags, initialMomentum, initialPosition );
    mcParticleTds->setTrackIndex( trackIndex );
    mcParticleTds->addStatusFlag( statusFlags );
 
    mcParticleTds->setVertexIndex0( vertexIndex0 );
    mcParticleTds->setVertexIndex1( vertexIndex1 );
    // mcParticleTds->finalize(finalMomentum, finalPosition);
    mcParticleTds->finalize( finalPosition );
 
    mcParticleTdsCol->push_back( mcParticleTds );
  }
 
  // Set Mother and DaughterList
  McParticleCol::iterator iter;
  int                     i = 0;
  for ( iter = mcParticleTdsCol->begin(); iter != mcParticleTdsCol->end(); iter++, i++ )
  {
    // cout<<"    *****   mothers[ "<<i<<"] = "<<mothers[i]<<endl;
    if ( mothers[i] != -99 )
    {
      McParticleCol::iterator mcParticleTds;
      for ( mcParticleTds = mcParticleTdsCol->begin();
            mcParticleTds != mcParticleTdsCol->end(); mcParticleTds++ )
      {
        int trackIndex = ( *mcParticleTds )->trackIndex();
        if ( trackIndex == mothers[i] )
        {
          ( *iter )->setMother( *mcParticleTds );
          ( *mcParticleTds )->addDaughter( *iter );
          break;
        }
      }
    }
    else { ( *iter )->setMother( *iter ); }
  }
  // m_mcParticleCol->Delete();  // wensp add 2005/12/30
  delete m_mcParticleCol;
  m_mcParticleCol = 0;
  return StatusCode::SUCCESS;
}

TObjectToDataObject() [2/3]

virtual StatusCode McParticleCnv::TObjectToDataObject ( DataObject *& obj )

virtual

transformation from root

Implements RootEventBaseCnv.

TObjectToDataObject() [3/3]

virtual StatusCode McParticleCnv::TObjectToDataObject ( DataObject *& obj )

virtual

transformation from root

Implements RootEventBaseCnv.

Friends And Related Symbol Documentation

CnvFactory< McParticleCnv >

friend class CnvFactory< McParticleCnv >

friend

Definition at line 9 of file Event/RootCnvSvc/include/RootCnvSvc/Mc/McParticleCnv.h.

---

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

• Event/RootCnvSvc/include/RootCnvSvc/Mc/McParticleCnv.h
• InstallArea/x86_64-el9-gcc13-dbg/include/RootCnvSvc/Mc/McParticleCnv.h
• InstallArea/x86_64-el9-gcc13-opt/include/RootCnvSvc/Mc/McParticleCnv.h
• McParticleCnv.cxx