modified: MakeVertex.C

MakeVertex.C

@@ -93,6 +93,7 @@ double qqqGain[MAX_QQQ][MAX_RING][MAX_WEDGE] = {{{0}}};
 bool qqqGainValid[MAX_QQQ][MAX_RING][MAX_WEDGE] = {{{false}}};
 double qqqCalib[MAX_QQQ][MAX_RING][MAX_WEDGE] = {{{0}}};
 bool qqqCalibValid[MAX_QQQ][MAX_RING][MAX_WEDGE] = {{{false}}};
+bool Diagnostics = false;
 
 double sx3BackGain[24][4][4] = {{{1.}}};
 double sx3FrontGain[24][4] = {{1.}};

run_sx3.sh

@@ -38,6 +38,7 @@ fi
 if [[ 1 -eq 1 ]] ; then
 export anode_offset=1
 root -q -b -x ../Remapped_files/27Al_Data/root_data/Run_015_mapped.root -e 'tree->Process("MakeVertex.C+O")'; mv Analyzer_SX3.root results_run15.root;
+exit
 root -q -b -x ../Remapped_files/27Al_Data/root_data/Run_017_mapped.root -e 'tree->Process("MakeVertex.C+O")'; mv Analyzer_SX3.root results_run17.root;
 root -q -b -x ../Remapped_files/27Al_Data/root_data/Run_018_mapped.root -e 'tree->Process("MakeVertex.C+O")'; mv Analyzer_SX3.root results_run18.root;
 root -q -b -x ../Remapped_files/27Al_Data/root_data/Run_019_mapped.root -e 'tree->Process("MakeVertex.C+O")'; mv Analyzer_SX3.root results_run19.root;

sx3cal/17F/backgains.dat

@@ -1,23 +1,55 @@
-1 front 1 back 2 4.4094
+1 front 0 back 3 3.68208
-1 front 1 back 1 4.4094
+1 front 0 back 0 3.89631
-1 front 2 back 2 4.4832
+1 front 1 back 2 4.44599
-1 front 3 back 2 4.52103
+1 front 1 back 3 3.81311
-1 front 3 back 1 4.5210
+1 front 1 back 0 3.98321
-3 front 1 back 2 3.63215
+1 front 2 back 2 4.52103
-3 front 1 back 1 3.70756
+1 front 2 back 3 3.68208
-3 front 2 back 2 3.68208
+1 front 2 back 0 3.95382
-3 front 2 back 1 3.86817
+1 front 3 back 2 4.55951
-3 front 3 back 2 3.7334
+1 front 3 back 3 3.68208
-3 front 3 back 1 3.70756
+1 front 3 back 0 3.89631
-7 front 1 back 2 3.60769
+3 front 1 back 2 3.81311
-7 front 1 back 1 3.46759
+3 front 1 back 1 3.92485
-7 front 2 back 2 3.58356
+3 front 1 back 3 4.13701
-7 front 2 back 1 3.49018
+3 front 1 back 0 4.26886
-7 front 3 back 2 3.60769
+3 front 2 back 2 3.86817
-7 front 3 back 1 3.46759
+3 front 2 back 1 4.04334
-9 front 1 back 2 3.58356
+3 front 2 back 3 3.98321
-9 front 1 back 1 3.44529
+3 front 2 back 0 4.23512
-9 front 2 back 2 3.58356
+3 front 3 back 2 3.92485
-9 front 2 back 1 3.46759
+3 front 3 back 1 3.92485
-9 front 3 back 2 3.63215
+3 front 3 back 3 3.95382
-9 front 3 back 1 3.46759
+3 front 3 back 0 4.13701
+7 front 0 back 2 3.65694
+7 front 0 back 1 3.53625
+7 front 0 back 3 4.23512
+7 front 0 back 0 3.60769
+7 front 1 back 2 3.65694
+7 front 1 back 1 3.53625
+7 front 1 back 3 4.23512
+7 front 1 back 0 3.58356
+7 front 2 back 2 3.68208
+7 front 2 back 1 3.55975
+7 front 2 back 3 4.13701
+7 front 2 back 0 3.60769
+7 front 3 back 2 3.68208
+7 front 3 back 1 3.53625
+7 front 3 back 3 4.23512
+7 front 3 back 0 3.58356
+9 front 0 back 2 3.68208
+9 front 0 back 1 3.53625
+9 front 0 back 3 3.51306
+9 front 0 back 0 3.42327
+9 front 1 back 2 3.70756
+9 front 1 back 1 3.49018
+9 front 1 back 3 3.63215
+9 front 1 back 0 3.46759
+9 front 2 back 2 3.70756
+9 front 2 back 1 3.55975
+9 front 2 back 3 3.53625
+9 front 2 back 0 3.44529
+9 front 3 back 2 3.7334
+9 front 3 back 1 3.51306
+9 front 3 back 3 3.53625
+9 front 3 back 0 3.42327

