Useful informations

Spectrometer settings

Usefuls runs

He 3 production

20 uA 1339 1340

10

40 uA 1343

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");

---

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>600")

---

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>600";

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");