TTools.C

Go to the documentation of this file.

// Class TTools
//
// Author     : F. Cassol Brunner
// Created    : 12/9/2002
// Last update: $Date: 2005/06/28 21:19:21 $
//          by: $Author: sauvan $
// Comment: class for general tool routines

#include "Marana/TTools.h"

#include "H1Mods/H1PartSelTrackArrayPtr.h"
#include "H1Mods/H1PartSelTrack.h"


// Import class
ClassImp(TTools);

TTools::TTools()
{
  GenPl = 27.598;
  GenPp = 919.971;  
  RunNumber=0;
  EventNumber=0;
  NewRun=kTRUE;
  MCFlag=0; 
  num=0;
  numMax=10000000;
}

TTools::TTools(Int_t Run, Int_t Event,Int_t MCFlag)
{
  GenPl = 27.598;
  GenPp = 919.971;  
  RunNumber=Run;
  EventNumber=Event;
  NewRun=kTRUE;
  MCFlag=MCFlag; 
  num=0;
  numMax=10000000;
}

Bool_t TTools::Next()
{
  if(num==numMax)return false;
  if(!gH1Tree->Next())return false;

  num++;
//...re-set run and event number
 if(RunNumber!=gH1Tree->GetRunNumber()) NewRun=kTRUE;
 else NewRun=kFALSE;

 RunNumber=gH1Tree->GetRunNumber();
 EventNumber=gH1Tree->GetEventNumber();

 static H1ShortPtr RunType("RunType");
 MCFlag=*RunType;


 return true;
}




void TTools::SetEventList(Char_t *ListFileName)
{
//     set an eventlist to the tree
  cout << "Set  eventlist from " << ListFileName <<endl;
  TFile *file=NULL;
  file = TFile::Open(ListFileName);
  TList *keys = file->GetListOfKeys();
  TIter next (keys);
  Bool_t haslist=kFALSE;
  
  while ( TKey *thekey = (TKey*)next() ) {
      if ( strcmp(thekey->GetClassName(),"H1EventList") == 0 ) {
          H1EventList *thelist = (H1EventList*)file->Get( thekey->GetName() );
          gH1Tree->SetEventList(thelist);
          haslist=kTRUE;
          break;
      }
  }

  if (! haslist ) {
      cerr << "Found no H1EventList on file! "<<endl;
  }

}



Double_t TTools::GetMod(Double_t x, Double_t y)
{
// Return modulo 
  return((x - ((Int_t) (x/y)) * y));
}

Double_t TTools::GetRealPhi(Double_t PhiIn)
{
// Get phi range of input angle PhiIn

  Double_t phi = PhiIn;
  if(phi <= TMath::Pi()){
    if(phi > -TMath::Pi()){ }//      return phi;  } 
    else if(phi > -3*TMath::Pi()){ phi += 2*TMath::Pi(); } 
    else { phi = -GetMod(TMath::Pi()-phi, 2*TMath::Pi()) + TMath::Pi(); }
  }
  else if(phi <= 3*TMath::Pi()) { phi -= 2*TMath::Pi(); }
  else { phi = GetMod(TMath::Pi()+phi, 2*TMath::Pi()) - TMath::Pi(); }

  if(fabs(phi) <= 1e-15) phi = 0.;
  return phi;
}

Double_t TTools::GetMergedPhi(Double_t PhiPart1,Double_t PtPart1,Double_t PhiPart2,Double_t PtPart2)
{
// Get phi angle between two particles

  Double_t WeightOfFirst = 0;
  Double_t PhiAngle = 0;  
  
  if (PtPart1>0 && PtPart2) {
    WeightOfFirst = GetRealPhi(PhiPart1-PhiPart2) * PtPart1/(PtPart1+PtPart2);
    PhiAngle = GetRealPhi(WeightOfFirst+PhiPart2);
  }
  return(PhiAngle);
}



Double_t TTools::GetProjPhi(Double_t Phi)
{
//     ProjPhi :
//      phi projection into -22.5-22.5 degrees to see cracks
//      Phi in rad.
Double_t PhiProjection=-99999.;
Double_t pi=3.14159265359;
 
PhiProjection=Phi;
 
        if(fabs(Phi*MYR2D)>22.5&&fabs(Phi*MYR2D)<=67.5) {
                if(Phi*MYR2D>0.) {
                        PhiProjection=Phi-pi/4.;
                } else {
                        PhiProjection=Phi+pi/4.;
                }
        }
        if(fabs(Phi*MYR2D)>67.5&&fabs(Phi*MYR2D)<=112.5) {
                if(Phi*MYR2D>0.) {
                        PhiProjection=Phi-pi/2.;
                } else {
                        PhiProjection=Phi+pi/2.;
                }
        }
        if(fabs(Phi*MYR2D)>112.5&&fabs(Phi*MYR2D)<=157.5) {
                        if(Phi*MYR2D>0.) {
                        PhiProjection=Phi-3.*pi/4.;
                } else {
                        PhiProjection=Phi+3.*pi/4.;
                }
        }
        if(fabs(Phi*MYR2D)>157.5&&fabs(Phi*MYR2D)<180.) {
//                 if(Phi*MYR2D>0.) {
//                         PhiProjection=Phi-9.*pi/8.;
//                 } else {
//                         PhiProjection=Phi+9.*pi/8.;
//                 }
                 if(Phi*MYR2D>0.) {
                        PhiProjection=Phi-pi;
                } else {
                        PhiProjection=Phi+pi;
                }
       }
        if(fabs(Phi*MYR2D)==180.) {
                if(Phi*MYR2D>0.) {
                        PhiProjection=Phi-pi;
                } else {
                        PhiProjection=Phi+pi;
                }
        }
 
return PhiProjection;
}    


Int_t TTools::GetMuonGrade(H1PartMuon *Muon)
{
//    decode and returns muon grade (best (==min.0 muon grade)
 Int_t iMinGr = 1000;
 for(Int_t iGr=1;iGr<=6;iGr++){
   if (Muon->TestGrade(iGr) == kTRUE)
     if (iGr < iMinGr) iMinGr = iGr;
 }
 return iMinGr;
}

Bool_t TTools::IsData2005(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
  if (MCFlag == 0 && RunNumber > 399101)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}

Bool_t TTools::IsData2004e(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
//  if (MCFlag == 0 && RunNumber >= 335423)
  if (MCFlag == 0 && RunNumber <= 398818 && RunNumber >= 396674)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}

Bool_t TTools::IsData0304R1(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
//   if (MCFlag == 0 && RunNumber <= 376562 && RunNumber >= 356241)
  if (MCFlag == 0 && RunNumber <= 376562 && RunNumber >= 357072) //--- duk ranges
    Selection = kTRUE;

  return(Selection);
}
Bool_t TTools::IsData0304L1(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
  if (MCFlag == 0 && RunNumber < 377533 && RunNumber > 376562)
    Selection = kTRUE;

  return(Selection);
}
Bool_t TTools::IsData0304L2(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
//  if (MCFlag == 0 && RunNumber < 387537 && RunNumber >= 377533)
  if (MCFlag == 0 && RunNumber <= 386696 && RunNumber >= 377533)//--- duk ranges
    Selection = kTRUE;

  return(Selection);
}
Bool_t TTools::IsData0304R2(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
//   if (MCFlag == 0 && RunNumber <= 392213 && RunNumber >= 387537)
  if (MCFlag == 0 && RunNumber <= 392213 && RunNumber >= 387942)//--- duk ranges
    Selection = kTRUE;

  return(Selection);
}

