--- PamelaLevel2/src/PamLevel2.cpp	2008/04/03 14:39:51	1.61
+++ PamelaLevel2/src/PamLevel2.cpp	2008/11/06 10:45:36	1.65
@@ -284,6 +284,7 @@
     trk_track  = 0;
     calo_track = 0;
     tof_track  = 0;
+    orb_track  = 0;
     candeleteobj = 0;
     pscore = 0;
     iscore = 0;
@@ -295,17 +296,19 @@
 /**
  * Constructor 
  */
-PamTrack::PamTrack(TrkTrack* t, CaloTrkVar* c, ToFTrkVar* o){ 
+PamTrack::PamTrack(TrkTrack* t, CaloTrkVar* c, ToFTrkVar* o, OrbitalInfoTrkVar *r){ 
 
     trk_track  = 0;
     calo_track = 0;
     tof_track  = 0;
+    orb_track  = 0;
 //     if(t)trk_track  = new TrkTrack(*t);
 //     if(c)calo_track = new CaloTrkVar(*c);
 //     if(o)tof_track  = new ToFTrkVar(*o);
     if(t)trk_track  = t;
     if(c)calo_track = c;
     if(o)tof_track  = o;
+    if(r)orb_track  = r;
     candeleteobj = 0;
 
 };
@@ -314,13 +317,16 @@
     TrkTrack   *t = track.trk_track;
     CaloTrkVar *c = track.calo_track;
     ToFTrkVar  *o = track.tof_track;
+    OrbitalInfoTrkVar  *r = track.orb_track;
 
     trk_track  = 0;
     calo_track = 0;
     tof_track  = 0;
+    orb_track  = 0;
     if(t)trk_track  = new TrkTrack(*t);
     if(c)calo_track = new CaloTrkVar(*c);
     if(o)tof_track  = new ToFTrkVar(*o);
+    if(r)orb_track  = new OrbitalInfoTrkVar(*r);
     candeleteobj = 1;
     
 }
@@ -331,10 +337,12 @@
 	if(trk_track)  trk_track->TrkTrack::Clear();
 	if(calo_track) calo_track->CaloTrkVar::Clear();//???
 	if(tof_track)  tof_track->ToFTrkVar::Clear();//???
+	if(orb_track)  orb_track->OrbitalInfoTrkVar::Clear();//???
     }else{
 	trk_track  = 0;
 	calo_track = 0;
 	tof_track  = 0;	
+	orb_track  = 0;	
     }
     pscore = 0;
     iscore = 0;
@@ -355,6 +363,10 @@
 	    tof_track->ToFTrkVar::Clear();//???
 	    delete tof_track;
 	}
+	if(orb_track){
+  	    orb_track->OrbitalInfoTrkVar::Clear();//???
+	    delete orb_track;
+	}
     }else{
 	Clear();
     }
