GSI Forum - RDF feed
type: channel
description: GSI Helmholtzzentrum für Schwerionenforschung
items: Field maps in Macros
items: Re: Field maps in Macros

Field maps in Macros
type: item
description: Hi, I have simplified the PndMultiField, it has now a new ctor: PndMultiField(TString Map) The Map can be: 1. Full which is equivalent to adding all maps in the macro i.e: PndDipoleMap *map_d1= new PndDipoleMap("DipoleMap1", "R"); PndDipoleMap *map_d2= new PndDipoleMap("DipoleMap2", "R"); PndSolenoidMap *map_s1= new PndSolenoidMap("SolenoidMap1", "R"); PndSolenoidMap *map_s2= new PndSolenoidMap("SolenoidMap2", "R"); PndSolenoidMap *map_s3= new PndSolenoidMap("SolenoidMap3", "R"); PndSolenoidMap *map_s4= new PndSolenoidMap("SolenoidMap4", "R"); AddField(map_t); AddField(map_d1); AddField(map_d2); AddField(map_s1); AddField(map_s2); AddField(map_s3); AddField(map_s4); 2. Dipole: equivalent to: PndDipoleMap *map_d1= new PndDipoleMap("DipoleMap1", "R"); PndDipoleMap *map_d2= new PndDipoleMap("DipoleMap2", "R"); AddField(map_d1); AddField(map_d2); 3. Solenoid: equivalent to PndSolenoidMap *map_s1= new PndSolenoidMap("SolenoidMap1", "R"); PndSolenoidMap *map_s2= new PndSolenoidMap("SolenoidMap2", "R"); PndSolenoidMap *map_s3= new PndSolenoidMap("SolenoidMap3", "R"); PndSolenoidMap *map_s4= new PndSolenoidMap("SolenoidMap4", "R"); AddField(map_s1); AddField(map_s2); AddField(map_s3); AddField(map_s4); so in the macro one can simply write: Quote: //---------------------Create and Set the Field(s)---------- PndMultiField *fField= new PndMultiField("FULL"); fRun->SetField(fField); also the parameter handling for the field is done now in the field itself so the lines: (PndMultiFieldPar*) rtdb->getContainer("PndMultiFieldPar"); if(fField) { fieldPar->SetParameters(fField); } fieldPar->setInputVersion(fRun->GetRunId(),1); fieldPar->setChanged(kTRUE); are not needed any more (see macro/run/run_sim.C ) regards Mohammad
creator: Mohammad Al-Turany
date: 2009-12-01T18:31:44-00:00

Re: Field maps in Macros
type: item
description: Cool!
creator: StefanoSpataro
date: 2009-12-01T18:44:57-00:00

Re: Field maps in Macros
type: item
description: Dear Mohammad, i have to use PndConstfield since i have to adapt the magnetic field value to low momenta particle. after the modifications of PndMultifield, the setting of PndConstField does not work anymore. please run at macro/hyp sim_pid.C hit_pid.C If i have missed something, please let me know how to proceed properly, thanks in advance ALicia. error message: root -l hit_pid.C >test.txt Warning in <TClassTable::Add>: class PndSttTrack already in TClassTable Warning in <TClassTable::Add>: class PndMvdPixel already in TClassTable Warning in <TClassTable::Add>: class PndMvdStrip already in TClassTable Warning in <TClassTable::Add>: class PndMvdMCPoint already in TClassTable Warning in <TClassTable::Add>: class PndMvdDigi already in TClassTable Warning in <TClassTable::Add>: class PndMvdDigiPixel already in TClassTable Warning in <TClassTable::Add>: class PndMvdApvHit already in TClassTable Warning in <TClassTable::Add>: class PndMvdPidCand already in TClassTable Warning in <TClassTable::Add>: class PndStringVector already in TClassTable Info in <TGeoManager::TGeoManager>: Geometry Geometry, default geometry created Warning in <TGeoManager::Init>: Deleting previous geometry: Geometry/default geometry Info in <TGeoManager::CloseGeometry>: Geometry loaded from file... Info in <TGeoManager::SetTopVolume>: Top volume is cave. Master volume is cave Info in <TGeoManager::Voxelize>: Voxelizing... Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100 Info in <TGeoManager::CloseGeometry>: 247 nodes/ 247 volume UID's in FAIR geometry Info in <TGeoManager::CloseGeometry>: ----------------modeler ready---------------- *** Break *** segmentation violation Generating stack trace... 0x00007f306a0dc95e in PndFieldCreator::createFairField() at /home/lorente/work/trunk/field/PndFieldCreator.cxx:89 from /home/lorente/work/build/lib/libField.so 0x00007f306ade1f6d in FairRunAna::Init() at /home/lorente/work/trunk/base/FairRunAna.cxx:200 from /home/lorente/work/build…
creator: asanchez
date: 2009-12-03T16:33:52-00:00

