Changeset 235 in Idarraga

Timestamp:

Sep 20, 2011, 1:42:11 PM (13 years ago)

Author:

idarraga

Message:

 

Location:

mafalda/ShallowAngleFEIX

Files:

• ShallowAngleFEIX.cpp (modified) (3 diffs)
• ShallowAngleFEIX.h (modified) (1 diff)
• configuration_algos_ShallowAngleFEIX.txt (modified) (1 diff)
• runNewAlgo.C (modified) (10 diffs)

mafalda/ShallowAngleFEIX/ShallowAngleFEIX.cpp

@@ -36 +36 @@
         getMyTree()->Branch("clusterSize", &m_clusterSize, "clusterSize/I");
         getMyTree()->Branch("mipLength", &m_mipLength);
+        getMyTree()->Branch("mipLength_X", &m_mipLength_X);
         getMyTree()->Branch("totPerSegment", &m_totPerSegment);
         getMyTree()->Branch("chargeWeights", &m_chargeWeights);
@@ -116 +117 @@
                                         limitPixels[__TAIL_INDX].second * m_pixSizeY);
                         m_mipLength.push_back( totalLength ); // mm
+                        m_mipLength_X.push_back( TMath::Abs ( limitPixels[__HEAD_INDX].first * m_pixSizeX
+                                        - limitPixels[__TAIL_INDX].first * m_pixSizeX ) );
                         m_totalMIPsLength += totalLength;
                         Log << MSG::DEBUG << "Integrated mips length = " << m_totalMIPsLength << " [mm]" << endreq;
@@ -204 +207 @@
                 }
 
-                // other clusters != MIP
+                // there might be more than 1 good mip in a frame
+                // Fill and clean up here
+
+                // Fill the output tree of this algorithm
+                getMyTree()->Fill();
+
+                // WARNING ! don't forget to clean up your variables for the next TTree::Fill call
+                m_distanceToCenter.clear();
+                m_totPerSegment.clear();
+                m_chargeWeights.clear();
+                m_mipLength.clear();
+                m_mipLength_X.clear();
 
         }
 
-        // Fill the output tree of this algorithm
-        getMyTree()->Fill();
-
-        // WARNING ! don't forget to clean up your variables for the next TTree::Fill call
-        m_distanceToCenter.clear();
-        m_totPerSegment.clear();
-        m_chargeWeights.clear();
 
 }

mafalda/ShallowAngleFEIX/ShallowAngleFEIX.h

@@ -47 +47 @@
   Int_t m_clusterSize;
   vector<Float_t> m_mipLength;
+  vector<Float_t> m_mipLength_X;
   Double_t m_totalMIPsLength;
   Int_t m_clusterHeight;

mafalda/ShallowAngleFEIX/configuration_algos_ShallowAngleFEIX.txt

@@ -1 @@
-FEI3Mips minMipLength 2 float
-FEI3Mips pixelSizeX 0.4 float
-FEI3Mips pixelSizeY 0.05 float
-FEI3Mips dRayBalSearch 1.2 float
 PRBasicSpecies circleToEllipseMin 1.4 float
 PRBasicSpecies circleToEllipseMax 5 float
 PRBasicSpecies mindistanceToLine 15 float
 PRBasicSpecies fractionOfPixelsAtMindistance 0.8 float
+ShallowAngleFEIX pixelSizeX 0.25 float
+ShallowAngleFEIX pixelSizeY 0.05 float
+ShallowAngleFEIX dRayBalSearch 1.4 float
 BlobsFinder border 0 int
 BlobsFinder discontinuity 0 int
 BlobsFinder lvl1cut -1 int
-FEI3Mips minNPixels 10 int
-FEI3Mips nDivisions 6 int
-FEI3Mips distanceCutHot 3 int
-FEI3Mips distancePerpHot 3 int
-FEI3Mips guardDistanceX 1 int
-FEI3Mips guardDistanceY 5 int
 PRBasicSpecies nInnerPixels 2 int
 PRBasicSpecies longGammaMax 3 int
 PRBasicSpecies minNPixelsMip 6 int
 PRBasicSpecies maxNPixelsCurly 50 int
+ShallowAngleFEIX minNPixels 5 int
+ShallowAngleFEIX minXWidth 6 int
+ShallowAngleFEIX guardDistanceX 1 int
+ShallowAngleFEIX guardDistanceY 5 int

mafalda/ShallowAngleFEIX/runNewAlgo.C

