GSI Forum - RDF feed
type: channel
description: GSI Helmholtzzentrum für Schwerionenforschung
items: Running the Kalman on the Mvd
items: Re: Running the Kalman on the Mvd
items: Re: Running the Kalman on the Mvd (solved dependency)

Running the Kalman on the Mvd
type: item
description: Hi all, I started to try out genfit on MvdHits. After some issues concerning the Mvd package itself I ran into some fitting problems which don't tell me much. Here we go: I checked in the stuff into rev. 2589 (and 2588). I run the macros in macro/mvd/Tracking The simulation does 1GeV protons with Mvd & Field The digi & hitreco runs MC Trackfinding runs Anyway somehow the Kalman cannot fit the tracks I make. Here you have some detailed output (I uncommented the info output in the Kalman class). Hopefully someone can help... Ralf. Toggle Spoiler PndMvdKalmanTask::Exec -I- PndMvdKalmanTask: contains 1 Tracks. Detailed Debug info on the tracks: TrackCand no. 0 has 6 hits. [ ihit | detid | index ] [ 0 | 3 | 0 ] [ 1 | 3 | 1 ] [ 2 | 3 | 2 ] [ 3 | 4 | 0 ] [ 4 | 4 | 1 ] [ 5 | 4 | 2 ] starting track0 -I- PndMvdRecoHit::PndMvdRecoHit(PndMvdHit*) called. Mvd hit in detector 1_1/34_0/119_1/118_2/98_14/92_1/90_1/89_1/ at (1.60909, 0.271733, 1.94354) cm with 46929.8 e, Cluster No. 0 o: 2.11 0.656733 1.94354 u: 0 1 0 v: -1 0 0 -I- PndMvdRecoHit::PndMvdRecoHit(PndMvdHit*) called. Mvd hit in detector 1_1/34_0/102_1/101_3/98_1/92_1/90_1/89_1/ at (4.16446, 0.653887, 5.0215) cm with 92495.9 e, Cluster No. 1 o: 4.08241 1.43459 4.5 u: -0.104528 0.994522 0 v: 0 0 1 -I- PndMvdRecoHit::PndMvdRecoHit(PndMvdHit*) called. Mvd hit in detector 1_1/34_0/66_2/65_12/59_1/57_1/56_1/ at (4.9599, 0.758694, 5.975) cm with 32364.3 e, Cluster No. 2 o: 4.10249 -1.00282 5.975 u: 0.707107 0.707107 0 v: -0.707107 0.707107 0 -I- PndMvdRecoHit::PndMvdRecoHit(PndMvdHit*) called. Mvd hit in detector 1_1/34_0/80_1/78_1/77_4/75_29/72_1/71_1/ at (7.89362, 1.08679, 9.48283) cm with 122466 e, Cluster No. 0 o: 7.67364 1.9824 10.65 u: 0 0 1 v: 0.238533 -0.971134 0 -I- PndMvdRecoHit::PndMvdRecoHit(PndMvdHit*) called. Mvd hit in detector 1_1/34_0/80_1/79_1/77_9/75_29/72_1/71_1/ at (7.44033, 1.04598, 8.9416) cm with 66616.9 e, Cluster No. 2 o: 7.4256 -0.00891832 10.65 u: 0 0 1 v: -0.0139622 -0.999903 0 -I- PndMvdRecoHit::P…
creator: Ralf Kliemt
date: 2008-04-24T09:54:20-00:00

Re: Running the Kalman on the Mvd
type: item
description: Hi Ralf! It looks to me as if you would still use the LSLTrackRep. For your exercise you should use GeaneTrackRep! Furthermore you have to make sure that you have reasonable starting values for your fit, also for the covariance matrix. How do you initialize the TrackRep? Cheers! Sebastian.
creator: Sebastian Neubert
date: 2008-04-24T13:51:10-00:00

