2009 Cascade Filter Proposal

Size: px

Start display at page:

Download "2009 Cascade Filter Proposal"

Avis Page
5 years ago
Views:

1 2009 Cascade Filter Proposal Eike Middell and Sebastian Panknin DESY Zeuthen and Humboldt University of Berlin / Group Meeting ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 1 / 17

2 Outline 1 Filter Proposal Time at Desy Design of the filter Organisation at Desy Feature Extractor Settings 2 Cascade Filter 2009 Filter Idea Optimized Settings ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 2 / 17

The Pole Filter Script JEB Joint Event Builder Script Design at the Pole Events may be interesting for many filters and some filters will use similar variables: all modules run over all events a

3 The Pole Filter Script JEB Joint Event Builder Script Design at the Pole Events may be interesting for many filters and some filters will use similar variables: all modules run over all events a module may refuse a specific event (precondition not fulfilled) the filter modules tag the events, if they pass the filter events without any tag are rejected different filters can use same modules (in different ways) it is helpful to have the whole filter script ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 3 / 17

4 Requirements for a Proposal TFT Board Request All filter proposal have to fulfill certain standards: show data, MC agreement for IC40 show cut variables for IC60 decide between different settings and variables calculate performance for IC60 bandwidth estimation same (retriggered) datasets for most groups similar script to calculate variables some settings are the same (Feature Extractor) ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 4 / 17

5 Organisation Cooperation Let s work together, because: same data sets needed same scripts can be used easier to compare results The OFU people are also interested in a online filter study. Eike Sebastian Robert F. and Anna F. processing retriggering cut variables study cascade filter FE threshold OFU ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 5 / 17

6 Batch System About 500 CPUh were used for variable calculation: Resource Usage High Disk Memory High space consumption of 7 TB: for one setting: about 500 GB datasets: about 100 GB output files: hdf files: 700 MB root files: 2 GB i3 files: 97 GB Space will be given free starting next week. very full i3 files better save than sorry ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 6 / 17

7 Batch System About 500 CPUh were used for variable calculation: Resource Usage High Disk Memory High space consumption of 7 TB: for one setting: about 500 GB datasets: about 100 GB output files: hdf files: 700 MB root files: 2 GB i3 files: 97 GB Space will be given free starting next week. very full i3 files better save than sorry ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 6 / 17

8 Batch System About 500 CPUh were used for variable calculation: Resource Usage High Disk Memory High space consumption of 7 TB: for one setting: about 500 GB datasets: about 100 GB output files: hdf files: 700 MB root files: 2 GB i3 files: 97 GB Space will be given free starting next week. very full i3 files better save than sorry ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 6 / 17

9 Batch System About 500 CPUh were used for variable calculation: Resource Usage High Disk Memory High space consumption of 7 TB: for one setting: about 500 GB datasets: about 100 GB output files: hdf files: 700 MB root files: 2 GB i3 files: 97 GB Space will be given free starting next week. very full i3 files better save than sorry ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 6 / 17

10 Batch System About 500 CPUh were used for variable calculation: Resource Usage High Disk Memory High space consumption of 7 TB: for one setting: about 500 GB datasets: about 100 GB output files: hdf files: 700 MB root files: 2 GB i3 files: 97 GB Space will be given free starting next week. very full i3 files better save than sorry ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 6 / 17

11 Batch System About 500 CPUh were used for variable calculation: Resource Usage High Disk Memory High space consumption of 7 TB: for one setting: about 500 GB datasets: about 100 GB output files: hdf files: 700 MB root files: 2 GB i3 files: 97 GB Space will be given free starting next week. very full i3 files better save than sorry ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 6 / 17

12 Batch System About 500 CPUh were used for variable calculation: Resource Usage High Disk Memory High space consumption of 7 TB: for one setting: about 500 GB datasets: about 100 GB output files: hdf files: 700 MB root files: 2 GB i3 files: 97 GB Space will be given free starting next week. very full i3 files better save than sorry ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 6 / 17

13 Used Datasets Datasets Minimal triggered IC80 simulation (V ): dataset (IC40): single muon corsika: 1516 double muon corsika: 1511 triple muon corsika: 1550 ν e MC: 1568 ν µ MC: 1560 each 100 up to 200 files available in Zeuthen retriggered to IC40 and IC60 ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 7 / 17

14 What is the ADCThreshold threshold is used for charge calculation for multi peak reconstruction threshold of 1 corresponds to the hardware discriminator old setting was 0.4 new (proposed) setting is 1.0 A scratch of two pulses: the 0.25 PE pulse defines the 1.0 FE threshold. ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 8 / 17

15 Effect on Pulse Charge Charge per Pulse: only for the 1.0 threshold clear peak at 1. many pulses are split with lower threshold. ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 9 / 17

16 Effect on N sth number of pulses N hit : very sensitive on setting total charge of event Q event : quite robust number of channels N chan : per definition (number of DOM launches) independent number of modules with pulses N OM : a small dependence on thresholds ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 10 / 17

