McParticleCnv Class Reference

Concrete converter for the McParticle collection. More...

#include <McParticleCnv.h>

Inheritance diagram for McParticleCnv:

Public Member Functions
virtual	~McParticleCnv ()
virtual	~McParticleCnv ()
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.
virtual StatusCode	DataObjectToTObject (DataObject *dat, RootAddress *addr)=0
	Do the concrete conversion from TDS to ROOT.
virtual StatusCode	TObjectToDataObject (DataObject *&dat)=0
	Do the concrete conversion from ROOT to TDS.
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.
virtual StatusCode	DataObjectToTObject (DataObject *dat, RootAddress *addr)=0
	Do the concrete conversion from TDS to ROOT.
virtual StatusCode	TObjectToDataObject (DataObject *&dat)=0
	Do the concrete conversion from ROOT to TDS.
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

Static Public Member Functions
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.

Protected Member Functions
	McParticleCnv (ISvcLocator *svc)
virtual StatusCode	DataObjectToTObject (DataObject *obj, RootAddress *addr)
	transformation to root
virtual StatusCode	TObjectToDataObject (DataObject *&obj)
	transformation from root
	McParticleCnv (ISvcLocator *svc)
virtual StatusCode	DataObjectToTObject (DataObject *obj, RootAddress *addr)
	transformation to root
virtual StatusCode	TObjectToDataObject (DataObject *&obj)
	transformation from root
Protected Member Functions inherited from Converter< Ty1, Ty2 >
virtual destination *	convert (const source &) const =0
virtual destination *	convert (const source &) const =0

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
Protected Attributes inherited from RootEventBaseCnv
RootCnvSvc *	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/RootCnvSvc-02-01-11/RootCnvSvc/Mc/McParticleCnv.h.

Constructor & Destructor Documentation

~McParticleCnv() [1/2]

virtual McParticleCnv::~McParticleCnv ( )

inlinevirtual

Definition at line 29 of file Event/RootCnvSvc/RootCnvSvc-02-01-11/RootCnvSvc/Mc/McParticleCnv.h.

{ };

McParticleCnv() [1/2]

McParticleCnv::McParticleCnv ( ISvcLocator *  svc )

protected

Definition at line 25 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() << endreq;
    m_rootBranchname ="m_mcParticleCol";
    //declareObject(EventModel::MC::McParticleCol, objType(), m_rootTreename, m_rootBranchname);
    m_adresses.push_back(&m_mcParticleCol);
    m_mcParticleCol=0;
}

~McParticleCnv() [2/2]

virtual McParticleCnv::~McParticleCnv ( )

inlinevirtual

Definition at line 29 of file InstallArea/include/RootCnvSvc/RootCnvSvc/Mc/McParticleCnv.h.

{ };

McParticleCnv() [2/2]

McParticleCnv::McParticleCnv ( ISvcLocator *  svc )

protected

Member Function Documentation

classID() [1/2]

static const CLID & McParticleCnv::classID ( )

inlinestatic

Definition at line 24 of file Event/RootCnvSvc/RootCnvSvc-02-01-11/RootCnvSvc/Mc/McParticleCnv.h.

  {
    return CLID_McParticleCol; 
  }

classID() [2/2]

static const CLID & McParticleCnv::classID ( )

inlinestatic

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

  {
    return CLID_McParticleCol; 
  }

DataObjectToTObject() [1/2]

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

protectedvirtual

transformation to root

Implements RootEventBaseCnv.

Definition at line 139 of file McParticleCnv.cxx.

                                                                                   {
 
  MsgStream log(msgSvc(), "McParticleCnv");
  log << MSG::DEBUG << "McParticleCnv::DataObjectToTObject" << endreq;
  StatusCode sc=StatusCode::SUCCESS;
 
  McParticleCol *mcParticleCnvTds=dynamic_cast<McParticleCol *> (obj);
  if (!mcParticleCnvTds) {
    log << MSG::ERROR << "Could not downcast to McParticleCol" << endreq;
    return StatusCode::FAILURE;
  }
 
  DataObject *evt;
  m_eds->findObject(EventModel::MC::Event,evt);
  if (evt==NULL) {
    log << MSG::ERROR << "Could not get McEvent in TDS "  << endreq;
    return StatusCode::FAILURE;
  }
  McEvent * devtTds=dynamic_cast<McEvent *> (evt);
  if (!devtTds) {
    log << MSG::ERROR << "McParticleCnv:Could not downcast to TDS McEvent" << endreq;
  }
  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/2]

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

protectedvirtual

transformation to root

Implements RootEventBaseCnv.

TObjectToDataObject() [1/2]

StatusCode McParticleCnv::TObjectToDataObject ( DataObject *&  obj )

protectedvirtual

transformation from root

Implements RootEventBaseCnv.

Definition at line 37 of file McParticleCnv.cxx.

                                                                     {
  // creation of TDS object from root object
 
    MsgStream log(msgSvc(), "McParticleCnv");
    log << MSG::DEBUG << "McParticleCnv::TObjectToDataObject" << endreq;
    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/2]

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

protectedvirtual

transformation from root

Implements RootEventBaseCnv.

Friends And Related Function Documentation

CnvFactory< McParticleCnv >

friend class CnvFactory< McParticleCnv >

friend

Definition at line 9 of file Event/RootCnvSvc/RootCnvSvc-02-01-11/RootCnvSvc/Mc/McParticleCnv.h.

---

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

• Event/RootCnvSvc/RootCnvSvc-02-01-11/RootCnvSvc/Mc/McParticleCnv.h
• InstallArea/include/RootCnvSvc/RootCnvSvc/Mc/McParticleCnv.h
• McParticleCnv.cxx