MUST2 SAMURAI_2018 S3-LEB-LPC SUPERNEMO MUGAST EXPAND SCALP GALATRON HiCARI VELO
  MUGAST  Not logged in MUGAST
Message ID: 27     Entry time: Fri Mar 1 12:07:44 2019
Author: Freddy 
ExpNbr: TEST 
Type: GRIT 
Category: General 
Status: Fixed 
Target-Source: 3-alpha 
Subject: Efficiency reconstruction 
blue: reconstructed efficiency using alpha run 40 (arbitratily normalized)
red: simulated efficiency

There is a difference in angle for trapezoids -> to be investigated

Remark: Simulation curve seems wrong (maybe a target in detector file not removed)
Attachment 1: comparison.pdf  21 kB  | Hide | Hide all
comparison.pdf
Attachment 2: comparison.C  20 kB  | Hide | Hide all
void comparison()
{
//=========Macro generated from canvas: c1/c1
//=========  (Fri Mar  1 12:50:07 2019) by ROOT version 6.16/00
   TCanvas *c1 = new TCanvas("c1", "c1",13,60,856,773);
   c1->ToggleEventStatus();
   c1->Range(1.5e+02,-0.11,1.7e+02,0.98);
   c1->SetFillColor(0);
   c1->SetBorderMode(0);
   c1->SetBorderSize(2);
   c1->SetFrameBorderMode(0);
   c1->SetFrameBorderMode(0);
   
   TH1F *h__1 = new TH1F("h__1","acos(PosZ/sqrt(PosX*PosX+PosY*PosY+PosZ*PosZ))*180./3.141592 {Mugast.TelescopeNumber && Mugast.DSSD_Y!=46}",180,0,1.8e+02);
   h__1->SetBinContent(103,0.0037);
   h__1->SetBinContent(104,0.047);
   h__1->SetBinContent(105,0.36);
   h__1->SetBinContent(106,0.46);
   h__1->SetBinContent(107,0.5);
   h__1->SetBinContent(108,0.49);
   h__1->SetBinContent(109,0.47);
   h__1->SetBinContent(110,0.5);
   h__1->SetBinContent(111,0.46);
   h__1->SetBinContent(112,0.45);
   h__1->SetBinContent(113,0.53);
   h__1->SetBinContent(114,0.54);
   h__1->SetBinContent(115,0.52);
   h__1->SetBinContent(116,0.47);
   h__1->SetBinContent(117,0.47);
   h__1->SetBinContent(118,0.46);
   h__1->SetBinContent(119,0.35);
   h__1->SetBinContent(120,0.44);
   h__1->SetBinContent(121,0.51);
   h__1->SetBinContent(122,0.48);
   h__1->SetBinContent(123,0.46);
   h__1->SetBinContent(124,0.51);
   h__1->SetBinContent(125,0.4);
   h__1->SetBinContent(126,0.44);
   h__1->SetBinContent(127,0.38);
   h__1->SetBinContent(128,0.41);
   h__1->SetBinContent(129,0.41);
   h__1->SetBinContent(130,0.41);
   h__1->SetBinContent(131,0.43);
   h__1->SetBinContent(132,0.49);
   h__1->SetBinContent(133,0.44);
   h__1->SetBinContent(134,0.5);
   h__1->SetBinContent(135,0.41);
   h__1->SetBinContent(136,0.48);
   h__1->SetBinContent(137,0.48);
   h__1->SetBinContent(138,0.42);
   h__1->SetBinContent(139,0.45);
   h__1->SetBinContent(140,0.41);
   h__1->SetBinContent(141,0.42);
   h__1->SetBinContent(142,0.41);
   h__1->SetBinContent(143,0.4);
   h__1->SetBinContent(144,0.42);
   h__1->SetBinContent(145,0.4);
   h__1->SetBinContent(146,0.41);
   h__1->SetBinContent(147,0.41);
   h__1->SetBinContent(148,0.3);
   h__1->SetBinContent(149,0.32);
   h__1->SetBinContent(150,0.28);
   h__1->SetBinContent(151,0.28);
   h__1->SetBinContent(152,0.2);
   h__1->SetBinContent(153,0.048);
   h__1->SetBinContent(159,0.0039);
   h__1->SetBinContent(160,0.6);
   h__1->SetBinContent(161,0.73);
   h__1->SetBinContent(162,0.7);
   h__1->SetBinContent(163,0.82);
   h__1->SetBinContent(164,0.79);
   h__1->SetBinContent(165,0.9);
   h__1->SetBinContent(166,0.92);
   h__1->SetBinContent(167,1.1);
   h__1->SetBinContent(168,0.97);
   h__1->SetBinContent(169,0.96);
   h__1->SetBinContent(170,0.13);
   h__1->SetBinError(103,0.00054);
   h__1->SetBinError(104,0.0019);
   h__1->SetBinError(105,0.0053);
   h__1->SetBinError(106,0.006);
   h__1->SetBinError(107,0.0062);
   h__1->SetBinError(108,0.0062);
   h__1->SetBinError(109,0.0061);
   h__1->SetBinError(110,0.0063);
   h__1->SetBinError(111,0.006);
   h__1->SetBinError(112,0.0059);
   h__1->SetBinError(113,0.0064);
   h__1->SetBinError(114,0.0065);
   h__1->SetBinError(115,0.0064);
   h__1->SetBinError(116,0.0061);
   h__1->SetBinError(117,0.0061);
   h__1->SetBinError(118,0.006);
   h__1->SetBinError(119,0.0052);
   h__1->SetBinError(120,0.0059);
   h__1->SetBinError(121,0.0063);
   h__1->SetBinError(122,0.0061);
   h__1->SetBinError(123,0.006);
   h__1->SetBinError(124,0.0063);
   h__1->SetBinError(125,0.0056);
   h__1->SetBinError(126,0.0059);
   h__1->SetBinError(127,0.0055);
   h__1->SetBinError(128,0.0057);
   h__1->SetBinError(129,0.0057);
   h__1->SetBinError(130,0.0056);
   h__1->SetBinError(131,0.0058);
   h__1->SetBinError(132,0.0062);
   h__1->SetBinError(133,0.0059);
   h__1->SetBinError(134,0.0062);
   h__1->SetBinError(135,0.0057);
   h__1->SetBinError(136,0.0061);
   h__1->SetBinError(137,0.0061);
   h__1->SetBinError(138,0.0058);
   h__1->SetBinError(139,0.0059);
   h__1->SetBinError(140,0.0057);
   h__1->SetBinError(141,0.0058);
   h__1->SetBinError(142,0.0057);
   h__1->SetBinError(143,0.0056);
   h__1->SetBinError(144,0.0057);
   h__1->SetBinError(145,0.0056);
   h__1->SetBinError(146,0.0057);
   h__1->SetBinError(147,0.0057);
   h__1->SetBinError(148,0.0048);
   h__1->SetBinError(149,0.005);
   h__1->SetBinError(150,0.0047);
   h__1->SetBinError(151,0.0047);
   h__1->SetBinError(152,0.004);
   h__1->SetBinError(153,0.0019);
   h__1->SetBinError(159,0.00055);
   h__1->SetBinError(160,0.0069);
   h__1->SetBinError(161,0.0076);
   h__1->SetBinError(162,0.0074);
   h__1->SetBinError(163,0.008);
   h__1->SetBinError(164,0.0079);
   h__1->SetBinError(165,0.0084);
   h__1->SetBinError(166,0.0085);
   h__1->SetBinError(167,0.0092);
   h__1->SetBinError(168,0.0087);
   h__1->SetBinError(169,0.0087);
   h__1->SetBinError(170,0.0031);
   h__1->SetEntries(3.7e+05);
   
   TPaveStats *ptstats = new TPaveStats(0.9,0.83,0.98,0.94,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(1);
   ptstats->SetFillColor(0);
   ptstats->SetTextAlign(12);
   ptstats->SetTextFont(42);
   TText *ptstats_LaTex = ptstats->AddText("h");
   ptstats_LaTex->SetTextSize(0.023);
   ptstats_LaTex = ptstats->AddText("Entries = 373600 ");
   ptstats_LaTex = ptstats->AddText("Mean  =  157.2");
   ptstats_LaTex = ptstats->AddText("Std Dev   =  5.644");
   ptstats->SetOptStat(1111);
   ptstats->SetOptFit(0);
   ptstats->Draw();
   h__1->GetListOfFunctions()->Add(ptstats);
   ptstats->SetParent(h__1);

   Int_t ci;      // for color index setting
   TColor *color; // for color definition with alpha
   ci = TColor::GetColor("#000099");
   h__1->SetLineColor(ci);
   h__1->GetXaxis()->SetRange(148,163);
   h__1->GetXaxis()->SetLabelFont(42);
   h__1->GetXaxis()->SetLabelSize(0.035);
   h__1->GetXaxis()->SetTitleSize(0.035);
   h__1->GetXaxis()->SetTitleOffset(1);
   h__1->GetXaxis()->SetTitleFont(42);
   h__1->GetYaxis()->SetLabelFont(42);
   h__1->GetYaxis()->SetLabelSize(0.035);
   h__1->GetYaxis()->SetTitleSize(0.035);
   h__1->GetYaxis()->SetTitleFont(42);
   h__1->GetZaxis()->SetLabelFont(42);
   h__1->GetZaxis()->SetLabelSize(0.035);
   h__1->GetZaxis()->SetTitleSize(0.035);
   h__1->GetZaxis()->SetTitleOffset(1);
   h__1->GetZaxis()->SetTitleFont(42);
   h__1->Draw("EHIST");
   
   TPaveText *pt = new TPaveText(0.15,0.93,0.85,1,"blNDC");
   pt->SetName("title");
   pt->SetBorderSize(0);
   pt->SetFillColor(0);
   pt->SetFillStyle(0);
   pt->SetTextFont(42);
   TText *pt_LaTex = pt->AddText("acos(PosZ/sqrt(PosX*PosX+PosY*PosY+PosZ*PosZ))*180./3.141592 {Mugast.TelescopeNumber && Mugast.DSSD_Y!=46}");
   pt->Draw();
   
   TH1F *hEfficiency__2 = new TH1F("hEfficiency__2","Efficiency",180,0,1.8e+02);
   hEfficiency__2->SetBinContent(5,0.0014);
   hEfficiency__2->SetBinContent(7,0.012);
   hEfficiency__2->SetBinContent(8,0.085);
   hEfficiency__2->SetBinContent(9,0.2);
   hEfficiency__2->SetBinContent(10,0.33);
   hEfficiency__2->SetBinContent(11,0.42);
   hEfficiency__2->SetBinContent(12,0.47);
   hEfficiency__2->SetBinContent(13,0.53);
   hEfficiency__2->SetBinContent(14,0.57);
   hEfficiency__2->SetBinContent(15,0.61);
   hEfficiency__2->SetBinContent(16,0.63);
   hEfficiency__2->SetBinContent(17,0.68);
   hEfficiency__2->SetBinContent(18,0.69);
   hEfficiency__2->SetBinContent(19,0.72);
   hEfficiency__2->SetBinContent(20,0.72);
   hEfficiency__2->SetBinContent(21,0.74);
   hEfficiency__2->SetBinContent(22,0.76);
   hEfficiency__2->SetBinContent(23,0.77);
   hEfficiency__2->SetBinContent(24,0.78);
   hEfficiency__2->SetBinContent(25,0.8);
   hEfficiency__2->SetBinContent(26,0.8);
   hEfficiency__2->SetBinContent(27,0.82);
   hEfficiency__2->SetBinContent(28,0.84);
   hEfficiency__2->SetBinContent(29,0.83);
   hEfficiency__2->SetBinContent(30,0.85);
   hEfficiency__2->SetBinContent(31,0.85);
   hEfficiency__2->SetBinContent(32,0.86);
   hEfficiency__2->SetBinContent(33,0.86);
   hEfficiency__2->SetBinContent(34,0.84);
   hEfficiency__2->SetBinContent(35,0.79);
   hEfficiency__2->SetBinContent(36,0.74);
   hEfficiency__2->SetBinContent(37,0.72);
   hEfficiency__2->SetBinContent(38,0.64);
   hEfficiency__2->SetBinContent(39,0.57);
   hEfficiency__2->SetBinContent(40,0.48);
   hEfficiency__2->SetBinContent(41,0.4);
   hEfficiency__2->SetBinContent(42,0.34);
   hEfficiency__2->SetBinContent(43,0.29);
   hEfficiency__2->SetBinContent(44,0.24);
   hEfficiency__2->SetBinContent(45,0.21);
   hEfficiency__2->SetBinContent(46,0.16);
   hEfficiency__2->SetBinContent(47,0.12);
   hEfficiency__2->SetBinContent(48,0.08);
   hEfficiency__2->SetBinContent(49,0.046);
   hEfficiency__2->SetBinContent(50,0.018);
   hEfficiency__2->SetBinContent(51,0.0071);
   hEfficiency__2->SetBinContent(52,0.00058);
   hEfficiency__2->SetBinContent(53,0.0014);
   hEfficiency__2->SetBinContent(54,0.0076);
   hEfficiency__2->SetBinContent(55,0.034);
   hEfficiency__2->SetBinContent(56,0.096);
   hEfficiency__2->SetBinContent(57,0.17);
   hEfficiency__2->SetBinContent(58,0.21);
   hEfficiency__2->SetBinContent(59,0.21);
   hEfficiency__2->SetBinContent(60,0.22);
   hEfficiency__2->SetBinContent(61,0.21);
   hEfficiency__2->SetBinContent(62,0.22);
   hEfficiency__2->SetBinContent(63,0.21);
   hEfficiency__2->SetBinContent(64,0.21);
   hEfficiency__2->SetBinContent(65,0.22);
   hEfficiency__2->SetBinContent(66,0.22);
   hEfficiency__2->SetBinContent(67,0.16);
   hEfficiency__2->SetBinContent(68,0.05);
   hEfficiency__2->SetBinContent(69,0.0067);
   hEfficiency__2->SetBinContent(70,0.00072);
   hEfficiency__2->SetBinContent(71,0.00037);
   hEfficiency__2->SetBinContent(72,0.00012);
   hEfficiency__2->SetBinContent(73,0.00024);
   hEfficiency__2->SetBinContent(84,0.00023);
   hEfficiency__2->SetBinContent(89,0.00011);
   hEfficiency__2->SetBinContent(90,0.00012);
   hEfficiency__2->SetBinContent(92,0.00011);
   hEfficiency__2->SetBinContent(93,0.00011);
   hEfficiency__2->SetBinContent(96,0.00011);
   hEfficiency__2->SetBinContent(101,0.00012);
   hEfficiency__2->SetBinContent(103,0.00012);
   hEfficiency__2->SetBinContent(104,0.00012);
   hEfficiency__2->SetBinContent(105,0.00012);
   hEfficiency__2->SetBinContent(106,0.00036);
   hEfficiency__2->SetBinContent(107,0.00024);
   hEfficiency__2->SetBinContent(108,0.00071);
   hEfficiency__2->SetBinContent(109,0.0013);
   hEfficiency__2->SetBinContent(110,0.006);
   hEfficiency__2->SetBinContent(111,0.049);
   hEfficiency__2->SetBinContent(112,0.25);
   hEfficiency__2->SetBinContent(113,0.47);
   hEfficiency__2->SetBinContent(114,0.53);
   hEfficiency__2->SetBinContent(115,0.54);
   hEfficiency__2->SetBinContent(116,0.53);
   hEfficiency__2->SetBinContent(117,0.54);
   hEfficiency__2->SetBinContent(118,0.54);
   hEfficiency__2->SetBinContent(119,0.55);
   hEfficiency__2->SetBinContent(120,0.53);
   hEfficiency__2->SetBinContent(121,0.55);
   hEfficiency__2->SetBinContent(122,0.54);
   hEfficiency__2->SetBinContent(123,0.53);
   hEfficiency__2->SetBinContent(124,0.53);
   hEfficiency__2->SetBinContent(125,0.52);
   hEfficiency__2->SetBinContent(126,0.53);
   hEfficiency__2->SetBinContent(127,0.53);
   hEfficiency__2->SetBinContent(128,0.53);
   hEfficiency__2->SetBinContent(129,0.52);
   hEfficiency__2->SetBinContent(130,0.53);
   hEfficiency__2->SetBinContent(131,0.53);
   hEfficiency__2->SetBinContent(132,0.51);
   hEfficiency__2->SetBinContent(133,0.52);
   hEfficiency__2->SetBinContent(134,0.51);
   hEfficiency__2->SetBinContent(135,0.51);
   hEfficiency__2->SetBinContent(136,0.5);
   hEfficiency__2->SetBinContent(137,0.52);
... 225 more lines ...
Attachment 3: alpha_eff.root  5 kB

ELOG Home Page