Bool_t TTools::IsData2004p(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
//  if (MCFlag == 0 && RunNumber >= 335423)
  if (MCFlag == 0 && RunNumber <= 392213 && RunNumber >= 367257)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsData2003(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
//  if (MCFlag == 0 && RunNumber >= 335423)
  if (MCFlag == 0 && RunNumber > 335000)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsDataHeraII(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
//  if (MCFlag == 0 && RunNumber >= 335423)
  if (MCFlag == 0 && RunNumber > 335000)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsData2002(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
//  if (MCFlag == 0 && RunNumber > 303906 && RunNumber <= 334964)
  if (MCFlag == 0 && RunNumber > 300000 && RunNumber < 335000)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}

Bool_t TTools::IsData2000(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
  if (MCFlag == 0 && RunNumber > 262144 && RunNumber < 279215)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsData1999p(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
  if (MCFlag == 0 && RunNumber >= 244968 && RunNumber <= 261338)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsData1999e(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
  if (MCFlag == 0 && RunNumber >= 231721 && RunNumber <= 241649)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsData1998(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
  if (MCFlag == 0 && RunNumber >= 218217 && RunNumber <= 230164)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsData1997(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
  if (MCFlag == 0 && RunNumber > 176421 && RunNumber <= 201519)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsData1996(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
  if (MCFlag == 0 && RunNumber >= 146416 && RunNumber <= 171156)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsData1995(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
  if (MCFlag == 0 && RunNumber >= 110946 && RunNumber <= 131045)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}
Bool_t TTools::IsData1994(void)
{
// Check Data set using runnumber
 
  Bool_t Selection = kFALSE;
 
  if (MCFlag == 0 && RunNumber > 0 && RunNumber < 110946)
    Selection = kTRUE;

  // Return selection flag
  return(Selection);
}


Bool_t TTools::IsData94(const UInt_t& irun)
{
   if(irun>=76137 && irun<=90419)return kTRUE;
   return kFALSE;
}
Bool_t TTools::IsData95S(const UInt_t& irun)
{
  if(irun>=110946 && irun<=131045)
    {
      if(irun>=130066 && irun<=130385)return kTRUE;
      else if (irun>=129473&& irun<=129933)return kTRUE;
      else return kFALSE;
    }
  return kFALSE;
}
Bool_t TTools::IsData95P(const UInt_t& irun)
{
  if(irun>=110946 && irun<=131045)
    {
      if(irun>=130066 && irun<=130385)return kFALSE;
      else if (irun>=129473&& irun<=129933)return kFALSE;
      else return kTRUE;
    }
  return kFALSE;
}
Bool_t TTools::IsData96P(const UInt_t& irun)
{
  if(irun>=146416 && irun<=171156) return kTRUE;
  return kFALSE;
}
Bool_t TTools::IsData97P(const UInt_t& irun)
{
  if(irun>=177920 && irun<=200407)return kTRUE;
  return kFALSE;
}
Bool_t TTools::IsData97S(const UInt_t& irun)
{
  if(irun>=200445 && irun<=201519)return kTRUE;
  return kFALSE;
}
Bool_t TTools::IsData98E(const UInt_t& irun)
{
  if(irun>=218217 && irun<=230164) return kTRUE;
  return kFALSE;
}

Bool_t TTools::IsData99S(const UInt_t& irun)
{
  if(irun>=259487 && irun<=261338)return kTRUE;
  return kFALSE;
}
Bool_t TTools::IsData99E(const UInt_t& irun)
{
  if(irun>=231721 && irun<=241649) return kTRUE;
  return kFALSE;
}

Bool_t TTools::IsData99P(const UInt_t& irun)
{
  if(irun>=244968 && irun<=261338)
    {
      if(irun>=244968 && irun<=259461) return kTRUE;
      else  if(irun>=259487 && irun<=261338)return kFALSE;
      else {cerr<<"Warning in <RTools::IsData99P>: Irun "<<irun<<" unassignable (year 1999)."<<endl;return kFALSE;}
    }
  return kFALSE;
}
Bool_t TTools::IsData00S(const UInt_t& irun)
{
  if(irun>=278691 && irun<=278978)return kTRUE;
  return kFALSE;
}
Bool_t TTools::IsData00P(const UInt_t& irun)
{
  if(irun>=262204 && irun<=279215)
    {
      if(irun>=278691 && irun<=278978)return kFALSE;
      else return kTRUE;
    }
  return kFALSE;
}








Bool_t TTools::IsHV(void)
{
// Are HV bits on??
  if(MCFlag) return kTRUE;      //--- valid only if DATA

//  if(IsData2005()||IsData2004e()) return IsHVCJCLoose();

   // High voltage bits
   static H1ShortPtr hv1("Mhv1");
   static H1ShortPtr hv2("Mhv2");
   static H1ShortPtr bbl3("Mbbl3");

   if(IsData2003()) {
     if ( (
         (hv1[0]||hv1[1]||hv2[0]||hv2[1]||hv2[2]||bbl3[0]>0) // CJC
         || (hv1[11]||bbl3[6]>0) // LAr
//       || (hv1[18]||hv1[19]||hv1[20]||hv2[17]||hv2[18]||hv2[19]||bbl3[12]>0) // lumi
//       || (hv1[20]||hv2[18]||bbl3[12]>0) // lumi
         || (hv1[20]||hv2[18]) // lumi
         || (hv1[14]||bbl3[9]>0) // ToF
         || bbl3[8]>0 // Spacal
         ))
         return(kFALSE);
      else 
         return(kTRUE); 
   }else{
     if ( (
         (hv1[0]||hv1[1]||hv2[0]||hv2[1]||hv2[2]||bbl3[0]>0) // CJC
         || (hv1[11]||bbl3[6]>0) // LAr
         || (hv1[18]||hv1[19]||hv1[20]||hv2[17]||hv2[18]||hv2[19]||bbl3[12]>0) // lumi
         || (hv1[14]||bbl3[9]>0) // ToF
         || bbl3[8]>0 // Spacal
         ))
         return(kFALSE);
      else 
         return(kTRUE); 
   }  

} 

Bool_t TTools::IsHVCJCLoose(void)
{
// Are HV bits on?? (with CJC1 or CJC2

  if(MCFlag) return kTRUE;      //--- valid only if DATA


   // High voltage bits
   static H1ShortPtr hv1("Mhv1");
   static H1ShortPtr hv2("Mhv2");
   static H1ShortPtr bbl3("Mbbl3");

   if(IsDataHeraII()) {
     if ( (
//          ((hv1[0]||hv2[0]>0)||
//          ((hv1[1]||hv2[1]||hv2[2]>0)&&bbl3[0]>0) // CJC1 or CJC2
         (hv1[1]||hv2[1]||hv2[2]>0) // CJC
         || (hv1[11]||bbl3[6]>0) // LAr
         || (hv1[20]||hv2[18]) // lumi
         || (hv1[14]||bbl3[9]>0) // ToF
         || bbl3[8]>0 // Spacal
         ))
         return(kFALSE);
      else 
         return(kTRUE); 
   }else{
     if ( (
         (hv1[0]||hv1[1]||hv2[0]||hv2[1]||hv2[2]||bbl3[0]>0) // CJC
         || (hv1[11]||bbl3[6]>0) // LAr
         || (hv1[18]||hv1[19]||hv1[20]||hv2[17]||hv2[18]||hv2[19]||bbl3[12]>0) // lumi
         || (hv1[14]||bbl3[9]>0) // ToF
         || bbl3[8]>0 // Spacal
         ))
         return(kFALSE);
      else 
         return(kTRUE); 
   }  

} 


Bool_t TTools::IsBdcOn()
{
  // 0/1 : CIP and COP off/on
// Are HV bits on??
  if(MCFlag) return kTRUE;      //--- valid only if DATA


   // High voltage bits
   static H1ShortPtr hv1("Mhv1");
   static H1ShortPtr hv2("Mhv2");
   static H1ShortPtr bbl3("Mbbl3");

  return hv1[21]||hv1[22]>0 ? 0 : 1;
}

Bool_t TTools::IsFmdOn()
{
  // 0/1 : CIP and COP off/on
// Are HV bits on??
  if(MCFlag) return kTRUE;      //--- valid only if DATA


   // High voltage bits
   static H1ShortPtr hv1("Mhv1");
   static H1ShortPtr hv2("Mhv2");
   static H1ShortPtr bbl3("Mbbl3");

  return hv1[31] || hv2[23]>0 ? 0 : 1;
}

Bool_t TTools::IsPlugOn()
{
  // 0/1 : CIP and COP off/on
// Are HV bits on??
  if(MCFlag) return kTRUE;      //--- valid only if DATA


   // High voltage bits
   static H1ShortPtr hv1("Mhv1");
   static H1ShortPtr hv2("Mhv2");
   static H1ShortPtr bbl3("Mbbl3");

  return hv1[12] || hv2[15] || bbl3[7]>0 ? 0 : 1;
}




Bool_t TTools::IsCipCopOn()
{
  // 0/1 : CIP and COP off/on
// Are HV bits on??
  if(MCFlag) return kTRUE;      //--- valid only if DATA


   // High voltage bits
   static H1ShortPtr hv1("Mhv1");
   static H1ShortPtr hv2("Mhv2");
   static H1ShortPtr bbl3("Mbbl3");

  return hv1[7]||hv1[8]||hv2[11]||hv2[12]||bbl3[14]>0 ? 0 : 1;
}

Bool_t TTools::IsCipOn()
{
  // 0/1 : CIP  off/on
// Are HV bits on??
  if(MCFlag) return kTRUE;      //--- valid only if DATA


   // High voltage bits
   static H1ShortPtr hv1("Mhv1");
   static H1ShortPtr hv2("Mhv2");
   static H1ShortPtr bbl3("Mbbl3");

  return  hv1[7]|| hv2[11]|| bbl3[14]>0 ? 0 : 1;
}





//------------------------------------------------------
Bool_t TTools::IsCjcLooseOn()  
{ 
  // 0/1 : CJC1 and CJC2 off / either CJC1 or CJC2 on
   // High voltage bits
   static H1ShortPtr hv1("Mhv1");
   static H1ShortPtr hv2("Mhv2");
   static H1ShortPtr bbl3("Mbbl3");

  if ((hv1[0]||hv2[0]>0)||
      (hv1[1]||hv2[1]||hv2[2]>0)&&bbl3[0]>0) return 0;
  else return 1;
}
//------------------------------------------------------



  
Bool_t TTools::IsZCrack(Float_t z)
{
// Is z Crack

  // selection is false at beg.
  Bool_t selection = kFALSE;

  // Define crack position in z
  const Int_t nCracks = 7;
//  Float_t zCrack[nCracks] = {-211.14,-151.43,-64.63,23.17,110.12,198.43,289.145};
  Float_t zCrack[nCracks] = {-211.14,-151.43,-64.63,23.17,110.12};
  // Define Min approach in z
  const Float_t ZAppr = 2.0;

  // Check all phi-cracks
  for(Int_t iCrack=0;iCrack<nCracks;iCrack++){
    if (fabs(z-zCrack[iCrack]) < ZAppr){
      selection = kTRUE;
      break;
    }
  }

  // Return selection 
  return(selection);

}

Bool_t TTools::IsPhiCrack(Float_t Phi)
{
// Check if Phi is phi-crack +- 2deg

  // selection is false at beg.
  Bool_t selection = kFALSE;

  // Define crack position in phi
  const Int_t nCracks = 9;
  Float_t PhiCrack[nCracks]={-180.0,-135.0,-90.0,-45.0,0.0,45.0,90.0,135.0,180.0};
  
  // Define Min approach in phi
  const Float_t PhiAppr = 2.0;

  // Phi in degree
  Float_t PhiDeg = Phi*180/3.1415;

  // Check all phi-cracks
  for(Int_t iCrack=0;iCrack<nCracks;iCrack++){
    if (fabs(PhiDeg-PhiCrack[iCrack]) < PhiAppr){
      selection = kTRUE;
      break;
    }
  }
  // Return selection 
  return(selection);
}


Bool_t TTools::IsLArCrack(Float_t z,Float_t phi)
{
  // init at false
  Bool_t selection=kFALSE;

  // exclude CB2-CB3 crack
  if(z<-180 || (z<=25 && z>=15) ) selection=kTRUE;

  // put Phi in 0-360 deg range 
  if(phi<0) phi+=2*TMath::Pi();
  phi*=MYR2D;
  Float_t phimod=GetMod(phi,45);

  // exlude cracks in barrrel part + 2 cracks in IF part
  if(   z <=290 && (phimod < 2 || phimod > 43)
        || z > 290 && ( phi<2 || (phi>178 && phi <180) || phi>358  )) selection=kTRUE;


  // Return selection 
  return(selection);
}



Bool_t TTools::IsLArCrack(Float_t z,Float_t zBBE,Float_t phi)
{
// Definition of LAr cracks taken into account BBE structure
// (Benjamin)

  // init at false
  Bool_t selection=kFALSE;

  // exclude CB2-CB3 crack (190)
  if(z<-190 || (z<=25 && z>=15) ){
    //    cout << " z= " << z;
        selection=kTRUE;
  } 
  if (z<=-150){
    if (zBBE <= -190) selection=kTRUE;
  }

 // if(z<=25 && z>=15) selection=kTRUE;

  // put Phi in 0-360 deg range 
  if(phi<0) phi+=2*TMath::Pi();
  phi*=MYR2D;
  Float_t phimod=GetMod(phi,45);

  // exlude cracks in barrrel part + 2 cracks in IF part
  if(   z <=292 && (phimod < 2 || phimod > 43)
        || z > 292 && ( phi<2 || (phi>178 && phi <180) || phi>358  )) selection=kTRUE;


  // Return selection 
  return(selection);
}


Bool_t TTools::IsHVCaloCJC(void)
{
// Are HV bits on??
  if(MCFlag) return kTRUE;      //--- valid only if DATA

   // High voltage bits
   static H1ShortPtr hv1("Mhv1");
   static H1ShortPtr hv2("Mhv2");
   static H1ShortPtr bbl3("Mbbl3");


   if ( (
         (hv1[0]||hv1[1]||hv2[0]||hv2[1]||hv2[2]||bbl3[0]>0) // CJC
         || (hv1[11]||bbl3[6]>0) // LAr
         || bbl3[8]>0 // Spacal
         ))
     return(kFALSE);
   else 
     return(kTRUE);   
} 




Bool_t TTools::IsCosmic(void)
{
// Check if QBGFMAR fired or not
 
  Bool_t selection = kFALSE;

  static H1IntPtr Ibg ("Ibg");     //standard array of background finders
//  static H1IntPtr Ibgam ("Ibgam"); // Add array of background finders

  // Check bits 0-8 of Ibgf
  for(Int_t iBit=0;iBit<9;iBit++)
    if (((*Ibg >> iBit) & 1) == 1){
          selection = kTRUE;
      break;
    }

  /*
  // Check bits 0-8 of Ibgfam
  for(Int_t iBit=0;iBit<7;iBit++)
    if ((*Ibgam & (1 << iBit)) == 1){
      selection = kTRUE;
      break;
    }
  */

  return(selection);

}

Bool_t TTools::IsIbgBit(Int_t iBit)
{
// Check if QBGFMAR fired or not
// bits ranging from 0 to 9 !
 
  Bool_t selection = kFALSE;
  static H1IntPtr Ibg ("Ibg");     //standard array of background finders

    if (((*Ibg >> iBit) & 1) == 1) selection = kTRUE;
  return(selection);

}


Bool_t TTools::IsL4Trigg(Int_t iTrigg)
{
// Has trigger iTrigg fired (L4 verified)
  H1BytePtr STrigg("Il4vst");

  // Has trigger iTrigg really fired
  if (STrigg[iTrigg] == kTRUE)
    return(kTRUE);
  else
    return(kFALSE);
}

Bool_t TTools::IsL1Trigg(Int_t iTrigg)
{
// Has trigger iTrigg fired (L1 actual subtrigger)
  H1BytePtr L1Trigg("Il1ac");

  // Has trigger iTrigg really fired
  if (L1Trigg[iTrigg] == kTRUE)
    return(kTRUE);
  else
    return(kFALSE);
}

Bool_t TTools::IsL2TTTrigg(Int_t iTrigg)
{
// Has trigger iTrigg fired (L1 actual subtrigger)
  H1BytePtr L2TTTrigg("Il2tte");

  // Has trigger iTrigg really fired
  if (L2TTTrigg[iTrigg] == kTRUE)
    return(kTRUE);
  else
    return(kFALSE);
}



Bool_t TTools::Is2004CalibratedTimingOK(void)
{
// Check timing of event, adapted to calibrated LAr timing in 2.7.7
  
  if(MCFlag) return kTRUE;      //--- valid only if DATA
  Bool_t selection = kFALSE;

  // Get timing from HAT
  static H1FloatPtr t0_CJC("T0_CJC");
  static H1FloatPtr t0_LAr("T0_LAr");

  // define offsets for timing determination of LAr
  const Float_t OffsetLAr0304 = -0.0097;

  // define offsets for timing determination of CJC
   Float_t OffsetCJC0304 = 458.0;
  
  // define cuts
  Float_t CJCCutWidth=50.0; // 50 ticks == 0.1 BC
//   Float_t LArCutWidth=0.54; // in BC
  Float_t LArCutWidth=0.7; // in BC

  if(RunNumber >=300000 && RunNumber <= 368000  ) {
        OffsetCJC0304 = 435.0;
        CJCCutWidth=75.0;
  }
  
  // Final timings values
  Float_t FinalLArTim=0.0; 
  Float_t FinalCJCTim=0.0;

  // Get present LAr /CJC timing
  FinalLArTim=*t0_LAr; 
  FinalCJCTim=*t0_CJC;
  FinalLArTim = FinalLArTim - OffsetLAr0304;
  FinalCJCTim = FinalCJCTim - OffsetCJC0304;


  Int_t NbBCLAr=-99;
  Int_t NbBCCJC=+99;

  if(fabs(FinalCJCTim)<CJCCutWidth) NbBCCJC=0;
  if(fabs(FinalCJCTim-500.)<CJCCutWidth) NbBCCJC=1;
  if(fabs(FinalCJCTim-1000.)<CJCCutWidth) NbBCCJC=2;
  if(fabs(FinalCJCTim+500.)<CJCCutWidth) NbBCCJC=-1;
  if(fabs(FinalCJCTim+1000.)<CJCCutWidth) NbBCCJC=-2;
  
  if(fabs(FinalLArTim)<LArCutWidth) NbBCLAr=0;
  if(fabs(FinalLArTim-1.)<LArCutWidth&&NbBCLAr!=0) NbBCLAr=1;
  if(fabs(FinalLArTim-2.)<LArCutWidth&&NbBCLAr!=1) NbBCLAr=2;
  if(fabs(FinalLArTim+1.)<LArCutWidth&&NbBCLAr!=0) NbBCLAr=-1;
  if(fabs(FinalLArTim+2.)<LArCutWidth&&NbBCLAr!=1) NbBCLAr=-2;
//   if(fabs(FinalLArTim-1.)<LArCutWidth&&NbBCCJC==1) NbBCLAr=1;
//   if(fabs(FinalLArTim-2.)<LArCutWidth&&NbBCCJC==2) NbBCLAr=2;
//   if(fabs(FinalLArTim+1.)<LArCutWidth&&NbBCCJC==-1) NbBCLAr=-1;
//   if(fabs(FinalLArTim+2.)<LArCutWidth&&NbBCCJC==-2) NbBCLAr=-2;


  Bool_t NbTrkOk=kFALSE;
  static H1ShortPtr ncent("NumCentralTracks");
  if((*ncent)>0.){
       NbTrkOk=kTRUE;
//        NbBCCJC=0; //---- if no tracks, CJC timing not used for cut !
  }// end of theta max cut

  if((NbBCLAr==NbBCCJC)||(NbTrkOk==kFALSE&&NbBCLAr==0)) selection=kTRUE;

  
  return(selection);   
}

 
 
Bool_t TTools::IsTimingOk(void)
{
// Check timing of event
  if(MCFlag) return kTRUE;      //--- valid only if DATA
  
  if(IsDataHeraII()) return Is2004CalibratedTimingOK();
  
  Bool_t selection = kTRUE;

  // define offsets for timing determination of CJC (1.guess)
  const Float_t Offset100000 = 500;
  const Float_t Offset209999 = 400;
  const Float_t Offset221499 = 437;
  const Float_t Offset244999 = 462;
  const Float_t Offset245000 = 444;
  const Float_t OffsetCJC03 = 450;

  // define offsets for timing determination of LAr
  const Float_t OffsetLAr94 = -0.0;
  const Float_t OffsetLAr95 = -0.055;
  const Float_t OffsetLAr96 = -0.022;
  const Float_t OffsetLAr97 = -0.012;
  const Float_t OffsetLAr9899e = 0.045;
  const Float_t OffsetLAr9900 = 0.077;
  const Float_t OffsetLAr03 = 0.125;

  // define offsets for timing determination of CJC
  const Float_t OffsetCJC94 = 0.0;
  const Float_t OffsetCJC95 = -1.48;
  const Float_t OffsetCJC96 = 10.61;
  const Float_t OffsetCJC97 = 1.61;
  const Float_t OffsetCJC9899e = 15.23;
  const Float_t OffsetCJC9900 = 13.41;

  // Final timings values
  Float_t FinalLArTim=0.0; 
  Float_t FinalCJCTim=0.0;


  // Get timing from HAT
  static H1FloatPtr t0_CJC("T0_CJC");
  static H1FloatPtr t0_LAr("T0_LAr");

  // Timing de CJC
  Float_t TCJC = 1000000;
  if (RunNumber < 100000){
    //    TCJC = *t0_CJC - Offset100000;
    TCJC = Offset100000;
  }
  else if (RunNumber >= 100000 && RunNumber <= 209999)
    TCJC = *t0_CJC - Offset209999;
  else if (RunNumber > 209999 && RunNumber <= 221499)
    TCJC = *t0_CJC - Offset221499;
  else if (RunNumber > 221499 && RunNumber <= 244999)
    TCJC = *t0_CJC - Offset244999;
  else if (RunNumber >= 245000 && RunNumber < 300000)
    TCJC = *t0_CJC - Offset245000;
  else if (IsData2003()){
    TCJC = *t0_CJC - OffsetCJC03;
  } else{
    printf("No timing available\n");
    return(kTRUE);
  }  

  // Get present LAr /CJC timing
  FinalLArTim=*t0_LAr; 
  FinalCJCTim=TCJC;

  // Correction finale 
  if (RunNumber >= 76345 && RunNumber <= 90417){
    //    FinalLArTim = FinalLArTim - OffsetLAr94;
    //    FinalCJCTim = FinalCJCTim - OffsetCJC94;
    FinalLArTim = OffsetLAr94;
    FinalCJCTim = OffsetCJC94;
  }
  else if (RunNumber >= 101902 && RunNumber <= 131045){
    FinalLArTim = FinalLArTim - OffsetLAr95;
    FinalCJCTim = FinalCJCTim - OffsetCJC95;
  }
  else if (RunNumber >= 144000 && RunNumber <= 171573){
    FinalLArTim = FinalLArTim - OffsetLAr96;
    FinalCJCTim = FinalCJCTim - OffsetCJC96;
  }
  else if (RunNumber >= 177920 && RunNumber <= 200407){
    FinalLArTim = FinalLArTim - OffsetLAr97;
    FinalCJCTim = FinalCJCTim - OffsetCJC97;
  }
  else if (RunNumber >= 218217 && RunNumber <= 241649){
    FinalLArTim = FinalLArTim - OffsetLAr9899e;
    FinalCJCTim = FinalCJCTim - OffsetCJC9899e;
  }
  else if (RunNumber >= 244968 && RunNumber <= 279215){
    FinalLArTim = FinalLArTim - OffsetLAr9900;
    FinalCJCTim = FinalCJCTim - OffsetCJC9900;
  }
  else if (IsData2003()){
    FinalLArTim = FinalLArTim - OffsetLAr03;
    FinalCJCTim = TCJC;
  }
  else {
    FinalLArTim = 0.0;
    FinalCJCTim = 0.0;
  }

//   // Fill timing histos
//   hEventtLAr->Fill(FinalLArTim);
//   hEventtCJC->Fill(FinalCJCTim);
//   hEventTTRTCAL->Fill(FinalLArTim,FinalCJCTim);


  // Check LAr timing

//  printf("LAr Timing = %f\n",TMath::Abs(FinalLArTim));
  if (TMath::Abs(FinalLArTim) > 0.7){
     return(kFALSE);
  }else if(IsData2003()&&(TMath::Abs(FinalLArTim) > 0.8)){
    return(kFALSE);
  }

  // Check absolute T_CJC 500 ticks = 1 BC
  
//  printf("CJC Timing = %f\n",TMath::Abs(FinalCJCTim));
  if (TMath::Abs(FinalCJCTim) > 60){
    return(kFALSE);
  }else if(IsData2003()&&(FinalCJCTim > 60 || FinalCJCTim < -100 )){
    return(kFALSE);
  }

  return(selection);   
}
 
 

// Bool_t TTools::IsTimingOk()
// {
// // Check timing of event
//   if(MCFlag) return kTRUE;   //--- valid only if DATA
//   
//   Bool_t selection = kTRUE;
// 
//   // define offsets for timing determination
//   const Float_t Offset100000 = 500;
//   const Float_t Offset209999 = 400;
//   const Float_t Offset221499 = 437;
//   const Float_t Offset244999 = 462;
//   const Float_t Offset245000 = 444;
// 
//   // Get timing from HAT
//   H1FloatPtr t0_CJC("T0_CJC");
//   H1FloatPtr t0_LAr("T0_LAr");
// 
//   // Timing de CJC
//   Float_t TCJC = 1000000;
//   if (RunNumber < 100000)
//     TCJC = *t0_CJC - Offset100000;
//   else if (RunNumber >= 100000 && RunNumber <= 209999)
//     TCJC = *t0_CJC - Offset209999;
//   else if (RunNumber > 209999 && RunNumber <= 221499)
//     TCJC = *t0_CJC - Offset221499;
//   else if (RunNumber > 221499 && RunNumber <= 244999)
//     TCJC = *t0_CJC - Offset244999;
//   else if (RunNumber >= 245000)
//     TCJC = *t0_CJC - Offset245000;
//   else{
// //    printf("No timing available\n");
//     return(kTRUE);
//   }  
// 
//   // Fill timing histos
//   const Float_t FinalDATAOffsetLAr = 0.077;
//   const Float_t FinalDATAOffsetCJC = 13.41;
//   Float_t FinalLArTim=*t0_LAr-FinalDATAOffsetLAr; 
//   Float_t FinalCJCTim=TCJC-FinalDATAOffsetCJC; 
// 
// 
// 
//   // Check LAr timing
//   if (TMath::Abs(FinalLArTim) > 0.7){
// //    printf("LAr Timing = %f\n",FinalLArTim);
//     return(kFALSE);
//   }
// 
//   // Check absolute T_CJC 500 ticks = 1 BC
//   if (TMath::Abs(FinalCJCTim) > 50){
// //    printf("CJC Timing = %f\n",FinalCJCTim);
//     return(kFALSE);
//   }
// 
//   return(selection);   
// }

 Bool_t TTools::IsCCCosmic(void)
{
// Is Special CC cosmic
// Test 24 bits of QBGFMAR

  Bool_t selection = kFALSE;

  static H1IntPtr Ibg ("Ibg");     //standard array of background finders
  static H1IntPtr Ibgam ("Ibgam"); // Add array of background finders
  static H1IntPtr Ibgfm ("Ibgfm"); // Add array of background finders

  // Check bits 0-9 of Ibg
  for(Int_t iBit=0;iBit<10;iBit++)
    if (((*Ibg >> iBit) & 1) == 1){
          selection = kTRUE;
      break;
    }

  // Check bits 0-5 and 8,9 of Ibgfm
  for(Int_t iBit=0;iBit<25;iBit++){
    if (((*Ibgfm >> iBit) & 1) == 1){
         if(iBit !=6 ) {//--- too inefficient bits disabled
                selection = kTRUE;
//              printf("CC Bg finder: Ibgfm bit %d fired !\n",iBit);
                break;
         }
    }
  }

  // Check bits 0-8 of Ibgfam
  for(Int_t iBit=0;iBit<9;iBit++)
    if (((*Ibgam >> iBit) & 1) == 1){
         if(iBit!=3 && iBit !=4 && iBit !=8 ) {
                selection = kTRUE; //--- too inefficient bits disabled
//              printf("UK CC Bg finder: Ibgfam bit %d fired !\n",iBit);
                break;
        }
    }

  // return selection
  return(selection);
}




Bool_t TTools::IsNonEP03(void)
{
// Is NOISE (electronic LAr noisy runs for 2003 oct-dec)

  Bool_t selection = kFALSE;
    
  if (RunNumber == 360178) return(kTRUE);
  else if(RunNumber == 360180) return(kTRUE);
  else if(RunNumber == 360181) return(kTRUE);
  else if(RunNumber == 360191) return(kTRUE);
  else if(RunNumber == 360213) return(kTRUE);
  else if(RunNumber == 360223) return(kTRUE);
  else if(RunNumber == 360224) return(kTRUE);
  else if(RunNumber == 360227) return(kTRUE);
  else if(RunNumber == 360229) return(kTRUE);
  else if(RunNumber == 360236) return(kTRUE);
//--- Runs with calo "OFF" ?? reported by Benji  
  else if(RunNumber >= 358730 && RunNumber <= 358735) return(kTRUE);
  else if(RunNumber == 367227) return(kTRUE);
  else if(RunNumber == 367523) return(kTRUE);
//--- jnp scanning generic analysis... (13/12/03)
  else if(RunNumber == 359925 && EventNumber ==  15639) return(kTRUE);
  else if(RunNumber == 360296 && EventNumber ==  14958) return(kTRUE);
  else if(RunNumber == 360296&& EventNumber == 137296 ) return(kTRUE);
  else if(RunNumber == 360298 && EventNumber == 58016 ) return(kTRUE);
  else if(RunNumber == 360409 && EventNumber == 648 ) return(kTRUE);
  else if(RunNumber == 361096 && EventNumber == 12481 ) return(kTRUE);
  else if(RunNumber ==  361096&& EventNumber == 37679 ) return(kTRUE);
  else if(RunNumber == 361175 && EventNumber == 21831 ) return(kTRUE);
  else if(RunNumber == 362266 && EventNumber == 4434 ) return(kTRUE);
  else if(RunNumber == 362268 && EventNumber == 12951 ) return(kTRUE);
  else if(RunNumber == 362412 && EventNumber ==  8810) return(kTRUE);
  else if(RunNumber == 362526 && EventNumber == 13957 ) return(kTRUE);
  else if(RunNumber == 363060 && EventNumber == 29468 ) return(kTRUE);
  else if(RunNumber == 363076 && EventNumber == 4538 ) return(kTRUE);
  else if(RunNumber == 363290 && EventNumber == 82572 ) return(kTRUE);
  else if(RunNumber == 363313 && EventNumber ==  29343) return(kTRUE);
  else if(RunNumber == 363411 && EventNumber == 22679 ) return(kTRUE);
  else if(RunNumber == 363411 && EventNumber ==  38055) return(kTRUE);
  else if(RunNumber ==  363455&& EventNumber ==  30769) return(kTRUE);
  else if(RunNumber == 363475 && EventNumber == 26888 ) return(kTRUE);
  else if(RunNumber == 363475 && EventNumber == 53887 ) return(kTRUE);
  else if(RunNumber == 364484 && EventNumber == 31657 ) return(kTRUE);
  else if(RunNumber == 365153 && EventNumber == 2943 ) return(kTRUE);
  else if(RunNumber == 365153 && EventNumber ==  32023) return(kTRUE);
  else if(RunNumber == 365414 && EventNumber == 31028) return(kTRUE);
  else if(RunNumber == 365512 && EventNumber == 29897 ) return(kTRUE);
  else if(RunNumber == 365520 && EventNumber == 49896 ) return(kTRUE);
//---- jjjnp scanning generic analysis ...(13/12/03)
  else if(RunNumber == 365164 && EventNumber == 84268 ) return(kTRUE);
//--- munp scanning generic analysis...(13/12/03)
  else if(RunNumber ==  358760&& EventNumber ==  37303) return(kTRUE);
  else if(RunNumber ==  360135&& EventNumber ==  6818) return(kTRUE);
  else if(RunNumber ==  360296&& EventNumber ==  8968) return(kTRUE);
  else if(RunNumber ==  360296&& EventNumber ==  119968) return(kTRUE);
  else if(RunNumber ==  361309&& EventNumber ==  1202) return(kTRUE);
  else if(RunNumber ==  362453&& EventNumber ==  98583) return(kTRUE);
  else if(RunNumber ==  363597&& EventNumber ==  19723) return(kTRUE);
  else if(RunNumber ==  363666&& EventNumber ==  6819) return(kTRUE);
  else if(RunNumber ==  363865&& EventNumber == 44606 ) return(kTRUE);
  else if(RunNumber ==  363874&& EventNumber ==  3430) return(kTRUE);
  else if(RunNumber ==  365083&& EventNumber ==  44578) return(kTRUE);
//--- CC Dq scanning ...(16/12/03)
  else if(RunNumber ==  360121&& EventNumber == 5463) return(kTRUE);
  else if(RunNumber ==  363284&& EventNumber == 102592) return(kTRUE);
  else if(RunNumber ==  363459&& EventNumber == 14932) return(kTRUE);
  else if(RunNumber ==  364484&& EventNumber == 10143) return(kTRUE);


  return(selection);
}

 

Bool_t TTools::IsNonEP()
{
// Is NOISE or cosmic
  if(MCFlag) return kFALSE;     //--- valid only if DATA

  Bool_t selection = kFALSE;

  // NOISE + EYE ID OF COSMICS RUNS/EVENTS
  // Check if actual event is noisy or cosmic

  // Eliminate cosmic events
  if (RunNumber == 265513 && EventNumber == 80488)
    return(kTRUE);
  else if (RunNumber == 274259 && EventNumber == 103462)
    return(kTRUE);
  else if (RunNumber == 262491 && EventNumber == 103171)
    return(kTRUE);
  else if (RunNumber == 252787 && EventNumber == 84899)
    return(kTRUE);
  else if (RunNumber == 256965 && EventNumber == 160670)
    return(kTRUE);
  else if (RunNumber == 258147 && EventNumber == 42406)
    return(kTRUE);
  else if (RunNumber == 258678 && EventNumber == 179589)
    return(kTRUE);
  else if (RunNumber == 257627 && EventNumber == 48985)
    return(kTRUE);
  else if (RunNumber == 258808 && EventNumber == 28416)
    return(kTRUE);
  else if (RunNumber == 259734 && EventNumber == 60542)
    return(kTRUE);
  else if (RunNumber == 261338 && EventNumber == 22402)
    return(kTRUE);
  // j-pho scanning
  else if (RunNumber == 253523 && EventNumber == 1603)
    return(kTRUE);
  else if (RunNumber == 256697 && EventNumber == 9599)
    return(kTRUE);
  else if (RunNumber == 256965 && EventNumber == 169670)
    return(kTRUE);
  else if (RunNumber == 258113 && EventNumber == 42406)
    return(kTRUE);
  else if (RunNumber == 258147 && EventNumber == 179589)
    return(kTRUE);
  else if (RunNumber == 258147 && EventNumber == 333546)
    return(kTRUE);
  else if (RunNumber == 258674 && EventNumber == 64200)
    return(kTRUE);
  else if (RunNumber == 258678 && EventNumber == 48985)
    return(kTRUE);
  else if (RunNumber == 257691 && EventNumber == 2219)
    return(kTRUE);
  else if (RunNumber == 258025 && EventNumber == 167887)
    return(kTRUE);
  else if (RunNumber == 258808 && EventNumber == 28601)
    return(kTRUE);
  else if (RunNumber == 259107 && EventNumber == 95002)
    return(kTRUE);
  else if (RunNumber == 260845 && EventNumber == 40143)
    return(kTRUE);
  else if (RunNumber == 260036 && EventNumber == 84418)
    return(kTRUE);
  else if (RunNumber == 260210 && EventNumber == 8222)
    return(kTRUE);
  else if (RunNumber == 260354 && EventNumber == 126164)
    return(kTRUE);
  // Eye scanning e-mu
  else if (RunNumber == 264246 && EventNumber == 12180)
    return(kTRUE);
  else if (RunNumber == 268319 && EventNumber == 28736)
    return(kTRUE);
  else if (RunNumber == 275086 && EventNumber == 35719)
    return(kTRUE);
  // Eye scanning mu-mu
  else if (RunNumber == 252377 && EventNumber == 93828)
    return(kTRUE);
  else if (RunNumber == 250996 && EventNumber == 48582)
    return(kTRUE);
  else if (RunNumber == 254433 && EventNumber == 21344)
    return(kTRUE);
  else if (RunNumber == 255109 && EventNumber == 75418)
    return(kTRUE);
  else if (RunNumber == 259107 && EventNumber == 66564)
    return(kTRUE);
  else if (RunNumber == 249271 && EventNumber == 114116)
    return(kTRUE);
  else if (RunNumber == 278496 && EventNumber == 21284)
    return(kTRUE);
  else if (RunNumber == 278541 && EventNumber == 17094)
    return(kTRUE);
  else if (RunNumber == 264171 && EventNumber == 155107)
    return(kTRUE);
  else if (RunNumber == 264257 && EventNumber == 9015)
    return(kTRUE);
  else if (RunNumber == 264442 && EventNumber == 57876)
    return(kTRUE);
  else if (RunNumber == 265270 && EventNumber == 13543)
    return(kTRUE);
  else if (RunNumber == 267206 && EventNumber == 205856)
    return(kTRUE);
  else if (RunNumber == 267207 && EventNumber == 112373)
    return(kTRUE);
  else if (RunNumber == 268338 && EventNumber == 57765)
    return(kTRUE);
  else if (RunNumber == 262534 && EventNumber == 5367)
    return(kTRUE);
  else if (RunNumber == 262752 && EventNumber == 8585)
    return(kTRUE);
  else if (RunNumber == 269754 && EventNumber == 1409)
    return(kTRUE);
  else if (RunNumber == 269633 && EventNumber == 17192)
    return(kTRUE);
  else if (RunNumber == 272932 && EventNumber == 58808)
    return(kTRUE);
  else if (RunNumber == 272037 && EventNumber == 77213)
    return(kTRUE);
  else if (RunNumber == 275081 && EventNumber == 1076)
    return(kTRUE);
  else if (RunNumber == 275591 && EventNumber == 143305)
    return(kTRUE);
  else if (RunNumber == 275615 && EventNumber == 12772)
    return(kTRUE);
  else if (RunNumber == 275722 && EventNumber == 69943)
    return(kTRUE);
  else if (RunNumber == 275843 && EventNumber == 31714)
    return(kTRUE);
  else if (RunNumber == 275993 && EventNumber == 80659)
    return(kTRUE);
  else if (RunNumber == 276036 && EventNumber == 120319)
    return(kTRUE);
  else if (RunNumber == 276662 && EventNumber == 82873)
    return(kTRUE);
  else if (RunNumber == 276839 && EventNumber == 3475)
    return(kTRUE);
  else if (RunNumber == 276965 && EventNumber == 64052)
    return(kTRUE);
  else if (RunNumber == 277996 && EventNumber == 37571)
    return(kTRUE);
  else if (RunNumber == 270343 && EventNumber == 15397)
    return(kTRUE);
  else if (RunNumber == 270321 && EventNumber == 40892)
    return(kTRUE);
  else if (RunNumber == 263675 && EventNumber == 170286)
    return(kTRUE);
  else if (RunNumber == 278496 && EventNumber == 21234)
    return(kTRUE);
  else if (RunNumber == 265372 && EventNumber == 89946)
    return(kTRUE);
  else if (RunNumber == 275709 && EventNumber == 63682)
    return(kTRUE);
  else if (RunNumber == 278078 && EventNumber == 4137)
    return(kTRUE);
  // Eye scanning mu-j
  else if (RunNumber == 273527 && EventNumber == 4799)
    return(kTRUE);
  else if (RunNumber == 276219 && EventNumber == 148101)
    return(kTRUE);
  else if (RunNumber == 271223 && EventNumber == 987)
    return(kTRUE);
  // Eye scanning mu-np
  else if (RunNumber == 249504 && EventNumber == 73715)
    return(kTRUE);
  else if (RunNumber == 252377 && EventNumber == 105453)
    return(kTRUE);
  else if (RunNumber == 254855 && EventNumber == 23864)
    return(kTRUE);
  else if (RunNumber == 258073 && EventNumber == 201294)
    return(kTRUE);
  else if (RunNumber == 257779 && EventNumber == 88729)
    return(kTRUE);
  else if (RunNumber == 259399 && EventNumber == 16567)
    return(kTRUE);
  else if (RunNumber == 258880 && EventNumber == 184942)
    return(kTRUE);
  else if (RunNumber == 259088 && EventNumber == 62683)
    return(kTRUE);
  else if (RunNumber == 259312 && EventNumber == 179641)
    return(kTRUE);
  else if (RunNumber == 278653 && EventNumber == 92192)
    return(kTRUE);
  else if (RunNumber == 268024 && EventNumber == 25327)
    return(kTRUE);
  else if (RunNumber == 267315 && EventNumber == 11062)
    return(kTRUE);
  else if (RunNumber == 267405 && EventNumber == 13005)
    return(kTRUE);
  else if (RunNumber == 272035 && EventNumber == 194608)
    return(kTRUE);
  else if (RunNumber == 273497 && EventNumber == 53854)
    return(kTRUE);
  else if (RunNumber == 274258 && EventNumber == 63331)
    return(kTRUE);
  else if (RunNumber == 276133 && EventNumber == 18337)
    return(kTRUE);
  else if (RunNumber == 276789 && EventNumber == 33060)
    return(kTRUE);
  else if (RunNumber == 277369 && EventNumber == 131891)
    return(kTRUE);
  else if (RunNumber == 276941 && EventNumber == 85829)
    return(kTRUE);
  else if (RunNumber == 270750 && EventNumber == 149391)
    return(kTRUE);
  else if (RunNumber == 270273 && EventNumber == 54901)
    return(kTRUE);
  // Eye scanning np-pho
  else if (RunNumber == 251667 && EventNumber == 9890)
    return(kTRUE);
  else if (RunNumber == 250414 && EventNumber == 15294)
    return(kTRUE);
  else if (RunNumber == 278460 && EventNumber == 119514)
    return(kTRUE);
  else if (RunNumber == 278849 && EventNumber == 26720)
    return(kTRUE);
  else if (RunNumber == 252642 && EventNumber == 51145)
    return(kTRUE);
  else if (RunNumber == 252725 && EventNumber == 119991)
    return(kTRUE);
  else if (RunNumber == 252792 && EventNumber == 82955)
    return(kTRUE);
  else if (RunNumber == 251866 && EventNumber == 86526)
    return(kTRUE);
  else if (RunNumber == 252380 && EventNumber == 112743)
    return(kTRUE);
  else if (RunNumber == 254229 && EventNumber == 20194)
    return(kTRUE);
  else if (RunNumber == 254385 && EventNumber == 37410)
    return(kTRUE);
  else if (RunNumber == 254818 && EventNumber == 178930)
    return(kTRUE);
  else if (RunNumber == 264403 && EventNumber == 7264)
    return(kTRUE);
  else if (RunNumber == 256562 && EventNumber == 10616)
    return(kTRUE);
  else if (RunNumber == 256425 && EventNumber == 38159)
    return(kTRUE);
  else if (RunNumber == 256496 && EventNumber == 10830)
    return(kTRUE);
  else if (RunNumber == 256551 && EventNumber == 31171)
    return(kTRUE);
  else if (RunNumber == 265528 && EventNumber == 15519)
    return(kTRUE);
  else if (RunNumber == 265528 && EventNumber == 92368)
    return(kTRUE);
  else if (RunNumber == 265019 && EventNumber == 18073)
    return(kTRUE);
  else if (RunNumber == 265180 && EventNumber == 71515)
    return(kTRUE);
  else if (RunNumber == 266281 && EventNumber == 159387)
    return(kTRUE);
  else if (RunNumber == 266171 && EventNumber == 61707)
    return(kTRUE);
  else if (RunNumber == 259561 && EventNumber == 1312)
    return(kTRUE);
  else if (RunNumber == 259347 && EventNumber == 16880)
    return(kTRUE);
  else if (RunNumber == 267771 && EventNumber == 50401)
    return(kTRUE);
  else if (RunNumber == 267827 && EventNumber == 149334)
    return(kTRUE);
  else if (RunNumber == 266936 && EventNumber == 3510)
    return(kTRUE);
  else if (RunNumber == 267477 && EventNumber == 5168)
    return(kTRUE);
  else if (RunNumber == 259733 && EventNumber == 15503)
    return(kTRUE);
  else if (RunNumber == 260155 && EventNumber == 9196)
    return(kTRUE);
  else if (RunNumber == 268338 && EventNumber == 69658)
    return(kTRUE);
  else if (RunNumber == 268652 && EventNumber == 82354)
    return(kTRUE);
  else if (RunNumber == 268661 && EventNumber == 3143)
    return(kTRUE);
  else if (RunNumber == 268072 && EventNumber == 88380)
    return(kTRUE);
  else if (RunNumber == 246631 && EventNumber == 20449)
    return(kTRUE);
  else if (RunNumber == 262817 && EventNumber == 24252)
    return(kTRUE);
  else if (RunNumber == 270024 && EventNumber == 26844)
    return(kTRUE);
  else if (RunNumber == 270819 && EventNumber == 51303)
    return(kTRUE);
  else if (RunNumber == 272793 && EventNumber == 132804)
    return(kTRUE);
  else if (RunNumber == 271156 && EventNumber == 6015)
    return(kTRUE);
  else if (RunNumber == 272319 && EventNumber == 133943)
    return(kTRUE);
  else if (RunNumber == 273242 && EventNumber == 56531)
    return(kTRUE);
  else if (RunNumber == 274918 && EventNumber == 76937)
    return(kTRUE);
  else if (RunNumber == 274480 && EventNumber == 38737)
    return(kTRUE);
  else if (RunNumber == 274627 && EventNumber == 105476)
    return(kTRUE);
  else if (RunNumber == 275242 && EventNumber == 16100)
    return(kTRUE);
  else if (RunNumber == 275353 && EventNumber == 111236)
    return(kTRUE);
  else if (RunNumber == 275381 && EventNumber == 50776)
    return(kTRUE);
  else if (RunNumber == 276133 && EventNumber == 37695)
    return(kTRUE);
  else if (RunNumber == 276839 && EventNumber == 97017)
    return(kTRUE);
  else if (RunNumber == 276459 && EventNumber == 130863)
    return(kTRUE);
  else if (RunNumber == 277452 && EventNumber == 30340)
    return(kTRUE);
  else if (RunNumber == 270464 && EventNumber == 169398)
    return(kTRUE);
  else if (RunNumber == 270509 && EventNumber == 14807)
    return(kTRUE);
  else if (RunNumber == 270099 && EventNumber == 44898)
    return(kTRUE);
  else if (RunNumber == 270137 && EventNumber == 56412)
    return(kTRUE);
  else if (RunNumber == 270147 && EventNumber == 195335)
    return(kTRUE);
  else if (RunNumber == 263754 && EventNumber == 14959)
    return(kTRUE);
  else if (RunNumber == 263846 && EventNumber == 129382)
    return(kTRUE);
  else if (RunNumber == 275616 && EventNumber == 47523)
    return(kTRUE);
  else if (RunNumber == 249108 && EventNumber == 24810)
    return(kTRUE);
  else if (RunNumber == 275715 && EventNumber == 3919)
    return(kTRUE);
  else if(RunNumber == 278414 && EventNumber == 21088)
    return(kTRUE);
  else if(RunNumber == 278999 && EventNumber == 6385)
    return(kTRUE);
  else if(RunNumber == 265490 && EventNumber == 146599)
    return(kTRUE);
  else if(RunNumber == 266171 && EventNumber == 61701)
    return(kTRUE);
  else if(RunNumber == 267949 && EventNumber == 54921)
    return(kTRUE);
  else if(RunNumber == 268211 && EventNumber == 8027)
    return(kTRUE);
  else if(RunNumber == 275071 && EventNumber == 7335)
    return(kTRUE);
  else if(RunNumber == 277996 && EventNumber == 43191)
    return(kTRUE);
  // Eliminate hot cells noise events 
  if (RunNumber == 274737 && EventNumber == 99795)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57052)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57162)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57236)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57447)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57488)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57497)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57543)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57748)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57794)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 58242)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 58257)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 58457)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 58697)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 58700)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 58783)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 58827)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 58945)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 58974)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 59113)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 59761)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 59852)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 59984)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 60304)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 60116)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 60450)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 60598)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 60618)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 60943)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 60988)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61043)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61100)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61186)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61206)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61321)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61600)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61645)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61702)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61780)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 61854)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 62039)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 62589)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 62690)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 62772)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 62859)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 62874)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 63060)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 63433)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 63610)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 63731)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 63797)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 63814)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 64018)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 64028)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 64311)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 64505)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 64628)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 64741)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 65027)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 65235)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 65366)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 65909)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 66073)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 66589)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 66606)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 66674)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 94597)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 94806)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 94925)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 94962)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 95072)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 95088)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 95225)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 95738)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 95924)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 96060)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 96130)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 96389)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 96760)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 96958)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 96991)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 97031)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 97103)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 97396)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 97443)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 97556)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 97827)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 98001)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 98233)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 98948)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 98964)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 99038)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 99153)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 99256)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 99369)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 99465)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 99597)
    return(kTRUE);
  else if (RunNumber == 269604 && EventNumber == 175)
    return(kTRUE);
  else if (RunNumber == 270856 && EventNumber == 75790)
    return(kTRUE);
  else if (RunNumber == 272435 && EventNumber == 25542)
    return(kTRUE);
  else if (RunNumber == 273711 && EventNumber == 45868)
    return(kTRUE);
  else if (RunNumber == 273711 && EventNumber == 45952)
    return(kTRUE);
  else if (RunNumber == 273711 && EventNumber == 45854)
    return(kTRUE);
  else if (RunNumber == 273744 && EventNumber == 4529)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 99765)
    return(kTRUE);
  else if (RunNumber == 275083 && EventNumber == 53806)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57167)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 57263)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 59994)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 60306)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 60416)
    return(kTRUE);
  else if (RunNumber == 274673 && EventNumber == 2066)
    return(kTRUE);
  else if (RunNumber == 274673 && EventNumber == 2123)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 55104)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56001)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56002)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56088)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56095)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56263)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56272)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56311)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56746)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56760)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56770)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56945)
    return(kTRUE);
  else if (RunNumber == 274737 && EventNumber == 56991)
    return(kTRUE);
  else if (RunNumber == 270132 && EventNumber == 118742)
    return(kTRUE);
  else if (RunNumber == 270132 && EventNumber == 118857)
    return(kTRUE);
  else if (RunNumber == 270132 && EventNumber == 118883)
    return(kTRUE);
  else if (RunNumber == 270132 && EventNumber == 119033)
    return(kTRUE);
  else if (RunNumber == 270270 && EventNumber == 103708)
    return(kTRUE);
  else if (RunNumber == 270270 && EventNumber == 111491)
    return(kTRUE);
  else if (RunNumber == 270270 && EventNumber == 111659)
    return(kTRUE);
  else if (RunNumber == 270270 && EventNumber == 111786)
    return(kTRUE);
  // mu-np scanning
  else if (RunNumber == 279003 && EventNumber == 20339)
    return(kTRUE);
  else if (RunNumber == 278612 && EventNumber == 71506)
    return(kTRUE);
  else if (RunNumber == 264735 && EventNumber == 83011)
    return(kTRUE);
  else if (RunNumber == 264403 && EventNumber == 117528)
    return(kTRUE);
  else if (RunNumber == 266830 && EventNumber == 61617)
    return(kTRUE);
  else if (RunNumber == 266000 && EventNumber == 216997)
    return(kTRUE);
  else if (RunNumber == 266926 && EventNumber == 150936)
    return(kTRUE);
  else if (RunNumber == 267665 && EventNumber == 82141)
    return(kTRUE);
  else if (RunNumber == 269049 && EventNumber == 14371)
    return(kTRUE);
  else if (RunNumber == 271908 && EventNumber == 62529)
    return(kTRUE);
  else if (RunNumber == 270899 && EventNumber == 2338)
    return(kTRUE);
  else if (RunNumber == 271115 && EventNumber == 139166)
    return(kTRUE);
  else if (RunNumber == 272398 && EventNumber == 16831)
    return(kTRUE);
  else if (RunNumber == 272218 && EventNumber == 12581)
    return(kTRUE);
  else if (RunNumber == 273348 && EventNumber == 25257)
    return(kTRUE);
  else if (RunNumber == 273497 && EventNumber == 116838)
    return(kTRUE);
  else if (RunNumber == 273105 && EventNumber == 231053)
    return(kTRUE);
  else if (RunNumber == 273230 && EventNumber == 26104)
    return(kTRUE);
  else if (RunNumber == 274361 && EventNumber == 95877)
    return(kTRUE);
  else if (RunNumber == 273694 && EventNumber == 27612)
    return(kTRUE);
  else if (RunNumber == 273738 && EventNumber == 37091)
    return(kTRUE);
  else if (RunNumber == 273744 && EventNumber == 2549)
    return(kTRUE);
  else if (RunNumber == 274736 && EventNumber == 17683)
    return(kTRUE);
  else if (RunNumber == 275256 && EventNumber == 54317)
    return(kTRUE);
  else if (RunNumber == 275401 && EventNumber == 31574)
    return(kTRUE);
  else if (RunNumber == 276305 && EventNumber == 142797)
    return(kTRUE);
  else if (RunNumber == 276132 && EventNumber == 98651)
    return(kTRUE);
  else if (RunNumber == 276756 && EventNumber == 67880)
    return(kTRUE);
  else if (RunNumber == 276462 && EventNumber == 18266)
    return(kTRUE);
  else if (RunNumber == 276974 && EventNumber == 6466)
    return(kTRUE);
  else if (RunNumber == 278206 && EventNumber == 27372)
    return(kTRUE);
  else if (RunNumber == 278307 && EventNumber == 13763)
    return(kTRUE);
  else if (RunNumber == 268008 && EventNumber == 214400)
    return(kTRUE);
  else if (RunNumber == 249371 && EventNumber == 71211)
    return(kTRUE);
  else if (RunNumber == 252276 && EventNumber == 118618)
    return(kTRUE);
  else if (RunNumber == 252428 && EventNumber == 31173)
    return(kTRUE);
  else if (RunNumber == 254195 && EventNumber == 1244)
    return(kTRUE);
  else if (RunNumber == 254466 && EventNumber == 117267)
    return(kTRUE);
  else if (RunNumber == 253567 && EventNumber == 1006)
    return(kTRUE);
  else if (RunNumber == 253973 && EventNumber == 10515)
    return(kTRUE);
  else if (RunNumber == 256916 && EventNumber == 193529)
    return(kTRUE);
  else if (RunNumber == 254853 && EventNumber == 26151)
    return(kTRUE);
  else if (RunNumber == 255077 && EventNumber == 27797)
    return(kTRUE);
  else if (RunNumber == 256343 && EventNumber == 7200)
    return(kTRUE);
  else if (RunNumber == 256553 && EventNumber == 29340)
    return(kTRUE);
  else if (RunNumber == 258708 && EventNumber == 125833)
    return(kTRUE);
  else if (RunNumber == 257772 && EventNumber == 50965)
    return(kTRUE);
  else if (RunNumber == 259688 && EventNumber == 19688)
    return(kTRUE);
  else if (RunNumber == 258812 && EventNumber == 67459)
    return(kTRUE);
  else if (RunNumber == 258867 && EventNumber == 175276)
    return(kTRUE);
  else if (RunNumber == 258917 && EventNumber == 119958)
    return(kTRUE);
  else if (RunNumber == 259090 && EventNumber == 69640)
    return(kTRUE);
  else if (RunNumber == 260035 && EventNumber == 94532)
    return(kTRUE);
  else if (RunNumber == 248744 && EventNumber == 20110)
    return(kTRUE);
  else if (RunNumber == 249089 && EventNumber == 7075)
    return(kTRUE);

  return(selection);
}

