--- PamelaLevel2/src/PamLevel2.cpp	2014/10/15 08:51:51	1.95
+++ PamelaLevel2/src/PamLevel2.cpp	2015/11/02 09:59:54	1.112
@@ -315,17 +315,34 @@
 
 //   candeleteobj = 0;
 
+//  cout << t<<"--"<<c<<"--"<<o<<"--"<<r<<"--"<< endl;
 
   if (t){
       trk_track = new TrkTrack(*t);
       trk_ext_track = new ExtTrack(*t);//NB!! ha dimensione 6 invece che 8
+  }else{
+      trk_track = new TrkTrack();
+      trk_ext_track = new ExtTrack();
+      
   }
+
   if (c)
     calo_track = new CaloTrkVar(*c);
+  else
+    calo_track = new CaloTrkVar();
+  
   if (o)
     tof_track = new ToFTrkVar(*o);
+  else
+    tof_track = new ToFTrkVar();
+
   if (r)
     orb_track = new OrbitalInfoTrkVar(*r);
+  else
+    orb_track = new OrbitalInfoTrkVar();
+
+//  cout << trk_track<<"--"<< calo_track <<"--"<<tof_track<<"--"<<orb_track<<"--"<< endl;
+
   candeleteobj = 1;
 
 }
@@ -335,6 +352,8 @@
  */
 PamTrack::PamTrack(ExtTrack* t, CaloTrkVar* c, ToFTrkVar* o, OrbitalInfoTrkVar *r) {
 
+    
+
   trk_ext_track = 0;
   trk_track = 0;
   calo_track = 0;
@@ -354,34 +373,54 @@
 
 //   candeleteobj = 0;
 
+//  cout << t<<"--"<<c<<"--"<<o<<"--"<<r<<"--"<< endl;
+
 
   if (t){
 ////      trk_track = new TrkTrack(*t);//in this case TrkTrack object remains null
       trk_ext_track = new ExtTrack(*t);
+  }else{
+      trk_ext_track = new ExtTrack();
   }
   if (c)
     calo_track = new CaloTrkVar(*c);
+  else
+      calo_track = new CaloTrkVar();
+
   if (o)
     tof_track = new ToFTrkVar(*o);
+  else
+    tof_track = new ToFTrkVar();
+
   if (r)
     orb_track = new OrbitalInfoTrkVar(*r);
+  else
+    orb_track = new OrbitalInfoTrkVar();
+
+//  cout << trk_track<<"--"<< calo_track <<"--"<<tof_track<<"--"<<orb_track<<"--"<< endl;
+
   candeleteobj = 1;
   pscore = 0;
   iscore = 0;
 
 }
 ;
+
 PamTrack::PamTrack(const PamTrack& track) {
 
   TrkTrack *t = track.trk_track;
+  ExtTrack *e = track.trk_ext_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;
+  trk_ext_track = 0;
+  trk_track     = 0;
+  calo_track    = 0;
+  tof_track     = 0;
+  orb_track     = 0;
+  if(e)
+      trk_ext_track = new ExtTrack(*e);
   if (t)
     trk_track = new TrkTrack(*t);
   if (c)
@@ -390,18 +429,36 @@
     tof_track = new ToFTrkVar(*o);
   if (r)
     orb_track = new OrbitalInfoTrkVar(*r);
+
   candeleteobj = 1;
   pscore = 0;
   iscore = 0;
 
 }
