--- quicklook/tracker/flight/macros/FTrkCalibQLook_EXPERT.cxx	2006/05/10 10:06:54	1.1
+++ quicklook/tracker/flight/macros/FTrkCalibQLook_EXPERT.cxx	2006/06/29 11:11:59	1.6
@@ -1,8 +1,8 @@
 /**
- * FTrkCalibQLookExpert.cpp
+ * FTrkCalibQLookExpert.cxx
  *
  * autor: D.Fedele
- * version 2.0
+ * version v1r04
  * Parameters:
  * 	file - the data file to analyze
  *      step - select =1 in order to analyze one event at time
@@ -18,6 +18,7 @@
 #include <TPaveText.h>
 #include <TLatex.h>
 #include <TCanvas.h>
+#include <TFile.h>
 #include <TTree.h>
 #include <TGraph.h>
 #include <TStyle.h>
@@ -29,19 +30,7 @@
 #include <CalibTrk2Event.h>
 //
 
-
-void stringcopy(TString& s1, const TString& s2, Int_t from=0, Int_t to=0){    
-    if ( to == 0 ){
-	Int_t t2length = s2.Length();
-	s1 = "";
-	to = t2length;
-    };
-    for (Int_t i = from; i<to; i++){
-	s1.Append(s2[i],1);
-    };
-};
-
- typedef struct caltrk_def{
+typedef struct caltrk_def{
   Int_t good0[2];
   Int_t daqmode[12];
   Int_t dspnum[12];
@@ -70,12 +59,11 @@
   // 
   //   obtain information about the data file and select the output dir
   const string filepath=file.Data(); 
-  Int_t dwpos = filepath.rfind("DW_");
-  Int_t dwpos1 = filepath.find(".root");
-  TString fpath=(filepath.c_str());
-  TString base,ffile;
-  stringcopy(base,fpath,0,dwpos);
-  stringcopy(ffile,fpath,dwpos,dwpos1);
+  Int_t dwpos = file.Last('/');
+  Int_t dwpos1 = file.Last('.');
+  TString base,ffile ;
+  ffile=file(dwpos+1,dwpos1-(dwpos+1));
+  if(dwpos>0) base=file(0,dwpos);
 
   TString out;
   if(outdir.Length()==0){
@@ -83,6 +71,7 @@
   }else{
     out = outdir;
   }
+  if(out.Last('/')+1<out.Length()) out+="/";
 
   //
   // inizialise the variables and open the file
@@ -306,13 +295,13 @@
     TH1F *histoasig[12];           //histos of sigma 
     TH1F *histoaped[12];          //histos of pedestals
     
-    Double_t posy = 0.95;    // up y-coord - top pads
-    Double_t hpad = 0.15;   // pad height
-    Double_t posx1=0;          // left  x-coord - pad column
-    Double_t posx2=0;          // right x-coord - pad olumn
-    Double_t posx0=0;          //       x-coord - column division
-    Double_t wrel = 0.6;     // relative x size of first sub-column
-    Double_t marg = 0.004;   // margin among pads
+    Float_t posy = 0.95;    // up y-coord - top pads
+    Float_t hpad = 0.15;   // pad height
+    Float_t posx1=0;          // left  x-coord - pad column
+    Float_t posx2=0;          // right x-coord - pad olumn
+    Float_t posx0=0;          //       x-coord - column division
+    Float_t wrel = 0.6;     // relative x size of first sub-column
+    Float_t marg = 0.004;   // margin among pads
     stringstream title;
     stringstream hid;
 
@@ -335,20 +324,20 @@
       /* -----------> HISTOGRAMS          */
       /* calibration parameters */
       title<<"DSP "<<n+1;
-      hid<<"h"<<n;
+      hid<<"h"<<n<<"i"<<i;
       histosig[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3072,0.5,3072.5);
       hid.str("");
-      hid<<"hh"<<n;
+      hid<<"hh"<<n<<"i"<<i;
       histoped[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3072,0.5,3072.5);
       hid.str("");
-      hid<<"hhh"<<n;
+      hid<<"hhh"<<n<<"i"<<i;
       title.str("");
       hid.str("");
       /* AVERAGE calibration parameters */    
-      hid<<"ah"<<n;
+      hid<<"ah"<<n<<"i"<<i;
       histoasig[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3,0.5,3072.5);
       hid.str("");
-      hid<<"ahh"<<n;
+      hid<<"ahh"<<n<<"i"<<i;
       histoaped[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3,0.5,3072.5);
       hid.str("");
     };                     //end loop on views
