Useful informations

From Hall A Wiki
Revision as of 00:19, 7 December 2006 by 129.57.164.26 (Talk) (d recoil vs phi target)

Jump to: navigation, search

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

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