Bool_t TTools::IsCJCDead(Float_t PhiTr )
{

   if(MCFlag>=1)
     return(kFALSE); 

   if ((MCFlag==0 &&RunNumber >= 257590 && RunNumber <= 262144)
        || (MCFlag>0 &&RunNumber >= 257590 && RunNumber <= 262144)) { //---- hole in 99 data
        if(PhiTr*57.3 < -80 && PhiTr*57.3 > -170) return(kTRUE);
        else return(kFALSE);
   }else if((MCFlag==0 && RunNumber<=201519)){  //--- 94-97 data 
        if(PhiTr*57.3 < -110 && PhiTr*57.3 > -130) {
          return(kTRUE);
        }
        else return(kFALSE);
   }
   return(kFALSE); // to make the compiler happy...
}

// Bool_t TMyEvent::IsCJCDead(Float_t PhiTr)
// {
// // Is CJC dead (Dead period in 1999 data)
// // Return kFALSE if CJC alive and kTRUE if CJC dead
// // Check MC or not
//   Float_t lumitot=65.85;
//   if(MCFlag>=3&&MCFlag <12) lumitot= 79.5405;
// 
//    if ((MCFlag==0 &&RunNumber >= 257590 && RunNumber <= 262144)
//         || (MCFlag>=12 && MCFlag<40 &&RunNumber >= 257637 && RunNumber <= 262144)) { //---- hole in 99 data
//    
//      if(PhiTr*57.3 < -80 && PhiTr*57.3 > -170) return(kTRUE);
//      else return(kFALSE);
//    }else if(MCFlag>0 && MCFlag<12) {//---- for old MC's
//      if((gRandom->Rndm()<8.0267/lumitot)&&(PhiTr*57.3 < -80 && PhiTr*57.3 > -170)) return(kTRUE);
//      else return(kFALSE);
//    }else if((MCFlag==0 && RunNumber<=201519)
//            ||(RunYear=TMarEvent::Y9497)){    //--- 94-97 data 
//      if(PhiTr*57.3 < -110 && PhiTr*57.3 > -130) {
//        //      cout <<"ICI\n";
//        return(kTRUE);
//      }
//      else return(kFALSE);
//         
//    }
//    // Sinon retourner kFALSE
//    return(kFALSE);
// }