@@ -710,7 +722,7 @@
 
 Bool_t PamLevel2::IsGood(){
     Bool_t goodev=true;
-    //  if(trk2_obj && trk2_obj->UnpackError() != 0 ) goodev = false; 
+    if(trk2_obj && trk2_obj->UnpackError() != 0 ) goodev = false; 
     if(calo2_obj && calo2_obj->good != 1) goodev = false;
     if(tof_obj && tof_obj->unpackError != 0) goodev = false;  
     if(trig_obj && trig_obj->unpackError != 0) goodev = false; 
@@ -744,7 +756,10 @@
 	return &trk2_obj;
     };
     if(!objname.CompareTo("TrkHough"))   {
-	if(!trkh_obj)  trkh_obj  = new TrkHough();
+	if(!trkh_obj)  {
+	    trkh_obj  = new TrkHough();
+	    trkh_obj->Set();
+	}
 	return &trkh_obj;
     };
     if(!objname.CompareTo("CaloLevel1")) {
@@ -866,7 +881,40 @@
 //
 //--------------------------------------
 /**
- * Give the pamela track associated to a tracker track, retrieving related calorimeter and tof track information. 
+ * Retrieves the OrbitalInfo track matching the seqno-th tracker stored track. 
+ * (If seqno = -1 retrieves the tracker-independent tof related track)
+ */
+OrbitalInfoTrkVar *PamLevel2::GetOrbitalInfoStoredTrack(int seqno){
+	 
+    if( !orb_obj )return 0;
+
+    if( orb_obj->OrbitalInfo::ntrk()==0 ){
+	cout << "PamLevel2::GetOrbitalInfoStoredTrack(int) : requested tracker SeqNo "<< seqno <<" but no OrbitalInfo tracks are stored"<<endl;
+	return NULL;
+    };
+     
+    OrbitalInfoTrkVar *c = 0;
+    Int_t it_tof=0;
+     
+    do{ 
+        c = orb_obj->OrbitalInfo::GetOrbitalInfoTrkVar(it_tof);
+	it_tof++; 
+    } while( c && seqno != c->trkseqno && it_tof < orb_obj->OrbitalInfo::ntrk());	
+     
+    if(!c || seqno != c->trkseqno){
+	c = 0;
+	if(seqno!=-1)cout << "PamLevel2::GetOrbitalInfoStoredTrack(int) : requested tracker SeqNo "<< seqno <<" does not match OrbitalInfo stored tracks"<<endl;
+    };
+    return c;
+     
+};
+
+//--------------------------------------
+//
+//
+//--------------------------------------
+/**
+ * Give the pamela track associated to a tracker track, retrieving related calorimeter, orbitalinfo and tof track information. 
  */
 PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t){
      
@@ -875,12 +923,14 @@
 
     CaloTrkVar *c = 0;
     ToFTrkVar  *o = 0;
+    OrbitalInfoTrkVar  *r = 0;
      
     if(CAL2) c = GetCaloStoredTrack(t->GetSeqNo());
     if(TOF) o = GetToFStoredTrack(t->GetSeqNo());
+    if(ORB) r = GetOrbitalInfoStoredTrack(t->GetSeqNo());
      
 //    if(t && c && o)track = new PamTrack(t,c,o);
-    PamTrack *track = new PamTrack(t,c,o);
+    PamTrack *track = new PamTrack(t,c,o,r);
      
     return track;
  
@@ -898,7 +948,7 @@
 PamTrack* PamLevel2::GetStoredTrack(Int_t itrk){
     
     cout <<"PamLevel2::GetStoredTrack(Int_t itrk) **to-be-updated** "<<endl;
-    cout <<"for the moment, better use separately the methods: TrkLevel2::GetStoredTrack(seqno) CaloLevel2::GetCaloTrkVar(Int_t notrack) ToFLevel2::GetToFTrkVar(Int_t notrack)"<<endl;
+    cout <<"for the moment, better use separately the methods: TrkLevel2::GetStoredTrack(seqno) CaloLevel2::GetCaloTrkVar(Int_t notrack) ToFLevel2::GetToFTrkVar(Int_t notrack) OrbitalInfo::GetOrbitalInfoTrkVar(Int_t notrack)"<<endl;
     cout <<"(if you use it, remember to delete the PamTrack object)"<<endl;
     PamTrack *track = 0;
     
@@ -1010,15 +1060,18 @@
 	TrkTrack   *ts = 0;
 	CaloTrkVar *cs = 0;
 	ToFTrkVar  *os = 0;
+	OrbitalInfoTrkVar  *rs = 0;
 	
 	// get tracker tracks
 	TrkTrack   *tp = trk2_obj->TrkLevel2::GetTrack(i); //tracker
 	CaloTrkVar *cp = GetCaloStoredTrack(tp->GetSeqNo());
 	ToFTrkVar  *op = GetToFStoredTrack(tp->GetSeqNo());
+	OrbitalInfoTrkVar  *rp = GetOrbitalInfoStoredTrack(tp->GetSeqNo());
 
 	TrkTrack   *ti = 0; //tracker (image)
 	CaloTrkVar *ci = 0;
 	ToFTrkVar  *oi = 0;
+	OrbitalInfoTrkVar  *ri = 0;
 	//	cout << "trk track n. "<<i << " "<<hex<< tp <<dec<< endl;
 	// if track has an image, check image selection
 
@@ -1032,6 +1085,7 @@
 	    ti = trk2_obj->TrkLevel2::GetTrackImage(i);              //tracker (image)
 	    ci = GetCaloStoredTrack(ti->GetSeqNo());
 	    oi = GetToFStoredTrack(ti->GetSeqNo());
+	    ri = GetOrbitalInfoStoredTrack(ti->GetSeqNo());
 	    
 	    //	    cout << "its image "<<i << " "<<hex<< ti <<dec<< endl;
 
@@ -1373,15 +1427,18 @@
 		ts = ti;//its image!!
 		cs = ci;
 		os = oi;
+		rs = ri;
 		Int_t totis = toti_score;
 
 		ti = tp;//its image!!
 		ci = cp;
 		oi = op;
+		ri = rp;
 
 		tp = ts;//its image!!
 		cp = cs;
 		op = os;
+		rp = rs;
 
 		toti_score = totp_score;
 		totp_score = totis;
@@ -1411,8 +1468,8 @@
 	//	cout<<"o "<<cp<<endl;
 	//	cout<<"o "<<op<<endl;
 
-	new(ttsorted[i]) PamTrack(tp,cp,op);
-	new(ttimage[i])  PamTrack(ti,ci,oi);
+	new(ttsorted[i]) PamTrack(tp,cp,op,rp);
+	new(ttimage[i])  PamTrack(ti,ci,oi,ri);
 
 	((PamTrack*)(ttsorted[i]))->SetPScore(totp_score);
 	((PamTrack*)(ttsorted[i]))->SetIScore(toti_score);
@@ -2532,26 +2589,10 @@
     if( GetOrbitalInfo() ){
 	abstime = GetOrbitalInfo()->absTime;
 	obt     = GetOrbitalInfo()->OBT;
-//	cout << " ABS >>> "<<abstime<<" OBT >>> "<<obt<<" "<<endl;
     }else{
   	abstime = GetRunInfo()->RUNHEADER_TIME;
  	obt     = GetRunInfo()->RUNHEADER_OBT;
     }
-//    cout << ISGP << " "<<abstime<<endl;
-    // ---------------------------------------------------------------
-    // the absolute time is necessary to relate the event with the run
-    // ---------------------------------------------------------------
-//     if( !GetOrbitalInfo() && !ISGP ){
-// 	cout << "Bool_t PamLevel2::UpdateRunInfo(Long64_t iev) -- ERROR -- missing OrbitalInfo "<<endl;
-// 	return(false);
-//     }
-    // -----------------------------------------------------------------------
-    // if it is simulation, assign abstime by hand (temporaneo!!!)
-    // -----------------------------------------------------------------------
-//     if(ISGP){
-//  	abstime = GetRunInfo()->RUNHEADER_TIME; // BARBATRUCCO
-//  	obt     = GetRunInfo()->RUNHEADER_OBT; // BARBATRUCCO
-//     }
     
 
     // +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-
@@ -2570,72 +2611,41 @@
 	    }
 	    totdltime[2]++;
 	}
