--- tof/flight/ToFReprocessing/ToFreproc.cpp	2008/11/28 14:58:50	1.1.1.1
+++ tof/flight/ToFReprocessing/ToFreproc.cpp	2009/02/05 09:12:35	1.2
@@ -72,8 +72,9 @@
   //  read input file/list
   //  --------------------
   pam_event = new PamLevel2(ddir,list,options);
+  pam_event->SetSELLI(2);
   TTree::SetMaxTreeSize(1000*Long64_t(2000000000)); 
- //
+  //
   TString outfile_t =0;
   outfile_t = OUTF;
   outfile_t.Append("b.root");
@@ -145,7 +146,7 @@
   for(ULong64_t iev=0; iev<nevents; iev++){      
 
 
-//==================================================================
+    //==================================================================
 
 
     //  for(ULong64_t iev=0; iev<50; iev++){      
@@ -182,104 +183,108 @@
       //
 
 
-//================================================================
-//==================================================================
-//---  Define absolute time
-       Float_t tabs=pam_event->GetOrbitalInfo()->absTime;
-
-//==================================================================
-Float_t dedx_corr_m[2000][48],dedx_corr[48];
-Double_t mtime[2000],t1,t2,tm;
-Float_t yhelp1,yhelp2,slope,inter,thelp1,thelp2;
-Float_t xmean1,xwidth1;
-
-Int_t ical,ii,j,jj;
-
-if (iev==0) {
-
-ical=0;  // counter set to zero if first-time reading
-
-//-----------------------------------------------------------
-// Here I  read the dEdx_korr parameters
-//-----------------------------------------------------------
-
-jj=0;
-
-ifstream fin("adcmonitor.7th.100k.cut.dat");
-
-while (! fin.eof()) {
-fin>>t1>>tm>>t2;
-//cout<<jj<<" "<<tm<<endl;
-cout << setiosflags(ios::fixed)  << setw(10) << setprecision(3) << tm << endl;
-mtime[jj]=tm;
-for (ii=0; ii<48;ii++) {
-fin>>j>>xmean1>>xwidth1;
-dedx_corr_m[jj][ii]=xmean1;
-                       }
-jj=jj+1;
-}
-
-fin.close();
-
-
-while (tabs > mtime[ical]) {
-ical = ical+1;
-          }
-ical = ical-1;
-cout<<"abs time "<<tabs<<" limit low "<<mtime[ical]<<" limit up "<<mtime[ical+1]<<" ical "<<ical<<endl;
-} // if iev==0...
-
-//==================================================================
-//==  End of first time reading & filling the array
-//==================================================================
-
-//==================================================================
-//==  if time is outside time limits: 
-//==================================================================
-
-if (tabs<mtime[ical] || tabs>mtime[ical+1]) {
-cout<<"Checking Time Limits!"<<endl;
-ical=0;
-while (tabs > mtime[ical]) {
-ical = ical+1;
-          }
-ical = ical-1;
-cout<<"abs time "<<tabs<<" limit low "<<mtime[ical]<<" limit up "<<mtime[ical+1]<<" ical "<<ical<<endl;
-                                           };
-
-//==================================================================
-//== interpolate betwen time limits
-//==================================================================
-
-       thelp1 = mtime[ical];
-       thelp2 = mtime[ical+1];
-
-       for (ii=0; ii<48;ii++) {
-       yhelp1 = dedx_corr_m[ical][ii];
-       yhelp2 = dedx_corr_m[ical+1][ii];
-       slope  = (yhelp2-yhelp1)/(thelp2-thelp1);
-       inter  = yhelp1 - slope*thelp1;
-       dedx_corr[ii] = slope*tabs + inter;
-//       if (ii==0) cout<<thelp1<<" "<<thelp2<<" "<<tabs<<" "<<yhelp1<<" "<<yhelp2<<" "<<dedx_corr[0]<<endl;
-                       }
+      //================================================================
+      //==================================================================
+      //---  Define absolute time
+      Float_t tabs=pam_event->GetOrbitalInfo()->absTime;
+
+      //==================================================================
+      Float_t dedx_corr_m[2000][48],dedx_corr[48];
+      Double_t mtime[2000],t1,t2,tm;
+      Float_t yhelp1,yhelp2,slope,inter,thelp1,thelp2;
+      Float_t xmean1,xwidth1;
+
+      Int_t ical,ii,j,jj;
+
+      if (iev==0) {
+
+	ical=0;  // counter set to zero if first-time reading
+
+	//-----------------------------------------------------------
+	// Here I  read the dEdx_korr parameters
+	//-----------------------------------------------------------
+
+	jj=0;
+
+	ifstream fin("adcmonitor.7th.100k.cut.dat");
+
+	while (! fin.eof()) {
+	  fin>>t1>>tm>>t2;
+	  //	  cout<<jj<<" "<<tm<<endl;
+	  cout << setiosflags(ios::fixed)  << setw(10) << setprecision(3) << tm << endl;
+	  mtime[jj]=tm;
+	  for (ii=0; ii<48;ii++) {
+	    fin>>j>>xmean1>>xwidth1;
+	    dedx_corr_m[jj][ii]=xmean1;
+	  }
+	  jj=jj+1;
+	  //	  printf(" kk %i \n",jj);
+	}
+	//	printf(" 1ical %i \n",ical);
+
+	fin.close();
+
+
+	while (tabs<mtime[ical] || tabs > mtime[ical+1]) {
+	  //	  printf(" ical %i \n",ical);
+	  ical = ical+1;
+	}
+	//	ical = ical-1;
+	cout<<"abs time "<<tabs<<" limit low "<<mtime[ical]<<" limit up "<<mtime[ical+1]<<" ical "<<ical<<endl;
+      } // if iev==0...
+
+      //==================================================================
+      //==  End of first time reading & filling the array
+      //==================================================================
+
+      //==================================================================
+      //==  if time is outside time limits: 
+      //==================================================================
+
+      if (tabs<mtime[ical] || tabs>mtime[ical+1]) {
+	cout<<"Checking Time Limits!"<<endl;
+	ical=0;
+	while (tabs > mtime[ical+1] || tabs<mtime[ical]) {
+	  ical = ical+1;
+	}
+	//	ical = ical-1;
+	cout<<"abs time "<<tabs<<" limit low "<<mtime[ical]<<" limit up "<<mtime[ical+1]<<" ical "<<ical<<endl;
+      };
 
-//================================================================
-//================================================================
+      //      printf(" 2ical %i \n",ical);
+      //==================================================================
+      //== interpolate betwen time limits
+      //==================================================================
+
+      thelp1 = mtime[ical];
+      thelp2 = mtime[ical+1];
+
+      for (ii=0; ii<48;ii++) {
+	yhelp1 = dedx_corr_m[ical][ii];
+	yhelp2 = dedx_corr_m[ical+1][ii];
+	slope  = (yhelp2-yhelp1)/(thelp2-thelp1);
+	inter  = yhelp1 - slope*thelp1;
+	dedx_corr[ii] = slope*tabs + inter;
+	//       if (ii==0) cout<<thelp1<<" "<<thelp2<<" "<<tabs<<" "<<yhelp1<<" "<<yhelp2<<" "<<dedx_corr[0]<<endl;
+      }
 
+      //================================================================
+      //================================================================
+      //      printf("cpippo \n");
 
       // process tof data
       //
       for (Int_t hh=0; hh<12;hh++){
         for (Int_t kk=0; kk<4;kk++){
-		adc[kk][hh] = 4095;
-		tdc[kk][hh] = 4095;
-		tdcc[kk][hh] = 4095.;
-		tofinput_.adc[hh][kk] = 4095;
-		tofinput_.tdc[hh][kk] = 4095;
+	  adc[kk][hh] = 4095;
+	  tdc[kk][hh] = 4095;
+	  tdcc[kk][hh] = 4095.;
+	  tofinput_.adc[hh][kk] = 4095;
+	  tofinput_.tdc[hh][kk] = 4095;
 	};
       };
-//      memset(adc, 0, 12*4*sizeof(Int_t));
-//      memset(tdc, 0, 12*4*sizeof(Int_t));
+      //      memset(adc, 0, 12*4*sizeof(Int_t));
+      //      memset(tdc, 0, 12*4*sizeof(Int_t));
       Int_t gg = 0;
       Int_t hh = 0;
       Int_t adcf[48];
@@ -289,10 +294,10 @@
       for (Int_t pm=0; pm < tofl2->ntrk() ; pm++){
 	ToFTrkVar *ttf = tofl2->GetToFTrkVar(pm);
 	for ( Int_t nc=0; nc < ttf->npmttdc; nc++){
-		if ( (ttf->tdcflag).At(nc) != 0 ) tdcf[(ttf->pmttdc).At(nc)] = 1;
+	  if ( (ttf->tdcflag).At(nc) != 0 ) tdcf[(ttf->pmttdc).At(nc)] = 1;
 	};
 	for ( Int_t nc=0; nc < ttf->npmtadc; nc++){
-		if ( (ttf->adcflag).At(nc) != 0 ) adcf[(ttf->pmtadc).At(nc)] = 1;
+	  if ( (ttf->adcflag).At(nc) != 0 ) adcf[(ttf->pmtadc).At(nc)] = 1;
 	};
       };
       //
@@ -300,12 +305,12 @@
 	ToFPMT *pmt = tofl2->GetToFPMT(pm);
 	tofl2->GetPMTIndex(pmt->pmt_id, gg, hh);
 	if ( adcf[pmt->pmt_id] == 0 ){
-		tofinput_.adc[gg][hh] = (int)pmt->adc;
-		adc[hh][gg] = (int)pmt->adc;
+	  tofinput_.adc[gg][hh] = (int)pmt->adc;
+	  adc[hh][gg] = (int)pmt->adc;
 	};
 	if ( tdcf[pmt->pmt_id] == 0 ){
-		tofinput_.tdc[gg][hh] = (int)pmt->tdc;
-		tdc[hh][gg] = (int)pmt->tdc;
+	  tofinput_.tdc[gg][hh] = (int)pmt->tdc;
+	  tdc[hh][gg] = (int)pmt->tdc;
 	};
 	tdcc[hh][gg] = (float)pmt->tdc_tw;
 	//        Int_t pppid = tof->GetPMTid(hh,gg);
@@ -313,14 +318,14 @@
       };
 
 
