Time Based Reconstruction

Transcription

1 Time Based Reconstruction Dominik Steinschaden

2 Content Introduction PndBranchBurstBuilder Time based reconstuction SttCellTrackFinder Status / Outlook 2/37

3 3/37

4 Perform online tracking and event building Software trigger... 4/37

5 Break continuois data stream into processable data packages Perform online tracking and event building Software trigger... 5/37

6 PndBranchBurstBuilder PndBranchBurstBuilder_eventBased PndBranchBurstBuilder_timeCut 6/37

7 class PndBranchBurstBuilder: public FairTask SetOutputPrefix(TString prefix) AddInputBranch(TString branchname) SetPersistence(Bool_t val) virtual InitStatus Init() virtual void Exec(Option_t* opt) (pure) virtual TClonesArray* GetBurstData(size_t inbranchnum) 7/37

8 PndBranchBurstBuilder PndBranchBurstBuilder_eventBased PndBranchBurstBuilder_timeCut 8/37

9 PndBranchBurstBuilder PndBranchBurstBuilder_eventBased Simulate burst data in event based simulation Combine defined number of events SetNEvents(int num) 9/37

10 PndBranchBurstBuilder PndBranchBurstBuilder_eventBased Simulate burst data in event based simulation Combine defined number of events SetNEvents(int num) /37

11 PndBranchBurstBuilder PndBranchBurstBuilder_eventBased Simulate burst data in event based simulation Combine defined number of events SetNEvents(int num) 11/37

12 PndBranchBurstBuilder PndBranchBurstBuilder_eventBased Not sorted in time! Simulate burst data in event based simulation Combine defined number of events SetNEvents(int num) 12/37

13 PndBranchBurstBuilder PndBranchBurstBuilder_timeCut 13/37

14 PndBranchBurstBuilder PndBranchBurstBuilder_timeCut Cut data stream periodically SetTimePeriod(int time) Only in timebased sim! Need sorted timestams! 14/37

15 PndBranchBurstBuilder PndBranchBurstBuilder_timeCut Cut data stream periodically SetTimePeriod(int time) Only in timebased sim! Need sorted timestams! HESR ~ 8% filled ~ 2 ns revolution time ~ 16 ns beam ~ 4 ns gap 15/37

16 PndBranchBurstBuilder PndBranchBurstBuilder_timeCut Cut data stream periodically SetTimePeriod(int time) Only in timebased sim! HESR ~ 8% filled ~ 2 ns revolution time ~ 16 ns beam ~ 4 ns gap Need sorted timestams! 16/37

17 PndBranchBurstBuilder PndBranchBurstBuilder_timeCut Cut data stream periodically SetTimePeriod(int time) Only in timebased sim! HESR ~ 8% filled ~ 2 ns revolution time ~ 16 ns beam ~ 4 ns gap Need sorted timestams! 17/37

18 PndBranchBurstBuilder PndBranchBurstBuilder_timeCut Cut data stream periodically SetTimePeriod(int time) Only in timebased sim! Need sorted timestams! 18/37

19 PndBranchBurstBuilder PndBranchBurstBuilder_eventBased PndBranchBurstBuilder_timeCut Upload to trunk ~ end of September Adjustment of Fairlinks missing (If needed) 19/37

20 Time based Reconstruction SttCellTrackFinder Status from PCM 2/217: Updated interface Works now with any input branch name PndSttHit You can specify an output prefix Works now also with bursts of events Evaluate the efficiency for time based data Work ongoing 2/37

21 Simulation Generator: DPM Momentum: 7 GeV/c events Digitization (event based) Digitization (time based) 2, 2 MHz event rate 16 ns beam 4 ns gap Reconstruction (SttCellTrackFinder) Event based Event based data bursts 1 6 events /burst Time based bursts 2, 2 MHz 21/37

22 Reco processing time processing time Burst reconstruction works! 16 Total processing time ~ burst size time [min] Burst size 22/37

23 Reco processing time processing time Burst reconstruction works! 16 Total processing time ~ burst size 14 Less packages More tracklets per package => more combinatorics 12 time [min] Burst size 23/37

24 Reco tracks and processing time processing time Burst reconstruction works! Total processing time ~ burst size Less packages More tracklets per package => more combinatorics Total # reconstructed tracks ~ burst size # MC Tracks constant => increased number of ghost tracks time [min] reconstructed tracks Burst size 24/37

