Skip to content

CaloPID line improvements (Run1)

Maximilien Chefdeville requested to merge chefdevi_calib_improvements into stripping21-patches

Modification of 2 calo lines. Before/after rate comparison done without re-running the calo reco on 2012 data (i.e. commenting out #importOptions("$STRIPPINGSELECTIONSROOT/tests/users/DV-RedoCaloPID-Stripping21.py") in user test script to allow running over more stat in a shorter time. Test with re-running the calo-reco will follow, hence the WIP).

Ds* -> Ds gamma from StrippingCaloPID.py had too high timing due to use of StdLooseDplus2KKPi. Now use StdVeryTightDsplus2KKPi which implement P & PID cut on daughters. These cuts were present at a later stage anyway in my selections, so no signal is lost.

Rate test over 50k with StdLooseDplus2KKPi:

StrippingReport                                                INFO Event 50000, Good event 50000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |!StrippingCaloPIDDsst2DsGammaLine                            |  0.0280|        14|  1.143|  15.584|
 |!StrippingCaloPIDDsst2DsGammaLine_TIMING                     |  0.0280|        14|  1.143|   0.036|

Rate test over 50k with StdVeryTightDsplus2KKPi:

StrippingReport                                                INFO Event 50000, Good event 50000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |!StrippingCaloPIDDsst2DsGammaLine                            |  0.0280|        14|  1.143|   2.769|
 |!StrippingCaloPIDDsst2DsGammaLine_TIMING                     |  0.0280|        14|  1.143|   0.035|

-> factor of 5 improvement on timing for same physics output.

D0 -> K pi gamma from StrippingD02KPiGamma.py. This line was heavily pre-scaled (5%) due to the large retention rate (despite a large 2 GeV PT cut on the photon). To reduce the rate and remove the prescale, the D* -> D0 pi+ decay was reconstructed by adding a slow pion. A new stripping line is created: StrippingDst2D0Pi.py, the previous D02KPiGamma line can be ignored in the upcoming stripping campaigns.

Rate test over 50k without reco of the slow pion and 5% prescale:

StrippingReport                                                INFO Event 50000, Good event 50000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |_StrippingGlobal_                                            |  0.0300|        15|       |   9.801|
 |_StrippingSequenceStreamCalibration_                         |  0.0300|        15|       |   9.787|
 |!StrippingD02KPiGamma                                        |  0.0140|         7|  1.143|   1.679|
 |!StrippingD02KPiGamma_TIMING                                 |  0.0160|         8|  1.250|   1.687|

Rate test over 50k with reco of the slow pion and no prescale:

StrippingReport                                                INFO Event 50000, Good event 50000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |_StrippingGlobal_                                            |  0.0040|         2|       |  15.793|
 |_StrippingSequenceStreamCalibration_                         |  0.0040|         2|       |  15.777|
 |!StrippingDst2D0Pi                                           |  0.0040|         2|  1.000|   9.276|
 |!StrippingDst2D0P_TIMING                                     |  0.0040|         2|  1.000|   0.040|

-> factor of 10 improvement on retention rate at the cost of increased timing (not sure why). If necessary a prescale of 0.5 could be introduced.

Rate tests with re-running the calo on 2011 & 2012 data will follow.

Rate tests re-running the calo reconstruction.

Reminder: those lines are for calo PID: they go to full DST and save the calo raw banks. They can exceed the limit of 0.01% retention rate. For more precise rate, see the results previously posted (over 50k + calo repro should not change the rates significantly).

CaloPID 2011

StrippingReport                                                INFO Event 10000, Good event 10000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |_StrippingGlobal_                                            |  0.0100|         1|       |   3.147|
 |_StrippingSequenceStreamCalibration_                         |  0.0100|         1|       |   3.130|
 |!StrippingCaloPIDDsst2DsGammaLine                            |  0.0100|         1|  1.000|   0.891|
 |!StrippingCaloPIDDsst2DsGammaLine_TIMING                     |  0.0100|         1|  1.000|   0.035|

CaloPID 2012

StrippingReport                                                INFO Event 10000, Good event 10000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |_StrippingGlobal_                                            |  0.0500|         5|       |   6.402|
 |_StrippingSequenceStreamCalibration_                         |  0.0500|         5|       |   6.384|
 |!StrippingCaloPIDDs2EtapPiLine                               |  0.0400|         4|  1.000|   2.198|
 |!StrippingCaloPIDDsst2DsGammaLine                            |  0.0100|         1|  1.000|   2.949|
 |!StrippingCaloPIDDs2EtapPiLine_TIMING                        |  0.0400|         4|  1.000|   0.047|
 |!StrippingCaloPIDDsst2DsGammaLine_TIMING                     |  0.0100|         1|  1.000|   0.038|

Dst2D0Pi 2011

StrippingReport                                                INFO Event 10000, Good event 10000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |_StrippingGlobal_                                            |  0.0100|         1|	   |   2.406|
 |_StrippingSequenceStreamCalibration_                         |  0.0100|         1|       |   2.387|
 |!StrippingDst2D0Pi                                           |  0.0100|         1|  1.000|   1.736|
 |!StrippingDst2D0P_TIMING                                     |  0.0100|         1|  1.000|   0.041|

Dst2D0Pi 2012

StrippingReport                                                INFO Event 10000, Good event 10000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |_StrippingGlobal_                                            |  0.0100|         1|       |   3.530|
 |_StrippingSequenceStreamCalibration_                         |  0.0100|         1|       |   3.510|
 |!StrippingDst2D0Pi                                           |  0.0100|         1|  1.000|   2.755|
 |!StrippingDst2D0P_TIMING                                     |  0.0100|         1|  1.000|   0.043|

Merge request reports