17 Effect on N sth number of pulses N hit : very sensitive on setting total charge of event Q event : quite robust number of channels N chan : per definition (number of DOM launches) independent number of modules with pulses N OM : a small dependence on thresholds ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 10 / 17

18 Effect on N sth number of pulses N hit : very sensitive on setting total charge of event Q event : quite robust number of channels N chan : per definition (number of DOM launches) independent number of modules with pulses N OM : a small dependence on thresholds ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 10 / 17

19 Effect on N sth number of pulses N hit : very sensitive on setting total charge of event Q event : quite robust number of channels N chan : per definition (number of DOM launches) independent number of modules with pulses N OM : a small dependence on thresholds ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 10 / 17

20 Cascade Filter 2009 ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 11 / 17

21 Overview Cascade Fiter Cascade Filter 2008: ToI EvalRatio > 0.12 and LineFit Velocity < 0.13 Sebastian's IC40 Test: Hz 17.7 Hz BG passing: 1.5%, nue E 2 passing: 70.3 % Cascade Filter 2009: 16 (+ 3) additional strings adjust filter to keep rate at about 17 Hz (eventually include some track fit zenith)

22 FeatureExtractor New Settings ADCThreshold (influence on NHits) NewDiscThreshold: True (use new routine to estimate discriminator threshold, data run has to be processed with False) LineFit: use only FirstLeadingEdge ToI: LineFit doesn't weight DOM with charge but with #pulses reduce influence of Nhits Loose Nchannel > 8 to consider StringTrigger events? Use coincified pulses in LineFit and ToI?

23 ADCThreshold no impact on nchannel small impact on charge huge influence on #pulses 1.0 setting physically preferred (histogram of pulse charges peaked at 1pe Sebastian's Wikipage) influence on final efficiency not totally understood

24 FilterVariables (2008 Filter Settings) ToI Evalratio Left: without DiscThreshold Right: with DiscThreshold

25 FilterVariables (2008 Filter Settings) Linefit Vel. Left: without DiscThreshold Right: with DiscThreshold

26 FilterVariables Gulliver TrackFit Zenith Left: without DiscThreshold Right: with DiscThreshold

27 Rates IC40 Last year: Hz Corsika: 1107 Hz Reco Success percentages (IC60) (preliminary): nue Data: 1076 Hz ToI: 78% IC60 Corsika: 1611 Hz Last year after filter: 17 Hz 2009 Goal: Suppress Background to 1 % Linefit 99 % Trackfit 78 % Corsika: ToI: 45 % Linefit: 97 % Trackfit: 53 % (Nchan > 10)

28 2009 Filter Optimization use max. LineFit velocity and min. ToI evalratio determine 1% contour line in corsika passing rate 2 plot on this search for maximum E passing rate along Nue E 2 passing corsika passing

29 IC60 - Old ToI and LineFit Settings (+coincify, +NewDiscThreshold, ADCTh 1.0) down to 1% background at 69% E 2 signal is feasible with old settings

30 New ToI and LineFit Settings (+coincify, +NewDiscThreshold, ADCTh 1.0) losen the nchannel cut in ToI more StringTrigger events to achieve the same BG supression only get 62.6% E 2 signal

31 New ToI and LineFit Settings (+no coincify, +NewDiscThreshold, ADCTh 1.0) Without coincify and with StringTrigger Events slightly better: 63.8% E 2 signal

32 Energy Efficiency Energy Dependence of the Cut: low efficiency below 1 TeV ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 12 / 17

33 Resulting numbers Optimal Cuts: This are the choosen cut values and the resultig efficiencies: LineFit vel ToI eval IC40(old) IC Resource requirements The following properties were calculated: bandwidth s = 6020 Mb d filter rate f = 16.1 Hz Sample Efficiency [%] BG MC 1 ν e E 1 89 ν e E 2 69 low calculation time (below 0.2 ms per event and module) ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 13 / 17

Summary 2009 Filter Proposals: good collaborations with OFU people high resource consumption Feature Extractor threshold of 1 physically preferred Cascade Filter LineFit and Tensor of Inertia as cut

34 Summary 2009 Filter Proposals: good collaborations with OFU people high resource consumption Feature Extractor threshold of 1 physically preferred Cascade Filter LineFit and Tensor of Inertia as cut variables 1 % background (16 Hz, 6 Gb d ) 69 % E 2 -signal efficiency Outlook fast clean up of unused data answering open questions on request ike Middell and Sebastian Panknin (DESY and HU) 2009 Cascade Filter Proposal Group Meeting 14 / 17

L2-Hitcleaning and other L2 tests. Matthias Danninger. - ed - RT -

Old FE NFE L2-Hitcleaning and other L2 tests - ed - RT - Overview NFE vs oldfe http://icecube.wisc.edu/~mda65/ic79/studies/nfe_fe_new.htm details about study different parameters (position & time & reconstruction)