GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » Fairroot » General Discussion » FairRunAna::SetGeomFile seems not to work
FairRunAna::SetGeomFile seems not to work [message #14842] Wed, 05 June 2013 22:08 Go to next message
Volker Friese is currently offline  Volker Friese
Messages: 365
Registered: April 2004
Location: GSI CBM
first-grade participant
From: *gsi.de
The method SetGeomFile of FairRunAna should enable to use a different geometry in the analysis run than was used in the simulation and stored in the parameter database, a feature which is needed e.g. for studying the effects of mis-alignment.

Unfortunately, the implementation seems not to be correct. Whatever geometry file is set in this method, what you get in the end is the geometry from the parameters. This can be seen from the screen output:

Quote:

[INFO ] Opening Geometry input file: data/000.geofile.root
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
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 <TGeoManager::CloseGeometry>: 2640 nodes/ 63 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
[INFO ] The input consists out of the following trees and files:
[INFO ] - cbmsim
[INFO ] - data/x200.mc.root
[INFO ] - FriendTree_1
[INFO ] - data/x200.raw.root
[INFO ] Parameter and input file are available, Assure that basic info is there for the run!
[INFO ] The number of entries in chain is 100

*************************************************************
initialisation for run id 1370261320
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet
Warning in <TGeoManager::Init>: Deleting previous geometry: FAIRGeom/FAIR 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 <TGeoManager::CloseGeometry>: 2640 nodes/ 63 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
Container FairBaseParSet initialized from ROOT file.


Evidently, the geometry is first read from the specified geometry file and then deleted and replaced with the one from the parameter container.

I tried to locate where that happens and followed from FairRunAna::Init() to FairRunTimeDb::initContainers to FairParSet::init() to FairParSet::init(FairParIo*) to FairParGenericSet::init(FairParIo*) to FairParDetIo::init(FairParGenericSet*).... but eventually gave up, being lost in this jungle.
Re: FairRunAna::SetGeomFile seems not to work [message #14844 is a reply to message #14842] Thu, 06 June 2013 09:43 Go to previous messageGo to next message
Mohammad Al-Turany is currently offline  Mohammad Al-Turany
Messages: 518
Registered: April 2004
Location: GSI, Germany
first-grade participant
From: *gsi.de
Hallo Volker,

Yes, this is true! This method was implemented at the time where the the geometry was saved in the data file, and as we move the geometry to the parameters we did not think about it (Sorry!!). Anyway, I created a new Parameter container for the the geometry now and this will be filled by the simulation and case you do not set your geometry manually it will be initialized. So nothing should changed for you except that it should work now! (Fairbase Trunk 20129)

I test it with our tutorials: runing run_digi macro in example/Tutorial3

Quote:


root run_digi.C -l
root [0]
Processing run_digi.C...
FairRootManager::OpenOutFile("data/testdigi.root")
[INFO ] The input consists out of the following trees and files:
[INFO ] - cbmsim
[INFO ] - data/testrun.root
[INFO ] Parameter and input file are available, Assure that basic info is there for the run!
[INFO ] The number of entries in chain is 10
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] No event Header was found!!!
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] Branch: EventHeader. not found in Tree

*************************************************************
initialisation for run id 1370503422
*************************************************************
-I- FairRunTimeDB::InitContainer() FairGeoParSet
Info in <TGeoManager::CloseGeometry>: Geometry loaded from file...
Info in <TGeoManager::SetTopVolume>: Top volume is cave. Master volume is cave
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
Info in <TGeoManager::Voxelize>: Voxelizing...
Info in <TGeoManager::CountLevels>: max level = 1, max placements = 8
Info in <TGeoManager::CloseGeometry>: 9 nodes/ 5 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
Container FairGeoParSet initialized from ROOT file.
-I- FairRunTimeDB::InitContainer() FairBaseParSet
Container FairBaseParSet initialized from ROOT file.

*************************************************************
initialisation for run id 1370503422
*************************************************************
-I- FairRunTimeDB::InitContainer() FairGeoParSet
-I- FairRunTimeDB::InitContainer() FairBaseParSet

*************************************************************
initialisation for run id 1370503422
*************************************************************
-I- FairRunTimeDB::InitContainer() FairGeoParSet
-I- FairRunTimeDB::InitContainer() FairBaseParSet
[INFO ] The number of entries in chain is 10


Macro finished successfully.
Output file is data/testdigi.root
Parameter file is data/testparams.root
Real time 0.039017 s, CPU time 0.03 s

Test passed
All ok
root [1]




and if you add the line fRun->SetGeomFile("data/geofile_full.root") to the macro you get:



Quote:

root run_digi.C -l
root [0]
Processing run_digi.C...
FairRootManager::OpenOutFile("data/testdigi.root")
[INFO ] Opening Geometry input file: data/geofile_full.root
Info in <TGeoManager::CloseGeometry>: Geometry loaded from file...
Info in <TGeoManager::SetTopVolume>: Top volume is cave. Master volume is cave
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
Info in <TGeoManager::Voxelize>: Voxelizing...
Info in <TGeoManager::CountLevels>: max level = 1, max placements = 8
Info in <TGeoManager::CloseGeometry>: 9 nodes/ 5 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
[INFO ] The input consists out of the following trees and files:
[INFO ] - cbmsim
[INFO ] - data/testrun.root
[INFO ] Parameter and input file are available, Assure that basic info is there for the run!
[INFO ] The number of entries in chain is 10
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] No event Header was found!!!
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] Branch: EventHeader. not found in Tree

*************************************************************
initialisation for run id 1370503422
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet
Container FairBaseParSet initialized from ROOT file.

*************************************************************
initialisation for run id 1370503422
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet

