/[PAMELA software]/DarthVader/ToFLevel2/src/ToFLevel2.cpp

Diff of /DarthVader/ToFLevel2/src/ToFLevel2.cpp

Parent Directory | Revision Log | View Patch Patch

-revision 1.26 by carbone,
Fri Nov 20 11:05:21 2009 UTC
+revision 1.42 by pam-fi,
Fri Oct 17 07:54:17 2014 UTC
 Line 21 
 ToFPMT::ToFPMT(){
    adc = 0.;
    tdc_tw = 0.;
    tdc = 0.;
+   l0flag_adc = 0.;
+   l0flag_tdc = 0.;
  }
  ToFPMT::ToFPMT(const ToFPMT &t){
-Line 111 
 void ToFLevel2::Set(){//ELENA
+Line 113 
 void ToFLevel2::Set(){//ELENA
      if(!PMT)PMT = new TClonesArray("ToFPMT",12); //ELENA
      if(!ToFTrk)ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA
  }//ELENA
+ //--------------------------------------
+ //
+ //
+ //--------------------------------------
+ void ToFLevel2::SetTrackArray(TClonesArray *track){//ELENA
+     if(track && strcmp(track->GetClass()->GetName(),"ToFTrkVar")==0){
+         if(ToFTrk)ToFTrk->Clear("C");
+         ToFTrk = track;
+     }
+ }
  void ToFLevel2::Clear(Option_t *t){
    //
-Line 118 
 void ToFLevel2::Clear(Option_t *t){
+Line 130 
 void ToFLevel2::Clear(Option_t *t){
    if(PMT)PMT->Delete(); //ELENA
    memset(tof_j_flag, 0, 6*sizeof(Int_t));
    unpackError = 0;
+   unpackWarning = 0;
    //
  };
-Line 148 
 ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t
+Line 161 
 ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t
    return toftrack;
  }
+ /**
+  * Retrieves the tof track matching the seqno-th tracker stored track.
+  *
+  */
+ ToFTrkVar *ToFLevel2::GetToFStoredTrack(int seqno){
+   if( ntrk()==0 ){
+     printf("ToFLevel2::GetToFStoredTrack(int) : requested tracker SeqNo %i but no ToFrimeter tracks are stored\n",seqno);
+     return NULL;
+   };
+   ToFTrkVar *c = 0;
+   Int_t it_tof=0;
+   do {
+     c = GetToFTrkVar(it_tof);
+     it_tof++;
+   } while( c && seqno != c->trkseqno && it_tof < ntrk());
+   if(!c || seqno != c->trkseqno){
+     c = 0;
+     if(seqno!=-1 ) printf("ToFLevel2::GetToFStoredTrack(int) : requested tracker SeqNo %i does not match ToFrimeter stored tracks\n",seqno);
+   };
+   return c;
+ }
  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit){
    //
    if(ihit >= npmt()){
-Line 221 
 Int_t ToFLevel2::GetNHitPaddles(Int_t pl
+Line 262 
 Int_t ToFLevel2::GetNHitPaddles(Int_t pl
      return npad;
  };
+ /**
+  * Method to get the number of hit paddles on a ToF plane.
+  * @param plane Plane ID (11, 12, 21, 22, 31, 32) or Plane index (0,1,2,3,4,5).
+  */
+ Int_t ToFLevel2::GetTrueNHitPaddles(Int_t plane){
+     Int_t npad=0;
+     TClonesArray* Pmt = this->PMT;
+     int paddle[24];
+     memset(paddle,0, 24*sizeof(int));
+     for(int i=0; i<Pmt->GetEntries(); i++) {  //loop per vedere quale TOF � colpito
+       ToFPMT* pmthit = (ToFPMT*)Pmt->At(i);
+       int pplane = -1;
+       int ppaddle = -1;
+       GetPMTPaddle(pmthit->pmt_id,pplane,ppaddle);
+       if ( pplane == plane ) paddle[ppaddle]++;
+     }
+     for(int i=0;i<24;i++) if ( paddle[i]>0 ) npad++;
+     return npad;
+ };
  //wm Nov 08
  //gf Apr 07
  /**
-Line 235 
 Int_t ToFLevel2::GetNHitPaddles(Int_t pl
+Line 297 
 Int_t ToFLevel2::GetNHitPaddles(Int_t pl
   * @param plane Plane index (0,1,2,3,4,5)
   * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
   */
- Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){
+ Float_t ToFLevel2::GetdEdx(ToFTrkVar *trk, Int_t plane, Int_t adcfl){
+   //  printf("fiffi\n");
+   Float_t dedx = 0.;
+   Float_t PadEdx =0.;
+   Int_t SatWarning;
+   Int_t pad=-1;
+   //
+   if(!trk) return 0; //ELENA
+   //
+   if ( trk->trkseqno == -1 ){ //standalone, only paddles along the track, or about...
+     //    printf("ciccio\n");
+     Float_t xleft=0;
+     Float_t xright=0;
+     Float_t yleft=0;
+     Float_t yright=0;
+     Float_t xtof_temp[6]={100.,100.,100.,100.,100.,100.};
+     Float_t ytof_temp[6]={100.,100.,100.,100.,100.,100.};
+     ToFTrkVar *t_tof = trk;
+     xtof_temp[0]=t_tof->xtofpos[0];
+     ytof_temp[0]=t_tof->ytofpos[0];
+     xtof_temp[1]=t_tof->xtofpos[0];
+     ytof_temp[1]=t_tof->ytofpos[0];
+     xtof_temp[2]=t_tof->xtofpos[1];
+     ytof_temp[2]=t_tof->ytofpos[1];
+     xtof_temp[3]=t_tof->xtofpos[1];
+     ytof_temp[3]=t_tof->ytofpos[1];
+     xtof_temp[4]=t_tof->xtofpos[2];
+     ytof_temp[4]=t_tof->ytofpos[2];
+     xtof_temp[5]=t_tof->xtofpos[2];
+     ytof_temp[5]=t_tof->ytofpos[2];
+     if(t_tof->xtofpos[0]<100. && t_tof->ytofpos[0]<100.){
+       xtof_temp[1]=t_tof->xtofpos[0];
+       ytof_temp[0]=t_tof->ytofpos[0];
+     }else if(t_tof->xtofpos[0]>=100. && t_tof->ytofpos[0]<100.){
+       ytof_temp[0]=t_tof->ytofpos[0];
+       this->GetPaddleGeometry(0,(Int_t)log2(this->tof_j_flag[0]),xleft, xright, yleft, yright);
+       xtof_temp[1]=xleft+2.55;
+     }else if(t_tof->ytofpos[0]>=100. && t_tof->xtofpos[0]<100.){
+       xtof_temp[1]=t_tof->xtofpos[0];
+       this->GetPaddleGeometry(1,(Int_t)log2(this->tof_j_flag[1]),xleft, xright, yleft, yright);
+       ytof_temp[0]=yleft+2.75;
+     }
+     if(t_tof->xtofpos[1]<100. && t_tof->ytofpos[1]<100.){
+       xtof_temp[2]=t_tof->xtofpos[1];
+       ytof_temp[3]=t_tof->ytofpos[1];
+     }else if(t_tof->xtofpos[1]>=100. && t_tof->ytofpos[1]<100.){
+       ytof_temp[3]=t_tof->ytofpos[1];
+       this->GetPaddleGeometry(3,(Int_t)log2(this->tof_j_flag[3]),xleft, xright, yleft, yright);
+       xtof_temp[2]=xleft+4.5;
+     }else if(t_tof->ytofpos[1]>=100. && t_tof->xtofpos[1]<100.){
+       xtof_temp[2]=t_tof->xtofpos[1];
+       this->GetPaddleGeometry(2,(Int_t)log2(this->tof_j_flag[2]),xleft, xright, yleft, yright);
+       ytof_temp[3]=yleft+3.75;
+     }
+     if(t_tof->xtofpos[2]<100. && t_tof->ytofpos[2]<100.){
+       xtof_temp[5]=t_tof->xtofpos[2];
+       ytof_temp[4]=t_tof->ytofpos[2];
+     }else if(t_tof->xtofpos[2]>=100. && t_tof->ytofpos[2]<100.){
+       ytof_temp[4]=t_tof->ytofpos[2];
+       this->GetPaddleGeometry(4,(Int_t)log2(this->tof_j_flag[4]),xleft, xright, yleft, yright);
+       xtof_temp[5]=xleft+3;
+     }else if(t_tof->ytofpos[2]>=100. && t_tof->xtofpos[2]<100.){
+       xtof_temp[5]=t_tof->xtofpos[2];
+       this->GetPaddleGeometry(5,(Int_t)log2(this->tof_j_flag[5]),xleft, xright, yleft, yright);
+       ytof_temp[4]=yleft+2.5;
+     }
+     if ( (xtof_temp[0])> 17.85 && fabs(xtof_temp[0])<17.85+4. ) xtof_temp[0] =17.84;
+     if ( (xtof_temp[1])> 17.85 && fabs(xtof_temp[1])<17.85+4. ) xtof_temp[1] =17.84;
+     if ( (ytof_temp[0])> 13.75 && fabs(ytof_temp[0])<13.75+4. ) ytof_temp[0] =13.74;
+     if ( (ytof_temp[1])> 13.75 && fabs(ytof_temp[1])<13.75+4. ) ytof_temp[1] =13.74;
+     if ( (xtof_temp[0])< -17.85 && fabs(xtof_temp[0])>-(17.85+4.) ) xtof_temp[0] =-17.84;
+     if ( (xtof_temp[1])< -17.85 && fabs(xtof_temp[1])>-(17.85+4.) ) xtof_temp[1] =-17.84;
+     if ( (ytof_temp[0])< -13.75 && fabs(ytof_temp[0])>-(13.75+4.) ) ytof_temp[0] =-13.74;
+     if ( (ytof_temp[1])< -13.75 && fabs(ytof_temp[1])>-(13.75+4.) ) ytof_temp[1] =-13.74;
+     if ( (xtof_temp[2])> 4.5 && (xtof_temp[2])<4.5+4. ) xtof_temp[2] =4.4;
+     if ( (xtof_temp[3])> 4.45 && (xtof_temp[3])<4.5+4. ) xtof_temp[3] =4.4;
+     if ( (ytof_temp[2])> 3.75 && (ytof_temp[2])<3.75+4. ) ytof_temp[2] =3.74;
+     if ( (ytof_temp[3])> 3.75 && (ytof_temp[3])<3.75+4. ) ytof_temp[3] =3.74;
+     if ( (xtof_temp[2])< -4.5 && (xtof_temp[2])>-(4.5+4.) ) xtof_temp[2] =-4.4;
+     if ( (xtof_temp[3])< -4.45 && (xtof_temp[3])>-(4.5+4.) ) xtof_temp[3] =-4.4;
+     if ( (ytof_temp[2])< -3.75 && (ytof_temp[2])>-(3.75+4.) ) ytof_temp[2] =-3.74;
+     if ( (ytof_temp[3])< -3.75 && (ytof_temp[3])>-(3.75+4.) ) ytof_temp[3] =-3.74;
+     if ( (xtof_temp[4])> 6. && (xtof_temp[4])<6.+4. ) xtof_temp[2] =5.9;
+     if ( (xtof_temp[5])> 6. && (xtof_temp[5])<6.+4. ) xtof_temp[3] =5.9;
+     if ( (ytof_temp[4])> 5. && (ytof_temp[4])<5.+4. ) ytof_temp[2] =4.9;
+     if ( (ytof_temp[5])> 5. && (ytof_temp[5])<5.+4. ) ytof_temp[3] =4.9;
+     if ( (xtof_temp[4])< -6. && (xtof_temp[4])>-(6.+4.) ) xtof_temp[2] =-5.9;
+     if ( (xtof_temp[5])< -6. && (xtof_temp[5])>-(6.+4.) ) xtof_temp[3] =-5.9;
+     if ( (ytof_temp[4])< -5. && (ytof_temp[4])>-(5.+4.) ) ytof_temp[2] =-4.9;
+     if ( (ytof_temp[5])< -5. && (ytof_temp[5])>-(5.+4.) ) ytof_temp[3] =-4.9;
+     for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
+       Int_t paddleid=ii;
+       pad = GetPaddleid(plane,paddleid);
+       Int_t IpaddleT=-1;
+       IpaddleT=this->GetPaddleIdOfTrack(xtof_temp[plane],ytof_temp[plane], plane,0.0);
+ //      printf("SSS %i %i %f %f  %f %f \n",IpaddleT,paddleid,xtof_temp[plane],ytof_temp[plane],t_tof->xtofpos[plane],t_tof->ytofpos[plane]);
+       //      if ( IpaddleT == paddleid || IpaddleT-1 == paddleid || IpaddleT+1 == paddleid ){
+       if ( IpaddleT == paddleid || GetTrueNHitPaddles(plane) == 1 ){
+         //IpaddleT-1 == paddleid || IpaddleT+1 == paddleid ){
+         GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
+         dedx += PadEdx;
+       }
+     };
+   } else {
+     for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
+       Int_t paddleid=ii;
+       pad = GetPaddleid(plane,paddleid);
+       GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
+       dedx += PadEdx;
+ //      printf("TTT %i %i %f\n",paddleid,plane,PadEdx);
+     };
+   }
+   //
+   return(dedx);
+ };
+ //wm Nov 08
+ //gf Apr 07
+ /**
+  * Method to get the mean dEdx from a ToF layer - ATTENTION:
+  * It will sum up the dEdx of all the paddles, but since by definition
+  * only the paddle hitted by the track gets a dEdx value and the other
+  * paddles are set to zero, the output is just the dEdx of the hitted
+  * paddle in each layer!
+  * The "adcfl" option is not very useful (an artificial dEdx is per
+  * definition= 1 mip and not a real measurement), anyway left in the code
+  * @param notrack Track Number
+  * @param plane Plane index (0,1,2,3,4,5)
+  * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
+  */
+ Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){
+   //  printf("fiffi\n");
    Float_t dedx = 0.;
    Float_t PadEdx =0.;
    Int_t SatWarning;
-Line 245 
 Float_t ToFLevel2::GetdEdx(Int_t notrack
+Line 449 
 Float_t ToFLevel2::GetdEdx(Int_t notrack
    ToFTrkVar *trk = GetToFTrkVar(notrack);
    if(!trk) return 0; //ELENA
    //
-   for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
+   if ( trk->trkseqno == -1 ){ //standalone, only paddles along the track, or about...
-     Int_t paddleid=ii;
+     //    printf("ciccio\n");
-     pad = GetPaddleid(plane,paddleid);
+     Float_t xleft=0;
-     GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);
+     Float_t xright=0;
-     dedx += PadEdx;
+     Float_t yleft=0;
-   };
+     Float_t yright=0;
+     Float_t xtof_temp[6]={100.,100.,100.,100.,100.,100.};
+     Float_t ytof_temp[6]={100.,100.,100.,100.,100.,100.};
+     ToFTrkVar *t_tof = trk;
+     xtof_temp[0]=t_tof->xtofpos[0];
+     ytof_temp[0]=t_tof->ytofpos[0];
+     xtof_temp[1]=t_tof->xtofpos[0];
+     ytof_temp[1]=t_tof->ytofpos[0];
+     xtof_temp[2]=t_tof->xtofpos[1];
+     ytof_temp[2]=t_tof->ytofpos[1];
+     xtof_temp[3]=t_tof->xtofpos[1];
+     ytof_temp[3]=t_tof->ytofpos[1];
+     xtof_temp[4]=t_tof->xtofpos[2];
+     ytof_temp[4]=t_tof->ytofpos[2];
+     xtof_temp[5]=t_tof->xtofpos[2];
+     ytof_temp[5]=t_tof->ytofpos[2];
+     if(t_tof->xtofpos[0]<100. && t_tof->ytofpos[0]<100.){
+       xtof_temp[1]=t_tof->xtofpos[0];
+       ytof_temp[0]=t_tof->ytofpos[0];
+     }else if(t_tof->xtofpos[0]>=100. && t_tof->ytofpos[0]<100.){
+       ytof_temp[0]=t_tof->ytofpos[0];
+       this->GetPaddleGeometry(0,(Int_t)log2(this->tof_j_flag[0]),xleft, xright, yleft, yright);
+       xtof_temp[1]=xleft+2.55;
+     }else if(t_tof->ytofpos[0]>=100. && t_tof->xtofpos[0]<100.){
+       xtof_temp[1]=t_tof->xtofpos[0];
+       this->GetPaddleGeometry(1,(Int_t)log2(this->tof_j_flag[1]),xleft, xright, yleft, yright);
+       ytof_temp[0]=yleft+2.75;
+     }
+     if(t_tof->xtofpos[1]<100. && t_tof->ytofpos[1]<100.){
+       xtof_temp[2]=t_tof->xtofpos[1];
+       ytof_temp[3]=t_tof->ytofpos[1];
+     }else if(t_tof->xtofpos[1]>=100. && t_tof->ytofpos[1]<100.){
+       ytof_temp[3]=t_tof->ytofpos[1];
+       this->GetPaddleGeometry(3,(Int_t)log2(this->tof_j_flag[3]),xleft, xright, yleft, yright);
+       xtof_temp[2]=xleft+4.5;
+     }else if(t_tof->ytofpos[1]>=100. && t_tof->xtofpos[1]<100.){
+       xtof_temp[2]=t_tof->xtofpos[1];
+       this->GetPaddleGeometry(2,(Int_t)log2(this->tof_j_flag[2]),xleft, xright, yleft, yright);
+       ytof_temp[3]=yleft+3.75;
+     }
+     if(t_tof->xtofpos[2]<100. && t_tof->ytofpos[2]<100.){
+       xtof_temp[5]=t_tof->xtofpos[2];
+       ytof_temp[4]=t_tof->ytofpos[2];
+     }else if(t_tof->xtofpos[2]>=100. && t_tof->ytofpos[2]<100.){
+       ytof_temp[4]=t_tof->ytofpos[2];
+       this->GetPaddleGeometry(4,(Int_t)log2(this->tof_j_flag[4]),xleft, xright, yleft, yright);
+       xtof_temp[5]=xleft+3;
+     }else if(t_tof->ytofpos[2]>=100. && t_tof->xtofpos[2]<100.){
+       xtof_temp[5]=t_tof->xtofpos[2];
+       this->GetPaddleGeometry(5,(Int_t)log2(this->tof_j_flag[5]),xleft, xright, yleft, yright);
+       ytof_temp[4]=yleft+2.5;
+     }
+     if ( (xtof_temp[0])> 17.85 && fabs(xtof_temp[0])<17.85+4. ) xtof_temp[0] =17.84;
+     if ( (xtof_temp[1])> 17.85 && fabs(xtof_temp[1])<17.85+4. ) xtof_temp[1] =17.84;
+     if ( (ytof_temp[0])> 13.75 && fabs(ytof_temp[0])<13.75+4. ) ytof_temp[0] =13.74;
+     if ( (ytof_temp[1])> 13.75 && fabs(ytof_temp[1])<13.75+4. ) ytof_temp[1] =13.74;
+     if ( (xtof_temp[0])< -17.85 && fabs(xtof_temp[0])>-(17.85+4.) ) xtof_temp[0] =-17.84;
+     if ( (xtof_temp[1])< -17.85 && fabs(xtof_temp[1])>-(17.85+4.) ) xtof_temp[1] =-17.84;
+     if ( (ytof_temp[0])< -13.75 && fabs(ytof_temp[0])>-(13.75+4.) ) ytof_temp[0] =-13.74;
+     if ( (ytof_temp[1])< -13.75 && fabs(ytof_temp[1])>-(13.75+4.) ) ytof_temp[1] =-13.74;
+     if ( (xtof_temp[2])> 4.5 && (xtof_temp[2])<4.5+4. ) xtof_temp[2] =4.4;
+     if ( (xtof_temp[3])> 4.45 && (xtof_temp[3])<4.5+4. ) xtof_temp[3] =4.4;
+     if ( (ytof_temp[2])> 3.75 && (ytof_temp[2])<3.75+4. ) ytof_temp[2] =3.74;
+     if ( (ytof_temp[3])> 3.75 && (ytof_temp[3])<3.75+4. ) ytof_temp[3] =3.74;
+     if ( (xtof_temp[2])< -4.5 && (xtof_temp[2])>-(4.5+4.) ) xtof_temp[2] =-4.4;
+     if ( (xtof_temp[3])< -4.45 && (xtof_temp[3])>-(4.5+4.) ) xtof_temp[3] =-4.4;
+     if ( (ytof_temp[2])< -3.75 && (ytof_temp[2])>-(3.75+4.) ) ytof_temp[2] =-3.74;
+     if ( (ytof_temp[3])< -3.75 && (ytof_temp[3])>-(3.75+4.) ) ytof_temp[3] =-3.74;
+     if ( (xtof_temp[4])> 6. && (xtof_temp[4])<6.+4. ) xtof_temp[2] =5.9;
+     if ( (xtof_temp[5])> 6. && (xtof_temp[5])<6.+4. ) xtof_temp[3] =5.9;
+     if ( (ytof_temp[4])> 5. && (ytof_temp[4])<5.+4. ) ytof_temp[2] =4.9;
+     if ( (ytof_temp[5])> 5. && (ytof_temp[5])<5.+4. ) ytof_temp[3] =4.9;
+     if ( (xtof_temp[4])< -6. && (xtof_temp[4])>-(6.+4.) ) xtof_temp[2] =-5.9;
+     if ( (xtof_temp[5])< -6. && (xtof_temp[5])>-(6.+4.) ) xtof_temp[3] =-5.9;
+     if ( (ytof_temp[4])< -5. && (ytof_temp[4])>-(5.+4.) ) ytof_temp[2] =-4.9;
+     if ( (ytof_temp[5])< -5. && (ytof_temp[5])>-(5.+4.) ) ytof_temp[3] =-4.9;
+     for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
+       Int_t paddleid=ii;
+       pad = GetPaddleid(plane,paddleid);
+       Int_t IpaddleT=-1;
+       IpaddleT=this->GetPaddleIdOfTrack(xtof_temp[plane],ytof_temp[plane], plane,0.0);
+ //      printf("SSS %i %i %f %f  %f %f \n",IpaddleT,paddleid,xtof_temp[plane],ytof_temp[plane],t_tof->xtofpos[plane],t_tof->ytofpos[plane]);
+       //      if ( IpaddleT == paddleid || IpaddleT-1 == paddleid || IpaddleT+1 == paddleid ){
+       if ( IpaddleT == paddleid || GetTrueNHitPaddles(plane) == 1 ){
+         //IpaddleT-1 == paddleid || IpaddleT+1 == paddleid ){
+         GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);
+         dedx += PadEdx;
+       }
+     };
+   } else {
+     for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
+       Int_t paddleid=ii;
+       pad = GetPaddleid(plane,paddleid);
+       GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);
+       dedx += PadEdx;
+ //      printf("TTT %i %i %f\n",paddleid,plane,PadEdx);
+     };
+   }
    //
    return(dedx);
  };
