FEEMC from ROOT file [message #5401] |
Wed, 14 November 2007 18:41 |
Aleksandra Biegun
Messages: 64 Registered: May 2007 Location: Groningen
|
continuous participant |
From: *KVI.nl
|
|
Hi All,
I just commit to svn a macro for creating a ROOT file for
forward end-cap of emc (FEEMC) "createRootGeoFileFwEnd.C"
which is in
/pandaroot/macro/emc/
directory. This end-cap geometry is created according to
design made by Herbert Loehner and Henk Smit at KVI
(discussed and accepted by PANDA collaboration at CM in Dubna).
I also adapted the PndEmc.cxx (*.h) code for it.
This endcap is also called "module 3",
so now one should not use both inputs for this endcap
at the same time.
If someone wants to run simulation with barrel and backward endcap (from "emc_module1245.dat", without module 3) and forward endcap (from "emc_module3.root") one should called twice PndEmc class:
CbmDetector *Emc = new PndEmc("EMC",kTRUE);
Emc->SetGeometryFileName("emc_module1245.dat");
Emc->SetVerboseLevel(3);
fRun->AddModule(Emc);
CbmDetector *Emc = new PndEmc("EMC",kTRUE);
Emc->SetGeometryFileName("emc_module3.root");
Emc->SetVerboseLevel(3);
fRun->AddModule(Emc);
In attachment you can find plots for Points when variables,
expecially energy of points, are compared for previous
(black) and new (ROOT) version (red) of forward end-cap.
If you have any comments or problems with running simulation
with the new file, please let me know.
Best regards,
Aleksandra.
Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
|
|
|
FEEMC from ROOT file [message #5618 is a reply to message #5401] |
Thu, 10 January 2008 17:24 |
Aleksandra Biegun
Messages: 64 Registered: May 2007 Location: Groningen
|
continuous participant |
From: *KVI.nl
|
|
Hi "Emc Friends",
I just submitted to svn changed macro for creation geometry for forward endcap (FwEndCap) of emc at
/pandaroot/macro/emc/ called createRootGeoFileFwEnd_changed.C (rev. 1996).
The way of creation each of 4 volumes - top, QuarterVol, SubunitVol, BoxVol (except for CrystalVol) is changed from TGeoVolume(name,shape,medium) to TGeoVolumeAssembly(name). The reason was a problem with G4, which could not recognize medium (even air) always for the most outside volume (G3 also had the same problem, but it went over the problem while G4 did not). By changing volumes into TGeoVolumeAssembly(name) we can also avoid overlaps.
There are also PndEmc.cxx (h) classes submitted (rev. 1995), where you can find solution for reading 2 types of geometry format. Changes are in:
a) ConstructGeometry(), where one can see the possibility of using "standard" geometry of FwEndCap of Emc, called module 3, and read together with other volumes of Emc from emc_module12345.dat and new version of FwEndCap of Emc read from emc_module3new.root;
b) ConstructRootGeometry() functions, where reading of the ROOT geometry is done.
Now, ROOT is exiting without any errors.
You can still use "standard" or "new" version of FwEndCap (of course not at the same moment!).
If somebody wants to use "standard" geometry, there is nothing changed, but if somebody wants to use geometry of Emc with NEW FwEndCap, in the sim_emc.C macro you should write:
PndEmc *Emc = new PndEmc("EMC",kTRUE);
Emc->SetGeometryFileNameDouble("emc_module1245.dat","emc_module3new.root ");
fRun->AddModule(Emc);
This changes were checked successfully at lxi008.gsi.de and at opteron cluster in Groningen and it works with G3, and _thanks_to_Mohammad_ with G4 as well.
If you have any problems, please let me know.
Best wishes,
Ola B.
Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
|
|
|
|
Re: FEEMC from ROOT file [message #6368 is a reply to message #6363] |
Fri, 11 April 2008 13:07 |
Aleksandra Biegun
Messages: 64 Registered: May 2007 Location: Groningen
|
continuous participant |
From: *KVI.nl
|
|
Hi Stefano,
now, both of these files with geometry are already at svn (emc_module3new.root you can also generate by yourself using createRootGeoFileFwEnd_changed.C, which is in the pandaroot/macro/emc/ directory).
If one wants to do simulation with both of the emc geometry, in the sim_emc.C macro one should use:
PndEmc *Emc = new PndEmc("EMC",kTRUE);
Emc->SetGeometryFileNameDouble("emc_module1245.dat","emc_module3new.root ");
fRun->AddModule(Emc);
With full_sim.C everything works fine for me, when I generated 1000 gammas with 1GeV or 10GeV energy, but when I tried to use each of the script separately (hit_emc.C, digi_emc.C and reco_emc.C), it crashed at digi level (there is still some strange detector ID produced in PndEmc.cxx, I am checking it).
Please run it to check if for you it works.
Cheers,
Ola.
Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
|
|
|
|
Re: FEEMC from ROOT file [message #6371 is a reply to message #6369] |
Fri, 11 April 2008 14:20 |
Aleksandra Biegun
Messages: 64 Registered: May 2007 Location: Groningen
|
continuous participant |
From: *KVI.nl
|
|
Hi,
have you used PndEmc:
PndEmc *Emc = new PndEmc("EMC",kTRUE);
instead of:
CbmDetector *Emc = new PndEmc("EMC",kTRUE);
?
Ola.
Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
|
|
|
|
Re: FEEMC from ROOT file [message #6373 is a reply to message #6372] |
Fri, 11 April 2008 14:36 |
Aleksandra Biegun
Messages: 64 Registered: May 2007 Location: Groningen
|
continuous participant |
From: *KVI.nl
|
|
Well...
It seams you do not have a ROOT file (emc_module3new.root)
in pandaroot/geometry/ directory, but I sent it and it's there (rev. 2475)
Have you downloaded it?
(for me it works!)
Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
|
|
|
|
|
|
Re: FEEMC from ROOT file [message #6377 is a reply to message #6376] |
Fri, 11 April 2008 15:14 |
Aleksandra Biegun
Messages: 64 Registered: May 2007 Location: Groningen
|
continuous participant |
From: *KVI.nl
|
|
I've done my sims at lxi009.gsi.de
with the new external package located at
/misc/cbmsoft/Debian3.1/mar08/fairsoft/
So, you do not have any warning when you use all macra, one by one, i.e.
sim_emc.C
hit_emc.C
digi_emc.C
reco_emc.C
or when you use
sim_emc.C + full_emc.C ?
Ola.
Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
|
|
|
|
Re: FEEMC from ROOT file [message #6631 is a reply to message #6595] |
Tue, 29 April 2008 12:33 |
StefanoSpataro
Messages: 2736 Registered: June 2005 Location: Torino
|
first-grade participant |
From: 140.181.10*
|
|
Is there a way to have the FWEEMC volumes with a different shape than TGeoArb8 ?
TGeoArb8 shape is not existing on Geant4, so if one tries to use Geant4 native navigation, the FWEEMC volumes are not recognised.
Is there a way to convert and to have them in the root file with a different shape, then?
[Updated on: Tue, 29 April 2008 12:33] Report message to a moderator
|
|
|
|
Re: FEEMC from ROOT file [message #6633 is a reply to message #6632] |
Tue, 29 April 2008 12:54 |
Aleksandra Biegun
Messages: 64 Registered: May 2007 Location: Groningen
|
continuous participant |
From: *KVI.nl
|
|
Hi,
I will try to change the shape and let you know...
Ola.
Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
|
|
|