-//	cout << setw(10)<<totdltime[0]<<setw(10)<<totdltime[1]<<setw(10)<<totdltime[2]<<endl;
 
-// 	// ---------------------------------------------------------------
-// 	// retrieve OBT and absolute time of the event
-// 	// ---------------------------------------------------------------
-// 	ULong64_t abstime = 0;
-// 	ULong64_t obt     = 0;
-// 	if( GetOrbitalInfo() ){
-// 	    abstime = GetOrbitalInfo()->absTime;
-// 	    obt     = GetOrbitalInfo()->OBT;
-// 	}
-// 	// ---------------------------------------------------------------
-// 	// the absolute time is necessary to relate the event with the run
-// 	// ---------------------------------------------------------------
-// 	if( !GetOrbitalInfo() && !ISGP ){
-// 	    cout << "Bool_t PamLevel2::UpdateRunInfo(Long64_t iev) -- ERROR -- missing OrbitalInfo "<<endl;
-// 	    return(false);
-// 	}
-      
-// 	// -----------------------------------------------------------------------
-// 	// the first time the routine is called, set run search from the beginning
-// 	// -----------------------------------------------------------------------
-// 	if ( irun < 0 ){
-// 	    irun = 0LL;
-// 	    irunt = 0LL;
-// 	    irunentry = 0;
-// 	    prevshift = 0;
-// 	    run_tree->GetEntry(irun);
-
-// 	    if( ISGP && run_tree->GetEntries()!=1 ){
-// 		cout << "** WARNING ** simulated files are assumed to have 1 single run, not "<< run_tree->GetEntries() << endl;
-// 		cout << "** WARNING ** run will not be updated"<<endl;
-// 	    }
-// 	};      
-
-// 	// -----------------------------------------------------------------------
-// 	// if it is simulation, assign abstime by hand (temporaneo!!!)
-// 	// -----------------------------------------------------------------------
-// 	if(ISGP){
-// 	    abstime = GetRunInfo()->RUNHEADER_TIME; // BARBATRUCCO
-// 	    obt     = GetRunInfo()->RUNHEADER_OBT; // BARBATRUCCO
-// 	}
 	//
 	bool fromfirst = true; // first loop over runs
 