Re: Running the Kalman on the Mvd
type: item
description: Hello again, Thanks for the hint to use GeaneTrackRep. Of course I still cannot run the Kalman, so I'll ask my question: When making a new CbmGeanePro() the program breaks. with: Toggle Spoiler PndFieldCreator::createCbmField() Program received signal SIGSEGV, Segmentation fault. [Switching to Thread -1224078560 (LWP 28569)] 0xb3f81377 in CbmGeanePro (this=0xa25b468) at /home/ralfk/Pandaroot/pandaroot/geane/CbmGeanePro.cxx:35 35 afErtrio=gMC3->fErtrio; (gdb) bt #0 0xb3f81377 in CbmGeanePro (this=0xa25b468) at /home/ralfk/Pandaroot/pandaroot/geane/CbmGeanePro.cxx:35 #1 0xb3d9b052 in PndMvdKalmanTask::Init (this=0x872ac08) at /home/ralfk/Pandaroot/pandaroot/mvd/MvdTracking/PndMvdKalmanTask.cxx:115 #2 0xb4bb72ef in CbmTask::InitTask (this=0x872ac08) at /home/ralfk/Pandaroot/pandaroot/base/CbmTask.cxx:36 #3 0xb4bb73d4 in CbmTask::InitTasks (this=0x8562078) at /home/ralfk/Pandaroot/pandaroot/base/CbmTask.cxx:101 #4 0xb4bb738a in CbmTask::InitTask (this=0x8562078) at /home/ralfk/Pandaroot/pandaroot/base/CbmTask.cxx:43 #5 0xb4bb35f3 in CbmRunAna::Init (this=0x8561fe0) at /home/ralfk/Pandaroot/pandaroot/base/CbmRunAna.cxx:114 #6 0xb4bf0cfb in G__CbmDict_531_0_4 (result7=0xbfd712ac, funcname=0x855f1b8 "\001", libp=0xbfd6d808, hash=0) at /home/ralfk/Pandaroot/build/base/CbmDict.cxx:9361 #7 0xb73fd5fb in Cint::G__ExceptionWrapper (funcp=0xb4bf0cd4 <G__CbmDict_531_0_4>, result7=0xbfd712ac, funcname=0x855f1b8 "\001", libp=0xbfd6d808, hash=0) at cint/src/Api.cxx:364 #8 0xb750d87a in G__call_cppfunc (result7=0xbfd712ac, libp=0xbfd6d808, ifunc=0x855f1b8, ifn=0) at cint/src/v6_newlink.cxx:512 #9 0xb74d4855 in G__interpret_func (result7=0xbfd712ac, funcname=0xbfd70dac "Init", libp=0xbfd6d808, hash=404, p_ifunc=0x855f1b8, funcmatch=1, memfunc_flag=1) at cint/src/v6_ifunc.cxx:5133 #10 0xb74c3e9e in G__getfunction (item=0xbfd74d86 "Init()", known3=0xbfd7334c, memfunc_flag=1) at cint/src/v6_func.cxx:2517 #11 0xb759b6d8 in G__getstructmem (store_var_type=112, varname=0xbfd72790 "'׿\025_\0…
creator: Ralf Kliemt
date: 2008-04-28T10:44:22-00:00

Re: Running the Kalman on the Mvd
type: item
description: Hm... That does not ring a bell... My only suggestion is, that you look into the demo to see what preparations are done there to set up Geane. Cheers! Sebastian.
creator: Sebastian Neubert
date: 2008-04-28T10:48:55-00:00

Re: Running the Kalman on the Mvd
type: item
description: I am also facing the same problem Ralf have mentioned here. It looks to me TGeant3 has pointer data member fErtrio, which is not handeled properly in the CbmGeanePro's Default constructor.
creator: Vanniarajan Suyam Jothi
date: 2008-04-28T14:35:22-00:00

Re: Running the Kalman on the Mvd (solved dependency)
type: item
description: Hi again, I solved my library problems to start the fitter: Add ${GEANT3_INCLUDE_DIR} ${CMAKE_SOURCE_DIR}/genfit ${CMAKE_SOURCE_DIR}/recotasks ${CMAKE_SOURCE_DIR}/trackbase ${CMAKE_SOURCE_DIR}/trackrep ${CMAKE_SOURCE_DIR}/geane to CMakeLists.txt Make a CbmGeane object in your macro Toggle Spoiler CbmGeane *Geane = new CbmGeane(MCFile); PndConstField *fMagField=new PndConstField(); fMagField->SetField(0, 0 ,20. ); // values are in kG fMagField->SetFieldRegion(-500, 500,-500, 500, -200, 200); Geane->SetField(fMagField); At my current stage I have to apply a correct handling of the MvdRecoHit. Greetings again, Ralf.
creator: Ralf Kliemt
date: 2008-04-28T15:24:09-00:00

