PWG 3 analysis barrel https twiki cern chtwikibinviewALICEPWG
PWG 3 analysis (barrel) https: //twiki. cern. ch/twiki/bin/view/ALICE/PWG 3 Hadron https: //twiki. cern. ch/twiki/bin/view/ALICE/PWG 3 Electron Andrea Dainese (INFN Padova) thanks to R. Bala, C. Di Giglio, S. Masciocchi, R. Romita, A. Rossi, C. Zampolli ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 1
Contents Heavy-flavour vertexing software PWG 3/vertexing. HF Visualization in alieve Production and analysis of candidates on grid (train) Status of event mixing Heavy-flavour electron software PWG 3/hfe (S. Masciocchi) ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 2
D 0 Kpi in alieve D. Caffarri Visualise AOD objects in top of ESD display works only with one-to-one file corresponce between ESD and AOD will we have a AOD-standalone visualisation? Use similar objects/classes as for V 0 s Ali. Eve. HF, Ali. Eve. HFEditor, Ali. Eve. HFList. Editor ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 3
D. Caffarri Changing the cuts on the variable less candidates are drawn. The ones fitting the requirements are checked. 4
D. Caffarri Focusing on a candidate we can see the D 0 line in yellow, the positive and the negative tracks 5
D. Caffarri 2 D display with the position of the secondary vertex. 3 D rotaiting display with only the selected HF decay. Some details about D 0 features. 6
Heavy-flavour vertexing: production of candidates Ali. Analysis. Vertexing. HF: main class to produce the candidated charm decays: 2 -prong, 3 -prong, 4 -prong, cascades. In one go. It is the “core” of Ali. Analysis. Task. SEVertexing. HF wagon of the Official Analysis Train Output written to AOD event (friend tree in “delta-AOD” file Ali. AOD. Vertexing. HF. root) Configured by Ali. Analysis. Vertxing. HF* Config. Vertexing. HF. C News since last OW: self-configuration macro Add. Task. Vertexing. HF. C provided possibility to run it in the same aliroot session as rest of train (using the AOD event as input), thanks to Andrei possibility to run it on CAF (support for creation of AOD+delta. AOD), thanks to Andrei ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 7
HF candidates in Ali. AODEvent aod. Tree in Ali. AOD. root (standard AOD, unchanged ) tracks Ali. AODTrack. . . V 0 vertices (from ESD) friend aod. Tree in Ali. AOD. Vertexing. HF. root Ali. AODVertex. . . vertices. HF Ali. AODVertex D 0 to. Kpi Dstar Ali. AODReco. Decay. . . Access to daghters: Ali. AODReco. Cascade. HF *c = arr. Dstar->At(i); Ali. AODTrack *pion = c->Get. Bachelor(); Ali. AODReco. Decay. HF 2 Prong *dzero = c->Get 2 Prong() ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 8
Heavy-flavour vertexing: Analysis of candidates from AOD Input: Ali. AODs. root + Ali. AOD. Vertexing. HF. root Prepared “D 2 H” analysis train: Steering macro Run. Analysis. AODVertexing. HF. C 10 wagons (Task, Task. SE) up to now Ø Compare. HF (vertex resolutions) Ø Dplus analysis Ø Like-sign bkg D 0 Kpi Ø Like-sign bkg Jpsi ee Ø Analysis for Jpsi from B Ø D 0 Kpi CORRFW task Ø 4 x Promt D 0 fraction (with 4 sets of cuts) runs on the grid with Alien plugin Ø added support for delta. AODs to Alien plugin (create AOD chain and friend AOD chain) to be adapted for running on CAF ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 9
Heavy-flavour vertexing: Analysis of candidates from AOD Train run by Renu Bala on AODs from LHC 09 a 4 and LHC 09 a 5 from Official Train (Mihaela) Use par files from trunk because many tasks came in the past week and had to apply a patch to Ali. Analysis. Alien. cxx (see later) Input: Ali. AODs. root + Ali. AOD. Vertexing. HF. root The two files have to be read in parallel Each pair should be on the same SE? (split by SE, 10 files per subjob) Currently, we could analyse 17 M out of 48 M for LHC 09 a 4 and 4 M out of 6 M for LHC 09 a 5 this is mainly due to the “unzip error” (that we reported at the beginning of May), now fixed in root trunk caused missing Ali. AOD. Vertexing. HF. root loss of matching of the two input chains we added a temporary (private) patch to Ali. Analysis. Alien. cxx to check existence of both files in the catalogue before adding to the chain for LHC 09 a 5 the unzip error was checked in the validation of the Official Train no orphan AOD files, but still loss of statistics ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 10
Heavy-flavour vertexing: Analysis of candidates from AOD Issues/requests: size of Ali. AOD. Vertexing. HF. root was too large initially (< v 4 -17 -03): 47% of Ali. AOD. Reduced with tighter cuts, now it is 13% of Ali. AOD Official Train to be re-run to overcome the unzip error from v 4 -17 -01 (several changes in manager/handlers), error messages appear when creating or reading the AODs: Ø Error in <TRef. Table: : Add>: Set. Parent must be called before adding uid=131300 (Andreas investigating) Ø apparently harmless (still there) but. . . from same moment, looping on candidates in Ali. AOD. Vertexing. HF. root has slown down by huge factor (~100): really slow! unknown reason (correlated with errors above? ) having a large statistics allows to debug/tune our code Ø will need Official Train to be rerun in cycles Ø would be good to analyse our AODs on CAF (need to add support for reading of delta. AODs, already in Andrei’s to-do list) ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 11
Some plots from LHC 09 a 4/a 5 4 M evts LHC 09 a 5 D+ R. Bala ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 12
Some plots from LHC 09 a 4/a 5 17 M evts LHC 09 a 4 D 0 A. Rossi ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 13
Some plots from LHC 09 a 4/a 5 like-sign unlike-sign ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 17 M evts LHC 09 a 4 Jpsi bkg C. Di. Giglio 14
Some plots from LHC 09 a 4/a 5 4 M evts LHC 09 a 5 D 0 reco eff, C. Zampolli gen acc trk rec trk sel D 0 rec D 0 sel pt y ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 15
Preparation of event mixing R. Romita Goal: possibility to use mixed events to study background Need to mix 2, 3 or 4 events Will use Ali. Analysis. Task. ME Dedicated class Ali. Mixed. Event : Ali. VEvent implemented by Andreas Event selection multiplicity (OK) vertex position in z (OK) event shape/jettiness (under way) Computation of mean primary vertex and translation of tracks to this common vertex (needed for vertexing) Integration in common code for candidates production is under way ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 16
ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 17
ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 18
ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 19
ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 20
Summary Implemented visualization of D 0 Kpi Visualization from AOD-standalone? Good to have high-stat MC We are running our analysis train on the AODs several wagons, more coming, easy to add new wagons easy to run with alien plugin In coming weeks/months several cycles of run-check-fixrerun foreseen (for Official Train and our trains) on grid (OK) and CAF (coming) would be good to have fixed schedule for trunk tags “fixed” schedule for Official Train runs ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 21
EXTRA SLIDES ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 22
Ali. Analysis. Vertexing. HF class single-track cuts using Ali. ESDtrack. Cuts [2 sets of cuts: “displaced track”, “soft pion (D*)”, can be extended e. g. electrons] build (+, -) pairs and like-sign with “displ. tracks” and compute secondary vtx for D+, Ds+, Lc+ for D 0 Kp, J/y ee create Ali. AODReco. Decay. HF 2 Prong apply reco cuts for D* D 0 p special D 0 cuts, loop on “soft pions” for D* candidates store D 0 store D*+ store J/y ALICE Offline Week, CERN, 24. 06. 09 loop on all tracks (+ & -): build triplets, create Ali. AODReco. Decay. HF 3 Prong, apply reco cuts (common for the 3 particles? ) for D 0 Kppp loop on all tracks (+ & -) … create Ali. AODReco. Decay. HF … store D+, Ds+, Lc Andrea Dainese store D 0 23
Ali. VEvent as input Motivation: AOD as well as ESD AODs are smaller in size, will be replicated in more SEs Strategy followed: AOD well suited for “kinematics” analysis, not for vertexing ESD well suited for vertexing (Ali. External. Track. Param+Ali. ESDVertex) Use “virtual” interfaces to read ESD or AOD with the same code Øin case of AOD input, convert to External. Track. Param and ESDVertex Do vertexing as from ESD ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 24
Use virtual interfaces Virtual Ali. VEvent Ali. VParticle Ali. VTrack Ali. VVertex interfaces: ESD world AOD world Ali. ESDEvent Ali. AODEvent Ali. ESDVertex Ali. AODVertex Ali. ESDtrack Ali. AODTrack Possible issues: no common interface for single-track cuts (Ali. ESDtrack. Cuts used, but a ESD track created from a AOD tracks lacks some info, e. g. n. TPC clusters) PID (e. g. electron ID in TRD will need momentum at the TRD) cuts used to create AOD from ESD have to be known!!! stored in alien together with AOD? database? ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 25
Storing the candidates in the AOD: Ali. AODReco. Decay family Ali. VParticle Ali. VTrack lib. AOD NEW Ali. AODReco. Decay Ali. AODv 0 lib. AOD Ali. AODReco. Decay. HF lib. PWG 3 Ali. AODReco. Decay. HF 2 Prong Ali. AODReco. Decay. HF 3 Prong D 0, J/y from B, . . . D+, Ds, Lc, . . . lib. PWG 3 Ali. AODReco. Decay. HF 4 Prong lib. PWG 3 D 0 ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 26
Ali. AODReco. Decay from Ali. VTrack Ali. AODReco. Decay derived from Ali. VTrack can use it for further vertexing use cases: Ø heavy-flavour cascades (e. g. D* D 0 p, B D 0 h, B D 0 e Ø calculate impact parameters (with error) to primary vertex (feed-down from beauty) Ø use reconstructed V 0 s (D mesons? ) in primary vertex fit Ali. VTrack can be used as input to construct a track Ali. External. Track. Param t(v. Track); // charged track Ali. Neutral. Track. Param t(v. Track); // neutral track (e. g. D 0, V 0) Ø Ali. Neutral. Track. Param: public Ali. External. Track. Param with J. Belikov Get. C(Double_t bz) {return 0. ; } Charge() {return 0; } Ali. AODReco. Decay: : Get. Covariance. XYZPx. Py. Pz() xyz from the sec. vertex pxpypz as the sum of the daughter tracks’ cov matrices ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 27
Inclusions of D*+ D 0 p+ New class Ali. AODReco. Cascade. HF: public Ali. AODReco. Decay. HF 2 Prong no members, only methods use the TRef of Ali. AODVertex to point to the 2 Prong (D 0) and to the bachelor (ps) Ø use TRef * instead of TRef. Array because the latter requires all TRefs to have the Process. ID (thanks to Andreas) ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 28
New class Ali. AODReco. Cascade. HF Ali. VParticle lib. AOD Ali. VTrack Ali. AODReco. Decay Ali. AODv 0 lib. AOD Ali. AODReco. Decay. HF lib. PWG 3 Ali. AODReco. Decay. HF 2 Prong Ali. AODReco. Decay. HF 3 Prong D 0, J/y from B, . . . D+, Ds, Lc, . . . lib. PWG 3 Ali. AODReco. Decay. HF 4 Prong lib. PWG 3 D 0 Ali. AODReco. Cascade. HF D* ALICE Offline Week, CERN, 24. 06. 09 Andrea Dainese 29
- Slides: 29