Data production models for the CDF experiment

The data production for the CDF experiment is conducted on a large Linux PC farm designed to meet the needs of data collection at a maximum rate of 40 MByte/sec. We present two data production models that exploits advances in computing and communication technology. The first production farm is a centralized system that has achieved a stable data processing rate of approximately 2 TByte per day. The recently upgraded farm is migrated to the SAM (Sequential Access to data via Metadata) data handling system. The software and hardware of the CDF production farms has been successful in providing large computing and data throughput capacity to the experiment.Comment: 8 pages, 9 figures; presented at HPC Asia2005, Beijing, China, Nov 30 - Dec 3, 200

Novel statistical approaches for non-normal censored immunological data: analysis of cytokine and gene expression data

Background: For several immune-mediated diseases, immunological analysis will become more complex in the future with datasets in which cytokine and gene expression data play a major role. These data have certain characteristics that require sophisticated statistical analysis such as strategies for non-normal distribution and censoring. Additionally, complex and multiple immunological relationships need to be adjusted for potential confounding and interaction effects. Objective: We aimed to introduce and apply different methods for statistical analysis of non-normal censored cytokine and gene expression data. Furthermore, we assessed the performance and accuracy of a novel regression approach in order to allow adjusting for covariates and potential confounding. Methods: For non-normally distributed censored data traditional means such as the Kaplan-Meier method or the generalized Wilcoxon test are described. In order to adjust for covariates the novel approach named Tobit regression on ranks was introduced. Its performance and accuracy for analysis of non-normal censored cytokine/gene expression data was evaluated by a simulation study and a statistical experiment applying permutation and bootstrapping. Results: If adjustment for covariates is not necessary traditional statistical methods are adequate for non-normal censored data. Comparable with these and appropriate if additional adjustment is required, Tobit regression on ranks is a valid method. Its power, type-I error rate and accuracy were comparable to the classical Tobit regression. Conclusion: Non-normally distributed censored immunological data require appropriate statistical methods. Tobit regression on ranks meets these requirements and can be used for adjustment for covariates and potential confounding in large and complex immunological datasets

Applied immuno-epidemiological research: an approach for integrating existing knowledge into the statistical analysis of multiple immune markers.

BACKGROUND: Immunologists often measure several correlated immunological markers, such as concentrations of different cytokines produced by different immune cells and/or measured under different conditions, to draw insights from complex immunological mechanisms. Although there have been recent methodological efforts to improve the statistical analysis of immunological data, a framework is still needed for the simultaneous analysis of multiple, often correlated, immune markers. This framework would allow the immunologists' hypotheses about the underlying biological mechanisms to be integrated. RESULTS: We present an analytical approach for statistical analysis of correlated immune markers, such as those commonly collected in modern immuno-epidemiological studies. We demonstrate i) how to deal with interdependencies among multiple measurements of the same immune marker, ii) how to analyse association patterns among different markers, iii) how to aggregate different measures and/or markers to immunological summary scores, iv) how to model the inter-relationships among these scores, and v) how to use these scores in epidemiological association analyses. We illustrate the application of our approach to multiple cytokine measurements from 818 children enrolled in a large immuno-epidemiological study (SCAALA Salvador), which aimed to quantify the major immunological mechanisms underlying atopic diseases or asthma. We demonstrate how to aggregate systematically the information captured in multiple cytokine measurements to immunological summary scores aimed at reflecting the presumed underlying immunological mechanisms (Th1/Th2 balance and immune regulatory network). We show how these aggregated immune scores can be used as predictors in regression models with outcomes of immunological studies (e.g. specific IgE) and compare the results to those obtained by a traditional multivariate regression approach. CONCLUSION: The proposed analytical approach may be especially useful to quantify complex immune responses in immuno-epidemiological studies, where investigators examine the relationship among epidemiological patterns, immune response, and disease outcomes

Data processing model for the CDF experiment

The data processing model for the CDF experiment is described. Data processing reconstructs events from parallel data streams taken with different combinations of physics event triggers and further splits the events into datasets of specialized physics datasets. The design of the processing control system faces strict requirements on bookkeeping records, which trace the status of data files and event contents during processing and storage. The computing architecture was updated to meet the mass data flow of the Run II data collection, recently upgraded to a maximum rate of 40 MByte/sec. The data processing facility consists of a large cluster of Linux computers with data movement managed by the CDF data handling system to a multi-petaByte Enstore tape library. The latest processing cycle has achieved a stable speed of 35 MByte/sec (3 TByte/day). It can be readily scaled by increasing CPU and data-handling capacity as required.Comment: 12 pages, 10 figures, submitted to IEEE-TN

Direct Search for Charged Higgs Bosons in Decays of Top Quarks

