JeffersonLab / halld_recon

Reconstruction for the GlueX Detector
7 stars 9 forks source link

BCAL Showers Peaking at 212 cms when using RandomBKG #498

Open karthik18495 opened 3 years ago

karthik18495 commented 3 years ago

As per the discussions form the CalWG meeting on 3rd Feb 2021. This issue is being opened

I did look into sample (Run 50715) with and without Random Bkg. There seems to be a spike at 212 cms (middle of the BCAL) when using a random Bkg.

Output.pdf contains the following.

The first page has plots of a Photon gun shot at theta = 40 and momentum = 1GeV.

The second page of the plots has the same set of plots but for a bggen sample with and without the Random BKG. Output.pdf

zihlmann commented 3 years ago

are the noise hits at 212cm associated with the same readout channels?

markdalton commented 3 years ago

I asked for this to be made as an issue because I suspect that it is somehow related to the timing of the random background. In order for points in the BCAL to have the position 212 cm, the time of the hits on each end must be exactly the same. I don't know the details but I understand that the timing of the random events must be adjusted when they are merged with simulation events. It seems to me that this timing adjustment is not being done quite right for the BCAL hits.

Karthik, it might be useful to point to hddm files with and without the random background in order for people to look more closely.

zihlmann commented 3 years ago

maybe a 2d plot with phi (azimuth angle) vs z might give some insight.

karthik18495 commented 3 years ago

No, I believe it is not. Here you can see in the 2D plot (Phi Vs Z) for the bggen with noise. It doesnt happen at a single channel.

Also the folders that has the hddm files are at

/work/halld/ksuresh/zPi0 Each folder is named as (particle)(with/without)Noise

PhiVsZ.pdf

zisis0 commented 3 years ago

This can be investigated by Karthik. But it is probably not just a few channels.

On Feb 3, 2021, at 1:07 PM, zihlmann notifications@github.com wrote:

are the noise hits at 212cm associated with the same readout channels?

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/JeffersonLab/halld_recon/issues/498#issuecomment-772747142, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADANEOSD42DHWKA73RMNDIDS5GNGLANCNFSM4XBMEROA.

Dr. Zisis Papandreou Head, Department of Physics Professor of Physics, Ph.D., P.Phys.

Department of Physics
University of Regina
3737 Wascana Parkway Regina, SK S4S 0A2 CANADA

Phone: (306) 585-5379 Fax: (306) 585-5659 Email: zisis@uregina.ca mailto:zisis@uregina.ca Website: http://www.uregina.ca/science/physics/people/faculty-research/zisis-papandreou/index.html http://www.uregina.ca/science/physics/people/faculty-research/zisis-papandreou/index.html

zihlmann commented 3 years ago

ah, now I notice, 212cm is exaclty the center of the BCAL. and the other two smaller spikes are at the start and end of it.

sdobbs commented 3 years ago

Looking at one of these events, the hits that cause the shower in the middle don't seem to have exactly the same time, though it's close - I extracted a dump below.

