McParticleCnv Class Reference

Concrete converter for the McParticle collection. More...

#include <McParticleCnv.h>

Inheritance diagram for McParticleCnv:

Public Member Functions
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
Public Member Functions inherited from Converter< Ty1, Ty2 >
destination *	operator (const source &) const

Static Public Member Functions
static const CLID &	classID ()
Static Public Member Functions inherited from RootEventBaseCnv
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
Protected Member Functions inherited from Converter< Ty1, Ty2 >
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
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 McParticleCnv.h.

Constructor & Destructor Documentation

~McParticleCnv()

virtual McParticleCnv::~McParticleCnv ( )

inlinevirtual

Definition at line 29 of file McParticleCnv.h.

{ };

McParticleCnv()

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

Member Function Documentation

classID()

static const CLID & McParticleCnv::classID ( )

inlinestatic

Definition at line 24 of file McParticleCnv.h.

  {
    return CLID_McParticleCol; 
  }

DataObjectToTObject()

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

TObjectToDataObject()

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

Friends And Related Function Documentation

CnvFactory< McParticleCnv >

friend class CnvFactory< McParticleCnv >

friend

Definition at line 9 of file McParticleCnv.h.

---

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

• McParticleCnv.h
• McParticleCnv.cxx