We present a search for charged Higgs bosons in decays of pair-produced top quarks in pbar p collisions at sqrt(s) = 1.8 TeV using 62.2 pb^-1 of data recorded by the D0 detector at the Fermilab Tevatron collider. No evidence is found for signal, and we exclude at 95% confidence most regions of the (M higgs, tan beta) parameter space where the decay t->H b has a branching fraction greater than 0.36 and B(H -> tau nu) is large.Comment: 11 pages, 4 figures, submitted to Phys. Rev. Let

Limits on Anomalous WWgamma and WWZ Couplings

Limits on the anomalous WWgamma and WWZ couplings are presented from a simultaneous fit to the data samples of three gauge boson pair final states in pbar-p collisions at sqrt(s)=1.8 TeV: Wgamma production with the W boson decaying to enu or munu, W boson pair production with both of the W bosons decaying to enu or munu, and WW or WZ production with one W boson decaying to enu and the other W boson or the Z boson decaying to two jets. Assuming identical WWgamma and WWZ couplings, 95 % C.L. limits on the anomalous couplings of -0.30<Delta kappa<0.43 (lambda = 0) and -0.20<lambda<0.20 (Delta kappa = 0) are obtained using a form factor scale Lambda = 2.0 TeV. Limits found under other assumptions on the relationship between the WWgamma and WWZ couplings are also presented.Comment: 13 pages, 3 figures, submitted to Physical Review

Search for First Generation Scalar Leptoquark Pairs in ppbar Collisions at sqrt(s)=1.8 TeV

We have searched for first generation scalar leptoquark (LQ) pairs in the enu+jets channel using ppbar collider data (integrated luminosity= 115 pb^-1) collected by the DZero experiment at the Fermilab Tevatron during 1992-96. The analysis yields no candidate events. We combine the results with those from the ee+jets and nunu+jets channels to obtain 95% confidence level (CL) upper limits on the LQ pair production cross section as a function of mass and of beta, the branching fraction to a charged lepton. Comparing with the next-to-leading order theory, we set 95% CL lower limits on the LQ mass of 225, 204, and 79 GeV/c^2 for beta=1, 1/2, and 0, respectively.Comment: 14 pages, 4 figures, submitted to Physical Review Letters Replaced to correct visitor addresse

Second Generation Leptoquark Search in p\bar{p} Collisions at s\sqrt{s} = 1.8 TeV

We report on a search for second generation leptoquarks with the D\O\ detector at the Fermilab Tevatron $p\bar{p}$ collider at $\sqrt{s}$ = 1.8 TeV. This search is based on 12.7 pb$^{-1}$ of data. Second generation leptoquarks are assumed to be produced in pairs and to decay into a muon and quark with branching ratio $\beta$ or to neutrino and quark with branching ratio $(1-\beta)$. We obtain cross section times branching ratio limits as a function of leptoquark mass and set a lower limit on the leptoquark mass of 111 GeV/c$^{2}$ for $\beta = 1$ and 89 GeV/c$^{2}$ for $\beta = 0.5$ at the 95%\ confidence level.Comment: 18 pages, FERMILAB-PUB-95/185-

Search for High Mass Photon Pairs in p-pbar --> gamma-gamma-jet-jet Events at sqrt(s)=1.8 TeV

A search has been carried out for events in the channel p-barp --> gamma gamma jet jet. Such a signature can characterize the production of a non-standard Higgs boson together with a W or Z boson. We refer to this non-standard Higgs, having standard model couplings to vector bosons but no coupling to fermions, as a "bosonic Higgs." With the requirement of two high transverse energy photons and two jets, the diphoton mass (m(gamma gamma)) distribution is consistent with expected background. A 90(95)% C.L. upper limit on the cross section as a function of mass is calculated, ranging from 0.60(0.80) pb for m(gamma gamma) = 65 GeV/c^2 to 0.26(0.34) pb for m(gamma gamma) = 150 GeV/c^2, corresponding to a 95% C.L. lower limit on the mass of a bosonic Higgs of 78.5 GeV/c^2.Comment: 9 pages, 3 figures. Replacement has new H->gamma gamma branching ratios and corresponding new mass limit

Search for Top Squark Pair Production in the Dielectron Channel

This report describes the first search for top squark pair production in the channel stop_1 stopbar_1 -> b bbar chargino_1 chargino_1 -> ee+jets+MEt using 74.9 +- 8.9 pb^-1 of data collected using the D0 detector. A 95% confidence level upper limit on sigma*B is presented. The limit is above the theoretical expectation for sigma*B for this process, but does show the sensitivity of the current D0 data set to a particular topology for new physics.Comment: Five pages, including three figures, submitted to PRD Brief Report