Bool_t TTools::IsCJCMCDead(Float_t PhiTr)
{
// Is CJC MC dead (Equivalent to Dead period in 1999 data)
// Return kFALSE if CJC alive and kTRUE if CJC dead

  // Check dead CJC segment 
  // => By eye analysis choice of Phi range
   if (PhiTr*57.3 > -180 && PhiTr*57.3<-150)
     return(kTRUE);
   else 
     return(kFALSE);
}


Bool_t TTools::IsEPzOK(Float_t Epz)
{
// Is E-Pz physical (E-Pz < 75 GeV)

  const Double_t EPZMAX = 75.0;

  // Check event variable Epz
  if (Epz < EPZMAX)
    return(kTRUE);
  else
    return(kFALSE);
}

Float_t TTools::qhadcor_98_00(Int_t ITRA,Int_t ITYPE,Float_t jetpos,Float_t ptin,Float_t qvar){
// ===========================================================================
//
// INPUT :
// -----
//              ITRA = 0  for calibration without track information
//                     1   "     "        with     "         "
//
//              ITYPE = 1  for DATA 98-99 (nrun<233404)
//                      2  for DATA 99    (nrun>233404)
//                      3  for NC MC 98 e-p
//                      4  for DATA 99-00
//                      5  for NC MC 99 e+p
//
//              jetpos = theta of the jet
//
//              ptin = transverse inmpulsion of the jet
//
//              qvar = energie or pt or px or py or pz
//
//
// OUTPUT :
// ------       qhadcor_98_00 = corrected varin
//
//                                                                             
  // No track info
  Float_t a_da9900[8]={0.868093,0.877024,0.921192,0.949944,0.977572,1.09807,1.17005,0.95448};
  Float_t b_da9900[8]={0.0,-1.08682,-1.70889,0.0,0.0,0.0,-0.175373,0.0};
  Float_t c_da9900[8]={5.65951,10.0505,17.0786,5.55913,5.84501,6.47199,8.74776,5.50889};

  // track info
  Float_t a_da9900_tr[8]={0.868376,0.879944,1.151,0.960589,0.987846,1.15057,1.76864,0.915961};
  Float_t b_da9900_tr[8]={0.0,-1.20845,-1.33753,-1.91097,-1.71629,-1.40015,-0.691473,0.0};
  Float_t c_da9900_tr[8]={5.60348,10.699,196.481,9.64211,9.27777,13.1271,66.3346,2.53625};

  // No track info


  // Define helper vars
  Int_t ith=0; 
  Float_t acorr=0,bcorr=0,ccorr=0,ptinter=0,scale=0;

  // Check validity of input parameters
  if (ITRA != 0 && ITRA != 1)
    return(qvar);
  if (ITYPE != 1 && ITYPE != 2 && ITYPE != 3 && ITYPE != 4 && ITYPE != 5) 
    return(qvar);
  if (ptin<0)
    return(qvar);

  if (jetpos>0 && jetpos<10 ) ith=1;
  if (jetpos>=10 && jetpos<15 ) ith=2;
  if (jetpos>=30 && jetpos<55 ) ith=4;
  if (jetpos>=5 && jetpos<80 ) ith=5;
  if (jetpos>=80 && jetpos<110) ith=6;
  if (jetpos>=110 && jetpos<135) ith=7;
  if (jetpos>=135 && jetpos<180) ith=8;                           

  // DATA 99/00
  if (ITYPE == 4){
    if (!ITRA){
      if (ith){
        acorr = a_da9900[ith];
        bcorr = b_da9900[ith];
        ccorr = c_da9900[ith];
      }
      else return(qvar);
    }
    else if (ITRA == 1){
      if (ith){
        acorr = a_da9900_tr[ith];
        bcorr = b_da9900_tr[ith];
        ccorr = c_da9900_tr[ith];
      }
      else return(qvar);
    }
  }
  // DATA 98/00, MC...

  // Without track
  if (!ITRA){
    // --- Constante for data 98 e-p : bin 8
    //               for data 99 e-p : bin 8
    //               for NC MC  e+p  : bin 8
    if ((ITYPE == 1 && ith==8)||(ITYPE==2 && ith ==8)||(ITYPE==5 && ith==8)){
     scale = acorr;
    }
    // Exponential function for all others
    else{
      scale = acorr*(1-TMath::Exp(bcorr-ptin/ccorr));
      ptinter = ptin/scale;
      scale = acorr*(1-TMath::Exp(bcorr-ptinter/ccorr));
    }
  }

  // With track
  if (ITRA == 1){
    // --- Constante for data 98 e-p  : if,of,cb3,cb2,cb1,bbe
    //               for data 99 e-p  : bbe
    //               for nc mc e+p    : cb3,bbe                                
    if ((ITYPE == 1 && (ith==5 || ith==7 || ith==8 || ith==2 || ith==3 || ith==6)) || (ITYPE==2 && ith==8) || (ITYPE==5 && (ith==5 || ith==8))){
      scale = acorr;
    }
    // Special function for nc mc e-p bbe
    else if ((ITYPE==3 && ith==8)){
      scale = acorr/ptin + ccorr;
      ptinter = ptin /scale;
      scale = acorr/ptinter + ccorr;
    }
    // Exponential function for all others
    else{
      scale = acorr*(1-TMath::Exp(bcorr-ptin/ccorr));
      ptinter = ptin/scale;
      scale = acorr*(1-TMath::Exp(bcorr-ptinter/ccorr));
    }
  }

  // Return correction of qvar
  return(qvar/scale);

}



 
  
