J/Psi picking
Goal is to pick those events containing a J/Psi to allow AA to be rerun
over the J/Psi sample to enhance the number of Ks, Lambda etc.
For p17 certification (contact Rick V.):
Summarized
location of raw J/psi samples (from Common Sample Group)
New
event/run lists for raw data skimming (including J/psi)
Data Samples:
set
1
set
2
set
3
set
4
set
5
set
6
set
7
set
8
set
9
set
10
set
11
set
12
set
13
set
14
set
15
set
16
set
17
set
18
set
19
set
20
set
21
set
22
set
23
set
24
Cuts used to find J/Psi
Dimuon trigger must fire
2 oppositely charged muons with nseg>0 and Pt>1.5
At least 2 muons with rank=1
For picking events dimuon mass between 2 and 4 GeV
(dimuon Mass between 2.6 and 3.5 for sets 11,12 ...)
J/Psi fit to Gaussian + linear background
Code can be found here
Results:
Plots of J/Psi mass:(Note for sets 1-4, multiple entries/event
are allowed, so the number of J/Psi in these events is an overestimate
due to duplicate muons. For set 5, only 1 entry/ event is allowed, so the
number of J/Psi in this sample will be underestimated slightly.
set
1 mass
set
2 mass
set
3 mass
set
4 mass
set
5 mass
set
6 mass
set
7 mass
set
8 mass
set
8 mass with nseg>=0
set
9 mass
set
10 mass
set
1 AA mass
set
2 AA mass
set
3 AA mass
set
4 AA mass
set
5n AA mass
set
11 mass
set
12 mass
set
13 mass
set
14 mass
set
15 mass
set
16 mass
set
17 mass
set
18 mass
set
19 mass
set
20 mass
set
21 mass
set
22 mass
set
23 mass
set
24 mass
Lists of run and event number:
set
1 run and event number
set
2 run and event number
set
3 run and event number
set
4 run and event number
set
5 run and event number
set
6 run and event number
set
7 run and event number
set
8 run and event number
set
9 run and event number
set
10 run and event number
set
11 run and event number
set
12 run and event number
set
13 run and event number
set
14 run and event number
set
15 run and event number
set
16 run and event number
set
17 run and event number
set
18 run and event number
set
19 run and event number
set
20 run and event number
set
21 run and event number
set
22 run and event number set
22 run and event number
with duplicates removed (use this one!!)
set
23 run and event number set
23 run and event number
with duplicates removed (use this one!!)
set
24 run and event number set
24 run and event number
with duplicates removed (use this one!!)
(instructions for removing duplicates if needed in the future:
/RunII/home/bayo/CHECKCOPY
--cDup is the executable file.
-- run it using ./cDup on the command line
-- the first prompt asks for the file used for checking
-- the second file prompts for the file you want to check
-- the third asks for the name of the file you want the output to go
into
-- the duplicate run and event numbers will be put in a file called
copy.out (so make sure there are isn't a file in your directory called
copy.out
Some statistics:
statistics
Instructions for picking raw events.
First prepare all of the necessary files for picking
** Updated June 12 2003 **
On d0mino make a directory for picking (i.e mkdir jpsi_picking_set5)
cd jpsi_picking_set5
cp /tmp_root/751/ok_2/babbott/jpsi_picking_set12/PickPrepare.py
cp /tmp_root/751/ok_2/babbott/jpsi_picking_set12/JobHelper3.py (JobHelper3.py
submits jobs every 5 minutes and resubmits jobs with less than 80 events)
copy a list of run and event numbers to your area. (i.e. cp run_event_numbers_set5_final.txt)
make the files to be run by typing python PickPrepare.py filename (i.e.
python PickPrepare.py run_event_numbers_set5_final.txt
Answer the questions: (for Example)
Input defname for this set :dimuon_pickevents-babbott-set5
Input groupname to work under :dzero
Input path for output (.(dot) for current path):.
Enter a base filename (ie. dimuons_set1.out :dimuons_set5.out
Enter a description for these events :Dimuon Events
This will prepare n files with 100 events each called event0.list, event1.list
... eventn.list, it will also create n files called run_pick_0.sh, run_pick_1.sh
... run_pick_n.list and a file called scripts which contains all of the
jobs you want to run
Next Submit the jobs for picking
setup D0RunII p14.03.01
setup sam
setup d0tools -t
python JobHelper3.py scripts
** when you run this program make sure you have no other pick_events job
running **
Here is a brief description of exactly what this does. For each job in
the list, the program submits ./run_pick_n.sh where run_pick_n.sh is the
shell script to run. It then records the jobID number and the proposed
output filename. The program then waits 10 minutes. Once the .out file
containing the raw data starts to be written a new job is submitted.
typically a job takes less than 1 hour to finish
do not delete the metadata files since they are needed to put the data
back into SAM
to see how many events were picked you can use: grep "Total events written"
reco*.out. A fully completed job should have 100 events written.
Location of Picked Events on Disk (scroll down for locations
in SAM):
set 1 : Rick Van Kooten : /prj_root/857/higgs_4/rickv/picked_p13_dimuons
and indiana-clued0:/rooms/living/rickv/pick_dimuons
set 2 : Burair sets 0-250 : /prj_root/877/ckm_1/burair/jpsi_picking_set2_new
set 2 : Daria sets 250-501 :/tmp_root/747/iu_1/daria/pick2
set 3 : Ike : in SAM (dataset definition bphysics-dimuons-picked-set3)
set 4 : Brad : In SAM (dataset definition bphysics-dimuons-picked-raw-set4)
set 5 : Heriberto/Pedro/Eduard : /prj_root/837/ckm_1/podesta2/pick_events_cinves_set5/
set 6 : Heriberto/Pedro/Eduard : /prj_root/837/ckm_1/podesta2/pick_events_cinves_set6/
set 7 : Paul : /work/diemen-clued0/balm/data/dimuons_set7
set 8 : Shashi: /prj_root/814/ckm_1/avdhesh/JPSIPICK/jpsi_picking_set8
set 9 : Ike: in SAM (dataset definition bphysics-dimuons-picked-set9)
set 10: Avdhesh: /prj_root/814/ckm_1/avdhesh/JPSIPICK/jpsi_picking_set10
set 10: Bayo: Second half: /prj_root/878/muid_2/bayo/set10_part2/
set 11: Bayo: /prj_root/878/muid_2/bayo/picked_set11/
set 12: Bayo/Brad: In SAM (dataset definition bphysics-dimuons-picked-raw-set12)
set 13: Bayo: /prj_root/878/muid_2/bayo/picked_set13/
set 14: Bayo: /prj_root/878/muid_2/bayo/picked_set14/
set 15: Abaz: /rooms/living/abaz/set15_raw/
set 16: Brad: In SAM (dataset definition bphysics-dimuons-picked-raw-set16)
set 17: Kin and Daria: in SAM (dataset definition bphysics-dimuons-picked-set17)
set 18: Jim Fast/Heriberto: First 300 files at: /prj_root/882/ckm_1/jfast/pickset18/
(in SAM)
set 18: Jim Fast/Heriberto: 300-399 at /prj_root/814/ckm_1/castilla/NPE_III/pe_File_*/
(in SAM)
set 18: Jim Fast/Heriberto: 400-635 at /prj_root/814/ckm_1/castilla/NPE_IV/pe_File_*/
(in SAM)
set 19: Jim Fast: /prj_root/882/ckm_1/jfast/pickset19/ (in SAM)
set 20: Jim Fast: /prj_root/882/ckm_1/jfast/pickset20/ (in SAM)
set 21: Jim Fast: /prj_root/882/ckm_1/jfast/pickset21/
set 22: Bayo, being picked into /tmp_root/747/iu_1/bayo/
set 23: Daria, being picked into /tmp_root/747/iu_1/daria/jpsi_picking_set23
set 24: RickV, being picked.into /prj_root/809/b_1/rickv1/jpsi_picking_set24
Post Fall (Nov. 2003) shutdown data:
Selection: 2MU skim (two loose muons), dimuon mass 2 - 4 GeV,
with obvious subdirectories
clued0:
/rooms/living/rickv/new_jspi_picking
/rooms/living/abaz/new_jpsi_picking
/rooms/foyer/pmcs/rickv/
For trigger studies (raw events, fine for reprocessing also)
indiana-clued0:/rooms/living/rickv/trigger_jpsis
Near Future Plans
Picking more than a million J/psi's is proving to
be unwieldy (plus we will also want skimmed raw or DST events for other
selections such as D mesons for reprocessing with htfaa-extended to increase
K_S and Lambda yields); therefore, work continues on implementing a general
system through the Common Sample Group of being able to skim DST events
from SAM in a systematic way using a given selection flag. The DST has
both RAW data and processed results (including the flag from the CSG indicating
the skim bit and classification of the event) that will allow this.
It is being worked on by Oleksiy Atramentov
(Iowa State), assisted by grad student Jason
Rieger (Indiana) as a general implementation. The first test case will
be operating over the AA or D skim, skimming out J/psi events identified
with the CS filter bit set: "SKIM_AA_JPSI" (J/psi & dimuon mass continuum
: selectJPsiFlag = true, see D-skim
documentation for details of selection).
Events Reprocessed with AA at:
/work/albion-clued0/bgv/data
TMB-trees of the events reprocessed with AA at:
/work/creek-clued0/tmb-trees/AA_dimuon/new_process/
TMB-trees of the events reprocessed with HTFAA at:
/rooms/living/chuluo/STREAM-Data/HtfAA/ (sets 1-10)
/rooms/living/abaz/tmb_dimuon_htfaa_set11
/rooms/living/abaz/tmb_dimuon_htfaa_set12
/rooms/living/abaz/tmb_dimuon_htfaa_set13
/rooms/living/abaz/tmb_dimuon_htfaa_set14
/rooms/living/abaz/tmb_dimuon_htfaa_set15
/rooms/living/abaz/tmb_dimuon_htfaa_set16
/rooms/living/abaz/tmb_dimuon_htfaa_set17
/rooms/living/abaz/tmb_dimuon_htfaa_set18
/rooms/living/abaz/tmb_dimuon_htfaa_set19
/rooms/living/abaz/tmb_dimuon_htfaa_set20
/rooms/living/abaz/tmb_dimuon_htfaa_set21
Storing the files in SAM
Kin has provided a description on how to put the raw files into SAM.
If you would like to free up disk space, you can put the files into SAM
and then delete them from disk. How
to put files into SAM
Some comments about storing files
It's not necessary that you do it on d0mino. You can do it also on flotsam-clued0.fnal.gov
(where SAM can be used properly)
The script does all of the renaming and storing so you only need to run
the one script
In set 4, the tier name is "reconstructed" in set 12, the tier name is
"reconstructed-bygroup". I am not sure when the tier name changed. Kin
wrote his script based on "reconstructed-bygroup" so if the tier name for
your set is "reconstructed" you will need to change the script slightly.
Dataset definitions available in SAM
set 3: bphysics-dimuons-picked-set3
set 4: bphysics-dimuons-picked-raw-set4
set 9: bphysics-dimuons-picked-set9
set 12: bphysics-dimuons-picked-raw-set12
set 16: bphysics-dimuons-picked-raw-set16
set 17: bphysics-dimuons-picked-set17
Last updated 19 July 2004, RvK