Changeset 124 in Idarraga

Timestamp:

Apr 12, 2011, 9:41:47 AM (13 years ago)

Author:

idarraga

Message:

 

Location:

mafalda

Files:

• AnalysisCore/MediPixAnalysisCoreLinkDef.h (modified) (1 diff)
• MAFTools/MAFTools.cpp (modified) (2 diffs)
• MAFTools/MAFTools.h (modified) (2 diffs)
• MPXAlgo/MediPixAlgo.h (modified) (1 diff)
• PPSResidual/PPSResidual.cpp (modified) (2 diffs)

mafalda/AnalysisCore/MediPixAnalysisCoreLinkDef.h

@@ -46 +46 @@
 #pragma link C++ class hola;
 #pragma link C++ class hola;
+#pragma link C++ class FEI4;
+#pragma link C++ class Test1;
+#pragma link C++ class FEI3Mips;
 
 #endif

mafalda/MAFTools/MAFTools.cpp

@@ -10 +10 @@
 
 #include "MAFTools.h"
+#include "MPXAlgo/Highlighter.h"
+#include "MPXAlgo/MediPixAlgo.h"
 
 #include <set>
 #include <iostream>
-  
+
 #include <TROOT.h>
 #include <TF1.h>
@@ -24 +26 @@
 
 
-  Bool_t LinearRegression(Float_t * slope, Float_t * cut,
-                          Float_t * chisquare_OverDof, 
-                          Float_t initslope, Float_t initycross,
-                          set< pair<Int_t, Int_t> > blobContent){
-    
-    set< pair<Int_t, Int_t> >::iterator setItr;
-    
-    // data -->N points
-    const Int_t N = (const Int_t) blobContent.size();
-    
-    // Fit to ax + b.  N has to be greater or equal to the number of
-    // linearly indepentend functions, in this case 2.
-    if(N < 2) {
-      // fit failed
-      *slope = 0;
-      *chisquare_OverDof = 0.;
-      return false;
-    }
-    
-    Float_t * xd = new Float_t[N];
-    Float_t * yd = new Float_t[N];
-    Int_t cntr = 0;
-    
-    for(setItr = blobContent.begin() ; setItr != blobContent.end() ; setItr++)
-      {
-        xd[cntr] = (Float_t)(*setItr).first;
-        yd[cntr] = (Float_t)(*setItr).second;
-        cntr++;
-      }
-    
-    // get the fit done, linear fit
-    TF1 * f1 = new TF1("f1","[0]*x + [1]",0. , 1.);
-    // enter the guess for the parameters
-    f1->SetParameters(initslope, initycross);
-    TGraph * g1 = new TGraph(N, xd, yd);
-    g1->Fit(f1,"NQ");
-    
-    *chisquare_OverDof = f1->GetChisquare()/(Float_t)f1->GetNDF();
-    *slope = f1->GetParameter(0);
-    *cut = f1->GetParameter(1);
-    
-    delete f1;
-    delete g1;
-    delete [] xd;
-    delete [] yd;
-    
-    return true;
-  }
+Bool_t LinearRegression(Float_t * slope, Float_t * cut,
+                Float_t * chisquare_OverDof,
+                Float_t initslope, Float_t initycross,
+                set< pair<Int_t, Int_t> > blobContent){
 
-  Int_t ConvertXYPositionToInteger(pair<Int_t, Int_t> position, Int_t dimX, Int_t dimY){
+        set< pair<Int_t, Int_t> >::iterator setItr;
 
-    return (position.first)%dimX + (position.second)*dimY;
+        // data -->N points
+        const Int_t N = (const Int_t) blobContent.size();
 
-  }
+        // Fit to ax + b.  N has to be greater or equal to the number of
+        // linearly indepentend functions, in this case 2.
+        if(N < 2) {
+                // fit failed
+                *slope = 0;
+                *chisquare_OverDof = 0.;
+                return false;
+        }
 
-  Int_t ConvertXYPositionToInteger(Int_t positionX, Int_t positionY, Int_t dimX, Int_t dimY){
+        Float_t * xd = new Float_t[N];
+        Float_t * yd = new Float_t[N];
+        Int_t cntr = 0;
 
-    pair<Int_t, Int_t> coordinates (positionX, positionY);
+        for(setItr = blobContent.begin() ; setItr != blobContent.end() ; setItr++)
+        {
+                xd[cntr] = (Float_t)(*setItr).first;
+                yd[cntr] = (Float_t)(*setItr).second;
+                cntr++;
+        }
 
-    return ConvertXYPositionToInteger(coordinates, dimX, dimY);
+        // get the fit done, linear fit
+        TF1 * f1 = new TF1("f1","[0]*x + [1]",0. , 1.);
+        // enter the guess for the parameters
+        f1->SetParameters(initslope, initycross);
+        TGraph * g1 = new TGraph(N, xd, yd);
+        g1->Fit(f1,"NQ");
 
-  }
-  
-  Double_t CalcDistance(Double_t x1, Double_t y1,
-                        Double_t x2, Double_t y2){
+        *chisquare_OverDof = f1->GetChisquare()/(Float_t)f1->GetNDF();
+        *slope = f1->GetParameter(0);
+        *cut = f1->GetParameter(1);
 
-    return TMath::Sqrt((y2 - y1)*(y2 - y1)
-                       +
-                       (x2 - x1)*(x2 - x1));
-    
-  }
+        delete f1;
+        delete g1;
+        delete [] xd;
+        delete [] yd;
 
-  Double_t CalcDistance(std::pair<Double_t, Double_t> p1,
-                        std::pair<Double_t, Double_t> p2){
-    
-    return CalcDistance(p1.first, p1.second,
+        return true;
+}
+
+Int_t ConvertXYPositionToInteger(pair<Int_t, Int_t> position, Int_t dimX, Int_t dimY){
+
+        return (position.first)%dimX + (position.second)*dimY;
+
+}
+
+Int_t ConvertXYPositionToInteger(Int_t positionX, Int_t positionY, Int_t dimX, Int_t dimY){
+
+        pair<Int_t, Int_t> coordinates (positionX, positionY);
+
+        return ConvertXYPositionToInteger(coordinates, dimX, dimY);
+
+}
+
+Double_t CalcDistance(Double_t x1, Double_t y1,
+                Double_t x2, Double_t y2){
+
+        return TMath::Sqrt((y2 - y1)*(y2 - y1)
+                        +
+                        (x2 - x1)*(x2 - x1));
+
+}
+
+Double_t CalcDistance(std::pair<Double_t, Double_t> p1,
+                std::pair<Double_t, Double_t> p2){
+
+        return CalcDistance(p1.first, p1.second,
                         p2.first, p2.second);
-  }
-  
+}
 
+Double_t CalcPerpDistanceToLine(Float_t slope, Float_t cut, pair<Int_t, Int_t> point){
+        return CalcPerpDistanceToLine(slope, cut, point.first, point.second);
+}
+
+Double_t CalcPerpDistanceToLine(Float_t slope, Float_t cut, Int_t px, Int_t py){
+
+        Float_t cutPrim = 0.;
+        Float_t slopePrim = -1./slope;
+
+        // Cut of the perpendicular line
+        cutPrim  = py;
+        cutPrim -= slopePrim * px;
+
+        // intersection of the two lines
+        Float_t interx = (cutPrim - cut)/(slope - slopePrim);
+        Float_t intery = interx * slope + cut;
+
+        return CalcDistance(px, py, interx, intery);
+}
+
+void DrawLine(MediPixAlgo * algo, Float_t x1, Float_t y1, Float_t x2, Float_t y2, Int_t wd = 2, Int_t st = 2, EColor color = kBlack) {
+
+        Highlighter * lineD = new Highlighter(x1, y1, x2, y2, "line", algo);
+        lineD->SetLineColor(color);
+        lineD->SetLineWidth(wd);
+        lineD->SetLineStyle(st);
+        algo->PullToStoreGateAccess(lineD, MPXDefs::DO_NOT_SERIALIZE_ME);
+
+}
+
+std::string DumpWordInBinary(UInt_t word) {
+
+        string wordS;
+
+        UInt_t oneBit = 0x0;
+        Int_t sizeOfWord = sizeof(word)*8; // number of bits
+
+        for (int i = 0 ; i < sizeOfWord ; i++) { // number of bits
+
+                oneBit = word & 0x1;
+
+                if(oneBit == 0x1) { wordS.append(1,'1'); }
+                else { wordS.append(1,'0'); }
+
+                word = word >> 1;
+        }
+
+        reverse(wordS.begin(),wordS.end());
+
+        return wordS;
+}
 
 }

