DHCAL Noise Filter
Embed Size (px)
description
Transcript of DHCAL Noise Filter

DHCAL Noise Filter

2
Analyze each ΔTS separately – assumption: noise is concentrated in 1TS Eliminate the noise TS and use the rest of the event
Cut if MaxBoardOcc>120
Board Noise
MC: DHCAL Oct10 setup, assume all hits in 1ΔTSe+: 2, 4, 6, 8, 10, 12, 16, 20, 25, 32 GeVpi+: 2, 4, 6, 8, 10, 12, 16, 20, 25, 32, 40, 50, 60 GeV
1
2 3
4
MC

Verify NoiseFilter1 2
In Time (using 2 ΔTS)Out Of Time

3 4
In Time (using 2 ΔTS)Out Of Time

MaxBoardFrameOcc>10 && MaxBoardFrameOcc==MaxBoardOcc next slide5
Board Frame Noise
Previously: Cut if MaxBoardFrameOcc>30Not so obvious now.
ASIC-related
For now, impose:Cut if (MaxBoardFrameOcc>10 && MaxBoardFrameOcc==CorrespondingBoardOcc)
Clear the ASIC-related noise first.

MaxBoardFrameOcc>10 && MaxBoardFrameOcc<20 && MaxBoardFrameOcc==MaxBoardOcc &&MaxBoardOcc<120
MaxBoardFrameOcc>20 && MaxBoardFrameOcc<30 && MaxBoardFrameOcc==MaxBoardOcc &&MaxBoardOcc<120
In Time (using 2 ΔTS)Out Of Time

In Time (using 2 ΔTS)Out Of Time
MaxBoardFrameOcc>30 && MaxBoardFrameOcc==MaxBoardOcc &&MaxBoardOcc<120

But, modify based on MC
Cut if (MaxBoardFrameOcc>15 && MaxBoardFrameOcc==CorrespondingBoardOcc)
MCMC - zoom

9
ASIC Noise
Cut if (MaxASICOcc≥45 && Nhits<300 or MaxASICOcc≥55 && Nhits≥300) && (ASICOcc≥45+0.8xMeanNeighLayerASICOcc)
Remark: No ASIC occupancy = 64 after MaxBoardOcc>120 cut – verified with 660210
MC
MC

10
After the cuts so far
20

11
21
22
23
Board borders aligned along the x-direction
Board borders aligned along the y-direction

20 21
Cut all (ASICOcc≥45+0.8xMeanNeighLayerASICOcc)

2223

ASIC Frame Noise
MC MC
Cut if (ASICFrameOcc>13 && ASICFrameOcc==CorrespondingASICOcc) (24)

24
Helper for the board frame noise which is not yet implemented.
Cut if (ASICFrameOcc>13 && ASICFrameOcc==CorrespondingASICOcc)

Board Frame Noise For the Cases 21, 22 and 23
21, 22, 23 all MC
MC
Implement two versions of cuts: • Loose: Require for each event, discard the event only if the noise is in time with the event.• Tight: Leave the implementation option to the analyzer. Cuts related to these plots on next slide
zoom

Loose Selection:
Noise Filter
Cut if MaxBoardOcc>120 or (MaxBoardFrameOcc>15 && MaxBoardFrameOcc==CorrespondingBoardOcc) or (ASICOcc≥45+0.8xMeanNeighLayerASICOcc) or (ASICFrameOcc>13 && ASICFrameOcc==CorrespondingASICOcc) or BoardBorderHits≥10 && MeanBorderLineNeigh≤1 for borders along X or BoardBorderHits≥10 && MeanBorderLineNeigh≤1 for borders along Y
Tight Selection:
Cut if MaxBoardOcc>120 or MaxBoardOcc>90 && Nhits<400 or (MaxBoardFrameOcc>15 && MaxBoardFrameOcc==CorrespondingBoardOcc) or (ASICOcc≥45+0.8xMeanNeighLayerASICOcc) or (ASICFrameOcc>13 && ASICFrameOcc==CorrespondingASICOcc) or BoardBorderHits≥10 for borders along X or BoardBorderHits≥15 for borders along Y or BoardBorderHits≥10 && MeanBorderLineNeigh≤1 for borders along Y

Noise Filter Statistics
Total number of events analyzed: 55600593 1
Events with bad timing (MPEvent≠MPRun2): 745428 (1.3%3)
Events after timing veto: 54855165
Events with in time noise4: 608909 (1.1%)
Events with out of time noise: 109040 (0.2%)
Events after noise veto: 54246256 (97.6% of grand total)
1 FNAL Runs ~ 23M, CERN Runs ~ 33M = ~ 56M.
2 The most popular (MP) ΔTS bin is not one of the two MPs of the run.
3 See the top left figure in slide 11 of ΔTS analysis.
4 In time with any of the selection patterns in slide 16 of ΔTS analysis.