25 Reco tracks and processing time processing time reconstructed tracks Burst reconstruction works! Total processing time ~ burst size Less packages More tracklets per package => more combinatorics Total # reconstructed tracks ~ burst size # MC Tracks constant => increased number of ghost tracks Event based = 1 event/burst reconstruction time [min] Burst size 25/37

26 Reco tracks and processing time processing time Burst reconstruction works! Total processing time ~ burst size Total # reconstructed tracks ~ burst size Event based = 1 event/burst reconstruction time [min] # reco tracks reconstructed tracks Burst size 26/37

27 Burst reconstruction works! Total processing time ~ burst size Total # reconstructed tracks ~ burst size Event based = 1 event/burst reconstruction pp ΛΛ pπ- pπ 27/37

28 Burst reconstruction works! Total processing time ~ burst size Total # reconstructed tracks ~ burst size Event based = 1 event/burst reconstruction Cern courier 31 may 212 Alice lead lead collision 28/37

29 Burst reconstruction works! Total processing time ~ burst size Total # reconstructed tracks ~ burst size Event based = 1 event/burst reconstruction Large burst size Less separated tracklets (seeds) Still high combinatoric (processing time) Less reconstructed tracks Cern courier 31 may 212 Alice lead lead collision 29/37

30 Reco tracks and processing time processing time Burst reconstruction works! Total processing time ~ burst size 16 Total # reconstructed tracks ~ burst size 5 Event based = 1 event/burst reconstruction 4 Large burst size Less separated tracklets (seeds) Still high combinatoric (processing time) Less reconstructed tracks # reco tracks 7 time [min] reconstructed tracks Time information not used Burst size 3/37

31 Reco tracks and processing time processing time Total processing time ~ burst size 16 Total # reconstructed tracks ~ burst size Event based = 1 event/burst reconstruction Large burst size 5 Less separated tracklets (seeds) Still high combinatoric (processing time) Less reconstructed tracks Time information not used Additional coordinate More separate tracklets Less combinatorics # reco tracks Burst reconstruction works! 7 time [min] reconstructed tracks Burst size 31/37

32 First full time based track reconstruction with pandaroot! Reco tracks and processing time processing time 7 Timebased reconstruction works! 2 MHz MHz 5 Total processing time ~ burst size Total # reconstructed tracks ~ burst size Event based = 1 event/burst reconstruction Time information not used Large burst size # reco tracks 2 6 time [min] reconstructed tracks Less separated tracklets (seeds) Still high combinatoric (processing time) Less reconstructed tracks Burst size 32/37

33 Reco tracks and processing time processing time Burst reconstruction works! 7 Timebased reconstruction works! 6 Total processing time ~ burst size 5 Total # reconstructed tracks ~ burst size Time information not used Loss of Ghosts or MC Tracks? # reco tracks 2 16 time [min] reconstructed tracks Burst size 33/37

34 Reco tracks and processing time processing time Burst reconstruction works! 7 Timebased reconstruction works! 6 Total processing time ~ burst size 5 Total # reconstructed tracks ~ burst size 4 Time information not used Loss of Ghosts or MC Tracks? We can not study that (now) # reco tracks 2 16 time [min] reconstructed tracks Burst size 34/37

35 Burst reconstruction works! Timebased reconstruction works! Tracking QA Task Total processing time ~ burst size Total # reconstructed tracks ~ burst size Works only for event based data! => we need a more general QA task Time information not used Loss of Ghosts or MC Tracks? Event based, burst, time based comparability We can not study that (now) 35/37

36 Time based digitization Detector MVD STT SciTil/Btof EMC Gem Ftof DIRC MDT FTS RICH DigiTask PndMvdDigiTask PndSttHitProducerRealFull PndSciTDigiTask Flight time Event time Time based FairLinks (time based) Pile up? PndGemDigitize PndFtofHitProducerIdeal???? 36/37

37 Time based digitization Time based digitization doesn t work out of the box Update time based example macros Evaluate status unclear detectors Detector MVD STT SciTil/Btof EMC Gem Ftof DIRC MDT FTS RICH DigiTask PndMvdDigiTask PndSttHitProducerRealFull PndSciTDigiTask Flight time Event time Time based FairLinks (time based) Pile up? PndGemDigitize PndFtofHitProducerIdeal???? 37/37