sx3cal/17F/frontgains.dat

@@ -1,12 +1,15 @@
-1 lengthcal front 1 1.5121 60.4839
+1 lengthcal front 0 0.302419 60.4839
-1 lengthcal front 2 -1.5625 62.5
+1 lengthcal front 1 0.907258 60.4839
-1 lengthcal front 3 2.72177 60.4839
+1 lengthcal front 2 -2.11694 60.4839
-3 lengthcal front 1 -0.595088 59.5088
+1 lengthcal front 3 3.32661 60.4839
-3 lengthcal front 2 -4.53935 58.5723
+3 lengthcal front 1 -0.302419 60.4839
-3 lengthcal front 3 4.08107 60.4603
+3 lengthcal front 2 -3.12146 59.4564
-7 lengthcal front 1 1.14329 45.7317
+3 lengthcal front 3 5.3066 62.4306
-7 lengthcal front 2 0.115661 46.2646
+7 lengthcal front 0 0.225805 45.1611
-7 lengthcal front 3 2.90179 44.6429
+7 lengthcal front 1 1.85157 46.2892
-9 lengthcal front 1 0.115732 46.2928
+7 lengthcal front 2 0.936622 46.8311
-9 lengthcal front 2 0.799176 45.6672
+7 lengthcal front 3 3.65527 45.6909
-9 lengthcal front 3 1.68159 48.0453
+9 lengthcal front 0 2.10806 46.8457
+9 lengthcal front 1 0.694127 46.2751
+9 lengthcal front 2 1.04095 46.2646
+9 lengthcal front 3 1.99048 46.8348

sx3cal/17F/rightgains.dat

@@ -1,16 +1,16 @@
-1 0 1678.38 1.0
+1 0 1239.16 0.649287
-1 1 1678.38 1.07163
+1 1 1853.11 1.08968
-1 2 1693.99 1.1035
+1 2 1852.65 1.10279
-1 3 1667.9 0.975015
+1 3 1821.01 0.965101
-3 0 1597.96 1.0
+3 0 1749.6 1.
-3 1 1597.96 1.02536
+3 1 1749.6 1.01204
-3 2 1821.48 1.29182
+3 2 1954.81 1.23687
-3 3 1607.52 0.928543
+3 3 1751.15 0.904329
-7 0 1773.34 1.
+9 0 1838.02 1.07361
-7 1 1773.34 1.14263
+9 1 1695.51 0.971866
-7 2 1573.79 1.06715
+9 2 1709.36 0.988575
-7 3 1542.41 0.956475
+9 3 1719.37 0.982937
-9 0 1555.08 1.
+7 0 1787.9 1.02266
-9 1 1555.08 0.998252
+7 1 1924.89 1.10649
-9 2 1559.36 1.00299
+7 2 1706.01 1.03175
-9 3 1585.79 1.01582
+7 3 1685.88 0.933633

sx3cal/EXFit.C

