Difference between revisions of "Some programs"

From New IAC Wiki
Jump to navigation Jump to search
(Created page 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 //to…')
 
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);
 +
 +
}

Revision as of 21:12, 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()


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

}

Retrieved from "https://wiki.iac.isu.edu/index.php?title=Some_programs&oldid=73945"