Difference between revisions of "Some programs"

From New IAC Wiki
Jump to navigation Jump to search
(Replaced content with '[http://wiki.iac.isu.edu/index.php/CODA go back] ==some Alex program from Sep 2011 run== //to see spectrum for all 9 detectors * plot22.C ** root *.root ** .x plot22.C …')
 
Line 20: Line 20:
 
** root -l NDpairStudy.C'()'
 
** root -l NDpairStudy.C'()'
 
** .x NDpairStudy.C()
 
** .x NDpairStudy.C()
−
 
−
 
−
{
 
−
gROOT->Reset();
 
−
gStyle->SetOptStat(1111111);
 
−
gStyle->SetOptDate(1);
 
−
 
−
        TH1F *TDC12=new TH1F("TDC12"," stop timing spectrum",1000,3,930);
 
−
TH1F *TDC11=new TH1F("TDC11"," stop timing spectrum",1000,3,930);
 
−
        TH1F *TDC15=new TH1F("TDC15","start by NaI(Tl) det, Cf-252 source",1000,3,930);
 
−
      // TH1F *Fsubtr=new TH1F("Fsubtr","Subtraction of the stops",1000,0,110);
 
−
 
−
  TH1F *Fsubtr=new TH1F("Fsubtr","Subtraction of the stops",1000,-930,930);
 
−
 
−
               
 
−
  TH2D *TDC11evt = new TH2D("TDC11evt", "TDC11evt", 900, 0., 90000., 4090, 0., 4096.);
 
−
  TH2D *TDC12evt = new TH2D("TDC12evt", "TDC12evt", 900, 0., 90000., 4096, 0., 4096.);
 
−
  TH2D *TDCsub = new TH2D("TDCsub", "Eveny-by-event subtraction", 1000, 0., 90000., 1000, -930., 930.);
 
−
 
−
TH2D *Subtr = new TH2D("Subtr", "Subtr", 900, 0., 900., 4096, 0., 4096.);
 
−
//TH1D *Subtr=new TH1D("Subtr","Subtraction of the stops",4096,0,4096);
 
−
TH1D *TDC11proj=new TH1D("TDC11proj","TDC11 projection",1000,0,4096);
 
−
TH1D *TDC12proj=new TH1D("TDC12proj","TDC12 projection",1000,0,4096);
 
−
 
−
TCanvas *c = new TCanvas("c", "TDC", 900, 800);
 
−
 
−
c->Divide(2,2);
 
−
c->Update();
 
−
c->cd(1);
 
−
//ntuple->Draw("500*(TDC12)/4096 >> TDC12");
 
−
//ntuple->Draw("TDC11*0.2272>>TDC11","(EventNumber >=750000 && EventNumber <=1180000)");
 
−
ntuple->Draw("TDC11*0.2272>>TDC11");
 
−
//ADC4->GetYaxis()->SetNdivisions(12);
 
−
        //ADC4->GetXaxis()->SetNdivisions(26);
 
−
c->cd(1)->SetLogy(0);
 
−
c->cd(1)->SetGrid(1,1);
 
−
TDC11->SetXTitle("Time [ns]");
 
−
TDC11->SetYTitle("Events/(1.075)");
 
−
TDC11->SetLineColor(4);
 
−
c->Update();
 
−
 
−
 
−
TLegend *leg = new TLegend(0.5,0.61,0.76,0.92);  //coordinates are fractions
 
−
                                    //of pad dimensions
 
−
 
−
  //the header (or "title") for the legend
 
−
  leg->SetHeader("Run# 3073");
 
−
  leg->AddEntry(TDC11,"Real stop","l");  // "l" means line
 
−
 
−
  leg->AddEntry(TDC11,"rate(Max) =  Hz", "0");            // use "0" for a blank text
 
−
  leg->AddEntry(TDC11,"events =  ", "0");
 
−
  leg->AddEntry(TDC11,"time =  min ", "0");
 
−
  leg->AddEntry(TDC11,"U(right) = -1250 V", "0");
 
−
  leg->AddEntry(TDC11,"U(left) = -1250 V", "0");
 
−
 
 
−
  // oops we forgot the header (or "title") for the legend
 
−
 
 
−
  leg->Draw();
 
−
 
−
c->cd(2);
 
−
//ntuple->Draw("500*(TDC12)/4096 >> TDC12");
 
−
ntuple->Draw("TDC15*0.2272>>TDC15");
 
−
//ADC4->GetYaxis()->SetNdivisions(12);
 
−
        //ADC4->GetXaxis()->SetNdivisions(26);
 
−
c->cd(2)->SetLogy(0);
 
−
c->cd(2)->SetGrid(1,1);
 
−
TDC15->SetXTitle("Time [ns]");
 
−
TDC15->SetYTitle("Events/(1.075)");
 
−
TDC15->SetLineColor(4);
 
−
c->Update();
 
−
 
−
 
−
TLegend *leg = new TLegend(0.5,0.61,0.76,0.92);  //coordinates are fractions
 
−
                                    //of pad dimensions
 
−
 
−
  //the header (or "title") for the legend
 
−
  leg->SetHeader("Run# 3073");
 
−
  leg->AddEntry(TDC15,"Real stop","l");  // "l" means line
 
−
 
−
  leg->AddEntry(TDC15,"rate(Max) =  Hz", "0");            // use "0" for a blank text
 
−
  leg->AddEntry(TDC15,"events =  ", "0");
 
−
  leg->AddEntry(TDC15,"time =  ", "0");
 
−
  leg->AddEntry(TDC15,"U(left) = -1250 V", "0");
 
−
  leg->AddEntry(TDC15,"U(right) = -1250 V", "0");
 
−
 
 
−
  // oops we forgot the header (or "title") for the legend
 
−
 
 
−
  leg->Draw();
 
−
 
−
c->cd(3);
 
−
//ntuple->Draw("500*(TDC12)/4096 >> TDC12");
 
−
//ntuple->Draw("TDC12*0.2272>>TDC12","(EventNumber >=750000 && EventNumber <=1181198)");
 
−
ntuple->Draw("TDC12*0.2272>>TDC12");
 
−
//ADC4->GetYaxis()->SetNdivisions(12);
 
−
        //ADC4->GetXaxis()->SetNdivisions(26);
 
−
c->cd(3)->SetLogy(0);
 
−
c->cd(3)->SetGrid(1,1);
 
−
TDC12->SetXTitle("Time [ns]");
 
−
TDC12->SetYTitle("Events/(1.075)");
 
−
TDC12->SetLineColor(4);
 
−
c->Update();
 
−
 
−
 
−
TLegend *leg = new TLegend(0.5,0.61,0.76,0.92);  //coordinates are fractions
 
−
                                    //of pad dimensions
 
−
 
−
  //the header (or "title") for the legend
 
−
  leg->SetHeader("Run# 3073");
 
−
  leg->AddEntry(TDC15,"Real stop","l");  // "l" means line
 
−
 
−
  leg->AddEntry(TDC12,"rate(Max) =  Hz", "0");            // use "0" for a blank text
 
−
  leg->AddEntry(TDC12,"events =  ", "0");
 
−
  leg->AddEntry(TDC12,"time =  min ", "0");
 
−
  leg->AddEntry(TDC12,"U(left) = -1250 V", "0");
 
−
  leg->AddEntry(TDC12,"U(right) = -1250 V", "0");
 
−
 
 
−
  // oops we forgot the header (or "title") for the legend
 
−
 
 
−
  leg->Draw();
 
−
 
−
c->cd(4);
 
−
c->cd(4)->SetLogy(0);
 
−
c->cd(4)->SetGrid(1,1);
 
−
//ntuple -> Draw("(((TDC11*0.2272)-(TDC12*0.2272))-103.1)/0.1097>>Fsubtr");
 
−
//ntuple -> Draw("(TDC11*0.2272)-(TDC12*0.2272)>>Fsubtr","(EventNumber >=750000 && EventNumber <=1181198)");
 
−
ntuple -> Draw("(TDC11*0.2272)-(TDC12*0.2272)>>Fsubtr");
 
−
Fsubtr->SetXTitle("X-coordinate, cm");
 
−
//TDC15->SetYTitle("Events/(1.075)");
 
−
TDC15->SetLineColor(4);
 
−
 
−
c->Update();
 
−
 
−
TCanvas *c2 = new TCanvas("c2", "TDC", 900, 800);
 
−
 
−
c2->Divide(2,2);
 
−
c2->Update();
 
−
 
−
c2->cd(1);
 
−
ntuple->Draw("(TDC12):EventNumber>>TDC12evt");
 
−
c2->Update();
 
−
 
−
c2->cd(2);
 
−
ntuple->Draw("(TDC11):EventNumber>>TDC11evt");
 
−
c2->Update();
 
−
 
−
c2->cd(3);
 
−
//TH1D *py = TDC11evt->ProjectionY("Y-projection TDC11", 0, 4096);
 
−
 
−
//TH1D *py = TDCsub->ProjectionY("Y-projection", -930, 930); // where firstXbin = 0 and lastXbin = 9
 
−
//py->Draw();
 
−
/*
 
−
Int_t nx11 = TDC11evt->GetNbinsX();
 
−
Int_t ny11 = TDC11evt->GetNbinsY();
 
−
 
−
for (int i=1;i<=nx11;i++) {
 
−
  for (int j=1;j<=ny11;j++) {
 
−
 
 
−
      Double_t c11 =TDC11evt->GetBinContent(i,j);
 
−
      Double_t c12 =TDC12evt->GetBinContent(i,j);
 
−
 
−
  if(c11) Int_t binN11 = j;
 
−
  if(c12) Int_t binN12 = j;
 
−
//if(c12) cout<<"c12: "<<c12<<endl;
 
−
    //if (c11 && c12 >  0.) {
 
−
Double_t a = binN11-binN12;
 
−
Subtr->SetBinContent(i,j,a);
 
−
binN11 = 0;
 
−
binN12 = 0;
 
−
//double a = c11;
 
−
       
 
−
//if(!a) Subtr->SetBinContent(0,0);
 
−
 
−
//a=0.;
 
−
      //    }
 
−
}
 
−
}
 
−
 
−
//TDC11evt->Draw();
 
−
//TDCsub->Draw("TDC11evt*0.2272-TDC12evt*0.2272");
 
−
Subtr->Draw();
 
−
*/
 
−
//TDC11evt->Add(TDC12evt,-1.);
 
−
 
−
//TDC11evt->Draw();
 
−
c2->Update();
 
−
 
−
 
−
c2->cd(4);
 
−
 
−
//TH1D *y2 = TDC12evt->ProjectionY("Y-projection TDC12", 0, 4096);
 
−
 
−
//TH1D *py = TDCsub->ProjectionY("Y-projection", -930, 930); // where firstXbin = 0 and lastXbin = 9
 
−
//y2->Draw();
 
−
 
−
 
−
c2->Update();
 
−
 
−
//TCanvas *c3 = new TCanvas("c3", "TDC", 900, 800);
 
−
 
−
//py->Add(y2,1);
 
−
//py->Draw();
 
−
 
−
char printName[256];
 
−
char printName2[256];
 
−
sprintf(printName,"/data/Oleksiy/Cf_252_test/5_2_2012/run3080.eps");
 
−
c->SaveAs(printName);
 
−
        sprintf(printName2,"/data/Oleksiy/Cf_252_test/5_2_2012/run3080.jpg");
 
−
c->SaveAs(printName2);
 
−
 
−
}
 

Latest revision as of 21:14, 3 May 2012

go back


some Alex program from Sep 2011 run

//to see spectrum for all 9 detectors

  • plot22.C
    • root *.root
    • .x plot22.C

//to asymmetry study

  • NDasymStudy3.C
  • ND_asymmetry.C
  • ND_integral.C
  • NDasymStudy2.C //D2O?

// to monitor beam current stability defined by pair spectrometer // from daq/data machine, Sep 21, 2011

  • NDpairStudy.C
  • ND_pair.C
    • root -l NDpairStudy.C'()'
    • .x NDpairStudy.C()
Retrieved from "https://wiki.iac.isu.edu/index.php?title=Some_programs&oldid=73946"