Re: Field maps in Macros
type: item
description: Hi, It should work now, I implemented the new handling for the multi field but for all the single field (Sorry!), anyway now it is also implemented for the const field. I also made a change in your macro, the parameter of the field are done internally now, you do not need to save them your self. regards Mohammad
creator: Mohammad Al-Turany
date: 2009-12-03T18:40:16-00:00

Re: Field maps in Macros
type: item
description: thank you very much,
creator: asanchez
date: 2009-12-03T21:29:52-00:00

Re: Field maps in Macros
type: item
description: Hi all, as a reminder to all the developers... please clean up your macros and apply the simplification introduced by mr. M. Greetings, Johan.
creator: Johan Messchendorp
date: 2010-01-13T19:38:25-00:00

Re: Field maps in Macros
type: item
description: Dear collegues, After the refresh the pandaroot v7467, /macro/run/sim_complete_tpc.C and /macro/emc/sim_emc.C are tested, but they simply doesn't work due to magnet field. The field definition have been already replaced as you recommended. fRun->SetBeamMom(15); PndMultiField *fField= new PndMultiField("FULL"); fRun->SetField(fField); If I exclude magnet field in the script, everything looks fine. Please check macro, specially the part of magnet accessing. Thank you in advance. Quote: ============== FairRunSim: Initialising simulation run ============== Info in <TGeoManager::TGeoManager>: Geometry FAIRGeom, FAIR geometry created -I- FairGeoMedia Read media -I container name PndEmcDigiPar ************************************************************* initialisation for run id 1890884962 ************************************************************* PndEmcDigiPar initialized from Ascii file -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/TransMap.1500.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/DipoleMap1.1500.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/DipoleMap2.1500.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap1.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap2.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap3.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap4.root create PndFieldPar container PndMultiFieldPar create PndFieldPar container PndTransPar create PndFieldPar container PndDipole1Par create PndFieldPar container PndDipole2Par create PndFieldPar container PndSolenoid1Par create PndFieldPar container PndSolenoid2Par create PndFieldPar container PndSolenoid3Par create PndFieldPar container PndSolenoid4Par *** Break *** segmentation …
creator: donghee
date: 2010-01-15T13:34:00-00:00

Re: Field maps in Macros
type: item
description: HI, I just test them again and all work! I also change the field setting in the sim_emc.C and emc_complete.C and it runs without any problems. Maybe your parameter files are corrupted, try to delete the old parameter root files from your run and emc directory, you can also use the gdb so that one can see where it crashes. i.e: gdb root.exe gdb> r root> .x sim_complete_tpc .C regards Mohammad
creator: Mohammad Al-Turany
date: 2010-01-15T13:53:06-00:00

