Simulation Reconstruction Study for Data Filtering

Comments/suggestions are appreciated!

Purpose of this page:

To find out filtering parameters for coming year (~10 strings).
Parameters which distinguish upgoing neutrino induced events from downgoing atmospheric background muon events.

Shown in this page:

1. 16 icetray MC sample information (NUSIM events and dCORSIKA background events)
2. hit distributions (NUSIM events and dCORSIKA background events)
3. Reconstruction results with first-guess methods (dipole, linefit and cfirst)
Analysis is done with IceRec current release (as of Aug14) V01-02-00.

Summary:

If you have no time to go through all just Jump to the first result plots
Observations:
1) Events look fine. But should we use anis instead of nusim? (I remember that was conclusion but not sure)
2) Hit distributions. To decide parameters for hit cleaning, is reasonable tray.SetParameter("clean", "Minimum Time", 16500.); tray.SetParameter("clean", "Maximum Time", 25000.);? (Not sure about parameters for conicidence distance and time, but need to decide very soon). No Amplitude nor TOT Cut.
3) Reconstruction. Linefit zenith looks good to distinguish atmospheric events from upgoing neutrino induced events. Not sure about the use of dipole fit for now (not so good reconstruction for neutrino events). Need more study on it.

We need cascade MC (where can we find?) to find out cfirst parameters
More fine tunining and statistics...

Comment:
I have encountered several segmentation breaks, assertion fails and PushBox fails! ...thus I had to modify several modules to run through.
If you want to see the errors I had, you can run for example testNUSIM.C on IceRec V01-02-00.

Relevant Parameters

Nusim comparison of zenith angle: mc truth and linefit results

Hit Cleaning and Coincidence Parameters

Cut 1

  tray.SetParameter("clean", "Output Name", "cleaned");
  tray.SetParameter("clean", "Minimum TOT", 125.);  //was 125
  tray.SetParameter("clean", "Maximum TOT", 2000.); //was 2000
  tray.SetParameter("clean", "Minimum Amplitude", 0.);//was 150
  tray.SetParameter("clean", "Maximum Amplitude", 200.); // was 400
  tray.SetParameter("clean", "Minimum Time", 16500.); // was 21000
  tray.SetParameter("clean", "Maximum Time", 25500.); //was 25500
  tray.SetParameter("iso", "Output Name", "coinc");
  tray.SetParameter("iso", "Input Selection", "cleaned");
  tray.SetParameter("iso", "Input Readout", "cleaned");
  tray.SetParameter("iso", "Isolation Time",     1000.);
  tray.SetParameter("iso", "Isolation Distance", 500.);

linefit parameter
  tray.SetParameter("linefit", "MinHits", 2);
dipole parameter
  tray.SetParameter("dipole",  "MinHits", 2);
cfirst parameter
  tray.SetParameter("c_first", "MinHits", 2);
  tray.SetParameter("c_first", "TriggerHits", 1);

Cut 2

  tray.SetParameter("clean", "Output Name", "cleaned");
  tray.SetParameter("clean", "Minimum TOT", 125.);  //was 125
  tray.SetParameter("clean", "Maximum TOT", 2000.); //was 2000
  tray.SetParameter("clean", "Minimum Amplitude", 0.);//was 150
  tray.SetParameter("clean", "Maximum Amplitude", 200.); // was 400
  tray.SetParameter("clean", "Minimum Time", 16500.); // was 21000
  tray.SetParameter("clean", "Maximum Time", 25500.); //was 25500
  tray.SetParameter("iso", "Isolation Time",     500.);
  tray.SetParameter("iso", "Isolation Distance", 70.);
linefit parameter
  tray.SetParameter("linefit", "MinHits", 2);
dipole parameter
  tray.SetParameter("dipole",  "MinHits", 2);
cfirst parameter
  tray.SetParameter("c_first", "MinHits", 2);
  tray.SetParameter("c_first", "TriggerHits", 1);

Cut 3

  tray.SetParameter("clean", "Output Name", "cleaned");
  tray.SetParameter("clean", "Minimum TOT", 125.);  //was 125
  tray.SetParameter("clean", "Maximum TOT", 2000.); //was 2000
  tray.SetParameter("clean", "Minimum Amplitude", 0.5);//was 150
  tray.SetParameter("clean", "Maximum Amplitude", 200.); // was 400
  tray.SetParameter("clean", "Minimum Time", 17000.); // was 21000
  tray.SetParameter("clean", "Maximum Time", 25000.); //was 25500
  tray.SetParameter("iso", "Isolation Time",     500.);
  tray.SetParameter("iso", "Isolation Distance", 30.);
linefit parameter
  tray.SetParameter("linefit", "MinHits", 4);
dipole parameter
  tray.SetParameter("dipole",  "MinHits", 4);
cfirst parameter
  tray.SetParameter("c_first", "MinHits", 5);
  tray.SetParameter("c_first", "TriggerHits", 4);

Muon and Cascade Reconstructions

other cfirst parameters

tray.SetParameter("c_first", "FirstLE", false);
tray.SetParameter("c_first", "AmpWeightPower", 1.0);
tray.SetParameter("c_first", "TriggerHits", 2);
tray.SetParameter("c_first", "DirectHitRadius", 100.0);
tray.SetParameter("c_first", "SmallClusterRadius", 30.0);
tray.SetParameter("c_first", "LargeClusterRadius", 60.0);
tray.SetParameter("c_first", "TriggerWindow", 200.0);
tray.SetParameter("c_first", "EarlyHitStart", -1500.0);
tray.SetParameter("c_first", "EarlyHitStop", -200.0);

dcorsika(black) and nusim(red) reconstruction with cut 3

Basic reconstruction information with Cut 3 dcorsika=background(black), nusim(red)

dcorsika(left) and nusim(right)

Basic reconstruction information with Cut 1 (dCORSIKA) input event # 2974, reconstructed event # 2974 In the lower right plot black is (reco zenith - true zenith) linefit and red is that of dipole fit	Basic reconstruction information with Cut 1(NUSIM) input event # 1509, reconstructed event # 1509 In the lower right plot black is (reco zenith - true zenith) linefit and red is that of dipole fit
Basic reconstruction information with Cut 2 (dCORSIKA) input event # 2974, reconstructed event # 2875 In the lower right plot black is (reco zenith - true zenith) linefit and red is that of dipole fit	Basic reconstruction information with Cut 2 (NUSIM) input event # 1509, reconstructed event # 1482 In the lower right plot black is (reco zenith - true zenith) linefit and red is that of dipole fit
Basic reconstruction information with Cut 3 (dCORSIKA) input event # 2974, reconstructed event # 1333 In the lower right plot black is (reco zenith - true zenith) linefit and red is that of dipole fit	Basic reconstruction information with Cut 3 (NUSIM) input event # 1509, reconstructed event #1070 In the lower right plot black is (reco zenith - true zenith) linefit and red is that of dipole fit