Bool_t TTools::IsMCQ2OK(void)
{
// Is Q2 OK (Combination of diff. MCs)

 static H1FloatPtr Q2gki("Q2Gki");

 // Check for RAPGAP NC Q2<100 GeV2
 if (MCFlag != MCRAPGAPNCID)
   return(kTRUE);
 else{
   printf("Found rapgap q2 = %#8.3f\n",*Q2gki);
   if (*Q2gki < 100)
     return(kTRUE);
   else{
     return(kFALSE);
   }
 }
}




Float_t TTools::abscor(Int_t icalfac,Float_t theta_rad,Float_t ptj,Float_t qvar){
// Final calib using qhadcor
// ONLY FORSEEN FOR DATA 99/00
 
  // Vars for calib
  Float_t ej,theta_deg;
  Int_t itype_new,ITRA;
 
  // Get energy
  ej = qvar;
  theta_deg = theta_rad*180/3.1415;

  // 99/00
  ITRA = 1; // (0/1) = (NO TRACK INFO/TRACK INFO)
  if (RunNumber>205000 && RunNumber < 233404 && icalfac == 1)
    itype_new = 1; // data 99 e- before oc
  else if (RunNumber>233404 && RunNumber <241650 && icalfac == 1)
    itype_new = 2; // data 99 e-, after oc 
  else if (RunNumber>241650 && icalfac == 1)
    itype_new = 4;
  else if (icalfac==2)
    itype_new = 5;

  if (RunNumber>205000 && icalfac == 1){
    return(qhadcor_98_00_trans(&ITRA,&itype_new,&theta_deg,&ptj,&ej));
  }
  else if (icalfac == 2){
    return(qhadcor_98_00_trans(&ITRA,&itype_new,&theta_deg,&ptj,&ej));
  }  
  else {
    printf("no calib\n");
    return(qvar);
  }
}