Re: Field maps in Macros
type: item
description: Dear collegues, The parameter root files are deleted from run and emc directory before execution of macro script, and running again in gdb. The output error messages for macro/run/sim_complete_tpc.C are still there. Quote: PndEmcDigiPar initialized from Ascii file -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/TransMap.1500.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/DipoleMap1.1500.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/DipoleMap2.1500.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap1.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap2.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap3.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap4.root create PndFieldPar container PndMultiFieldPar create PndFieldPar container PndTransPar create PndFieldPar container PndDipole1Par create PndFieldPar container PndDipole2Par create PndFieldPar container PndSolenoid1Par create PndFieldPar container PndSolenoid2Par create PndFieldPar container PndSolenoid3Par create PndFieldPar container PndSolenoid4Par Program received signal SIGSEGV, Segmentation fault. [Switching to Thread 0xb6b1a6c0 (LWP 10442)] 0xb3d79ac4 in ROOT::Init68 () from /home/donghee/GSI/buildpanda/lib/libField.so Detail backtrace of all stack frames shows failure of field map reading at below. In ordet to do very simple test, pandaroot trunk has been downloaded, installed and directly example macros are excuted without any modification. Every script have same problem in my current OS, in which I'm using now Ubuntu 32bit. Do you have some more idea for this problem? Thank you for your help! Donghee Quote: (gdb) bt #0 0xb3d79ac4 in ROOT::Init68 () from /home/donghee/GSI/buildpanda/lib/libF…
creator: donghee
date: 2010-01-15T20:52:14-00:00

Re: Field maps in Macros
type: item
description: Dear collegues, I have tested the pandaroot with GSI etch 64bit. It works fine, but only in ubuntu-32bit macro script have trouble. Is there a simple solution? Thanks...
creator: donghee
date: 2010-01-17T12:03:02-00:00

Re: Field maps in Macros
type: item
description: Dear Donghee, Unfortunately, I cannot test it, since my MAC machine with various linux Virtual Machines died last year and I am waiting for a new one. Does anyone else has a Ubuntu 32-bit (Lenny or Etch, please specify???) and would be able to test this, or even better, place it on the Dashboard? Greetings, Johan.
creator: Johan Messchendorp
date: 2010-01-17T12:32:20-00:00

Re: Field maps in Macros
type: item
description: In my SL4.7 32bit I am able to run fine the macro. The only thing is that I have thousands of messages: PndTpcDetector::ProcessHits: EXIT q==0 Is it normal? I think it would be better to remove these messages.
creator: StefanoSpataro
date: 2010-01-18T08:39:52-00:00

Re: Field maps in Macros
type: item
description: Felix, what to do with these TPC message. Could you remove them, or are they serious and do they need to be solved. Please your comments, Johan.
creator: Johan Messchendorp
date: 2010-01-18T08:50:01-00:00

Re: Field maps in Macros
type: item
description: Johan Messchendorp wrote on Sun, 17 January 2010 13:32 Dear Donghee, Unfortunately, I cannot test it, since my MAC machine with various linux Virtual Machines died last year and I am waiting for a new one. Does anyone else has a Ubuntu 32-bit (Lenny or Etch, please specify???) and would be able to test this, or even better, place it on the Dashboard? Greetings, Johan. Hi Johan, I have got a 64-bit Debian 5.x (Lenny). Is nightly enough? I can set a 32-bit version available, from the next week, or is it too late? Greetings,
creator: M.Babai
date: 2010-01-18T13:35:50-00:00

Re: Field maps in Macros
type: item
description: Sorry for replying so late (I was on vacation and now I am ill). The message you see is a harmless warning. Everything is working fine. I removed this debug output.
creator: Felix Boehmer
date: 2010-01-21T12:29:24-00:00

Re: Field maps in Macros
type: item
description: Dear Panda Developers, Since 2 months ago I have still trouble to run all scripts in Pandaroot under Ubuntu 8.04.2 hardy 32bit. During the MC simulation, error message looks following. Could somebody test it in virtual machine. If you can give some comment ,any suggestion, and any semi-solution, then would be happy at this moment. I guess that there are some kind of mishandling for memory in the function of SetFieldContainer or too large/short memory allocation, especailly with 32bit OS. The correct issues and indications for this problem are shown at the end of error message, but unfortunately I couldn't reach to get a solution with my poor knowledge of OS system. I know already a simple solution. Please, don't say something like this! "Donghee, You can simply use 64bit or upgrade 8.04 to 8.10, that is already tested and properly, nicely working!, let's thinking about it!" I know this that is a final solution. not right now! I'm really sorry for my stupid question! Best wishes, Donghee Quote: ************************************************************* initialisation for run id 1838244001 ************************************************************* FairRunSim::Init() create visualisation manager -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/TransMap.1500.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/DipoleMap1.1500.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/DipoleMap2.1500.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap1.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap2.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap3.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap4.root create PndFieldPar container PndMultiFieldPar create PndField…
creator: donghee
date: 2010-02-18T16:07:17-00:00

