Main Page | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members | Related Pages

McParticleCnv Class Reference

Concrete converter for the McParticle collection. More...

#include <McParticleCnv.h>

Inheritance diagram for McParticleCnv:

RootEventBaseCnv RootEventBaseCnv Converter< Ty1, Ty2 > Converter< Ty1, Ty2 > Converter< Ty1, Ty2 > Converter< Ty1, Ty2 > List of all members.

Public Types

typedef Ty2 destination
typedef Ty2 destination
typedef Ty1 source
typedef Ty1 source

Public Member Functions

virtual StatusCode createObj (IOpaqueAddress *addr, DataObject *&dat)
 Convert the persistent object to transient.
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 createRep (DataObject *pObject, IOpaqueAddress *&refpAddress)
 Convert the transient object to the requested representation.
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.
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 fillObjRefs (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 fillRepRefs (IOpaqueAddress *pAddress, DataObject *pObject)
 Resolve the references of the converted object.
virtual StatusCode fillRepRefs (IOpaqueAddress *pAddress, DataObject *pObject)
 Resolve the references of the converted object.
virtual StatusCode finalize ()
virtual StatusCode finalize ()
TObject * getReadObject () const
 get the object to be read
TObject * getReadObject () const
 get the object to be read
virtual StatusCode initialize ()
virtual StatusCode initialize ()
destination * operator (const source &) const
destination * operator (const source &) const
virtual long repSvcType () const
virtual long repSvcType () const
virtual ~McParticleCnv ()
virtual ~McParticleCnv ()

Static Public Member Functions

const CLID & classID ()
const CLID & classID ()
const unsigned char storageType ()
 Storage type and class ID.
const unsigned char storageType ()
 Storage type and class ID.

Protected Member Functions

virtual destination * convert (const source &) const =0
virtual destination * convert (const source &) const =0
virtual StatusCode DataObjectToTObject (DataObject *obj, RootAddress *addr)
 transformation to root
virtual StatusCode DataObjectToTObject (DataObject *obj, RootAddress *addr)
 transformation to root
 McParticleCnv (ISvcLocator *svc)
 McParticleCnv (ISvcLocator *svc)
virtual StatusCode TObjectToDataObject (DataObject *&obj)
 transformation from root
virtual StatusCode TObjectToDataObject (DataObject *&obj)
 transformation from root

Protected Attributes

CLID CLID_top
 the CLID of the upper converter if any
std::vector< void * > m_adresses
 each converter knows the corresponding adresses
std::vector< void * > m_adresses
 each converter knows the corresponding adresses
int m_branchNr
 the branchNr of this converter for writing
int m_branchNrDst
int m_branchNrEvtHeader
int m_branchNrEvtRec
int m_branchNrMc
int m_branchNrRecon
TArrayS * m_branchNumbers
 array with number of branches for reading
TArrayS * m_branchNumbers
 array with number of branches for reading
RootCnvSvc * m_cnvSvc
RootCnvSvc * m_cnvSvc
std::string m_currentFileName
IDataProviderSvc * m_eds
 pointer to eventdataservice
IDataProviderSvc * m_eds
 pointer to eventdataservice
RootEvtSelector * m_evtsel
RootEvtSelector * m_evtsel
std::vector< RootCnvSvc::Leaf > m_leaves
std::vector< RootCnvSvc::Leaf > m_leaves
TObject * m_objRead
 the object that was read
TObject * m_objRead
 the object that was read
std::string m_rootBranchname
 root branchname (may be concatenated of severals)
RootInterface * m_rootInterface
 pointer to the RootInterface
RootInterface * m_rootInterface
 pointer to the RootInterface
std::string m_rootTreename
 each converter knows it's treename

Private Attributes

commonData m_common
 relational maps
TObjArray * m_mcParticleCol
 root object to be read
TObjArray * m_mcParticleCol
 root object to be read

Friends

class CnvFactory<McParticleCnv>

Detailed Description

Concrete converter for the McParticle collection.

Member Typedef Documentation

template<class Ty1, class Ty2>
typedef Ty2 Converter< Ty1, Ty2 >::destination [inherited]
 

template<class Ty1, class Ty2>
typedef Ty2 Converter< Ty1, Ty2 >::destination [inherited]
 

template<class Ty1, class Ty2>
typedef Ty1 Converter< Ty1, Ty2 >::source [inherited]
 

template<class Ty1, class Ty2>
typedef Ty1 Converter< Ty1, Ty2 >::source [inherited]
 

Constructor & Destructor Documentation

virtual McParticleCnv::~McParticleCnv (   )  [inline, virtual]
  

00029 { };

McParticleCnv::McParticleCnv (  ISvcLocator *  svc  )  [protected]
  

00026 : RootEventBaseCnv(classID(), svc)
00027 {
00028     // Here we associate this converter with the /Event path on the TDS.
00029     MsgStream log(msgSvc(), "McParticleCnv");
00030     //log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
00031     m_rootBranchname ="m_mcParticleCol";
00032     //declareObject(EventModel::MC::McParticleCol, objType(), m_rootTreename, m_rootBranchname);
00033     m_adresses.push_back(&m_mcParticleCol);
00034     m_mcParticleCol=0;
00035 }

virtual McParticleCnv::~McParticleCnv (   )  [inline, virtual]
  

00029 { };

McParticleCnv::McParticleCnv (  ISvcLocator *  svc  )  [protected]
 

Member Function Documentation

const CLID& McParticleCnv::classID (   )  [inline, static]
  

00025   {
00026     return CLID_McParticleCol; 
00027   }

const CLID& McParticleCnv::classID (   )  [inline, static]
  

00025   {
00026     return CLID_McParticleCol; 
00027   }

template<class Ty1, class Ty2>
virtual destination* Converter< Ty1, Ty2 >::convert (  const source &   )  const [protected, pure virtual, inherited]
 

template<class Ty1, class Ty2>
virtual destination* Converter< Ty1, Ty2 >::convert (  const source &   )  const [protected, pure virtual, inherited]
 

virtual StatusCode RootEventBaseCnv::createObj (  IOpaqueAddress *  addr,
DataObject *&  dat
)  [virtual, inherited]
 

Convert the persistent object to transient.

Reimplemented in EventCnv, and EventCnv.

StatusCode RootEventBaseCnv::createObj (  IOpaqueAddress *  addr,
DataObject *&  dat
)  [virtual, inherited]
 

Convert the persistent object to transient.

Reimplemented in EventCnv, and EventCnv.

00157                                                            {
00158   // transform ROOT object to TDS object
00159   MsgStream log(msgSvc(), "RootEventBaseCnv");
00160   log << MSG::DEBUG << "RootEventBaseCnv::createObj with clid " <<addr->clID()<< endreq;
00161   StatusCode sc;
00162 
00163   // add 2005-11-29
00164   // log<<MSG::INFO<<"######### RootEventBaseCnv ::createObj begin of createObj:  m_branchNumbers "<<m_branchNumbers->GetSize()<<"###############"<<endreq;
00165 
00166   RootAddress *raddr=dynamic_cast<RootAddress *>(addr);
00167   if (!raddr) {
00168     log << MSG::ERROR << "Could not downcast to Root address" << endreq;
00169     return StatusCode::FAILURE;
00170   }
00171  
00172   static int temp =0;   //control the begin of each files 2005-12-01
00173   static int entryN =0; //control the event number of each files 2005-21-01
00174   static int brN =0; //control munber of branch of the tree;
00175   int lastBrn = brN;
00176   //lastBrn = brN; 
00177   static int branchN=0;
00178   static bool isSet=true;
00179 
00180   static int entryBefore = 0;
00181   static bool addEntryEachFile = true;
00182  
00183  
00184    if(m_rootInterface->getENDFILE() || (temp >0 && temp < branchN)){ // if the file has get the end:y the go to next file to create a new tree
00185     
00186     if(m_rootInterface->getENDFILE() ) {
00187       entryN = 0;
00188     }
00189     
00190     temp++;
00191     
00192     delete m_branchNumbers;
00193     m_branchNumbers = new TArrayS(0);
00194    
00195     if(temp == branchN) {
00196       temp =0;
00197     }
00198    }
00199 
00200    if(m_rootInterface->getENDFILE()) addEntryEachFile = true;
00201    // the 2nd method 
00202    if(m_evtsel->getRecId() - entryBefore == 0) { // first event in this file
00203      delete m_branchNumbers;
00204      m_branchNumbers = new TArrayS(0);
00205    }
00206 
00207    //new method to initialize the branchNumber 
00208    if(m_currentFileName=="") m_currentFileName = m_rootInterface->getCurrentFileName();
00209    if(!(m_currentFileName == m_rootInterface->getCurrentFileName())){ 
00210            m_currentFileName = m_rootInterface->getCurrentFileName();
00211            delete m_branchNumbers;
00212            m_branchNumbers = new TArrayS(0);
00213    }   
00214    //---------------------------------------- 
00215 
00216    
00217   if (m_branchNumbers->GetSize()<=0) {
00218     if(isSet) brN++;
00219     int branchNumber;
00220     for (int nb=0;nb<raddr->getNrBranches();nb++) {
00221       sc=m_rootInterface->setBranchAddress(raddr->getTreename().c_str(),raddr->getBranchname(nb).c_str(),m_adresses[nb],branchNumber);
00222       if (!sc.isSuccess()) 
00223       {
00224         if(isSet) brN--;    //liangyt: if fail to retrieve this branch, this will be not a effective branch.
00225         //entryN++;   //liangyt: this is the second method
00226         if(temp>0) temp--;  //temp > 0 means recording effective branch number.
00227         return sc;
00228       }
00229       m_branchNumbers->Set(nb+1);
00230       m_branchNumbers->AddAt(branchNumber,nb);
00231      
00232     }
00233   }
00234 
00236   if(addEntryEachFile&&(m_evtsel->getRecId()>entryBefore)){   // for a new file, add entry for ONLY one time.
00237     entryBefore += m_rootInterface->getEntries();
00238     addEntryEachFile = false;
00239   }
00240 
00241   if(lastBrn == brN && isSet ){
00242     branchN = brN;
00243     isSet=false;    
00244   }
00245 
00246   if(isSet==false) log << MSG::INFO <<" 1st method set event as : "<<int(entryN/branchN)<<endreq;
00247   if(isSet==false) raddr->setEntryNr(int(entryN/branchN));//former method, keep it to be backup.
00248   if(m_evtsel) log << MSG::INFO <<" event id = "<<m_evtsel->getRecId()<<endreq;
00249 
00250   int eventID = 0;
00251   if(entryBefore > m_evtsel->getRecId())
00252     eventID = m_evtsel->getRecId() + m_rootInterface->getEntries() - entryBefore;
00253   else if(entryBefore == m_evtsel->getRecId()) eventID = 0;
00254   else log << MSG::ERROR <<"eventId error!!!"<<endreq;
00255   
00256   log << MSG::INFO <<" 2nd  method set event as : "<<eventID<<endreq;
00257   if(m_evtsel) raddr->setEntryNr(eventID);
00258   // read branch
00259   
00260   if (m_branchNumbers->GetSize()>0) {
00261     int nbtot=0,nb;
00262     for (int ib=0;ib<m_branchNumbers->GetSize();ib++) {
00263       //sc=m_rootInterface->getBranchEntry(m_branchNumbers->At(ib),raddr->getEntryNr(),nb);
00264       //change to get branch entry with addr(set address for each entry) liangyt
00265       sc=m_rootInterface->getBranchEntry(m_branchNumbers->At(ib),raddr->getEntryNr(),m_adresses[ib],nb);      
00266       if (sc.isFailure()) {
00267         log << MSG::ERROR << "Could not read branch " << raddr->getBranchname(nb) << endreq;
00268         return sc;
00269       }
00270       nbtot+=nb;
00271     }
00272   }
00273      
00274   else {   // get ROOT object
00275     if (CLID_top) {
00276       IConverter *p=conversionSvc()->converter(CLID_top);
00277       RootEventBaseCnv *cnv=dynamic_cast<RootEventBaseCnv *>(p);
00278       if (!cnv) {
00279         log << MSG::ERROR << "Could not downcast to RootEventBaseCnv "  << endreq;
00280         return StatusCode::FAILURE;
00281       }
00282       m_objRead=cnv->getReadObject();
00283     }
00284   }
00285 
00286   //do concrete transformation in derived converter
00287   sc = TObjectToDataObject(refpObject);
00288   if (sc.isFailure()) {
00289     log << MSG::ERROR << "Could not transform object" << endreq;
00290     return sc;
00291   }
00292 
00293   // verify if we have to register
00294   IRegistry* ent = addr->registry();
00295   if ( ent == 0)   {
00296        sc=m_eds->registerObject(raddr->getPath(),refpObject);
00297        if (sc.isFailure()) {
00298          log << MSG::ERROR << "Could not register object " << raddr->getPath()<<" status "<<sc.getCode()<<endreq;
00299        }
00300        //    }
00301   }
00302 
00303   entryN++;
00304   return StatusCode::SUCCESS;
00305 }

virtual StatusCode RootEventBaseCnv::createRep (  DataObject *  pObject,
IOpaqueAddress *&  refpAddress
)  [virtual, inherited]
 

Convert the transient object to the requested representation.

StatusCode RootEventBaseCnv::createRep (  DataObject *  pObject,
IOpaqueAddress *&  refpAddress
)  [virtual, inherited]
 

Convert the transient object to the requested representation. 

00078                                                            {
00079   // Purpose and Method: Convert the transient object to ROOT
00080 
00081   MsgStream log(msgSvc(), "RootEventBaseCnv");
00082  
00083   StatusCode sc= StatusCode::SUCCESS;
00084   // get the corresponding address 
00085   RootAddress *rootaddr;
00086   sc=m_cnvSvc->createAddress(obj,addr);
00087 
00088   rootaddr = dynamic_cast<RootAddress *>(addr);
00089 
00090   if (sc.isFailure() || !rootaddr  ) { 
00091     log << MSG::ERROR << "Could not create address for clid " <<obj->clID()<<", objname "<<obj->name()<<endreq;
00092     return StatusCode::FAILURE;
00093   }
00094 
00095   //  do the real conversion in the derived converter
00096   sc = DataObjectToTObject(obj,rootaddr);
00097 
00098   delete addr;
00099   addr = NULL;
00100 
00101   if (sc.isFailure()) {
00102     log << MSG::ERROR << "Could not transform object" << endreq;
00103     return sc;
00104   }
00105 
00106   return StatusCode::SUCCESS;
00107 }

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

transformation to root

Implements RootEventBaseCnv.

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

transformation to root

Implements RootEventBaseCnv.

00139                                                                                    {
00140 
00141   MsgStream log(msgSvc(), "McParticleCnv");
00142   log << MSG::DEBUG << "McParticleCnv::DataObjectToTObject" << endreq;
00143   StatusCode sc=StatusCode::SUCCESS;
00144  
00145   McParticleCol *mcParticleCnvTds=dynamic_cast<McParticleCol *> (obj);
00146   if (!mcParticleCnvTds) {
00147     log << MSG::ERROR << "Could not downcast to McParticleCol" << endreq;
00148     return StatusCode::FAILURE;
00149   }
00150  
00151   DataObject *evt;
00152   m_eds->findObject(EventModel::MC::Event,evt);
00153   if (evt==NULL) {
00154     log << MSG::ERROR << "Could not get McEvent in TDS "  << endreq;
00155     return StatusCode::FAILURE;
00156   }
00157   McEvent * devtTds=dynamic_cast<McEvent *> (evt);
00158   if (!devtTds) {
00159     log << MSG::ERROR << "McParticleCnv:Could not downcast to TDS McEvent" << endreq;
00160   }
00161   IOpaqueAddress *addr;
00162 
00163   m_cnvSvc->getMcCnv()->createRep(evt,addr); 
00164   TMcEvent *McEvt=m_cnvSvc->getMcCnv()->getWriteObject();
00165 
00166   const TObjArray *m_mcParticleCol = McEvt->getMcParticleCol();
00167   if (!m_mcParticleCol) { return sc; }
00168   McEvt->clearMcParticleCol(); //necessary in case there is I/O at the same time since array is static
00169   
00170   McParticleCol::const_iterator mcParticleTds;
00171 
00172   for (mcParticleTds = mcParticleCnvTds->begin(); mcParticleTds != mcParticleCnvTds->end(); mcParticleTds++) {
00173     Int_t particleID  = (*mcParticleTds)->particleProperty();
00174     Int_t trackIndex  = (*mcParticleTds)->trackIndex();
00175 
00176     Int_t vertexIndex0 = (*mcParticleTds)->vertexIndex0();
00177     Int_t vertexIndex1 = (*mcParticleTds)->vertexIndex1();
00178     Int_t statusFlags = (*mcParticleTds)->statusFlags();
00179 
00180     HepLorentzVector initialPosition = (*mcParticleTds)->initialPosition();
00181     HepLorentzVector finalPosition   = (*mcParticleTds)->finalPosition();
00182     HepLorentzVector initialFourMomentum = (*mcParticleTds)->initialFourMomentum();
00183     //HepLorentzVector finalFourMomentum   = (*mcParticleTds)->finalFourMomentum();
00184 
00185     //cout << " HepID " << (*mcParticleTds)->particleProperty() << endl;
00186     //cout << " init pos " << initialPosition.x() << "  " << initialPosition.y() << "  " << initialPosition.z() << endl;
00187     Int_t mother = -99;
00188     if(!(*mcParticleTds)->primaryParticle()) mother = ((*mcParticleTds)->mother()).trackIndex();
00189  
00190    vector<Int_t> daughters;
00191     SmartRefVector<McParticle> daughterList = (*mcParticleTds)->daughterList();
00192     //SmartRefVector<McParticle>::iterator iter; 
00193     for (int iPar = 0; iPar < daughterList.size(); iPar++) {
00194       //cout <<"daughter Index   " <<daughterList[iPar]->getTrackIndex()<<endl; 
00195       daughters.push_back(daughterList[iPar]->trackIndex());
00196     }
00197     /*    
00198     for (int idau = 0; idau<daughters.size();idau++){
00199       cout <<"daughter Index      " <<daughters[idau]<<endl;
00200     }
00201     
00202     cout<<"###############################"<<endl;
00203     */
00204     TMcParticle *mcParticleRoot = new TMcParticle();
00205     //m_common.m_mcParticleMap[(*mcParticleTds)] = mcParticleRoot;
00206 
00207     mcParticleRoot->setParticleID(particleID);
00208     mcParticleRoot->setTrackIndex(trackIndex);
00209 
00210     mcParticleRoot->setVertexIndex0(vertexIndex0);
00211     mcParticleRoot->setVertexIndex1(vertexIndex1);
00212 
00213     mcParticleRoot->setStatusFlags(statusFlags);
00214 
00215     mcParticleRoot->setInitialPositionX(initialPosition.x());
00216     mcParticleRoot->setInitialPositionY(initialPosition.y());
00217     mcParticleRoot->setInitialPositionZ(initialPosition.z());
00218     mcParticleRoot->setInitialPositionT(initialPosition.t());
00219 
00220     mcParticleRoot->setFinalPositionX(finalPosition.x());
00221     mcParticleRoot->setFinalPositionY(finalPosition.y());
00222     mcParticleRoot->setFinalPositionZ(finalPosition.z());
00223     mcParticleRoot->setFinalPositionT(finalPosition.t());
00224 
00225     mcParticleRoot->setInitialMomentumX(initialFourMomentum.x());
00226     mcParticleRoot->setInitialMomentumY(initialFourMomentum.y());
00227     mcParticleRoot->setInitialMomentumZ(initialFourMomentum.z());
00228     mcParticleRoot->setInitialMomentumE(initialFourMomentum.e());
00229 
00230     //mcParticleRoot->setFinalMomentumX(finalFourMomentum.x());
00231     //mcParticleRoot->setFinalMomentumY(finalFourMomentum.y());
00232     //mcParticleRoot->setFinalMomentumZ(finalFourMomentum.z());
00233     //mcParticleRoot->setFinalMomentumE(finalFourMomentum.e());
00234 
00235     mcParticleRoot->setMother(mother);
00236     mcParticleRoot->setDaughters(daughters);
00237     /*
00238     vector<int> dau = mcParticleRoot->getDaughters();       
00239 
00240     for (int idau = 0; idau<daughters.size();idau++){
00241       cout <<"daughter Index root      " <<dau[idau]<<endl;
00242     } 
00243     */
00244     
00245     McEvt->addMcParticle(mcParticleRoot);
00246   }
00247   return StatusCode::SUCCESS;
00248 }

void RootEventBaseCnv::declareObject (  const std::string &  fullPath,
const CLID &  clid,
const std::string &  treename,
const std::string &  branchname
)  [inherited]
 

Store TDS path to link a particular converter to an object on the TDS.

void RootEventBaseCnv::declareObject (  const std::string &  fullPath,
const CLID &  clid,
const std::string &  treename,
const std::string &  branchname
)  [inherited]
 

Store TDS path to link a particular converter to an object on the TDS. 

00150                                                                                       {
00151     // Purpose and Method:  Save the path on the TDS, treename, pathname in the m_leaves vector, 
00152     //   corresponding to the DataObject that the converter handles.
00153    m_leaves.push_back(RootCnvSvc::Leaf(path, cl,treename,branchname));
00154 }

virtual StatusCode RootEventBaseCnv::fillObjRefs (  IOpaqueAddress *  pAddress,
DataObject *  pObject
)  [virtual, inherited]
 

Resolve the references of the converted object.

StatusCode RootEventBaseCnv::fillObjRefs (  IOpaqueAddress *  pAddress,
DataObject *  pObject
)  [virtual, inherited]
 

Resolve the references of the converted object. 

00117                                                  {
00118     // Purpose and Method:  Resolve the references of the converted object.
00119     //     It is expected that derived classes will override this method.
00120     MsgStream log(msgSvc(), "RootEventBaseCnv");
00121     return StatusCode::SUCCESS;
00122 }

virtual StatusCode RootEventBaseCnv::fillRepRefs (  IOpaqueAddress *  pAddress,
DataObject *  pObject
)  [virtual, inherited]
 

Resolve the references of the converted object.

StatusCode RootEventBaseCnv::fillRepRefs (  IOpaqueAddress *  pAddress,
DataObject *  pObject
)  [virtual, inherited]
 

Resolve the references of the converted object. 

00109                                                  {
00110     // Purpose and Method:  Resolve the references of the converted object.
00111     //     It is expected that derived classes will override this method.
00112     MsgStream log(msgSvc(), "RootEventBaseCnv");
00113     return StatusCode::SUCCESS;
00114 }

virtual StatusCode RootEventBaseCnv::finalize (   )  [virtual, inherited]
 

StatusCode RootEventBaseCnv::finalize (   )  [virtual, inherited]
  

00141                                         {
00142     if ( m_cnvSvc )     {
00143         m_cnvSvc->release();
00144         m_cnvSvc=0;
00145     }
00146     return Converter::finalize();
00147 }

TObject* RootEventBaseCnv::getReadObject (   )  const [inline, inherited]
 

get the object to be read 

00124 { return m_objRead;}

TObject* RootEventBaseCnv::getReadObject (   )  const [inline, inherited]
 

get the object to be read 

00124 { return m_objRead;}

virtual StatusCode RootEventBaseCnv::initialize (   )  [virtual, inherited]
 

Reimplemented in DigiCnv, DstCnv, EvtHeaderCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, TrigCnv, DigiCnv, DstCnv, EvtHeaderCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, and TrigCnv.

StatusCode RootEventBaseCnv::initialize (   )  [virtual, inherited]
 

Reimplemented in DigiCnv, DstCnv, EvtHeaderCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, TrigCnv, DigiCnv, DstCnv, EvtHeaderCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, and TrigCnv.

00125                                           {
00126  
00127   StatusCode status = Converter::initialize();
00128 
00129   if ( status.isSuccess() )   {
00130     IService* isvc = 0;
00131     status = serviceLocator()->service("RootCnvSvc", isvc, false);
00132     if ( !status.isSuccess() )   status = serviceLocator()->service("EventCnvSvc", isvc, true);
00133     if ( status.isSuccess() )   {
00134       status = isvc->queryInterface(IID_IRootCnvSvc, (void**)&m_cnvSvc);
00135     }
00136   }
00137   
00138   return status;
00139 }

template<class Ty1, class Ty2>
destination* Converter< Ty1, Ty2 >::operator (  const source &   )  const [inline, inherited]
 

template<class Ty1, class Ty2>
destination* Converter< Ty1, Ty2 >::operator (  const source &   )  const [inline, inherited]
 

virtual long RootEventBaseCnv::repSvcType (   )  const [inline, virtual, inherited]
  

00087                                     { 
00088       return ROOT_StorageType; 
00089     }

virtual long RootEventBaseCnv::repSvcType (   )  const [inline, virtual, inherited]
  

00087                                     { 
00088       return ROOT_StorageType; 
00089     }

const unsigned char RootEventBaseCnv::storageType (   )  [inline, static, inherited]
 

Storage type and class ID. 

00083                                                   {
00084       return ROOT_StorageType;
00085     }

const unsigned char RootEventBaseCnv::storageType (   )  [inline, static, inherited]
 

Storage type and class ID. 

00083                                                   {
00084       return ROOT_StorageType;
00085     }

virtual StatusCode McParticleCnv::TObjectToDataObject (  DataObject *&  obj  )  [protected, virtual]
 

transformation from root

Implements RootEventBaseCnv.

StatusCode McParticleCnv::TObjectToDataObject (  DataObject *&  obj  )  [protected, virtual]
 

transformation from root

Implements RootEventBaseCnv.

00037                                                                      {
00038   // creation of TDS object from root object
00039 
00040     MsgStream log(msgSvc(), "McParticleCnv");
00041     log << MSG::DEBUG << "McParticleCnv::TObjectToDataObject" << endreq;
00042     StatusCode sc=StatusCode::SUCCESS;
00043 
00044     // create the TDS location for the MdcMc Collection
00045     McParticleCol* mcParticleTdsCol = new McParticleCol;
00046     refpObject=mcParticleTdsCol;
00047 
00048 
00049     // now convert
00050     if (!m_mcParticleCol) return sc;
00051     TIter mcParticleIter(m_mcParticleCol);
00052     TMcParticle *mcParticleRoot = 0;
00053     
00054     vector<int> mothers;
00055     // vector<vector<int> > daughtList;
00056 
00057     while ((mcParticleRoot = (TMcParticle*)mcParticleIter.Next())) {
00058         unsigned int particleID = mcParticleRoot ->getParticleID();
00059         unsigned int trackIndex = mcParticleRoot ->getTrackIndex();
00060 
00061         int vertexIndex0 = mcParticleRoot ->getVertexIndex0();
00062         int vertexIndex1 = mcParticleRoot ->getVertexIndex1();
00063         unsigned int  statusFlags = mcParticleRoot ->getStatusFlags();
00064         
00065         double xInitialPosition = mcParticleRoot->getInitialPositionX();
00066         double yInitialPosition = mcParticleRoot->getInitialPositionY();
00067         double zInitialPosition = mcParticleRoot->getInitialPositionZ();
00068         double tInitialPosition = mcParticleRoot->getInitialPositionT();
00069         
00070         double xFinalPosition =   mcParticleRoot->getFinalPositionX();
00071         double yFinalPosition =   mcParticleRoot->getFinalPositionY();
00072         double zFinalPosition =   mcParticleRoot->getFinalPositionZ();
00073         double tFinalPosition =   mcParticleRoot->getFinalPositionT();
00074         
00075         double xInitialMomentum = mcParticleRoot->getInitialMomentumX();
00076         double yInitialMomentum = mcParticleRoot->getInitialMomentumY();
00077         double zInitialMomentum = mcParticleRoot->getInitialMomentumZ();
00078         double eInitialMomentum = mcParticleRoot->getInitialMomentumE();
00079 
00080         //double xFinalMomentum =   mcParticleRoot->getFinalMomentumX();
00081         //double yFinalMomentum =   mcParticleRoot->getFinalMomentumY();
00082         //double zFinalMomentum =   mcParticleRoot->getFinalMomentumZ();
00083         //double eFinalMomentum =   mcParticleRoot->getFinalMomentumE();
00084 
00085         HepLorentzVector initialMomentum(xInitialMomentum, yInitialMomentum, zInitialMomentum, eInitialMomentum);
00086         HepLorentzVector initialPosition(xInitialPosition, yInitialPosition, zInitialPosition, tInitialPosition);
00087         //HepLorentzVector finalMomentum(xFinalMomentum, yFinalMomentum, zFinalMomentum, eFinalMomentum);
00088         HepLorentzVector finalPosition(xFinalPosition, yFinalPosition, zFinalPosition, tFinalPosition);
00089         
00090         
00091         int mother =  mcParticleRoot->getMother();
00092         vector<int> daughters = mcParticleRoot->getDaughters();
00093 
00094         mothers.push_back(mother);
00095         //daughtList.push_back(daughters);
00096         
00097         McParticle *mcParticleTds = new McParticle;
00098         m_common.m_rootMcParticleMap[mcParticleRoot] = mcParticleTds;
00099         
00100         
00101         mcParticleTds->initialize(particleID, statusFlags, initialMomentum, initialPosition);
00102         mcParticleTds->setTrackIndex(trackIndex);
00103         mcParticleTds->addStatusFlag(statusFlags);
00104 
00105         mcParticleTds->setVertexIndex0(vertexIndex0);
00106         mcParticleTds->setVertexIndex1(vertexIndex1);
00107         //mcParticleTds->finalize(finalMomentum, finalPosition);
00108   mcParticleTds->finalize(finalPosition);
00109 
00110         mcParticleTdsCol->push_back(mcParticleTds);
00111      }
00112         
00113     //Set Mother and DaughterList 
00114     McParticleCol::iterator iter;
00115     int i =0;
00116     for (iter = mcParticleTdsCol->begin(); iter != mcParticleTdsCol->end(); iter++,i++){
00117       //cout<<"    *****   mothers[ "<<i<<"] = "<<mothers[i]<<endl;
00118       if(mothers[i] != -99 ){
00119               McParticleCol::iterator mcParticleTds;
00120         for (mcParticleTds = mcParticleTdsCol->begin(); mcParticleTds != mcParticleTdsCol->end(); mcParticleTds++) {
00121                 int trackIndex = (*mcParticleTds)->trackIndex();
00122                 if( trackIndex == mothers[i] ){ 
00123              (*iter)->setMother(*mcParticleTds);
00124                    (*mcParticleTds)->addDaughter(*iter);
00125                    break;
00126           }
00127                 }
00128       }
00129       else{
00130         (*iter)->setMother(*iter);
00131       }
00132     }
00133     //m_mcParticleCol->Delete();  // wensp add 2005/12/30
00134     delete m_mcParticleCol;
00135     m_mcParticleCol = 0;
00136    return StatusCode::SUCCESS;
00137 }

Friends And Related Function Documentation

CnvFactory<McParticleCnv> [friend]
 

Member Data Documentation

CLID RootEventBaseCnv::CLID_top [protected, inherited]
 

the CLID of the upper converter if any

std::vector<void *> RootEventBaseCnv::m_adresses [protected, inherited]
 

each converter knows the corresponding adresses

std::vector<void *> RootEventBaseCnv::m_adresses [protected, inherited]
 

each converter knows the corresponding adresses

int RootEventBaseCnv::m_branchNr [protected, inherited]
 

the branchNr of this converter for writing

int RootEventBaseCnv::m_branchNrDst [protected, inherited]
 

int RootEventBaseCnv::m_branchNrEvtHeader [protected, inherited]
 

int RootEventBaseCnv::m_branchNrEvtRec [protected, inherited]
 

int RootEventBaseCnv::m_branchNrMc [protected, inherited]
 

int RootEventBaseCnv::m_branchNrRecon [protected, inherited]
 

TArrayS* RootEventBaseCnv::m_branchNumbers [protected, inherited]
 

array with number of branches for reading

TArrayS* RootEventBaseCnv::m_branchNumbers [protected, inherited]
 

array with number of branches for reading

RootCnvSvc* RootEventBaseCnv::m_cnvSvc [protected, inherited]
 

RootCnvSvc* RootEventBaseCnv::m_cnvSvc [protected, inherited]
 

commonData McParticleCnv::m_common [private]
 

relational maps

std::string RootEventBaseCnv::m_currentFileName [protected, inherited]
 

IDataProviderSvc* RootEventBaseCnv::m_eds [protected, inherited]
 

pointer to eventdataservice

IDataProviderSvc* RootEventBaseCnv::m_eds [protected, inherited]
 

pointer to eventdataservice

RootEvtSelector* RootEventBaseCnv::m_evtsel [protected, inherited]
 

RootEvtSelector* RootEventBaseCnv::m_evtsel [protected, inherited]
 

std::vector<RootCnvSvc::Leaf> RootEventBaseCnv::m_leaves [protected, inherited]
 

std::vector<RootCnvSvc::Leaf> RootEventBaseCnv::m_leaves [protected, inherited]
 

TObjArray* McParticleCnv::m_mcParticleCol [private]
 

root object to be read

TObjArray* McParticleCnv::m_mcParticleCol [private]
 

root object to be read

TObject* RootEventBaseCnv::m_objRead [protected, inherited]
 

the object that was read

TObject* RootEventBaseCnv::m_objRead [protected, inherited]
 

the object that was read

std::string RootEventBaseCnv::m_rootBranchname [protected, inherited]
 

root branchname (may be concatenated of severals)

RootInterface* RootEventBaseCnv::m_rootInterface [protected, inherited]
 

pointer to the RootInterface

RootInterface* RootEventBaseCnv::m_rootInterface [protected, inherited]
 

pointer to the RootInterface

std::string RootEventBaseCnv::m_rootTreename [protected, inherited]
 

each converter knows it's treename

The documentation for this class was generated from the following files:
Generated on Wed Feb 2 16:22:35 2011 for BOSS6.5.5 by  doxygen 1.3.9.1