-void PamTrack::Clear() {
 
-  //    cout << "PamTrack::Clear() "<<candeleteobj<<endl;
+void PamTrack::Copy( PamTrack& track) const {
+
+    track.trk_track = trk_track;
+    track.trk_ext_track = trk_ext_track;
+    track.calo_track = calo_track;
+    track.tof_track = tof_track;
+    track.orb_track = orb_track;
+    
+    track.candeleteobj = candeleteobj;
+    track.pscore = pscore;
+    track.iscore = iscore;
+    
+}
+
+
+
+void PamTrack::Clear(Option_t *option) {
+
+//    cout << "PamTrack::Clear( "<<option<<" ) "<<candeleteobj<<endl;
   if (candeleteobj) {
       
     if (trk_ext_track)
-      trk_ext_track->ExtTrack::Clear();
+      trk_ext_track->ExtTrack::Clear(option);
     if (trk_track)
       trk_track->TrkTrack::Clear();
     if (calo_track)
@@ -426,7 +483,8 @@
   //    cout << "PamTrack::Delete() "<<candeleteobj<<endl;
   if (candeleteobj) {
     if (trk_ext_track) {
-      trk_ext_track->ExtTrack::Clear();
+//	trk_ext_track->ExtTrack::Clear("C");//Clear is called for all the array elements 
+	trk_ext_track->ExtTrack::Clear("C+C");//Clear is called for all the array elements passing option 'C'
       delete trk_ext_track;
     }
     if (trk_track) {
@@ -561,6 +619,7 @@
 
   run_obj = 0;//new GL_RUN();
   soft_obj = 0;// Emiliano
+  proc_obj = 0;// Emiliano
   irun = -1LL;
   irunt = -1LL;
   totrunentry = 0LL;
@@ -571,6 +630,7 @@
   gltsync = 0; // Emiliano
   fUpdateRunInfo = true; // Emiliano
   fUseDBinRunInfo = true; // Emiliano
+  fDiscarded = false; //EM
   isSync = false; // by default assume that the level2 file(s) is(are) not sinchronized between L0/DB and L2, that is we miss some packets in L2 due to nested/DV-skipped events
   il0entry = 0LL;
   //  hasL0EE = true;
@@ -586,6 +646,10 @@
 
   run_tree = NULL;
   run_tree_clone = NULL;
+
+  proc_tree = NULL;
+  proc_tree_clone = NULL;
+
   sel_tree = NULL;
   sel_tree_clone = NULL;
 
@@ -617,6 +681,8 @@
   if (strcmp(pamdbpsw, ""))
     psw = pamdbpsw;
 
+  customString = "";
+
   //    sorted_tracks = 0;//new TRefArray();
 
   CAL0 = false;
@@ -633,11 +699,12 @@
   ND = true;
   AC = true;
   ORB = true;
+  PROC = true;
   GP = false;
 
   EXT = false;
   NUC = false;
-  trkAlg = "";//default tracking algorythm
+  trkAlg = "STD";//default tracking algorythm
 
   RUN = true;
 
@@ -675,6 +742,8 @@
     delete run_obj;
   if (soft_obj)
     delete soft_obj; //Emiliano
+  if (proc_obj)
+    delete proc_obj; //Emiliano
 
   //    cout << "void PamLevel2::Clear()"<<endl;
   if (h0_obj)
@@ -875,6 +944,8 @@
     orb2_obj->Clear();
   if (gp_obj)
     gp_obj->Clear();
+  if (proc_obj)
+    proc_obj->Clear();
 
   //    if(sorted_tracks)sorted_tracks->Clear();
   //    sorted_tracks.Clear();
@@ -955,6 +1026,10 @@
   if (sel_tree)
     sel_tree->Delete();;
   sel_tree = NULL;
+
+  if (proc_tree)
+    proc_tree->Delete();
+  proc_tree = NULL;
   //
   // Close file
   //
@@ -978,6 +1053,7 @@
   ac_obj = 0;
   orb2_obj = 0;
   gp_obj = 0;
+  proc_obj = 0;
 
   trk_ext_obj     = 0;
   trk_ext_nuc_obj = 0;
@@ -1009,6 +1085,7 @@
   //
   run_obj = 0;//new GL_RUN();
   soft_obj = 0;// Emiliano
+  proc_obj = 0;// Emiliano
   irun = -1;
   irunt = -1;
   totrunentry = 0LL;
@@ -1183,6 +1260,12 @@
   if (!objname.CompareTo("SoftInfo"))
     return &soft_obj; // Emiliano
 
+  if (!objname.CompareTo("ProcInfo")){
+    if (!proc_obj) 
+      proc_obj = new ProcInfo();    
+    return &proc_obj; // Emiliano
+  }
+
   return NULL;
 }
 ;
@@ -1200,9 +1283,11 @@
     return 0;
 
   if (calo2_obj->CaloLevel2::ntrk() == 0) {
-    cout << "PamLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo " << seqno
-        << " but no Calorimeter tracks are stored" << endl;
-    return NULL;
+      if( seqno >=0 ){
+	  cout << "PamLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo " << seqno;
+	  cout << " but no Calorimeter tracks are stored" << endl;
+      }
+      return NULL;
   };
 
   CaloTrkVar *c = 0;
@@ -1215,7 +1300,7 @@
 
   if (!c || seqno != c->trkseqno) {
     c = 0;
-    if (seqno != -1)
+    if (seqno != -1 && seqno>=0) 
       cout << "PamLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo " << seqno
           << " does not match Calorimeter stored tracks" << endl;
   };
@@ -1334,65 +1419,65 @@
 //
 //
 //--------------------------------------
-/**
- * Give the pamela track associated to a tracker track, retrieving related calorimeter, orbitalinfo and tof track information.
- */
-PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t) {
-
-  cout << "PamLevel2::GetPamTrackAlong(TrkTrack* t) **obsolete** " << endl;
-  cout << "(if you use it, remember to delete the PamTrack object)" << endl;
-
-  CaloTrkVar *c = 0;
-  ToFTrkVar *o = 0;
-  OrbitalInfoTrkVar *r = 0;
+// /**
+//  * Give the pamela track associated to a tracker track, retrieving related calorimeter, orbitalinfo and tof track information.
+//  */
+// PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t) {
+
+//   cout << "PamLevel2::GetPamTrackAlong(TrkTrack* t) **obsolete** " << endl;
+//   cout << "(if you use it, remember to delete the PamTrack object)" << endl;
+
+//   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 (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, r);
 
-  //    if(t && c && o)track = new PamTrack(t,c,o);
-  PamTrack *track = new PamTrack(t, c, o, r);
+//   return track;
 
-  return track;
-
-}
-;
+// }
+// ;
 //--------------------------------------
 //
 //
 //--------------------------------------
-/**
- * Retrieves the it-th stored track.
- * It override TrkLevel2::GetTrack(int it).
- * @param itrk Track number, ranging from 0 to GetNTracks().
- */
+// /**
+//  * Retrieves the it-th stored track.
+//  * It override TrkLevel2::GetTrack(int it).
+//  * @param itrk Track number, ranging from 0 to GetNTracks().
+//  */
+
+// 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) OrbitalInfo::GetOrbitalInfoTrkVar(Int_t notrack)"
+//       << endl;
+//   cout << "(if you use it, remember to delete the PamTrack object)" << endl;
+//   PamTrack *track = 0;
 
-PamTrack* PamLevel2::GetStoredTrack(Int_t itrk) {
+//   if (itrk >= 0 && itrk < trk2_obj->TrkLevel2::ntrk()) {
 
-  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) OrbitalInfo::GetOrbitalInfoTrkVar(Int_t notrack)"
-      << endl;
-  cout << "(if you use it, remember to delete the PamTrack object)" << endl;
-  PamTrack *track = 0;
+//     TrkTrack *t = trk2_obj->TrkLevel2::GetStoredTrack(itrk);
+//     track = GetPamTrackAlong(t);
 
-  if (itrk >= 0 && itrk < trk2_obj->TrkLevel2::ntrk()) {
-
-    TrkTrack *t = trk2_obj->TrkLevel2::GetStoredTrack(itrk);
-    track = GetPamTrackAlong(t);
-
-  }
-  else {
-    cout << "PamLevel2::GetStoredTrack(int) : tracker track SeqNo " << itrk << " does not exist (GetNTracks() = "
-        << trk2_obj->TrkLevel2::GetNTracks() << ")" << endl;
-  };
+//   }
+//   else {
+//     cout << "PamLevel2::GetStoredTrack(int) : tracker track SeqNo " << itrk << " does not exist (GetNTracks() = "
+//         << trk2_obj->TrkLevel2::GetNTracks() << ")" << endl;
+//   };
 
-  return track;
+//   return track;
 
-}
+// }
 //--------------------------------------
 //
 
@@ -1446,14 +1531,20 @@
   Int_t ObjectNumber = TProcessID::GetObjectCount();
 
   // create TCloneArrays to store tracks and its images
-  if (!tsorted)
-    tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
-  tsorted->Delete();
-  TClonesArray &ttsorted = *tsorted;
+//   if (!tsorted)
+//     tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
+//   tsorted->Clear("C+C");//Delete();
+//   if (!timage)
+//     timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
+//   timage->Clear("C+C");//Delete();
+
+  if(tsorted)delete tsorted;
+  if(timage) delete timage;
+  tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
+  timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
 
-  if (!timage)
-    timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
-  timage->Delete();
+
+  TClonesArray &ttsorted = *tsorted;
   TClonesArray &ttimage = *timage;
 
 
@@ -1895,6 +1986,9 @@
     //	cout<<"o "<<cp<<endl;
     //	cout<<"o "<<op<<endl;
 
+//     cout <<"old p:"<< tp<<"--"<<cp<<"--"<<op<<"--"<<rp<<"--"<< endl;
+//     cout <<"old i:"<< ti<<"--"<<ci<<"--"<<oi<<"--"<<ri<<"--"<< endl;
+
     new (ttsorted[i]) PamTrack(tp, cp, op, rp);
     new (ttimage[i]) PamTrack(ti, ci, oi, ri);
 
@@ -1922,6 +2016,8 @@
   issorted = true;
   lastsorted = GetReadEntry();
 
+//  cout <<" SortTracks() -- end"<<endl;
+
 }
 ;
 //
@@ -1989,36 +2085,55 @@
   //-----------------------------------------------------------
   // create/reset TCloneArrays to store tracks and their images
   //-----------------------------------------------------------
+
+//  cout << " PamLevel2::SortTracksNew() --- Clear TClonesArray objects"<<endl;
   
   // main tracks from standard alg
-  if (!tsorted)
-      tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
-  tsorted->Delete();
-  // track images from standard alg
-  if (!timage)
-      timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
-  timage->Delete();
-  // tracks from extended algorythm
-  if(EXT && !text)
-      text = new TClonesArray("PamTrack",trk_ext_obj->GetEntries()); 
-  if(text)text->Delete();
+//   if (!tsorted)
+//       tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
+//   tsorted->Clear("C+C");//Delete();
+//   // track images from standard alg
+//   if (!timage)
+//       timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
+//   timage->Clear("C+C");//Delete();
+//   // tracks from extended algorythm
+//   if(EXT && !text)
+//       text = new TClonesArray("PamTrack",trk_ext_obj->GetEntries()); 
+//   if(text)text->Clear("C+C");//Delete();
+
+  if(tsorted)delete tsorted;
+  if(timage) delete timage;
+  if(text)   delete text;
+  tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
+  timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
+  text = new TClonesArray("PamTrack",trk_ext_obj->GetEntries()); 
 
   //-----------------------------------------------------------
   // create/reset TCloneArrays to store tracks and their images
   //-----------------------------------------------------------
   if(NUC){
+
+      
+      if(tsorted_nuc)delete tsorted_nuc;
+      if(timage_nuc) delete timage_nuc;
+      if(text_nuc)   delete text_nuc;
+      tsorted_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
+      timage_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
+      text_nuc = new TClonesArray("PamTrack",trk_ext_nuc_obj->GetEntries()); 
+
   // main tracks from standard alg
-      if (!tsorted_nuc)
-	  tsorted_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
-      tsorted_nuc->Delete();
-      // track images from standard alg
-      if (!timage_nuc)
-	  timage_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
-      timage_nuc->Delete();
-      // tracks from extended algorythm
-      if(EXT && !text_nuc)
-	  text_nuc = new TClonesArray("PamTrack",trk_ext_nuc_obj->GetEntries()); 
-      if(text_nuc)text_nuc->Delete();
+//       if (!tsorted_nuc)
+// 	  tsorted_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
+//       tsorted_nuc->Clear("C+C");//Delete();
+//       // track images from standard alg
+//       if (!timage_nuc)
+// 	  timage_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
+//       timage_nuc->Clear("C+C");//Delete();
+//       // tracks from extended algorythm
+//       if(EXT && !text_nuc)
+// 	  text_nuc = new TClonesArray("PamTrack",trk_ext_nuc_obj->GetEntries()); 
+//       if(text_nuc)text_nuc->Clear("C+C");//Delete();
+
   }
   //--------------------------------------------------
   // retrieve sorting method
@@ -2066,34 +2181,30 @@
 
       if(doit == 0){
 
-  	  trk2  = trk2_obj;
-	  calo2 = calo2_obj;
-	  tof2  = tof2_obj;
-	  orb2  = orb2_obj;
-      
-	  trkext = trk_ext_obj;
-	  caloext = calo_ext_obj;
-	  tofext = tof_ext_obj;
-	  orbext = orb_ext_obj;
-    
-      
-      
+  	  trk2  = (TRK2 ? trk2_obj: NULL);;
+	  calo2 = (CAL2 ? calo2_obj: NULL);;
+	  tof2  = (TOF  ?tof2_obj: NULL);;
+	  orb2  = (ORB ? orb2_obj: NULL);;
       
+	  trkext = (TRK2 ? trk_ext_obj: NULL);;
+	  caloext = (CAL2 ? calo_ext_obj: NULL);;
+	  tofext = (TOF  ?tof_ext_obj: NULL);;
+	  orbext = (ORB ? orb_ext_obj: NULL);;
 
       }else if (doit == 1){
 
 	  if(!NUC)break;
 
       
-	  trk2  = trk2_nuc_obj;
-	  calo2 = calo2_nuc_obj;
-	  tof2  = tof2_nuc_obj;
-	  orb2  = orb2_nuc_obj;
+	  trk2  = (TRK2 ?trk2_nuc_obj: NULL);;
+	  calo2 = (CAL2 ? calo2_nuc_obj: NULL);;
+	  tof2  = (TOF  ?tof2_nuc_obj: NULL);;
+	  orb2  = (ORB ? orb2_nuc_obj: NULL);;
       
-	  trkext = trk_ext_nuc_obj;
-	  caloext = calo_ext_nuc_obj;
-	  tofext = tof_ext_nuc_obj;
-	  orbext = orb_ext_nuc_obj;
+	  trkext = (TRK2 ?trk_ext_nuc_obj: NULL);;
+	  caloext = (CAL2 ? calo_ext_nuc_obj: NULL);;
+	  tofext = (TOF  ?tof_ext_nuc_obj: NULL);;
+	  orbext = (ORB ? orb_ext_nuc_obj: NULL);;
       
       
       
@@ -2126,10 +2237,12 @@
 	  OrbitalInfoTrkVar *rs = 0;
 
 	  // get tracker tracks
-	  TrkTrack *tp = trk2->GetTrack(i); //tracker
-	  CaloTrkVar *cp = calo2->GetCaloStoredTrack(tp->GetSeqNo());
-	  ToFTrkVar *op = tof2->GetToFStoredTrack(tp->GetSeqNo());
-	  OrbitalInfoTrkVar *rp = orb2->GetOrbitalInfoStoredTrack(tp->GetSeqNo());
+	  TrkTrack *tp          = (TRK2 ? trk2->GetTrack(i): NULL); //tracker
+	  CaloTrkVar *cp        = (CAL2 ? calo2->GetCaloStoredTrack(tp->GetSeqNo()) : NULL);
+	  ToFTrkVar *op         = (TOF  ? tof2->GetToFStoredTrack(tp->GetSeqNo())   : NULL);
+	  OrbitalInfoTrkVar *rp = (ORB ? orb2->GetOrbitalInfoStoredTrack(tp->GetSeqNo()) : NULL);
+
+//	  cout << "ORB="<<ORB<<" rp="<<rp<<endl;
 
 	  TrkTrack *ti = 0; //tracker (image)
 	  CaloTrkVar *ci = 0;
@@ -2145,10 +2258,10 @@
 
 	  if (tp->HasImage()) {
 
-	      ti = trk2->GetTrackImage(i); //tracker (image)
-	      ci = calo2->GetCaloStoredTrack(ti->GetSeqNo());
-	      oi = tof2->GetToFStoredTrack(ti->GetSeqNo());
-	      ri = orb2->GetOrbitalInfoStoredTrack(ti->GetSeqNo());
+	      ti = (TRK2 ? trk2->GetTrackImage(i) : NULL); //tracker (image)
+	      ci = (CAL2 ? calo2->GetCaloStoredTrack(ti->GetSeqNo()): NULL);
+	      oi = (TOF  ? tof2->GetToFStoredTrack(ti->GetSeqNo()): NULL);
+	      ri = (ORB ? orb2->GetOrbitalInfoStoredTrack(ti->GetSeqNo()): NULL);
 
 	      //	    cout << "its image "<<i << " "<<hex<< ti <<dec<< endl;
 
@@ -2516,6 +2629,11 @@
 	  //	cout<<"o "<<cp<<endl;
 	  //	cout<<"o "<<op<<endl;
 
+// 	  cout <<"p:"<< tp<<"--"<<cp<<"--"<<op<<"--"<<rp<<"--"<< endl;
+// 	  cout <<"i:"<< ti<<"--"<<ci<<"--"<<oi<<"--"<<ri<<"--"<< endl;
+
+
+
 	  new (ttsorted[i]) PamTrack(tp, cp, op, rp);
 	  new (ttimage[i]) PamTrack(ti, ci, oi, ri);
 
@@ -2539,7 +2657,7 @@
 //   }
 
 
-//       cout << "text "<<doit<<endl;
+      //      cout << "text "<<orbext<<endl;
 
 
       //--------------------------------------------------
@@ -2548,11 +2666,30 @@
       if(EXT){
 	  for(int it=0; it<trkext->GetEntries(); it++){
 	      
-	      new (ttext[it]) PamTrack((ExtTrack*)(*trkext)[it], (CaloTrkVar*)(*caloext)[it], (ToFTrkVar*)(*tofext)[it], (OrbitalInfoTrkVar*)(*orbext)[it]);
+//	      if( 
+// 		  it < caloext->GetEntries() &&
+// 		  it < tofext->GetEntries() &&
+// 		  it < orbext->GetEntries() &&
+//		  true){
+	      
+	      ExtTrack *t = (trkext ?(ExtTrack*)(*trkext)[it]:NULL ); 
+	      CaloTrkVar *c =(caloext ?(CaloTrkVar*)(*caloext)[it]:NULL ); 
+	      ToFTrkVar *o = (tofext ?(ToFTrkVar*)(*tofext)[it]:NULL ); 
+	      OrbitalInfoTrkVar *r =(orbext ?(OrbitalInfoTrkVar*)(*orbext)[it]:NULL ); 
+//	      cout <<"ext:"<< t<<"--"<<c<<"--"<<o<<"--"<<r<<"--"<< endl;
 	      
+	      new (ttext[it]) PamTrack(t, c, o, r);
+// 	      }else{
+// 		  cout << " PamLevel2::SortTracksNew() --> ORRORE E RACCAPRICCIO!!!"<<endl;
+// 		  cout << " trk  ext-tracks = "<<trkext->GetEntries()<<endl;
+// 		  cout << " calo ext-tracks = "<<caloext->GetEntries()<<endl;
+// 		  cout << " tof  ext-tracks = "<<tofext->GetEntries()<<endl;
+// 		  cout << " orb  ext-tracks = "<<orbext->GetEntries()<<endl;
+// 	      }
 	  }
       }
       
+//  cout <<" SortTracksNew() -- end"<<endl;
 
 
   };
@@ -2580,7 +2717,7 @@
 /**
  * This method overrides TrkLevel2::GetTracks(), where sorting is done by decreasing number of fit points and increasing chi^2.
  * PamLevel2::GetTracks() keeps the same track order given by TrkLevel2::GetTracks(), but checks image selection by using calorimeter and ToF tracking information.
- /*
+ */
 
 // TRefArray *PamLevel2::GetTracks(){
 
@@ -2598,12 +2735,13 @@
 //
 //
 //--------------------------------------
+
 /**
  * Retrieves the it-th Pamela "physical" track.
  * It override TrkLevel2::GetTrack(int it).
  * @param it Track number, ranging from 0 to GetNTracks().
  */
-PamTrack *PamLevel2::GetTrack(int it) {
+PamTrack *PamLevel2::GetTrackOld(int it) {
 
   PamTrack *track = NULL;
 
@@ -2612,15 +2750,22 @@
   //  *-*-*-*-*-*-*-*-*-*-*-*-*
   if (!tsorted)
     return track;
-  if (!tsorted->GetEntries())
+//   if (!tsorted->GetEntries())
+//     return track;
+  if (tsorted->GetEntries()==0)
     return track;
 
 
-  if (it >= 0 && it < trk2_obj->TrkLevel2::GetNTracks()) {
-    track = (PamTrack*)((*tsorted)[it]);
+  if (
+      tsorted && 
+      it >= 0 && 
+      it < trk2_obj->TrkLevel2::GetNTracks() && 
+      it < tsorted->GetEntries() && 
+      true) {
+      track = (PamTrack*)((*tsorted)[it]);
   }
   else {
-    cout << "PamLevel2::GetTrack(int) : tracker track SeqNo " << it << " does not exist (GetNTracks() = "
+    cout << "PamLevel2::GetTrackOld(int) : tracker track SeqNo " << it << " does not exist (GetNTracks() = "
         << trk2_obj->TrkLevel2::GetNTracks() << ")" << endl;
   };
 
@@ -2636,10 +2781,10 @@
  * @param it Track number, ranging from 0 to GetNTracks().
  * @param alg Algorythm, see SetTrakingAlgorythm(char *alg) for explanation.
  */
-PamTrack *PamLevel2::GetTrack(int it, char* alg) {
+PamTrack *PamLevel2::GetTrack(int it, const char* alg) {
 
     TString s(alg);
-    if(!s.CompareTo("") ||!s.CompareTo("STD") )return GetTrack(it); //old algorythm
+    if(!s.CompareTo("") ||!s.CompareTo("STD") )return GetTrackOld(it); //old algorythm
 
 
     SortTracksNew();
@@ -2687,6 +2832,7 @@
     };
     
     
+    cout << "PamTrack *PamLevel2::GetTrack("<<it<<","<<alg<<") -- wrong track number or unrecognised algorithm"<<endl;
     
     return NULL;
 
@@ -2702,7 +2848,7 @@
 
 }
 ;
-Int_t  PamLevel2::GetNTracks(char* alg) {
+Int_t  PamLevel2::GetNTracks(const char* alg) {
 
 
 //    cout << " trk_nuc_obj->GetEntries() "<<trk_nuc_obj->GetEntries()<<" trk2_nuc_obj->GetNTracks() "<<trk2_nuc_obj->GetNTracks()<<endl;
@@ -2723,6 +2869,8 @@
     if(s.Contains("NUC", TString::kIgnoreCase) && NUC )
 	return trk2_nuc_obj->TrkLevel2::GetNTracks(); 
     
+    cout << "Int_t  PamLevel2::GetNTracks("<<alg<<") -- unrecognised algorithm"<<endl;
+
     return 0;
 
 }
@@ -2736,7 +2884,7 @@
  * Retrieves (if present) the image of the it-th Pamela "physical" track, sorted by the method PamLevel2::SortTracks().
  * @param it Track number, ranging from 0 to GetNTracks().
  */
-PamTrack *PamLevel2::GetTrackImage(int it) {
+PamTrack *PamLevel2::GetTrackImageOld(int it) {
 
   //  *-*-*-*-*-*-*-*-*-*-*-*-*
   SortTracks();
@@ -2760,7 +2908,7 @@
     };
   }
   else {
-    cout << "PamLevel2::GetTrackImage(int) : Tracker track SeqNo " << it << " does not exist (GetNTracks() = "
+    cout << "PamLevel2::GetTrackImageOld(int) : Tracker track SeqNo " << it << " does not exist (GetNTracks() = "
         << trk2_obj->TrkLevel2::GetNTracks() << ")" << endl;
   };
 
@@ -2771,10 +2919,10 @@
  * @param it Track number, ranging from 0 to GetNTracks().
  * @param alg Algorythm, see SetTrakingAlgorythm(char *alg) for explanation.
  */
-PamTrack *PamLevel2::GetTrackImage(int it, char* alg) {
+PamTrack *PamLevel2::GetTrackImage(int it, const char* alg) {
 
     TString s(alg);
-    if(!s.CompareTo("") || !s.CompareTo("STD"))return GetTrack(it); //old algorythm
+    if(!s.CompareTo("") || !s.CompareTo("STD"))return GetTrackImageOld(it); //old algorythm
 
 
     SortTracksNew();
@@ -2792,6 +2940,8 @@
 		PamTrack *temp = (PamTrack*) t[it];
 		if (temp->GetTrkTrack()->HasImage()) {
 		    return (PamTrack*)((*timage_nuc)[it]);                 //ok return the track
+		}else{
+		    return NULL;
 		}
 		
 	    }
@@ -2805,11 +2955,15 @@
 		PamTrack *temp = (PamTrack*) t[it];
 		if (temp->GetTrkTrack()->HasImage()) {		    
 		    return (PamTrack*)((*timage)[it]);                 //ok return the track
+		}else{
+		    return NULL;
 		}
 	    }
 	}
 
     } 
+
+//    cout << "PamTrack *PamLevel2::GetTrackImage("<<it<<","<<alg<<") -- wrong track number or unrecognised algorithm"<<endl;
       
     return NULL;
 
@@ -3269,6 +3423,7 @@
   TChain *G = 0;
 
   TChain *L = 0;
+  TChain *P = 0;
 
   if (TRK2 || TRK1 || TRKh)
     T = new TChain("Tracker");
@@ -3288,6 +3443,8 @@
     B = new TChain("OrbitalInfo");
   if (GP)
     G = new TChain("h20");
+  if (PROC)
+    P = new TChain("ProcessingInfo");
   L = new TChain("SelectionList");
 
   // loop over files and create chains
@@ -3315,6 +3472,8 @@
         B->Add(name);
       if (GP)
         G->Add(name);
+      if (P)
+        P->Add(name);
       if (SELLI == 1)
         L->Add(name);
      };
@@ -3440,12 +3599,25 @@
       L->Delete();
   };
 
+  //ProcessingInfo EM
+  if ( P && P->GetEntries() ){
+    cout << "----------------------------------------------------" << endl;
+    cout << ">>> Found ProcessingInfo <<<" << endl;
+    //	L->SetBranchAddress("RunEntry",&irun);
+    P->SetBranchAddress("ProcInfo", &proc_obj);//NEWNEW
+    proc_tree = P;
+  } else {
+    //    proc_tree = new TChain("ProcessingInfo","Log of data processing");
+    //    proc_tree->Branch("ProcInfo", "ProcInfo", &proc_obj);
+    cout << "----------------------------------------------------" << endl;
+    cout << ">>> ProcessingInfo not found, PROC set to false and continue (not 10RED files?)<<<" << endl;
+    PROC = false;
+  }
   //  --------------------------------------------
   //  return the pamela chain with all the friends
   //  --------------------------------------------
 
   pam_tree = Trout;
-
   return Trout;
 }
 
@@ -3767,36 +3939,36 @@
  * "NUCEXT" --> as "EXT", but for nuclei 
  * "NUCEXTF" --> as "EXTF", but for nuclei
  */
-void PamLevel2::SetTrakingAlgorythm(char *alg){
+// void PamLevel2::SetTrackingAlgorythm(const char *alg){
  
 
-    TString s(alg);
-    if(s.Contains("NUC", TString::kIgnoreCase) && !NUC)
-	cout << "Warning! NUC algorythm requested, but branches are missing"<<endl;
-    if(s.Contains("EXT", TString::kIgnoreCase) && !EXT)
-	cout << "Warning! EXT algorythm requested, but branches are missing"<<endl;;
+//     TString s(alg);
+//     if(s.Contains("NUC", TString::kIgnoreCase) && !NUC)
+// 	cout << "Warning! NUC algorythm requested, but branches are missing"<<endl;
+//     if(s.Contains("EXT", TString::kIgnoreCase) && !EXT)
+// 	cout << "Warning! EXT algorythm requested, but branches are missing"<<endl;;
 	    
-    trkAlg = alg;
+//     trkAlg = alg;
 
-};
-char* PamLevel2::GetTrakingAlgorythm(){
+// };
+// const char* PamLevel2::GetTrackingAlgorythm(){
 
 
-    cout<<endl<<" Implemented tracking algorythm: ";
-    cout<<endl<<"  \"\" or \"STD\"  --> take the output of the standard tracking algorythm";
-    cout<<endl<<"  \"NUC\"      --> take the output of the standard tracking algorythm for nuclei cluster selection";
-    cout<<endl<<"  \"EXT\"      --> in case the standard tracking algorythm has not found any track,";
-    cout<<endl<<"                 take the output of the extended one";
-    cout<<endl<<"  \"EXTF\"     --> force the extended tracking algorythm";
-    cout<<endl<<"  \"NUCEXT\"   --> as \"EXT\", but for nuclei ";
-    cout<<endl<<"  \"NUCEXTF\"  --> as \"EXTF\", but for nuclei";
-
-    cout<<endl;
-    cout<<" <<Currently set algorythm>> "<<trkAlg<<endl;
-    cout<<endl;
+//     cout<<endl<<" Implemented tracking algorythm: ";
+//     cout<<endl<<"  \"\" or \"STD\"  --> take the output of the standard tracking algorythm";
+//     cout<<endl<<"  \"NUC\"      --> take the output of the standard tracking algorythm for nuclei cluster selection";
+//     cout<<endl<<"  \"EXT\"      --> in case the standard tracking algorythm has not found any track,";
+//     cout<<endl<<"                 take the output of the extended one";
+//     cout<<endl<<"  \"EXTF\"     --> force the extended tracking algorythm";
+//     cout<<endl<<"  \"NUCEXT\"   --> as \"EXT\", but for nuclei ";
+//     cout<<endl<<"  \"NUCEXTF\"  --> as \"EXTF\", but for nuclei";
+
+//     cout<<endl;
+//     cout<<" <<Currently set algorythm>> "<<trkAlg<<endl;
+//     cout<<endl;
 
-    return trkAlg;
-};
+//     return trkAlg;
+// };
 
 
 
@@ -3955,7 +4127,7 @@
     isFragment = false;
     run_tree->GetEntry(irun);
     if (!GetOrbitalInfo())
-      cout << "** WARNING ** missing OrbitalInfo ---> run info might be not correctly updated " << endl;
+	cout << "PamLevel2::UpdateRunInfo(Long64_t "<<iev<<") ** WARNING ** missing OrbitalInfo    ORB="<<ORB << endl;
     if ( fUseDBinRunInfo ){
       if (gltsync)
         delete gltsync; //Emiliano
@@ -4537,6 +4709,7 @@
   //    cout << "Checking file: "<<f->GetName()<<endl;
   if (!f || f->IsZombie()) {
     cout << "File: " << f->GetName() << " Non valid root file" << endl;
+    fDiscarded = true;
     return;
   }
 
@@ -4743,6 +4916,7 @@
   TFile *f = new TFile(name.Data());
   if (!f || f->IsZombie()) {
     cout << "File: " << f->GetName() << " discarded ---- Non valid root file" << endl;
+    fDiscarded = true;
     return false;
   }
   //    cout << "Get list of keys: "<<f<<endl;
@@ -4781,6 +4955,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- SelectionList tree has " << nevt
               << " events instead of " << nev << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4795,6 +4970,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- Trigger tree has " << nevt << " events instead of "
               << nev << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4808,6 +4984,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- ToF tree has " << nevt << " events instead of " << nev
               << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4821,6 +4998,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- S4 tree has " << nevt << " events instead of " << nev
               << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4835,6 +5013,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- NeutronD tree has " << nevt << " events instead of "
               << nev << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4848,6 +5027,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- Anticounter tree has " << nevt << " events instead of "
               << nev << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4861,6 +5041,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- OrbitalInfo tree has " << nevt << " events instead of "
               << nev << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4874,6 +5055,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- Tracker tree has " << nevt << " events instead of "
               << nev << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4897,6 +5079,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- Calorimeter tree has " << nevt << " events instead of "
               << nev << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4919,6 +5102,7 @@
         if (nev && nevt != nev) {
           cout << "File: " << f->GetName() << " discarded ---- h20 tree has " << nevt << " events instead of " << nev
               << endl;
+          fDiscarded = true;
           return false;
         }
         nev = nevt;
@@ -4931,10 +5115,12 @@
     SELLI = (Int_t) SELLI__ok;
   if (SELLI == 0 && SELLI__ok) {
     cout << "File: " << f->GetName() << " discarded ---- found SelectionList (it is not a full-event file)" << endl;
+    fDiscarded = true;
     return false;
   }
   if (SELLI == 1 && !SELLI__ok) {
     cout << "File: " << f->GetName() << " discarded ---- SelectionList missing" << endl;
+    fDiscarded = true;
     return false;
   }
 
@@ -4973,51 +5159,63 @@
 
   if (CAL1 && !CAL1__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing CaloLevel1 branch" << endl;
+    fDiscarded = true;
     return false;
   };
   if (CAL2 && !CAL2__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing CaloLevel2 branch" << endl;
+    fDiscarded = true;
     return false;
   };
   if (TRK2 && !TRK2__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing TrkLevel2 branch" << endl;
+    fDiscarded = true;
     return false;
   };
   if (TRK1 && !TRK1__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing TrkLevel1 branch" << endl;
+    fDiscarded = true;
     return false;
   };
   if (TRKh && !TRKh__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing TrkHough branch" << endl;
+    fDiscarded = true;
     return false;
   };
   if (ORB && !ORB__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing ORB tree" << endl;
+    fDiscarded = true;
     return false;
   };
   if (AC && !AC__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing AC tree" << endl;
+    fDiscarded = true;
     return false;
   };
   if (S4 && !S4__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing S4 tree" << endl;
+    fDiscarded = true;
     return false;
   };
   if (TOF && !TOF__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing ToF tree" << endl;
+    fDiscarded = true;
     return false;
   };
 
   if (ND && !ND__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing ND tree" << endl;
+    fDiscarded = true;
     return false;
   };
   if (TRG && !TRG__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing Trigger tree" << endl;
+    fDiscarded = true;
     return false;
   };
   if (GP && !GP__ok) {
     cout << "File: " << f->GetName() << " discarded ---- Missing h20 tree" << endl;
+    fDiscarded = true;
     return false;
   };
 
@@ -5117,8 +5315,39 @@
     //	cout << i<< "\t | "<<GetRunInfo()->ID<<"\t "<<GetRunInfo()->NEVENTS<< "\t "<<GetRunInfo()->RUNHEADER_TIME<<" <---> "<<GetRunInfo()->RUNTRAILER_TIME<<endl;
     run_tree_clone->Fill();
   }
-  cout << "----------------------------------------------------" << endl;
-
+  //
+  // replicate processinginfo tree
+  //
+  if ( PROC ){ // EMEMEM
+    proc_tree_clone = new TTree("ProcessingInfo","Log of data processing");
+    proc_tree_clone->Branch("ProcInfo", "ProcInfo", GetPointerTo("ProcInfo"));
+    cout << "ProcessingInfo: branch ProcessingInfo" << endl;  
+    // ------------------
+    // replicate processinginfo tree
+    // ------------------
+    //    cout << "----------------------------------------------------"<<endl;
+    //    cout << "irun\t | RUN\t NEVENTS\t absolute time"<<endl;
+    for (Int_t i = 0; i < proc_tree->GetEntries(); i++) {
+      proc_tree->GetEntry(i);
+      //      cout << i<< "\t | "<<endl;
+      proc_tree_clone->Fill();
+    }
+    if ( SELLI != 2 ){      
+      proc_obj->runID = 0;
+      TTimeStamp *dt = new TTimeStamp();
+      proc_obj->date = dt->AsString();
+      delete dt;
+      proc_obj->commandLine = Form("PamelaLevel2 was called: CAL2 %i CAL1 %i CAL0 %i TRK2 %i TRK1 %i TRKh %i TRK0 %i TOF %i TOF0 %i TRG %i \n                                            S4 %i ND %i AC %i ORB %i GP %i EXT %i NUC %i RUN %i ISGP %i SELLI %i \n Custom string   = %s",CAL2,CAL1,CAL0,TRK2,TRK1,TRKh,TRK0,TOF,TOF0,TRG,S4,ND,AC,ORB,GP,EXT,NUC,RUN,ISGP,SELLI,customString.Data());
+      proc_obj->outputFilename = ofile->GetName();
+      proc_obj->localDir = gSystem->WorkingDirectory();
+      proc_obj->uname = gSystem->GetFromPipe("uname -a");
+      if (!dbc || (dbc && !dbc->IsConnected())) SetDBConnection(); 
+      proc_obj->DB = Form("mysql://%s/%s",dbc->GetHost(),dbc->GetDB());
+      dbc->Close();
+      proc_tree_clone->Fill();
+    }
+    cout << "----------------------------------------------------" << endl;
+  }
   // ------------------------------------
   // add branch with dead and live times
   // ------------------------------------
@@ -5309,6 +5538,11 @@
     if (!name.CompareTo(na))
       return sel_tree_clone;
   }
+  if (proc_tree_clone && PROC) {
+    TString na = proc_tree_clone->GetName();
+    if (!name.CompareTo(na))
+      return proc_tree_clone;
+  }
   return NULL;
 
 }
@@ -5328,9 +5562,13 @@
   for (Int_t i = 0; i < NCLONES; i++) {
     if (pam_tree_clone[i]) {
       cout << pam_tree_clone[i]->GetName() << endl;
-      pam_tree_clone[i]->Write();
+      pam_tree_clone[i]->Write(pam_tree_clone[i]->GetName(),TObject::kOverwrite);
     };
   }
+  
+  if ( PROC ){//EMEMEMEM
+    proc_tree_clone->Write("ProcessingInfo",TObject::kOverwrite);
+  }
   cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" << endl;
 
 }