--- quicklook/tracker/flight/macros/FTrkCalibQLook_BASIC.cxx	2006/06/01 15:23:12	1.2
+++ quicklook/tracker/flight/macros/FTrkCalibQLook_BASIC.cxx	2006/09/19 15:15:42	1.8
@@ -2,7 +2,7 @@
  * FTrkCalibQLook_BASIC.cxx
  *
  * autor: D.Fedele
- * version 2.0
+ * version v1r10
  * Parameters:
  * 	file - the data file to analyze
  *      fromevent - first event to analyze
@@ -15,6 +15,7 @@
 #include <TLatex.h>
 #include <TCanvas.h>
 #include <TGraph.h>
+#include <TFile.h>
 #include <TTree.h>
 #include <TStyle.h>
 #include <TString.h>
@@ -25,18 +26,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];
@@ -64,14 +54,11 @@
 {
   // 
   //   obtain information about the data file and select the output dir
-  const string filepath=file.Data(); 
-  Int_t dwpos = filepath.rfind("/");
-  Int_t dwpos1 = filepath.rfind(".root");
-  TString fpath=(filepath.c_str());
-  TString base,ffile;
-  stringcopy(base,fpath,0,dwpos);
-  stringcopy(ffile,fpath,dwpos+1,dwpos1);
-  if(dwpos>0) base+="/";
+  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){
@@ -79,7 +66,8 @@
   }else{
     out = outdir;
   };
-  
+  if(out.Last('/')+1<out.Length()) out+="/";
+
   //
   // inizialise the variables and open the file
   struct caltrk_def ctrk;
@@ -175,7 +163,7 @@
     } else if (toevent > nevents) {	   
       maxevent = nevents;
     } else {	   
-      maxevent = fromevent;
+      maxevent = nevents;
     };
   };
 