@@ -4 +4 @@
 
 //#define __NDIV          8
-#define __THICKNESS 0.320 // [mm]
+//#define __THICKNESS 0.320 // [mm]  KEK FE-I3
+#define __THICKNESS 0.150 // [mm]    KEK FE-I4 run 50734 to 50744
+
+#define __CHARGE  10
+#define __ENTRIES 20
 
 void runNewAlgo(){
@@ -10 +14 @@
         gROOT->ProcessLine(".x ~/styles/AtlasStyle.C");
 
-
         TChain * T1 = new TChain("ShallowAngleFEIX");
-        T1->Add("/home/idarraga/analysis/mafalda/MAFOutput_MPXNtuple_USBPIX_10_dep.root");
+        T1->Add("/storage1/idarraga/Testbeam_July_CERN_2011/KEK_100V_4deg/MAFOutput/MAFOutput_MPXNtuple_USBPIXI4_20_100V.root");
+        TString titleT = "KEK FE-I4 ~4deg DUT3";
+        //T1->Add("/home/idarraga/analysis/mafalda/MAFOutput_MPXNtuple_USBPIX_10_dep.root");
         //T->Add("/home/idarraga/analysis/mafalda/MAFOutput_KEK_run050620.root");
         //T->Add("/home/idarraga/analysis/mafalda/MAFOutput_KEK_run050621.root");
         //T->Add("/home/idarraga/analysis/mafalda/MAFOutput_KEK_run050622.root");
 
+        Int_t plotType = __CHARGE;
+        Int_t ndiv = 10;
+        Int_t fulldiv = 10;
+        Int_t rangeYMax = 14;
+
         TCanvas * c10 = new TCanvas();
 
-        TGraphErrors * g1 = GetGraph(T1, "dep", 7, 7);
+        TGraphErrors * g1 = GetGraph(T1, "V100V", ndiv, fulldiv, plotType);
         c10->cd();
         g1->Draw("A*");
         g1->GetXaxis()->SetTitle("depth [mm]");
-        g1->GetYaxis()->SetTitle("mean charge [TOT]");
-        g1->GetYaxis()->SetRangeUser(0., 20);
-        g1->SetMarkerStyle(22);
+        g1->GetYaxis()->SetTitle("Q [TOT]");
+        g1->GetYaxis()->SetRangeUser(0., rangeYMax);
+        g1->SetMarkerStyle(26);
+        g1->SetMarkerSize(2.0);
         g1->SetMarkerColor(kBlue);
 
+        TLatex * lt1 = new TLatex();
+        lt1->DrawLatex(0.050, 9, "backside");
+        lt1->DrawLatex(0.050, 13, "frontside");
+        lt1->DrawLatex(0.06, 4, titleT);
+
+        TLatex * lt2 = new TLatex();
+        lt2->SetTextSize(0.03);
+        lt2->DrawLatex(0.06, 1, "MAFalda, idarraga@cern.ch");
+
+        /////////////////////////////////////////////////////////////////////////////////////
 
         TChain * T2 = new TChain("ShallowAngleFEIX");
-        T2->Add("/home/idarraga/analysis/mafalda/MAFOutput_MPXNtuple_USBPIX_10_f.root");
-
-        TGraphErrors * g2 = GetGraph(T2, "udep", 6, 7);
+        T2->Add("/storage1/idarraga/Testbeam_July_CERN_2011/KEK_15V_4deg/MAFOutput/MAFOutput_MPXNtuple_USBPIXI4_21_15V.root");
+
+        TGraphErrors * g2 = GetGraph(T2, "15V", 9, fulldiv, plotType);
         // fix last point
         c10->cd();
@@ -39 +60 @@
         g2->SetMarkerColor(kRed);
 
-
-
+        /////////////////////////////////////////////////////////////////////////////////////
+
+        TChain * T3 = new TChain("ShallowAngleFEIX");
+        T3->Add("/storage1/idarraga/Testbeam_July_CERN_2011/KEK_7V_4deg/MAFOutput/MAFOutput_MPXNtuple_USBPIXI4_21_7V.root");
+
+        TGraphErrors * g3 = GetGraph(T3, "7V", 9, fulldiv, plotType);
+        // fix last point
+        c10->cd();
+        g3->Draw("*");
+        g3->SetMarkerStyle(20);
+        g3->SetMarkerColor(kGreen);
+
+        /////////////////////////////////////////////////////////////////////////////////////
         TLegend * l1 = new TLegend(0.5, 0.7, 0.9, 0.9);
-        l1->AddEntry(g1, "fully depleted", "P");
-        l1->AddEntry(g2, "under depleted", "P");
+        l1->AddEntry(g1, "100V run50734", "P");
+        l1->AddEntry(g2, "15V run50739", "P");
+        l1->AddEntry(g3, "7V run50743", "P");
         l1->Draw();
+        l1->SetBorderSize(1);
+        l1->SetFillColor(kWhite);
 
 
 }
 