-//	Bool_t hasfrag = false;
-//	if( GetRunInfo()->ID_RUN_FRAG!=0 && GetRunInfo()->ID_RUN_FRAG != GetRunInfo()->ID)hasfrag=true;
-//	ULong64_t fragid = GetRunInfo()->ID_RUN_FRAG;
 
 	// ------------------------------------------------------
 	// loop over runs to find the one that contains the event 
 	// ------------------------------------------------------
 	while ( 
 	    (
-	        (
-		    !(abstime >= GetRunInfo()->RUNHEADER_TIME &&    // check on absolute time (s)
-		      abstime <= GetRunInfo()->RUNTRAILER_TIME) && 
-		    !(obt >= GetRunInfo()->RUNHEADER_OBT &&         // additional check on OBT (ms)
-		      obt <= GetRunInfo()->RUNTRAILER_OBT)
-		    ) 
-  		|| GetRunInfo()->NEVENTS==0 
+// 	        (
+// 		    !(abstime >= GetRunInfo()->RUNHEADER_TIME &&    // check on absolute time (s)
+// 		      abstime <= GetRunInfo()->RUNTRAILER_TIME) && 
+// 		    !(obt >= GetRunInfo()->RUNHEADER_OBT &&         // additional check on OBT (ms)
+// 		      obt <= GetRunInfo()->RUNTRAILER_OBT)
+// 		    ) 
+	        
+		!(abstime >= GetRunInfo()->RUNHEADER_TIME &&    // check on absolute time (s)
+		  abstime <= GetRunInfo()->RUNTRAILER_TIME) 
+		||
+		!(obt >= GetRunInfo()->RUNHEADER_OBT &&         // additional check on OBT (ms)
+		  obt <= GetRunInfo()->RUNTRAILER_OBT)		 
+  		|| 
+		GetRunInfo()->NEVENTS==0 
 //		|| !(irunentry < GetRunInfo()->NEVENTS-1-prevshift) // ERRORE!!! fa saltare i run con 1 evento
-		|| !(irunentry <= GetRunInfo()->NEVENTS-1-prevshift)
+		|| 
+		!(irunentry <= GetRunInfo()->NEVENTS-1-prevshift)
 		)
 	    && irun < run_tree->GetEntries() ){
 
+
+// 	    if( !(abstime >= GetRunInfo()->RUNHEADER_TIME &&abstime <= GetRunInfo()->RUNTRAILER_TIME)  )cout << "ABS TIME "<<abstime << " " <<GetRunInfo()->RUNTRAILER_TIME <<endl;
+// 	    if( !(obt >= GetRunInfo()->RUNHEADER_OBT && obt <= GetRunInfo()->RUNTRAILER_OBT) )cout << "OBT TIME "<< obt <<" "<< GetRunInfo()->RUNTRAILER_OBT << endl;
+// 	    if( GetRunInfo()->NEVENTS==0  )cout <<"GetRunInfo()->NEVENTS==0 "<<endl;
+// 	    if( !(irunentry <= GetRunInfo()->NEVENTS-1-prevshift) ) cout << "irunentry > "<<GetRunInfo()->NEVENTS-1-prevshift << endl;
 	    // - - - - - - - - - - - - -
 	    // irunentry = position of current entry, relative to the run
             // prevshift = shift needed to synchronize l0 and l2 data (nested events)
@@ -2645,35 +2655,35 @@
 	    // store dead and live-time of previous run
 	    // -----------------------------------------
 //	    if(SELLI==0){
-		if(fromfirst){
-		    if(oldrun==irun){
-			/// decrement counters
-			if( GetTrigLevel2()){
-			    totdltime[0]-=GetTrigLevel2()->dltime[0];//live-time
-			    totdltime[1]-=GetTrigLevel2()->dltime[1];//dead-time
-			}
-			totdltime[2]--;                              //event counter
-			cout << endl;
-			cout << "n.events     : "<<totdltime[2]<<endl;
-			cout << "RUN LIVE-TIME: "<<totdltime[0]*0.16<<" ms"<<endl;
-			cout << "RUN DEAD-TIME: "<<totdltime[1]*0.01<<" ms"<<endl;	    
-		    }else{
-			totdltime[0]=0;//live-time
-			totdltime[1]=0;//dead-time
-			totdltime[2]=0;                             //event counter
-			cout << " *** JUMP RUN *** irun "<<irun<<endl; 
+	    if(fromfirst){
+		if(oldrun==irun){
+		    /// decrement counters
+		    if( GetTrigLevel2()){
+			totdltime[0]-=GetTrigLevel2()->dltime[0];//live-time
+			totdltime[1]-=GetTrigLevel2()->dltime[1];//dead-time
 		    }
-		    /// add an entry
-		    if(run_tree_clone)
-			if(run_tree_clone->GetBranch("DeadLiveTime")->GetEntries() < run_tree->GetEntries())
-			    run_tree_clone->GetBranch("DeadLiveTime")->Fill();
-		    /// reset counters
-		    if( GetTrigLevel2() ){
-			totdltime[0]=GetTrigLevel2()->dltime[0];//live-time
-			totdltime[1]=0;                         //dead-time
-		    }
-		    totdltime[2]=1;                             //event counter
+		    totdltime[2]--;                              //event counter
+		    cout << endl;
+		    cout << "n.events     : "<<totdltime[2]<<endl;
+		    cout << "RUN LIVE-TIME: "<<totdltime[0]*0.16<<" ms"<<endl;
+		    cout << "RUN DEAD-TIME: "<<totdltime[1]*0.01<<" ms"<<endl;	    
+		}else{
+		    totdltime[0]=0;//live-time
+		    totdltime[1]=0;//dead-time
+		    totdltime[2]=0;                             //event counter
+		    cout << " *** JUMP RUN *** irun "<<irun<<endl; 
+		}
+		/// add an entry
+		if(run_tree_clone)
+		    if(run_tree_clone->GetBranch("DeadLiveTime")->GetEntries() < run_tree->GetEntries())
+			run_tree_clone->GetBranch("DeadLiveTime")->Fill();
+		/// reset counters
+		if( GetTrigLevel2() ){
+		    totdltime[0]=GetTrigLevel2()->dltime[0];//live-time
+		    totdltime[1]=0;                         //dead-time
 		}
+		totdltime[2]=1;                             //event counter
+	    }
 //	    }
 
 	    irun++;	   
@@ -2704,7 +2714,7 @@
 	    if(irun>0)runfirstentry += (GetRunInfo()->NEVENTS)-prevshift;
 	    irunentry = 0;
 	    prevshift = 0;	    
-	    run_tree->GetEntry(irun);
+	    run_tree->GetEntry(irun);//update runinfo
 	    irunt = irun - irunoffset[run_tree->GetTreeNumber()];
 	    if(GetRunInfo()->RUNHEADER_OBT>GetRunInfo()->RUNTRAILER_OBT ){
 		cout << "Bool_t PamLevel2::UpdateRunInfo(Long64_t iev) -- WARNING -- irun "<<irun<<"  has RUNHEADER_OBT>=RUNTRAILER_OBT " <<endl;
@@ -2727,12 +2737,18 @@
 	// --------------------------------------
 
 
-
-
 	// --------------------------------------
 	// ---> exit with TRUE
 	// --------------------------------------
-	cout << endl << " ))))) UPDATE RUN INFO (((((  @iev "<<iev<<" run "<<GetRunInfo()->ID<<" irun "<<irun<<endl;
+	if(DBG)cout << endl << " ))))) UPDATE RUN INFO (((((  @iev "<<iev<<" run "<<GetRunInfo()->ID<<" irun "<<irun<<endl;
+	// ----------------------------------------------------
+	// update the tracker parameters 
+	// (non ho trovato nessun altro modo sicuro di farlo...)
+	// ----------------------------------------------------
+	if(!dbc || (dbc && !dbc->IsConnected()) )SetDBConnection();
+	TrkParams::Set(GetRunInfo() ,dbc);
+	if(dbc)dbc->Close();
+
 	// ----------------------------------------------------
 	// then check if the run has a fragment 
 	// in this case we have to switch to the next fragment
@@ -2741,6 +2757,7 @@
 	if( 
 	    GetRunInfo()->ID_RUN_FRAG != 0 && 
 //	    GetRunInfo()->ID_RUN_FRAG != GetRunInfo()->ID && 
+	    DBG &&
 	    true ){
 	    cout << "* fragment *"<<endl;		
 	} 
@@ -2760,7 +2777,7 @@
 	    // check if the time is ok (with merged files it is not...)
 	    // if not loop over run and look for the proper entry
 	    bool SECONDO_GIRO=false;
-	    Long64_t irun_start   = irun;
+	    //	    Long64_t irun_start   = irun;
 	    int      offset_start = irunoffset[sel_tree->GetTreeNumber()];
 	    while ( 
 		(
@@ -2786,8 +2803,6 @@
 		    cout << " obt "<<obt;
 		    cout <<" >> "<<GetRunInfo()->RUNHEADER_OBT<<" "<<GetRunInfo()->RUNTRAILER_OBT;
 		    cout << " *** JUMP RUN *** irun "<<irun;
-// 		if(!SECONDO_GIRO)cout << " (don't worry)";
-// 		else             cout << " (start worring...)";
 		    cout << endl; 
 		}
 //		irun++;
@@ -2819,8 +2834,15 @@
 		cout << " obt "<<obt;
 		cout <<" >> "<<GetRunInfo()->RUNHEADER_OBT<<" "<<GetRunInfo()->RUNTRAILER_OBT;
 	    }
-	    cout << endl; 
-	    cout << endl << " ))))) UPDATE RUN INFO (((((  @iev "<<iev<<" run "<<GetRunInfo()->ID<<" (run-entry "<<irun<<")"<<endl;
+	    if(DBG)cout << endl; 
+	    if(DBG)cout << endl << " ))))) UPDATE RUN INFO (((((  @iev "<<iev<<" run "<<GetRunInfo()->ID<<" (run-entry "<<irun<<")"<<endl;
+	    // ----------------------------------------------------
+	    // update the tracker parameters 
+	    // (non ho trovato nessun altro modo sicuro di farlo...)
+	    // ----------------------------------------------------
+	    if(!dbc || (dbc && !dbc->IsConnected()) )SetDBConnection();
+	    TrkParams::Set(GetRunInfo() ,dbc);
+	    if(dbc)dbc->Close();
 //	    cout << endl; 
 	    prevshift = 0;
 	    return true;
@@ -3892,17 +3914,17 @@
 //     }
 
     if ( TRK0 ){
-	TrkParams::Load(6);
-	if( !TrkParams::IsLoaded(6) ){
-	    cout << " TTree* PamLevel2::GetYodaTree( ) -- WARNING -- VK-mask not loaded"<<endl;
-	};
-	if(!dbc || (dbc && !dbc->IsConnected()))SetDBConnection();
+// 	TrkParams::Load(6);
+// 	if( !TrkParams::IsLoaded(6) ){
+// 	    cout << " TTree* PamLevel2::GetYodaTree( ) -- WARNING -- VK-mask not loaded"<<endl;
+// 	};
+//	if(!dbc || (dbc && !dbc->IsConnected()))SetDBConnection();
 	TrkParams::SetCalib(run_obj,dbc);
 	TrkParams::LoadCalib( );
 	if( !TrkParams::CalibIsLoaded() ){
 	    cout << " TTree* PamLevel2::GetYodaTree( ) -- WARNING -- Calibration not loaded"<<endl;
 	};
-	dbc->Close(); // EMILIANO, do not leave open connections, open only when needed
+// 	if(dbc)dbc->Close(); // EMILIANO, do not leave open connections, open only when needed
     }
 
 
@@ -3971,8 +3993,8 @@
     do{
 	if(shift>0){	 
 	    cout << " PKTNUM  L2 --- "<< pktn << " --- L0 --- "<< GetEventHeader()->GetPscuHeader()->GetCounter()<<endl;
-	    cout << "         RUN: ID "<< GetRunInfo()->ID << " ID_ROOT_L0 "<<run_obj->ID_ROOT_L0<<" ID_RUN_FRAG "<<GetRunInfo()->ID_RUN_FRAG << " EV_FROM "<<GetRunInfo()->EV_FROM  <<endl;
-	    cout << "         L2 <--> L0 mismatch ( irun "<<irun<<" irunentry "<<irunentry<<" shift "<<shift<<" prevshift "<<prevshift<<" )"<<endl;
+	    if(DBG)cout << "         RUN: ID "<< GetRunInfo()->ID << " ID_ROOT_L0 "<<run_obj->ID_ROOT_L0<<" ID_RUN_FRAG "<<GetRunInfo()->ID_RUN_FRAG << " EV_FROM "<<GetRunInfo()->EV_FROM  <<endl;
+	    if(DBG)cout << "         L2 <--> L0 mismatch ( irun "<<irun<<" irunentry "<<irunentry<<" shift "<<shift<<" prevshift "<<prevshift<<" )"<<endl;
 	}
 	answer = l0_tree->GetEntry(quellagiusta+(Long64_t)shift+(Long64_t)prevshift);
 	shift++;
@@ -4020,6 +4042,7 @@
  */
 Bool_t PamLevel2::SetDBConnection(){
 
+//    cout << "PamLevel2::SetDBConnection()" << endl;
     if(DBG){
 	cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<<endl; 
 	cout<<"Connecting to DB"<<endl;