Float_t TTools::qhadcor_98_00_trans(Int_t *itra,Int_t *itype,Float_t *jetpos, Float_t *ptin, Float_t *qvar)
{
// Translated qhadcor

    /* Initialized data */

    Float_t a_da98__[8] = { 0.879904,0.896318,0.906617,0.952872,0.955324,1.11211,1.14113,0.877522 };
    Float_t a_da99_tr__[8] = { 0.909524,0.964727,
            1.53204,0.980311,1.00269,1.14101,
            1.14013,0.984392 };
    Float_t b_da99_tr__[8] = { 0.,-1.56025,-.830946,-2.10923,0.,0.,0.,0. };
    Float_t c_da99_tr__[8] = { 5.48443,28.4622,
            406.365,10.539,3.07687,2.63095,
            3.14387,0. };
    Float_t a_ncem__[8] = { 0.924608,0.936066,1.01124,0.992175,0.958481,1.06982,1.09196,0.951766 };
    Float_t b_ncem__[8] = { 0.,-1.0471,-1.74525,-1.83088,0.,0.,0.,0. };
    Float_t c_ncem__[8] = { 5.60523,11.141,34.3981,19.3739,4.75546,5.05682,5.48667,4.65611 };
    Float_t a_ncem_tr__[8] = { 0.923356,0.935513,1.0688,1.12909,0.971356,1.09247,1.10269,1.4647 };
    Float_t b_ncem_tr__[8] = { 0.,-1.21237,-1.79435,-1.81835,0.,0.,0.,0. };
    Float_t c_ncem_tr__[8] = { 5.4633,11.7644,
            80.2606,186.33,2.62928,3.13916,
            3.11973,0.803558 };
    Float_t a_da9900__[8] = { 0.868093,0.877024,
            .921192,0.949944,0.977572,1.09807,
            1.17005,0.95448 };
    Float_t b_da98__[8] = { 0.,0.,0.,0.,0.,0.,0.,0. };
    Float_t b_da9900__[8] = { 0.,-1.08682,-1.70889,0.,0.,0.,-.175373,0. };
    Float_t c_da9900__[8] = { 5.65951,10.0507,
            17.0786,5.55913,5.84501,6.47199,
            8.74776,5.50889 };
    Float_t a_da9900_tr__[8] = { 0.868376,0.879944,
            1.151,0.960589,0.987846,1.15057,1.76864,
            0.915961 };
    Float_t b_da9900_tr__[8] = { 0.,-1.20845,
            -1.33753,-1.91097,-1.71629,-1.40015,
            -.691473,0. };
    Float_t c_da9900_tr__[8] = { 5.60348,10.699,
            196.481,9.64211,9.27777,13.1271,
            66.3346,2.53625 };
    Float_t a_ncep__[8] = { 0.89243,0.937066,0.949565,0.96273,0.934171,1.05762,1.06794,0.807267 };
    Float_t b_ncep__[8] = { 0.,-1.72168,-2.03856,-1.49543,0.,0.,0.,0. };
    Float_t c_ncep__[8] = { 4.95748,23.913,18.111,11.2861,4.50758,5.62712,5.91176,0.0 };
    Float_t a_ncep_tr__[8] = { 0.892234,0.909234,
            1.13882,0.981473,0.943926,1.07693,
            1.10456,0.877618 };
    Float_t b_ncep_tr__[8] = { 0.,-1.40548,-1.56322,-2.71634,0.,0.,-1.52819,0. };
    Float_t c_da98__[8] = { 4.47972,6.72565,5.85448,5.28593,3.44614,6.03909,6.58739,0. };
    Float_t c_ncep_tr__[8] = { 4.8741,11.5655,
            246.091,25.2448,0.,3.20811,9.36505,0. };
    Float_t a_da98_tr__[8] = { 0.878587,0.857764,
            0.885084,0.957587,0.979401,1.0804,
            1.05171,0.906565 };
    Float_t b_da98_tr__[8] = { 0.,0.,0.,0.,0.,0.,0.,0. };
    Float_t c_da98_tr__[8] = { 4.39525,0.,0.,
            2.63626,0.,0.,0.,0. };
    Float_t a_da99__[8] = { 0.911278,0.954436,0.993601,0.964861,1.00978,1.13403,1.15267,0.939532 };
    Float_t b_da99__[8] = { 0.,-1.4355,-1.65515,0.,0.,0.,0.,0. };
    Float_t c_da99__[8] = { 5.66306,22.6958,39.045,5.45322,6.05385,5.75632,6.33544,0. };

    /* System generated locals */
    Float_t ret_val;

    /* Local variables */
    Float_t scale=0, acorr=0, bcorr=0, ccorr=0;
    Int_t ith;
    Float_t ptinter;

/*=======================================================================
====*/

/* INPUT : */
/* ----- */
/*              ITRA = 0  for calibration without track information */
/*                     1   "     "        with     "         " */

/*              ITYPE = 1  for DATA 98-99 (nrun<233404) */
/*                      2  for DATA 99    (nrun>233404) */
/*                      3  for NC MC 98 e-p */
/*                      4  for DATA 99-00 */
/*                      5  for NC MC 99 e+p */

/*              jetpos = theta of the jet */

/*              ptin = transverse inmpulsion of the jet */

/*              qvar = energie or pt or px or py or pz */


/* OUTPUT : */
/* ------       qhadcor_98_00 = corrected varin */


/*===========================================================================*/
/* correction factors for  DATA 98/99 e-p (before online calibration change)*/
/*===========================================================================*/
    if (*itra != 0 && *itra != 1) {
        goto L101;
    }
    if (*itype != 1 && *itype != 2 && *itype != 3 && *itype != 4 && *itype != 
            5) {
        goto L101;
    }
    if (*ptin <= 0.) {
        goto L101;
    }
    ith = 0;
    if (*jetpos >= 0. && *jetpos < 10.) {
        ith = 1;
    }
    if (*jetpos >= 10. && *jetpos < 15.) {
        ith = 2;
    }
    if (*jetpos >= 15. && *jetpos < 30.) {
        ith = 3;
    }
    if (*jetpos >= 30. && *jetpos < 55.) {
        ith = 4;
    }
    if (*jetpos >= 55. && *jetpos < 80.) {
        ith = 5;
    }
    if (*jetpos >= 80. && *jetpos < 110.) {
        ith = 6;
    }
    if (*jetpos >= 110. && *jetpos < 135.) {
        ith = 7;
    }
    if (*jetpos >= 135. && *jetpos < 180.) {
        ith = 8;
    }
/* -------------------------------------------------------------- */
    if (*itype == 1) {
/* data 98/99 */
        if (*itra == 0) {
            if (ith != 0) {
                acorr = a_da98__[ith - 1];
                bcorr = b_da98__[ith - 1];
                ccorr = c_da98__[ith - 1];
            } else {
                goto L101;
            }
        }
        if (*itra == 1) {
            if (ith != 0) {
                acorr = a_da98_tr__[ith - 1];
                bcorr = b_da98_tr__[ith - 1];
                ccorr = c_da98_tr__[ith - 1];
            } else {
                goto L101;
            }
        }
    }
/* -------------------------------------------------------------- */
    if (*itype == 2) {
/* data 99 */
        if (*itra == 0) {
            if (ith != 0) {
                acorr = a_da99__[ith - 1];
                bcorr = b_da99__[ith - 1];
                ccorr = c_da99__[ith - 1];
            } else {
                goto L101;
            }
        }
        if (*itra == 1) {
            if (ith != 0) {
                acorr = a_da99_tr__[ith - 1];
                bcorr = b_da99_tr__[ith - 1];
                ccorr = c_da99_tr__[ith - 1];
            } else {
                goto L101;
            }
        }
    }
/* -------------------------------------------------------------- */
    if (*itype == 4) {
/* data 99/00 */
        if (*itra == 0) {
            if (ith != 0) {
                acorr = a_da9900__[ith - 1];
                bcorr = b_da9900__[ith - 1];
                ccorr = c_da9900__[ith - 1];
            } else {
                goto L101;
            }
        }
        if (*itra == 1) {
            if (ith != 0) {
                acorr = a_da9900_tr__[ith - 1];
                bcorr = b_da9900_tr__[ith - 1];
                ccorr = c_da9900_tr__[ith - 1];
            } else {
                goto L101;
            }
        }
    }
/* -------------------------------------------------------------- */
    if (*itype == 3) {
/* NC MC e-p */
        if (*itra == 0) {
            if (ith != 0) {
                acorr = a_ncem__[ith - 1];
                bcorr = b_ncem__[ith - 1];
                ccorr = c_ncem__[ith - 1];
            } else {
                goto L101;
            }
        }
        if (*itra == 1) {
            if (ith != 0) {
                acorr = a_ncem_tr__[ith - 1];
                bcorr = b_ncem_tr__[ith - 1];
                ccorr = c_ncem_tr__[ith - 1];
            } else {
                goto L101;
            }
        }
    }
/* -------------------------------------------------------------- */
    if (*itype == 5) {
/* NC MC e+p */
        if (*itra == 0) {
            if (ith != 0) {
                acorr = a_ncep__[ith - 1];
                bcorr = b_ncep__[ith - 1];
                ccorr = c_ncep__[ith - 1];
            } else {
                goto L101;
            }
        }
        if (*itra == 1) {
            if (ith != 0) {
                acorr = a_ncep_tr__[ith - 1];
                bcorr = b_ncep_tr__[ith - 1];
                ccorr = c_ncep_tr__[ith - 1];
            } else {
                goto L101;
            }
        }
    }
/* ================================================================ */
/*                    without track */
/* ================================================================ */
    if (*itra == 0) {

/* --- Constante for data 98 e-p : bin 8 */
/*               for data 99 e-p : bin 8 */
/*               for NC MC  e+p  : bin 8 */

        if (*itype == 1 && ith == 8 || *itype == 2 && ith == 8 || *itype == 5 
                && ith == 8) {
            scale = acorr;
            goto L99;
        }

/* --- Exponential function for all others -------------------------- 
*/

        scale = acorr * (1. - exp(bcorr - *ptin / ccorr));
        ptinter = *ptin / scale;
        scale = acorr * (1. - exp(bcorr - ptinter / ccorr));
    }
/* ================================================================ */
/*                 with track */
/* ================================================================ */
    if (*itra == 1) {

/* --- Constante for data 98 e-p  : if,of,cb3,cb2,cb1,bbe */
/*               for data 99 e-p  : bbe */
/*               for nc mc e+p    : cb3,bbe */
        if (*itype == 1 && (ith == 5 || ith == 7 || ith <= 8 || ith == 2 || 
                ith == 3 || ith == 6) || *itype == 2 && ith == 8 || *itype == 
                5 && (ith == 5 || ith == 8)) {
            scale = acorr;
            goto L99;
        }

/* --- Special function for nc mc e-p bbe */

        if (*itype == 3 && ith == 8) {
            scale = acorr / *ptin + ccorr;
            ptinter = *ptin / scale;
            scale = acorr / ptinter + ccorr;
            goto L99;
        }

/* --- Exponential function for all others -------------------------- 
*/

        scale = acorr * (1. - exp(bcorr - *ptin / ccorr));
        ptinter = *ptin / scale;
        scale = acorr * (1. - exp(bcorr - ptinter / ccorr));
    }
/* ================================================================ */
L99:
    ret_val = *qvar / scale;
    goto L100;
L101:
    ret_val = *qvar;
L100:
    return ret_val;
} /* qhadcor_98_00__ */



 
Bool_t TTools::ReadNoiseRuns(Char_t *InputFile)
{
// Read in noise runs

  // Def file pointer
  FILE *fpIn;

  // Def char. to read
  Char_t Ch; 
  Int_t cCounter = 0;
  Int_t EvNumber,RnNumber;
  Int_t MAXNOISEEVTS=1497;
  Int_t NoiseRuns[MAXNOISEEVTS];   // Noise Runs
  Int_t NoiseEvents[MAXNOISEEVTS]; // Noise events 

  // Check if input file is existing
  if ((fpIn = fopen(InputFile,"r")) == NULL) {
    printf("File %s not found\n",InputFile);
    return(kFALSE);
  }

  // Loop over all input numbers
  while ((Ch = fscanf(fpIn,"%d %d",&RnNumber,&EvNumber)) != EOF) {

    // For each line fill in event.NoiseRuns[]/event.NoiseEvents
    NoiseRuns[cCounter] = RnNumber;
    NoiseEvents[cCounter] = EvNumber; 

    // Increase counter
    if (cCounter < MAXNOISEEVTS)
      cCounter++;
    else {
      printf("ARRAY of NoiseRuns out of bounds => Please check validity\n");
      fclose(fpIn);
      return(kFALSE);
    }
  }
  // Close file
  fclose(fpIn);

  return(kTRUE);

}