Re: Field maps in Macros
type: item
description: For this thing I am out of the business, I hope somebody else will comment and solve it. Sorry.
creator: StefanoSpataro
date: 2010-02-18T16:48:52-00:00

Re: Field maps in Macros
type: item
description: Dear Donghee, For this, I would refer to the following thread on the forum: http://forum.gsi.de/index.php?t=tree&goto=10217&rid=981&S=94 6a02b8347d89e91707581eef92e617#msg_10217 Mohammad B., can we make Ubuntu 8.04, 32-bit available as a virtual machine and give Donghee access? Johan.
creator: Johan Messchendorp
date: 2010-02-18T23:58:54-00:00

Re: Field maps in Macros
type: item
description: Hi Donghee Could you please try with an older version of the external packages, e.g. jul09. I had already a problem where Root crashes with some similar crazy error messages and this happens only with the new external packages. This will not solve the general problem but it may help you to go on with your work. Ciao Florian
creator: Florian Uhlig
date: 2010-02-19T07:33:05-00:00

Re: Field maps in Macros
type: item
description: Is it possible that the field root files should be created with the new root version?
creator: StefanoSpataro
date: 2010-02-19T09:36:14-00:00

Re: Field maps in Macros
type: item
description: Johan Messchendorp wrote on Fri, 19 February 2010 00:58 Dear Donghee, For this, I would refer to the following thread on the forum: http://forum.gsi.de/index.php?t=tree&goto=10217&rid=981&S=94 6a02b8347d89e91707581eef92e617#msg_10217 Mohammad B., can we make Ubuntu 8.04, 32-bit available as a virtual machine and give Donghee access? Johan. No prob., today, I will install a 32-bit Ubuntu 8.10. Donghee, I will contact you asa the installation is ready to use. Best wishes, /MB
creator: M.Babai
date: 2010-02-19T10:18:12-00:00

Re: Field maps in Macros
type: item
description: Dear Donghee, We could reproduce your problem related to reading the field maps on a virtual machine running Ubuntu 8.04. However, so far we could not find a quick solution for your problem. It might well be that we find a solution in the course of time, but we will not put a very high priority to focus on solving these problems, in particular, since fairroot and pandaroot is already supported on various platforms, and there are so many other issues with high priority to tackle. For a list of supported sites, please take a look at the fairroot website (http://fairroot.gsi.de). You have to understand that it always remains a nightmare (and a lot of effort) to make it compatible for everything, although we (well, in particular the fairroot-team) are trying hard. Therefore, we will continue to find the origin of the problem, but it will be on a very low priority. As a general remark, I would advise those who wants to run the code "out-of-the-box" should to try to stick to the sites which are "officially" supported. This does not mean that we are not interested in the compatibility with "other" platforms etc, but a quick solution cannot be expected. Best wishes, Johan.
creator: Johan Messchendorp
date: 2010-02-21T16:11:33-00:00

Re: Field maps in Macros
type: item
description: Dear Johan, I can fully understand all situation, and agree on your comment about the issue of high or low priority. I'm already satisfied for any kind help and all efforts of developers. One more remark is that I have tried to update ubuntu 8.04 to 8.10(32bit). It doesn't help as I expected, so one of simple solution is to use new OS platform with 64bit, specially ubuntu case. Best regards, Donghee Kang
creator: donghee
date: 2010-02-22T09:53:54-00:00