-//      for (Int_t pm=0; pm <48 ; pm++){
-//	tof->GetPMTIndex(pm, gg, hh);
-//	tofinput_.tdc[hh][gg] = (int)500.;
-//	tofinput_.adc[hh][gg] = (int)500.;
-//	tdc[hh][gg] = (int)500.;
-//	adc[hh][gg] = (int)500.;
-//	//	printf(" hh %i gg %i tdc %f adc %f \n",hh,gg,pmt->tdc,pmt->adc);
-//    };
+      //      for (Int_t pm=0; pm <48 ; pm++){
+      //	tof->GetPMTIndex(pm, gg, hh);
+      //	tofinput_.tdc[hh][gg] = (int)500.;
+      //	tofinput_.adc[hh][gg] = (int)500.;
+      //	tdc[hh][gg] = (int)500.;
+      //	adc[hh][gg] = (int)500.;
+      //	//	printf(" hh %i gg %i tdc %f adc %f \n",hh,gg,pmt->tdc,pmt->adc);
+      //    };
       //
       for (Int_t hh=0; hh<5;hh++){
 	tofinput_.patterntrig[hh]=pam_event->GetTrigLevel2()->patterntrig[hh];
@@ -397,8 +402,8 @@
       //
       for (Int_t hh=0; hh<12;hh++){
 	for (Int_t kk=0; kk<4;kk++){
-         // new WM
-//	  if ( tofoutput_.tdc_c[hh][kk] < 4095 || adc[kk][hh] < 4095  || tdc[kk][hh] < 4095 ){	   
+	  // new WM
+	  //	  if ( tofoutput_.tdc_c[hh][kk] < 4095 || adc[kk][hh] < 4095  || tdc[kk][hh] < 4095 ){	   
 	  if ( tdcc[kk][hh] < 4095. || adc[kk][hh] < 4095  || tdc[kk][hh] < 4095 ){	   
 	    //
 	    t_pmt->pmt_id = tof->GetPMTid(kk,hh);
@@ -435,7 +440,7 @@
 	  // Get tracker related variables for this track
 	  //
 	  toftrk();
-//	  toftrk(thelp);
+	  //	  toftrk(thelp);
 	  //
 	  // Copy values in the class from the structure (we need to use a temporary class to store variables).
 	  //
@@ -458,7 +463,7 @@
 	  for (Int_t hh=0; hh<12;hh++){
 	    for (Int_t kk=0; kk<4;kk++){
 	      if ( tofoutput_.adc_c[hh][kk] < 1000 ){
-  	     //	    t_tof->dedx.AddAt(tofoutput_.adc_c[hh][kk],t_tof->npmtadc); // EMILIANO
+		//	    t_tof->dedx.AddAt(tofoutput_.adc_c[hh][kk],t_tof->npmtadc); // EMILIANO
 	        pmt_id = tof->GetPMTid(kk,hh); 
   	        t_tof->dedx.AddAt((tofoutput_.adc_c[hh][kk]*4./dedx_corr[pmt_id]),t_tof->npmtadc); // EMILIANO
 		t_tof->pmtadc.AddAt(pmt_id,t_tof->npmtadc);
@@ -526,87 +531,87 @@
 
 void usage(){
 
-    cout << "------------------------------------------------------------"<<endl;
-    cout << "Loop over events (on one or more Level2-files), applying some selection cuts (defined in My-Selection.cpp), \n";
-    cout << "creating output histograms (defined in My-Histos.cpp) and/or trees with selected events. \n \n ";
-    cout << "USAGE:"<<endl;
-    cout << "-processDir  DIR     -  Level2 data directory \n";
-    cout << "-processList LIST    -  list  of files (.txt) or single file (.root) to be analysed \n";
-    cout << "-outputFile  PATH    -  name of the output file \n";
-    cout << "-NumEvents   XXX     -  number of events to be analysed \n";
-    cout << "--debug, -g          -  debug mode \n";
-    cout << "--help, -h           -  print this help \n";
-    cout << "-options [ options ] -  options: \n";
-    cout << "                        fillHistos --> create an output file with histograms  \n";
-    cout << "                        fillTree   --> create an output file with trees storing the selected events \n ";
-    cout << "                        +(-)ALL    --> inlcude(exclude) all trees and branches \n "   ;
-    cout << "                        +(-)TRK1 +(-)TRK2 +(-)CAL1 +(-)CAL2 +(-)TOF +(-)TRG +(-)ND +(-)S4 +(-)ORB --> inlcude(exclude) trees and branches  \n"  ; 
-    cout << "------------------------------------------------------------"<<endl;
+  cout << "------------------------------------------------------------"<<endl;
+  cout << "Loop over events (on one or more Level2-files), applying some selection cuts (defined in My-Selection.cpp), \n";
+  cout << "creating output histograms (defined in My-Histos.cpp) and/or trees with selected events. \n \n ";
+  cout << "USAGE:"<<endl;
+  cout << "-processDir  DIR     -  Level2 data directory \n";
+  cout << "-processList LIST    -  list  of files (.txt) or single file (.root) to be analysed \n";
+  cout << "-outputFile  PATH    -  name of the output file \n";
+  cout << "-NumEvents   XXX     -  number of events to be analysed \n";
+  cout << "--debug, -g          -  debug mode \n";
+  cout << "--help, -h           -  print this help \n";
+  cout << "-options [ options ] -  options: \n";
+  cout << "                        fillHistos --> create an output file with histograms  \n";
+  cout << "                        fillTree   --> create an output file with trees storing the selected events \n ";
+  cout << "                        +(-)ALL    --> inlcude(exclude) all trees and branches \n "   ;
+  cout << "                        +(-)TRK1 +(-)TRK2 +(-)CAL1 +(-)CAL2 +(-)TOF +(-)TRG +(-)ND +(-)S4 +(-)ORB --> inlcude(exclude) trees and branches  \n"  ; 
+  cout << "------------------------------------------------------------"<<endl;
 }
 //
 int HandleInputPar(int argc, char **argv){
 
-    if(argc>1){
+  if(argc>1){
 		
-	if(!strcmp(argv[1], "-h") || !strcmp(argv[1], "--help") ){
-	    usage();
-	    return(1);
-	};
-// -----------------------
-// Read input parameters
-// -----------------------  
-	DIR     = gSystem->WorkingDirectory();
-	LIST    = "";
-	OUTF    = "myfile";
-	OPTIONS = "+AUTO -TOF";
+    if(!strcmp(argv[1], "-h") || !strcmp(argv[1], "--help") ){
+      usage();
+      return(1);
+    };
+    // -----------------------
+    // Read input parameters
+    // -----------------------  
+    DIR     = gSystem->WorkingDirectory();
+    LIST    = "";
+    OUTF    = "myfile";
+    OPTIONS = "+AUTO -TOF";
 		
-	for (int i = 1; i < argc; i++){		
-	    // -----------------------------------------------------//
-	    if (!strcmp(argv[i], "-processDir")){
-		if (++i >= argc) throw -1;
-		DIR = argv[i];
-		cout << "processDir "<<DIR<<endl;
-		continue;
-	    }
-	    // -----------------------------------------------------//
-	    else if (!strcmp(argv[i], "-processList")){
-		if (++i >= argc) throw -1;
-		LIST = argv[i];
-		cout << "processList "<<LIST<<endl;
-		continue;
-	    }  
-	    // -----------------------------------------------------//
-	    else if (!strcmp(argv[i], "-outputFile")){
-		if (++i >= argc) throw -1;
-		OUTF = argv[i];
-		cout << "outputFile "<<OUTF<<endl;
-		continue;
-	    }  
-	    // -----------------------------------------------------//
-	    else if (!strcmp(argv[i], "-options")){
-		if (++i >= argc) throw -1;
-		OPTIONS = argv[i];
-		if( OPTIONS.Contains("[") ){
-		    do{
-			if (++i >= argc) throw -1;
-			OPTIONS.Append(argv[i]);
-		    }while(!OPTIONS.Contains("]"));
-		}else cout << "wrong option format --> ignoring " << endl;
-	    }
-	    else{
-		cout << "Unidentified input parameter. Ignored."<< endl;
-	    };
-	};
-    }else{
-	usage();
-	return(1);
+    for (int i = 1; i < argc; i++){		
+      // -----------------------------------------------------//
+      if (!strcmp(argv[i], "-processDir")){
+	if (++i >= argc) throw -1;
+	DIR = argv[i];
+	cout << "processDir "<<DIR<<endl;
+	continue;
+      }
+      // -----------------------------------------------------//
+      else if (!strcmp(argv[i], "-processList")){
+	if (++i >= argc) throw -1;
+	LIST = argv[i];
+	cout << "processList "<<LIST<<endl;
+	continue;
+      }  
+      // -----------------------------------------------------//
+      else if (!strcmp(argv[i], "-outputFile")){
+	if (++i >= argc) throw -1;
+	OUTF = argv[i];
+	cout << "outputFile "<<OUTF<<endl;
+	continue;
+      }  
+      // -----------------------------------------------------//
+      else if (!strcmp(argv[i], "-options")){
+	if (++i >= argc) throw -1;
+	OPTIONS = argv[i];
+	if( OPTIONS.Contains("[") ){
+	  do{
+	    if (++i >= argc) throw -1;
+	    OPTIONS.Append(argv[i]);
+	  }while(!OPTIONS.Contains("]"));
+	}else cout << "wrong option format --> ignoring " << endl;
+      }
+      else{
+	cout << "Unidentified input parameter. Ignored."<< endl;
+      };
     };
-// -----------------------
-// Check input parameters
-// -----------------------  
+  }else{
+    usage();
+    return(1);
+  };
+  // -----------------------
+  // Check input parameters
+  // -----------------------  
 	
 		
-    return(0);
+  return(0);
 	
 };
 // 
@@ -614,14 +619,14 @@
 int main(int argc, char **argv)
 {
 	
-    if( HandleInputPar(argc,argv) )return(1);
+  if( HandleInputPar(argc,argv) )return(1);
 	
-    cout << "OPTIONS "<<OPTIONS<<endl;
-    Loop(DIR,LIST,OPTIONS);
+  cout << "OPTIONS "<<OPTIONS<<endl;
+  Loop(DIR,LIST,OPTIONS);
     
-    cout << "Back to main - end"<<endl;
+  cout << "Back to main - end"<<endl;
     
-    return 0;
+  return 0;
 	
 };