Fork: 0
ibezshyi / TB_Chris
Transfer to URL with SHA Find file
Newer
Older
TB_Chris / TbUT / scripts / .svn / text-base / ClusterWithTrackAna.C.svn-base
@iaro iaro on 9 May 2016 47 KB first attempt of automated anal
Raw Blame History 
  1. #define ClusterWithTrackAna_cxx
  2. //#include "ClusterWithTrackAna_Inputs.h"
  3. #include "ClusterWithTrackAna.h"
  4. #include <TH2.h>
  5. #include <TStyle.h>
  6. #include <TLatex.h>
  7. #include <TCanvas.h>
  8. #include <iostream>
  9. #include <fstream>
  10.  
  11. void addGraphics(TH1 *h, int iCol = 1, TString XTitle="", TString YTitle="")
  12. {
  13.   h->SetXTitle(XTitle);
  14.   h->SetYTitle(YTitle);
  15.   h->SetStats(kFALSE);
  16.   h->SetLineColor(iCol);
  17.   h->SetMarkerColor(iCol);
  18.   h->SetMinimum(0.0);
  19.   h->SetMaximum(1.2*h->GetMaximum());
  20.   h->SetTitleSize(0.1);
  21.   
  22.   //h->SetLineColor(kBlack);
  23.   h->SetMarkerSize(0.7);
  24.   h->SetMarkerStyle(20);
  25.   h->GetXaxis()->SetTitleOffset(1.0);  
  26.   h->GetYaxis()->SetTitleOffset(1.2);
  27.   h->GetXaxis()->SetTitleSize(0.045);  
  28.   h->GetYaxis()->SetTitleSize(0.045);
  29.   h->GetXaxis()->SetLabelSize(0.04);  
  30.   h->GetYaxis()->SetLabelSize(0.04);  
  31.   h->SetNdivisions(505,"X");
  32.   h->SetNdivisions(505,"Y");
  33.   h->SetLineWidth(2);
  34. }
  35.  
  36. void addGraphics(TH2 *h, int iCol = 1, TString XTitle="", TString YTitle="")
  37. {
  38.   //float bw = h->GetBinWidth(1);
  39.   h->SetXTitle(XTitle);
  40.   h->SetYTitle(YTitle);
  41.   h->SetStats(kFALSE);
  42.   h->SetLineColor(iCol);
  43.   h->SetMarkerColor(iCol);
  44.   h->SetMinimum(0.0);
  45.   h->SetMaximum(1.2*h->GetMaximum());
  46.   h->SetTitleSize(0.1);
  47.   
  48.   //h->SetLineColor(kBlack);
  49.   //h->SetMarkerSize(0.7);
  50.   //h->SetMarkerStyle(20);
  51.   h->GetXaxis()->SetTitleOffset(1.0);  
  52.   h->GetYaxis()->SetTitleOffset(1.2);
  53.   h->GetXaxis()->SetTitleSize(0.045);  
  54.   h->GetYaxis()->SetTitleSize(0.045);
  55.   h->GetXaxis()->SetLabelSize(0.04);  
  56.   h->GetYaxis()->SetLabelSize(0.04);  
  57.   h->SetNdivisions(505,"X");
  58.   h->SetNdivisions(505,"Y");
  59.   h->SetLineWidth(2);
  60. }
  61.  
  62.  
  63. void ClusterWithTrackAna::Loop()
  64. {
  65. //   In a ROOT session, you can do:
  66. //      Root > .L ClusterWithTrackAna.C
  67. //      Root > ClusterWithTrackAna t
  68. //      Root > t.GetEntry(12); // Fill t data members with entry number 12
  69. //      Root > t.Show();       // Show values of entry 12
  70. //      Root > t.Show(16);     // Read and show values of entry 16
  71. //      Root > t.Loop();       // Loop on all entries
  72. //
  73.  
  74. //     This is the loop skeleton where:
  75. //    jentry is the global entry number in the chain
  76. //    ientry is the entry number in the current Tree
  77. //  Note that the argument to GetEntry must be:
  78. //    jentry for TChain::GetEntry
  79. //    ientry for TTree::GetEntry and TBranch::GetEntry
  80. //
  81. //       To read only selected branches, Insert statements like:
  82. // METHOD1:
  83. //    fChain->SetBranchStatus("*",0);  // disable all branches
  84. //    fChain->SetBranchStatus("branchname",1);  // activate branchname
  85. // METHOD2: replace line
  86. //    fChain->GetEntry(jentry);       //read all branches
  87. //by  b_branchname->GetEntry(ientry); //read only this branch
  88.    if (fChain == 0) return;
  89.  
  90.    Int_t nentries = fChain->GetEntriesFast();
  91.  
  92.    TString m_board2 = m_board;
  93.    m_board2 = m_board2.ReplaceAll("_All","");
  94.    m_board2 = m_board2.ReplaceAll("_Full","");
  95.    m_board2 = m_board2.ReplaceAll("_v7","");
  96.    TString f_out = m_fileOutdir + plotdir + "/AnalysisOutput_" + m_board2 + "_" + m_bias + "_" + m_sector + ".root";
  97.    if(m_angle != "0"){
  98.      f_out			= m_fileOutdir + plotdir + "/AnalysisOutput_" + m_board2 + "_" + m_bias + "_" + m_sector + "_" + m_angle + ".root";
  99.    }
  100.    cout << "Will write out file: " << f_out << endl;
  101.    
  102.  
  103.  
  104.    fout = new TFile(f_out,"RECREATE");
  105.  
  106.    TH1F* h0 = new TH1F("h0","#DeltaX between strip hit and track projection (strips)",101,-50.5,50.5);
  107.    TH1F* h1 = new TH1F("h1","#DeltaX",800,-2.0,2.0);
  108.    TH2F* h1vsx = new TH2F("h1vsx","#DeltaX vs X",50,-5,5,100,-0.2,0.2);
  109.    TH1F* h1mpa = new TH1F("h1mpa","#DeltaX",400,-20.0,20.0);
  110.    TH1F* h1mbpa = new TH1F("h1mbpa","#DeltaX",400,-20.0,20.0);
  111.    TH1F* h1fpa = new TH1F("h1fpa","#DeltaX",400,-20.0,20.0);
  112.    TH1F* h1fbpa = new TH1F("h1fbpa","#DeltaX",400,-20.0,20.0);
  113.  
  114.    TH1F* h1mpa1 = new TH1F("h1mpa1","#DeltaX",800,-20.0,20.0);
  115.    TH1F* h1mpa2 = new TH1F("h1mpa2","#DeltaX",800,-20.0,20.0);
  116.    TH1F* h1mpa3 = new TH1F("h1mpa3","#DeltaX",800,-20.0,20.0);
  117.    TH1F* h1mpa4 = new TH1F("h1mpa4","#DeltaX",800,-20.0,20.0);
  118.    TH1F* h1mpa5 = new TH1F("h1mpa5","#DeltaX",800,-20.0,20.0);
  119.  
  120.    TH1F* h1mpaL = new TH1F("h1mpaL","Strip# of Missed Hit in Lower PA region",512,0,512);
  121.    TH1F* h1mpaU = new TH1F("h1mpaU","Strip# of Missed Hit in Upper PA region",512,0,512);
  122.  
  123.  
  124.    TH1F* h1s = new TH1F("h1s","Seed strip of cluster",512,0.0,512.0);
  125.    TH1F* h1a = new TH1F("h1a","#DeltaX, 1 strip",800,-2.0,2.0);
  126.    TH1F* h1b = new TH1F("h1b","#DeltaX, 2 strip",800,-2.0,2.0);
  127.    TH1F* h1z = new TH1F("h1z","#DeltaX, low Charge",800,-2.0,2.0);
  128.    TH1F* h1w = new TH1F("h1w","#DeltaX",20000,-100.0,100.0);
  129.    TH1F* h1wY = new TH1F("h1wY","#DeltaX",20000,-100.0,100.0);
  130.    TH2F* h2 = new TH2F("h2","X_{DUT} vs X_{trk}",100,-8.0,8.0,100,-8,8);
  131.    TH2F* h3 = new TH2F("h3","Y_{trk} vs X_{trk}, with cluster",640,-8,8.0,640,-8,8);
  132.    TH2F* h3a = new TH2F("h3a","Y_{trk} vs X_{trk} in fiducial",640,-8,8.0,640,-8,8);
  133.    TH2F* h3b = new TH2F("h3b","Y_{trk} vs X_{trk}, with missed cluster",640,-8,8.0,640,-8,8);
  134.    TH2F* h3c = new TH2F("h3c","Y_{trk} vs X_{trk}, with found cluster",640,-8,8.0,640,-8,8);
  135.    TH1F* h4 = new TH1F("h4","Detector strip # of cluster with track",512,0.0,512);
  136.    TH1F* h4a = new TH1F("h4a","Electonic Strip # of cluster with track",512,0.0,512);
  137.    TH1F* h4b = new TH1F("h4b","Strip # of cluster with track",1024,0.0,1024);
  138.    TH1F* h4c = new TH1F("h4c","Strip # of cluster with track",1024,0.0,1024);
  139.    TH1F* h5 = new TH1F("h5","#theta_{X}",500,-5.0,5.0);
  140.    TH1F* h6 = new TH1F("h6","#theta_{Y}",500,-5.0,5.0);
  141.    TH1F* h5c = new TH1F("h5c","#theta_{X}",500,-5.0,5.0);
  142.    TH1F* h6c = new TH1F("h6c","#theta_{Y}",500,-5.0,5.0);
  143.    TH1F* h5a = new TH1F("h5a","X position of track",400,-10.0,10.0);
  144.    TH1F* h6a = new TH1F("h6a","Y position of track",400,-10.0,10.0);
  145.    TH1F* h5b = new TH1F("h5b","X position of matched cluster",400,-10.0,10.0);
  146.    TH1F* h6b = new TH1F("h6b","Y position of matched cluster",400,-10.0,10.0);
  147.  
  148.    TProfile *h8 = new TProfile("h8","#DeltaX vs #theta_{trk}",100,-5,5,-1.0,1.0);
  149.    TProfile *h9 = new TProfile("h9","#DeltaX vs Y_{trk} at DUT",1600,-8,8,-1.0,1.0);
  150.    TProfile *h9a = new TProfile("h9a","#DeltaX vs X_{trk} at DUT",20,-5,5,-1.0,1.0);
  151.  
  152.    TProfile *h10a = new TProfile("h10a","<ADC> vs strip",512,0,512,0.0,1000.0);
  153.    TProfile *h10b = new TProfile("h10b","<ADC> vs strip",512,0,512,0.0,1000.0);
  154.    TProfile *h10c = new TProfile("h10c","<ADC> vs strip",512,0,512,0.0,1000.0);
  155.    TProfile *h10d = new TProfile("h10d","<ADC> vs Y_{trk}",100,-8,8,0.0,1000.0);
  156.    TProfile *h10e = new TProfile("h10e","<ADC> vs X_{trk}",100,-8,8,0.0,1000.0);
  157.    TH1F* h11n = new TH1F("h11n","Strip # of matched cluster",512,0.0,512);
  158.    TH1F* h11d = new TH1F("h11d","Strip # of track",512,0.0,512);
  159.    h11n->Sumw2();
  160.    h11d->Sumw2();
  161.    
  162.    TH1F* h12 = new TH1F("h12","Y position of matched cluster",400,-10.0,10.0); h12->Sumw2();
  163.    TH1F* h12c = new TH1F("h12c","X position of matched cluster",400,-10.0,10.0);h12c->Sumw2();
  164.    TH1F* h12a = new TH1F("h12a","Y position of track",1600,-10.0,10.0); h12a->Sumw2();
  165.    TH1F* h12b = new TH1F("h12b","Y position of track",1600,-10.0,10.0); h12b->Sumw2();
  166.  
  167.    TH1F* h12dn = new TH1F("h12dn","X position of track",200,-10.0,10.0); h12dn->Sumw2();
  168.    TH1F* h12en = new TH1F("h12en","X position of track",200,-10.0,10.0);h12en->Sumw2();
  169.    TH1F* h12fn = new TH1F("h12fn","X position of track",200,-10.0,10.0);h12fn->Sumw2();
  170.    TH1F* h12gn = new TH1F("h12gn","X position of track",200,-10.0,10.0);h12gn->Sumw2();
  171.    TH1F* h12dd = new TH1F("h12dd","X position of track",200,-10.0,10.0);h12dd->Sumw2();
  172.    TH1F* h12ed = new TH1F("h12ed","X position of track",200,-10.0,10.0);h12ed->Sumw2();
  173.    TH1F* h12fd = new TH1F("h12fd","X position of track",200,-10.0,10.0);h12fd->Sumw2();
  174.    TH1F* h12gd = new TH1F("h12gd","X position of track",200,-10.0,10.0);h12gd->Sumw2();
  175.  
  176.    TH1F* h12hn = new TH1F("h12hn","X position of track",50,-0.5,0.5);h12hn->Sumw2();
  177.    TH1F* h12in = new TH1F("h12in","X position of track",50,-0.5,0.5);h12in->Sumw2();
  178.    TH1F* h12jn = new TH1F("h12jn","X position of track",50,-0.5,0.5);h12jn->Sumw2();
  179.    TH1F* h12kn = new TH1F("h12kn","X position of track",50,-0.5,0.5);h12kn->Sumw2();
  180.    TH1F* h12hd = new TH1F("h12hd","X position of track",50,-0.5,0.5);h12hd->Sumw2();
  181.    TH1F* h12id = new TH1F("h12id","X position of track",50,-0.5,0.5);h12id->Sumw2();
  182.    TH1F* h12jd = new TH1F("h12jd","X position of track",50,-0.5,0.5);h12jd->Sumw2();
  183.    TH1F* h12kd = new TH1F("h12kd","X position of track",50,-0.5,0.5);h12kd->Sumw2();
  184.  
  185.    TProfile *h12m = new TProfile("h12m","<ADC> vs interstrip pos",50,-0.5,0.5,0.0,1000.0);
  186.    TProfile *h12n = new TProfile("h12n","<ClusterSize> vs interstrip pos",50,-0.5,0.5,0.0,1000.0);
  187.  
  188.    TH1F* h12on = new TH1F("h12on","dist of track to cutout",700,-2.0,5.0);h12on->Sumw2();
  189.    TH1F* h12od = new TH1F("h12od","dist of track to cutout",700,-2.0,5.0);h12od->Sumw2();
  190.  
  191.    TH1F *h13mpa = new TH1F("h13mpa","#Clusters/Ev",10,-0.5,9.5);
  192.    TH1F *h13mbpa = new TH1F("h13mbpa","#Clusters/Ev",10,-0.5,9.5);
  193.    TH1F *h13fpa = new TH1F("h13fpa","#Clusters/Ev",10,-0.5,9.5);
  194.    TH1F *h13fbpa = new TH1F("h13fbpa","#Clusters/Ev",10,-0.5,9.5);   
  195.  
  196.    TH1F* h14a = new TH1F("h14a","Nom Strip, no matched cluster",512,0.0,512);
  197.    TH1F* h14b = new TH1F("h14b","Nom Strip, matched cluster",512,0.0,512);
  198.  
  199.    TH1F *h15a = new TH1F("h15a","Chisq of tracks with matched DUT hit",200,0.0,200.0);
  200.    TH1F *h15b = new TH1F("h15b","TDC time with matched DUT hit",12,0.0,12.0);
  201.    TH1F *h15c = new TH1F("h15c","Time Between Trigger and Track",5,0.0,5.0);
  202.    TH1F *h15d = new TH1F("h15d","TDC time of all clusters",12,0.0,12.0);
  203.    TH1F *h15e = new TH1F("h15e","TDC time of good clusters",12,0.0,12.0);
  204.  
  205.    TH1F *h16a = new TH1F("h16a","Chisq of tracks without matched DUT hit",200,0.0,200.0);
  206.    TH1F *h16b = new TH1F("h16b","TDC time without matched DUT hit",12,0.0,12.0);
  207.    TH1F *h16c = new TH1F("h16c","Time Between Trigger and Track",5,0.0,5.0);
  208.    h15c->Sumw2();
  209.    h16c->Sumw2();
  210.    
  211.    TH1F* h17 = new TH1F("h17","Q_{L} / (Q_{L}+Q_{R})",50,0.0,1.0);
  212.    TH2F* h17a = new TH2F("h17a","Q_{L} / (Q_{L}+Q_{R}) vs Interstrip Pos",20,-0.5,0.5,120,-0.1,1.1);
  213.    TProfile* h17b = new TProfile("h17b","Q_{R} / (Q_{L}+Q_{R}) vs Interstrip Pos",20,-0.5,0.5,-0.1,1.1);
  214.  
  215.    TH1F* hcmpa = new TH1F("hcmpa","Cluster charge",100,0.0,1000.0);
  216.    TH1F* hcmbpa = new TH1F("hcmbpa","Cluster charge",100,0.0,1000.0);
  217.    TH1F* hcfpa = new TH1F("hcfpa","Cluster charge",100,0.0,1000.0);
  218.    TH1F* hcfbpa = new TH1F("hcfbpa","Cluster charge",100,0.0,1000.0);
  219.  
  220.    TH1F* h18a = new TH1F("h18a","Seed/Total charge",100,0.0,2.0);
  221.  
  222.    TH2F* h18b = new TH2F("h18b","#DeltaX vs Interstrip Pos (1 strip Clu)",100,-0.5,0.5,100,-0.2,0.2);
  223.    TH2F* h18c = new TH2F("h18c","#DeltaX vs Interstrip Pos (2 strip Clu)",100,-0.5,0.5,100,-0.2,0.2);
  224.    TH2F* h18d = new TH2F("h18d","Seed/Total charge vs Interstrip Pos",100,-0.5,0.5,100,0.0,1.1);
  225.  
  226.    TH1F* hcTrkXY[16];
  227.    for(int i=0; i<4; i++){
  228.      for(int j=0; j<4; j++){
  229.        int iblock = 4*i+j; 
  230.        TString hname = Form("hcTrkXY_%d",iblock+1);     
  231.        TString htitle = Form("Cluster charge, block %d",iblock+1);     
  232.        hcTrkXY[iblock]= new TH1F(hname,htitle,100,0.0,1000.0);
  233.      }
  234.    }
  235.    
  236.    TH1F *hlandau[512];
  237.    for(int i=0; i<nChan; i++){
  238.      hlandau[i] = new TH1F(Form("hlandau_%d",i),"Cluster charge",100,0.0,1000.0);
  239.    }
  240.    
  241.  
  242.  
  243.    TH1F* hcAll = new TH1F("hcAll","Cluster charge",100,0.0,1000.0);
  244.    TH1F* hcTrk = new TH1F("hcTrk","Cluster charge",100,0.0,1000.0);
  245.    TH1F* hcTrkCorr = new TH1F("hcTrkCorr","Cluster charge",100,0.0,1000.0);
  246.    TH1F* hcTrkSeed = new TH1F("hcTrkSeed","Cluster charge",100,0.0,1000.0);
  247.    TH1F* hcTrkSeed1 = new TH1F("hcTrkSeed1","Cluster charge",100,0.0,1000.0);
  248.    TH1F* hcTrkSeed2 = new TH1F("hcTrkSeed2","Cluster charge",100,0.0,1000.0);
  249.    //TH1F* hcTrkY1 = new TH1F("hcTrkY1","Cluster charge",100,0.0,1000.0);
  250.    //TH1F* hcTrkY2 = new TH1F("hcTrkY2","Cluster charge",100,0.0,1000.0);
  251.    //TH1F* hcTrkY3 = new TH1F("hcTrkY3","Cluster charge",100,0.0,1000.0);
  252.    //TH1F* hcTrkY4 = new TH1F("hcTrkY4","Cluster charge",100,0.0,1000.0);
  253.    TH1F* hcTrkNoTop = new TH1F("hcTrkNoTop","Cluster charge",100,0.0,1000.0);
  254.    TH1F* hcTrkLowX = new TH1F("hcTrkLowX","Cluster charge",100,0.0,1000.0);
  255.    TH1F* hcTrkHighX = new TH1F("hcTrkHighX","Cluster charge",100,0.0,1000.0);
  256.    TH1F* hcTrka = new TH1F("hcTrka","Cluster charge, limited region",100,0.0,1000.0);
  257.    TH1F* hcTrk1 = new TH1F("hcTrk1","Cluster charge, Size = 1",100,0.0,1000.0);
  258.    TH1F* hcTrk2 = new TH1F("hcTrk2","Cluster charge, Size = 2",100,0.0,1000.0);
  259.    TProfile *h2p = new TProfile("h2p","Cluster Charge vs TDC time",12,0,12,100,1000);
  260.  
  261.  
  262.    TH2F* h31a = new TH2F("h31a","Y_{trk} vs X_{trk}, with cluster",64,-8,8.0,64,-8,8);
  263.    TH2F* h31b = new TH2F("h31b","Y_{trk} vs X_{trk}, with cluster, low ADC",64,-8,8.0,64,-8,8);
  264.    TH1F* h32a = new TH1F("h32a","Strip # of cluster with track",512,0.0,512);
  265.    TH1F* h32b = new TH1F("h32b","Strip # of cluster with track & low ADC",512,0.0,512);
  266.    TH1F* h33a = new TH1F("h33a","Y position of matched cluster",200,-10.0,10.0);
  267.    TH1F* h33b = new TH1F("h33b","X position of matched cluster",200,-10.0,10.0);
  268.    TH1F* h34 = new TH1F("h34","Y position of matched cluster & low ADC",200,-10.0,10.0);
  269.  
  270.    TH1F *hnoise = new TH1F("hnoise","Noise in connected channels",100,-200,200);
  271.    TH1F *hnoiseChan = new TH1F("hnoiseChan","Noise in connected channels",200,0,200);
  272.    TH1F *hnoisePerChannel = new TH1F("hnoisePerChannel","Noise",512,0,512);
  273.  
  274.    TH1F* h35 = new TH1F("h35","No. clusters / event",50,0.0,50.0);
  275.    
  276.    TH1F *h41[10];
  277.    TH1F *h42[10];
  278.    TH1F *h43[10];
  279.    TH1F *h44[10];
  280.    for(int i=0; i<10; i++){
  281.      h41[i] = new TH1F(Form("h41_%d",i),Form("ADC_{L}-ADC_{R}, Bin %d, Odd Ch",i),80,-400,400);
  282.      h42[i] = new TH1F(Form("h42_%d",i),Form("ADC_{L}-ADC_{R}, Bin %d, Odd Ch",i),80,-400,400);
  283.      h43[i] = new TH1F(Form("h43_%d",i),Form("ADC_{L}-ADC_{R}, Bin %d, Odd Ch",i),160,-400,400);
  284.      h44[i] = new TH1F(Form("h44_%d",i),Form("ADC_{L}-ADC_{R}, Bin %d, Odd Ch",i),160,-400,400);
  285.    }
  286.  
  287.    for(int i=0; i<nChan; i++){
  288.      hnoisePerChannel->Fill(i+0.5,noise[i]);
  289.    }
  290.    
  291.  
  292.  
  293.    //------------------------------------------------------------------------------------------------------
  294.    
  295.    //-----------------------------
  296.    // Prepare DUT (Alignment, etc)
  297.    //-----------------------------
  298.  
  299.    PrepareDUT();
  300.    //return;
  301.    float biasVal = atof(m_bias);
  302.    cout << "chargeCorrSlopeOdd, chargeCorrSlopeEven = " << chargeCorrSlopeOdd << " " << chargeCorrSlopeEven << " " 
  303.         << m_bias << " " << biasVal << endl;
  304.    
  305.  
  306.    std::ofstream myfile;
  307.    if(writeEventsWithMissinhHitsToFile){
  308.      myfile.open("MissingDUTHits.dat");
  309.    }
  310.  
  311.    int iChan = nChan;
  312.    double nomStrip = 0, detStrip = 0;
  313.    int nPrint = 0;
  314.    double dxh[10];
  315.    int njump = 10000;
  316.    if(nentries > 100000) njump = 50000;
  317.    Long64_t nbytes = 0, nb = 0;
  318.    cout << "Begin loop over " << nentries << " events" << endl;
  319.    for (Long64_t jentry=0; jentry<nentries;jentry++) {
  320.       Long64_t ientry = LoadTree(jentry);
  321.       if (ientry < 0) break;
  322.       nb = fChain->GetEntry(jentry);   nbytes += nb;
  323.       if(jentry%njump==0) cout << "====> At entry = " << jentry << endl;
  324.       if(n_tp3_tracks != 1) continue;
  325.  
  326.       //==========================================
  327.       // Example here of acccessing CMS strip data
  328.       //==========================================
  329.       if(fCMS){
  330.         for(int k=0;k<iChan;k++){
  331.           if(fCMS->cmsData[k]!=0) hnoise->Fill(fCMS->cmsData[k]);
  332.         }      
  333.       }
  334.  
  335.       h35->Fill(clusterNumberPerEvent);
  336.       if(clusterNumberPerEvent>10) continue;
  337.       
  338.       
  339.       // Loop over TPIX tracks in event
  340.       for(int k=0; k<n_tp3_tracks; k++){
  341.         if(dtime > trackTriggerTimeDiffCut) continue;
  342.         double x_trk = vec_trk_tx->at(k)*z_DUT+vec_trk_x->at(k);
  343.         double y_trk = vec_trk_ty->at(k)*z_DUT+vec_trk_y->at(k);
  344.         
  345.         transformTrackToDUTFrame(k, x_trk, y_trk, nomStrip, detStrip);
  346.  
  347.         //if(isInCutoutRegion(x_trk, y_trk)) continue;
  348.         double distToCutout = DistToCutoutRegion(x_trk, y_trk);
  349.         bool awayFromCutout = distToCutout > minDistFromHole;
  350.  
  351.         int closestStrip = nomStrip;
  352.         double fracStrip = nomStrip - closestStrip - 0.5;
  353.  
  354.         double tx = 1000*vec_trk_tx->at(k);
  355.         double ty = 1000*vec_trk_ty->at(k);
  356.  
  357.         h5->Fill(tx);
  358.         h6->Fill(ty);
  359.  
  360.         // Test cuts here, e.g.  to reject regions
  361.         //if(nomStrip<422 || nomStrip>442) continue;
  362.         //if(y_trk>2.4) continue;
  363.         //if(nomStrip<420 || nomStrip>445) continue;
  364.         
  365.  
  366.         bool goodRegion = true;
  367.         for(int id = 0; id<nDeadRegion; id++){
  368.           if(x_trk>=deadRegionLo[id]  && x_trk<=deadRegionHi[id]) goodRegion = false;  
  369.         }
  370.         if(!goodRegion) continue;
  371.  
  372.         h5a->Fill(x_trk);
  373.         h6a->Fill(y_trk);
  374.         bool goodTrack = false;
  375.         bool inFiducialX = false;
  376.         bool inFiducialY = false;
  377.  
  378.         if(x_trk>xMin && x_trk<xMax) inFiducialX = true;          
  379.         if(y_trk>yMin && y_trk<yMax) inFiducialY = true;          
  380.         bool inFiducial = inFiducialX && inFiducialY;
  381.         inFiducial = inFiducial && (x_trk<xLeftHole || x_trk>xRightHole);
  382.         
  383.         if(tx>txMin && tx<txMax && ty>tyMin && ty<tyMax) goodTrack = true;        
  384.         bool goodTime =  (clustersTDC >= tdcLo && clustersTDC <= tdcHi);
  385.  
  386.         if(goodTrack && goodTime && inFiducial) h12od->Fill(distToCutout);
  387.  
  388.         if(goodTrack && goodTime && inFiducial && awayFromCutout) {
  389.           h3a->Fill(x_trk,y_trk);
  390.           h12a->Fill(y_trk);
  391.           h12dd->Fill(x_trk);
  392.           if(y_trk>yInt1[0]&&y_trk<yInt1[1]) h12ed->Fill(x_trk);
  393.           if(y_trk>yInt2[0]&&y_trk<yInt2[1]) h12fd->Fill(x_trk);
  394.           if(y_trk>yInt3[0]&&y_trk<yInt3[1]) h12gd->Fill(x_trk);
  395.           h12hd->Fill(fracStrip);
  396.           if(y_trk>yInt1[0]&&y_trk<yInt1[1] && x_trk>-3.5&&x_trk<2.5) h12id->Fill(fracStrip);
  397.           if(y_trk>yInt2[0]&&y_trk<yInt2[1] && x_trk>-3.5&&x_trk<2.5) h12jd->Fill(fracStrip);
  398.           if(y_trk>yInt3[0]&&y_trk<yInt3[1] && x_trk>-3.5&&x_trk<2.5) h12kd->Fill(fracStrip);
  399.         }
  400.         
  401.         bool foundHit = false;
  402.         bool foundHitNoFid = false;
  403.         double x_trk0 = x_trk;
  404.         for(int j=0; j<min(clusterNumberPerEvent,10); j++){
  405.           dxh[j] = -999;
  406.           if(clustersPosition[j] < 0.1) continue;
  407.           if(polarity*clustersCharge[j] < 0.2*kClusterChargeMin) continue;
  408.           double x_dut = getDUTHitPosition(j);
  409.           x_trk = x_trk0;
  410.        
  411.           int iPeak = 1;
  412.           if(clustersSeedPosition[j]%2==0) iPeak = 0;
  413.           
  414.           h1s->Fill(clustersSeedPosition[j]);
  415.           if(clustersSeedPosition[j]<iLo-5 || clustersSeedPosition[j]>iHi+5) continue;
  416.  
  417.           double dx = x_dut - x_trk;
  418.           dxh[j] = dx;
  419.  
  420.           if(inFiducial && goodTime && awayFromCutout) h12->Fill(y_trk); 
  421.           if(goodTrack && inFiducial && fabs(dx)<dxWin && awayFromCutout) {
  422.             h2p->Fill(clustersTDC+0.1,polarity*clustersCharge[j]);
  423.             h15d->Fill(clustersTDC);          
  424.             if(goodTime) h15e->Fill(clustersTDC);          
  425.           }
  426.           
  427.  
  428.           if(goodTrack && inFiducial && goodTime && fabs(dx)<dxWin){
  429.             h3->Fill(x_trk,y_trk);
  430.           }
  431.           
  432.           if(goodTrack && inFiducial && goodTime && fabs(dx)<dxWin) foundHitNoFid = true;
  433.  
  434.  
  435.           if(goodTrack && inFiducial && goodTime && awayFromCutout) {
  436.             hcAll->Fill(polarity*clustersCharge[j]);
  437.             h2->Fill(x_trk, x_dut);
  438.             h1->Fill(dx);
  439.             if(clustersSize[j]==1) h1a->Fill(dx);
  440.             if(clustersSize[j]==2) h1b->Fill(dx);
  441.             h1w->Fill(dx);
  442.             if(y_trk>2.5) h1wY->Fill(dx);
  443.             if(polarity*clustersCharge[j] < 250) h1z->Fill(dx);
  444.             h11d->Fill(detStrip);        
  445.             
  446.             if(fabs(dx)<dxWin) {
  447.               int ichan = clustersSeedPosition[j];
  448.               h4c->Fill(clustersSeedPosition[j]);
  449.               if(ichan>=0 && ichan<=511){
  450.                 hlandau[ichan]->Fill(polarity*clustersCharge[j]);
  451.               }
  452.               
  453.               hnoiseChan->Fill(noise[ichan]);
  454.               h18a->Fill(clustersSeedCharge[j]/clustersCharge[j]);
  455.               foundHit = true;
  456.               if(clustersSize[j]==1) h18b->Fill(fracStrip,dx);
  457.               if(clustersSize[j]==2) h18c->Fill(fracStrip,dx);
  458.               h18d->Fill(fracStrip,clustersSeedCharge[j]/clustersCharge[j]);
  459.               
  460.               h12m->Fill(fracStrip,polarity*clustersCharge[j]);
  461.               h12n->Fill(fracStrip,clustersSize[j]);
  462.               h1vsx->Fill(x_trk,dx);
  463.               if(y_trk>yMid&&y_trk<yMax) h10a->Fill(clustersPosition[j],polarity*clustersCharge[j]);
  464.               if(y_trk>yMin&&y_trk<yMid) h10b->Fill(clustersPosition[j],polarity*clustersCharge[j]);
  465.               if(y_trk>yHi2&&y_trk<yMax) h10c->Fill(clustersPosition[j],polarity*clustersCharge[j]);
  466.               h10d->Fill(y_trk,polarity*clustersCharge[j]);
  467.               h10e->Fill(x_trk,polarity*clustersCharge[j]);
  468.               double chleft = polarity*clustersCharge1StripLeft[j];
  469.               double chright = polarity*clustersCharge1StripRight[j];
  470.               double rc = -999;
  471.               if(detStrip <= clustersSeedPosition[j]) {
  472.                 if(chleft>0) rc = (chleft/(chleft+clustersSeedCharge[j]));              
  473.               }else{ 
  474.                 if(chright>0) rc = (clustersSeedCharge[j]/(chright+clustersSeedCharge[j]));
  475.               }
  476.               if( rc>=0 ) {
  477.                 h17->Fill(rc);
  478.                 h17a->Fill(fracStrip,rc);
  479.                 h17b->Fill(fracStrip,rc);
  480.               }
  481.               h8->Fill(1000*vec_trk_tx->at(k),dx);    
  482.               h9->Fill(y_trk,dx);
  483.               h9a->Fill(x_trk,dx);
  484.               h5b->Fill(x_trk);
  485.               h6b->Fill(y_trk);
  486.               h5c->Fill(tx);
  487.               h6c->Fill(ty);
  488.  
  489.               h11n->Fill(detStrip);
  490.               hcTrk->Fill(polarity*clustersCharge[j]);
  491.               if(clustersSize[j]==1) hcTrkCorr->Fill(polarity*clustersCharge[j]);
  492.               if(clustersSize[j]==2 && iPeak==1) hcTrkCorr->Fill(polarity*clustersCharge[j]*(1.0-chargeCorrSlopeOdd));
  493.               if(clustersSize[j]==2 && iPeak==0) hcTrkCorr->Fill(polarity*clustersCharge[j]*(1.0-chargeCorrSlopeEven));
  494.               hcTrkSeed->Fill(polarity*clustersSeedCharge[j]);
  495.               if(clustersSize[j]==1) hcTrkSeed1->Fill(polarity*clustersSeedCharge[j]);
  496.               if(clustersSize[j]==2) hcTrkSeed2->Fill(polarity*clustersSeedCharge[j]);
  497.               if(y_trk - yMax < 5.0) hcTrkNoTop->Fill(polarity*clustersCharge[j]);
  498.               if(x_trk < (xMin+xMax)/2.) hcTrkLowX->Fill(polarity*clustersCharge[j]);
  499.               if(x_trk >= (xMin+xMax)/2.) hcTrkHighX->Fill(polarity*clustersCharge[j]);
  500.               int ix = 4*(x_trk - xMin - 0.001) / (xMax-xMin);
  501.               int iy = 4*(y_trk - yMin - 0.001) / (yMax-yMin);
  502.               int ihist = 4*ix + iy;
  503.               hcTrkXY[ihist]->Fill(polarity*clustersCharge[j]);              
  504.               //if(y_trk>yMin && y_trk<=(yMin+0.25*yRange)) hcTrkY1->Fill(polarity*clustersCharge[j]);
  505.               //if(y_trk>(yMin+0.25*yRange) && y_trk<(yMin+0.50*yRange)) hcTrkY2->Fill(polarity*clustersCharge[j]);
  506.               //if(y_trk>(yMin+0.50*yRange) && y_trk<(yMin+0.75*yRange)) hcTrkY3->Fill(polarity*clustersCharge[j]);
  507.               //if(y_trk>(yMin+0.75*yRange) && y_trk<(yMin+1.00*yRange)) hcTrkY4->Fill(polarity*clustersCharge[j]);
  508.               //if(y_trk>yMin&&y_trk<yMid+1.0) hcTrkY2->Fill(polarity*clustersCharge[j]);
  509.  
  510.               if(clustersSize[j]==1) hcTrk1->Fill(polarity*clustersCharge[j]);
  511.               if(clustersSize[j]==2) hcTrk2->Fill(polarity*clustersCharge[j]);
  512.               if(clustersPosition[j]>170  && clustersPosition[j]<190) hcTrka->Fill(polarity*clustersCharge[j]);
  513.               h31a->Fill(x_trk,y_trk);
  514.               h32a->Fill(clustersPosition[j]);
  515.               h33a->Fill(y_trk);
  516.               h33b->Fill(x_trk);
  517.               if(polarity*clustersCharge[j] < 250) {
  518.                 h31b->Fill(x_trk,y_trk);
  519.                 h32b->Fill(clustersPosition[j]);
  520.                 h34->Fill(y_trk);
  521.               }
  522.               double chr = clustersCharge1StripRight[j]*polarity;
  523.               double chl = clustersCharge1StripLeft[j]*polarity;
  524.               double chr2 = clustersCharge2StripRight[j]*polarity;
  525.               double chl2 = clustersCharge2StripLeft[j]*polarity;
  526.               double pch = polarity*clustersSeedCharge[j];
  527.               int ic = pch/50.;
  528.               if(ic>=0 && ic<10 && clustersSize[j]<=2 ){
  529.                 if(iPeak==1) h41[ic]->Fill(chl-chr);
  530.                 if(iPeak==0) h42[ic]->Fill(chl-chr);
  531.                 if(chr2!=0 and chl2!=0){
  532.                   if(iPeak==1) h43[ic]->Fill(chl2-chr2);
  533.                   if(iPeak==0) h44[ic]->Fill(chl2-chr2);
  534.                 }
  535.                 
  536.               } 
  537.             }
  538.  
  539.             double clstrip = getCorrChannel(clustersPosition[j]);
  540.             h4->Fill(clstrip);
  541.             h4b->Fill(clstrip + channelOffset);
  542.             h4a->Fill(clustersPosition[j]);
  543.             h0->Fill(detStrip - clstrip);
  544.           }
  545.         }
  546.  
  547.         if(inFiducial && goodTrack && goodTime && foundHitNoFid) {
  548.           h12on->Fill(distToCutout);
  549.         }
  550.         
  551.         if(inFiducial && goodTrack && goodTime && awayFromCutout) {
  552.           h16c->Fill(dtime);
  553.           if(foundHit) {
  554.             h3c->Fill(x_trk,y_trk);
  555.             for(int j=0; j<min(clusterNumberPerEvent,10); j++){
  556.               hcfpa->Fill(polarity*clustersCharge[j]);
  557.               if(polarity*clustersCharge[j]>120) h1fpa->Fill(dxh[j]);
  558.             }        
  559.             h15c->Fill(dtime);
  560.             h15b->Fill(clustersTDC+0.1);
  561.             h15a->Fill(vec_trk_chi2ndf->at(k));
  562.             h12b->Fill(y_trk);
  563.             h12dn->Fill(x_trk); 
  564.             if(y_trk>yInt1[0]&&y_trk<yInt1[1]) h12en->Fill(x_trk);
  565.             if(y_trk>yInt2[0]&&y_trk<yInt2[1]) h12fn->Fill(x_trk);
  566.             if(y_trk>yInt3[0]&&y_trk<yInt3[1]) h12gn->Fill(x_trk);
  567.             h12hn->Fill(fracStrip);
  568.             if(y_trk>yInt1[0]&&y_trk<yInt1[1] && x_trk>-3.5&&x_trk<2.5) h12in->Fill(fracStrip);
  569.             if(y_trk>yInt2[0]&&y_trk<yInt2[1] && x_trk>-3.5&&x_trk<2.5) h12jn->Fill(fracStrip);
  570.             if(y_trk>yInt3[0]&&y_trk<yInt3[1] && x_trk>-3.5&&x_trk<2.5) h12kn->Fill(fracStrip);
  571.             if(y_trk>yInt1[0] && y_trk < yInt2[1] && x_trk>-3.5&&x_trk<-2.0) {
  572.               h14b->Fill(detStrip);
  573.               h13fpa->Fill(clusterNumberPerEvent);
  574.             }else if(y_trk>yInt3[0] && y_trk < yInt3[1]) {
  575.               h13fbpa->Fill(clusterNumberPerEvent);
  576.               for(int j=0; j<min(clusterNumberPerEvent,10); j++){
  577.                 h1fbpa->Fill(dxh[j]);
  578.                 hcfbpa->Fill(polarity*clustersCharge[j]);
  579.               }
  580.             }                
  581.           }else{     
  582.             h3b->Fill(x_trk,y_trk);          
  583.             h13mpa->Fill(clusterNumberPerEvent);
  584.             h16a->Fill(vec_trk_chi2ndf->at(k));
  585.             h16b->Fill(clustersTDC+0.1);
  586.             for(int j=0; j<min(clusterNumberPerEvent,10); j++){
  587.               hcmpa->Fill(polarity*clustersCharge[j]);
  588.               if(polarity*clustersCharge[j]>120) h1mpa->Fill(dxh[j]);
  589.             }            
  590.             nPrint++;
  591.             if(y_trk>yInt1[0] && y_trk < yInt2[1]) {
  592.               h14a->Fill(detStrip);
  593.               if(writeEventsWithMissinhHitsToFile) 
  594.                 myfile << jentry << " " << detStrip << " " << x_trk << " " << y_trk << endl;              
  595.               //if(nPrint < 100) cout << "Missed hit, event, #clu, nomStrip =  " << jentry << " " 
  596.               //                      << clusterNumberPerEvent << " " << nomStrip << endl;
  597.               for(int j=0; j<min(clusterNumberPerEvent,10); j++){
  598.                 if(x_trk>-3.5&&x_trk<-2.0) h1mpa1->Fill(dxh[j]);
  599.                 if(x_trk>-2.0&&x_trk<-1.0) h1mpa2->Fill(dxh[j]);
  600.                 if(x_trk>-1.0&&x_trk<0.0) h1mpa3->Fill(dxh[j]);
  601.                 if(x_trk>0.0&&x_trk<1.0)  h1mpa4->Fill(dxh[j]);
  602.                 if(x_trk>1.5&&x_trk<3.5)  h1mpa5->Fill(dxh[j]);
  603.                 if(y_trk<yInt1[1]) h1mpaL->Fill(clustersPosition[j]);                         
  604.                 if(y_trk>=yInt1[1]) h1mpaU->Fill(clustersPosition[j]);
  605.               }
  606.             }else if(y_trk>yInt3[0] && y_trk < yInt3[1]) {
  607.               h13mbpa->Fill(clusterNumberPerEvent);
  608.               for(int j=0; j<min(clusterNumberPerEvent,10); j++){              
  609.                 h1mbpa->Fill(dxh[j]);
  610.                 hcmbpa->Fill(polarity*clustersCharge[j]);
  611.               }    
  612.             }
  613.           }
  614.         } 
  615.       }
  616.    }
  617.    if(writeEventsWithMissinhHitsToFile) myfile.close();
  618.  
  619.    int i1 = h1->FindBin(-0.3);
  620.    int i2 = h1->FindBin(0.3);
  621.    cout << "Number of track - DUT hit matchs: " << h1->Integral(i1,i2) << endl;
  622.    cout << "Ry = " << Ry << endl;
  623.    
  624.    
  625.    TString r1 = Form("%6.3f < Y < %6.3f mm",yInt1[0], yInt1[1]);
  626.    TString r2 = Form("%6.3f < Y < %6.3f mm ",yInt2[0], yInt2[1]);
  627.    TString r3 = Form("%6.3f < Y < %6.3f mm",yInt3[0], yInt3[1]);   
  628.  
  629.  
  630.    TF1 *funchole = new TF1("funchole","[0]+[1]*x+[2]*x*x",-8.0,8.0);
  631.    funchole->SetParameters(holeQuadPar[0],holeQuadPar[1],holeQuadPar[2]);
  632.    funchole->SetLineColor(4);
  633.    funchole->SetLineWidth(2);
  634.    cout << "Hole Parameters: " << holeQuadPar[0] << " " << holeQuadPar[1] << " " << holeQuadPar[2] << endl;
  635.    
  636.    gStyle->SetOptStat(1000000001);
  637.    gStyle->SetOptFit(0011);
  638.    gStyle->SetStatH(0.05);
  639.    gStyle->SetStatW(0.2);
  640.  
  641.    addGraphics(h3b, 1, "X_{trk} [mm]", "Y_{trk} [mm] ");
  642.    h3b->GetXaxis()->SetRangeUser(-5,5);
  643.    h3b->GetYaxis()->SetRangeUser(-5,5);
  644.  
  645.    TCanvas *c = new TCanvas("c","Residuals",1500,1000);
  646.    c->Divide(4,3);
  647.    c->cd(1)->SetLeftMargin(0.13);
  648.    addGraphics(h2, 1, "X_{trk} [mm]", "X_{DUT} [mm]");
  649.    double dM = yMax;
  650.    if(xMax>dM) dM = xMax;
  651.    h2->GetXaxis()->SetRangeUser(-1.0*dM-1,dM+1);
  652.    h2->GetYaxis()->SetRangeUser(-1.0*dM-1,dM+1);
  653.    h2->SetMinimum(1);
  654.    h2->Draw();
  655.  
  656.    c->cd(2)->SetLeftMargin(0.13);;
  657.    addGraphics(h1, 1, "#DeltaX [mm]", "");
  658.    addGraphics(h1z, 2, "#DeltaX [mm]", "");
  659.    addGraphics(h1a, 3, "#DeltaX [mm]", "");
  660.    addGraphics(h1b, 6, "#DeltaX [mm]", "");
  661.    h1->GetXaxis()->SetRangeUser(-0.3,0.3);
  662.    if(m_board.Contains("A1")) h1->GetXaxis()->SetRangeUser(-1.5,1.5);
  663.    h1->SetMaximum(1.2*h1->GetMaximum());
  664.    h1->Draw();   
  665.    //h1z->Draw("same");
  666.    h1a->Draw("same");
  667.    h1b->Draw("same");
  668.    TLine *l1 = new TLine(-stripPitch/2.,0,-stripPitch/2.,0.5*h1->GetMaximum());
  669.    TLine *l2 = new TLine(stripPitch/2,0,stripPitch/2,0.5*h1->GetMaximum());
  670.    TLine *l1a = new TLine(-dxWin,0,-dxWin,0.5*h1->GetMaximum());
  671.    TLine *l2a = new TLine(dxWin,0,dxWin,0.5*h1->GetMaximum());
  672.    l1->SetLineColor(kRed); l2->SetLineColor(kRed);
  673.    l1a->SetLineColor(kBlue); l2a->SetLineColor(kBlue); l1a->SetLineWidth(2); l2a->SetLineWidth(2);
  674.    l1->Draw(); l2->Draw();
  675.    //l1a->Draw(); l2a->Draw();
  676.  
  677.    TLegend* legend3 = new TLegend(0.15,0.70,0.40,0.89);
  678.    legend3->SetFillStyle(0);
  679.    legend3->SetBorderSize(0);
  680.    legend3->SetFillStyle(0);
  681.    legend3->SetTextSize(0.045);
  682.    
  683.    legend3->AddEntry(h1,"All","L"); 
  684.    //legend3->AddEntry(h1z,"ADC < 250","L"); 
  685.    legend3->AddEntry(h1a,"1-strip","L"); 
  686.    legend3->AddEntry(h1b,"2-strip","L"); 
  687.    legend3->Draw();
  688.    
  689.    c->cd(3)->SetLeftMargin(0.13);;
  690.    addGraphics(h8, 1, "#theta_{X}^{trk} [mrad]", "#DeltaX [mm]");
  691.    h8->GetXaxis()->SetRangeUser(txMin-0.1,txMax+0.1);   
  692.    h8->GetYaxis()->SetRangeUser(-0.1,0.1);   
  693.    h8->Draw();
  694.    c->cd(4)->SetLeftMargin(0.13);;
  695.    addGraphics(h9, 1, "Y^{trk} [mm]", "#DeltaX [mm]");
  696.    h9->GetXaxis()->SetRangeUser(yMin-1,yMax+1);   
  697.    h9->GetYaxis()->SetRangeUser(-0.1,0.1);   
  698.    h9->Draw();
  699.    c->cd(5)->SetLeftMargin(0.13);;
  700.    addGraphics(h9a, 1, "X^{trk} [mm]", "#DeltaX [mm]");
  701.    h9a->GetXaxis()->SetRangeUser(-4,4);   
  702.    h9a->GetYaxis()->SetRangeUser(-0.1,0.1);   
  703.    h9a->Draw();
  704.  
  705.    c->cd(6)->SetLeftMargin(0.13);;
  706.    addGraphics(h3, 1, "X_{trk} [mm]", "Y_{trk} [mm]");
  707.    h3->GetXaxis()->SetRangeUser(xMin-1,xMax+1);   
  708.    h3->GetYaxis()->SetRangeUser(yMin-1,yMax+1);   
  709.    h3->Draw();
  710.    if(holeQuadPar[0]!=0) funchole->Draw("same");
  711.    if(fabs(xLeftHole)<900 && fabs(xRightHole)<900){
  712.      TBox *b1 = new TBox(xMin,yMin,xLeftHole,yMax);
  713.      TBox *b2 = new TBox(xRightHole,yMin,xMax,yMax);
  714.      b1->SetLineColor(kBlue); b1->SetLineWidth(2);b1->SetFillStyle(0);   
  715.      b1->Draw();
  716.      b2->SetLineColor(kBlue); b2->SetLineWidth(2);b2->SetFillStyle(0);   
  717.      b2->Draw();
  718.    }else{   
  719.      TBox *b = new TBox(xMin,yMin,xMax,yMax);
  720.      b->SetLineColor(kBlue); b->SetLineWidth(2);b->SetFillStyle(0);   
  721.      b->Draw();
  722.    }
  723.    
  724.  
  725.    c->cd(7)->SetLeftMargin(0.13);;
  726.    addGraphics(h5, 1, "#theta_{X}^{trk} [mrad]", "");
  727.    h5->GetXaxis()->SetRangeUser(txMin-0.1,txMax+0.1);   
  728.    h5->Draw();
  729.    h5c->SetLineColor(2); h5c->Draw("same");
  730.    TLine *l1b = new TLine(txMin,0,txMin,h5->GetMaximum());
  731.    TLine *l2b = new TLine(txMax,0,txMax,h5->GetMaximum());
  732.    l1b->SetLineColor(kBlue); l2b->SetLineColor(kBlue); l1b->SetLineWidth(2); l2b->SetLineWidth(2);
  733.    l1b->Draw(); l2b->Draw();
  734.  
  735.    c->cd(8)->SetLeftMargin(0.13);
  736.    addGraphics(h6, 1, "#theta_{Y}^{trk} [mrad]", "");
  737.    h6->GetXaxis()->SetRangeUser(tyMin-0.1,tyMax+0.1);   
  738.    h6->Draw();
  739.    h6c->SetLineColor(2); h6c->Draw("same");
  740.    TLine *l1c = new TLine(tyMin,0,tyMin,h6->GetMaximum());
  741.    TLine *l2c = new TLine(tyMax,0,tyMax,h6->GetMaximum());
  742.    l1c->SetLineColor(kBlue); l2c->SetLineColor(kBlue); l1c->SetLineWidth(2); l2c->SetLineWidth(2);   
  743.    l1c->Draw(); l2c->Draw();
  744.    
  745.    c->cd(9)->SetLeftMargin(0.13);;
  746.    addGraphics(h4, 1, "Strip # with cluster", "");
  747.    addGraphics(h4a, 2, "Strip # with cluster", "");
  748.    h4->GetXaxis()->SetRangeUser(max(iLo-60.0,1.0),min(iHi+60.0,512.0));
  749.    h4a->GetXaxis()->SetRangeUser(max(iLo-60.0,1.0),min(iHi+60.0,512.0));
  750.    double vmax = h4a->GetMaximum();
  751.    //for(int ii=iLo; ii<=iHi;ii++){
  752.    //  if(h4a->GetBinContent(ii)>vmax) vmax = h4a->GetBinContent(ii);
  753.    //}
  754.    h4->SetMaximum(1.25*vmax);
  755.    TLine *l1d = new TLine(iLo,0,iLo,0.75*h4a->GetMaximum());
  756.    TLine *l2d = new TLine(iHi,0,iHi,0.75*h4a->GetMaximum());
  757.    l1d->SetLineColor(kBlue); l2d->SetLineColor(kBlue); l1d->SetLineWidth(2); l2d->SetLineWidth(2);   
  758.    h4->Draw();
  759.    h4a->Draw("same");
  760.    l1d->Draw(); l2d->Draw();
  761.    
  762.    TLegend* legend4 = new TLegend(0.15,0.70,0.94,0.89);
  763.    legend4->SetFillStyle(0);
  764.    legend4->SetBorderSize(0);
  765.    legend4->SetFillStyle(0);
  766.    legend4->SetTextSize(0.045);
  767.    
  768.    legend4->AddEntry(h4,"Sensor strip ch #","L"); 
  769.    legend4->AddEntry(h4a,"Electr. ch #","L");
  770.  
  771.    legend4->Draw();
  772.  
  773.  
  774.    c->cd(10)->SetLeftMargin(0.13);;
  775.    addGraphics(hcAll, 1, "Cluster charge [ADC]", "");
  776.    addGraphics(hcTrk1, 3, "Cluster charge [ADC]", "");
  777.    addGraphics(hcTrk2, 6, "Cluster charge [ADC]", "");
  778.    hcAll->SetMaximum(1.25*hcAll->GetMaximum());
  779.    hcAll->Draw();
  780.    hcTrk->SetLineColor(kBlue); hcTrk->SetLineWidth(2);
  781.    hcTrk->Draw("same"); 
  782.    hcTrk1->Draw("same"); 
  783.    hcTrk2->Draw("same"); 
  784.   
  785.    TLegend* legend2 = new TLegend(0.15,0.70,0.94,0.89);
  786.    legend2->SetFillStyle(0);
  787.    legend2->SetBorderSize(0);
  788.    legend2->SetFillStyle(0);
  789.    legend2->SetTextSize(0.045);
  790.    
  791.    legend2->AddEntry(hcAll,"All clusters, trk in Fid","L"); 
  792.    legend2->AddEntry(hcTrk,"Clusters, |#DeltaX|<200 #mum, trk in Fid.","L");
  793.    legend2->AddEntry(hcTrk1,"1-strip Clusters","L");
  794.    legend2->AddEntry(hcTrk2,"2-strip Clusters","L");
  795.  
  796.    legend2->Draw();
  797.  
  798.    
  799.  
  800.    c->cd(11)->SetLeftMargin(0.13);;
  801.    addGraphics(h10a, 1, "Strip #", "<ADC>");
  802.    addGraphics(h10b, 2, "Strip #", "<ADC>");
  803.    addGraphics(h10c, 4, "Strip #", "<ADC>");
  804.    int ixl = iLo;
  805.    int ixh = iHi;
  806.    int yh = hcTrk->GetMean()-50;
  807.    h10a->GetXaxis()->SetRangeUser(ixl-2,ixh+2);   
  808.    h10a->GetYaxis()->SetRangeUser(yh-200,yh+300);
  809.    /*
  810.    h10a->Draw();
  811.    h10b->SetLineColor(2);h10b->SetMarkerColor(2);
  812.    h10b->Draw("same");
  813.    h10c->SetLineColor(4);h10c->SetMarkerColor(4);
  814.    h10c->Draw("same");
  815.  
  816.    TLegend* legend1 = new TLegend(0.25,0.75,0.94,0.89);
  817.    legend1->SetFillStyle(0);
  818.    legend1->SetBorderSize(0);
  819.    legend1->SetFillStyle(0);
  820.    legend1->SetTextSize(0.045);
  821.    
  822.    legend1->AddEntry(h10a,"Top half in Y","LEP"); 
  823.    legend1->AddEntry(h10b,"Bottom half in Y","LEP"); 
  824.    legend1->AddEntry(h10c,"Top 2 mm in Y","LEP"); 
  825.    legend1->Draw();
  826.    */
  827.    addGraphics(h35, 1, "#DUT clusters", "Entries");   
  828.    h35->GetYaxis()->SetTitleOffset(1.2);
  829.    h35->Draw();
  830.  
  831.    c->cd(12)->SetLeftMargin(0.13);;
  832.    addGraphics(h2p, 1, "TDC time / 2.5 ns", "<ADC>");
  833.    h2p->Draw();
  834.    TLine *l1e = new TLine(tdcLo,0,tdcLo,h2p->GetMaximum());
  835.    TLine *l2e = new TLine(tdcHi,0,tdcHi,h2p->GetMaximum());
  836.    l1e->SetLineColor(kBlue); l2e->SetLineColor(kBlue); l1e->SetLineWidth(2); l2e->SetLineWidth(2);   
  837.    l1e->Draw(); l2e->Draw();
  838.  
  839.    c->Print("Plots/plot_"+m_board+"_s"+m_sector+"_vb"+m_bias+".png");
  840.  
  841.    //return;
  842.    
  843.  
  844.  
  845.    TCanvas *c4 = new TCanvas("c4","Plot 4",800,600);
  846.    c4->Divide(2,2);
  847.    c4->cd(1)->SetLeftMargin(0.13);
  848.    addGraphics(h12en, 1, "X_{trk} [mm]", "");
  849.    addGraphics(h12fn, 1, "X_{trk} [mm]", "");
  850.  
  851.    TH1F *hepa1 = (TH1F*)h12en->Clone("hepa1");
  852.    TH1F *hepa2 = (TH1F*)h12fn->Clone("hepa2");
  853.    TH1F *henpa = (TH1F*)h12gn->Clone("henpa");
  854.    hepa1->Divide(h12en,h12ed,1.0,1.0,"B");
  855.    hepa2->Divide(h12fn,h12fd,1.0,1.0,"B");
  856.    henpa->Divide(h12gn,h12gd,1.0,1.0,"B");
  857.    addGraphics(hepa1, 1, "X_{trk} [mm]", "#Good DUT hit / # Track ");
  858.    addGraphics(hepa2, 1, "X_{trk} [mm]", "#Good DUT hit / # Track ");
  859.    h3b->GetYaxis()->SetRangeUser(-5,5);
  860.    h3b->Draw();
  861.  
  862.    c4->cd(2)->SetLeftMargin(0.13);
  863.    hepa1->GetXaxis()->SetRangeUser(xMin-0.5,xMax+0.5);
  864.    hepa1->SetMinimum(0.0);
  865.    hepa1->SetMaximum(1.2);
  866.    hepa1->SetLineWidth(1);
  867.    hepa1->SetTitle(r1);
  868.    hepa1->Draw("hist");
  869.    c4->cd(3)->SetLeftMargin(0.13);
  870.    hepa2->GetXaxis()->SetRangeUser(xMin-0.5,xMax+0.5);
  871.    hepa2->SetMinimum(0.0);
  872.    hepa2->SetMaximum(1.2);
  873.    hepa2->SetLineWidth(1);
  874.    hepa2->SetTitle(r2);
  875.    hepa2->Draw("hist");
  876.  
  877.    c4->cd(4)->SetLeftMargin(0.13);
  878.    addGraphics(henpa, 1, "X_{trk} [mm]", "#Good DUT hit / # Track ");
  879.    henpa->GetXaxis()->SetRangeUser(xMin-0.5,xMax+0.5);
  880.    henpa->SetMinimum(0.0);
  881.    henpa->SetMaximum(1.2);
  882.    henpa->SetLineWidth(1);
  883.    henpa->SetTitle(r3);
  884.    henpa->Draw("hist");   
  885.  
  886.  
  887.    TCanvas *c5 = new TCanvas("c5","Plot 5",800,600);
  888.    c5->Divide(2,2);
  889.    c5->cd(1)->SetLeftMargin(0.13);
  890.    addGraphics(h12hn, 1, "Rel. Strip Pos", "");
  891.    addGraphics(h12in, 1, "Rel. Strip Pos", "");
  892.    addGraphics(h12jn, 1, "Rel. Strip Pos", "");
  893.    addGraphics(h12kn, 1, "Rel. Strip Pos", "");
  894.  
  895.    TH1F *hepas1 = (TH1F*)h12hn->Clone("hepas1");
  896.    TH1F *hepas2 = (TH1F*)h12in->Clone("hepas2");
  897.    TH1F *hepas3 = (TH1F*)h12jn->Clone("hepas3");
  898.    TH1F *hepas4 = (TH1F*)h12kn->Clone("hepas4");
  899.    hepas1->Divide(h12hn,h12hd,1.0,1.0,"B");
  900.    hepas2->Divide(h12in,h12id,1.0,1.0,"B");
  901.    hepas3->Divide(h12jn,h12jd,1.0,1.0,"B");
  902.    hepas4->Divide(h12kn,h12kd,1.0,1.0,"B");
  903.    addGraphics(hepas1, 1, "Rel. Strip Pos.", "#Good DUT hit / # Track ");
  904.    addGraphics(hepas2, 1, "Rel. Strip Pos.", "#Good DUT hit / # Track ");
  905.    addGraphics(hepas3, 1, "Rel. Strip Pos.", "#Good DUT hit / # Track ");
  906.    addGraphics(hepas4, 1, "Rel. Strip Pos.", "#Good DUT hit / # Track ");
  907.    hepas1->SetTitle("Full Beam Spot");
  908.    hepas1->SetMinimum(0.8);
  909.    hepas1->SetMaximum(1.2);
  910.    hepas1->SetLineWidth(1);
  911.    hepas1->Draw("hist");
  912.    c5->cd(2)->SetLeftMargin(0.13);
  913.    hepas2->SetTitle(r1);
  914.    hepas2->SetMinimum(0.8);
  915.    hepas2->SetMaximum(1.2);
  916.    hepas2->SetLineWidth(1);
  917.    hepas2->Draw("hist");
  918.    c5->cd(3)->SetLeftMargin(0.13);
  919.    hepas3->SetTitle(r2);
  920.    hepas3->SetMinimum(0.8);
  921.    hepas3->SetMaximum(1.2);
  922.    hepas3->SetLineWidth(1);
  923.    hepas3->Draw("hist");
  924.    c5->cd(4)->SetLeftMargin(0.13);
  925.    hepas4->SetTitle(r3);
  926.    hepas4->SetMinimum(0.8);
  927.    hepas4->SetMaximum(1.2);
  928.    hepas4->SetLineWidth(1);
  929.    hepas4->Draw("hist");
  930.  
  931.    c5->Print("Plots/plot5_"+m_board+"_s"+m_sector+"_vb"+m_bias+".png");
  932.    //return;
  933.    
  934.  
  935.    TCanvas *c1 = new TCanvas("c1","Plot 1",1600,800);
  936.    c1->Divide(3,2);
  937.  
  938.    c1->cd(1)->SetLeftMargin(0.13);
  939.    addGraphics(h10d, 2, "Y_{trk} [mm]", "<ADC>");
  940.    h10d->GetXaxis()->SetRangeUser(yMin-0.5,yMax+0.5);
  941.    h10d->Draw();
  942.  
  943.    //c1->cd(2)->SetLeftMargin(0.13);
  944.    //h3b->Draw();   
  945.  
  946.    c1->cd(2)->SetLeftMargin(0.13);
  947.    addGraphics(h12b, 2, "Y_{trk} [mm]", "");
  948.    addGraphics(h12a, 1, "Y_{trk} [mm]", "");
  949.    h12a->SetMaximum(1.2*h12a->GetMaximum());
  950.    h12a->GetXaxis()->SetRangeUser(yMin-0.5,yMax+0.5);
  951.    h12a->Draw("hist");
  952.    h12b->Draw("hist,same");
  953.  
  954.    TLegend* legend0 = new TLegend(0.15,0.75,0.94,0.89);
  955.    legend0->SetFillStyle(0);
  956.    legend0->SetBorderSize(0);
  957.    legend0->SetFillStyle(0);
  958.    legend0->SetTextSize(0.045);   
  959.    legend0->AddEntry(h12a,"Track Y pos","LEP"); 
  960.    legend0->AddEntry(h12b,"Track Y pos with good cluster","LEP");
  961.    legend0->Draw();
  962.  
  963.    c1->cd(3)->SetLeftMargin(0.13);
  964.    TH1F *he = (TH1F*)h12b->Clone("he");
  965.    he->Divide(h12b,h12a,1.0,1.0,"B");
  966.    float bw = 1000*h12b->GetBinWidth(1);
  967.    TString yt = Form("(#Good DUT hit / # Track) / %3.0f #mum",bw);
  968.    addGraphics(he, 1, "Y_{trk} [mm]", yt);
  969.    he->SetTitle("DUT Efficiency vs Y_{trk}");
  970.    he->GetXaxis()->SetRangeUser(yMin-0.5,yMax+0.5);
  971.    he->SetMinimum(0.8);
  972.    he->SetMaximum(1.1);
  973.    he->Draw("e");
  974.  
  975.    c1->cd(4)->SetLeftMargin(0.13);
  976.    addGraphics(h10e, 2, "X_{trk} [mm]", "<ADC>");
  977.    h10e->GetXaxis()->SetRangeUser(xMin-0.5,xMax+0.5);
  978.    h10e->Draw();
  979.  
  980.    c1->cd(5)->SetLeftMargin(0.13);
  981.    addGraphics(h12dd, 2, "X_{trk} [mm]", "");
  982.    addGraphics(h12dn, 1, "X_{trk} [mm]", "");
  983.    h12dd->SetMaximum(1.2*h12dd->GetMaximum());
  984.    h12dd->GetXaxis()->SetRangeUser(xMin-0.5,xMax+0.5);
  985.    h12dd->Draw("hist");
  986.    h12dn->Draw("hist,same");
  987.  
  988.    c1->cd(6)->SetLeftMargin(0.13);
  989.    TH1F *he2 = (TH1F*)h12dn->Clone("he");
  990.    he2->Divide(h12dn,h12dd,1.0,1.0,"B");
  991.    bw = 1000*h12dd->GetBinWidth(1);
  992.    yt = Form("(#Good DUT hit / # Track) / %3.0f #mum",bw);
  993.    addGraphics(he2, 1, "X_{trk} [mm]", yt);
  994.    he2->SetTitle("DUT Efficiency vs X_{trk}");
  995.    he2->GetXaxis()->SetRangeUser(xMin-0.5,xMax+0.5);
  996.    he2->SetMinimum(0.8);
  997.    he2->SetMaximum(1.1);
  998.    he2->Draw("e");
  999.    c1->Print("Plots/plot1_"+m_board+"_s"+m_sector+"_vb"+m_bias+".png");
  1000.  
  1001.    TCanvas *c8 = new TCanvas("c8","Plot 8",800,400);
  1002.    c8->Divide(2,1);
  1003.    TH1F *he4 = (TH1F*)h15c->Clone("Eff_trigger-DUT_time");
  1004.    he4->Divide(h15c,h16c,1.0,1.0,"B");
  1005.    yt = "#Good DUT hit / # Track";
  1006.    addGraphics(h16c, 1, "DUT time - Track Time (ns)", "Entries");
  1007.    addGraphics(he4, 1, "DUT time - Track Time (ns)", "Efficiency");
  1008.    c8->cd(1)->SetLeftMargin(0.15);;
  1009.    c8->cd(1)->SetRightMargin(0.05);;
  1010.    h16c->Draw();
  1011.    c8->cd(2)->SetLeftMargin(0.13);;
  1012.    c8->cd(2)->SetRightMargin(0.05);;
  1013.    he4->Draw();
  1014.    c8->Print("Plots/plot8_"+m_board+"_s"+m_sector+"_vb"+m_bias+".png");   
  1015.  
  1016.    if(holeSector){
  1017.      TCanvas *c7 = new TCanvas("c7","Plot 7",800,800);
  1018.      c7->SetLeftMargin(0.13);
  1019.      TH1F *he3 = (TH1F*)h12on->Clone("Erf()_fit");
  1020.      he3->Divide(h12on,h12od,1.0,1.0,"B");
  1021.      bw = 1000*h12od->GetBinWidth(1);
  1022.      yt = Form("(#Good DUT hit / # Track) / %3.0f #mum",bw);
  1023.      addGraphics(he3, 1, "Dist to cutout [mm]", yt);
  1024.      he3->SetStats(kTRUE);
  1025.      he3->SetTitle("DUT Efficiency vs Dist to Cutout");
  1026.      he3->GetXaxis()->SetRangeUser(-0.2,2.0);
  1027.      he3->GetYaxis()->SetRangeUser(0.0,1.3);
  1028.      he3->Draw("e");
  1029.      TLatex *myLatex = new TLatex();
  1030.      myLatex->SetTextFont(42); myLatex->SetTextColor(1); 
  1031.      myLatex->SetTextAlign(12); myLatex->SetNDC(kTRUE); myLatex->SetTextSize(0.047);
  1032.      TString text = "Sector "+m_sector;
  1033.      c7->cd();
  1034.      myLatex->DrawLatex(0.18,0.8,text);
  1035.      
  1036.      TF1 *f4a = new TF1("f4a","0.5*[0]*(1+TMath::Erf((x-[1])/[2]))",-0.1,2.0);
  1037.      f4a->SetParameters(0.5,0.0,0.04);
  1038.      f4a->SetParNames("Const","Mean","Sigma");
  1039.      he3->Fit("f4a","R");
  1040.      TPad *inset = new TPad("inset","y vx x",0.44,0.1,0.89,0.55);
  1041.      inset->SetLeftMargin(0.13);
  1042.      inset->SetRightMargin(0.05);
  1043.      h3->SetTitle("");
  1044.      inset->Draw();
  1045.      inset->cd();
  1046.      h3->Draw();
  1047.      h3b->SetMarkerSize(0.2);h3b->SetMarkerStyle(20);h3b->SetMarkerColor(kRed); 
  1048.      h3b->Draw("same");
  1049.      if(holeQuadPar[0]!=0) funchole->Draw("same");   
  1050.      
  1051.      c7->Print("Plots/plot7_"+m_board+"_s"+m_sector+"_vb"+m_bias+".png");   
  1052.    }
  1053.    
  1054.    /*
  1055.    TCanvas *c2 = new TCanvas("c2","Plot 2",1600,600);
  1056.    c2->Divide(5,2);
  1057.    addGraphics(h1mpa, 2, "X_{clu}-X_{trk} [mm]", ""); h1mpa->SetTitle("2.8<Y_{trk}<4.1 mm, All X, No Match Clu");
  1058.    addGraphics(h1mbpa, 4, "X_{clu}-X_{trk} [mm]", ""); h1mbpa->SetTitle("1.0<Y_{trk}<2.3 mm, All X, No Match Clu");
  1059.    addGraphics(h1fpa, 1, "X_{clu}-X_{trk} [mm]", ""); h1fpa->SetTitle("2.8<Y_{trk}<4.1 mm, All X, Match Clu");
  1060.    addGraphics(h1fbpa, 1, "X_{clu}-X_{trk} [mm]", ""); h1fbpa->SetTitle("1.0<Y_{trk}<2.3 mm, All X, Match Clu");
  1061.  
  1062.    addGraphics(h1mpa1, 2, "X_{clu}-X_{trk} [mm]", ""); h1mpa1->SetTitle("2.8<Y_{trk}<4.1 mm, -3.5<X_{trk}<-2.0 mm, No Match Clu");
  1063.    addGraphics(h1mpa2, 2, "X_{clu}-X_{trk} [mm]", ""); h1mpa2->SetTitle("2.8<Y_{trk}<4.1 mm, -2.0<X_{trk}<-1.0 mm, No Match Clu");
  1064.    addGraphics(h1mpa3, 2, "X_{clu}-X_{trk} [mm]", ""); h1mpa3->SetTitle("2.8<Y_{trk}<4.1 mm, -1.0<X_{trk}<0.0 mm, No Match Clu");
  1065.    addGraphics(h1mpa4, 2, "X_{clu}-X_{trk} [mm]", ""); h1mpa4->SetTitle("2.8<Y_{trk}<4.1 mm,  0.0<X_{trk}<1.0 mm, No Match Clu");
  1066.    addGraphics(h1mpa5, 2, "X_{clu}-X_{trk} [mm]", ""); h1mpa5->SetTitle("2.8<Y_{trk}<4.1 mm,  1.5<X_{trk}<3.5 mm, No Match Clu");
  1067.    h1mpa->SetLineWidth(1); 
  1068.    h1mpa1->SetLineWidth(1); 
  1069.    h1mpa2->SetLineWidth(1); 
  1070.    h1mpa3->SetLineWidth(1); 
  1071.    h1mpa4->SetLineWidth(1); 
  1072.    h1mpa5->SetLineWidth(1); 
  1073.    h1fpa->SetLineWidth(1); 
  1074.    h1fbpa->SetLineWidth(1); 
  1075.  
  1076.    
  1077.    h1fpa->GetXaxis()->SetRangeUser(-3,9); 
  1078.    h1fbpa->GetXaxis()->SetRangeUser(-3,9); 
  1079.    h1mbpa->GetXaxis()->SetRangeUser(-3,9); 
  1080.    h1mpa->GetXaxis()->SetRangeUser(-3,9); 
  1081.    h1mpa1->GetXaxis()->SetRangeUser(-3,9); 
  1082.    h1mpa2->GetXaxis()->SetRangeUser(-3,9); 
  1083.    h1mpa3->GetXaxis()->SetRangeUser(-3,9); 
  1084.    h1mpa4->GetXaxis()->SetRangeUser(-3,9); 
  1085.    h1mpa5->GetXaxis()->SetRangeUser(-3,9); 
  1086.  
  1087.    c2->cd(1);
  1088.    h3b->Draw();   
  1089.    double xl[5] = {-3.5, -2.0, -1.0, 0.0, 1.5};
  1090.    double xh[5] = {-2.0, -1.0,  0.0, 1.0, 3.5};
  1091.    
  1092.    
  1093.    for(int i=0;i<5;i++){
  1094.      TBox *b = new TBox(xl[i],yInt1[0],xh[i],yInt2[1]);
  1095.      b->SetLineColor(2); b->SetFillStyle(0);b->SetLineWidth(2);
  1096.      b->Draw();
  1097.    }
  1098.    TBox *b = new TBox(xl[0],yInt3[0],xh[4],yInt3[1]);
  1099.    b->SetLineColor(4);b->SetFillStyle(0);b->SetLineWidth(2);
  1100.    b->Draw();
  1101.    
  1102.  
  1103.  
  1104.    c2->cd(2);
  1105.    h1mpa->Draw();   
  1106.    c2->cd(3);
  1107.    h1fpa->GetYaxis()->SetRangeUser(0,30);
  1108.    h1fpa->Draw();   
  1109.    c2->cd(4);
  1110.    h1mbpa->Draw();   
  1111.    c2->cd(5);
  1112.    h1fbpa->GetYaxis()->SetRangeUser(0,30);
  1113.    h1fbpa->Draw();   
  1114.    c2->cd(6);   
  1115.    h1mpa1->Draw();   
  1116.    c2->cd(7);
  1117.    h1mpa2->Draw();   
  1118.    c2->cd(8);
  1119.    h1mpa3->Draw();   
  1120.    c2->cd(9);
  1121.    h1mpa4->Draw();   
  1122.    c2->cd(10);
  1123.    h1mpa5->Draw();
  1124.    c2->Print("Plots/plot2_"+m_board+"_s"+m_sector+"_vb"+m_bias+".png");
  1125.    */
  1126.  
  1127.    /*
  1128.    TCanvas *c3 = new TCanvas("c3","Plot 3",800,600);
  1129.    c3->Divide(1,2);
  1130.    c3->cd(1);
  1131.    h1mpaU->GetXaxis()->SetRangeUser(iLo-10,iHi+10);
  1132.    h1mpaU->Draw();
  1133.    c3->cd(2);
  1134.    h1mpaL->GetXaxis()->SetRangeUser(iLo-10,iHi+10);
  1135.    h1mpaL->Draw();   
  1136.    c3->Print("Plots/plot3_"+m_board+"_s"+m_sector+"_vb"+m_bias+".png");
  1137.    */
  1138.  
  1139.    addGraphics(h12m,1,"Interstrip Pos","<ADC>");
  1140.    addGraphics(h12n,1,"Interstrip Pos","<Cluster Size>");
  1141.    TCanvas *c6 = new TCanvas("c6","Plot 6",1200,400);
  1142.    c6->Divide(3,1);
  1143.    c6->cd(1);
  1144.    hepas1->SetMaximum(1.2);
  1145.    hepas1->SetMinimum(0.5);
  1146.    hepas1->Draw();
  1147.    c6->cd(2);
  1148.    h12m->SetMinimum(0);
  1149.    h12m->SetMaximum(1.2*h12m->GetMaximum());
  1150.    h12m->Draw();
  1151.    c6->cd(3);
  1152.    h12n->SetMinimum(0.5);
  1153.    h12n->SetMaximum(2.0);
  1154.    h12n->Draw();
  1155.    
  1156.    c6->Print("Plots/plot6_"+m_board+"_s"+m_sector+"_vb"+m_bias+".png");
  1157.  
  1158.  
  1159.    fout->Write();
  1160.    
  1161.  
  1162. }