@@ -453,7 +442,10 @@
       li.SetLineStyle(3);
       li.SetLineWidth(2);
       
-      
+      Float_t maxhist=0;
+      TBox b;
+      b.SetFillColor(6);
+      b.SetFillStyle(3945);
       /* plot PEDESTAL */
       c1->cd();
       trkpadtext[nn]->Draw();
@@ -466,25 +458,42 @@
       histoped[nn]->SetLineWidth(1);
       histoped[nn]->GetYaxis()->SetTitle("PED (ADC channels)");
       histoped[nn]->GetYaxis()->CenterTitle();
-      if((nn+1)%2==1)  histoped[nn]->GetYaxis()->SetRangeUser(2500,3200);
-      if((nn+1)%2==0)  histoped[nn]->GetYaxis()->SetRangeUser(1000,1700);
+      if((nn+1)%2==1)  histoped[nn]->GetYaxis()->SetRangeUser(2200,3200);
+      if((nn+1)%2==0)  histoped[nn]->GetYaxis()->SetRangeUser(700,1700);
       histoaped[nn]->SetLineColor(5);
       histoaped[nn]->SetLineWidth(1);
-      if(ctrk.good0[0]==1 && ctrk.good0[1]==1) histoped[nn]->Draw("b");
+      if(ctrk.good0[0]==1 && ctrk.good0[1]==1){
+	histoped[nn]->Draw("b");
+	if(nn==1){
+	  maxhist=histoped[nn]->GetMaximum();
+	  b.DrawBox(2816.,700.,3060.,maxhist);
+	}
+	else if(nn==6){
+	  maxhist=histoped[nn]->GetMaximum();
+	  b.DrawBox(2560.,2200.,2816.,maxhist);
+	  b.DrawBox(512.,2200.,768.,maxhist);
+	  b.DrawBox(1024.,2200.,1792.,maxhist);
+	}
+	else if(nn==11){
+	  maxhist=histoped[nn]->GetMaximum();
+	  b.DrawBox(768.,700.,1024.,maxhist);
+	}
+      }
       else histoped[nn]->Draw("axis");
       histoaped[nn]->Draw("same");
       if((nn+1)%2==1) {
-	li.DrawLine(1024.5,2500,1024.5,3200);
-	li.DrawLine(2048.5,2500,2048.5,3200);
+	li.DrawLine(1024.5,2200,1024.5,3200);
+	li.DrawLine(2048.5,2200,2048.5,3200);
       }
       if((nn+1)%2==0) {
-	li.DrawLine(1024.5,1000,1024.5,1700);
-	li.DrawLine(2048.5,1000,2048.5,1700);
+	li.DrawLine(1024.5,700,1024.5,1700);
+	li.DrawLine(2048.5,700,2048.5,1700);
       }
       
+ 
       
       /* plot SIGMA */
-      Double_t max=500.;
+      Float_t max=500.;
       c2->cd();
       trkpadtext[nn]->Draw();
       trkpad2[nn]->SetLogy();
@@ -501,12 +510,30 @@
       histosig[nn]->GetYaxis()->CenterTitle();
       histoasig[nn]->SetLineColor(5);
       histoasig[nn]->SetLineWidth(1);
-      if(ctrk.good0[0]==1 && ctrk.good0[1]==1) histosig[nn]->Draw("b"); 
+      if(ctrk.good0[0]==1 && ctrk.good0[1]==1){
+	histosig[nn]->Draw("b");
+	if(nn==1){
+	  maxhist=histosig[nn]->GetMaximum();
+	  b.DrawBox(2816.,0.,3060.,maxhist);
+	}
+	else if(nn==6){
+	  maxhist=histosig[nn]->GetMaximum();
+	  b.DrawBox(2560.,0.,2816.,maxhist);
+	  b.DrawBox(512.,0.,768.,maxhist);
+	  b.DrawBox(1024.,0.,1792.,maxhist);
+	}
+	else if(nn==11){
+	  maxhist=histosig[nn]->GetMaximum();
+	  b.DrawBox(768.,0.,1024.,maxhist);
+	}
+      }
       else histosig[nn]->Draw("axis");
       histoasig[nn]->Draw("same");
       li.DrawLine(1024.5,0,1024.5,max);
       li.DrawLine(2048.5,0,2048.5,max);
-        
+       
+
+
     };//end loop on views
     
     c1->Update();//draw pads in canvas