RecEmcShowerCnv Class Reference

#include <RecEmcShowerCnv.h>

Inheritance diagram for RecEmcShowerCnv:

Public Member Functions
virtual	~RecEmcShowerCnv ()
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
	RecEmcShowerCnv (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< RecEmcShowerCnv >

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

Definition at line 15 of file RecEmcShowerCnv.h.

Constructor & Destructor Documentation

~RecEmcShowerCnv()

virtual RecEmcShowerCnv::~RecEmcShowerCnv ( )

inlinevirtual

Definition at line 25 of file RecEmcShowerCnv.h.

{ };

RecEmcShowerCnv()

RecEmcShowerCnv::RecEmcShowerCnv ( ISvcLocator *  svc )

protected

Definition at line 34 of file RecEmcShowerCnv.cxx.

: RootEventBaseCnv(classID(), svc)
{
 
  RecEmcShower *aShower = new RecEmcShower();
  //cout<<"in RecEmcShowerCnv::constructor:: showid= "<<aShower->clID()<<endl;
  delete aShower;
 
  // Here we associate this converter with the /Event path on the TDS.
  MsgStream log(msgSvc(), "RecEmcShowerCnv");
  //log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
  //m_rootTreename ="Rec";
  m_rootBranchname ="m_recEmcShowerCol";
  //declareObject(EventModel::Recon::RecEmcShowerCol, objType(), m_rootTreename, m_rootBranchname);
  m_adresses.push_back(&m_recEmcShowerCol);
  m_recEmcShowerCol=0;
}

Member Function Documentation

classID()

static const CLID & RecEmcShowerCnv::classID ( )

inlinestatic

Definition at line 20 of file RecEmcShowerCnv.h.

  {
       return CLID_RecEmcShowerCol; 
  }

DataObjectToTObject()

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

protectedvirtual

transformation to root

Implements RootEventBaseCnv.

Definition at line 223 of file RecEmcShowerCnv.cxx.

{
  MsgStream log(msgSvc(), "RecEmcShowerCnv");
  log << MSG::DEBUG << "RecEmcShowerCnv::DataObjectToTObject" << endreq;
  //cout << "RecEmcShowerCnv::DataObjectToTObject" << endl;
  StatusCode sc=StatusCode::SUCCESS;
 
  RecEmcShowerCol * emcShowerColTds=dynamic_cast<RecEmcShowerCol *> (obj);
  if (!emcShowerColTds) {
    log << MSG::ERROR << "Could not downcast to RecEmcShowerCol" << endreq;
    return StatusCode::FAILURE;
  }
 
  DataObject *evt;
  m_eds->findObject(EventModel::Recon::Event,evt);
  if (evt==NULL) {
    log << MSG::ERROR << "Could not get RecEvent in TDS "  << endreq;
    return StatusCode::FAILURE;
  }
  ReconEvent * devtTds=dynamic_cast<ReconEvent *> (evt);
  if (!devtTds) {
    log << MSG::ERROR << "RecEmcShowerCnv:Could not downcast to TDS DstEvent" << endreq;
  }
  IOpaqueAddress *addr;
 
  m_cnvSvc->getRecTrackCnv()->createRep(evt,addr); 
  TRecTrackEvent *recEvt=m_cnvSvc->getRecTrackCnv()->getWriteObject();
 
  const TObjArray *m_emcShowerCol = recEvt->getEmcShowerCol();
  if (!m_emcShowerCol) return sc;
  recEvt->clearEmcShowerCol(); //necessary in case there is I/O at the same time since array is static
  RecEmcShowerCol::const_iterator emcShowerTds;
 
  for (emcShowerTds = emcShowerColTds->begin(); emcShowerTds != emcShowerColTds->end(); emcShowerTds++) {
    Int_t       trackId = (*emcShowerTds)->trackId();
    Int_t       numHits = (*emcShowerTds)->numHits(); 
    Int_t       status =(*emcShowerTds)->status(); 
    Int_t       cellId =(*emcShowerTds)->getShowerId(); 
    Int_t       clusterId = (*emcShowerTds)->getClusterId(); 
    Int_t       module =(*emcShowerTds)->module(); 
    Double_t    x = (*emcShowerTds)->x();      
    Double_t    y = (*emcShowerTds)->y();      
    Double_t    z = (*emcShowerTds)->z();      
    Double_t    theta = (*emcShowerTds)->theta();  
    Double_t    dtheta =(*emcShowerTds)->dtheta(); 
    Double_t    phi    =(*emcShowerTds)->phi();    
    Double_t    dphi   =(*emcShowerTds)->dphi();   
    Double_t    energy =  (*emcShowerTds)->energy();  
    Double_t    dE    = (*emcShowerTds)->dE();     
    Double_t    eSeed  = (*emcShowerTds)->eSeed();
    Double_t    e3x3  = (*emcShowerTds)->e3x3();
    Double_t    e5x5  = (*emcShowerTds)->e5x5();
    Double_t    eall  = (*emcShowerTds)->getEAll();
    Double_t    elepton  = (*emcShowerTds)->getELepton();
    Double_t    time  = (Int_t)(*emcShowerTds)->time();
    Double_t    secondMoment = (*emcShowerTds)->secondMoment();
    Double_t    latMoment = (*emcShowerTds)->latMoment();
    Double_t    a20Moment = (*emcShowerTds)->a20Moment();
    Double_t    a42Moment = (*emcShowerTds)->a42Moment();
 
    HepSymMatrix matrix   = (*emcShowerTds)->errorMatrix();
    //cout<<"matrix: \n"<<matrix<<endl;
    Double_t    err[6];
    err[0] = matrix[0][0];
    err[1] = matrix[1][1];
    err[2] = matrix[2][2];
    err[3] = matrix[0][1];
    err[4] = matrix[0][2];
    err[5] = matrix[1][2];
 
    map<Int_t, Double_t> cellIdMap;
    RecEmcFractionMap::const_iterator pFractionMap;
 
    for(pFractionMap=(*emcShowerTds)->Begin();
        pFractionMap!=(*emcShowerTds)->End();
        pFractionMap++) {
      Int_t id = pFractionMap->first;
      Double_t frac = pFractionMap->second.getFraction();
      cellIdMap[id]=frac;
    }
 
    vector<Int_t> cellId3x3;
    RecEmcIDVector id3x3 = (*emcShowerTds)->getCellId3x3();
    ci_RecEmcIDVector ciCellId;
    for(ciCellId=id3x3.begin();
        ciCellId!=id3x3.end();
        ciCellId++) {
      Int_t id = *ciCellId;
      cellId3x3.push_back(id);
    }
 
    vector<Int_t> cellId5x5;
    RecEmcIDVector id5x5 = (*emcShowerTds)->getCellId5x5();
    for(ciCellId=id5x5.begin();
        ciCellId!=id5x5.end();
        ciCellId++) {
      Int_t id = *ciCellId;
      cellId5x5.push_back(id);
    }
 
    TRecEmcShower *emcShowerRoot = new TRecEmcShower();
    //m_common.m_recEmcShowerMap[(*emcShowerTds)] = emcShowerRoot;
 
    //emcShowerRoot->initialize(numHits ,status ,cellId ,module ,x ,dx,y ,dy,z ,dz,theta ,dtheta ,phi,dphi ,energy ,dE,cosx ,cosy ,cosz);
    emcShowerRoot->setTrackId(trackId );
    emcShowerRoot->setNumHits(numHits );    
    emcShowerRoot->setStatus(status );     
    emcShowerRoot->setCellId( cellId );     
    emcShowerRoot->setClusterId( clusterId );     
    emcShowerRoot->setModule( module);     
    emcShowerRoot->setX( x    );     
    emcShowerRoot->setY( y    );     
    emcShowerRoot->setZ( z    );     
    emcShowerRoot->setEnergy( energy  );    
    emcShowerRoot->setDE( dE  );     
    emcShowerRoot->setTheta( theta  );  
    emcShowerRoot->setDtheta( dtheta     );  
    emcShowerRoot->setPhi( phi       );  
    emcShowerRoot->setDphi( dphi     );  
    emcShowerRoot->setESeed( eSeed     );
    emcShowerRoot->setE3x3( e3x3    );
    emcShowerRoot->setE5x5( e5x5     );  
    emcShowerRoot->setEAll( eall     );  
    emcShowerRoot->setELepton( elepton     );  
    emcShowerRoot->setTime( time );
    emcShowerRoot->setSecondMoment(secondMoment); 
    emcShowerRoot->setLatMoment(latMoment); 
    emcShowerRoot->setA20Moment(a20Moment); 
    emcShowerRoot->setA42Moment(a42Moment); 
    emcShowerRoot->setErr( err );   
    //emcShowerRoot->setVecHits(vecHits);
 
    emcShowerRoot->setCellIdMap(cellIdMap);
    emcShowerRoot->setCellId3x3(cellId3x3);
    emcShowerRoot->setCellId5x5(cellId5x5);
 
    recEvt->addEmcShower(emcShowerRoot);
  }
  return StatusCode::SUCCESS;
}

TObjectToDataObject()

StatusCode RecEmcShowerCnv::TObjectToDataObject ( DataObject *&  obj )

protectedvirtual

transformation from root

Implements RootEventBaseCnv.

Definition at line 52 of file RecEmcShowerCnv.cxx.

                                                                       {
  // creation of TDS object from root object
 
  MsgStream log(msgSvc(), "RecEmcShowerCnv");
  log << MSG::DEBUG << "RecEmcShowerCnv::TObjectToDataObject" << endreq;
  //cout << "RecEmcShowerCnv::TObjectToDataObject" << endl;
  StatusCode sc=StatusCode::SUCCESS;
 
  IDataProviderSvc* eventSvc;
  Gaudi::svcLocator()->service("EventDataSvc", eventSvc);
  
  SmartDataPtr<RecEmcHitCol> emcRecHitCol(eventSvc,
      EventModel::Recon::RecEmcHitCol);
  if(!emcRecHitCol) log << MSG::INFO << "can't retrieve RecEmcHitCol" << endreq;
  
  SmartDataPtr<RecEmcClusterCol> emcRecClusterCol(eventSvc,
      EventModel::Recon::RecEmcClusterCol);
  if(!emcRecClusterCol) log << MSG::INFO << "can't retrieve RecEmcClusterCol" << endreq;
  
  // create the TDS location for the RecEmcShower Collection
  RecEmcShowerCol* emcShowerTdsCol = new RecEmcShowerCol;
  refpObject=emcShowerTdsCol;
 
  // now convert
  if (!m_recEmcShowerCol) return sc;
  TIter emcShowerIter(m_recEmcShowerCol);
  TRecEmcShower *emcShowerRoot = 0;
  while ((emcShowerRoot = (TRecEmcShower*)emcShowerIter.Next())) {
    int       trackId = emcShowerRoot->trackId();
    int       numHits = emcShowerRoot->numHits(); 
    int       status =  emcShowerRoot->status(); 
    int       cellId =  emcShowerRoot->cellId(); 
    int       module =  emcShowerRoot->module(); 
    double    x =   emcShowerRoot->x();      
    double    y =   emcShowerRoot->y();      
    double    z =   emcShowerRoot->z();      
    double    theta =   emcShowerRoot->theta();  
    double    dtheta =  emcShowerRoot->dtheta(); 
    double    phi    =  emcShowerRoot->phi();    
    double    dphi   =  emcShowerRoot->dphi();   
    double    energy =  emcShowerRoot->energy();  
    double    dE    =   emcShowerRoot->dE();     
    double  eSeed =  emcShowerRoot->eSeed();
    double  e3x3 =  emcShowerRoot->e3x3();
    double  e5x5 =  emcShowerRoot->e5x5();
    double  eall =  emcShowerRoot->eAll();
    double  elepton =  emcShowerRoot->eLepton();
    double  time = emcShowerRoot->time();
    double  secondMoment = emcShowerRoot->secondMoment();
    double  latMoment = emcShowerRoot->latMoment();
    double  a20Moment = emcShowerRoot->a20Moment();
    double  a42Moment = emcShowerRoot->a42Moment();
 
    HepSymMatrix matrix(3);
    matrix[0][0]=emcShowerRoot->err(0);
    matrix[1][1]=emcShowerRoot->err(1);
    matrix[2][2]=emcShowerRoot->err(2);
    matrix[0][1]=emcShowerRoot->err(3);
    matrix[0][2]=emcShowerRoot->err(4);
    matrix[1][2]=emcShowerRoot->err(5);
 
    RecEmcID showerId(cellId);
    RecEmcID clusterId(emcShowerRoot->clusterId());
    
    RecEmcIDVector id3x3;
    vector<Int_t> cellId3x3 = emcShowerRoot->cellId3x3();
    vector<Int_t>::iterator iCellId;
    //cout<<"id3x3"<<endl;
    for(iCellId=cellId3x3.begin();
        iCellId!=cellId3x3.end();
        iCellId++) {
      RecEmcID id(*iCellId);
      //cout<<"id="<<id<<endl;
      id3x3.push_back(id);
    }
 
    RecEmcIDVector id5x5;
    vector<Int_t> cellId5x5 = emcShowerRoot->cellId5x5();
    //cout<<"id5x5"<<endl;
    for(iCellId=cellId5x5.begin();
        iCellId!=cellId5x5.end();
        iCellId++) {
      RecEmcID id(*iCellId);
      //cout<<"id="<<id<<endl;
      id5x5.push_back(id);
    }
 
    RecEmcShower *emcShowerTds = new RecEmcShower();
    m_common.m_rootRecEmcShowerMap[emcShowerRoot] = emcShowerTds;
 
    emcShowerTds->setTrackId(trackId );
    emcShowerTds->setNumHits(numHits );    
    emcShowerTds->setStatus(status );     
    emcShowerTds->setCellId( cellId );     
    emcShowerTds->setModule( module);     
 
    HepPoint3D pos(x,y,z);
    emcShowerTds->setPosition(pos);
    emcShowerTds->setEnergy( energy  );    
    emcShowerTds->setDE( dE  );     
    emcShowerTds->setDtheta( dtheta     );  
    emcShowerTds->setDphi( dphi     );  
    emcShowerTds->setESeed(eSeed);
    emcShowerTds->setE3x3(e3x3);
    emcShowerTds->setE5x5(e5x5);
    emcShowerTds->EAll(eall);
    emcShowerTds->ELepton(elepton);
    emcShowerTds->setTime(time);
    emcShowerTds->setSecondMoment(secondMoment); 
    emcShowerTds->setLatMoment(latMoment); 
    emcShowerTds->setA20Moment(a20Moment); 
    emcShowerTds->setA42Moment(a42Moment); 
    emcShowerTds->setErrorMatrix( matrix );   
 
    emcShowerTds->ShowerId(showerId);     
    emcShowerTds->ClusterId(clusterId);     
    emcShowerTds->CellId3x3(id3x3);
    emcShowerTds->CellId5x5(id5x5);
    
    //put fraction map into shower
    if(emcRecHitCol) {
      map<Int_t, Double_t> cellIdMap=emcShowerRoot->cellIdMap();
      map<Int_t, Double_t>::iterator iCellIdMap;
      for(iCellIdMap=cellIdMap.begin();
          iCellIdMap!=cellIdMap.end();
          iCellIdMap++) {
        
        RecEmcID id(iCellIdMap->first);
        
        RecEmcHitCol::iterator iHit;
        for(iHit=emcRecHitCol->begin();
            iHit!=emcRecHitCol->end();
            iHit++) {
          
          RecEmcID idHit((*iHit)->getCellId());
 
          if(id==idHit) {
            RecEmcFraction frac(*(*iHit));
            frac.Fraction(iCellIdMap->second);
            emcShowerTds->Insert(frac);
            break;
          }
        } //RecEmcHitCol
      } //CellIdMap
    }
 
    if(emcRecClusterCol) {
      RecEmcClusterCol::iterator iCluster;
      for(iCluster=emcRecClusterCol->begin();
          iCluster!=emcRecClusterCol->end();
          iCluster++) {
        if(clusterId==(*iCluster)->getClusterId()) {
          emcShowerTds->Cluster(*iCluster);
          break;
          //emcRecClusterCol->InsertShower(*iCluster);
        }
      }
    }
 
    emcShowerTdsCol->push_back(emcShowerTds);
    //     delete emcShowerTds;
    // emcShowerTds = NULL;
  }
 
  //m_recEmcShowerCol->Delete();  // wensp add 2005/12/30
  delete m_recEmcShowerCol;
  m_recEmcShowerCol = 0;  
 
  return StatusCode::SUCCESS;
}

Friends And Related Function Documentation

CnvFactory< RecEmcShowerCnv >

friend class CnvFactory< RecEmcShowerCnv >

friend

Definition at line 9 of file RecEmcShowerCnv.h.

---

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

• RecEmcShowerCnv.h
• RecEmcShowerCnv.cxx