bdrum / cern-physics

My analysis workflow
1 stars 1 forks source link

Comparison with Valeri data #15

Closed bdrum closed 4 years ago

bdrum commented 4 years ago

Now me and Valeri have discrepancies in the data.

He has less number of processed runs, but more statistics than me.

What is the reason of such discrepancies?

He'll provide to me file with event info and we can find out the difference...

bdrum commented 4 years ago

pt < 150MeV Valeri : 2317 vs my: 1807

There are list of runs, that doesn't contains events for me, but does for Valeri: {245705, 245952, 246042, 246115, 246178, 246182, 246276, 246495, 246759, 246808, 246945}

Valeri has 1408 unique events (I don't have them) I have 898 unique events (Valeri doesn't have them)

The accuracy of pt between intersection of events is 5.018930435141122e-07.

Run by Run comparison

download

bdrum commented 4 years ago

Comparison have a sense only for the same runs i.e. I have to check only Valeri's run list.

Also we've seen that my conditions for dca is too hard and seems I have to change it to just abs(dca[01])>3 continue;

bdrum commented 4 years ago

We have found out that the reason in the condition:

if (i >= 200) break;

where i is the number of track. It means that events with total track number greater than 200 I skipped that is wrong.

Now the condition looks like

if (nTracks >= 200) break;

i.e. in case of event has more than 200 good tracks I just break loop by tracks and save event.

I've ran scipts to the grid with such condition and share result when it will be dobe.

bdrum commented 4 years ago

After latest corrections:

Valeri good events number: 2317 Boris good events number: 2364

img