/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/lcio/src/cpp/src/TESTS/test_trackerhit.cc File Reference

#include "tutil.h"
#include "lcio.h"
#include "EVENT/LCIO.h"
#include "IO/LCReader.h"
#include "IO/LCWriter.h"
#include "IMPL/LCEventImpl.h"
#include "IMPL/LCCollectionVec.h"
#include "IMPL/TrackerHitImpl.h"
#include "IMPL/LCFlagImpl.h"

Include dependency graph for test_trackerhit.cc:

Go to the source code of this file.

Functions
int	main (int argc, char **argv)
Variables
static const int	NEVENT = 10
static const int	NHITS = 1000
static string	FILEN = "trackerhits.slcio"
static const string	testname = "test_trackerhit"

---

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 32 of file test_trackerhit.cc.

                                {

    // this should be the first line in your test
    TEST MYTEST=TEST( testname, std::cout );

    try{


        MYTEST.LOG( " writing TrackerHits " );

        // create sio writer
        LCWriter* lcWrt = LCFactory::getInstance()->createLCWriter()  ;

        lcWrt->open( FILEN , LCIO::WRITE_NEW ) ;

        // EventLoop - create some events and write them to the file
        for(int i=0;i<NEVENT;i++){

            // we need to use the implementation classes here 
            LCEventImpl*  evt = new LCEventImpl() ;


            evt->setRunNumber( 4711  ) ;
            evt->setEventNumber( i ) ;

            LCCollectionVec* trkHits = new LCCollectionVec( LCIO::TRACKERHIT )  ;

            for(int j=0;j<NHITS;j++){
                TrackerHitImpl* trkHit = new TrackerHitImpl ;
                //trkHit->setEDep( i*j*117. ) ;
                trkHit->setdEdx( i*j*117. ) ;
                trkHit->setEDepError( (i+j)*.3 ) ;
                double pos[3] = { i, j, i*j } ;
                trkHit->setPosition( pos ) ;
                trkHits->addElement( trkHit ) ;
            }
            evt->addCollection( trkHits , "TrackerHits") ;

            lcWrt->writeEvent(evt) ;

            delete evt ;
        }


        lcWrt->close() ;

        MYTEST.LOG(" reading back TrackerHits from file " ) ;

        // create sio reader
        LCReader* lcRdr = LCFactory::getInstance()->createLCReader()  ;

        lcRdr->open( FILEN ) ;

        for(int i=0;i<NEVENT;i++){

            //std::cout << " testing event " << i << std::endl ;

            LCEvent*  evt =  lcRdr->readNextEvent() ;

            MYTEST( evt->getRunNumber() ,  4711 , " run number " )  ;

            MYTEST(  evt->getEventNumber() ,  i , " event number "  ) ;

            LCCollection* trkHits = evt->getCollection( "TrackerHits") ;

            for(int j=0;j<NHITS;j++) {

                //std::cout << " testing hit " << j << std::endl ;

                TrackerHit* trkHit = dynamic_cast<TrackerHit*>(trkHits->getElementAt(j)) ;

                //MYTEST( trkHit->getEDep() ,  i*j*117. , "EDep" ) ;
                MYTEST( trkHit->getdEdx() ,  i*j*117. , "dEdx" ) ;
                // remove float converstion and check what happens ;)
                MYTEST( trkHit->getEDepError() ,  float((i+j)*.3) , "EDepError" ) ;
                //MYTEST( trkHit->getEDepError() ,  (i+j)*.3 , "EDepError" ) ;

                const double* pos = trkHit->getPosition() ;

                MYTEST( pos[0] , i , " pos[0] " ) ;
                MYTEST( pos[1] , j , " pos[1] " ) ;
                MYTEST( pos[2] , i*j , " pos[2] " ) ;

            }
        }
        lcRdr->close() ;


    } catch( Exception &e ){
        MYTEST.FAILED( e.what() );
    }

    return 0;
}

---

Variable Documentation

const int NEVENT = 10 [static]

Definition at line 22 of file test_trackerhit.cc.

const int NHITS = 1000 [static]

Definition at line 23 of file test_trackerhit.cc.

string FILEN = "trackerhits.slcio" [static]

Definition at line 25 of file test_trackerhit.cc.

const string testname = "test_trackerhit" [static]

Definition at line 28 of file test_trackerhit.cc.

---