-TGraphErrors * GetGraph(TChain * T, TString sufx, Int_t __NDIV, Int_t fullDiv){
+TGraphErrors * GetGraph(TChain * T, TString sufx, Int_t __NDIV, Int_t fullDiv, Int_t graphType){
 
         gROOT->ProcessLine(".x ~/styles/AtlasStyle.C");
@@ -60 +95 @@
         cout << "Entries = " << T->GetEntries() << endl;
 
-        TCanvas * c1 = new TCanvas("c_"+sufx);
+        TCanvas * c1 = new TCanvas("c_"+sufx, "c_"+sufx);
         c1->Divide(__NDIV/2  + 1 , 2);
 
@@ -67 +102 @@
         vector<double> Qp;
         vector<double> errQp;
+
+        vector<double> entriesV;
+
         double maxQp = 0;
 
@@ -75 +113 @@
 
         int partitionCntr = 0;
-        for (indx = __NDIV-1 ; indx >= 0 ; indx--) {
+        //for (indx = __NDIV-1 ; indx >= 0 ; indx--) {
+        for (indx = 0 ; indx < __NDIV ; indx++) {
 
                 wordDraw = "totPerSegment[";
@@ -83 +122 @@
                 wordDraw += "_";
                 wordDraw += indx;
-                wordDraw += "(100, 0, 100)";
+                wordDraw += "(20, 0, 20)";
 
                 histo = "h_";
@@ -92 +131 @@
                 cout << "word --> " << wordDraw << endl;
                 c1->cd(indx+1);
-                T->Draw(wordDraw, "clusterHeight <= 3");
+                TString cutS = "clusterHeight < 3 && ";
+                cutS += "clusterSize == ";
+                cutS += __NDIV;
+                TCut cutC = cutS;
+                T->Draw(wordDraw, cutC);
                 //"nDeltaRaySoft == 0 && nDeltaRayExplicit == 0");
                 //T->Draw(wordDraw, "nDeltaRayExplicit == 0");
                 TH1F * h = gROOT->FindObject(histo);
-                h->GetXaxis()->SetTitle("mean charge [TOT]");
+                h->GetXaxis()->SetTitle("TOT");
+                h->GetXaxis()->SetTitleSize(0.15);
+                h->GetXaxis()->SetTitleOffset(0.5);
+                h->GetXaxis()->SetLabelSize(0.06);
+
 
                 double meanQp = h->GetMean();
@@ -109 +156 @@
                 a += indx; a += "]=";
                 a +     = TString::Format("%.1f", meanQp);
-                Tl->DrawLatex(meanQp*1.2, h->GetMaximum()/2., a);
+                Tl->DrawLatex(meanQp, h->GetMaximum()*0.9, a);
 
                 TLatex * Tb = new TLatex;
                 Tb->SetTextAlign(12);
-                Tb->SetTextSize(0.1);
+                Tb->SetTextSize(0.15);
                 Int_t nentries = h->GetEntries();
                 cout << "entries in histo = " << nentries << endl;
-                TString b1 = "entries: ";
+                TString b1 = "ent ";
+                entriesV.push_back( nentries );
                 b1 +    = TString::Format("%d", nentries);
-                Tb->DrawLatex(meanQp*1.2, h->GetMaximum()/3., b1);
+                Tb->DrawLatex(meanQp*0.5, h->GetMaximum()/3., b1);
 
                 // for final plot
@@ -144 +192 @@
         c2->cd();
 
-        TGraphErrors * g = new TGraphErrors((int)Qp.size(), &depth[0], &Qp[0], &errDepth[0], &errQp[0]);
+        TGraphErrors * g;
+        if(graphType == __CHARGE){
+                g = new TGraphErrors((int)Qp.size(), &depth[0], &Qp[0], &errDepth[0], &errQp[0]);
+        }
+        else if(graphType == __ENTRIES){
+                g = new TGraphErrors((int)Qp.size(), &depth[0], &entriesV[0], 0x0, 0x0);
+        }
+
 
         return g;