Difference between revisions of "Useful informations"

From Hall A Wiki
Jump to: navigation, search
 
(→‎How to replay data)
 
(22 intermediate revisions by 5 users not shown)
Line 1: Line 1:
−
[[http://www.jlab.org/~lerose/instructions/CSR_magnets_inst.htm]]
+
[http://www.jlab.org/~lerose/instructions/CSR_magnets_inst.htm Spectrometer settings]
 +
==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==
 
==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

Latest revision as of 15:52, 1 March 2007

Spectrometer settings

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
Retrieved from "https://hallaweb.jlab.org/wiki/index.php?title=Useful_informations&oldid=1584"