mafalda/MAFTools/MAFTools.h

@@ -11 +11 @@
 #include <set>
 #include <iostream>
+#include <vector>
 
 #include <TROOT.h>
+
+class MediPixAlgo;
+
+using namespace std;
 
 namespace MAFTools {
@@ -30 +35 @@
                         Double_t, Double_t);
   
+  Double_t CalcPerpDistanceToLine(Float_t, Float_t, pair<Int_t, Int_t>);
+  Double_t CalcPerpDistanceToLine(Float_t, Float_t, Int_t, Int_t);
+
+
+  /* Drawing */
+  void DrawLine(MediPixAlgo *, Float_t, Float_t, Float_t, Float_t, Int_t, Int_t, EColor);
+
+  /* Word handling */
+  std::string DumpWordInBinary(UInt_t word);
 
 }

mafalda/MPXAlgo/MediPixAlgo.h

@@ -125 +125 @@
   /* Other Gets */
   TApplication * GetApplication();
-  // get the tree corresponding to this algo
+  // Get the TTree corresponding to this algo
   TTree * getMyTree();
+  TTree * GetMyTree(){return getMyTree();};
+  // Get the ROOTFile;
   TFile * getMyROOTFile();
+  TFile * GetMyROOTFile(){return getMyROOTFile();};
+  // output
   MSG::Endreq endreq;
 

mafalda/PPSResidual/PPSResidual.cpp

@@ -7 +7 @@
 #include "MAFTools/MAFTools.h"
 #include "TVector3.h"
+#include "TF1.h"
 
 #include <vector>
@@ -17 +18 @@
 PPSResidual::PPSResidual(){
 
+        TF1 * f1 = new TF1("hola","x",0., 0.);
 
 }