Useful informations
Contents
Calibrations macros
In the onlana directory, there are a couple of script to generate histograms on particular detectors.
- DrawLS2.C
- DrawS2m.C
- DrawRS1.C
- DrawCer.C
this are scripts being worked on to plots and find the peak ofthe events. Most of the time you will have to edit the ranges to get sensible fit.
also available :
/adaqfs/home/adaq/e04018/onlana/Cali
- id_pedestals.C
- Determine_VDC_T0.C
Robert Feuerbach also has a couple of macros to calibrate the scintillator timing :
/adaqfs/home/adaq/e04018/rob
- time_s1_paddles.C
- time_s2_paddles.C
Usefuls runs
He 3 production
20 uA 1339 1340
10
40 uA 1343
He 3 production no pedestal subtraction run 1453 http://www.jlab.org/~adaq/halog/html/0612_archive/061211134332.html
How to replay data
Log in on any adaq computer as adaq. Type goonlana to go in the replay directory. Different kind of replays are available,
- replay.C usual replay it create a root file on /adaql3/work1/e04018/onlroot/onlana_XXXX.root
- replayFull.C same as replay.C but with output on /adaql10/work1/e04018/onlroot/onlan_XXXX.root for long replay when short replay for detector check is running
- replay_split.C takes all the split files and only decode the coincidence events with output on /adaql10/work1/e04018/onlroot/onlan_T5_XXXX.root
To replay just type : analyzer <name of script above> and answer the questions ( run number, target type and number of events to be replayed )
To look at predefined histograms after the replay is done. Use the script online.C
- For the standard replay
- .x online.C+("detectorL",runnumber);
- .x online.C+("detectorR",runnumber);
- .x online.C+("coin",runnumber);
- For the coincidence runs
- .x online.C+("detectorLT5",runnumber);
- .x online.C+("detectorRT5",runnumber);
- .x online.C+("coinT5",runnumber);
How to look at data
First open the file most useful one are the long replay on coincidence for example.
TFile * f=new TFile("/adaql10/work1/e04018/onlroot/onlana_T5_1365.root")
Then you can look at the tree ( ntuple ) named T.
To look at coincidence time :
T->Draw("D.ctimel","D.edtpl==0")
T->Draw("D.ctimel","D.edtpl==0&&L.cer.asum_c>600")
How to analyze many files. To chain many root files together one uses the TChain class :
ch = new TChain("T"); ch->Add("/adaql10/work1/e04018/onlroot/onlana_T5_1365.root"); ch->Add("/adaql10/work1/e04018/onlroot/onlana_T5_1366.root");
ch = new TChain("T"); ch->Add("/adaql10/work1/e04018/onlroot/onlana_T5_1346.root"); ch->Add("/adaql10/work1/e04018/onlroot/onlana_T5_1347.root");
T->Draw("D.ctimel","R.s1.la_p>400&&R.s1.la_p<800&&R.s1.lt>840&&R.s1.lt<1500") T->Draw("D.ctimel","R.s1.la_p>400&&R.s1.la_p<800&&R.s1.lt>840&&R.s1.lt<1500&&L.cer.asum_c>700")
Definition of the cut on s1 amplitude.
TCut s1cut = "R.s1.la_p>400&&R.s1.la_p<800&&R.s1.lt>840&&R.s1.lt<1500&&L.cer.asum_c>700";
How to print
gPad->Print("temp.ps");gSystem->Exec("pdq -Pcha2hp temp.ps");
d recoil vs phi target
T->Draw("R.gold.dp:R.tr.tg_ph","R.s1.la_p>400&&R.s1.la_p<800&&R.s1.lt>840&&R.s1.lt<1500&&L.cer.asum_c>600");
Kinematic run list
- ep calibration
1272
1352
- Kin 4
1353 1364
- Kin 5
1365 1375