@@ -284,53 +272,61 @@
     for (Int_t m = 0; m < 6; m++){
       ph1 = eh1->GetPscuHeader();
       OBT[0]= ph1->GetOrbitalTime();
-      ctrk.daqmode[trk1->DSPnumber[m]-1]=trk1->DAQmode[m];
-      ctrk.dspnum[trk1->DSPnumber[m]-1]=trk1->DSPnumber[m];
-      ctrk.calibnum[trk1->DSPnumber[m]-1]=trk1->calibnumber[m];
-      ctrk.ncalev[trk1->DSPnumber[m]-1]=trk1->ncalib_event[m];
-      ctrk.ped1[trk1->DSPnumber[m]-1]=trk1->ped_l1[m];
-      ctrk.ped2[trk1->DSPnumber[m]-1]=trk1->ped_l2[m];
-      ctrk.ped3[trk1->DSPnumber[m]-1]=trk1->ped_l3[m];
-      ctrk.sig1[trk1->DSPnumber[m]-1]=trk1->sig_l1[m];
-      ctrk.sig2[trk1->DSPnumber[m]-1]=trk1->sig_l2[m];
-      ctrk.sig3[trk1->DSPnumber[m]-1]=trk1->sig_l3[m];
-      ctrk.nbad1[trk1->DSPnumber[m]-1]=trk1->nbad_l1[m];
-      ctrk.nbad2[trk1->DSPnumber[m]-1]=trk1->nbad_l2[m];
-      ctrk.nbad3[trk1->DSPnumber[m]-1]=trk1->nbad_l3[m];
-      ctrk.calfl[trk1->DSPnumber[m]-1]=trk1->cal_flag[m];
-      ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk1->crc_cal[m][0];
-      ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk1->crc_cal[m][1];
-      ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk1->crc_cal[m][2];
-      ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk1->crc_hcal[m];
-      for (Int_t j = 0; j < 3072; j++){
-	ctrk.dspped[trk1->DSPnumber[m]-1][j]=trk1->DSPped_par[m][j];
-	ctrk.dspsig[trk1->DSPnumber[m]-1][j]=trk1->DSPsig_par[m][j];
-	ctrk.dspbad[trk1->DSPnumber[m]-1][j]=trk1->DSPbad_par[m][j];
+      if(trk1->DSPnumber[m]>0){
+	if(trk1->DSPnumber[m]<13){
+	  ctrk.daqmode[trk1->DSPnumber[m]-1]=trk1->DAQmode[m];
+	  ctrk.dspnum[trk1->DSPnumber[m]-1]=trk1->DSPnumber[m];
+	  ctrk.calibnum[trk1->DSPnumber[m]-1]=trk1->calibnumber[m];
+	  ctrk.ncalev[trk1->DSPnumber[m]-1]=trk1->ncalib_event[m];
+	  ctrk.ped1[trk1->DSPnumber[m]-1]=trk1->ped_l1[m];
+	  ctrk.ped2[trk1->DSPnumber[m]-1]=trk1->ped_l2[m];
+	  ctrk.ped3[trk1->DSPnumber[m]-1]=trk1->ped_l3[m];
+	  ctrk.sig1[trk1->DSPnumber[m]-1]=trk1->sig_l1[m];
+	  ctrk.sig2[trk1->DSPnumber[m]-1]=trk1->sig_l2[m];
+	  ctrk.sig3[trk1->DSPnumber[m]-1]=trk1->sig_l3[m];
+	  ctrk.nbad1[trk1->DSPnumber[m]-1]=trk1->nbad_l1[m];
+	  ctrk.nbad2[trk1->DSPnumber[m]-1]=trk1->nbad_l2[m];
+	  ctrk.nbad3[trk1->DSPnumber[m]-1]=trk1->nbad_l3[m];
+	  ctrk.calfl[trk1->DSPnumber[m]-1]=trk1->cal_flag[m];
+	  ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk1->crc_cal[m][0];
+	  ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk1->crc_cal[m][1];
+	  ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk1->crc_cal[m][2];
+	  ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk1->crc_hcal[m];
+	  for (Int_t j = 0; j < 3072; j++){
+	    ctrk.dspped[trk1->DSPnumber[m]-1][j]=trk1->DSPped_par[m][j];
+	    ctrk.dspsig[trk1->DSPnumber[m]-1][j]=trk1->DSPsig_par[m][j];
+	    ctrk.dspbad[trk1->DSPnumber[m]-1][j]=trk1->DSPbad_par[m][j];
+	  }
+	}
       }
       ph2 = eh2->GetPscuHeader();
       OBT[1]= ph2->GetOrbitalTime();
-      ctrk.daqmode[trk2->DSPnumber[m]-1]=trk2->DAQmode[m];
-      ctrk.dspnum[trk2->DSPnumber[m]-1]=trk2->DSPnumber[m];
-      ctrk.calibnum[trk2->DSPnumber[m]-1]=trk2->calibnumber[m];
-      ctrk.ncalev[trk2->DSPnumber[m]-1]=trk2->ncalib_event[m];
-      ctrk.ped1[trk2->DSPnumber[m]-1]=trk2->ped_l1[m];
-      ctrk.ped2[trk2->DSPnumber[m]-1]=trk2->ped_l2[m];
-      ctrk.ped3[trk2->DSPnumber[m]-1]=trk2->ped_l3[m];
-      ctrk.sig1[trk2->DSPnumber[m]-1]=trk2->sig_l1[m];
-      ctrk.sig2[trk2->DSPnumber[m]-1]=trk2->sig_l2[m];
-      ctrk.sig3[trk2->DSPnumber[m]-1]=trk2->sig_l3[m];
-      ctrk.nbad1[trk2->DSPnumber[m]-1]=trk2->nbad_l1[m];
-      ctrk.nbad2[trk2->DSPnumber[m]-1]=trk2->nbad_l2[m];
-      ctrk.nbad3[trk2->DSPnumber[m]-1]=trk2->nbad_l3[m];
-      ctrk.calfl[trk2->DSPnumber[m]-1]=trk2->cal_flag[m];
-      ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk2->crc_cal[m][0];
-      ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk2->crc_cal[m][1];
-      ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk2->crc_cal[m][2];
-      ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk2->crc_hcal[m];
-      for (Int_t j = 0; j < 3072; j++){
-	ctrk.dspped[trk2->DSPnumber[m]-1][j]=trk2->DSPped_par[m][j];
-	ctrk.dspsig[trk2->DSPnumber[m]-1][j]=trk2->DSPsig_par[m][j];
-	ctrk.dspbad[trk2->DSPnumber[m]-1][j]=trk2->DSPbad_par[m][j];
+      if(trk2->DSPnumber[m]>0){
+	if(trk2->DSPnumber[m]<13){
+	  ctrk.daqmode[trk2->DSPnumber[m]-1]=trk2->DAQmode[m];
+	  ctrk.dspnum[trk2->DSPnumber[m]-1]=trk2->DSPnumber[m];
+	  ctrk.calibnum[trk2->DSPnumber[m]-1]=trk2->calibnumber[m];
+	  ctrk.ncalev[trk2->DSPnumber[m]-1]=trk2->ncalib_event[m];
+	  ctrk.ped1[trk2->DSPnumber[m]-1]=trk2->ped_l1[m];
+	  ctrk.ped2[trk2->DSPnumber[m]-1]=trk2->ped_l2[m];
+	  ctrk.ped3[trk2->DSPnumber[m]-1]=trk2->ped_l3[m];
+	  ctrk.sig1[trk2->DSPnumber[m]-1]=trk2->sig_l1[m];
+	  ctrk.sig2[trk2->DSPnumber[m]-1]=trk2->sig_l2[m];
+	  ctrk.sig3[trk2->DSPnumber[m]-1]=trk2->sig_l3[m];
+	  ctrk.nbad1[trk2->DSPnumber[m]-1]=trk2->nbad_l1[m];
+	  ctrk.nbad2[trk2->DSPnumber[m]-1]=trk2->nbad_l2[m];
+	  ctrk.nbad3[trk2->DSPnumber[m]-1]=trk2->nbad_l3[m];
+	  ctrk.calfl[trk2->DSPnumber[m]-1]=trk2->cal_flag[m];
+	  ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk2->crc_cal[m][0];
+	  ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk2->crc_cal[m][1];
+	  ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk2->crc_cal[m][2];
+	  ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk2->crc_hcal[m];
+	  for (Int_t j = 0; j < 3072; j++){
+	    ctrk.dspped[trk2->DSPnumber[m]-1][j]=trk2->DSPped_par[m][j];
+	    ctrk.dspsig[trk2->DSPnumber[m]-1][j]=trk2->DSPsig_par[m][j];
+	    ctrk.dspbad[trk2->DSPnumber[m]-1][j]=trk2->DSPbad_par[m][j];
+	  }
+	}
       }
     }