Int_t TTools::OctogonIntersection(TVector3 ClusStart,Float_t Theta, Float_t Phi,
        Float_t RCurv,Float_t Zmin,Float_t Zmax,Float_t RBarr,TVector3 &XOut,Float_t &PhiOut)
{
// Calculer intersection entre la trace et le calorimetre 
Int_t IFlag=0;

Float_t Ph=0,DPh=0;
Float_t L8=0;
Float_t XOUT[2],X[2];

Float_t Dz=0;
Int_t ifl;

Float_t CosU;
Float_t CosV;
Float_t SinU;
Float_t Psi;
Float_t ZDelta;
Float_t RR;
Int_t N2Pi;


//--- Outer radius of octogon
 Float_t ROctog=RBarr/cos(TMath::Pi()/8.);

//--- Start poInt_tshould be inside octogone
if(ClusStart.Z()<=Zmin|| ClusStart.Z()>=Zmax) return IFlag;
Float_t R=sqrt(pow(ClusStart.X(),2)+pow(ClusStart.Y(),2));
if(R>=RBarr) {
        Ph=TMath::ATan2(ClusStart.Y(),ClusStart.X());
        if(Ph<0) Ph+=2.*TMath::Pi();
        L8=Ph/(TMath::Pi()/4.);
        DPh=Ph-L8*TMath::Pi()/4.;
        DPh=fabs(DPh-TMath::Pi()/8.);
        if(R*cos(DPh)>=RBarr) return IFlag;
}

//--- curvature radius
 Float_t RC=fabs(RCurv);

//--- sign of curvature defined to be 1 if phi increases with particle movement (otherwise -1.)
 Int_t IC=(Int_t) TMath::Sign(1,(Int_t)RCurv);

 Float_t CosT=cos(Theta);
 Float_t SinT=sin(Theta);

//--- unit vector from center of helix circle to start point
 Float_t UC[2],XC[2],VC[2],XX[2];
 UC[0]=cos(Phi-IC*TMath::Pi()/2.);
 UC[1]=sin(Phi-IC*TMath::Pi()/2.);
//--- center of helix circle
 XC[0]=ClusStart.X()-RC*UC[0];
 XC[1]=ClusStart.Y()-RC*UC[1];
//--- unit vector in direction of particle at start point
 VC[0]=cos(Phi);
 VC[1]=sin(Phi);

 if(R+2*RC<RBarr) goto jump100;
 if(fabs(SinT)<0.001) goto jump100;

//--- intersection with octogon in XY plane
 X[0]=ClusStart.X();
 X[1]=ClusStart.Y();
 OctgXY(X,Phi,IC,RC,ROctog,XOUT,ifl);
 XOut.SetX(XOUT[0]);
 XOut.SetY(XOUT[1]);
//  printf("Bad flag ... %d\n",ifl);
 if(ifl!=1) goto jump100;
//--- get turn angle of a particle (psi)
 XX[0]=XOut.X()-XC[0];
 XX[1]=XOut.Y()-XC[1];
 RR=sqrt(XX[0]*XX[0]+X[1]*XX[1]);
 XX[0]/=RR;
 XX[1]/=RR;
 CosU=UC[0]*XX[0]+UC[1]*XX[1];
 CosV=VC[0]*XX[0]+VC[1]*XX[1];
 SinU=CosV;
 Psi=TMath::ATan2(SinU,CosU);
//--- Range 0 -2*pi
 if(Psi<0) Psi+=2*TMath::Pi();

//--- distance to outer z plane
 if(CosT>=0) Dz=Zmax-ClusStart.Z();
 else Dz=ClusStart.Z()-Zmin;

//--- check where intersection occurs
 ZDelta=RC*Psi*CosT/SinT;
 if(fabs(ZDelta)>=Dz) goto jump100;
 XOut.SetZ(ClusStart.Z()+ZDelta);
 IFlag=1;
 goto jump200;


//--- intersection with z planes
jump100:
 if(CosT>0.0001) XOut.SetZ(Zmax);
 else if(CosT<-0.0001) XOut.SetZ(Zmin);
 else return IFlag;

//--- get turn angle of a particle
 if(RC>0.0001) {
        Dz=XOut.Z()-ClusStart.Z();
        Psi=Dz*SinT/(RC*CosT);
        Psi=fabs(Psi);
        N2Pi=(Int_t) (Psi/(2*TMath::Pi()));
        Psi=Psi-N2Pi*2*TMath::Pi();
 }else{
       Psi=0.;
 }

//--- leaving poInt_tXY plane
 XOut.SetX(XC[0]+RC*(cos(Psi)*UC[0]+sin(Psi)*VC[0]));
 XOut.SetY(XC[1]+RC*(cos(Psi)*UC[1]+sin(Psi)*VC[1]));

 IFlag=2;

jump200:
//--- range 0-4*Pi
 PhiOut=Phi+IC*Psi+(1-IC)*TMath::Pi();
//--- range 0-2*Pi
 if(PhiOut>2*TMath::Pi()) PhiOut-=2*TMath::Pi();


return IFlag;
}