Re: Running the Kalman on the Mvd (solved dependency)
type: item
description: Hello Ralf, Can you tell me in which CMakeList you have made this Changes. Cheers, Vanni
creator: Vanniarajan Suyam Jothi
date: 2008-05-02T12:22:19-00:00

Re: Running the Kalman on the Mvd (solved dependency)
type: item
description: Hi Vanni, I made the additions to the pandaroot/mvd/CMakeLists.txt. It should be connected to the library where the kalman task is located. Regards, Ralf. Toggle Spoiler # Create a library called "" which includes the source files given in # the array . # The extension is already found. Any number of sources could be listed here. set(INCLUDE_DIRECTORIES ${ROOT_INCLUDE_DIR} ${GEANT3_INCLUDE_DIR} ${CLHEP_INCLUDE_DIR} ${CMAKE_SOURCE_DIR}/geobase ${CMAKE_SOURCE_DIR}/parbase ${CMAKE_SOURCE_DIR}/base ${CMAKE_SOURCE_DIR}/field ${CMAKE_SOURCE_DIR}/mcstack ${CMAKE_SOURCE_DIR}/passive ${CMAKE_SOURCE_DIR}/mvd/MvdMC ${CMAKE_SOURCE_DIR}/mvd/MvdData ${CMAKE_SOURCE_DIR}/mvd/MvdDigi ${CMAKE_SOURCE_DIR}/mvd/PndMvdReco ${CMAKE_SOURCE_DIR}/mvd/MvdTools ${CMAKE_SOURCE_DIR}/mvd/transport ${CMAKE_SOURCE_DIR}/mvd ${CMAKE_SOURCE_DIR}/genfit ${CMAKE_SOURCE_DIR}/recotasks ${CMAKE_SOURCE_DIR}/trackbase ${CMAKE_SOURCE_DIR}/trackrep ${CMAKE_SOURCE_DIR}/geane ) include_directories( ${INCLUDE_DIRECTORIES}) set(LINK_DIRECTORIES ${ROOT_LIBRARY_DIR} ) link_directories( ${LINK_DIRECTORIES}) ########################### create libMvd ###################################### set(MVD_SRCS MvdData/PndMvdDigi.cxx MvdData/PndMvdDigiPixel.cxx MvdData/PndMvdDigiStrip.cxx MvdData/PndMvdMCPoint.cxx MvdDigi/PndMvdCalcFePixel.cxx MvdDigi/PndMvdCalcPixel.cxx MvdDigi/PndMvdCalcStrip.cxx MvdDigi/PndMvdDigiAna.cxx MvdDigi/PndMvdHybridHitProducer.cxx MvdDigi/PndMvdStripHitProducer.cxx MvdDigi/PndMvdHitProducerIdeal.cxx MvdDigi/PndMvdDigiTask.cxx MvdMC/PndMvdContFact.cxx MvdMC/PndMvdDetector.cxx MvdMC/PndMvdGeo.cxx MvdMC/PndMvdGeoPar.cxx MvdMC/PndMvdGeoMappingPar.cxx MvdMC/PndMvdMaterialDet.cxx MvdTools/PndMvdPixelDigiPar.cxx MvdTools/PndMvdStripDigiPar.cxx MvdTools/PndStringVector.cxx MvdTools/PndMvdGeoHandling.cxx MvdTools/PndMvdEventAna.cxx MvdTools/PndMvdMCEventAna.cxx MvdTools/PndMvdDigiEventAna.cxx MvdTools/PndMvdAllDataEventAna.cxx MvdTools/PndMvdFileNameCreator.cxx MvdTools/PndGeoHitList.cxx MvdTools/PndEventDisp…
creator: Ralf Kliemt
date: 2008-05-02T12:47:42-00:00