 |
|
|
| Home | Computing | DAQ | Notes | Meetings | Subsystems | Search |
|
|
| |||||||||||
|
This page last edited by PK on July 18, 2007 . |
This page explains how to analyse stripped BB events from the DC04v1 and v2 strippings. It assumes that you have already have designed a "final" selection (we'll call "your DC04 selection") using signal data, and that this selection selects a subset of the data accepted by a preselection run in the stripping (we'll call it "your preselection").
How to use the SETC or FETC in DaVinci: For quick access to your events
How to use the FETC in ROOT: To measure the downscaling factor
Statistics for DC04v2 (obsolete)
For the first round of the stripping, we run the preselections (and the TDR-selections too) using DaVinci v12r4 on BB events generated during DC04 v1. Remember that these events are affected by the random seed bug, which makes it difficult to use to assess B/S limits.
The second round called DC04v2 has been stripped using DaVinci v12r11 and there is no known random seeds problem. This is now obsolete. It has been replaced by the results of a third round done with DaVinci v12r16.
Practically a stripping job consists of three steps for each set of n DST files (n~80):
The selections are run on each of the n files, producing n temporary stripped DST and n temporary event tag collection files.
The n event tag collection files are merged into one, the FETC.
The n DSTs are merged by running DaVinci on them. Too noisy algorithm (h>1/1000) are down-scaled. This produces the output DST and another event tag collection, the SETC.
Warning:
The DC04v1 stripping was made using
this patch to correct the default
value of the ExclusiveSelection option in the
CombinedParticleMaker.
| Input DSTs: | The n DSTs used as input. |
| Stripped DST: |
The DST containing the events accepted by the preselections. |
| FETC: |
Full Event tag Collection. It is the result of the merging of the event tag collections. It contains the links to the original input DSTs used in the stripping. |
| SETC: | Stripped Event Tag Collection. The ETC file is the one produced in the same step as the output DST. It contains links to this DST only. |
|
|
An xml file that links the opaque ID that allows to link to the DST from any event in the ETC to the physical location of the file. |
To get stripped DSTs, go to the
bookkeeping database,
select DC04v1 data, inclusive bb, DST1 datatype and set "step 1" to "DaVinci
v12r4". Submit, click on the Gaudi logo and paste the output into your job
options (we recommend to replace all rfio: by castor:).
WARNING: Some of the stripping DSTs in the bookkeeping are corrupted (same events appears several times). They will be removed someday. For some other files, there is not the DST and all ETCs available. A coherent set of DC04v1 DSTs can be found here. It corresponds to 18060757 BB events.
The same applies for DC04v2. I maintain on my public area (~pkoppenb/public/stripping/DSC04v2.opts)
a set of DST files which have their corresponding ETC files. They are the only
ones for which one can determine the number of initial BB events.
Make sure that you run your DC04 selection only on events accepted by your preselection. To do this put your selection code in a GaudiSequencer like this:
ApplicationMgr.TopAlg += { "GaudiSequencer/AnaBu2LLKSeq"
};
AnaBu2LLKSeq.Members += { "CheckSelResult/CheckPreselBu2LLK"
, "PrintHeader/PrintOfficialPreselBu2LLK"
, "GaudiSequencer/SeqPreselBu2LLK" };
CheckPreselBu2LLK.Algorithms = { "PreselBu2LLK" };
where "PreselBu2LLK" has to be replaced by the name
of your preselection algorithm and "GaudiSequencer/SeqPreselBu2LLK"
would typically be your DC04 selection sequence. The
CheckSelResult
algorithm checks in the DST what decision has been made by your preselection in
this event. This is more than just a trick:
It ensures that you only consider events selected by your preselection,
For the noisy preselections, it also skips events that your preselection may have selected but which the downscaling decided not to keep. It may happen that such an event is saved by another (not downscaled) preselection. It is OK to look at these events, but be aware that they are very biased. Only very dirty (multi-collision for instance) events tend to be accepted by several preselections. You certainly don't want your efficiency to depend on others' preselections.
If for any reason you decide to re-run your preselection on
these events, just make sure
CheckSelResult
is executed before your preselection, or it would read the result of your
algorithm, not the stripping.
There is an example option file DVReadStripping.opts
in DaVinci v12r11 (re-tagged. You may need to cvs update).
How to access data quicker using the stripped or full ETC
This does not work with DC04v1. It should work with v2 but has not been tested. An attempt of how-to is here.
How to estimate your downscaling ratio using the full ETC
The full ETC allows to get back to the original input DSTs. It also allows to look at the efficiency of selections and their correlations simply in root. An analysis using this data has the advantage of allowing to access the events that have been discarded for too noisy algorithms. But this is strongly discouraged as it implies that you need to stage all available statistics. If you do this, the stripping was pointless.
Full ETC, when used as a simple root file, allow you to recover the efficiency of the downscaling for your preselection. If it matched the 1/1000 requirement safely (i.e. including statistical fluctuations in 40000/1000=40 events), there is no need to do this. Else you will have to open all FETC files in root and count the number of events that have been accepted by your preselection before the merging.
Fortunately, there is a small program (many thanks to Luis) that
does this for you. Presently you can find it on ~pkoppenb/public/stripping.
I will put it in a release soon. simply do
> source pkoppenb/public/stripping/setup_slc3_root40200.csh
> ~pkoppenb/public/stripping/stripping myETC.txt [ MyPresel ]
where myETC.txt is a file containing all ETC (FETC
or SETC) corresponding to your data sample. The format has to be the one generated by the bookkeeping. MyPresel is your
preselection. It is an optional argument. If you don't give it you'll get all
preselections including yours. There is more information in a README file in the
same directory. Note that the compiled version works for the DC04v2 stripping.
For DC04v1, you need to change a define in the code and recompile.
WARNING: There are some FETC in the bookkeeping that do not contain the expected number of events. You'd better avoid them. The full list of ETC that are OK can be downloaded from here: FETC. SETC. DST (DC04v1), FETC. SETC. DST. LFN. (DC04v2r3).
Measured Downscaling Ratios for DC04v1
Using this program and the list of files described above, I measured the downscaling factors listed below. You can safely use these factors for your preselection even if you did not run on the full sample.
Statistics based on 18060757 events in 464 files. Average = 38924 Selection Rate Saved Selected Lost Downscaling =============================================================================== PreselMuon 15.2379 o/oo. 18246 of 275209 (256963) ->15.0832 PreselBs2DsH 9.0050 o/oo. 18246 of 162638 (144392) -> 8.9136 PreselBd2RhoRho 6.2010 o/oo. 18246 of 111995 (93749) -> 6.1380 PreselBx2HHpi0Gamma 3.9646 o/oo. 18246 of 71605 (53359) -> 3.9244 PreselBx2nHGamma 3.5257 o/oo. 18246 of 63677 (45431) -> 3.4899 PreselBu2KsPi 2.2175 o/oo. 18243 of 40050 (21807) -> 2.1953 PreselBu2RhoRho0 1.8760 o/oo. 17990 of 33883 (15893) -> 1.8834 PreselB2DstarX2D02hh 1.6618 o/oo. 18242 of 30014 (11772) -> 1.6453 PreselBs2KsKs 1.5705 o/oo. 17527 of 28365 (10838) -> 1.6183 PreselBs2PhiPhi 1.5688 o/oo. 18241 of 28335 (10094) -> 1.5533 PreselBs2KstarKstar 1.4851 o/oo. 17749 of 26823 (9074) -> 1.5112 PreselBd2MuMuKstar 1.3485 o/oo. 18207 of 24356 (6149) -> 1.3377 PreselJpsi2MuMu 1.3117 o/oo. 18186 of 23692 (5506) -> 1.3027 PreselBd2PhiKs 1.1757 o/oo. 17979 of 21235 (3256) -> 1.1811 PreselBd2a1Dstar 1.1087 o/oo. 17795 of 20024 (2229) -> 1.1252 PreselD02MuMu 0.9747 o/oo. 16768 of 17604 (836) -> 1.0498 PreselB2hh 0.9118 o/oo. 15867 of 16468 (601) -> 1.0378 PreselBs2MuMu 0.8609 o/oo. 15343 of 15549 (206) -> 1.0134 PreselBd2Rho0Rho0 0.8328 o/oo. 14744 of 15042 (298) -> 1.0202 PreselJpsi2ee 0.7841 o/oo. 14062 of 14162 (100) -> 1.0071 PreselB2KstarGamma2ee 0.7405 o/oo. 13339 of 13375 (36) -> 1.0026 PreselBd2PiPiPi0 0.7295 o/oo. 13158 of 13177 (19) -> 1.0014 PreselBd2DstarPi 0.5939 o/oo. 10727 of 10728 (1) -> 1.0000 PreselBd2DstarPiInclusive 0.5483 o/oo. 9899 of 9903 (4) -> 1.0004 PreselBs2Phi2KK_Etac 0.5157 o/oo. 9314 of 9314 PreselBd2D02HH_Kstar2KPi 0.5031 o/oo. 9087 of 9088 (1) -> 1.0001 PreselB2PiPiPi02GammaGamma2ee 0.4111 o/oo. 7426 of 7426 PreselB2Pi0Pi02GammaGamma2ee 0.3737 o/oo. 6750 of 6750 PreselBs2DsDs 0.1782 o/oo. 3220 of 3220 PreselBu2LLK 0.1553 o/oo. 2805 of 2805 PreselB2Pi02GammaGamma2eePi02 0.0805 o/oo. 1454 of 1454
TDRselBd2Jpsi2MuMu_Ks2PiPi 0.0091 o/oo. 166 of 166 TDRselBd2DstarPi 0.0073 o/oo. 133 of 133 TDRselBs2DsPi 0.0046 o/oo. 84 of 84 TDRselBd2KPi 0.0042 o/oo. 76 of 76 TDRselBd2PiPi 0.0039 o/oo. 72 of 72 TDRselBd2Jpsi2ee_Ks2PiPi 0.0032 o/oo. 58 of 58 TDRselBs2PiK 0.0021 o/oo. 38 of 38 TDRselBd2D02KKKstar 0.0018 o/oo. 34 of 34 TDRselBs2KK 0.0011 o/oo. 21 of 21 TDRselBs2DsK 0.0011 o/oo. 21 of 21 TDRselBs2Phi2KK_Etac2PiPiKK 0.0010 o/oo. 19 of 19 TDRselBd2D02pipiKstar 0.0008 o/oo. 15 of 15 TDRselBd2D02KpiKstar 0.0007 o/oo. 13 of 13 TDRselBs2PhiGamma 0.0005 o/oo. 10 of 10 TDRselBd2KstarGamma 0.0005 o/oo. 10 of 10 TDRselBs2Phi2KK_Etac2PiPiPiPi 0.0003 o/oo. 6 of 6 TDRselBs2Phi2KK_Etac2KKKK 0.0002 o/oo. 4 of 4 TDRselBd2MuMuKstar 0.0002 o/oo. 4 of 4 TDRselBs2Jpsi2MuMu_Eta2Gamma 0.0001 o/oo. 3 of 3 TDRselBs2Jpsi2ee_Eta2GammaGa 0.0001 o/oo. 3 of 3 TDRselBs2PhiPhi 0.0001 o/oo. 2 of 2
Statistics for DC04v2 - OBSOLETE
Statistics based on 27291931 events in 764 files. Average = 35722 Total: 680528 events saved from 3325998 selected (122.86 o/oo) =================================================================================================== PreselTriggerSel: 27404 events saved from 2551081 selected (93.4738 o/oo) PreselMuon: 27404 events saved from 416636 selected (15.2659 o/oo) PreselBx2HHpi0Gamma: 27403 events saved from 102078 selected (3.74023 o/oo) PreselBx2nHGamma: 27402 events saved from 96841 selected (3.54834 o/oo) PreselBd2OmegaGamma: 27400 events saved from 77412 selected (2.83644 o/oo) PreselBd2KstarGamma: 27400 events saved from 67388 selected (2.46915 o/oo) PreselBu2KsPi: 27397 events saved from 63542 selected (2.32823 o/oo) PreselBs2PhiGamma: 27382 events saved from 48295 selected (1.76957 o/oo) PreselB2DstarX2D02hh: 27392 events saved from 46026 selected (1.68643 o/oo) PreselBs2PhiPhi: 27385 events saved from 42780 selected (1.5675 o/oo) PreselBc2RhoPi: 27379 events saved from 41151 selected (1.50781 o/oo) PreselBu2D0K: 27366 events saved from 38851 selected (1.42353 o/oo) PreselBd2MuMuKstar: 27326 events saved from 37039 selected (1.35714 o/oo) PreselJpsi2MuMu: 36291 events saved from 36291 selected (1.32973 o/oo) (not vetoed) PreselBu2RhoRho0: 25268 events saved from 26583 selected (0.974024 o/oo) PreselD02MuMu: 25139 events saved from 26440 selected (0.968785 o/oo) PreselBd2PhiKs: 23848 events saved from 24455 selected (0.896052 o/oo) PreselBs2DsH: 23249 events saved from 23679 selected (0.867619 o/oo) PreselBs2KsKs: 23231 events saved from 23626 selected (0.865677 o/oo) PreselB2hh: 23145 events saved from 23557 selected (0.863149 o/oo) PreselBs2MuMu: 23043 events saved from 23444 selected (0.859008 o/oo) PreselBd2RhoRho: 21945 events saved from 22171 selected (0.812365 o/oo) PreselJpsi2ee: 21366 events saved from 21522 selected (0.788585 o/oo) PreselBd2PiPiPi0: 21027 events saved from 21157 selected (0.775211 o/oo) PreselB2KstarGamma2ee: 20377 events saved from 20442 selected (0.749013 o/oo) PreselBs2KstarKstar: 20239 events saved from 20313 selected (0.744286 o/oo) PreselBd2DstarPi: 16316 events saved from 16318 selected (0.597906 o/oo) PreselBu2LLK: 15838 events saved from 15843 selected (0.580501 o/oo) PreselBs2PhiEtac: 15552 events saved from 15568 selected (0.570425 o/oo) PreselBd2DstarPiInclusive: 15071 events saved from 15073 selected (0.552288 o/oo) PreselBd2Rho0Rho0: 14323 events saved from 14325 selected (0.52488 o/oo) PreselBd2D02HH_Kstar2KPi: 13648 events saved from 13652 selected (0.500221 o/oo) PreselBd2DPi: 12551 events saved from 12553 selected (0.459953 o/oo) PreselB2PiPiPi02GammaGamma2ee: 11163 events saved from 11164 selected (0.409059 o/oo) PreselBd2a1Dstar: 10098 events saved from 10099 selected (0.370036 o/oo) PreselB2Pi0Pi02GammaGamma2ee: 9951 events saved from 9951 selected (0.364613 o/oo) PreselBs2DsDs: 4763 events saved from 4763 selected (0.17452 o/oo) PreselB2Pi02GammaGamma2eePi02GammaGamma 2057 events saved from 2057 selected (0.0753703 o/oo) TDRselBd2Jpsi2MuMu_Ks2PiPi: 244 events saved from 244 selected (0.00894037 o/oo) TDRselBd2DstarPi: 164 events saved from 164 selected (0.0060091 o/oo) TDRselBd2PiPi: 148 events saved from 148 selected (0.00542285 o/oo) TDRselBd2KPi: 137 events saved from 137 selected (0.0050198 o/oo) TDRselBs2DsPi: 103 events saved from 103 selected (0.00377401 o/oo) TDRselBd2Jpsi2ee_Ks2PiPi: 84 events saved from 84 selected (0.00307783 o/oo) TDRselBs2PiK: 64 events saved from 64 selected (0.00234502 o/oo) TDRselBd2D02KKKstar: 56 events saved from 56 selected (0.00205189 o/oo) TDRselBs2KK: 38 events saved from 38 selected (0.00139235 o/oo) TDRselBd2D02KpiKstar: 35 events saved from 35 selected (0.00128243 o/oo) TDRselBd2D02pipiKstar: 27 events saved from 27 selected (0.000989303 o/oo) TDRselBs2PhiGamma: 22 events saved from 22 selected (0.000806099 o/oo) TDRselBd2KstarGamma: 19 events saved from 19 selected (0.000696176 o/oo) TDRselBs2DsK: 17 events saved from 17 selected (0.000622895 o/oo) TDRselBs2PhiPhi: 12 events saved from 12 selected (0.00043969 o/oo) TDRselBd2MuMuKstar: 11 events saved from 11 selected (0.00040305 o/oo) TDRselBs2MuMu: 1 events saved from 1 selected (3.66409e-05 o/oo)
DC04v2r3 : 33 926 781 BB events analysed in 965 files. Average = 35157.3 Total: events saved from 3698307 selected (109.009 o/oo) =================================================================================================== PreselTriggerSel: 34300 events saved from 2802704 selected (82.6104 o/oo) PreselMuon: 34276 events saved from 490700 selected (14.4635 o/oo) PreselBx2HHpi0Gamma: 34259 events saved from 126385 selected (3.72523 o/oo) PreselBx2nHGamma: 34260 events saved from 120268 selected (3.54493 o/oo) PreselBd2OmegaGamma: 34254 events saved from 96088 selected (2.83222 o/oo) PreselBu2KsPi: 34231 events saved from 77639 selected (2.28843 o/oo) PreselHltDstarSel: 34231 events saved from 73790 selected (2.17498 o/oo) PreselB2DstarX2D02hh: 34225 events saved from 57035 selected (1.68112 o/oo) PreselBu2RhoRho0: 31431 events saved from 33133 selected (0.976603 o/oo) PreselBs2DsH: 29016 events saved from 29552 selected (0.871052 o/oo) PreselB2DiMuon: 28956 events saved from 29505 selected (0.869667 o/oo) PreselB2hh: 28694 events saved from 29264 selected (0.862563 o/oo) PreselBs2KsKs: 28671 events saved from 29231 selected (0.861591 o/oo) PreselB2DiElectron: 27776 events saved from 28115 selected (0.828696 o/oo) PreselBd2RhoRho: 27263 events saved from 27618 selected (0.814047 o/oo) PreselJpsi2MuMu: 27192 events saved from 27241 selected (0.802935 o/oo) PreselJpsi2ee: 26551 events saved from 26758 selected (0.788698 o/oo) PreselBd2PiPiPi0: 26167 events saved from 26348 selected (0.776614 o/oo) PreselLb2Lambda1670Gamma: 26049 events saved from 26243 selected (0.773519 o/oo) PreselBs2MuMu_II: 25433 events saved from 25607 selected (0.754772 o/oo) PreselB2KstarGamma2ee: 25400 events saved from 25513 selected (0.752002 o/oo) PreselLb2Lambda1115Gamma: 25294 events saved from 25444 selected (0.749968 o/oo) PreselBs2KstarKstar: 25062 events saved from 25206 selected (0.742953 o/oo) PreselBd2RhoGamma: 24662 events saved from 24794 selected (0.730809 o/oo) PreselBd2KstarGamma: 21260 events saved from 21314 selected (0.628235 o/oo) PreselB2KstarD0viaK0S: 20497 events saved from 20541 selected (0.605451 o/oo) PreselBd2DstarPi: 20262 events saved from 20309 selected (0.598613 o/oo) PreselBs2PhiEtac: 19115 events saved from 19160 selected (0.564746 o/oo) PreselBd2MuMuKstarNew: 18235 events saved from 18278 selected (0.538748 o/oo) PreselD02MuMu: 18117 events saved from 18150 selected (0.534976 o/oo) PreselBc2RhoPi: 18106 events saved from 18139 selected (0.534651 o/oo) PreselBd2Rho0Rho0: 17712 events saved from 17743 selected (0.522979 o/oo) PreselBd2DstarPiInclusive: 17532 events saved from 17561 selected (0.517615 o/oo) PreselBu2D0Knew: 17080 events saved from 17130 selected (0.504911 o/oo) PreselBu2LLK: 17044 events saved from 17083 selected (0.503526 o/oo) PreselBd2D02HH_Kstar2KPi: 16970 events saved from 17007 selected (0.501285 o/oo) PreselHltExclusiveSel: 16645 events saved from 16675 selected (0.4915 o/oo) PreselBd2DstarMu: 16232 events saved from 16264 selected (0.479385 o/oo) PreselBd2DPi: 15536 events saved from 15580 selected (0.459224 o/oo) PreselBs2MuMu: 15527 events saved from 15558 selected (0.458576 o/oo) PreselBs2PhiGamma: 15091 events saved from 15120 selected (0.445666 o/oo) PreselBd2a1Dstar: 14329 events saved from 14352 selected (0.423029 o/oo) PreselB2PiPiPi02GammaGamma2ee: 13801 events saved from 13833 selected (0.407731 o/oo) PreselBd2MuMuKstar: 13573 events saved from 13600 selected (0.400863 o/oo) PreselB2Pi0Pi02GammaGamma2ee: 12303 events saved from 12329 selected (0.3634 o/oo) PreselBs2DsRho: 12105 events saved from 12126 selected (0.357417 o/oo) PreselBu2D0K_KsPiPi: 10698 events saved from 10719 selected (0.315945 o/oo) PreselBs2PhiPhi: 10348 events saved from 10363 selected (0.305452 o/oo) PreselBs2DsHjb: 9953 events saved from 9967 selected (0.29378 o/oo) PreselBs2DsMuX: 7352 events saved from 7363 selected (0.217026 o/oo) PreselBd2EtapKs: 7173 events saved from 7184 selected (0.21175 o/oo) PreselBs2DsDs: 5974 events saved from 5988 selected (0.176498 o/oo) PreselB2DstarH2D02KsHH: 5675 events saved from 5688 selected (0.167655 o/oo) PreselBu2D0K_4hnew: 4554 events saved from 4560 selected (0.134407 o/oo) PreselBu2D0Mu: 4021 events saved from 4025 selected (0.118638 o/oo) PreselBu2D0K_KsKK: 3770 events saved from 3777 selected (0.111328 o/oo) PreselBd2KstarMuMu: 2676 events saved from 2684 selected (0.0791115 o/oo) PreselB2Pi02GammaGamma2eePi02GammaGamm: 2570 events saved from 2574 selected (0.0758693 o/oo) PreselBu2D0Pi: 2545 events saved from 2554 selected (0.0752798 o/oo) PreselBu2D0K: 2363 events saved from 2370 selected (0.0698563 o/oo) PreselBs2PhiMuMu: 2073 events saved from 2074 selected (0.0611316 o/oo) PreselBd2PhiKs: 1213 events saved from 1217 selected (0.0358714 o/oo) DC04selBu2JpsiK: 346 events saved from 348 selected (0.0102574 o/oo) TDRselBd2Jpsi2MuMu_Ks2PiPi: 242 events saved from 243 selected (0.00716248 o/oo) DC04selBu2LLK: 213 events saved from 214 selected (0.0063077 o/oo) TDRselBd2DstarPi: 213 events saved from 213 selected (0.00627823 o/oo) TDRselBd2PiPi: 197 events saved from 197 selected (0.00580662 o/oo) TDRselBd2Jpsi2MuMu_Kstar2KPi: 193 events saved from 193 selected (0.00568872 o/oo) TDRselBd2KPi: 178 events saved from 178 selected (0.00524659 o/oo) TDRselBs2DsPi: 130 events saved from 130 selected (0.00383178 o/oo) DC04selBu2D0Knew: 128 events saved from 128 selected (0.00377283 o/oo) TDRselBd2Jpsi2ee_Ks2PiPi: 123 events saved from 123 selected (0.00362545 o/oo) TDRselBs2PiK: 86 events saved from 86 selected (0.00253487 o/oo) TDRselBd2D02KKKstar: 70 events saved from 70 selected (0.00206327 o/oo) TDRselBs2KK: 52 events saved from 52 selected (0.00153271 o/oo) TDRselBd2D02KpiKstar: 41 events saved from 41 selected (0.00120848 o/oo) TDRselBd2D02pipiKstar: 33 events saved from 33 selected (0.000972683 o/oo) TDRselBd2KstarGamma: 25 events saved from 25 selected (0.000736881 o/oo) TDRselBs2Jpsi2ee_Phi2KK: 24 events saved from 24 selected (0.000707406 o/oo) TDRselBs2PhiGamma: 23 events saved from 23 selected (0.000677931 o/oo) TDRselBs2DsK: 21 events saved from 21 selected (0.00061898 o/oo) TDRselBs2PhiPhi: 15 events saved from 15 selected (0.000442129 o/oo) TDRselBd2MuMuKstar: 9 events saved from 9 selected (0.000265277 o/oo)
The full correlation table is here.
Any comment is very welcome. Flame me here.