@@ -1,6 +1,6 @@
 {
     int index = 1;
-    TFile *f = new TFile("../results_run19.root");
+    TFile *f = new TFile("../results_run01.root");
     TH2F *h2=NULL;
     TH1F *h1x=NULL, *h1y=NULL;
     //f->cd("evsx");
@@ -11,7 +11,7 @@
 
     std::ofstream ofb(Form("backgains%d.dat",index));
     std::ofstream off(Form("frontgains%d.dat",index));
-    for(int i=1; i<4; i++) {
+    for(int i=0; i<4; i++) {
         //do it for pad#2
         int backnum=2;
         h2 = (TH2F*)(f->Get(Form("evsx/be_vs_x_sx3_id_%d_f%d_b%d",index,i,backnum)));

sx3cal/EXFit2.C

@@ -0,0 +1,79 @@
+{
+    int index = 9;
+    TFile *f = new TFile("../results_run05.root");
+    TH2F *h2=NULL;
+    TH1F *h1x=NULL, *h1y=NULL;
+    //f->cd("evsx");
+    //f->ls();
+
+    double known_xpos[4][2] = {{0,18.75},{-18.75,0}};
+    std::vector<double> xpos, xposkn; //first = x = known position, second = y = unknown position
+
+    std::ofstream ofb(Form("backgains%d.dat",index));
+    std::ofstream off(Form("frontgains%d.dat",index));
+    for(int i=0; i<4; i++) {
+        //do it for pad#2
+        int backnum=2;
+        h2 = (TH2F*)(f->Get(Form("evsx/be_vs_x_sx3_id_%d_f%d_b%d",index,i,backnum)));
+        auto macro = [&]() {
+            h1x = (TH1F*)(h2->ProjectionX("_px"));
+            double xleft = h1x->GetBinCenter(h1x->FindFirstBinAbove(h1x->GetMaximum()*0.4));
+            double xright = h1x->GetBinCenter(h1x->FindLastBinAbove(h1x->GetMaximum()*0.4));
+            //h1x->GetXaxis()->SetRangeUser(4*xleft, xright*4);
+            h1x->Draw();
+            TLine L1(xleft,0,xleft,h1x->GetMaximum()); L1.SetLineColor(kRed); L1.Draw("SAME");
+            TLine L2(xright,0,xright,h1x->GetMaximum()); L2.SetLineColor(kRed); L2.Draw("SAME");
+            gPad->Modified();
+            gPad->Update();
+            if(backnum==1 || backnum==2) xpos.push_back(xleft); xposkn.push_back(known_xpos[backnum-1][0]);
+            if(backnum==1 || backnum==2) xpos.push_back(xright); xposkn.push_back(known_xpos[backnum-1][1]);
+            while(gPad->WaitPrimitive());
+
+            h1y = (TH1F*)(h2->ProjectionY("_py"));
+            double ycenter = h1y->GetBinCenter(h1y->GetMaximumBin());
+//            std::cout << "front " << i << " back " << backnum << " "  << xleft << " " << xright << " " << ycenter << " " << 5486/ycenter << std::endl;
+            ofb << index <<" front " << i << " back " << backnum << " "  << 5486/ycenter << std::endl;
+            h1y->GetXaxis()->SetRangeUser(ycenter-200,ycenter+200);
+            h1y->Draw();
+            TLine L3(ycenter,0,ycenter,h1y->GetMaximum()*1.1); L3.SetLineColor(kRed); L3.Draw("SAME");
+
+            gPad->Modified();
+            gPad->Update();
+            while(gPad->WaitPrimitive());
+        };
+        if(h2)
+            macro();
+
+        //repeat for pad#1
+        backnum=1;
+        h2 = (TH2F*)(f->Get(Form("evsx/be_vs_x_sx3_id_%d_f%d_b%d",index,i,backnum)));
+        if(h2)
+            macro();
+
+        //repeat for pad#2
+        backnum=3;
+        h2 = (TH2F*)(f->Get(Form("evsx/be_vs_x_sx3_id_%d_f%d_b%d",index,i,backnum)));
+        if(h2)
+            macro();
+
+        //repeat for pad3
+        backnum=0;
+        h2 = (TH2F*)(f->Get(Form("evsx/be_vs_x_sx3_id_%d_f%d_b%d",index,i,backnum)));
+        if(h2)
+            macro();
+
+        double xtofit[] = {xpos[0],xpos[3]};
+        double xktofit[] = {xposkn[0],xposkn[3]};
+        TGraph G1(xpos.size(),xpos.data(),xposkn.data());
+        G1.Draw("APL*");
+        G1.Fit("pol1","Q");
+        off << index<<" lengthcal front " << i << " " << G1.GetFunction("pol1")->GetParameter(0) << " " << G1.GetFunction("pol1")->GetParameter(1) << std::endl;
+        gPad->Modified(); gPad->Update();
+        while(gPad->WaitPrimitive());
+        xpos.clear();
+        xposkn.clear();
+    }
+    ofb.close();
+    off.close();
+    f->Close();
+}

sx3cal/LRFit.C

@@ -1,10 +1,10 @@
 {
-    TFile *f = new TFile("../results_run19.root");
+    TFile *f = new TFile("../results_run01.root");
     f->cd("l_vs_r");
     gDirectory->ls();
-    int clkpos = 13;
+    int clkpos = 1;
     std::ofstream ofile(Form("rightgains%d.dat",clkpos));
-    for(int i=1; i<4; i++) {
+    for(int i=0; i<4; i++) {
         TH2F h2(*(TH2F*)(f->Get(Form("l_vs_r/l_vs_r_sx3_id_%d_f%d",clkpos,i))));
         h2.Draw();
 

sx3cal/backgains1.dat

@@ -1,3 +1,12 @@
-1 front 1 back 2 4.4094
+1 front 0 back 2 4.44599
-1 front 2 back 2 4.4832
+1 front 0 back 3 3.68208
-1 front 3 back 2 4.52103
+1 front 0 back 0 3.89631
+1 front 1 back 2 4.44599
+1 front 1 back 3 3.81311
+1 front 1 back 0 3.98321
+1 front 2 back 2 4.52103
+1 front 2 back 3 3.68208
+1 front 2 back 0 3.95382
+1 front 3 back 2 4.55951
+1 front 3 back 3 3.68208
+1 front 3 back 0 3.89631

sx3cal/backgains3.dat

@@ -0,0 +1,12 @@
+3 front 1 back 2 3.81311
+3 front 1 back 1 3.92485
+3 front 1 back 3 4.13701
+3 front 1 back 0 4.26886
+3 front 2 back 2 3.86817
+3 front 2 back 1 4.04334
+3 front 2 back 3 3.98321
+3 front 2 back 0 4.23512
+3 front 3 back 2 3.92485
+3 front 3 back 1 3.92485
+3 front 3 back 3 3.95382
+3 front 3 back 0 4.13701

sx3cal/backgains7.dat

@@ -0,0 +1,16 @@
+7 front 0 back 2 3.65694
+7 front 0 back 1 3.53625
+7 front 0 back 3 4.23512
+7 front 0 back 0 3.60769
+7 front 1 back 2 3.65694
+7 front 1 back 1 3.53625
+7 front 1 back 3 4.23512
+7 front 1 back 0 3.58356
+7 front 2 back 2 3.68208
+7 front 2 back 1 3.55975
+7 front 2 back 3 4.13701
+7 front 2 back 0 3.60769
+7 front 3 back 2 3.68208
+7 front 3 back 1 3.53625
+7 front 3 back 3 4.23512
+7 front 3 back 0 3.58356

sx3cal/backgains9.dat

@@ -1,6 +1,16 @@
-9 front 1 back 2 3.58356
+9 front 0 back 2 3.68208
-9 front 1 back 1 3.44529
+9 front 0 back 1 3.53625
-9 front 2 back 2 3.58356
+9 front 0 back 3 3.51306
-9 front 2 back 1 3.46759
+9 front 0 back 0 3.42327
-9 front 3 back 2 3.63215
+9 front 1 back 2 3.70756
-9 front 3 back 1 3.46759
+9 front 1 back 1 3.49018
+9 front 1 back 3 3.63215
+9 front 1 back 0 3.46759
+9 front 2 back 2 3.70756
+9 front 2 back 1 3.55975
+9 front 2 back 3 3.53625
+9 front 2 back 0 3.44529
+9 front 3 back 2 3.7334
+9 front 3 back 1 3.51306
+9 front 3 back 3 3.53625
+9 front 3 back 0 3.42327

sx3cal/frontgains1.dat

@@ -1,3 +1,4 @@
-1 lengthcal front 1 1.5121 60.4839
+1 lengthcal front 0 -0.292969 58.5938
-1 lengthcal front 2 -1.5625 62.5
+1 lengthcal front 1 0.302419 60.4839
+1 lengthcal front 2 -2.55682 56.8182
 1 lengthcal front 3 2.72177 60.4839

sx3cal/frontgains3.dat

@@ -0,0 +1,3 @@
+3 lengthcal front 1 0.468555 62.474
+3 lengthcal front 2 -3.17087 60.3975
+3 lengthcal front 3 5.24289 63.5502

sx3cal/frontgains7.dat

@@ -0,0 +1,4 @@
+7 lengthcal front 0 0 44.0628
+7 lengthcal front 1 1.68159 48.0453
+7 lengthcal front 2 0.80963 46.2646
+7 lengthcal front 3 3.5855 46.2646

sx3cal/frontgains9.dat

@@ -1,3 +1,4 @@
-9 lengthcal front 1 0.115732 46.2928
+9 lengthcal front 0 2.10937 46.875
-9 lengthcal front 2 0.799176 45.6672
+9 lengthcal front 1 0.57866 46.2928
-9 lengthcal front 3 1.68159 48.0453
+9 lengthcal front 2 1.17187 46.875
+9 lengthcal front 3 1.6345 46.6999

sx3cal/rightgains1.dat

@@ -1,3 +1,4 @@
-1 1 1678.38 1.07163
+1 0 1229.18 1.00572
-1 2 1693.99 1.1035
+1 1 1765.87 0.968534
-1 3 1667.9 0.975015
+1 2 1778.36 0.877985
+1 3 1843.61 0.985871

sx3cal/rightgains3.dat

@@ -0,0 +1,3 @@
+3 1 1749.6 1.01204
+3 2 1954.81 1.23687
+3 3 1751.15 0.904329

sx3cal/rightgains7.dat

@@ -0,0 +1,4 @@
+7 0 1787.9 1.02266
+7 1 1924.89 1.10649
+7 2 1706.01 1.03175
+7 3 1685.88 0.933633

sx3cal/rightgains9.dat

@@ -1,3 +1,4 @@
-9 1 1555.08 0.998252
+9 0 1838.02 1.07361
-9 2 1559.36 1.00299
+9 1 1695.51 0.971866
-9 3 1585.79 1.01582
+9 2 1709.36 0.988575
+9 3 1719.37 0.982937