Home » PANDA » PandaRoot » Bugs, Fixes, Releases » PndSttHitProducerRealFast() depends on event type?
| Re: PndSttHitProducerRealFast() depends on event type? [message #8675 is a reply to message #8622] |
Tue, 26 May 2009 17:03  |
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
|
first-grade participant |
From: *47-151.net24.it
|
|
Hi,
I ran your macros: they both run to the end without problems.
Did you see also some crash of the code or is it just a matter of the results?
Concerning the ouput, here is what I get:
Event Number 0
-I- PndSttTrackFinderIdeal: STTTrack 0 created from MCTrack 3 (3 STTPoints)
-I- PndSttTrackFinderIdeal: STTTrack 1 created from MCTrack 4 (3 STTPoints)
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 53
MCTracks: total 5, accepted 2, reconstructable: 2
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
-------------------------------------------------------
Event Number 1
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 5
MCTracks: total 6, accepted 0, reconstructable: 0
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
-------------------------------------------------------
Event Number 2
-I- PndSttTrackFinderIdeal: STTTrack 0 created from MCTrack 0 (3 STTPoints)
-I- PndSttTrackFinderIdeal: STTTrack 1 created from MCTrack 3 (3 STTPoints)
-I- PndSttTrackFinderIdeal: STTTrack 2 created from MCTrack 4 (3 STTPoints)
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 56
MCTracks: total 5, accepted 3, reconstructable: 3
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
-------------------------------------------------------
Event Number 3
-I- PndSttTrackFinderIdeal: STTTrack 0 created from MCTrack 0 (3 STTPoints)
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 26
MCTracks: total 3, accepted 1, reconstructable: 1
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
-------------------------------------------------------
Event Number 4
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 7
MCTracks: total 8, accepted 0, reconstructable: 0
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
-------------------------------------------------------
Event Number 5
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 1
MCTracks: total 3, accepted 0, reconstructable: 0
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
-------------------------------------------------------
Event Number 6
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 0
MCTracks: total 3, accepted 0, reconstructable: 0
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
-------------------------------------------------------
Event Number 7
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 0
MCTracks: total 3, accepted 0, reconstructable: 0
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
-------------------------------------------------------
Event Number 8
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 5
MCTracks: total 5, accepted 0, reconstructable: 0
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
-------------------------------------------------------
Event Number 9
-------------------------------------------------------
-I- Ideal STT track finding -I-
Hits: 1
MCTracks: total 3, accepted 0, reconstructable: 0
SttHits not found : 0
SttPoints not found : 0
MCTracks not found : 0
SttTracks not found : 0
So basically I see it can reconstruct: 2 tracks in evt 0, 3 tracks in evt 2 and 1 track in evt 3.
I checked the number of hits per track in the simulation stage and what I get is:
evt 0 track 4 hits 27
evt 0 track 3 hits 26
evt 1 track 0 hits 3
evt 1 track 5 hits 2
evt 2 track 4 hits 27
evt 2 track 3 hits 25
evt 2 track 0 hits 4
evt 3 track 0 hits 26
evt 4 track 7 hits 7
evt 5 track 0 hits 1
evt 8 track 4 hits 3
evt 8 track 3 hits 2
evt 9 track 0 hits 1
I looked into the PndSttTrackFinderIdeal and I found that there is actually a cut for tracks with less than 3 hits, so we expect the finder to reconstruct only the blue colored lines in the above list, which correspond to the events and number of tracks the track finder was actually able to reconstruct, except for the event number 4 (I have to understand it, because I see several hits in the same tube and this is strange, since we register an hit each time the particles leaves the volume, so a track should leave only one hit in each tube, with an entrance and en exit position... I will investigate this).
The problem I see is that there are too few hits in the stt to construct a track (even if we throw away the cut on the number of hits, the helix construction will fail) and I think this is due to the particular kind of events (this would also explain why everything runs fine with the box generator).
Can you try to generate particles only in the transverse plane for example, in such a way that they enter the stt and leave more hits? In this case everything should run fine.
Please let me know if there are still problems or if I didn' t answer properly to your questions.
Best regards,
Lia.
-
Attachment: numhits.C
(Size: 1.47KB, Downloaded 246 times)
|
|
|
|
 |
|
PndSttHitProducerRealFast() depends on event type?
By: donghee on Wed, 15 April 2009 16:23 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: Anonymous Poster on Wed, 15 April 2009 18:52
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: donghee on Thu, 16 April 2009 10:27 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: Anonymous Poster on Thu, 16 April 2009 12:30
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: donghee on Thu, 16 April 2009 14:20 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: donghee on Thu, 16 April 2009 14:53 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: donghee on Thu, 16 April 2009 16:19 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: donghee on Thu, 16 April 2009 23:12 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: donghee on Thu, 16 April 2009 23:39 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: asanchez on Fri, 17 April 2009 11:33 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: donghee on Thu, 16 April 2009 10:28 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
By: donghee on Sat, 23 May 2009 16:41 |
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
|
|
Re: PndSttHitProducerRealFast() depends on event type?
|
Goto Forum:
Current Time: Tue Sep 17 18:42:35 CEST 2024
Total time taken to generate the page: 0.00898 seconds
|
GSI Forum
Search
Help
Members
Login
Home
