DC2 Kickoff

• Full and realistic detector simulation, i.e., realistic in that simulations are imperfect, including things like dead strips, deadtime, etc.

• Full reconstruction chain from low level electronics quantities, to high level event information.

• Updated and detailed background model, including all the things that are not gamma-rays from celestial sources; also includes orbit variations.

• Rich description of the gamma-ray sky, including variable sources, pulsing sources and Gamma-Ray Bursts (GRBs).

• Revamped event classification and background rejection analysis, i.e., updated instrument performance plots.

This increased detail and improvements of the detector simulation, event classification analysis, and the sky results in a much broader range of analysis options and data types than have been made available in previous data challenges and science tools checkouts.

Links to Key Information: (linked from DC2 navbars???)

• Data Challenge II Workshop

• Getting Ready for DC2

• Science Tools Tutorials

• DC2 Users Forum (Confluence wiki pages)

• DC2 Analysis Results (Confluence wiki pages)

Note: If unable to see (or edit) a Confluence DC2 page, make sure you are logged in.

Data Challenge II Workshop
March 1-3, 2006
Stanford Linear Accelerator Center

Install DC2 Code Release

To install the current DC2 code release for your platform (Red Hat Enterprise Linux 3, Scientific Linux 4, Macintosh, or Windows), go to:

DC2 Code Releases

DATA PRODUCTS:

Need buttons for: FT1, FT2 (Pointing Livetime

For an introduction to the types of LAT data, see LAT Data Products and FITS Files by David Band.

Definition of LAT Event Summary File (FT1) (for DC2) -need to copy current?? (for archive)??? -David, Masa???

Glast Data Server's LAT Data - DC2 Summary Ntuple - same data different formats for those who want to see the full set of events or who want to view them in an event display or want to do their analysis in ROOT.

GLAST Science Support Center - to access the following "precut" FT1 data, go to
GSSC's DC2 Data Access site:

• LAT Photon data
• All sky data files (data from the entire sky in ~1 week chunks)
• LAT Pointing History and Livetime data file - (popup to anchored frame in David's ???doc. Seth???Julie ???)
• Livetime cubes
• GBM Burst data
• Source Catalog
• DC2 Pulsar Ephemeris
• Galactic Diffuse Emission Model

Note: You can access either TCut EventClass A or TCut EventClass B data from this site.

DC2 Data Servers

• Glast Data Server

Also see: Using the Glast Data Server

• Glast Science Support Center (GSSC) - Data Challenge 2 Data Access

??? DC2 Data Portal (Preliminary). Welcome to the DC2 data portal. DC2 data represents 55 days of simulated data. The easiest way to get started is to use the pre generated data sets.

DC2 Data Portal (Preliminary)

Welcome to the DC2 data portal. DC2 data represents 55 days of simulated data. The easiest way to get started is to use the pre generated data sets.

The following servers allow you to create your own custom datasets, with selected events, cuts, columns. They also allow you to access the full root-trees (digi, recon) and browse the data using the event display.

Astro Data Server (allows to select events by time, region in sky, etc).

Pruner (allows you to filter the Merit tuples by applying cuts, and selecting selected columns)

Peeler (allows you to select full root trees for specific events)

EVENT CLASSIFICATIONS (ask Bill to comment)???

Note: EventClassA and EventClassB are the Quality flags mentioned on the Data Server page. … users choose ‘A’ or ‘B’ (what are they? and what do I do with that knowledge?)

BACKTRACK on the following: http://www-glast.slac.stanford.edu/software/AnaGroup/atwood-2006Jan9-DC2.pdf

• TCut EventClassA = GoodEvent3;

• TCut EventClassB = GoodEvent1&&!GoodEvent3;

TCuts

You can download a copy of these TCut definitions from here (so no need to cut and paste from this page)

// required cuts for all event classes

TCut DC2Trigger="(GltWord&10)>0&&(GltWord!=35)";

TCut DC2Filter="FilterStatus_HI==0";

TCut DC2PrefilterCal="CalEnergyRaw>5&&CalCsIRLn>4";

TCut DC2AcdVeto="(AcdCornerDoca>-5&&AcdCornerDoca<50&&CTBTkrLATEdge<100)||((AcdActiveDist3D>0 || AcdRibbonActDist>0)&&Tkr1SSDVeto<2)";

// filter out high energy electrons (old definition)
// TCut DC2ElectronVeto="((min(abs(Tkr1XDir),abs(Tkr1YDir)) < .01 &&
// Tkr1DieEdge < 10 && AcdActiveDist3D > 0 && AcdActDistTileEnergy > .2) ||
//(Tkr1SSDVeto < 2 && AcdActiveDist3D > -3 && AcdActDistTileEnergy > .3) ||
//( AcdActiveDist3D >(-30 + 30*(Tkr1FirstLayer-2)))) && CTBBestEnergy > 5000";

// filter out high energy electrons
TCut DC2ElectronVeto="((min(abs(Tkr1XDir),abs(Tkr1YDir)) < .01 && Tkr1DieEdge < 10 && AcdActiveDist3D > 0 && AcdActDistTileEnergy > .2) || (Tkr1SSDVeto < 7 && AcdActiveDist3D > -3 && AcdActDistTileEnergy > .15) || ( AcdActiveDist3D >(-30 + 30*(Tkr1FirstLayer-2)))) && (CTBGAM+0.17*CTBBestLogEnergy)<1.75";

//filter out some events at low-med energy where the Track 2 starts higher up
//than Track 1.
TCut DC2AnotherVeto="(Tkr1FirstLayer - Tkr2FirstLayer) < 0 && Tkr2FirstLayer > 2 && Tkr2TkrHDoca>10 && (CTBGAM+0.16*CTBBestLogEnergy)<1.32 ";

//Heavy Ion Filter
TCut HeavyIonVeto = "CTBBestEnergy>1000 && (((CalTransRms-1.5)*Tkr1ToTTrAve)<5)&&CTBGAM>0.5";

//Anti-correlated filter
TCut AntiCorrVeto = "CTBBestEnergy<500&&((CalCsIRLn+2.5*Tkr1CoreHC/Tkr1Hits)<8 || (Tkr1CoreHC/Tkr1Hits)<0.03)";

//Cosmic proton filter
TCut ProtonVeto = "Tkr1FirstLayer<6&&AcdActiveDist3D>-80 && ((AcdActDistTileEnergy + AcdActiveDist3D/100)>1)";

//Global Ribbon Extension and AcdCornerDoca Extension
TCut GlobalRibbonVeto = "(AcdRibbonEnergy>0.2 && AcdRibbonActDist > -10) || (AcdCornerDoca >-5 && AcdCornerDoca<50 &&CTBTkrLATEdge<200)";

TCut DC2Vetos = DC2AcdVeto||DC2ElectronVeto||DC2AnotherVeto||HeavyIonVeto||AntiCorrVeto||ProtonVeto||GlobalRibbonVeto;

//------------------------------------------------------------------

//These are handy to use to make pruned tuples.

TCut Basic="CTBCORE>0.1&&CTBBestEnergyProb>0.1&&CTBGAM>0.";

TCut ratecut="CTBBestZDir<-0.3&&CTBBestEnergy>100";

//-------------------------------------------------------------

// Use one of the following to select an event class.

TCut DC2Base1="CTBCORE>0.1&&CTBBestEnergyProb>0.3&&CTBGAM>0.35";
TCut DC2Base2="CTBCORE>0.1&&CTBBestEnergyProb>0.1&&CTBGAM>0.55";
TCut DC2Base3="CTBCORE>0.35&&CTBBestEnergyProb>0.35&&CTBGAM>0.50";

// Final Analysis Classes

TCut GoodEvent1=(DC2Base1&&DC2Trigger&&DC2Filter&&DC2PrefilterCal)&&!DC2Vetos;
TCut GoodEvent2=(DC2Base2&&DC2Trigger&&DC2Filter&&DC2PrefilterCal)&&!DC2Vetos;
TCut GoodEvent3=(DC2Base3&&DC2Trigger&&DC2Filter&&DC2PrefilterCal)&&!DC2Vetos;

// For DC2 we propose using the GoodEvent1 and GoodEvent3 analysis classes.
// This allows the data to be divided into a set of nested event classes.

TCut EventClassA = GoodEvent3;
TCut EventClassB = GoodEvent1&&!GoodEvent3;

Trigger and OnboardFilter settings for DC2

The effect of the trigger and filter settings on trigger and downlink rates is discussed in this confluence space. For DC2 we have decided to use configuration (2) for the trigger setting. This requires either a tracker (3-in a row) or a cal hi trigger and also requires for Tkr triggers, that either the acd throttle bit is not set or cal-lo is set. For the onboard filter we require that FilterStatus_HI==0.

Useful References (???)

• High-Level Analysis of LAT Data: Use Cases - presentation by Seth Digel prepared for the SSC-LAT Working Group, June 2002

• DC1 Kickoff and DC1 Closeout, February 2004

• DC2 Software Workshop, June 2005

Meetings: (archived planning info???).

Mailing List; Confluence Forum;

General Information DC2 Planning Outline Steve Ritz's DC1-DC2 Plan Steering Committee Meetings Schedule	Workshops GSFC June 27-29, 2005 DC II Kickoff March 1-3, 2006
DC-2 Data Current merit ntuple (GR v7r0p2) Pipeline Status FTP access to MC runs GLAST data server	Communications DC-2 e-mail list archives (need GLAST user ID/password) Confluence pages
Documentation LAT SAS: Webpage, Workbook LAT Science Tools:Webpage, Workbook Merit Ntuple Description	Useful Information LAT Analysis/Calibration Group

Note: Test link: for page review purposes only.

http://glast-ground.slac.stanford.edu/workbook/science-tools/pages/DC2/dc2_Kickoff.htm