We should check to see if these showers show up in the random trigger data themselves - maybe there really is some background in there? I pulled some files from tape, Karthik could you please run your code over them, when they show up in this location: `/cache/halld/RunPeriod-2018-08/recon/ver02/converted_random/merged


================================================================
Event: 27
DBCALCluster:
     r:   phi:      z: theta:      t:     E: E_preshower:  E_L2:  E_L3:  E_L4: N_cell: charge: t_rms:
-----------------------------------------------------------------------------------------------------
 68.95   1.58  144.84   0.71    3.31   1.06         0.05   0.64   0.31   0.05      17       0   0.39 
 68.30   4.60  232.39   0.39  -32.59   0.13         0.00   0.13   0.00   0.00       1       0   0.00 
 65.30   3.55  218.18   0.40  -44.48   0.12         0.12   0.00   0.00   0.00       1       0   0.00 

DBCALDigiHit:
 module: layer: sector:        end: pulse_integral: pulse_peak: pulse_time: pedestal: QF: nsamples_integral: nsamples_pedestal:
-------------------------------------------------------------------------------------------------------------------------------
      9      1       3  downstream            1035           0        5678         1   1                  1                  1 
     13      1       3    upstream            1404         109        1784         1   1                  1                  1 
     13      1       3  downstream            1080          84        1910         1   1                  1                  1 
     13      2       1    upstream             428          33        1775         1   1                  1                  1 
     13      2       1  downstream             460          35        1906         1   1                  1                  1 
     13      2       2    upstream            3332         260        1788         1   1                  1                  1 
     13      2       2  downstream            1920         150        1905         1   1                  1                  1 
     13      2       3    upstream           12328         963        1777         1   1                  1                  1 
     13      2       3  downstream            9096         710        1907         1   1                  1                  1 
     13      2       4    upstream             732          57        1782         1   1                  1                  1 
     13      2       4  downstream             720          56        1912         1   1                  1                  1 
     13      3       1    upstream             480          37        1781         1   1                  1                  1 
     13      3       1  downstream             316          24        1923         1   1                  1                  1 
     13      3       2    upstream            1980         154        1782         1   1                  1                  1 
     13      3       2  downstream            1968         153        1912         1   1                  1                  1 
     13      3       3    upstream            4464         348        1787         1   1                  1                  1 
     13      3       3  downstream            3360         262        1908         1   1                  1                  1 
     13      3       4    upstream             332          25        1787         1   1                  1                  1 
     13      3       4  downstream             496          38        1907         1   1                  1                  1 
     13      4       1    upstream             200          15        1810         1   1                  1                  1 
     13      4       1  downstream             148          11        1903         1   1                  1                  1 
     13      4       2    upstream             280          21        1806         1   1                  1                  1 
     13      4       2  downstream             576          45        1910         1   1                  1                  1 
     13      4       3    upstream             188          14        1792         1   1                  1                  1 
     13      4       3  downstream             252          19        1884         1   1                  1                  1 
     13      4       4    upstream             188          14        1781         1   1                  1                  1 
     13      4       4  downstream             232          18        1901         1   1                  1                  1 
     14      3       3    upstream             204          15        1801         1   1                  1                  1 
     14      3       3  downstream             284          22        1907         1   1                  1                  1 
     26      1       2    upstream            1352           0        5422         1   1                  1                  1 
     28      1       3    upstream            2836           0        1082         1   1                  1                  1 
     28      1       3  downstream            2807           0        1070         1   1                  1                  1 
     12      3       2    upstream             168          13        1806         1   1                  1                  1 
     12      3       2  downstream             116           9        1915         1   1                  1                  1 
     12      3       4    upstream              96           7        1792         1   1                  1                  1 
     12      3       4  downstream              84           6        1887         1   1                  1                  1 
     14      2       3    upstream             216          16        1789         1   1                  1                  1 
     14      2       3  downstream             116           9        1903         1   1                  1                  1 
     29      4       1    upstream             132          10        1873         1   1                  1                  1 
     36      2       3    upstream            2931           0        1290         1   1                  1                  1 
     36      2       3  downstream            2900           0        1251         1   1                  1                  1 
     44      1       1    upstream             939           0        5630         1   1                  1                  1 

DBCALHit:
 module: layer: sector:        end: E(GeV):  t(ns):
---------------------------------------------------
      9      1       3  downstream   0.031  254.88 
     13      1       3    upstream   0.043   11.50 
     13      1       3  downstream   0.033   19.38 
     13      2       1    upstream   0.013   10.94 
     13      2       1  downstream   0.014   19.12 
     13      2       2    upstream   0.101   11.75 
     13      2       2  downstream   0.058   19.06 
     13      2       3    upstream   0.379   11.06 
     13      2       3  downstream   0.278   19.19 
     13      2       4    upstream   0.022   11.38 
     13      2       4  downstream   0.022   19.50 
     13      3       1    upstream   0.015   11.31 
     13      3       1  downstream   0.010   20.19 
     13      3       2    upstream   0.060   11.38 
     13      3       2  downstream   0.060   19.50 
     13      3       3    upstream   0.136   11.69 
     13      3       3  downstream   0.102   19.25 
     13      3       4    upstream   0.010   11.69 
     13      3       4  downstream   0.015   19.19 
     13      4       1    upstream   0.006   13.12 
     13      4       1  downstream   0.004   18.94 
     13      4       2    upstream   0.008   12.88 
     13      4       2  downstream   0.017   19.38 
     13      4       3    upstream   0.006   12.00 
     13      4       3  downstream   0.008   17.75 
     13      4       4    upstream   0.006   11.31 
     13      4       4  downstream   0.007   18.81 
     14      3       3    upstream   0.006   12.56 
     14      3       3  downstream   0.009   19.19 
     26      1       2    upstream   0.041  238.88 
     28      1       3    upstream   0.086  -32.38 
     28      1       3  downstream   0.086  -33.12 
     12      3       2    upstream   0.005   12.88 
     12      3       2  downstream   0.003   19.69 
     12      3       4    upstream   0.003   12.00 
     12      3       4  downstream   0.003   17.94 
     14      2       3    upstream   0.007   11.81 
     14      2       3  downstream   0.003   18.94 
     29      4       1    upstream   0.004   17.06 
     36      2       3    upstream   0.089  -19.38 
     36      2       3  downstream   0.088  -21.81 
     44      1       1    upstream   0.028  251.88 

DBCALPoint:
 E(GeV): t(ns): z(cm): r(cm):   phi: module: layer: sector:
-----------------------------------------------------------
  0.006    4.6   88.7   73.3  1.424      12      3       2 
  0.004    3.3   96.1   73.3  1.489      12      3       4 
  0.054    3.8   81.9   65.3  1.587      13      1       3 
  0.020    3.4   78.4   68.3  1.522      13      2       1 
  0.113    3.8   85.7   68.3  1.554      13      2       2 
  0.471    3.5   79.0   68.3  1.587      13      2       3 
  0.032    3.8   79.0   68.3  1.620      13      2       4 
  0.017    4.1   71.2   73.3  1.522      13      3       1 
  0.088    3.8   77.6   73.3  1.554      13      3       2 
  0.171    3.8   82.3   73.3  1.587      13      3       3 
  0.019    3.8   82.9   73.3  1.620      13      3       4 
  0.008    4.4   95.9   81.2  1.522      13      4       1 
  0.020    4.5   89.9   81.2  1.554      13      4       2 
  0.010    3.2   96.4   81.2  1.587      13      4       3 
  0.009    3.4   81.2   81.2  1.620      13      4       4 
  0.007    3.7   87.3   68.3  1.718      14      2       3 
  0.011    4.2   90.3   73.3  1.718      14      3       3 
  0.125  -44.4  153.0   65.3  3.551      28      1       3 
  0.129  -32.2  167.2   68.3  4.598      36      2       3 

DBCALShower:
     E:      x:      y:     z:     t:     r:    phi: E_preshower:  E_L2:  E_L3:  E_L4: N_cell: Q:    dE:    dx:    dy:    dz:    dt: EXcorr: EYcorr: EZcorr: ETcorr: XYcorr: XZcorr: XTcorr: YZcorr: YTcorr: ZTcorr: sigLong: sigTrans: sigTheta: rmsTime:
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 1.065   -0.70   68.95  144.8    3.5   69.0   1.581        0.054  0.642  0.315  0.046      17  0  0.047  0.370  1.384   2.01  0.083  -0.004  -0.004  -0.005  -0.018  -0.030  -0.026  -0.026   0.849   0.681   0.565    1.673     0.008     0.011    0.390 
 0.133   -7.81  -67.85  232.4  -32.2   68.3  -1.685        0.000  0.129  0.000  0.000       1  0  0.019  0.633  1.237   4.28  0.180  -0.045  -0.196   0.150   0.079   0.168  -0.148  -0.071  -0.716  -0.380   0.354    2.873     0.009     0.009    0.000 
 0.129  -59.91  -25.97  218.2  -44.4   65.3  -2.733        0.125  0.000  0.000  0.000       1  0  0.019  1.196  0.770   4.25  0.182  -0.188  -0.135   0.141   0.081   0.504  -0.690  -0.363  -0.480  -0.256   0.351    2.874     0.009     0.008    0.000 
karthik18495 commented 3 years ago

Here I had run over 50k events over the file at /cache/halld/RunPeriod-2018-08/recon/ver02/converted_random/merged/converted_random_051302.hddm

The plots now show peaking at 2 places one at 212cm and another at around 175 cms OutputData.pdf

sdobbs commented 3 years ago

So, it looks like these features exist in the random triggered data that we are using to simulate backgrounds. We could look at data taken with other triggers, to see if it's there, maybe look at some other run periods. Does anyone have a hypothesis for why we might have these correlated hits? @markdalton ?

zihlmann commented 3 years ago

maybe?: random hit means equally possible hits for any time within the readout window. However hits close to either end of this window on both ends of the module will result in a time difference of either zero or +- full length of the module. Not sure how to go about this.

On 2/5/21 10:54 AM, Sean Dobbs wrote:

So, it looks like these features exist in the random triggered data that we are using to simulate backgrounds. We could look at data taken with other triggers, to see if it's there, maybe look at some other run periods. Does anyone have a hypothesis for why we might have these correlated hits? @markdalton https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_markdalton&d=DwMCaQ&c=CJqEzB1piLOyyvZjb8YUQw&r=Hy7ijcc6pcMoP-QxZxtQH4-vodW_VGkrA9xiBc7InXk&m=zIb2eOjCQMIHuuLlYtMc0KgukXTtMYWnVsIUTBZr6x0&s=CP02KFlii-lnRAT86wjNmFs_QDRqsDVh_RYIBEuRx5M&e= ?

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_JeffersonLab_halld-5Frecon_issues_498-23issuecomment-2D774118598&d=DwMCaQ&c=CJqEzB1piLOyyvZjb8YUQw&r=Hy7ijcc6pcMoP-QxZxtQH4-vodW_VGkrA9xiBc7InXk&m=zIb2eOjCQMIHuuLlYtMc0KgukXTtMYWnVsIUTBZr6x0&s=tZ4ipqLBhUZAh-n-hjxIL41l26vbg7ys_RDPqH2SwNg&e=, or unsubscribe https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_ADF7ACY3JNS2CXCW65XOSNDS5QIEDANCNFSM4XBMEROA&d=DwMCaQ&c=CJqEzB1piLOyyvZjb8YUQw&r=Hy7ijcc6pcMoP-QxZxtQH4-vodW_VGkrA9xiBc7InXk&m=zIb2eOjCQMIHuuLlYtMc0KgukXTtMYWnVsIUTBZr6x0&s=WENQkB4G_PtmbFqhbhyT5aUPRvG9xTAayCFtKTr58-w&e=.