GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung
Home » PANDA » PandaRoot » Bugs, Fixes, Releases » [FIXED] Problem with FTS Mapper PndFtsMapCreator
[FIXED] Problem with FTS Mapper PndFtsMapCreator [message #15797] Mon, 20 January 2014 11:49 Go to next message
Denis Bertini is currently offline  Denis Bertini
Messages: 6
Registered: February 2004
Location: GSI
occasional visitor
 From: *gsi.de
Hi
I found some inconsistencies with the geometry information corresponding the
the straw tubes in the FS.
All FS straw tubes (PndFtsTube Objects ) are stored in a TClonesArray
and indexed vis the TubeID.

A dump of values shows that for the 2nd layer the tube position is wrong ( notably wrong Z ) :


-I- PndFtsClusterFinder:: ReadFtsHits()
-I- Digits Data: fHits[0] is branchID 5 name# FTSHit entries# 49
-I- PndFtsClusterFinder ChamberID# 1 LayerID# 1 TubeID# 51
-I- PndFtsClusterFinder Tube position X: 13.635 Y:0 Z:294.895 R_meas: 0.208595
-I- PndFtsClusterFinder wire direction
wx# 0 wy# 1 wz# 0
0
-I- PndFtsClusterFinder ChamberID# 1 LayerID# 2 TubeID# 192
-I- PndFtsClusterFinder Tube position X: -95.95 Y:0 Z:423.765 R_meas: 0.330289
-I- PndFtsClusterFinder wire direction
wx# 0 wy# 1 wz# 0

The Z position should be instead Z=295.75 ?

Any idea ?

[Updated on: Mon, 27 January 2014 09:26] by Moderator

Report message to a moderator

Re: Problem with FTS Mapper PndFtsMapCreator [message #15798 is a reply to message #15797] Mon, 20 January 2014 16:01 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant
From: *to.infn.it
The new mapper has some bug, Isabella is taking a look.

Report message to a moderator

Re: Problem with FTS Mapper PndFtsMapCreator [message #15801 is a reply to message #15797] Tue, 21 January 2014 16:09 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
 From: *physik.uni-giessen.de
I think it makes sense to post this information here in case other people will run into the same problem with the FTS digitization:

Previously, it was discovered that the digitization has problems for tracks which are close to the beam pipe.

As a temporary workaround one can use the MCPoint positions. To do that PndFtsHitProducerRealFast.cxx has to be modified like this:

1) add comment to the line 225: // TVector3 position = tube->GetPosition();
2) uncomment line 224: TVector3 position(point->GetX(), point->GetY(), point->GetZ()); //point info

Report message to a moderator

Re: Problem with FTS Mapper PndFtsMapCreator [message #15802 is a reply to message #15801] Tue, 21 January 2014 16:16 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
 From: *physik.uni-giessen.de
In a new post (so that it won't be overlooked) my new problem:

I get 
Info in <PndFtsMapCreator::PndFtsMapCreator>: Geometry -1 not supported by map

when I try to run my reco task in a macro. I have the following in the Init method of my class:

FairRuntimeDb *rtdb= FairRun::Instance()->GetRuntimeDb();
fFtsParameters=(PndGeoFtsPar*)(rtdb->getContainer("PndGeoFtsPar"));

// FTS Tube Array
PndFtsMapCreator *mapperFts = new PndFtsMapCreator(fFtsParameters);
fFtsTubeArray = mapperFts->FillTubeArray();


I have declared the variables as private as follows:
PndGeoFtsPar *fFtsParameters;
TClonesArray *fFtsTubeArray;


In my .h file I put
class TClonesArray;
class PndGeoFtsPar;

and in the .cxx 
#include "PndGeoFtsPar.h"
#include "PndFtsMapCreator.h"


Am I doing something wrong?

PandaRoot Revision: 23483, external packages apr13, Ubuntu 12.04 64 bit

[Updated on: Tue, 21 January 2014 16:17]

Report message to a moderator

Re: Problem with FTS Mapper PndFtsMapCreator [message #15803 is a reply to message #15797] Tue, 21 January 2014 16:26 Go to previous messageGo to next message
Denis Bertini is currently offline  Denis Bertini
Messages: 6
Registered: February 2004
Location: GSI
occasional visitor
 From: *gsi.de
You should implement the SetParContainer() virtual function of the FaitTasks:

void PndMyTask::SetParContainers() {

FairRuntimeDb* rtdb = FairRunAna::Instance()->GetRuntimeDb();
fFtsParameters = (PndGeoFtsPar*) rtdb->getContainer("PndGeoFtsPar");
}

in order to get the container at the right time in the initialisation.

Report message to a moderator

Re: Problem with FTS Mapper PndFtsMapCreator [message #15804 is a reply to message #15803] Tue, 21 January 2014 16:32 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
 From: *physik.uni-giessen.de
Thank you, that worked.

Report message to a moderator

Re: Problem with FTS Mapper PndFtsMapCreator [message #15810 is a reply to message #15804] Thu, 23 January 2014 17:02 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant
From: *to.infn.it
I have jsut committed an update in the trunk which should fix the mapper, thanks to Isabella.
Please update your fts and geometry folder, and try again. Hopefully the problem is fixed now.

Report message to a moderator

Re: Problem with FTS Mapper PndFtsMapCreator [message #15811 is a reply to message #15810] Fri, 24 January 2014 09:44 Go to previous messageGo to next message
Denis Bertini is currently offline  Denis Bertini
Messages: 6
Registered: February 2004
Location: GSI
occasional visitor
 From: *gsi.de
Unfortunately after updating both:

fts.geo ( geometry directory)

PndFtsMapCreator ( fts directory)

The issue with the tube mapping mentioned previously still
exists.
For my part, i solved the problem for the FTS digitizer by taking the tube
geometrical information directly from the TGeoManager nodes and i am
storing the tube geometry directly in the PndFtsHit structure ( for now )

When i run the FTS digitizer using both the new tube mapping ( with the updated
fts.geo ) i am getting now a systematic X coordinate mismatch i.e :
In the text dump that follows, one has to compare geometrical information

- from the MC points ( Geant point )
- from the mapper ( Tube from mapper )
- from the TGeoManager directly ( Tgeo Tube position ) :


-I- TubeId -------------> 5758
-I- Geant Point -------------> X Y Z 69.0942 : 11.3914 : 607.028
-I- Tube From Mapper -------------> X Y Z 75.245 : 0 : 606.995

-I- Tube from TGeo: Rmin Rmax dz 0 : 0.505 : 59
-I- TGeo Translation i= 0 : 69.185
-I- TGeo Translation i= 1 : 0
-I- TGeo Translation i= 2 : 606.995
-I- TGeo Rotation:
i= 0 : 1 i= 1 : 0 i= 2 : 0 i= 3 : 0 i= 4 : 0 i= 5 : 1 i= 6 : 0 i= 7 : 1 i= 8 : 0
-I- TGeo: Tube position X: 69.185 Y:0 Z:606.995 R_meas: 0.121848

-I- TubeId -------------> 6170
-I- Geant Point -------------> X Y Z 69.4256 : 11.412 : 607.964
-I- Tube From Mapper -------------> X Y Z 87.87 : 0 : 607.87

-I- Tube from TGeo: Rmin Rmax dz 0 : 0.505 : 59
-I- TGeo Translation i= 0 : 69.69
-I- TGeo Translation i= 1 : 0
-I- TGeo Translation i= 2 : 607.87
-I- TGeo Rotation:
i= 0 : 1 i= 1 : 0 i= 2 : 0 i= 3 : 0 i= 4 : 0 i= 5 : 1 i= 6 : 0 i= 7 : 1 i= 8 : 0
-I- TGeo: Tube position X: 69.69 Y:0 Z:607.87 R_meas: 0.277454

-I- TubeId -------------> 6581
-I- Geant Point -------------> X Y Z 70.6534 : 11.4883 : 611.426
-I- Tube From Mapper -------------> X Y Z 101.893 : 0 : 611.49

-I- Tube from TGeo: Rmin Rmax dz 0 : 0.505 : 59
-I- TGeo Translation i= 0 : 71.477
-I- TGeo Translation i= 1 : 0
-I- TGeo Translation i= 2 : 611.49
-I- TGeo Rotation:
i= 0 : 0.996195 i= 1 : 0 i= 2 : -0.0871557 i= 3 : 0.0871557 i= 4 : 0 i= 5 : 0.996195 i= 6 : 0 i= 7 : 1 i= 8 : 0
-I- TGeo: Tube position X: 71.477 Y:0 Z:611.49 R_meas: 0.190151

-I- TubeId -------------> 6992
-I- Geant Point -------------> X Y Z 70.986 : 11.5091 : 612.362
-I- Tube From Mapper -------------> X Y Z 114.566 : 0 : 612.365

-I- Tube from TGeo: Rmin Rmax dz 0 : 0.505 : 59
-I- TGeo Translation i= 0 : 71.9839
-I- TGeo Translation i= 1 : 0
-I- TGeo Translation i= 2 : 612.365
-I- TGeo Rotation:
i= 0 : 0.996195 i= 1 : 0 i= 2 : -0.0871557 i= 3 : 0.0871557 i= 4 : 0 i= 5 : 0.996195 i= 6 : 0 i= 7 : 1 i= 8 : 0
-I- TGeo: Tube position X: 71.9839 Y:0 Z:612.365 R_meas: 0.0175759

Report message to a moderator

Re: Problem with FTS Mapper PndFtsMapCreator [message #15812 is a reply to message #15811] Fri, 24 January 2014 18:39 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant
From: *to.infn.it
Dear Denis,
maybe it would be better to copy here the code you are using, since Isabella (who cannot read the forum) is not able to reproduce your problem. Just send your macros, sim+digi+check.

By the way, could it be that you forgot to remove the parameter file after updating the geometry, and that you had an old geometry definition there screwing the mapper?

Report message to a moderator

Re: Problem with FTS Mapper PndFtsMapCreator [message #15814 is a reply to message #15812] Mon, 27 January 2014 08:28 Go to previous message
Denis Bertini is currently offline  Denis Bertini
Messages: 6
Registered: February 2004
Location: GSI
occasional visitor
 From: *gsi.de
Hi Stefano,
You were right, the new fts geometry was not in sync with the stored
params.root.
Updating the parameter file solved then the problem.
Thanks !

Report message to a moderator

Previous Topic: [FIXED] TGeant4 Error with pbarpSystem
Next Topic: [FIXED] Double tracks/double MC match issue
Goto Forum:
  

[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Sat Dec 07 01:07:53 CET 2024

Total time taken to generate the page: 0.00568 seconds
.:: Contact :: Home :: Legal Notice :: Data Protection Policy ::.

Powered by: FUDforum 3.1.3.
Copyright ©2001-2023 FUDforum Bulletin Board Software