*************************************************************
initialisation for run id 1370503422
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet
[INFO ] The number of entries in chain is 10


Macro finished successfully.
Output file is data/testdigi.root
Parameter file is data/testparams.root
Real time 0.0395799 s, CPU time 0.04 s

Test passed
All ok







Cheers,

Mohammad



Still some problems [message #14857 is a reply to message #14844] Mon, 10 June 2013 19:06 Go to previous messageGo to next message
Volker Friese is currently offline  Volker Friese
Messages: 365
Registered: April 2004
Location: GSI CBM
first-grade participant
From: *gsi.de
Hi Mohammad,

thanks for the quick help! I tried with the fairbase trunk, but still have some problems.

When I now set a geometry file to run_ana, it seems to work (see log_with_geofile), although there is a strange error message in the initialisation of FairBaseParSet:

Quote:

*************************************************************
initialisation for run id 1370260005
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet
Error in <TBufferFile::ReadClass>: got wrong class: TGeoManager
Error in <TBufferFile::ReadObject>: got object of wrong class! requested TObjArray but got TGeoManager
Error in <TBufferFile::CheckByteCount>: object of class FairBaseParSet read too few bytes: 69750 instead of 69833
Container FairBaseParSet initialized from ROOT file.


When, however, I try to run without explicitely specifying a geometry file, the program crashes (see log_without_geofile.txt).
Re: Still some problems [message #14858 is a reply to message #14857] Mon, 10 June 2013 21:53 Go to previous message
Mohammad Al-Turany is currently offline  Mohammad Al-Turany
Messages: 518
Registered: April 2004
Location: GSI, Germany
first-grade participant
From: *dip0.t-ipconnect.de
Hallo Volker,

I just try it with the standard macros in cbmroot/macro/run and the trunk of today, it works fine for me. Could it be that something went wrong with switching cbmroot to use the trunk of FairRoot in your local repository!

cheers,

Mohammad


Quote:

run_reco.C...
[INFO ] Opening Geometry input file: geo.root
FairRootManager::OpenOutFile("data/test.eds.root")
Starting CbmMvdDigitizeL::CbmMvdDigitizeL()
CbmTofHitProducerNew instantiated with verbose = 1
Info in <TGeoManager::CloseGeometry>: Geometry loaded from file...
Info in <TGeoManager::SetTopVolume>: Top volume is cave. Master volume is cave
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
Info in <TGeoManager::Voxelize>: Voxelizing...
Error in <TGeoVoxelFinder::SortAll>: Cannot voxelize pole_alu :less than 2 boundaries on Z
Info in <TGeoManager::CountLevels>: max level = 6, max placements = 286
Info in <TGeoManager::CloseGeometry>: 419024 nodes/ 184 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
[INFO ] The input consists out of the following trees and files:
[INFO ] - cbmsim
[INFO ] - data/test.mc.root
[INFO ] Parameter and input file are available, Assure that basic info is there for the run!
[INFO ] The number of entries in chain is 2
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] No event Header was found!!!
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] Branch: EventHeader. not found in Tree

*************************************************************
initialisation for run id 1370892023
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet
Container FairBaseParSet initialized from ROOT file.
-I container name CbmMvdGeoPar
-I container name CbmGeoStsPar
-I container name CbmStsDigiPar
-I container name CbmGeoMuchPar
-I container name CbmGeoTofPar
Error in <FairRuntimeDb::getContainer(Text_t*)>: Container CbmGeoSttPar not created!
create CbmFieldPar container CbmFieldPar
[INFO ] Get the digi parameters for tof
-I container name CbmTofDigiPar






and without geometry file:


Quote:



Processing run_reco.C...
FairRootManager::OpenOutFile("data/test.eds.root")
Starting CbmMvdDigitizeL::CbmMvdDigitizeL()
CbmTofHitProducerNew instantiated with verbose = 1
[INFO ] The input consists out of the following trees and files:
[INFO ] - cbmsim
[INFO ] - data/test.mc.root
[INFO ] Parameter and input file are available, Assure that basic info is there for the run!
[INFO ] The number of entries in chain is 2
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] No event Header was found!!!
[INFO ] Branch: EventHeader. not found in Tree
[INFO ] Branch: EventHeader. not found in Tree

*************************************************************
initialisation for run id 1370892023
*************************************************************
-I- FairRunTimeDB::InitContainer() FairGeoParSet
Info in <TGeoManager::CloseGeometry>: Geometry loaded from file...
Info in <TGeoManager::SetTopVolume>: Top volume is cave. Master volume is cave
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
Info in <TGeoManager::Voxelize>: Voxelizing...
Error in <TGeoVoxelFinder::SortAll>: Cannot voxelize pole_alu :less than 2 boundaries on Z
Info in <TGeoManager::CountLevels>: max level = 6, max placements = 286
Info in <TGeoManager::CloseGeometry>: 419024 nodes/ 184 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
Container FairGeoParSet initialized from ROOT file.
-I- FairRunTimeDB::InitContainer() FairBaseParSet
Container FairBaseParSet initialized from ROOT file.
-I container name CbmMvdGeoPar
-I container name CbmGeoStsPar
-I container name CbmStsDigiPar
-I container name CbmGeoMuchPar
-I container name CbmGeoTofPar
Error in <FairRuntimeDb::getContainer(Text_t*)>: Container CbmGeoSttPar not created!
create CbmFieldPar container CbmFieldPar
[INFO ] Get the digi parameters for tof
-I container name CbmTofDigiPar

*************************************************************
initialisation for run id 1370892023
*************************************************************



Previous Topic: trouble running event display
Next Topic: FairBaseParSet not initialized
Goto Forum:
  


Current Time: Mon Oct 07 22:07:44 CEST 2024

Total time taken to generate the page: 0.00688 seconds