-Line 404 
 void ToFLevel2::GetPMTIndex(Int_t ind, I
+Line 716 
 void ToFLevel2::GetPMTIndex(Int_t ind, I
   */
  void ToFLevel2::GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){
- /*
+   /*
- Float_t  PMTsat[48] = {
+     Float_t  PMTsat[48] = {
 .14, 3165.48, 3153.85, 3085.73, 3089.65, 3107.64, 3097.52, 3078.37,
 .05, 3087.07, 3112.22, 3102.92, 3080.58, 3092.55, 3087.94, 3125.03,
 .09, 3143.16, 3125.51, 3181.27, 3092.09, 3124.98, 3069.3, 3095.53,
 .11, 3133.53, 3114.73, 3113.01, 3091.19, 3097.99, 3033.84, 3134.98,
 .37, 3111.04, 3066.77, 3108.17, 3133, 3111.06, 3052.52, 3140.66,
 .33, 3094.85, 3150.85, 3118.8, 3096.24, 3118.47,3111.36, 3117.11 } ;
- */
+   */
- // new values from Napoli dec 2008
+   // new values from Napoli dec 2008
- Float_t  PMTsat[48] = {
+   Float_t  PMTsat[48] = {
 .35,3178.19,3167.38,3099.73,3117.00,3126.29,3111.44,3092.27,
 .48,3094.41,3132.13,3115.37,3099.32,3110.97,3111.80,3143.14,
 .72,3153.44,3136.00,3188.96,3104.73,3140.45,3073.18,3106.62,
 .48,3146.92,3127.24,3136.52,3109.59,3112.89,3045.15,3147.26,
 .92,3121.05,3083.25,3123.62,3150.92,3125.30,3067.60,3160.18,
 .36,3108.92,3164.77,3133.64,3111.47,3131.98,3128.87,3135.56 };
- for (Int_t i=0; i<48;i++) PMTsat[i] = PMTsat[i] - 5.;  // safety margin
+   for (Int_t i=0; i<48;i++) PMTsat[i] = PMTsat[i] - 5.;  // safety margin
    PadEdx = 0.;
- //  SatWarning = 1000;
+   //  SatWarning = 1000;
    SatWarning = 0;   // 0=good, increase for each bad PMT
    Float_t dEdx[48] = {0};
-Line 474 
 for (Int_t i=0; i<48;i++) PMTsat[i] = PM
+Line 786 
 for (Int_t i=0; i<48;i++) PMTsat[i] = PM
    }
- //  if( adcraw[pmtleft] >3000 || adcraw[pmtright] >3000)SatWarning=1;  //old version
+   //  if( adcraw[pmtleft] >3000 || adcraw[pmtright] >3000)SatWarning=1;  //old version
- // Increase SatWarning Counter for each PMT>Sat
+   // Increase SatWarning Counter for each PMT>Sat
    if( adcraw[pmtleft] > PMTsat[pmtleft])SatWarning++;
    if( adcraw[pmtright] > PMTsat[pmtright])SatWarning++;
- // if ADC  > sat set dEdx=1000
+   // if ADC  > sat set dEdx=1000
    if( adcraw[pmtleft] > PMTsat[pmtleft]) dEdx[pmtleft] = 1000.;
    if( adcraw[pmtright] > PMTsat[pmtright]) dEdx[pmtright] = 1000. ;
- // if two PMT are good, take mean dEdx, otherwise only the good dEdx
+   // if two PMT are good, take mean dEdx, otherwise only the good dEdx
    if(dEdx[pmtleft]<1000 && dEdx[pmtright]<1000) PadEdx = (dEdx[pmtleft]+dEdx[pmtright])*0.5;
    if(dEdx[pmtleft]==1000 && dEdx[pmtright]<1000) PadEdx = dEdx[pmtright];
    if(dEdx[pmtleft]<1000 && dEdx[pmtright]==1000) PadEdx = dEdx[pmtleft];
  };
  //
+ //  wm Nov 08 revision - saturation values included
+ /// gf Apr 07
+ /**
+  * Method to get the dEdx from a given ToF paddle.
+  * If two PMTs are good, the mean dEdx of both PMTs is taken, otherwise
+  * just the dEdx of the "good" PMT. If both PMTs are above saturation => dEdx=1000
+  * @param notrack Track Number
+  * @param Paddle index (0,1,...,23).
+  * @param adcflag in the paddle (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
+  * @param PadEdx dEdx from a given ToF paddle
+  * @param SatWarning 1 if the PMT ios near saturation region (adcraw ~3000)
+  */
+ void ToFLevel2::GetdEdxPaddle(ToFTrkVar *trk, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){
+   /*
+     Float_t  PMTsat[48] = {
+.14, 3165.48, 3153.85, 3085.73, 3089.65, 3107.64, 3097.52, 3078.37,
+.05, 3087.07, 3112.22, 3102.92, 3080.58, 3092.55, 3087.94, 3125.03,
+.09, 3143.16, 3125.51, 3181.27, 3092.09, 3124.98, 3069.3, 3095.53,
+.11, 3133.53, 3114.73, 3113.01, 3091.19, 3097.99, 3033.84, 3134.98,
+.37, 3111.04, 3066.77, 3108.17, 3133, 3111.06, 3052.52, 3140.66,
+.33, 3094.85, 3150.85, 3118.8, 3096.24, 3118.47,3111.36, 3117.11 } ;
+   */
+   // new values from Napoli dec 2008
+   Float_t  PMTsat[48] = {
+.35,3178.19,3167.38,3099.73,3117.00,3126.29,3111.44,3092.27,
+.48,3094.41,3132.13,3115.37,3099.32,3110.97,3111.80,3143.14,
+.72,3153.44,3136.00,3188.96,3104.73,3140.45,3073.18,3106.62,
+.48,3146.92,3127.24,3136.52,3109.59,3112.89,3045.15,3147.26,
+.92,3121.05,3083.25,3123.62,3150.92,3125.30,3067.60,3160.18,
+.36,3108.92,3164.77,3133.64,3111.47,3131.98,3128.87,3135.56 };
+   for (Int_t i=0; i<48;i++) PMTsat[i] = PMTsat[i] - 5.;  // safety margin
+   PadEdx = 0.;
+   //  SatWarning = 1000;
+   SatWarning = 0;   // 0=good, increase for each bad PMT
+   Float_t dEdx[48] = {0};
+   Int_t pmt_id = -1;
+   Float_t adcraw[48];
+   //
+   if(!trk) return; //ELENA
+   //
+   Int_t pmtleft=-1;
+   Int_t pmtright=-1;
+   GetPaddlePMT(paddleid, pmtleft, pmtright);
+   adcraw[pmtleft] = 4095;
+   adcraw[pmtright] = 4095;
+   for (Int_t jj=0; jj<npmt(); jj++){
+     ToFPMT *pmt = GetToFPMT(jj);
+     if(!pmt)break; //ELENA
+     pmt_id = pmt->pmt_id;
+     if(pmt_id==pmtleft){
+       adcraw[pmtleft] = pmt->adc;
+     }
+     if(pmt_id==pmtright){
+       adcraw[pmtright] = pmt->adc;
+     }
+   }
+   for (Int_t i=0; i<trk->npmtadc; i++){
+     if((trk->adcflag).At(i)==0 || adcfl==100){
+       if((trk->pmtadc).At(i) == pmtleft)dEdx[pmtleft] = (trk->dedx).At(i);
+       if((trk->pmtadc).At(i) == pmtright)dEdx[pmtright] = (trk->dedx).At(i);
+     }else{
+       if((trk->pmtadc).At(i) == pmtleft)dEdx[pmtleft] = 0.;
+       if((trk->pmtadc).At(i) == pmtright)dEdx[pmtright] = 0.;
+     }
+   }
+   //  if( adcraw[pmtleft] >3000 || adcraw[pmtright] >3000)SatWarning=1;  //old version
+   // Increase SatWarning Counter for each PMT>Sat
+   if( adcraw[pmtleft] > PMTsat[pmtleft])SatWarning++;
+   if( adcraw[pmtright] > PMTsat[pmtright])SatWarning++;
+   // if ADC  > sat set dEdx=1000
+   if( adcraw[pmtleft] > PMTsat[pmtleft]) dEdx[pmtleft] = 1000.;
+   if( adcraw[pmtright] > PMTsat[pmtright]) dEdx[pmtright] = 1000. ;
+   // if two PMT are good, take mean dEdx, otherwise only the good dEdx
+   if(dEdx[pmtleft]<1000 && dEdx[pmtright]<1000) PadEdx = (dEdx[pmtleft]+dEdx[pmtright])*0.5;
+   if(dEdx[pmtleft]==1000 && dEdx[pmtright]<1000) PadEdx = dEdx[pmtright];
+   if(dEdx[pmtleft]<1000 && dEdx[pmtright]==1000) PadEdx = dEdx[pmtleft];
+ };
  // gf Apr 07
  /**
-Line 1149 
 Int_t ToFLevel2::Process(TrkLevel2 *trk,
+Line 1560 
 Int_t ToFLevel2::Process(TrkLevel2 *trk,
    //
    // Copiare qui qualcosa di simile a calonuclei per evitare di riprocessare sempre tutto
    //
+   printf("\n\n\n ERROR: NOT IMPLEMENTED ANYMORE, write Emiliano if you need this method (Emiliano.Mocchiutti@ts.infn.it) \n\n\n");
+   return(-1);
+   //   //
+   //   // structures to communicate with F77
-   //
+   //   //
-   // structures to communicate with F77
+   //   extern struct ToFInput  tofinput_;
-   //
+ //   extern struct ToFOutput tofoutput_;
-   extern struct ToFInput  tofinput_;
+ //   //
-   extern struct ToFOutput tofoutput_;
+ //   // DB connection
-   //
+ //   //
-   // DB connection
+ //   TString host;
-   //
+ //   TString user;
-   TString host;
+ //   TString psw;
-   TString user;
+ //   const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");
-   TString psw;
+ //   const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");
-   const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");
+ //   const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");
-   const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");
+ //   if ( !pamdbhost ) pamdbhost = "";
-   const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");
+ //   if ( !pamdbuser ) pamdbuser = "";
-   if ( !pamdbhost ) pamdbhost = "";
+ //   if ( !pamdbpsw ) pamdbpsw = "";
-   if ( !pamdbuser ) pamdbuser = "";
+ //   if ( strcmp(pamdbhost,"") ) host = pamdbhost;
-   if ( !pamdbpsw ) pamdbpsw = "";
+ //   if ( strcmp(pamdbuser,"") ) user = pamdbuser;
-   if ( strcmp(pamdbhost,"") ) host = pamdbhost;
+ //   if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;
-   if ( strcmp(pamdbuser,"") ) user = pamdbuser;
+ //   //
-   if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;
+ //   //
-   //
+ //   TSQLServer *dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());
-   //
+ //   if ( !dbc->IsConnected() ) return 1;
-   TSQLServer *dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());
+ //   stringstream myquery;
-   if ( !dbc->IsConnected() ) return 1;
+ //   myquery.str("");
-   stringstream myquery;
+ //   myquery << "SET time_zone='+0:00';";
-   myquery.str("");
+ //   dbc->Query(myquery.str().c_str());
-   myquery << "SET time_zone='+0:00'";
+ //   delete dbc->Query("SET sql_mode = 'NO_UNSIGNED_SUBTRACTION';");
-   dbc->Query(myquery.str().c_str());
+ //   GL_PARAM *glparam = new GL_PARAM();
-   GL_PARAM *glparam = new GL_PARAM();
+ //   glparam->Query_GL_PARAM(1,1,dbc); // parameters stored in DB in GL_PRAM table
-   glparam->Query_GL_PARAM(1,1,dbc); // parameters stored in DB in GL_PRAM table
+ //   trk->LoadField(glparam->PATH+glparam->NAME);
-   trk->LoadField(glparam->PATH+glparam->NAME);
+ //   //
-   //
+ //   Bool_t defcal = true;
-   Bool_t defcal = true;
+ //   Int_t error=glparam->Query_GL_PARAM(run->RUNHEADER_TIME,201,dbc); // parameters stored in DB in GL_PRAM table
-   Int_t error=glparam->Query_GL_PARAM(run->RUNHEADER_TIME,201,dbc); // parameters stored in DB in GL_PRAM table
+ //   if ( error<0 ) {
-   if ( error<0 ) {
+ //     return(1);
-     return(1);
+ //   };
-   };
+ //   printf(" Reading ToF parameter file: %s \n",(glparam->PATH+glparam->NAME).Data());
-   printf(" Reading ToF parameter file: %s \n",(glparam->PATH+glparam->NAME).Data());
+ //   if ( (UInt_t)glparam->TO_TIME != (UInt_t)4294967295UL ) defcal = false;
-   if ( (UInt_t)glparam->TO_TIME != (UInt_t)4294967295UL ) defcal = false;
+ //   //
-   //
+ //   Int_t nlen = (Int_t)(glparam->PATH+glparam->NAME).Length();
-   Int_t nlen = (Int_t)(glparam->PATH+glparam->NAME).Length();
+ //   rdtofcal((char *)(glparam->PATH+glparam->NAME).Data(),&nlen);
-   rdtofcal((char *)(glparam->PATH+glparam->NAME).Data(),&nlen);
+ //   //
-   //
+ //   Int_t adc[4][12];
-   Int_t adc[4][12];
+ //   Int_t tdc[4][12];
-   Int_t tdc[4][12];
+ //   Float_t tdcc[4][12];
-   Float_t tdcc[4][12];
+ //   //
-   //
+ //   // process tof data
-   // process tof data
+ //   //
-   //
+ //   for (Int_t hh=0; hh<12;hh++){
-   for (Int_t hh=0; hh<12;hh++){
+ //     for (Int_t kk=0; kk<4;kk++){
-     for (Int_t kk=0; kk<4;kk++){
+ //            adc[kk][hh] = 4095;
-            adc[kk][hh] = 4095;
+ //            tdc[kk][hh] = 4095;
-            tdc[kk][hh] = 4095;
+ //            tdcc[kk][hh] = 4095.;
-            tdcc[kk][hh] = 4095.;
+ //            tofinput_.adc[hh][kk] = 4095;
-            tofinput_.adc[hh][kk] = 4095;
+ //            tofinput_.tdc[hh][kk] = 4095;
-            tofinput_.tdc[hh][kk] = 4095;
+ //     };
-     };
+ //   };
-   };
+ //   Int_t ntrkentry = 0;
-   Int_t ntrkentry = 0;
+ //   Int_t npmtentry = 0;
-   Int_t npmtentry = 0;
+ //   Int_t gg = 0;
-   Int_t gg = 0;
+ //   Int_t hh = 0;
-   Int_t hh = 0;
+ //   Int_t adcf[48];
-   Int_t adcf[48];
+ //   memset(adcf, 0, 48*sizeof(Int_t));
-   memset(adcf, 0, 48*sizeof(Int_t));
+ //   Int_t tdcf[48];
-   Int_t tdcf[48];
+ //   memset(tdcf, 0, 48*sizeof(Int_t));
-   memset(tdcf, 0, 48*sizeof(Int_t));
+ //   for (Int_t pm=0; pm < this->ntrk() ; pm++){
-   for (Int_t pm=0; pm < this->ntrk() ; pm++){
+ //      ToFTrkVar *ttf = this->GetToFTrkVar(pm);
-      ToFTrkVar *ttf = this->GetToFTrkVar(pm);
+ //      for ( Int_t nc=0; nc < ttf->npmttdc; nc++){
-      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++){
-      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;
+ //      };
-      };
+ //   };
-   };
+ //   //
-   //
+ //   for (Int_t pm=0; pm < this->npmt() ; pm++){
-   for (Int_t pm=0; pm < this->npmt() ; pm++){
+ //      ToFPMT *pmt = this->GetToFPMT(pm);
-      ToFPMT *pmt = this->GetToFPMT(pm);
+ //      this->GetPMTIndex(pmt->pmt_id, gg, hh);
-      this->GetPMTIndex(pmt->pmt_id, gg, hh);
+ //      if ( adcf[pmt->pmt_id] == 0 ){
-      if ( adcf[pmt->pmt_id] == 0 ){
+ //              tofinput_.adc[gg][hh] = (int)pmt->adc;
-              tofinput_.adc[gg][hh] = (int)pmt->adc;
+ //              adc[hh][gg] = (int)pmt->adc;
-              adc[hh][gg] = (int)pmt->adc;
+ //      };
-      };
+ //      if ( tdcf[pmt->pmt_id] == 0 ){
-      if ( tdcf[pmt->pmt_id] == 0 ){
+ //              tofinput_.tdc[gg][hh] = (int)pmt->tdc;
-              tofinput_.tdc[gg][hh] = (int)pmt->tdc;
+ //              tdc[hh][gg] = (int)pmt->tdc;
-              tdc[hh][gg] = (int)pmt->tdc;
+ //      };
-      };
+ //      tdcc[hh][gg] = (float)pmt->tdc_tw;
-      tdcc[hh][gg] = (float)pmt->tdc_tw;
+ //      // Int_t pppid = this->GetPMTid(hh,gg);
-      // Int_t pppid = this->GetPMTid(hh,gg);
+ //      //      printf(" pm %i pmt_id %i pppid %i hh %i gg %i tdcc %f tdc %f adc %f \n",pm,pmt->pmt_id,pppid,hh,gg,pmt->tdc_tw,pmt->tdc,pmt->adc);
-      //      printf(" pm %i pmt_id %i pppid %i hh %i gg %i tdcc %f tdc %f adc %f \n",pm,pmt->pmt_id,pppid,hh,gg,pmt->tdc_tw,pmt->tdc,pmt->adc);
+ //   };
-   };
+ //   //
-   //
+ //   Int_t unpackError = this->unpackError;
-   Int_t unpackError = this->unpackError;
+ //   //
-   //
+ //   for (Int_t hh=0; hh<5;hh++){
-   for (Int_t hh=0; hh<5;hh++){
+ //      tofinput_.patterntrig[hh]=trg->patterntrig[hh];
-      tofinput_.patterntrig[hh]=trg->patterntrig[hh];
+ //   };
-   };
+ //   //
-   //
+ //   this->Clear();
-   this->Clear();
+ //   //
-   //
+ //       Int_t pmt_id = 0;
-       Int_t pmt_id = 0;
+ //       ToFPMT *t_pmt = new ToFPMT();
-       ToFPMT *t_pmt = new ToFPMT();
+ //       if(!(this->PMT)) this->PMT = new TClonesArray("ToFPMT",12); //ELENA
-       if(!(this->PMT)) this->PMT = new TClonesArray("ToFPMT",12); //ELENA
+ //       TClonesArray &tpmt = *this->PMT;
-       TClonesArray &tpmt = *this->PMT;
+ //       ToFTrkVar *t_tof = new ToFTrkVar();
-       ToFTrkVar *t_tof = new ToFTrkVar();
+ //       if(!(this->ToFTrk)) this->ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA
-       if(!(this->ToFTrk)) this->ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA
+ //       TClonesArray &t = *this->ToFTrk;
-       TClonesArray &t = *this->ToFTrk;
+ //       //
-       //
+ //       //
-       //
+ //       // Here we have calibrated data, ready to be passed to the FORTRAN routine which will extract common and track-related  variables.
-       // Here we have calibrated data, ready to be passed to the FORTRAN routine which will extract common and track-related  variables.
+ //       //
-       //
+ //       npmtentry = 0;
-       npmtentry = 0;
+ //       //
-       //
+ //       ntrkentry = 0;
-       ntrkentry = 0;
+ //       //
-       //
+ //       // Calculate tracks informations from ToF alone
-       // Calculate tracks informations from ToF alone
+ //       //
-       //
+ //       tofl2com();
-       tofl2com();
+ //       //
-       //
+ //       memcpy(this->tof_j_flag,tofoutput_.tof_j_flag,6*sizeof(Int_t));
-       memcpy(this->tof_j_flag,tofoutput_.tof_j_flag,6*sizeof(Int_t));
+ //       //
-       //
+ //       t_tof->trkseqno = -1;
-       t_tof->trkseqno = -1;
+ //       //
-       //
+ //       // and now we must copy from the output structure to the level2 class:
-       // and now we must copy from the output structure to the level2 class:
+ //       //
-       //
+ //       t_tof->npmttdc = 0;
-       t_tof->npmttdc = 0;
+ //       //
-       //
+ //       for (Int_t hh=0; hh<12;hh++){
-       for (Int_t hh=0; hh<12;hh++){
+ //         for (Int_t kk=0; kk<4;kk++){
-         for (Int_t kk=0; kk<4;kk++){
+ //           if ( tofoutput_.tofmask[hh][kk] != 0 ){
-           if ( tofoutput_.tofmask[hh][kk] != 0 ){
+ //             pmt_id = this->GetPMTid(kk,hh);
-             pmt_id = this->GetPMTid(kk,hh);
+ //             t_tof->pmttdc.AddAt(pmt_id,t_tof->npmttdc);
-             t_tof->pmttdc.AddAt(pmt_id,t_tof->npmttdc);
+ //             t_tof->tdcflag.AddAt(tofoutput_.tdcflagtof[hh][kk],t_tof->npmttdc); // gf: Jan 09/07
-             t_tof->tdcflag.AddAt(tofoutput_.tdcflagtof[hh][kk],t_tof->npmttdc); // gf: Jan 09/07
+ //             t_tof->npmttdc++;
-             t_tof->npmttdc++;
+ //           };
-           };
+ //         };
-         };
+ //       };
-       };
+ //       for (Int_t kk=0; kk<13;kk++){
-       for (Int_t kk=0; kk<13;kk++){
+ //         t_tof->beta[kk] = tofoutput_.betatof_a[kk];
-         t_tof->beta[kk] = tofoutput_.betatof_a[kk];
+ //       }
-       }
+ //       //
-       //
+ //       t_tof->npmtadc = 0;
-       t_tof->npmtadc = 0;
+ //       for (Int_t hh=0; hh<12;hh++){
-       for (Int_t hh=0; hh<12;hh++){
+ //         for (Int_t kk=0; kk<4;kk++){
-         for (Int_t kk=0; kk<4;kk++){
+ //           if ( tofoutput_.adctof_c[hh][kk] < 1000 ){
-           if ( tofoutput_.adctof_c[hh][kk] < 1000 ){
+ //             t_tof->dedx.AddAt(tofoutput_.adctof_c[hh][kk],t_tof->npmtadc);
-             t_tof->dedx.AddAt(tofoutput_.adctof_c[hh][kk],t_tof->npmtadc);
+ //             pmt_id = this->GetPMTid(kk,hh);
-             pmt_id = this->GetPMTid(kk,hh);
+ //             t_tof->pmtadc.AddAt(pmt_id,t_tof->npmtadc);
-             t_tof->pmtadc.AddAt(pmt_id,t_tof->npmtadc);
+ //             t_tof->adcflag.AddAt(tofoutput_.adcflagtof[hh][kk],t_tof->npmtadc); // gf: Jan 09/07
-             t_tof->adcflag.AddAt(tofoutput_.adcflagtof[hh][kk],t_tof->npmtadc); // gf: Jan 09/07
+ //             t_tof->npmtadc++;
-             t_tof->npmtadc++;
+ //           };
-           };
+ //         };
-         };
+ //       };
-       };
+ //       //
-       //
+ //       memcpy(t_tof->xtofpos,tofoutput_.xtofpos,sizeof(t_tof->xtofpos));
-       memcpy(t_tof->xtofpos,tofoutput_.xtofpos,sizeof(t_tof->xtofpos));
+ //       memcpy(t_tof->ytofpos,tofoutput_.ytofpos,sizeof(t_tof->ytofpos));
-       memcpy(t_tof->ytofpos,tofoutput_.ytofpos,sizeof(t_tof->ytofpos));
+ //       memcpy(t_tof->xtr_tof,tofoutput_.xtr_tof,sizeof(t_tof->xtr_tof));
-       memcpy(t_tof->xtr_tof,tofoutput_.xtr_tof,sizeof(t_tof->xtr_tof));
+ //       memcpy(t_tof->ytr_tof,tofoutput_.ytr_tof,sizeof(t_tof->ytr_tof));
-       memcpy(t_tof->ytr_tof,tofoutput_.ytr_tof,sizeof(t_tof->ytr_tof));
+ //       //
-       //
+ //       new(t[ntrkentry]) ToFTrkVar(*t_tof);
-       new(t[ntrkentry]) ToFTrkVar(*t_tof);
+ //       ntrkentry++;
-       ntrkentry++;
+ //       t_tof->Clear();
-       t_tof->Clear();
+ //       //
-       //
+ //       //
-       //
+ //       //
-       //
+ //       t_pmt->Clear();
-       t_pmt->Clear();
+ //       //
-       //
+ //       for (Int_t hh=0; hh<12;hh++){
-       for (Int_t hh=0; hh<12;hh++){
+ //         for (Int_t kk=0; kk<4;kk++){
-         for (Int_t kk=0; kk<4;kk++){
+ //          // new WM
-          // new WM
+ //           if ( tofoutput_.tdc_c[hh][kk] < 4095 || adc[kk][hh] < 4095  || tdc[kk][hh] < 4095 ){
-           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 ){
- //          if ( tdcc[kk][hh] < 4095. || adc[kk][hh] < 4095  || tdc[kk][hh] < 4095 ){
+ //             //
-             //
+ //             t_pmt->pmt_id = this->GetPMTid(kk,hh);
-             t_pmt->pmt_id = this->GetPMTid(kk,hh);
+ //             t_pmt->tdc_tw = tofoutput_.tdc_c[hh][kk];
-             t_pmt->tdc_tw = tofoutput_.tdc_c[hh][kk];
+ //             t_pmt->adc = (Float_t)adc[kk][hh];
-             t_pmt->adc = (Float_t)adc[kk][hh];
+ //             t_pmt->tdc = (Float_t)tdc[kk][hh];
-             t_pmt->tdc = (Float_t)tdc[kk][hh];
+ //             //
-             //
+ //             new(tpmt[npmtentry]) ToFPMT(*t_pmt);
-             new(tpmt[npmtentry]) ToFPMT(*t_pmt);
+ //             npmtentry++;
-             npmtentry++;
+ //             t_pmt->Clear();
-             t_pmt->Clear();
+ //           };
-           };
+ //         };
-         };
+ //       };
-       };
+ //       //
-       //
+ //       // Calculate track-related variables
-       // Calculate track-related variables
+ //       //
-       //
+ //       if ( trk->ntrk() > 0 ){
-       if ( trk->ntrk() > 0 ){
+ //         //
-         //
+ //         // We have at least one track
-         // We have at least one track
+ //         //
-         //
+ //         //
-         //
+ //         // Run over tracks
-         // Run over tracks
+ //         //
-         //
+ //         for(Int_t nt=0; nt < trk->ntrk(); nt++){
-         for(Int_t nt=0; nt < trk->ntrk(); nt++){
+ //           //
-           //
+ //           TrkTrack *ptt = trk->GetStoredTrack(nt);
-           TrkTrack *ptt = trk->GetStoredTrack(nt);
+ //           //
-           //
+ //           // Copy the alpha vector in the input structure
-           // Copy the alpha vector in the input structure
+ //           //
-           //
+ //           for (Int_t e = 0; e < 5 ; e++){
-           for (Int_t e = 0; e < 5 ; e++){
+ //             tofinput_.al_pp[e] = ptt->al[e];
-             tofinput_.al_pp[e] = ptt->al[e];
+ //           };
-           };
+ //           //
-           //
+ //           // Get tracker related variables for this track
-           // Get tracker related variables for this track
+ //           //
-           //
+ //           toftrk();
-           toftrk();
+ //           //
-           //
+ //           // Copy values in the class from the structure (we need to use a temporary class to store variables).
-           // Copy values in the class from the structure (we need to use a temporary class to store variables).
+ //           //
-           //
+ //           t_tof->npmttdc = 0;
-           t_tof->npmttdc = 0;
+ //           for (Int_t hh=0; hh<12;hh++){
-           for (Int_t hh=0; hh<12;hh++){
+ //             for (Int_t kk=0; kk<4;kk++){
-             for (Int_t kk=0; kk<4;kk++){
+ //               if ( tofoutput_.tofmask[hh][kk] != 0 ){
-               if ( tofoutput_.tofmask[hh][kk] != 0 ){
+ //                 pmt_id = this->GetPMTid(kk,hh);
-                 pmt_id = this->GetPMTid(kk,hh);
+ //                 t_tof->pmttdc.AddAt(pmt_id,t_tof->npmttdc);
-                 t_tof->pmttdc.AddAt(pmt_id,t_tof->npmttdc);
+ //                 t_tof->tdcflag.AddAt(tofoutput_.tdcflag[hh][kk],t_tof->npmttdc); // gf: Jan 09/07
-                 t_tof->tdcflag.AddAt(tofoutput_.tdcflag[hh][kk],t_tof->npmttdc); // gf: Jan 09/07
+ //                 t_tof->npmttdc++;
-                 t_tof->npmttdc++;
+ //               };
-               };
+ //             };
-             };
+ //           };
-           };
+ //           for (Int_t kk=0; kk<13;kk++){
-           for (Int_t kk=0; kk<13;kk++){
+ //             t_tof->beta[kk] = tofoutput_.beta_a[kk];
-             t_tof->beta[kk] = tofoutput_.beta_a[kk];
+ //           };
-           };
+ //           //
-           //
+ //           t_tof->npmtadc = 0;
-           t_tof->npmtadc = 0;
+ //           for (Int_t hh=0; hh<12;hh++){
-           for (Int_t hh=0; hh<12;hh++){
+ //             for (Int_t kk=0; kk<4;kk++){
-             for (Int_t kk=0; kk<4;kk++){
+ //               if ( tofoutput_.adc_c[hh][kk] < 1000 ){
-               if ( tofoutput_.adc_c[hh][kk] < 1000 ){
+ //                 t_tof->dedx.AddAt(tofoutput_.adc_c[hh][kk],t_tof->npmtadc);
-                 t_tof->dedx.AddAt(tofoutput_.adc_c[hh][kk],t_tof->npmtadc);
+ //                 pmt_id = this->GetPMTid(kk,hh);
-                 pmt_id = this->GetPMTid(kk,hh);
+ //                 t_tof->pmtadc.AddAt(pmt_id,t_tof->npmtadc);
-                 t_tof->pmtadc.AddAt(pmt_id,t_tof->npmtadc);
+ //                 t_tof->adcflag.AddAt(tofoutput_.adcflag[hh][kk],t_tof->npmtadc); // gf: Jan 09/07
-                 t_tof->adcflag.AddAt(tofoutput_.adcflag[hh][kk],t_tof->npmtadc); // gf: Jan 09/07
+ //                 t_tof->npmtadc++;
-                 t_tof->npmtadc++;
+ //               };
-               };
+ //             };
-             };
+ //           };
-           };
+ //           //
-           //
+ //           memcpy(t_tof->xtofpos,tofoutput_.xtofpos,sizeof(t_tof->xtofpos));
-           memcpy(t_tof->xtofpos,tofoutput_.xtofpos,sizeof(t_tof->xtofpos));
+ //           memcpy(t_tof->ytofpos,tofoutput_.ytofpos,sizeof(t_tof->ytofpos));
-           memcpy(t_tof->ytofpos,tofoutput_.ytofpos,sizeof(t_tof->ytofpos));
+ //           memcpy(t_tof->xtr_tof,tofoutput_.xtr_tof,sizeof(t_tof->xtr_tof));
-           memcpy(t_tof->xtr_tof,tofoutput_.xtr_tof,sizeof(t_tof->xtr_tof));
+ //           memcpy(t_tof->ytr_tof,tofoutput_.ytr_tof,sizeof(t_tof->ytr_tof));
-           memcpy(t_tof->ytr_tof,tofoutput_.ytr_tof,sizeof(t_tof->ytr_tof));
+ //           //
-           //
+ //           // Store the tracker track number in order to be sure to have shyncronized data during analysis
-           // Store the tracker track number in order to be sure to have shyncronized data during analysis
+ //           //
-           //
+ //           t_tof->trkseqno = nt;
-           t_tof->trkseqno = nt;
+ //           //
-           //
+ //           // create a new object for this event with track-related variables
-           // create a new object for this event with track-related variables
+ //           //
-           //
+ //           new(t[ntrkentry]) ToFTrkVar(*t_tof);
-           new(t[ntrkentry]) ToFTrkVar(*t_tof);
+ //           ntrkentry++;
-           ntrkentry++;
+ //           t_tof->Clear();
-           t_tof->Clear();
+ //           //
-           //
+ //         }; // loop on all the tracks
-         }; // loop on all the tracks
+ //       //
-       //
+ //       this->unpackError = unpackError;
-       this->unpackError = unpackError;
+ //       if ( defcal ){
-       if ( defcal ){
+ //         this->default_calib = 1;
-         this->default_calib = 1;
+ //       } else {
-       } else {
+ //         this->default_calib = 0;
-         this->default_calib = 0;
+ //       };
-       };
+ //};
-  };
+ //  return(0);
+ }
+ bool ToFLevel2::bit(int decimal, char pos){
-   return(0);
+   return( (decimal>>pos)%2 );
+ }
+ bool ToFLevel2::checkPMT(TString givenpmt){
+   TClonesArray* Pmt = this->PMT;
+   //  printf(" ou %s entries %i \n",givenpmt.Data(),Pmt->GetEntries());
+   for(int i=0; i<Pmt->GetEntries(); i++) {
+     ToFPMT* pmthit = (ToFPMT*)Pmt->At(i);
+     TString pmtname = this->GetPMTName(pmthit->pmt_id);
+     //    printf(" name %s \n",pmtname.Data());
+     if ( !strcmp(pmtname.Data(),givenpmt.Data()) )
+       return true;
+   }
+   //  printf(" PMT %s missing \n",givenpmt.Data());
+   return false;
+ }
+ bool ToFLevel2::checkPMTpatternPMThit(TrigLevel2 *trg, int &pmtpattern, int &pmtnosignal){
+   UInt_t *patterntrig = trg->patterntrig;
+   pmtpattern = 0;
+   pmtnosignal = 0;
+   bool good = true;
+   //S3
+   if ( this->bit(patterntrig[2],0) ){ pmtpattern++;  if ( !this->checkPMT("S31_1A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],1) ){ pmtpattern++;  if ( !this->checkPMT("S31_2A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],2) ){ pmtpattern++;  if ( !this->checkPMT("S31_3A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],3) ){ pmtpattern++;  if ( !this->checkPMT("S31_1B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],4) ){ pmtpattern++;  if ( !this->checkPMT("S31_2B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],5) ){ pmtpattern++;  if ( !this->checkPMT("S31_3B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],6) ){ pmtpattern++;  if ( !this->checkPMT("S32_1A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],7) ){ pmtpattern++;  if ( !this->checkPMT("S32_2A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],8) ){ pmtpattern++;  if ( !this->checkPMT("S32_3A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],9) ){ pmtpattern++;  if ( !this->checkPMT("S32_1B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],10) ){ pmtpattern++;  if ( !this->checkPMT("S32_2B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[2],11) ){ pmtpattern++;  if ( !this->checkPMT("S32_3B")){ pmtnosignal++; good = false;}}
+   //S2
+   if ( this->bit(patterntrig[3],0) ){ pmtpattern++;  if ( !this->checkPMT("S21_1A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[3],1) ){ pmtpattern++;  if ( !this->checkPMT("S21_2A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[3],2) ){ pmtpattern++;  if ( !this->checkPMT("S21_1B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[3],3) ){ pmtpattern++;  if ( !this->checkPMT("S21_2B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[3],4) ){ pmtpattern++;  if ( !this->checkPMT("S22_1A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[3],5) ){ pmtpattern++;  if ( !this->checkPMT("S22_2A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[3],6) ){ pmtpattern++;  if ( !this->checkPMT("S22_1B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[3],7) ){ pmtpattern++;  if ( !this->checkPMT("S22_2B")){ pmtnosignal++; good = false;}}
+   //S12
+   if ( this->bit(patterntrig[4],0) ){ pmtpattern++;  if ( !this->checkPMT("S12_1A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],1) ){ pmtpattern++;  if ( !this->checkPMT("S12_2A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],2) ){ pmtpattern++;  if ( !this->checkPMT("S12_3A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],3) ){ pmtpattern++;  if ( !this->checkPMT("S12_4A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],4) ){ pmtpattern++;  if ( !this->checkPMT("S12_5A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],5) ){ pmtpattern++;  if ( !this->checkPMT("S12_6A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],6) ){ pmtpattern++;  if ( !this->checkPMT("S12_1A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],7) ){ pmtpattern++;  if ( !this->checkPMT("S12_2A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],8) ){ pmtpattern++;  if ( !this->checkPMT("S12_3A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],9) ){ pmtpattern++;  if ( !this->checkPMT("S12_4B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],10) ){ pmtpattern++; if ( !this->checkPMT("S12_5B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[4],11) ){ pmtpattern++; if ( !this->checkPMT("S12_6B")){ pmtnosignal++; good = false;}}
+   //S11
+   if ( this->bit(patterntrig[5],0) ){ pmtpattern++;  if ( !this->checkPMT("S11_1A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],1) ){ pmtpattern++;  if ( !this->checkPMT("S11_2A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],2) ){ pmtpattern++;  if ( !this->checkPMT("S11_3A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],3) ){ pmtpattern++;  if ( !this->checkPMT("S11_4A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],4) ){ pmtpattern++;  if ( !this->checkPMT("S11_5A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],5) ){ pmtpattern++;  if ( !this->checkPMT("S11_6A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],6) ){ pmtpattern++;  if ( !this->checkPMT("S11_7A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],7) ){ pmtpattern++;  if ( !this->checkPMT("S11_8A")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],8) ){ pmtpattern++;  if ( !this->checkPMT("S11_1B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],9) ){ pmtpattern++;  if ( !this->checkPMT("S11_2B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],10) ){ pmtpattern++; if ( !this->checkPMT("S11_3B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],11) ){ pmtpattern++; if ( !this->checkPMT("S11_4B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],12) ){ pmtpattern++; if ( !this->checkPMT("S11_5B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],13) ){ pmtpattern++; if ( !this->checkPMT("S11_6B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],14) ){ pmtpattern++; if ( !this->checkPMT("S11_7B")){ pmtnosignal++; good = false;}}
+   if ( this->bit(patterntrig[5],15) ){ pmtpattern++; if ( !this->checkPMT("S11_8B")){ pmtnosignal++; good = false;}}
+   return good;
+ }
+ bool ToFLevel2::checkPMTpmttrig(TrigLevel2 *trg){
+   //  UInt_t *patterntrig = trg->patterntrig;
+   int rS11 = 0;
+   int rS12 = 0;
+   int rS21 = 0;
+   int rS22 = 0;
+   int rS31 = 0;
+   int rS32 = 0;
+   // trigger configuration for the event from saved pmts
+   TClonesArray* Pmt = this->PMT;
+   for(int i=0; i<Pmt->GetEntries(); i++) {
+     ToFPMT* pmthit = (ToFPMT*)Pmt->At(i);
+     TString pmtname = this->GetPMTName(pmthit->pmt_id);
+     if ( pmtname.Contains("S11") ) rS11++;
+     if ( pmtname.Contains("S12") ) rS12++;
+     if ( pmtname.Contains("S21") ) rS21++;
+     if ( pmtname.Contains("S22") ) rS22++;
+     if ( pmtname.Contains("S31") ) rS31++;
+     if ( pmtname.Contains("S32") ) rS32++;
+   }
+   int rTOF1 = (rS11 + rS12) * (rS21 + rS22) * (rS31 + rS32);
+   int rTOF2 = (rS11 * rS12) * (rS21 * rS22) * (rS31 * rS32);
+   int rTOF3 = (rS21 + rS22) * (rS31 + rS32);
+   int rTOF4 = (rS21 * rS22) * (rS31 * rS32);
+   int rTOF5 = rS12 * (rS21 * rS22);
+   int rTOF6 = (rS11 + rS12) * (rS31 + rS32);
+   int rTOF7 = (rS11 * rS12) * (rS31 * rS32);
+   // trigger configuration of the run
+   bool TCTOF1 = false;
+   bool TCTOF2 = false;
+   bool TCTOF3 = false;
+   bool TCTOF4 = false;
+   bool TCTOF5 = false;
+   bool TCTOF6 = false;
+   bool TCTOF7 = false;
+   if ( trg->trigconf & (1<<0) ) TCTOF1 = true;
+   if ( trg->trigconf & (1<<1) ) TCTOF2 = true;
+   if ( trg->trigconf & (1<<2) ) TCTOF3 = true;
+   if ( trg->trigconf & (1<<3) ) TCTOF4 = true;
+   if ( trg->trigconf & (1<<4) ) TCTOF5 = true;
+   if ( trg->trigconf & (1<<5) ) TCTOF6 = true;
+   if ( trg->trigconf & (1<<6) ) TCTOF7 = true;
+   // do patterntrig pmts match the trigger configuration?
+   bool pmtsconf_trigconf_match = true;
+   if ( rTOF1 == 0 && TCTOF1 ) pmtsconf_trigconf_match = false;
+   if ( rTOF2 == 0 && TCTOF2 ) pmtsconf_trigconf_match = false;
+   if ( rTOF3 == 0 && TCTOF3 ) pmtsconf_trigconf_match = false;
+   if ( rTOF4 == 0 && TCTOF4 ) pmtsconf_trigconf_match = false;
+   if ( rTOF5 == 0 && TCTOF5 ) pmtsconf_trigconf_match = false;
+   if ( rTOF6 == 0 && TCTOF6 ) pmtsconf_trigconf_match = false;
+   if ( rTOF7 == 0 && TCTOF7 ) pmtsconf_trigconf_match = false;
+   return pmtsconf_trigconf_match;
+ }
+ void ToFLevel2::printPMT(){
+   TClonesArray* Pmt = this->PMT;
+   for(int i=0; i<Pmt->GetEntries(); i++) {
+     ToFPMT* pmthit = (ToFPMT*)Pmt->At(i);
+     TString pmtname = this->GetPMTName(pmthit->pmt_id);
+     printf(" PMT hit: %s \n",pmtname.Data());
+   }
  }
  ToFdEdx::ToFdEdx()
  {
    memset(conn,0,12*sizeof(Bool_t));
    memset(ts,0,12*sizeof(UInt_t));
    memset(te,0,12*sizeof(UInt_t));
+   eDEDXpmt = new TArrayF(48);
    Define_PMTsat();
    Clear();
  }
+ ToFdEdx::~ToFdEdx(){
+   Clear();
+   Delete();
+ }
+ void ToFdEdx::Delete(Option_t *option){
+   if ( eDEDXpmt ){
+     eDEDXpmt->Set(0);
+     if ( eDEDXpmt) delete eDEDXpmt;
+   }
+ }
  //------------------------------------------------------------------------
  void ToFdEdx::CheckConnectors(UInt_t atime, GL_PARAM *glparam, TSQLServer *dbc)
  {
-Line 1493 
 void ToFdEdx::Clear(Option_t *option)
+Line 2042 
 void ToFdEdx::Clear(Option_t *option)
  {
    //
    // Set arrays and initialize structure
+   //  eDEDXpmt.Set(48);    eDEDXpmt.Reset(-1);   // Set array size  and reset structure
-   eDEDXpmt.Set(48);    eDEDXpmt.Reset(-1);   // Set array size  and reset structure
+   eDEDXpmt->Set(48);    eDEDXpmt->Reset(-1);   // Set array size  and reset structure
-   eZpmt.Set(48);       eZpmt.Reset(-1);
-   eDEDXpad.Set(24);    eDEDXpad.Reset(-1);
-   eZpad.Set(24);       eZpad.Reset(-1);
-   eDEDXlayer.Set(6);   eDEDXlayer.Reset(-1);
-   eZlayer.Set(6);      eZlayer.Reset(-1);
-   eDEDXplane.Set(3);   eDEDXplane.Reset(-1);
-   eZplane.Set(3);      eZplane.Reset(-1);
-   INFOpmt.Set(48);     INFOpmt.Reset(0);
-   INFOlayer.Set(6);    INFOlayer.Reset(0);
    //
  };
-Line 1515 
 void ToFdEdx::Print(Option_t *option)
+Line 2055 
 void ToFdEdx::Print(Option_t *option)
  };
  //------------------------------------------------------------------------
- // void ToFdEdx::InitPar(TString parname, TString parfile)
+ void ToFdEdx::Init(pamela::tof::TofEvent *tofl0)
- // {
+ {
- //   // expensive function - call it once/run
+   //
+   ToFLevel2 tf;
+   for (Int_t gg=0; gg<4;gg++){
+     for (Int_t hh=0; hh<12;hh++){
- //   ReadParAtt(            Form("%s/attenuation.txt"              , pardir) );
+       //          tofinput_.tdc[hh][gg]=tofEvent->tdc[gg][hh];
- //   ReadParPos(            Form("%s/desaturation_position.txt"    , pardir) );
+       int mm = tf.GetPMTid(gg,hh);
- //   ReadParBBneg(          Form("%s/BetheBloch.txt"               , pardir) );
+       adc[mm]= (0xFFF & tofl0->adc[gg][hh]); // EM, exclude warning bits
- //   ReadParBBpos(          Form("%s/BetheBloch_betagt1.txt"       , pardir) );
+     };
- //   ReadParDesatBB(        Form("%s/desaturation_beta.txt"        , pardir) );
+   };
- // };
+ };
  //------------------------------------------------------------------------
- void ToFdEdx::Process(UInt_t atime, Float_t betamean, Float_t *xtr_tof, Float_t *ytr_tof, pamela::tof::TofEvent *tofl0 )
+ void ToFdEdx::Init(Int_t gg, Int_t hh, Float_t adce)
+ {
+   //
+   ToFLevel2 tf;
+   //  for (Int_t gg=0; gg<4;gg++){
+   //    for (Int_t hh=0; hh<12;hh++){
+   int mm = tf.GetPMTid(gg,hh);
+   adc[mm]=adce;
+ };
+ //------------------------------------------------------------------------
+ void ToFdEdx::Process(UInt_t atime, Float_t betamean, Float_t *xtr_tof, Float_t *ytr_tof, Int_t exitat)
  {
+   bool debug = false;
+   if ( debug ) printf(" INSIDE TOFDEDX PROCESS \n");
    // the parameters should be already initialised by InitPar()
+   //  printf(" in process \n");
    Clear();
-   //  Float_t betamean = fabs(trackTRK->GetToFTrack()->beta[12]);
-   if(betamean<0.05 || betamean>2){
-     for(int i=0;i<48;i++)INFOpmt[i]=1;
-   }
   // define angle:
    double dx   = xtr_tof[1] - xtr_tof[5];
    double dy   = ytr_tof[0] - ytr_tof[4];
    double dr   = sqrt(dx*dx+dy*dy);
    double theta=atan(dr/76.81);
+   //
+   if ( xtr_tof[1] > 99. ||  xtr_tof[5] > 99. || ytr_tof[0] > 99. ||  ytr_tof[4] > 99. ) theta = 0.;
+   for (Int_t ii=0; ii<6; ii++){
-   //  TArrayF adc;
+     if ( xtr_tof[ii] > 99. ) xtr_tof[ii] = 0.;
-   Float_t adc[48];
+     if ( ytr_tof[ii] > 99. ) ytr_tof[ii] = 0.;
-   ToFLevel2 tf;
-   for (Int_t gg=0; gg<4;gg++){
-     for (Int_t hh=0; hh<12;hh++){
-       //          tofinput_.tdc[hh][gg]=tofEvent->tdc[gg][hh];
-       int mm = tf.GetPMTid(gg,hh);
-       adc[mm]=tofl0->adc[gg][hh];
-     };
    };
+   //
+   if ( debug ) printf(" theta %f \n",theta);
+   if ( debug ) printf(" xtr_tof %.1f %.1f %.1f %.1f %.1f %.1f \n",xtr_tof[0],xtr_tof[1],xtr_tof[2],xtr_tof[3],xtr_tof[4],xtr_tof[5]);
+   if ( debug ) printf(" ytr_tof %.1f %.1f %.1f %.1f %.1f %.1f \n",ytr_tof[0],ytr_tof[1],ytr_tof[2],ytr_tof[3],ytr_tof[4],ytr_tof[5]);
+   //--------------------- TABLE OF PERIODS WITH HV PROBLEMS ----------------------------
+   int Aconn=conn[0];    // PMT 0,20,22,24
+   int Bconn=conn[1];    // PMT 6,12,26,34
+   int Cconn=conn[2];    // PMT 4,14,28,32
+   int Dconn=conn[3];    // PMT 2,8,10,30
+   int Econn=conn[4];    // PMT 42,43,44,47
+   int Fconn=conn[5];    // PMT 7,19,23,27
+   int Gconn=conn[6];    // PMT 3,11,25,33
+   int Hconn=conn[7];    // PMT 1,9,13,21
+   int Iconn=conn[8];    // PMT 5,29,31,35
+   int Lconn=conn[9];    // PMT 37,40,45,46
+   int Mconn=conn[10];    // PMT 15,16,17,18
+   int Nconn=conn[11];    // PMT 36,38,39,41
+   if( false ) cout << Gconn << Iconn << Lconn <<endl; // to avoid compilation warnings
+   //  printf(" size %i \n",eDEDXpmt.GetSize());
    for( int ii=0; ii<48; ii++ ) {
-     if( adc[ii] >= PMTsat[ii]-5 )  continue;
+     //
-     if( adc[ii] <= 0. )            continue;
+     //    eDEDXpmt.SetAt(-1.,ii);
+     //    printf(" ii %i beta %f atime %u xtr 1 %f ytr 1 %f adc %f \n",ii,betamean,atime,xtr_tof[0],ytr_tof[0],adc[ii]);
-     double adcpC   = f_adcPC( adc[ii] );    // - adc conversion in pC
+     if ( debug ) printf("II %i adc %f \n",ii,adc[ii]);
-     double adccorr = adcpC*fabs(cos(theta));
+     if( adc[ii] >= 4095. ){
-          if(adccorr<=0.)           continue;
+       //      eDEDXpmt[ii] = 0.;
+       eDEDXpmt->AddAt(0.,ii);
+       if ( debug ) printf(" %i adc>4095 \n",ii);
+       continue; // EMILIANO
+     };
-     //--------------------- TABLE OF PERIODS WITH HV PROBLEMS ----------------------------
+     if( adc[ii] >= (PMTsat[ii]-5.) && adc[ii] < 4095. ){
+       eDEDXpmt->AddAt(1000.,ii);
+       if ( debug ) printf(" %i adc> pmtsat && adc<4095 \n",ii);
+       continue; // EMILIANO
+     };
-     int Aconn=conn[0];    // PMT 0,20,22,24
+     if( adc[ii] <= 0. ) {
-     int Bconn=conn[1];    // PMT 6,12,26,34
+       eDEDXpmt->AddAt(1500.,ii);
-     int Cconn=conn[2];    // PMT 4,14,28,32
+       if ( debug ) printf(" %i adc<=0 \n",ii);
-     int Dconn=conn[3];    // PMT 2,8,10,30
+       continue;
-     int Econn=conn[4];    // PMT 42,43,44,47
+     };
-     int Fconn=conn[5];    // PMT 7,19,23,27
+     //
-     int Gconn=conn[6];    // PMT 3,11,25,33
+     double adcpC   = f_adcPC( adc[ii] );    // - adc conversion in pC
-     int Hconn=conn[7];    // PMT 1,9,13,21
+     if ( exitat == 0 ){
-     int Iconn=conn[8];    // PMT 5,29,31,35
+       eDEDXpmt->AddAt((Float_t)adcpC,ii);
-     int Lconn=conn[9];    // PMT 37,40,45,46
+       continue;
-     int Mconn=conn[10];    // PMT 15,16,17,18
+     }
-     int Nconn=conn[11];    // PMT 36,38,39,41
+     //    printf(" e qua? \n");
+     double adccorr = adcpC*fabs(cos(theta));
+     if ( debug ) printf(" adccorr %f \n",adccorr);
+     if(adccorr<=0.){
+       if ( debug ) printf(" %i adccorr<=0 \n",ii);
+       //      eDEDXpmt->AddAt((Float_t)adcpC,ii);//?
+       continue;
+     }
+     if ( exitat == 1 ){
+       eDEDXpmt->AddAt((Float_t)adccorr,ii);
+       continue;
+     }
+     //    printf(" e quo? \n");
      //    int standard=0;
-     if( false ) cout << Gconn << Iconn << Lconn <<endl;
      int S115B_ok=0;
      int S115B_break=0;
- //   if(atime>=1153660001 && atime<=1154375000)Dconn=1;
- //     else if(atime>=1155850001 && atime<=1156280000){
- //       Hconn=1;
- //       Nconn=1;
- //     }
- //  else if(atime>=1168490001 && atime<=1168940000)Dconn=1;
- //     else if(atime>=1168940001 && atime<=1169580000){
- //       Fconn=1;
- //       Mconn=1;
- //     }
- //  else if(atime>=1174665001 && atime<=1175000000)Bconn=1;
- //     else if(atime>=1176120001 && atime<=1176800000)Hconn=1;
- //     else if(atime>=1176800001 && atime<=1178330000)Econn=1;
- //     else if(atime>=1178330001 && atime<=1181322000)Hconn=1;
- //     else if(atime>=1182100001 && atime<=1183030000)Aconn=1;
- //     else if(atime>=1184000001 && atime<=1184570000)Hconn=1;
- //     else if(atime>=1185090001 && atime<=1185212000)Dconn=1;
- //     else if(atime>=1191100001 && atime<=1191940000)Dconn=1;
- //     else if(atime>=1196230001 && atime<=1196280000)Hconn=1;
- //     else if(atime>=1206100001 && atime<=1206375600)Cconn=1;
- //     else if(atime>=1217989201 && atime<=1218547800)Econn=1;
- //     else if(atime>=1225789201 && atime<=1226566800)Econn=1;
- //     else if(atime>=1229400901 && atime<=1229700000)Econn=1;
- //     else if(atime>=1230318001 && atime<=1230415200)Econn=1;
- //     else {
- //       standard=1;
- //     }
      if(atime<1158720000)S115B_ok=1;
      else S115B_break=1;
-  //------------------------------------------------------------------------
+     //------------------------------------------------------------------------
+     //    printf(" e qui? \n");
- //---------------------------------------------------- Z reconstruction
+     //---------------------------------------------------- Z reconstruction
- double adcHe, adcnorm, adclin, dEdx, Zeta;
+     double adcHe, adcnorm, adclin, dEdx;//, Zeta; // EM GCC4.7
-  adcHe=-2;
+     adcHe=-2;
-  adcnorm=-2;
+     adcnorm=-2;
-  adclin=-2;
+     adclin=-2;
-  dEdx=-2;
+     dEdx=-2;
-  Zeta=-2;
+     //    Zeta=-2;//EM GCC4.7
+     Double_t correction = 1.;
- //  float ZetaH=-2;
- //  float dEdxH=-2;
- //  double day = (atime-1150000000)/84600;
      if(Aconn==1 && (ii==0 || ii==20 || ii==22 || ii==24)){
-        adcHe   = (Get_adc_he(ii, xtr_tof, ytr_tof))/1.675;
+       correction = 1.675;
      }
      else if(Bconn==1 && (ii==6 || ii==12 || ii==26 || ii==34)){
-        adcHe   = (Get_adc_he(ii, xtr_tof, ytr_tof))/2.482;
+       correction = 2.482;
      }
      else if(Cconn==1 && (ii==4 || ii==14 || ii==28 || ii==32)){
-       adcHe   = (Get_adc_he(ii, xtr_tof, ytr_tof))/1.464;
+       correction = 1.464;
      }
      else if(Dconn==1 && (ii==2 || ii==8 || ii==10 || ii==30)){
-        adcHe   = (Get_adc_he(ii, xtr_tof, ytr_tof))/1.995;
+       correction = 1.995;
      }
      else if(Econn==1 && (ii==42 || ii==43 || ii==44 || ii==47)){
-        adcHe   = (Get_adc_he(ii, xtr_tof, ytr_tof))/1.273;
+       correction = 1.273;
      }
      else if(Fconn==1 && (ii==7 || ii==19 || ii==23 || ii==27)){
-        adcHe   = (Get_adc_he(ii, xtr_tof, ytr_tof))/1.565;
+       correction = 1.565;
      }
      else if(Mconn==1 && (ii==15 || ii==16 || ii==17 || ii==18)){
-        adcHe   = (Get_adc_he(ii, xtr_tof, ytr_tof))/1.565;
+       correction = 1.565;
      }
      else if(Nconn==1 && (ii==36 || ii==38 || ii==39 || ii==41)){
-        adcHe   = (Get_adc_he(ii, xtr_tof, ytr_tof))/1.018;
+       correction = 1.018;
      }
      else if(Hconn==1 && (ii==1 || ii==13 || ii==21 || (ii==9&&S115B_ok==1))){
-       adcHe   = (Get_adc_he(ii, xtr_tof, ytr_tof))/1.84;
+       correction = 1.84;
-     }
-     else if(S115B_break==1 && ii==9 && Hconn==0){
-        adcHe   = f_att5B( ytr_tof[0] );   //N.B.: this function refers to the Carbon!!!
      }
      else if(S115B_break==1 && ii==9 && Hconn==1){
-        adcHe   = (f_att5B( ytr_tof[0] ))/1.64;
+       correction = 1.64;
+     }
+     else correction = 1.;
+     if( ii==9 && S115B_break==1 ){
+       adcHe   = f_att5B( ytr_tof[0] )/correction;
+     } else {
+       adcHe   = Get_adc_he(ii, xtr_tof, ytr_tof)/correction;
+     };
+     if(adcHe<=0){
+       if ( debug ) printf(" %i adcHe<=0 \n",ii);
+       //      eDEDXpmt->AddAt((Float_t)adccorr,ii); //?
+       continue;
+     }
+     if ( exitat == 2 ){
+       if(ii==9 && S115B_break==1)  eDEDXpmt->AddAt(36.*(Float_t)adccorr/adcHe,ii);
+       else  adclin  = 4.*(Float_t)adccorr/adcHe;
+       continue;
      }
-     else  adcHe   = Get_adc_he(ii, xtr_tof, ytr_tof);
-     if(adcHe<=0)   continue;
      if(ii==9 && S115B_break==1)  adcnorm = f_pos5B(adccorr);
      else adcnorm = f_pos( (parPos[ii]), adccorr);
+     if(adcnorm<=0){
-     if(adcnorm<=0) continue;
+       if ( debug ) printf(" %i adcnorm<=0 \n",ii);
+       //      eDEDXpmt->AddAt((Float_t)adccorr,ii);//?
+       continue;
+     }
+     if ( debug ) printf(" adcnorm %f \n",adcnorm);
      if(ii==9 && S115B_break==1)  adclin  = 36.*adcnorm/adcHe;
      else  adclin  = 4.*adcnorm/adcHe;
+     if ( debug ) printf(" adclin %f \n",adclin);
-     if(adclin<=0)  continue;
+     if(adclin<=0){
+       if ( debug ) printf(" %i adclin<=0 \n",ii);
+       //      eDEDXpmt->AddAt((Float_t)adccorr,ii);//?
+       continue;
+     }
+     if ( exitat == 3 ){
+       if(ii==9 && S115B_break==1)  eDEDXpmt->AddAt((Float_t)adclin,ii);
+       else  eDEDXpmt->AddAt((Float_t)adclin,ii);
+       continue;
+     }
+     //
+     if ( betamean > 99. ){
+       //      eDEDXpmt.AddAt((Float_t)adclin,ii);
+       eDEDXpmt->AddAt((Float_t)adclin,ii);
+       //      printf(" AAPMT IS %i dedx is %f vector is %f \n",ii,adclin,eDEDXpmt[ii]);
+       if ( debug ) printf(" %i betamean > 99 \n",ii);
+       continue;
+     };
+     //
      double dEdxHe=-2;
      if(ii==9 && S115B_break==1){
        if( betamean <1. ) dEdxHe = f_BB5B( betamean );
-Line 1709 
 double adcHe, adcnorm, adclin, dEdx, Zet
+Line 2275 
 double adcHe, adcnorm, adclin, dEdx, Zet
        if( betamean <1. ) dEdxHe = f_BB( (parBBneg[ii]), betamean );
        else                       dEdxHe = parBBpos[ii];
      }
-     if(dEdxHe<=0)  continue;
-     if(ii==9 && S115B_break==1)  dEdx = f_desatBB5B( adclin );
-     else  dEdx = f_desatBB((parDesatBB[ii]), adclin );
-     if(dEdx<=0)    continue;
-     if(ii==9 && S115B_break==1)  Zeta = sqrt(36.*(dEdx/dEdxHe));
-     else  Zeta = sqrt(4.*(dEdx/dEdxHe));
-     if(Zeta<=0)    continue;
- //--------------------- TIME DEPENDENCE ----------------------------------
- //      TArrayF &binx = TDx[ii];
- //      TArrayF &biny = TDy[ii];
- //      for(int k=0; k<200; k++) {
- //        if (day > binx[k]-2.5 && day<=binx[k]+2.5 && biny[k]>0)  {
- //       ZetaH=Zeta/biny[k]*6;
- //       dEdxH=dEdx/(pow(biny[k],2))*36;
- //        }
- //      }
- //      if(ZetaH!=-2)eZpmt[ii]=(Float_t)ZetaH;
- //      else eZpmt[ii]=(Float_t)Zeta;
- //      if(dEdxH!=-2)eDEDXpmt[ii]=(Float_t)dEdxH;
- //      else eDEDXpmt[ii]=(Float_t)dEdx;
- //      printf("%5d %8.2f %8.2f %8.2f  %8.2f %8.2f  %8.2f %5.4f \n",               ii, adcpC,  adccorr, adcHe, dEdxHe, dEdx, Zeta, betamean );
-     eZpmt[ii]=(Float_t)Zeta;
-     eDEDXpmt[ii]=(Float_t)dEdx;
-  }  //end loop on 48 PMT
- //---------------------------------------------------  paddle + layer --------------------
-   for(int j=0;j<48;j++){
-     int k=100;
-     if(j%2==0 || j==0)k=j/2;
-     double zpdl=-1;
-     if((j%2==0 || j==0) && eZpmt[j]!=-1 && eZpmt[j+1]!=-1){
-       zpdl=0.5*(eZpmt[j]+eZpmt[j+1]);
-     }else if((j%2==0 || j==0) && eZpmt[j]!=-1 && eZpmt[j+1]==-1){
-       zpdl=eZpmt[j];
-     }else if((j%2==0 || j==0) && eZpmt[j]==-1 && eZpmt[j+1]!=-1){
-       zpdl=eZpmt[j+1];
-     }
-     if(j%2==0 || j==0)eZpad[k]= (Float_t)zpdl;
-     if((j%2==0 || j==0)&&eZpad[k]!=-1){
-       if(k>=0&&k<8)eZlayer[0]=eZpad[k];
-       if(k>=8&&k<14)eZlayer[1]=eZpad[k];
-       if(k>=14&&k<16)eZlayer[2]=eZpad[k];
-       if(k>=16&&k<18)eZlayer[3]=eZpad[k];
-       if(k>=18&&k<21)eZlayer[4]=eZpad[k];
-       if(k>=21)eZlayer[5]=eZpad[k];
-     }
-     if(eZlayer[0]!=-1&&eZlayer[1]!=-1&&fabs(eZlayer[0]-eZlayer[1])<1.5)eZplane[0]=0.5*(eZlayer[0]+eZlayer[1]);
-     else if(eZlayer[0]!=-1&&eZlayer[1]==-1)eZplane[0]=eZlayer[0];
-     else if(eZlayer[1]!=-1&&eZlayer[0]==-1)eZplane[0]=eZlayer[1];
-     if(eZlayer[2]!=-1&&eZlayer[3]!=-1&&fabs(eZlayer[2]-eZlayer[3])<1.5)eZplane[1]=0.5*(eZlayer[2]+eZlayer[3]);
-     else if(eZlayer[2]!=-1&&eZlayer[3]==-1)eZplane[1]=eZlayer[2];
-     else if(eZlayer[3]!=-1&&eZlayer[2]==-1)eZplane[1]=eZlayer[3];
-     if(eZlayer[4]!=-1&&eZlayer[5]!=-1&&fabs(eZlayer[4]-eZlayer[5])<1.5)eZplane[2]=0.5*(eZlayer[4]+eZlayer[5]);
-     else if(eZlayer[4]!=-1&&eZlayer[5]==-1)eZplane[2]=eZlayer[4];
-     else if(eZlayer[5]!=-1&&eZlayer[4]==-1)eZplane[2]=eZlayer[5];
-   }
-   for(int jj=0;jj<48;jj++){
-     int k=100;
-     if(jj%2==0 || jj==0)k=jj/2;
-     double dedxpdl=-1;
-     if((jj%2==0 || jj==0) && eDEDXpmt[jj]!=-1 && eDEDXpmt[jj+1]!=-1){
-       dedxpdl=0.5*(eDEDXpmt[jj]+eDEDXpmt[jj+1]);
-     }else if((jj%2==0 || jj==0) && eDEDXpmt[jj]!=-1 && eDEDXpmt[jj+1]==-1){
-       dedxpdl=eDEDXpmt[jj];
-     }else if((jj%2==0 || jj==0) && eDEDXpmt[jj]==-1 && eDEDXpmt[jj+1]!=-1){
-       dedxpdl=eDEDXpmt[jj+1];
-     }
-     if(jj%2==0 || jj==0)eDEDXpad[k]= (Float_t)dedxpdl;
+     if ( debug ) printf(" dEdxHe %f \n",dEdxHe);
-     if((jj%2==0 || jj==0)&&eDEDXpad[k]!=-1){
+     if(dEdxHe<=0){
-       if(k>=0&&k<8)eDEDXlayer[0]=eDEDXpad[k];
+       eDEDXpmt->AddAt((Float_t)adclin,ii);
-       if(k>=8&&k<14)eDEDXlayer[1]=eDEDXpad[k];
+       if ( debug ) printf(" %i dEdxHe<=0 \n",ii);
-       if(k>=14&&k<16)eDEDXlayer[2]=eDEDXpad[k];
+       continue;
-       if(k>=16&&k<18)eDEDXlayer[3]=eDEDXpad[k];
+     };
-       if(k>=18&&k<21)eDEDXlayer[4]=eDEDXpad[k];
-       if(k>=21)eDEDXlayer[5]=eDEDXpad[k];
-     }
-     if(eDEDXlayer[0]!=-1&&eDEDXlayer[1]!=-1&&fabs(eDEDXlayer[0]-eDEDXlayer[1])<10)eDEDXplane[0]=0.5*(eDEDXlayer[0]+eDEDXlayer[1]);
+     if(ii==9 && S115B_break==1)  dEdx = f_desatBB5B( adclin );
-     else if(eDEDXlayer[0]!=-1&&eDEDXlayer[1]==-1)eDEDXplane[0]=eDEDXlayer[0];
+     else  dEdx = f_desatBB((parDesatBB[ii]), adclin );
-     else if(eDEDXlayer[1]!=-1&&eDEDXlayer[0]==-1)eDEDXplane[0]=eDEDXlayer[1];
-     if(eDEDXlayer[2]!=-1&&eDEDXlayer[3]!=-1&&fabs(eDEDXlayer[2]-eDEDXlayer[3])<10)eDEDXplane[1]=0.5*(eDEDXlayer[2]+eDEDXlayer[3]);
+     if(dEdx<=0){
-     else if(eDEDXlayer[2]!=-1&&eDEDXlayer[3]==-1)eDEDXplane[1]=eDEDXlayer[2];
+       eDEDXpmt->AddAt((Float_t)adclin,ii);
-     else if(eDEDXlayer[3]!=-1&&eDEDXlayer[2]==-1)eDEDXplane[1]=eDEDXlayer[3];
+       if ( debug ) printf(" %i dEdx<=0 \n",ii);
+       continue;
+     };
-     if(eDEDXlayer[4]!=-1&&eDEDXlayer[5]!=-1&&fabs(eDEDXlayer[4]-eDEDXlayer[5])<10)eDEDXplane[2]=0.5*(eDEDXlayer[4]+eDEDXlayer[5]);
+     if ( debug ) printf(" dEdx %f \n",dEdx);
-     else if(eDEDXlayer[4]!=-1&&eDEDXlayer[5]==-1)eDEDXplane[2]=eDEDXlayer[4];
+     eDEDXpmt->AddAt((Float_t)dEdx,ii);
-     else if(eDEDXlayer[5]!=-1&&eDEDXlayer[4]==-1)eDEDXplane[2]=eDEDXlayer[5];
+     //    eDEDXpmt.AddAt((Float_t)dEdx,ii);
-   }
+     //    printf(" PMT IS %i dedx is %f vector is %f \n",ii,dEdx,eDEDXpmt[ii]);
+   }  //end loop on 48 PMT
  };
  //------------------------------------------------------------------------
- void ToFdEdx::PrintTD()
- {
-   for(int i=0; i<48; i++) {
-     TArrayF &binx = TDx[i];
-     TArrayF &biny = TDy[i];
-     for(int k=0; k<200; k++) {  // bin temporali
-       printf("%d %d %f %f", i,k, binx[k], biny[k]);
-     }
-   }
- }
- //------------------------------------------------------------------------
  void ToFdEdx::Define_PMTsat()
  {
    Float_t  sat[48] = {
-Line 1861 
 void ToFdEdx::Define_PMTsat()
+Line 2318 
 void ToFdEdx::Define_PMTsat()
  }
  //------------------------------------------------------------------------
- // void ToFdEdx::ReadParTD( Int_t ipmt, const char *fname )
- // {
- //   printf("read %s\n",fname);
- //   if(ipmt<0)  return;
- //   if(ipmt>47) return;
- //   FILE *fattin = fopen( fname , "r" );
- //   Float_t yTD[200],xTD[200];
- //   for(int j=0;j<200;j++){
- //     float x,y,ym,e;
- //     if(fscanf(fattin,"%f %f %f %f",
- //            &x, &y, &ym, &e )!=4) break;
- //     xTD[j]=x;
- //     if(ym>0&&fabs(y-ym)>1)  yTD[j]=ym;
- //     else                    yTD[j]=y;
- //   }
- //   TDx[ipmt].Set(200,xTD);
- //   TDy[ipmt].Set(200,yTD);
- //   fclose(fattin);
- // }
- //------------------------------------------------------------------------
  void ToFdEdx::ReadParBBpos( const char *fname )
  {
-   printf("read %s\n",fname);
+   //  printf("read %s\n",fname);
    parBBpos.Set(48);
    FILE *fattin = fopen( fname , "r" );
    for (int i=0; i<48; i++) {
-Line 1900 
 void ToFdEdx::ReadParBBpos( const char *
+Line 2336 
 void ToFdEdx::ReadParBBpos( const char *
  //------------------------------------------------------------------------
  void ToFdEdx::ReadParDesatBB( const char *fname )
  {
-   printf("read %s\n",fname);
+   //  printf("read %s\n",fname);
    FILE *fattin = fopen( fname , "r" );
    for (int i=0; i<48; i++) {
      int   tid=0;
-Line 1917 
 void ToFdEdx::ReadParDesatBB( const char
+Line 2353 
 void ToFdEdx::ReadParDesatBB( const char
  void ToFdEdx::ReadParBBneg( const char *fname )
  {
-   printf("read %s\n",fname);
+   //  printf("read %s\n",fname);
    FILE *fattin = fopen( fname , "r" );
    for (int i=0; i<48; i++) {
      int   tid=0;
-Line 1932 
 void ToFdEdx::ReadParBBneg( const char *
+Line 2368 
 void ToFdEdx::ReadParBBneg( const char *
  //------------------------------------------------------------------------
  void ToFdEdx::ReadParPos( const char *fname )
  {
-   printf("read %s\n",fname);
+   //  printf("read %s\n",fname);
    FILE *fattin = fopen( fname , "r" );
    for (int i=0; i<48; i++) {
      int   tid=0;
-Line 1947 
 void ToFdEdx::ReadParPos( const char *fn
+Line 2383 
 void ToFdEdx::ReadParPos( const char *fn
  //------------------------------------------------------------------------
  void ToFdEdx::ReadParAtt( const char *fname )
  {
-   printf("read %s\n",fname);
+   //  printf("read %s\n",fname);
    FILE *fattin = fopen( fname , "r" );
    for (int i=0; i<48; i++) {
      int   tid=0;
-Line 2020 
 float ToFdEdx::Get_adc_he( int id, float
+Line 2456 
 float ToFdEdx::Get_adc_he( int id, float
    //             pl_x - coord x of the tof plane
    //             pl_y - coord y
-    adc_he = 0;
+   adc_he = 0;
    if( eGeom.GetXY(id)==1 )  adc_he = f_att( (parAtt[id]), pl_x[eGeom.GetPlane(id)] );
    if( eGeom.GetXY(id)==2 )  adc_he = f_att( (parAtt[id]), pl_y[eGeom.GetPlane(id)] );
    return adc_he;
-Line 2055 
 double ToFdEdx::f_desatBB5B( float x )
+Line 2491 
 double ToFdEdx::f_desatBB5B( float x )
 .009*x*x;
  }

 Legend:



Removed from v.1.26
 


changed lines


 
Added in v.1.42
 Legend:



Removed from v.1.26
 


changed lines


 
Added in v.1.42
-Removed from v.1.26
+Added in v.1.42

	ViewVC Help
Powered by ViewVC 1.1.23