ePIC pfRICH software mailing list

[Youqi Song <youqi.song AT yale.edu>]

Hi Kolja, all,

This indeed doesn't look trivial. I tried a different approach, where I printed the number of particles that fall into each bin (so not dividing by the total number of truth particles per row), with their sqrt(N) errors, like in the screenshot attached. I guess this can give us a rough estimate of the relative uncertainty of the LUT that I showed last Monday. Do you think this is fine? Or should I still try to look into rigorously propagating the errors with TEfficiency?

Best,

Youqi

On Mon, Jun 10, 2024 at 12:35 PM Kolja Kauder <kkauder AT gmail.com> wrote:

Thanks.

Well, it's not trivial. I turned on Sumw2, and I added "error" to the draw command, but the problem is that these histos are generated using

h[pbin][tbin][phibin]->SetBinContent(1, 4, fee[pbin][tbin][phibin]);

I.e., there's no error propagation.

Fixing that will require a bit of work and maths. I found this class which may help, https://root.cern.ch/doc/master/classTEfficiency.html 

Best,

Kolja

On Mon, Jun 10, 2024 at 12:12 PM Youqi Song <youqi.song AT yale.edu> wrote:

Hi Kolja,

I did p0t0phi0->Draw("colz, text")

Best,

Youqi

On Mon, Jun 10, 2024 at 11:45 AM Kolja Kauder <kkauder AT gmail.com> wrote:

Thanks, Youqi! I can run it, what draw commands did you use for the actual plots?

On Mon, Jun 10, 2024 at 11:29 AM Youqi Song <youqi.song AT yale.edu> wrote:

Hi Kolja, all,

I generated the ROOT file for spot checks with the same code as https://github.com/eic/pfRICH/blob/main/readTree.C. I ran the code with

 

root -l 'readTree.C("merged0606.root", "all_out_0606.root", "LUT_0606.txt")'

You can find the input and output files at /gpfs02/eic/youqi/pfRICH/standaloneSimu/pfRICH.

Best,

Youqi

On Wed, Jun 5, 2024 at 5:14 PM Youqi Song <youqi.song AT yale.edu> wrote:

Hi Kolja,

That makes sense. I will submit a batch of jobs to run overnight, so that we can look at the output tomorrow.

Best,

Youqi

On Wed, Jun 5, 2024 at 12:33 PM Kolja Kauder <kkauder AT gmail.com> wrote:

I mean, ultimately we want to regenerate the LUTs anyway, there's
nothing prohibiting us from just sending out the jobs for the full
set.

On Wed, Jun 5, 2024 at 9:47 AM Youqi Song <youqi.song AT yale.edu> wrote:
>
> Hi all,
>
> Perhaps I missed it during our Monday meeting, but what spot checks should we do exactly? Do we want to redo the exercise from the May campaign by running events throughout the acceptance, but with lower statistics, and then look at the migration matrices? Or do we want to generate events in some specific regions in theta and phi, and just look at the matrices for those specific regions?
>
> Best,
> Youqi
>
> On Tue, Jun 4, 2024 at 10:21 AM Kolja Kauder via Eic-projdet-pfrich-software-l <eic-projdet-pfrich-software-l AT lists.bnl.gov> wrote:
>>
>> Dear All,
>>
>> The new field map has been added to the standalone repository. To use it,
>> optionally switch a pre-existing repository to the new url if you
>> haven't done so.
>> git remote set-url origin git AT github.com:eic/pfRICH.git
>> git pull
>>
>> To use the new map, simply switch to the branch:
>> git checkout newmap
>>
>> I've also created a pull request but let's first run some checks
>> before merging it into main as the default.
>>
>> Thanks,
>> Kolja
>>
>>
>>
>> --
>> ________________________
>> Kolja Kauder, Ph.D.
>> NPPS, EIC
>> Brookhaven National Lab, Upton, NY
>> +1 (631) 344-5935
>> he/him/his
>> ________________________
>> --
>> Eic-projdet-pfrich-software-l mailing list
>> Eic-projdet-pfrich-software-l AT lists.bnl.gov
>> https://lists.bnl.gov/mailman/listinfo/eic-projdet-pfrich-software-l



--
________________________
Kolja Kauder, Ph.D.
NPPS, EIC
Brookhaven National Lab, Upton, NY
+1 (631) 344-5935
he/him/his
________________________

--

________________________
Kolja Kauder, Ph.D.
NPPS, EIC
Brookhaven National Lab, Upton, NY
+1 (631) 344-5935

he/him/his
________________________

--

________________________
Kolja Kauder, Ph.D.
NPPS, EIC
Brookhaven National Lab, Upton, NY
+1 (631) 344-5935

he/him/his
________________________