void TTools::OctgXY(Float_t *xl,Float_t phi,Int_t ic,Float_t rc,
                                Float_t roctog,Float_t *xout,Int_t &iflag)
{
// Calc. X,Y

    Int_t i1, i2;


    /* Local variables */
    Float_t phil;
    Int_t l;
    Float_t r;
    Int_t l1, l2;
    Float_t xylin[2];
    Float_t philin, xycirc[2];
    Int_t ll1, ll2;
    Float_t xy1[2], xy2[2];
    Int_t ifl = 0;
    Float_t phi1, phi2;

/* *********************************************************************** */
/* PURPOSE: Intersection of a particle with an octogon in xy plane      * */
/*                                                                      * */
/* INPUT:   XL(2)  - start poInt_tinside octogon                         * */
/*          PHI    - gives direction of a particle at start poInt_t(rad.)* */
/*          IC     - sign of curvature defined to be 1 if phi increases * */
/*                   with particle movement (otherwise -1.)             * */
/*          RC     - curvature radius                             (cm)  * */
/*          ROCTOG - outer radius of octogon                            * */
/*                                                                      * */
/* OUTPUT:  XOUT(2)  - intersection (leaving) poInt_t                    * */
/*          IFLAG  = 1 if particle intersects octogon, otherwise =0     * */
/*                                                                      * */
/* AUTHOR:  V. Shekelyan                                at 14.11.91     * */
/* *********************************************************************** */

    iflag = 0;
/* find intersection of a line defined by XL(2) and PHI with an outer */
/* circle of octogon */
    xycirc[0] = 0.;
    xycirc[1] = 0.;
    //    printf("xl[0]=%f xl[1]=%f phi = %f circ[0] = %f circ[1] = %f r = %f xy1[0] = %f xy1[1] = %f xy2[0] = %f xy2[1] = %f ifl = %d\n",xl[0],xl[1],phi,xycirc[0],xycirc[1],roctog,xy1[0],xy1[1],xy2[0],xy2[1],ifl);
    LinCir(xl, phi, xycirc, roctog, xy1, xy2, ifl);
    if (ifl != 1) {
//     printf("First return \n");
    return;
    }
//     printf("test octxy ....... \n");

/* octant numbers of intersection of straight line from X with octogon */
/* ----- in particle direction */
    phi1 = TMath::ATan2(xy2[1], xy2[0]);
    if (phi1 < 0.) phi1 += 2*TMath::Pi();
    ll1 = (Int_t) (phi1 / (TMath::Pi()/4));
    if (ll1 > 7) ll1 = 7;

/* ----- in opposite direction */
    phi2 = TMath::ATan2(xy1[1], xy1[0]);
    if (phi2 < 0.) phi2 += 2*TMath::Pi();
    ll2 = (Int_t) (phi2 / (TMath::Pi()/4));
    if (ll2 > 7) ll2 = 7;

/* loop boundaries over sides of octogon from LL1 to LL2 */
/* in direction defined by a sign of curvature */
    l1 = ll1;
    l2 = ll2;
    if (ic > 0) {
        if (l2 < l1) {
            l2 = ll2 + 8;
        }
    } else {
        if (l2 > l1) {
            l1 = ll1 + 8;
        }
    }
/* center of particle circle */
    xycirc[0] = xl[0] + rc * cos(phi + ic * TMath::Pi()/2);
    xycirc[1] = xl[1] + rc * sin(phi + ic * TMath::Pi()/2);
/* try sides of octogon (first acceptable intersection is a solution) */
    i1 = l2;
    i2 = ic;
//     printf("l1=%d l2=%d %d %d\n",l1,l2,ll1,ll2);
    for (l = l1; i2 < 0 ? l >= i1 : l <= i1; l += i2) {
//   for (l = l1; l<=l2; l += ic) {
        phil = l * TMath::Pi()/4;
/* ------ define some poInt_tat octogon side */
        xylin[0] = roctog * cos(phil);
        xylin[1] = roctog * sin(phil);
/* ------ define direction of a side line */
        philin = phil + ic * TMath::Pi()/2 + TMath::Pi()/8;
        LinCir(xylin, philin, xycirc, rc, xout, xy2, ifl);
//      printf("LinCir : xout(%f,%f) ifl=%d\n",xout[0],xout[1],ifl);
        if (ifl != 1) {
//      printf("Badflag LinCir \n");
        continue;
        }
/* ------ check first intersection along side line direction */
/*        (with care about accuracy of calculations) */
        r = sqrt(xout[0] * xout[0] + xout[1] * xout[1]);
//      printf("%f roctog= %f \n",r,roctog);
        if (r < roctog + 0.001) goto L200;
    }
    return;

L200:
    iflag = 1;
    return;
}


void TTools::LinCir(Float_t *xl,Float_t phi,Float_t *xc,Float_t rc,
                        Float_t *xy1, Float_t *xy2, Int_t &iflag)
{
 Float_t cosp, sinp, a, b;
 Int_t isign;
 Float_t aa, bb, cc, det;
 Int_t iax1, iax2;

/* *********************************************************************** */
/* PURPOSE: Intersection of of a line with a circle in xy plane         * */
/*          equations:      (y-yc)**2+(x-xc)**2=rc**2                   * */
/*                          ax+by+c=0.                                  * */
/* INPUT:   XL(2),PHI - poInt_tat line and direction given by phi        * */
/*          XC(2),RC  - center of circle and radius                     * */
/*                                                                      * */
/* OUTPUT:  XY1(2),XY2(2) - intersection points (first and second one)  * */
/*                          ordered along line direction                * */
/*          IFLAG = 1 if there is intersection, otherwise =0            * */
/*                                                                      * */
/* AUTHOR:  V. Shekelyan                                at 14.11.91     * */
/* *********************************************************************** */

    iflag = 0;
    sinp = sin(phi);
    cosp = cos(phi);
/* define equation for a line */
    isign = 1;
    if (fabs(cosp) > fabs(sinp)) {
/* ------ y=ax+b */
        iax1 = 1;
        iax2 = 2;
        a = sinp / cosp;
        if (cosp < 0.) isign = -1;
    } else {
/* ------ x=ay+b */
        iax1 = 2;
        iax2 = 1;
        a = cosp / sinp;
        if (sinp < 0.) isign = -1;
    }
    b = xl[iax2-1] - a * xl[iax1-1];
/* solve equation aa*x(IAX1)**2+2.*bb*x(IAX1)+cc=0. */
/* which gives intersection points */
    aa = a * a + 1.;
    bb = a * (b - xc[iax2-1]) - xc[iax1-1];
    cc = (b -xc[iax2-1])*(b-xc[iax2-1])+xc[iax1-1]*xc[iax1-1]-rc*rc;
    det = bb*bb-aa*cc;
    if (det < 0.) return ;

    det = sqrt(det);
/* get intersection points ordered along line direction */
    xy1[iax1-1] = (-bb - isign * det) / aa;
    xy2[iax1-1] = (-bb + isign * det) / aa;
    xy1[iax2-1] = a * xy1[iax1-1] + b;
    xy2[iax2-1] = a * xy2[iax1-1] + b;

    iflag = 1;
    return ;
}




void TTools::CalcExtraPolTrack(TVector3 XStart,H1PartSelTrack *Track,Float_t &XOut,Float_t &YOut,Float_t &ZOut,Float_t &ThTrOut,Float_t &PhTrOut)
{
// Calculate track extrapolation into calo using octogon function
// Variables: - Track     => Track object to consider
//            - XOut      => X position extrapolated
//            - PhTrOut   => Extrapolated phi of track

//--- octogon extrapolation test
   TVector3 VecOut(0,0,0);
// 
//    // Get vertex of event
//    TVector3 XStart(0.,0.,0.);
//    if(Vertex[2] != -999999) {
//      XStart.SetX(Vertex[0]);
//      XStart.SetY(Vertex[1]);
//      XStart.SetZ(Vertex[2]);
//    }

   // Get phi of track
   TVector3 TrackEndPosition(0,0,0);
   Float_t PhiOut = 0.0;
   Float_t PhiTrk=Track->GetPhi();
   if(PhiTrk<0) PhiTrk+=2*TMath::Pi();

   // Calc curvature of track with nominal field (CONST !!!!!)
   Float_t RCurv= -1.0*Track->GetCharge()
        /((1/Track->GetPt())*BFIELDCONST * BFIELDNOM);

   // Calc intersection with octogon (if any)
   //   printf("x(z)=%f th=%f ph=%f r=%f z=%f z2=%f r=%f\n",XStart.Z(),Track->GetTheta(),PhiTrk,RCurv,zEDBBE,zIF1,rCalo);
      Int_t iflag=OctogonIntersection(XStart,Track->GetTheta(),PhiTrk,RCurv,zEDBBE,zIF1,rCalo,VecOut,PhiOut);

//  printf("OctoExt : (%f %f %f) flag=%d \n",XOut.X(),XOut.Y(),XOut.Z(),iflag);

  if(VecOut.Z()<=zEDBBE) {
        iflag=OctogonIntersection(XStart,Track->GetTheta(),PhiTrk,RCurv,zBACKDET,zIF1,rCalo,VecOut,PhiOut);
//   printf("OctoExt2 : (%f %f %f) flag=%d \n",VecOut.X(),VecOut.Y(),VecOut.Z(),iflag);
  }
//  printf("CylinderExt : (%f %f %f) \n",TrackEndPosition.X(),TrackEndPosition.Y(),TrackEndPosition.Z());


  // Final values
  TrackEndPosition=VecOut;
//  Double_t trackPhi   = TrackEndPosition.Phi(); 
  XOut = VecOut.X();  
  YOut = VecOut.Y();  
  ZOut = VecOut.Z();  
  PhTrOut = PhiOut;
}

---