//*-- Author : Jochen Markert 18.07.2007

#include "hparticletracksorter.h"
#include "hparticletool.h"
#include "hades.h"
#include "hruntimedb.h"
#include "hevent.h"
#include "hparticlecand.h"
#include "hemccluster.h"
#include "hparticledef.h"
#include "emcdef.h"
#include "htool.h"

#include <algorithm>
#include <iostream>
#include <iomanip>

using namespace std;
#define DEBUG 0

//_HADES_CLASS_DESCRIPTION
////////////////////////////////////////////////////////////////////////////////
//
//
// HParticleTrackSorter
//
// This Class loops over the HParticleCand category and
// analyze the objects. If the ECAL detector is used a loop over
// HEmcCluster is performed too to mark the clusters used by selected
// candidates. The Purpose of the procedure is to categorize
// the objects in the case of multiple usage of the single detector
// hits (RICH, InnerMDC, OuterMDC, META) as Double_t hit and to provide
// a hint which of the candidates is the best with respect to certain
// criteria (described by the bit flags below). If a HParticleCand
// object has been marked before by an other task with HParticleCand::kIsUsed it
// will be ignored by the sort process and stay unchanged. For the description
// of the differnt set/get functions of the HParticleCand objects see the
// documentation.
// To use HParticleTrackSorter without changing the flags of the HParticleCand objects
// (non persitent way inside a macro or reconstructor) one can call the
// backupFlags() function before and modification of the flags and
// restore the flags by calling restoreFlags() after the privat work has been
// done.
//
// For debugging purpose one cand set screen printouts to check the sorting
// (setDebug() + setPrintLevel(Int_t level)(level = 1-3)).
//
// The analysis of the HParticleCand objects is done in several steps:
//
// 1. Create the HParticleTrackSorter object and call the HParticleTrackSorter::Init()
//    To initialze the sorter object with category pointers etc ...
// 2. The internal structures are filled per event from the HParticleCand
//    category by calling the
//    Int_t HParticleTrackSorter::fill(Bool_t (*function)(HParticleCand* ))
//    function which expects a function pointer to a function taking an
//    HParticleCand pointer as argument and returning a Bool_t if the conditions on
//    the HParticleCand object are fullfilled (like for example in
//    Bool_t HParticleTrackSorter::selectHadrons(HParticleCand* pcand). For objects
//    failing the test the kIsRejected bit is set. The user can provide
//    his own selection functions to replace the build in selectLeptons(.....)
//    and selectHadrons(....) select functions.
//    The function pointer can be a pointer to a global function or member function
//    of an Object for example:
//
//    Bool_t select(HParticleCand* cand){             // global function
//             if ( put all selection criteria here ) return kTRUE;
//             else                                   return kFALSE;
//    } or
//
//    static Bool_t dummy::select(HParticleCand* cand){ // member function of object dummy
//                                                      // needs to be declared static !
//             if ( put all selection criteria here ) return kTRUE;
//             else                                   return kFALSE;
//    }
//    would be called in the code like
//
//    dummy d;
//    HParticleCand* p= new HParticleCand() // just to get an object
//    HParticleTrackSorter sorter;               // create the sorter object
//    // now we have 3 possibilities to call the
//    // test function
//    1. sorter.fill(select)        // global function
//    2. sorter.fill(dummy::select) // member function of object dummy (static call without object creation)
//    3. sorter.fill(d.select)      // member function of real object dummy
// 3. For all accepted objects the quality criteria and Double_t hit flags are
//    evaluated in within the selected sample and flagged to the HParticleCand objects.
//    The objects which fullfill the given minimum criteria the kIsAcceptedXXX bit is
//    flagged. For the RICH, InnerMDC, OuterMDC and META the kIsDoubleHitXXX bit are set
//    if the same detector hit index has been used more than once in the sample.
//    For the double hit rejection of Meta hits only the selected Meta hits
//    (ParticleCand::getSelectedMeta()) are taken into account. In the case of RPCClst+SHOWERHIT
//    The shower hit might be used by more than one candidate.
// 4. With the call of
//    Int_t HParticleTrackSorter::selectBest(Particle::ESwitch byQuality, Int_t byParticle)
//    the final selection on the objects is performed. The by the fill() procedure
//    selected objects are sorted by the selected quality criteria. Starting from the
//    best object one goes down the list and marks all objects with the kIsUsed and
//    kIsBestXXX bit if they do not reuse an already used detector hit index
//    (Double_t hit rejection).
//    For leptons the kIsLepton bit is flagged in addition. Individual detector hits
//    can be ignored in the selection procedure with the setIgnoreXXX(kTRUE) functions
//    (by default all hits are taken into account).
//    Particle::ESwitch byQuality can be used as quality criteria switch for
//
//          Particle::kIsBestHitRichSorter          (by number of pads)
//          Particle::kIsBestHitInnerMdcSorter      (by chi2, non fitted segments last)
//	    Particle::kIsBestHitOuterMdcSorter      (by chi2, non fitted segments last)
//          Particle::kIsBestHitMetaSorter          (by RK metamatch quality (or metaMatch quality if no RK , lower proirity))
//          Particle::kIsBestRKSorter               (by RK chi2, none fitted outer segments with lower priority)
//          Particle::kIsBestRKRKMETASorter         (by RK chi2 * RK META match quality)
//          Particle::kIsBestRKRKMETARadiusSorter   (by RK chi2 * RK META match radius)
//          Particle::kIsBestUserSorter             (by user provided sort function, no default function!)
//
//    user provide sort function:
//    static void   setUserSort(Bool_t    (*function)(candidateSort*, candidateSort*)
//    The user has to provide a function for sorting of the candidates:
//
//    Bool_t mysort(candidateSort* a, candidateSort* b) {
//      // this example does the same as kIsBestRKRKMETA
//      // sort in ascending order by RK_Chi2 * RK_META_match_Quality
//      // the function can use all data members of struct candidateSort defined in hparticletracksorter.h
//      return (a->RK_Chi2 * a->RK_META_match_Quality)  < (b->RK_Chi2 * b->RK_META_match_Quality);
//    }
//
//
// 5. The scheme allows for multiple selection following each other as
//    the already used objects are excluded from the next steps (for example first
//    select leptons and in the second iteration hadrons). The detector hit indicies
//    of the previous selection of kIsUsed marked objects are added to list of used
//    indicies by default not allowing to use a detector hit twice. To skip that
//    memory one can set setIgnorePreviousIndex(kTRUE).
// 6. After the analysis call HParticleTrackSorter::finalize(). After the selection
//    the user has to take into account only HParticleCand objects which  have
//    been flagged with kIsUsed. The objects slected by the Lepton selection will
//    be marked in addition with kIsLepton to identify them easily.
// 7. For lepton selections the internal selectLepton() function cand be manipulated.
//    With setRICHMatching(Particle::EMatch ,Float_t window) the matching between
//    inner MDC and RICH hit can be selected (default: HParticleTrackSorter::kUseRKRICHWindow,4 degree ).
// 8. For Better suppression of background the build in select functions can apply
//    a spacial match to the meta cell in y coordinate. setUseYMatching(kTRUE) (default)
//    makes use of a match in y in within 3 mm. ( see HParticleTool::isGoodMetaCell(HParticleCand* c,Float_t bound))
// 9. For Better suppression of background the build in select functions can require a valid beta measurement (beta>0)
//    setUseBeta(kTRUE) (default), for leptons in addition a low cut in beta can be set setBetaLetonCut(0.9) (default)
//10. For Better suppression of background the build in select functions can suppress fake measurement
//    setUseFakeRejection(kTRUE) (default)
//
//
//    A selection of unique lepton and hadron tracks could look like
//    the following example :
//
// void myMacro(Int_t nEvents,TString inputDir, TString inFile)
// {
//    Hades* myHades = new Hades;
//    gHades->setQuietMode(2);
//
//    HRootSource* source = new HRootSource();
//    source->setDirectory(((Text_t *)inputDir.Data()));
//    source->addFile((Text_t *)inFile.Data());
//
//    myHades->setDataSource(source);
//
//    if(!myHades->init()){
//	cout<<"Hades Init() failed!"<<endl;
//	exit(1);
//    }
//
//
//    //------------------------------------------------------------------------
//    // get category ointers and iterators
//    HCategory* pidCat  = (HCategory*)gHades->getCurrentEvent()->getCategory(catParticleCand);
//    HIterator* iter    = NULL;
//    if(pidCat)
//    {
//         iter = (HIterator *)pidCat->MakeIterator("native");
//    } else {
//        cout<<"Category Pointer is NULL!"<<endl;
//	exit(1);
//    }
//    //------------------------------------------------------------------------
//
//    --------------------------------------------------------------------------------------------
//    //--------------------------CONFIGURATION---------------------------------------------------
//    //At begin of the program (outside the event loop)
//    HParticleTrackSorter sorter;
//    //sorter.setDebug();                                            // for debug
//    //sorter.setPrintLevel(3);                                      // max prints
//    //sorter.setRICHMatching(Particle::kUseRKRICHWindowSorter,4.); // select matching RICH-MDC for selectLeptons() function
//    //sorter.setUseYMatching(kTRUE);                                // use meta cell build in select functions
//    //sorter.setMetaBoundary(4.0);                                   // match in 4 mm to meta cell
//    //sorter.setUseBeta(kTRUE);                                     // require beta >0 in build in select functions
//    //sorter.setTOFMAXCut( 60.);                                    // maximum tof [ns] allowd for match ( if beta is used)
//    //sorter.setBetaLeptonCut(0.9);                                 // lower beta cut for lepton select
//    //sorter.setUseFakeRejection(kTRUE);                            // reject fakes in  build in select functions (default kTRUE)
//    //sorter.setUseMETAQA(kTRUE);                                   // metaqa in  build in select functions (default kTRUE)
//    //sorter.setMETAQACut(3.0);                                     // metaqa cut for ahdrons+leptons
//    //sorter.setIgnoreInnerMDC();                                   // do not reject Double_t inner MDC hits
//    //sorter.setIgnoreOuterMDC();                                   // do not reject Double_t outer MDC hits
//    //sorter.setIgnoreMETA();                                       // do not reject Double_t META hits
//    //sorter.setIgnorePreviousIndex();                              // do not reject indices from previous selctions
//    sorter.init();                                                  // get catgegory pointers etc...
//    --------------------------------------------------------------------------------------------
//
//
//    //------------------------------------------------------------------------
//    // Event loop
//    for(Int_t i = 1;i < nEvents; i ++)
//    {
//        //----------break if last event is reached-------------
//        if(!gHades->eventLoop(1)) break;
//
//        //------------------------------------------------------------------------
//        // run your selection and set all flags
//
//        //------------------------------------------------------------------------
//        // clean vectors and index arrays
//        sorter.cleanUp();
//        //------------------------------------------------------------------------
//
//        // if you do not want to modify the flags
//        // for the following tasks etc (needs call of
//        // restoreFlags() at the end of the eventloop)
//        //sorter.backupFlags();
//
//        sorter.resetFlags(kTRUE,kTRUE,kTRUE,kTRUE);     // reset all flags for flags (0-28) ,reject,used,lepton
//        Int_t nCandLep     = sorter.fill(HParticleTrackSorter::selectLeptons);   // fill only good leptons
//        Int_t nCandLepBest = sorter.selectBest(Particle::kIsBestRKSorter,HParticleTrackSorter::kIsLeptonSorter);
//        Int_t nCandHad     = sorter.fill(HParticleTrackSorter::selectHadrons);   // fill only good hadrons (already marked good leptons will be skipped)
//        Int_t nCandHadBest = sorter.selectBest(Particle::kIsBestRKSorter,HParticleTrackSorter::kIsHadronSorter);
//        //------------------------------------------------------------------------
//
//        //------------------------------------------------------------------------
//        // analylize marked tracks (kIsUsed) ....
//    	  iter->Reset();
//        HParticleCand* pcand;
//        while((pcand = (HParticleCand*)iter->Next()) != NULL)
//	  {
//           // accept only objects marked as good
//	     if(!pcand->isFlagBit(HParticleCand::kIsUsed) ) continue;
//
//           // do other selection......
//        }
//        //------------------------------------------------------------------------
//
//        // if you do not want to modify the flags
//        // for the following tasks etc (needs call of
//        // backupFlags() at the beginning of the eventloop)
//        // sorter.restoreFlags();
//
//
//    }
//    //------------------------------------------------------------------------
//    //At the end of program
//    sorter.finalize();
// }
////////////////////////////////////////////////////////////////////////////////

ClassImp(sorter_setup)
ClassImp(HParticleTrackSorter)

Bool_t  HParticleTrackSorter::kDebug               = kFALSE;
Int_t   HParticleTrackSorter::printLevel           = 1;

Int_t   HParticleTrackSorter::kSwitchIndex         = 0;
Int_t   HParticleTrackSorter::kSwitchQuality       = 0;
Int_t   HParticleTrackSorter::kSwitchParticle      = 0;
Int_t   HParticleTrackSorter::kSwitchRICHMatching  = Particle::kUseRKRICHWindowSorter;
Float_t HParticleTrackSorter::fRICHMDCWindow       = 4.;

Bool_t  HParticleTrackSorter::kIsFieldOn           = kTRUE;
Bool_t  HParticleTrackSorter::kIgnoreRICH          = kFALSE;
Bool_t  HParticleTrackSorter::kIgnoreInnerMDC      = kFALSE;
Bool_t  HParticleTrackSorter::kIgnoreOuterMDC      = kFALSE;
Bool_t  HParticleTrackSorter::kIgnoreMETA          = kFALSE;
Bool_t  HParticleTrackSorter::kIgnorePreviousIndex = kFALSE;
Bool_t  HParticleTrackSorter::kUseYMatching        = kTRUE;
Bool_t  HParticleTrackSorter::kUseYMatchingScaling = kTRUE;
Float_t HParticleTrackSorter::fMetaBoundary        = 4.0;
Bool_t  HParticleTrackSorter::kUseBeta             = kTRUE;
Float_t HParticleTrackSorter::fBetaLepCut          = 0.9;
Bool_t  HParticleTrackSorter::kUseFakeRejection    = kTRUE;
Bool_t  HParticleTrackSorter::kUseMETAQA           = kTRUE;
Float_t HParticleTrackSorter::fMETAQACut           = 3.0;
Float_t HParticleTrackSorter::fTOFMAXCut           = 60.0;
Float_t HParticleTrackSorter::fDeltaSegCut         = 2.0;

Bool_t (*HParticleTrackSorter::pUserSort)(candidateSort*, candidateSort*) = NULL;

HParticleTrackSorter::HParticleTrackSorter(void)
    : TNamed("ParticleSorter","PID track candidate sorter")
{
    clear();
    setup.setDefault();
}

HParticleTrackSorter::HParticleTrackSorter(TString name, TString title)
: TNamed(name, title)
{
    clear();
    setup.setDefault();
}

HParticleTrackSorter::~HParticleTrackSorter(void)
{
    if(nameIndex)           delete [] nameIndex;
    if(nameQuality)         delete [] nameQuality;
    if(iterParticleCandCat) delete iterParticleCandCat;
}

void HParticleTrackSorter::clear()
{
    // initialize the member variables
    nt                   = NULL;
    isSimulation         = kFALSE;
    nameIndex            = 0;
    nameQuality          = 0;
    pParticleCandCat     =  NULL;
    pEmcClusterCat       =  NULL;
    iterParticleCandCat  =  NULL;
    fill_Iteration       =  0;
    selectBest_Iteration =  0;
    currentEvent         = 100000000;
    fEvent               = NULL;
}

Bool_t HParticleTrackSorter::init(HRecEvent* evt)
{
    // get the pointer to the HParticleCand categegory and create the
    // corresponding iterator. HRecEvent* evt == 0 (default) the current event
    // will be used otherise the event provided.

    if(evt) { fEvent = evt; }
    else    { fEvent = (HRecEvent*)gHades->getCurrentEvent(); }

    pParticleCandCat = fEvent->getCategory(catParticleCand);
    if (!pParticleCandCat) {
	Error("HParticleTrackSorter::init()"," No HParticleCand Input -> Switching HParticleTrackSorter OFF");
    }
    else {
        if(iterParticleCandCat) delete iterParticleCandCat;
	iterParticleCandCat = (HIterator *) pParticleCandCat->MakeIterator();
    }
    pEmcClusterCat   = fEvent->getCategory(catEmcCluster);

    HCategory* catKine = fEvent->getCategory(catGeantKine);
    if(catKine) isSimulation = kTRUE;

    if(kDebug) {

	Int_t size;
	nameIndex   = HTool::parseString("kIsIndexRICH,kIsIndexInnerMDC,kIsIndexOuterMDC,"
					 "kIsIndexTOF,kIsIndexSHOWER,kIsIndexRPC,kIsIndexMETA",size,",",kFALSE);
        cout<<"-------------------------------------------"<<endl;
	cout<<"HParticleTrackSorter::init() :"<<endl;
	cout<<"index :: " <<endl;
	for(Int_t i = 0; i < size; i ++){
	    cout<<i <<" "<<nameIndex[i].Data()<<endl;
	}
	nameQuality = HTool::parseString("kIsBestHitRICH,kIsBestHitInnerMDC,kIsBestHitOuterMDC,kIsBestHitMETA,"
					 "kIsBestRK,kIsBestRKRKMETA,kIsBestRKRKMETARadius",size,",",kFALSE);
        for(Int_t i = 0; i < size; i ++){
	    cout<<i + SORTER_NUMBER_OF_INDEX<<" "<<nameQuality[i].Data()<<endl;
	}
	cout<<"-------------------------------------------"<<endl;
    }

     return kTRUE;
}

Bool_t HParticleTrackSorter::finalize(void)
{
    // clean up the memory and
    // write ntuple to file (if it exist)
    cleanUp(kTRUE);
    return kTRUE ;
}

Bool_t  HParticleTrackSorter::isDeltaSegOk(HParticleCand* cand)
{
    // returns kTRUE if inner and outer segment
    // agree in theta, phi inside fDeltaSegCut

    if(cand->getTheta2() ==-1) return kFALSE;

    if((fabs(cand->getPhi()   - cand->getPhi2())   < fDeltaSegCut) &&
       (fabs(cand->getTheta() - cand->getTheta2()) < fDeltaSegCut)
      ) return kTRUE;

    return kFALSE;

}


Int_t HParticleTrackSorter::clearVector( vector<candidateSort*>& all_candidateSorts)
{
    // delete all candidate structs in vector
    // from heap and call the clear function of
    // vector afterwards.
    for(UInt_t i = 0; i < all_candidates.size() ; ++ i)
    {
	candidateSort* cand = all_candidates[i];
        delete cand;
    }
    all_candidates.clear();
    return 0;
}


Bool_t HParticleTrackSorter::cmpIndex(candidateSort* a, candidateSort* b)
{
    // sorts the candidate list ascending by the index
    // of the detector hit of interest.

    switch (kSwitchIndex) {
    case Particle::kIsIndexRICHSorter:
	// sort in ascending order by Rich index
	return a->ind_RICH     < b->ind_RICH;
    case Particle::kIsIndexInnerMDCSorter:
	// sort in ascending order by inner MDC segment index
	return a->ind_innerMDC < b->ind_innerMDC;
    case Particle::kIsIndexOuterMDCSorter:
	// sort in ascending order by outer MDC segment index
	return a->ind_outerMDC < b->ind_outerMDC;
    case Particle::kIsIndexTOFSorter:
	// sort in ascending order by TOF index
	return a->ind_TOF      < b->ind_TOF;
    case Particle::kIsIndexSHOWERSorter:
	// sort in ascending order by SHOWER index
	return a->ind_SHOWER   < b->ind_SHOWER;
    case Particle::kIsIndexRPCSorter:
	// sort in ascending order by RPC index
	return a->ind_RPC      < b->ind_RPC;

    default:
	::Error("HParticleTrackSorter::CmpIndex()","Unknown switch statement %i!",kSwitchIndex);
	return kFALSE;

    }
}

Bool_t HParticleTrackSorter::cmpQuality(candidateSort* a, candidateSort* b)
{
    // sorts the candidate list by the quality critera
    // of interest (ascending or decending depending on
    // the quantity). The best object is always the first.

    switch (kSwitchQuality) {
    case Particle::kIsBestHitRICHSorter:
	// sort in ascending order by inner MDC chi2
	return a->nPadsRich             < b->nPadsRich;
    case Particle::kIsBestHitInnerMDCSorter:
	// sort in ascending order by inner MDC chi2
	return a->innerMDCChi2          < b->innerMDCChi2;
    case Particle::kIsBestHitOuterMDCSorter:
	// sort in ascending order by outer MDC chi2
	return a->outerMDCChi2          < b->outerMDCChi2;
    case Particle::kIsBestHitMETASorter:
	// sort in ascending order by RK METAMATCH QUALITY
	return a->RK_META_match_Quality < b->RK_META_match_Quality;
    case Particle::kIsBestRKSorter:
        // sort in ascending order by RK chi2
	return a->RK_Chi2               < b->RK_Chi2;
    case Particle::kIsBestRKRKMETASorter:
	// sort in ascending order by RK_Chi2 * RK_META_match_Quality
        return (a->RK_Chi2 * a->RK_META_match_Quality)  < (b->RK_Chi2 * b->RK_META_match_Quality);
    case Particle::kIsBestRKRKMETARadiusSorter:
	// sort in ascending order by RK_Chi2 * RK_META_match_Quality
        return (a->RK_Chi2 * a->RK_META_match_Radius)  < (b->RK_Chi2 * b->RK_META_match_Radius);
    case Particle::kIsBestUserSorter:
	// sort in ascending order by user function
	if(pUserSort == NULL) {
            ::Error("HParticleTrackSorter::CmpQuality()","kIsBestUser is set but user function pointer i NULL!");
	    return kTRUE;
	}
        else return (*pUserSort)(a,b);
    default:
	::Error("HParticleTrackSorter::CmpQuality()","Unknown switch statement %i!",kSwitchQuality);
        return kFALSE;

    }
}


Bool_t HParticleTrackSorter::rejectIndex(candidateSort* a, Particle::ESwitch byIndex, Int_t& index)
{
    // if function return kTRUE the Object will be rejected
    // rejects candidates from the actual list, where the
    // detector hit of interest is not defined.

    switch (byIndex) {
    case Particle::kIsIndexRICHSorter:
	// NO RICH
	return ( index = a->ind_RICH )     < 0;
    case Particle::kIsIndexInnerMDCSorter:
	// NO inner MDC segment
	return ( index = a->ind_innerMDC ) < 0;
    case Particle::kIsIndexOuterMDCSorter:
	// NO outer MDC segment
	return ( index = a->ind_outerMDC ) < 0;
    case Particle::kIsIndexTOFSorter:
	// NO TOF
	return ( index = a->ind_TOF )      < 0;
    case Particle::kIsIndexSHOWERSorter:
	// NO SHOWER
	return ( index = a->ind_SHOWER )   < 0;
    case Particle::kIsIndexRPCSorter:
	// NO RPC
	return ( index = a->ind_RPC )      < 0;
    default:
	::Error("HParticleTrackSorter::rejectIndex()","Unknown switch statement %i!",byIndex);
        return kFALSE;
    }
}
Bool_t HParticleTrackSorter::rejectQuality(candidateSort* a, Particle::ESwitch byQuality)
{
    // if function return kTRUE the Object will be rejected
    // rejects the candidates from the actual list if the quality
    // of interest does not match the mimnum requirement.

    switch (byQuality) {
    case Particle::kIsBestHitRICHSorter:
	// RICH ring
	return a->ind_RICH < 0;
    case Particle::kIsBestHitInnerMDCSorter:
	// NO fitted inner MDC segments
	return a->innerMDCChi2 < 0;
    case Particle::kIsBestHitOuterMDCSorter:
	// NO fitted outer MDC segments
	return (a->outerMDCChi2 < 0 || a->ind_outerMDC < 0);
    case Particle::kIsBestHitMETASorter:
	// NO META
	return (a->ind_TOF < 0 && a->ind_SHOWER < 0 && a->ind_RPC < 0);
    case Particle::kIsBestRKSorter:
	// NO outer MDC segment or NO RK
	return (a->RK_Chi2 < 0);
    case Particle::kIsBestRKRKMETASorter:
	// NO META or NO RK
        return (a->ind_TOF < 0 && a->ind_SHOWER < 0 && a->ind_RPC < 0) || (a->RK_Chi2 < 0);
    case Particle::kIsBestRKRKMETARadiusSorter:
	// NO META or NO RK
        return (a->ind_TOF < 0 && a->ind_SHOWER < 0 && a->ind_RPC < 0) || (a->RK_Chi2 < 0);
    case Particle::kIsBestUserSorter:
	// NO META or NO RK
        return (a->ind_TOF < 0 && a->ind_SHOWER < 0 && a->ind_RPC < 0) || (a->RK_Chi2 < 0);

    default:
	::Error("HParticleTrackSorter::rejectQuality()","Unknown switch statement %i!",byQuality);
        return kFALSE;

    }
}

Int_t HParticleTrackSorter::flagAccepted(vector <candidateSort*>& all_candidates, Particle::ESwitch byQuality)
{
    // flag all candidates which does satisfy the quality
    // criteria from input with kIsAcceptedXXX bit. If kDebug
    // is set and printLevel >=2 the statistics for the accepted
    // objects is printed.


    Int_t ctAll           = 0;
    Int_t ctAcceptQuality = 0;
    HParticleCand* pcand;
    for(UInt_t i = 0; i < all_candidates.size(); ++ i)
    {
        ++ ctAll;
	candidateSort* cand = all_candidates[i];
	if(rejectQuality(cand, byQuality))      continue;
	++ ctAcceptQuality;

	pcand = dynamic_cast<HParticleCand*>(pParticleCandCat->getObject(cand->ind_Cand));
	if(!pcand) {
	    Error("HParticleTrackSorter::flagAccepted()","NULL pointer retrieved for HParticleCand at index %i",cand->ind_Cand);
	    continue;
	}
	switch (byQuality) {
	case Particle::kIsBestHitRICHSorter:
	    pcand -> setFlagBit(Particle::kIsAcceptedHitRICH);
            //------------------------------------------------
            // check the ring correlation
	    if(pcand->getChi2() >= 0){
		if( pcand->getRingCorr() == Particle::kIsRICHRK ||
		    pcand->getRingCorr() == (Particle::kIsRICHRK|Particle::kIsRICHMDC)
		  ) pcand->setFlagBit(Particle::kIsAcceptedRKRICH);
	    }
	    if( pcand->getRingCorr() == Particle::kIsRICHMDC ||
	        pcand->getRingCorr() == (Particle::kIsRICHRK|Particle::kIsRICHMDC)
	      ) pcand->setFlagBit(Particle::kIsAcceptedHitRICHMDC);
            //------------------------------------------------
	    break;
	case Particle::kIsBestHitInnerMDCSorter:
	    pcand -> setFlagBit(Particle::kIsAcceptedHitInnerMDC);
	    break;
	case Particle::kIsBestHitOuterMDCSorter:
	    pcand -> setFlagBit(Particle::kIsAcceptedHitOuterMDC);
	    break;
	case Particle::kIsBestHitMETASorter:
	    pcand -> setFlagBit(Particle::kIsAcceptedHitMETA);
            //------------------------------------------------
            // check if it was accepted for RK
	    if(pcand->getChi2() >= 0) {
		pcand-> setFlagBit(Particle::kIsAcceptedRKMETA);
	    }
            //------------------------------------------------
	    break;
	case Particle::kIsBestRKSorter:
	    pcand -> setFlagBit(Particle::kIsAcceptedRK);
	    break;
	case Particle::kIsBestRKRKMETASorter:
	    pcand -> setFlagBit(Particle::kIsAcceptedRKRKMETA);
	    break;
	case Particle::kIsBestRKRKMETARadiusSorter:
	    pcand -> setFlagBit(Particle::kIsAcceptedRKRKMETARadius);
	    break;
	case Particle::kIsBestUserSorter:
	    pcand -> setFlagBit(Particle::kIsAcceptedUser);
	    break;
	default :

	    Error("HParticleTrackSorter::flagAccepted()","Unknown switch statement : %i!",byQuality);
	    exit(1);
	    break;
	}

    }
    if(kDebug && printLevel >= 2){
	cout<<"flagAccepted(): "
	    <<"s: "                  <<selectBest_Iteration
	    <<", f: "                <<fill_Iteration
	    <<", kSwitchParticle: "  <<kSwitchParticle
	    <<", quality "           <<nameQuality[byQuality-SORTER_NUMBER_OF_INDEX].Data()
	    <<", candidate= "        <<ctAll
	    <<", acc. quality = "    <<ctAcceptQuality
	    <<endl;

    }
    return 0;
}

Int_t HParticleTrackSorter::flagDouble(vector <candidateSort*>& all_candidates, Particle::ESwitch byIndex)
{
    // sort vector in ascending order by the particular index byindex. Loops
    // over the vector and extract objects with same detetector index and marks them
    // as Double_t hits in HParticleCand (kIsDoubleHitXXX bit). Non valid detecor
    // hits (index -1) are not taken into account. If kDebug
    // is set and printLevel >=2 the statistics for the accepted
    // objects is printed.

    kSwitchIndex   = byIndex;
    sort(all_candidates.begin(), all_candidates.end(), this->cmpIndex);


    Int_t index;
    Int_t old_index =-99;
    vector <candidateSort*> daughter;

    Int_t ctAll          = 0;
    Int_t ctAcceptIndex  = 0;
    Int_t ctAcceptDoubles= 0;
    Int_t ctCheck        = 0;

    for(UInt_t i = 0; i < all_candidates.size(); ++ i)
    {
        ++ ctAll;
	candidateSort* cand = all_candidates[i];
	if(rejectIndex  (cand, byIndex, index)) continue;
        ++ ctAcceptIndex;
	if(index == old_index || old_index == -99)
	{ // first time  or  still same index
	    daughter.push_back(cand);
	} else {
            // finalize list
	    if(!daughter.empty())
	    {
                setFlagsDouble(daughter,byIndex);
		if(daughter.size() >= 2){
		    ctAcceptDoubles += daughter.size();
		}
		ctCheck += daughter.size();
		daughter.clear();
	    }
	    daughter.push_back(cand);
	}
	old_index = index;
    }
    // finalize list last group
    if(!daughter.empty())
    {
	setFlagsDouble(daughter, byIndex);
	if(daughter.size() >= 2){
	    ctAcceptDoubles += daughter.size();
	}
	ctCheck += daughter.size();
    }
    if(kDebug && printLevel >= 2){
	cout<<"flagDouble(): "
	    <<"s: "                <<selectBest_Iteration
	    <<", f: "              <<fill_Iteration
	    <<", kSwitchParticle: "<<kSwitchParticle
	    <<", index "           <<nameIndex  [byIndex].Data()
	    <<", candidate= "      <<ctAll
	    <<", acc. index = "    <<ctAcceptIndex
            <<", acc. doubles = "  <<ctAcceptDoubles
            <<", check = "         <<ctCheck
	    <<endl;

    }
    return 0;

}

Int_t HParticleTrackSorter::setFlagsDouble(vector<candidateSort*>& daughter, Particle::ESwitch byIndex)
{
    // flag the HParticleCand objects with the rank of the Double_t hit criteria as output
    // form the sorting algorithms. Expects as input a vector of candidates with the
    // same hit index.

    if(daughter.empty()   ) return 0;    // nothing to do
    if(daughter.size() < 2) return 0;    // no Double_t hit


    HParticleCand* cand;
    for(UInt_t i = 0; i < daughter.size(); ++ i)
    {  // loop over list and set flags
	cand = dynamic_cast<HParticleCand*>(pParticleCandCat->getObject(daughter[i]->ind_Cand));
	if(!cand) {
	    Error("HParticleTrackSorter::setFlagsDouble()","NULL pointer retrieved for HParticleCand at index %i",daughter[i]->ind_Cand);
	    continue;
	}
	switch (byIndex)
	{
	case Particle::kIsIndexRICHSorter:
	    cand -> setFlagBit(Particle::kIsDoubleHitRICH);
	    break;
	case Particle::kIsIndexInnerMDCSorter:
	    cand -> setFlagBit(Particle::kIsDoubleHitInnerMDC);
	    break;
	case Particle::kIsIndexOuterMDCSorter:
	    cand -> setFlagBit(Particle::kIsDoubleHitOuterMDC);
	    break;
	case Particle::kIsIndexTOFSorter:
	    cand -> setFlagBit(Particle::kIsDoubleHitMETA);
	    break;
	case Particle::kIsIndexSHOWERSorter:
	    cand -> setFlagBit(Particle::kIsDoubleHitMETA);
	    break;
	case Particle::kIsIndexRPCSorter:
	    cand -> setFlagBit(Particle::kIsDoubleHitMETA);
	    break;
	default :

	    Error("HParticleTrackSorter::setFlagsDouble()","Unknown switch statement : %i!",byIndex);
	    exit(1);
	    break;

	}
    }

    return 0;
}

Bool_t HParticleTrackSorter::flagEmcClusters()
{
    // loop over EmcCluster cat and reset IsUsedInParticleCand.
    // loop over all candidates and flag emc clusters used in
    // selected candidates only
    if(pEmcClusterCat&&pParticleCandCat)
    {
        HEmcCluster* emc = 0;
	Int_t n = pEmcClusterCat->getEntries();
	for(Int_t i = 0; i < n; i ++){
            emc = (HEmcCluster*) pEmcClusterCat->getObject(i);
            if(emc) emc -> unsetIsUsedInParticleCand();
	}

        HParticleCand* cand = 0;
	n = pParticleCandCat->getEntries();
	for(Int_t i = 0; i < n; i ++){
	    cand = dynamic_cast<HParticleCand*>(pParticleCandCat->getObject(i));
	    if(cand->isFlagBit(Particle::kIsUsed) ){
		Int_t ind = cand->getEmcInd();
		if(ind>=0){
		    emc = (HEmcCluster*) pEmcClusterCat->getObject(ind);
                    if(emc)  emc->setIsUsedInParticleCand();
		}
	    }
	}
    }
     return kTRUE;
}


Int_t HParticleTrackSorter::fillInput(vector < candidateSort* >& all_candidates)
{
    // loop over HParticleCand category and fill a vector of temporary objects
    // which will be used for sorting and flagging of the objects in the end.
    // The candidate objects are initialzied with values of interest in a well
    // defined way, which makes the sorting and rejecting of objects later on easier.

    HParticleCand* pCand ;

    candidateSort* cand;
    iterParticleCandCat->Reset();
    while ((pCand = dynamic_cast<HParticleCand*>(iterParticleCandCat->Next())) != 0 )
    {

	if(pCand->isFlagBit(Particle::kIsUsed))
	{
	    if(!kIgnorePreviousIndex && fill_Iteration == 2)
	    {
		//-------------------------------------------------------------------
		// add detector indices of already used objects to
		// the current index list
		if(pCand->getRichInd()     >= 0 &&
		   (find(index_RICH.begin(),index_RICH.end()        ,pCand->getRichInd())     == index_RICH.end()) ){
		    index_RICH.push_back(pCand->getRichInd());
		}
		if(pCand->getInnerSegInd() >= 0 &&
		   (find(index_InnerMDC.begin(),index_InnerMDC.end(),pCand->getInnerSegInd()) == index_InnerMDC.end()) ){
		    index_InnerMDC.push_back(pCand->getInnerSegInd());
		}
		if(pCand->getOuterSegInd() >= 0 &&
		   (find(index_OuterMDC.begin(),index_OuterMDC.end(),pCand->getOuterSegInd()) == index_OuterMDC.end()) ){
		    index_OuterMDC.push_back(pCand->getOuterSegInd());
		}

		//-------------------------------------------------------------------
                // TOF Clst + TOF hit special case
		if(pCand->getSelectedMeta() == Particle::kTofClst && pCand->getTofClstInd()      >= 0 &&
		   (find(index_TOFClst.begin(),index_TOFClst.end()  ,pCand->getTofClstInd())  == index_TOFClst.end()) ){
		    index_TOFClst.push_back(pCand->getTofClstInd());
		}
		if((pCand->getSelectedMeta() == Particle::kTofHit1 || pCand->getSelectedMeta() == Particle::kTofHit2) &&
		   pCand->getTofHitInd()       >= 0 &&
		   (find(index_TOFHit.begin(),index_TOFHit.end()    ,pCand->getTofHitInd()+1000)  == index_TOFHit.end()) ){
		    index_TOFHit.push_back(pCand->getTofHitInd()+1000);
		}
		//-------------------------------------------------------------------
		if(pCand->getSelectedMeta() == Particle::kShowerHit && pCand->getShowerInd()   >= 0 &&
		   (find(index_SHOWER.begin(),index_SHOWER.end()    ,pCand->getShowerInd())   == index_SHOWER.end()) ){
		    index_SHOWER.push_back(pCand->getShowerInd());
		}

		if(pCand->getSelectedMeta() == Particle::kRpcClst && pCand->getRpcInd()      >= 0 &&
		   (find(index_RPC.begin(),index_RPC.end()          ,pCand->getRpcInd())      == index_RPC.end()) ){
		    index_RPC.push_back(pCand->getRpcInd());
		}

	    }
	    //-------------------------------------------------------------------

	    continue; // skip objects which are marked
	}
	if(pCand->isFlagBit(Particle::kIsRejected))  continue; // skip objects which are marked

	for(Int_t i = 0; i < 29; i ++) { pCand->unsetFlagBit(i);  }  // clean condition, only reject/used/lepton bit stays!

        //-------------------------------------------------------------------
	// fill the structur objects with needed
        // information

	cand = new candidateSort;
        cand->cand                  = pCand;
	cand->ind_Cand              = pCand->getIndex();

	cand->ind_RICH              = pCand->getRichInd();
	cand->ind_innerMDC          = pCand->getInnerSegInd();
	cand->ind_outerMDC          = pCand->getOuterSegInd();
	cand->ind_TOF               = (pCand->getTofClstInd() >= 0)? pCand->getTofClstInd() : (pCand->getTofHitInd() >= 0)? pCand->getTofHitInd() + 1000 : -1 ; // special case tof
	cand->ind_SHOWER            = pCand->getShowerInd();
	cand->ind_RPC               = pCand->getRpcInd();
	cand->selectedMeta          = pCand->getSelectedMeta();
        cand->nPadsRich             = pCand->getRingNumPads();
	cand->innerMDCChi2          = pCand->getInnerSegmentChi2();
	cand->outerMDCChi2          = pCand->getOuterSegmentChi2();
                                      
	cand->sec                   =  pCand->getSector();
	cand->beta                  =  pCand->getBeta();
	cand->mom                   =  pCand->getMomentum();
	cand->MDC_dEdx              =  pCand->getMdcdEdx();
	cand->nMdcLayerInnerSeg     =  pCand->getNLayer(0);
        cand->nMdcLayerOuterSeg     =  pCand->getNLayer(1);
	cand->RK_Chi2               = (pCand->getOuterSegmentChi2() < 0 && pCand->getChi2() >= 0 )? pCand->getChi2() + 100 : pCand->getChi2(); // lower priority sorting
        cand->RICH_RK_Corr          =  pCand->getRingCorr();
	cand->RK_RICH_match_Quality = (pCand->getRichInd() >= 0)?  pCand->getRichMatchingQuality(): -1;
	cand->RK_META_match_Quality = (finite(pCand->getMetaMatchQuality()))? pCand->getMetaMatchQuality() : -1;
	cand->RK_META_match_Radius  = (finite(pCand->getMetaMatchRadius()))? pCand->getMetaMatchRadius() : -1;
	cand->RK_META_dx            = (finite(pCand->getRkMetaDx()))? pCand->getRkMetaDx() : -1;
	cand->RK_META_dy            = (finite(pCand->getRkMetaDy()))? pCand->getRkMetaDy() : -1;
	cand->RK_META_match_Shr_Quality = pCand->getMetaMatchQualityShower();
        cand->innerMDCTheta = pCand->getTheta();
        cand->innerMDCPhi   = pCand->getPhi();
        cand->outerMDCTheta = pCand->getTheta2();
        cand->outerMDCPhi   = pCand->getPhi2();

	cand->userSort1=0;
	cand->userSort2=0;
	cand->userSort3=0;


	if(pCand->getSelectedMeta() == kTofClst ||
	   pCand->getSelectedMeta() == kTofHit1 ||
	   pCand->getSelectedMeta() == kRpcClst
	  ) {
	    cand->MetaModule = pCand->getMetaModule(0);
	    cand->MetaCell   = pCand->getMetaCell(0);
	} else if (pCand->getSelectedMeta() == kTofHit2){
	    cand->MetaModule = pCand->getMetaModule(1);
	    cand->MetaCell   = pCand->getMetaCell(1);
	} else {
	    cand->MetaModule = -1;
            cand->MetaCell   = -1;
	}

	if(cand->RK_Chi2 < 0){
	    cand->RK_META_match_Quality += 100; // lower priority in sorting
	    cand->RK_META_match_Radius  += 1000; // lower priority in sorting
	}

	//-------------------------------------------------------------------


	all_candidates.push_back(cand);
    }
    return 0;
}

void HParticleTrackSorter::printCand(candidateSort* cand, Int_t i,TString spacer)
{
    // print one candidate object and the flags which have been
    // already set to HParticleCand
    cout<<spacer.Data()<<i<<" --------------------------------------------"<<endl;
    cout<<spacer.Data()<<"cand " <<cand->ind_Cand
	<<", R: "  <<cand->ind_RICH
	<<", iM: " <<cand->ind_innerMDC
	<<", oM: " <<cand->ind_outerMDC
	<<", T: "  <<cand->ind_TOF
	<<", S: "  <<cand->ind_SHOWER
	<<", RPC: "<<cand->ind_RPC
	<<", selMeta: "<<cand->selectedMeta
	<<endl;
    cout<<spacer.Data()
	<<"Npads:"    <<cand->nPadsRich
	<<", ichi2: " <<cand->innerMDCChi2
	<<", ochi2: " <<cand->outerMDCChi2
	<<", RKM: "   <<cand->RK_META_match_Quality
	<<", RKMSHR: "<<cand->RK_META_match_Shr_Quality
	<<", RK: "    <<cand->RK_Chi2
	<<", RKRKMETA: "  <<(cand->RK_Chi2 != -1 && cand->RK_META_match_Quality != -1 ? cand->RK_Chi2*cand->RK_META_match_Quality : -1)
	<<", RKRKMETARadius: "  <<(cand->RK_Chi2 != -1 && cand->RK_META_match_Radius != -1 ? cand->RK_Chi2*cand->RK_META_match_Radius : -1)
	<<", isRingRK: "<<(Int_t) cand->RICH_RK_Corr
	<<endl;
    HParticleCand* pcand = dynamic_cast<HParticleCand*>(pParticleCandCat->getObject(cand->ind_Cand));
    if(cand)
    {
	cout<<spacer.Data()<<flush;
	pcand->printFlags();
    } else {
	::Error("HParticleTrackSorter::printCand()","NULL pointer retrieved for HParticleCand Index %i",cand->ind_Cand);

    }

}
void HParticleTrackSorter::printEvent(TString comment)
{
    Int_t ctcand = 0;

    cout<<"----------------------------------------"<<endl;
    cout<<"Event kSwitchParticle: "<<kSwitchParticle<<" "
 	<<setw(6)                  <<fEvent->getHeader()->getEventSeqNumber()
	<<"s: "                    <<selectBest_Iteration
	<<", f: "                  <<fill_Iteration
	<<" "                      <<comment.Data()
	<<endl;
    HParticleCand* pcand;
    iterParticleCandCat->Reset();
    while ((pcand = dynamic_cast<HParticleCand*>(iterParticleCandCat->Next())) != 0 )
    {
	cout <<setw(6)<<ctcand<<" "<<flush;
	pcand->print();
	++ ctcand;
    }
}
void HParticleTrackSorter::printSetup()
{
    // print seeting of this sorter:
    // column 1 : backup
    // column 2 : setup
    // column 3 : actual

    cout<<"--------------------------------------------------------------------"<<endl;
    cout<<"HParticleTrackSorter: evt "<<gHades->getEventCounter()<<endl;
    cout<<"                        backup  setup  actual"<<endl;
    cout<<"kSwitchRICHMatching  "<<setup_backup.kSwitchRICHMatching <<" : "<<setup.kSwitchRICHMatching <<" : "<<kSwitchRICHMatching <<endl;
    cout<<"fRICHMDCWindow       "<<setup_backup.fRICHMDCWindow      <<" : "<<setup.fRICHMDCWindow      <<" : "<<fRICHMDCWindow      <<endl;
    cout<<"kIsFieldOn           "<<setup_backup.kIsFieldOn          <<" : "<<setup.kIsFieldOn          <<" : "<<kIsFieldOn          <<endl;
    cout<<"kIgnoreRICH          "<<setup_backup.kIgnoreRICH         <<" : "<<setup.kIgnoreRICH         <<" : "<<kIgnoreRICH         <<endl;
    cout<<"kIgnoreInnerMDC      "<<setup_backup.kIgnoreInnerMDC     <<" : "<<setup.kIgnoreInnerMDC     <<" : "<<kIgnoreInnerMDC     <<endl;
    cout<<"kIgnoreOuterMDCB     "<<setup_backup.kIgnoreOuterMDC     <<" : "<<setup.kIgnoreOuterMDC     <<" : "<<kIgnoreOuterMDC     <<endl;
    cout<<"kIgnoreMETA          "<<setup_backup.kIgnoreMETA         <<" : "<<setup.kIgnoreMETA         <<" : "<<kIgnoreMETA         <<endl;
    cout<<"kIgnorePreviousIndex "<<setup_backup.kIgnorePreviousIndex<<" : "<<setup.kIgnorePreviousIndex<<" : "<<kIgnorePreviousIndex<<endl;
    cout<<"pUserSort            "<<setup_backup.pUserSort           <<" : "<<setup.pUserSort           <<" : "<<pUserSort           <<endl;
    cout<<"kUseYMatching        "<<setup_backup.kUseYMatching       <<" : "<<setup.kUseYMatching       <<" : "<<kUseYMatching       <<endl;
    cout<<"kUseYMatchingScaling "<<setup_backup.kUseYMatchingScaling<<" : "<<setup.kUseYMatchingScaling<<" : "<<kUseYMatchingScaling<<endl;
    cout<<"fMetaBoundary        "<<setup_backup.fMetaBoundary       <<" : "<<setup.fMetaBoundary       <<" : "<<fMetaBoundary       <<endl;
    cout<<"kUseBeta             "<<setup_backup.kUseBeta            <<" : "<<setup.kUseBeta            <<" : "<<kUseBeta            <<endl;
    cout<<"fBetaLepCut          "<<setup_backup.fBetaLepCut         <<" : "<<setup.fBetaLepCut         <<" : "<<fBetaLepCut         <<endl;
    cout<<"kUseFakeRejection    "<<setup_backup.kUseFakeRejection   <<" : "<<setup.kUseFakeRejection   <<" : "<<kUseFakeRejection   <<endl;
    cout<<"kUseMETAQA           "<<setup_backup.kUseMETAQA          <<" : "<<setup.kUseMETAQA          <<" : "<<kUseMETAQA          <<endl;
    cout<<"fMetaBoundary        "<<setup_backup.fMetaBoundary       <<" : "<<setup.fMetaBoundary       <<" : "<<fMetaBoundary       <<endl;
    cout<<"fMETAQACut           "<<setup_backup.fMETAQACut          <<" : "<<setup.fMETAQACut          <<" : "<<fMETAQACut          <<endl;
    cout<<"fTOFMAXCut           "<<setup_backup.fTOFMAXCut          <<" : "<<setup.fTOFMAXCut          <<" : "<<fTOFMAXCut          <<endl;
    cout<<"fDeltaSegCut         "<<setup_backup.fDeltaSegCut        <<" : "<<setup.fDeltaSegCut        <<" : "<<fDeltaSegCut        <<endl;
    cout<<"--------------------------------------------------------------------"<<endl;

}

void HParticleTrackSorter::resetFlags(Bool_t flag, Bool_t reject, Bool_t used, Bool_t lepton)
{
    // reset all flags of HParticleCand by looping over the
    // HParticleCand category. The flags kIsRejected, kIsUsed
    // and kIsLepton are special for selection and handles independend
    // from the other flags. The lower flags of objects marked
    // with kIsUsed cand be only reste if used == kRTUE.

    if(!iterParticleCandCat) return;
    HParticleCand* pCand;
    iterParticleCandCat->Reset();
    while ((pCand = dynamic_cast<HParticleCand*>(iterParticleCandCat->Next())) != 0) // begin of ParticleCand iterator
    {
	if(flag){
	    if((pCand->isFlagBit(Particle::kIsUsed) && used)  ||
	       !pCand->isFlagBit(Particle::kIsUsed) ){
		// we want to keep the flags from previous selection
		// clean condition, only reject,used and lepton bit stays!
		for(Int_t i=0;i<Particle::kIsLepton;i++) { pCand->unsetFlagBit(i);  }  
	    }
	}
	if(reject) pCand->unsetFlagBit(Particle::kIsRejected);
	if(used)   pCand->unsetFlagBit(Particle::kIsUsed);
	if(lepton) pCand->unsetFlagBit(Particle::kIsLepton);
    }
}
void HParticleTrackSorter::selection(Bool_t (*function)(HParticleCand* ))
{
    // loops over the HParticleCand category and marks all objects
    // which does not fullfill the selection criteria with the reject flag.
    // Objects marked as kIsUsed are skipped.
    if(!iterParticleCandCat) return;
    HParticleCand* pCand;
    iterParticleCandCat->Reset();
    while ((pCand = dynamic_cast<HParticleCand*>(iterParticleCandCat->Next())) != 0) // begin of ParticleCand iterator
    {
	if( pCand->isFlagBit(Particle::kIsUsed)) { continue;} // do not touch already used objects!
	if(!pCand->select(function)) {pCand->setFlagBit(Particle::kIsRejected); }
    }
}
void HParticleTrackSorter::cleanUp(Bool_t final)
{
    // delete temporary candidate objects from
    // heap and clear the candidate vector. clear
    // all index vectors.

    clearVector(all_candidates);

    if(kIgnorePreviousIndex || final){
	index_RICH    .clear();
	index_InnerMDC.clear();
	index_OuterMDC.clear();
	index_SHOWER  .clear();
	index_TOFHit  .clear();
	index_TOFClst .clear();
	index_RPC     .clear();
    }
}

Int_t  HParticleTrackSorter::fillAndSetFlags()
{
    // clean up the memory and index vectors and fill the
    // candidate vector new. Flag the kIsAcceptedXXX criteria
    // and kIsDoubleHitXXX criteria of the HParticleCand objects of
    // the selection (skipp objects which are marked kIsUsed and
    // kIsRejected) returns the number of candidates.


    if(!pParticleCandCat) return 0; // nothing to do
    fill_Iteration++;

    cleanUp(kFALSE);
    //------------------------------------------------------------------
    // fill HParticleCand Info to temorary objects
    // which will be jused for sorting
    // all flags in HParticleCand will be reset (besides kIsRejected)
    fillInput(all_candidates);


    // list of thing to to
    // for each Double_t hit
    Particle::ESwitch myIndex  [] = {
	Particle::kIsIndexRICHSorter,
	Particle::kIsIndexInnerMDCSorter,
	Particle::kIsIndexOuterMDCSorter,
	Particle::kIsIndexTOFSorter,
	Particle::kIsIndexSHOWERSorter,
        Particle::kIsIndexRPCSorter
    };

    // and for different quality criteria
    Particle::ESwitch myQuality[] = {
        Particle::kIsBestHitRICHSorter,
        Particle::kIsBestHitInnerMDCSorter,
        Particle::kIsBestHitOuterMDCSorter,
        Particle::kIsBestHitMETASorter,
        Particle::kIsBestRKSorter,
	Particle::kIsBestRKRKMETASorter,
	Particle::kIsBestRKRKMETARadiusSorter,
        Particle::kIsBestUserSorter
   };

    kSwitchIndex   = Particle::kIsIndexRICHSorter;
    kSwitchQuality = Particle::kIsBestHitRICHSorter;

    // check all objects for validity with respect to
    // quality criteria and flag them in HParticleCand
    for(UInt_t i = 0; i < (sizeof(myQuality) / sizeof(Int_t)); i ++){
	flagAccepted(all_candidates,myQuality[i]);
    }

    // vector is resorted by index! flag all Double_t used
    // detector hits in HParticleCand
    for(UInt_t i = 0; i < (sizeof(myIndex) / sizeof(Int_t)); i ++){
        flagDouble(all_candidates,myIndex[i]);
    }

    return all_candidates.empty() ? 0 : all_candidates.size();
}

void  HParticleTrackSorter::backupFlags(Bool_t onlyFlags)
{
    // backup the flags of HParticleCand objects. The flags cand be
    // restored to the objects by calling restoreFlags. The flags
    // which have been stored previously will be cleared.

    if(!onlyFlags){
	setup_backup.copyFromStatic();
    }
    old_flags.clear();
    old_flags_emc.clear();

    if(!iterParticleCandCat) return;

    HParticleCand* pCand ;
    iterParticleCandCat->Reset();

    while ((pCand = dynamic_cast<HParticleCand*>(iterParticleCandCat->Next())) != 0 )
    {
         old_flags.push_back(pCand->getFlagField());
    }
    if(pEmcClusterCat){
	HEmcCluster* emc = 0;
	Int_t n = pEmcClusterCat->getEntries();
	for(Int_t i = 0; i < n; i ++){
	    emc = (HEmcCluster*) pEmcClusterCat->getObject(i);
	    if(emc) {
		Int_t fl = emc -> isUsedInParticleCand();
		old_flags_emc.push_back(fl);
	    }
	}
    }

}

Bool_t  HParticleTrackSorter::restoreFlags(Bool_t onlyFlags)
{
    // restore the flags of HParticleCand objects since the last call of
    // backupFlags(). Returns kFALSE if size of vector old_flags does
    // not match the number of objects in HParticleCand category

    if(!onlyFlags){
	setup_backup.copyToStatic();
    }

    if(!iterParticleCandCat) return kFALSE;

    HParticleCand* pCand ;
    iterParticleCandCat->Reset();

    vector<Int_t>::iterator iter;
    iter = old_flags.begin();
    while ((pCand = dynamic_cast<HParticleCand*>(iterParticleCandCat->Next())) != 0 )
    {
	if(iter != old_flags.end()){
	    pCand->setFlagField(*iter);
	    iter ++;
	} else {
            Error("restoreFlags()","Size of backuped flags does not match number of HParticleCand objects!");
	    return kFALSE;
	}
    }

    if(pEmcClusterCat){
	HEmcCluster* emc = 0;
	Int_t n = pEmcClusterCat->getEntries();
	for(Int_t i = 0; i < n; i ++){
	    emc = (HEmcCluster*) pEmcClusterCat->getObject(i);
	    if(emc) {
		emc -> unsetIsUsedInParticleCand();
                if(old_flags_emc[i] == 1 ) emc -> isUsedInParticleCand();
	    }
	}
    }

    return kTRUE;
}
void  HParticleTrackSorter::backupSetup()
{
    // backup the setup of HParticleTrackSorter. The setup cand be
    // restored to the objects by calling restoreSetup.

    setup_backup.copyFromStatic();
}

Bool_t  HParticleTrackSorter::restoreSetup()
{
    // restore the setup of HParticleTrackSorter since the last call of
    // backupSetup().

    setup_backup.copyToStatic();

    return kTRUE;
}

Int_t HParticleTrackSorter::fill(Bool_t (*function)(HParticleCand* ))
{
    // fill the input candidate vector. The kIsRejected flag and
    // lower flags (below kIsLepton) is reseted for all objects
    // not marked as kIsUsed. The filling is done twice once for
    // all objects and second after the application of the user
    // provide test function. Returns the number of candidates after
    // the selection. If kDebug is set and printLevel >= 1 the number
    // of candidates before and after selection will be printed.

    resetFlags(kTRUE,kTRUE,kFALSE,kFALSE);  // need to reset reject/double hit bits (already used objects will stay)!
    Int_t first  = fillAndSetFlags();       // fill all candidates
    selection(function);                    // flag candidates rejected if the they fail the test
    Int_t second = fillAndSetFlags();       // fill only good candidates
    if(kDebug && printLevel >= 1) {
	cout<<"fill(): "
	    <<"s: "                <<selectBest_Iteration
	    <<", f: "              <<fill_Iteration
	    <<", kSwitchParticle: "<<kSwitchParticle
	    <<", nCands first: "   <<first
	    <<", nCands second :"  <<second
	    <<endl;
    }
    return all_candidates.empty() ?  0 : all_candidates.size();
}

Int_t HParticleTrackSorter::selectBest(Particle::ESwitch byQuality, Int_t byParticle)
{
    // perfomes the selection of the best tracks (with respect to the chosen quality
    // criteria), which makes shure, that each detector hit is used only once. If
    // In case Int_t byParticle == kIsLepton the RICH index
    // is taken into account, otherwise not. For leptons the kIsLepton bit is flagged
    // in addition. Individual detector hits can be ignored in the selection procedure with
    // the setIgnoreXXX(kTRUE) functions (by default all hits are taken into account).
    // Returns the number of selected best HParticleCand objects.
    // if kDebug is set and  printLevel >= 1 the number of sected objects and
    // the statistics of the rejection on Double_t hits i the snhgle detectors is printed.
    // if printLevel >=3 the properties of the candidates (accepted and rejected) are
    // printed.

    if(fEvent->getHeader()->getEventSeqNumber() != currentEvent){
	currentEvent         = fEvent->getHeader()->getEventSeqNumber();
        selectBest_Iteration = 0;
    }
    selectBest_Iteration++;

    kSwitchQuality  = byQuality;
    kSwitchParticle = byParticle;
    Int_t ct        = 0;
    Int_t ctdouble[4] = {0};

    if(kSwitchQuality==Particle::kIsBestUserSorter && pUserSort==NULL){
	::Error ("selectBest()","qualtiy is kIsBestUser, but pUserSort==NULL. No sort done.");
    } else {
	sort(all_candidates.begin(),all_candidates.end(),this->cmpQuality);
    }
    for(UInt_t i = 0; i < all_candidates.size(); ++ i)
    {
        candidateSort* cand = all_candidates[i];

        //-------------------------------------------------------
	// skip the candidates which contain a single detector hit
        // which has been used before
        Int_t isDoubleHit = 0;
        /*
	// disabled :: HParticleCand objects are not created for all comb.
	// of innerMDC  and RICH ! This caused problems with close pairs
	// (2 good tracks matched to same RICH ring).
	if( (kSwitchParticle == Particle::kIsLeptonSorter && !kIgnoreRICH ) ) {
            // look to RICH hits only in case of leptons
	    if(cand->ind_RICH     >= 0 &&
	       (find(index_RICH.begin(),index_RICH.end()        ,cand->ind_RICH)     != index_RICH.end()) ){
		isDoubleHit |= 0x1 << 0;
		ctdouble[0] ++;
	    }
	} */
	if( !kIgnoreInnerMDC ) {
	    if(cand->ind_innerMDC >= 0 &&
	       (find(index_InnerMDC.begin(),index_InnerMDC.end(),cand->ind_innerMDC) != index_InnerMDC.end()) ){
		isDoubleHit |= 0x1 << 1;
		ctdouble[1] ++;
	    }
	}
	if( !kIgnoreOuterMDC ) {
	    if(cand->ind_outerMDC >= 0 &&
	       (find(index_OuterMDC.begin(),index_OuterMDC.end(),cand->ind_outerMDC) != index_OuterMDC.end()) ){
		isDoubleHit |= 0x1 << 2;
		ctdouble[2] ++;
	    }
	}
	if( !kIgnoreMETA ) {
	    if( (cand->selectedMeta == Particle::kTofHit1 || cand->selectedMeta == Particle::kTofHit2)  &&
	       cand->ind_TOF      >= 0 &&
	       (find(index_TOFHit.begin(),index_TOFHit.end()    ,cand->ind_TOF)      != index_TOFHit.end()) ){
		isDoubleHit |= 0x1 << 3;
		ctdouble[3] ++;

	    }

	    if(cand->selectedMeta == Particle::kTofClst && cand->ind_TOF      >= 0 &&
	       (find(index_TOFClst.begin(),index_TOFClst.end()   ,cand->ind_TOF)     != index_TOFClst.end()) ){
		isDoubleHit |= 0x1 << 3;
		ctdouble[3] ++;

	    }
            /*
	    // disabled :: HParticleCand objects are not created for all comb.
	    // of outerMDC  and SHOWER!
	    if(cand->selectedMeta == Particle::kShowerHit && cand->ind_SHOWER   >= 0 &&
	       (find(index_SHOWER.begin(),index_SHOWER.end()    ,cand->ind_SHOWER)   != index_SHOWER.end()) ){
		isDoubleHit |= 0x1 << 3;
		ctdouble[3] ++;
	    }
            */
	    if(cand->selectedMeta == Particle::kRpcClst   && cand->ind_RPC      >= 0 &&
	       (find(index_RPC.begin(),index_RPC.end()          ,cand->ind_RPC)      != index_RPC.end()) ){
		isDoubleHit |= 0x1 << 3;
		ctdouble[3] ++;
	    }
	}
	if(isDoubleHit){
            if(kDebug && printLevel >= 1){
		TString out="";
		for(Int_t j = 3; j >= 0; j --){
		    out += ( (isDoubleHit>>j) & 0x1 );
		}
		cout<<"--selectBest(): skip Double_t hit with pattern : "<<out.Data()<<endl;
		printCand(cand,i);
	    }
	    continue;
	}
        //-------------------------------------------------------

        //-------------------------------------------------------
	// mark the candidate as best and add the single detector hit indices
	// to the lookup

	/*
	// disabled :: HParticleCand objects are not created for all comb.
	// of innerMDC and RICH ! This caused problems with close pairs
	// (2 good tracks matched to same RICH ring).
	if(cand->ind_RICH     >= 0 ){
	    index_RICH.push_back(cand->ind_RICH);
	}
	*/

	if(cand->ind_innerMDC >= 0 ){
	    index_InnerMDC.push_back(cand->ind_innerMDC);
	}
	if(cand->ind_outerMDC >= 0 ){
	    index_OuterMDC.push_back(cand->ind_outerMDC);
	}
	if( (cand->selectedMeta == Particle::kTofHit1 || cand->selectedMeta == Particle::kTofHit2) &&
	   cand->ind_TOF      >= 0 ){
	    index_TOFHit.push_back(cand->ind_TOF);
	}
	if(cand->selectedMeta == Particle::kTofClst &&
	   cand->ind_TOF      >= 0 ){
	    index_TOFClst.push_back(cand->ind_TOF);
	}
        /*
        // disabled :: HParticleCand objects are not created for all comb.
	// of outerMDC and SHOWER !
	if(cand->selectedMeta == Particle::kShowerHit &&
	   cand->ind_SHOWER   >= 0 ){
	    index_SHOWER.push_back(cand->ind_SHOWER);
	}
	*/
	if(cand->selectedMeta == Particle::kRpcClst &&
	   cand->ind_RPC      >= 0 ){
	    index_RPC.push_back(cand->ind_RPC);
	}

	HParticleCand* pcand = dynamic_cast<HParticleCand*>(pParticleCandCat->getObject(cand->ind_Cand));
	if(!pcand) {
	    Error("HParticleTrackSorter::selectBest()","NULL pointer retrieved for HParticleCand at index %i",cand->ind_Cand);
	    continue;
	}
	switch (byQuality)
	{
	case Particle::kIsBestHitRICHSorter:
	    pcand -> setFlagBit(Particle::kIsBestHitRICH);
	    break;
	case Particle::kIsBestHitInnerMDCSorter:
	    pcand -> setFlagBit(Particle::kIsBestHitInnerMDC);
	    break;
	case Particle::kIsBestHitOuterMDCSorter:
	    pcand -> setFlagBit(Particle::kIsBestHitOuterMDC);
	    break;
	case Particle::kIsBestHitMETASorter:
	    pcand -> setFlagBit(Particle::kIsBestHitMETA);
	    break;
	case Particle::kIsBestRKSorter:
	    pcand -> setFlagBit(Particle::kIsBestRK);
	    break;
	case Particle::kIsBestRKRKMETASorter:
	    pcand -> setFlagBit(Particle::kIsBestRKRKMETA);
	    break;
	case Particle::kIsBestRKRKMETARadiusSorter:
	    pcand -> setFlagBit(Particle::kIsBestRKRKMETARadius);
	    break;
	case Particle::kIsBestUserSorter:
	    pcand -> setFlagBit(Particle::kIsBestUser);
	    break;
	default :

	    Error("HParticleTrackSorter::selectBest()","Unknown switch statement : %i!",byQuality);
	    exit(1);
	    break;

	}
        ++ ct;
	pcand ->setFlagBit(Particle::kIsUsed);
	if(byParticle == Particle::kIsLeptonSorter){
	    pcand ->setFlagBit(Particle::kIsLepton);
	}
	if(kDebug && printLevel >= 1){
	    cout<<"++selectBest(): select hit : "<<endl;
	    printCand(cand,i);
	}
	//-------------------------------------------------------

    }
    if(kDebug && printLevel >= 1)
    {
	cout<<"selectBest(): "
	    <<"s: "                  <<selectBest_Iteration
	    <<", f: "                <<fill_Iteration
	    <<", kSwitchParticle: "  <<kSwitchParticle
	    <<", nBest cand: "       <<ct
	    <<", Double_t RICH: "    <<ctdouble[0]
            <<", InnerMDC: "         <<ctdouble[1]
            <<", Double_t OuterMDC: "<<ctdouble[2]
            <<", Double_t META: "    <<ctdouble[3]
	    <<endl;
    }
    fill_Iteration = 0;
    flagEmcClusters();
  return ct;
}

void HParticleTrackSorter::setRICHMatching(Particle::ERichMatch match, Float_t window)
{
    // switch between the matching cut RICH-MDC inside the selectLeptons() function. The method
    // is selected by the options defined in Particle::ERichMatch.
	kSwitchRICHMatching = match;
	fRICHMDCWindow      = window;
	setup.kSwitchRICHMatching = match;
	setup.fRICHMDCWindow      = window;
}


Bool_t HParticleTrackSorter::selectLeptons(HParticleCand* pcand)
{
    // build in selection function for lepton candidates.
    // Requires besides an RICH hit, RK + META and fitted
    // inner+outer segment. The RICH-MDC matching cut depends on the
    // settings of the matching method and window (set by
    // HParticleTrackSorter::setRICHMatching . default is
    // kUseRKRICHWindow, window = 4. degree)


    if(kUseFakeRejection&&pcand->isFakeRejected()) return kFALSE;
    Bool_t test = kFALSE;
    if( pcand->isFlagAND(5,
			 Particle::kIsAcceptedHitRICH,
			 Particle::kIsAcceptedHitInnerMDC,
			 Particle::kIsAcceptedHitOuterMDC,
			 Particle::kIsAcceptedHitMETA,
			 Particle::kIsAcceptedRK
			)
       &&
       pcand->getInnerSegmentChi2() > 0
       &&
       pcand->getChi2()             < 1000      // RK
      ) test = kTRUE;

    if(!test) return kFALSE;

    if(kUseMETAQA && (pcand->getMetaMatchQuality() > fMETAQACut || pcand->getMetaMatchQuality() < 0)) return kFALSE;
    if(kUseYMatching && !HParticleTool::isGoodMetaCell(pcand,fMetaBoundary,kUseYMatchingScaling)) return kFALSE;
    if(kUseBeta && pcand->getBeta()<fBetaLepCut) return kFALSE;
    if(!kIsFieldOn && !isDeltaSegOk(pcand)) return kFALSE;

    if(test)
    {   // rich match only!
	if( kSwitchRICHMatching == Particle::kUseRICHIndexSorter)
	{
	    return kTRUE;
	}
	else if ( kSwitchRICHMatching == Particle::kUseRKRICHCorrelationSorter	&&  pcand -> getRingCorr() )
	{
	    return kTRUE;
	}
	else if ( kSwitchRICHMatching == Particle::kUseRKRICHWindowSorter
		 &&
		 fabs(pcand->getDeltaTheta()) < fRICHMDCWindow
		 &&
		 fabs(pcand->getDeltaPhi())   < fRICHMDCWindow
		)
	{
	    return kTRUE;
	}
	else return kFALSE;

    }
    else return kFALSE;
}
Bool_t HParticleTrackSorter::selectLeptonsNoOuterFit(HParticleCand* pcand)
{
    // build in selection function for lepton candidates.
    // Requires besides an RICH hit, RK + META and fitted
    // inner. Fit for outer segment not reuired. The
    // RICH-MDC matching cut depends on the
    // settings of the matching method and window (set by
    // HParticleTrackSorter::setRICHMatching . default is
    // kUseRKRICHWindow, window = 4. degree)

    if(kUseFakeRejection&&pcand->isFakeRejected()) return kFALSE;
    Bool_t test = kFALSE;
    if( pcand->isFlagAND(4,
			 Particle::kIsAcceptedHitRICH,
			 Particle::kIsAcceptedHitInnerMDC,
			 Particle::kIsAcceptedHitMETA,
			 Particle::kIsAcceptedRK
			)
       &&
       pcand->getInnerSegmentChi2() > 0
       &&
       pcand->getOuterSegInd()      > 0
       &&
       pcand->getChi2()             < 1000      // RK
      ) test = kTRUE;

    if(!test) return kFALSE;

    if(kUseMETAQA && (pcand->getMetaMatchQuality() > fMETAQACut || pcand->getMetaMatchQuality() < 0)) return kFALSE;
    if(kUseYMatching && !HParticleTool::isGoodMetaCell(pcand,fMetaBoundary,kUseYMatchingScaling)) return kFALSE;
    if(kUseBeta && pcand->getBeta()<fBetaLepCut) return kFALSE;
    if(!kIsFieldOn && !isDeltaSegOk(pcand)) return kFALSE;

    if(test)
    {   // rich match only!
	if( kSwitchRICHMatching == Particle::kUseRICHIndexSorter)
	{
	    return kTRUE;
	}
	else if ( kSwitchRICHMatching == Particle::kUseRKRICHCorrelationSorter	&&  pcand -> getRingCorr() )
	{
	    return kTRUE;
	}
	else if ( kSwitchRICHMatching == Particle::kUseRKRICHWindowSorter
		 &&
		 fabs(pcand->getDeltaTheta()) < fRICHMDCWindow
		 &&
		 fabs(pcand->getDeltaPhi())   < fRICHMDCWindow
		)
	{
	    return kTRUE;
	}
	else return kFALSE;
    }
    else return kFALSE;
}
Bool_t HParticleTrackSorter::selectHadrons(HParticleCand* pcand )
{
    // build in selection function for hadron candidates.
    // Requires besides RK + META and fitted inner+outer MDC.

    if(kUseFakeRejection&&pcand->isFakeRejected()) return kFALSE;
    Bool_t test=kFALSE;
    if( pcand->isFlagAND(4,
			 Particle::kIsAcceptedHitInnerMDC,
			 Particle::kIsAcceptedHitOuterMDC,
			 Particle::kIsAcceptedHitMETA,
			 Particle::kIsAcceptedRK
			)
       &&
       pcand->getInnerSegmentChi2() > 0
       &&
       pcand->getChi2()             < 1000
      ) test=kTRUE;
    if(!test) return kFALSE;

    if(kUseMETAQA && (pcand->getMetaMatchQuality() > fMETAQACut || pcand->getMetaMatchQuality() < 0)) return kFALSE;
    if(kUseYMatching && !HParticleTool::isGoodMetaCell(pcand,fMetaBoundary,kUseYMatchingScaling)) return kFALSE;
    if(kUseBeta && (pcand->getBeta()<0 || pcand->getTof()>fTOFMAXCut)) return kFALSE;
    if(!kIsFieldOn && !isDeltaSegOk(pcand)) return kFALSE;

    return test;
}
Bool_t HParticleTrackSorter::selectHadronsNoOuterFit(HParticleCand* pcand )
{
    // build in selection function for hadron candidates.
    // Requires besides RK + META and fitted inner MDC. Fit
    // for outer MDC not required.

    if(kUseFakeRejection&&pcand->isFakeRejected()) return kFALSE;
    Bool_t test=kFALSE;
    if( pcand->isFlagAND(3,
			 Particle::kIsAcceptedHitInnerMDC,
			 Particle::kIsAcceptedHitMETA,
			 Particle::kIsAcceptedRK
			)
       &&
       pcand->getInnerSegmentChi2() > 0
       &&
       pcand->getOuterSegInd()      > 0
       &&
       pcand->getChi2()             < 1000
      ) test=kTRUE;

    if(!test) return kFALSE;

    if(kUseMETAQA && (pcand->getMetaMatchQuality() > fMETAQACut || pcand->getMetaMatchQuality() < 0)) return kFALSE;
    if(kUseYMatching && !HParticleTool::isGoodMetaCell(pcand,fMetaBoundary,kUseYMatchingScaling)) return kFALSE;
    if(kUseBeta && (pcand->getBeta()<0 || pcand->getTof()>fTOFMAXCut)) return kFALSE;
    if(!kIsFieldOn && !isDeltaSegOk(pcand)) return kFALSE;

    return test;
}



void sorter_setup::copyFromStatic()
{

	kSwitchRICHMatching  = HParticleTrackSorter::kSwitchRICHMatching;
	fRICHMDCWindow       = HParticleTrackSorter::fRICHMDCWindow;
	kIsFieldOn           = HParticleTrackSorter::kIsFieldOn;
	kIgnoreRICH          = HParticleTrackSorter::kIgnoreRICH;
	kIgnoreInnerMDC      = HParticleTrackSorter::kIgnoreInnerMDC;
	kIgnoreOuterMDC      = HParticleTrackSorter::kIgnoreOuterMDC;
	kIgnoreMETA          = HParticleTrackSorter::kIgnoreMETA;
	kIgnorePreviousIndex = HParticleTrackSorter::kIgnorePreviousIndex;
	pUserSort            = HParticleTrackSorter::pUserSort;
	kUseYMatching        = HParticleTrackSorter::kUseYMatching;
	kUseYMatchingScaling = HParticleTrackSorter::kUseYMatchingScaling;
	fMetaBoundary        = HParticleTrackSorter::fMetaBoundary;
	kUseBeta             = HParticleTrackSorter::kUseBeta;
	fBetaLepCut          = HParticleTrackSorter::fBetaLepCut;
	kUseFakeRejection    = HParticleTrackSorter::kUseFakeRejection;
	kUseMETAQA           = HParticleTrackSorter::kUseMETAQA;
	fMETAQACut           = HParticleTrackSorter::fMETAQACut;
	fTOFMAXCut           = HParticleTrackSorter::fTOFMAXCut;
	fDeltaSegCut         = HParticleTrackSorter::fDeltaSegCut;
};

void sorter_setup::copyToStatic()
{

	HParticleTrackSorter::kSwitchRICHMatching  = kSwitchRICHMatching ;
	HParticleTrackSorter::fRICHMDCWindow       = fRICHMDCWindow;
	HParticleTrackSorter::kIsFieldOn           = kIsFieldOn;
	HParticleTrackSorter::kIgnoreRICH          = kIgnoreRICH;
	HParticleTrackSorter::kIgnoreInnerMDC      = kIgnoreInnerMDC;
	HParticleTrackSorter::kIgnoreOuterMDC      = kIgnoreOuterMDC;
	HParticleTrackSorter::kIgnoreMETA          = kIgnoreMETA;
	HParticleTrackSorter::kIgnorePreviousIndex = kIgnorePreviousIndex;
	HParticleTrackSorter::pUserSort            = pUserSort;
	HParticleTrackSorter::kUseYMatching        = kUseYMatching;
	HParticleTrackSorter::kUseYMatchingScaling = kUseYMatchingScaling;
	HParticleTrackSorter::fMetaBoundary        = fMetaBoundary;
	HParticleTrackSorter::kUseBeta             = kUseBeta;
	HParticleTrackSorter::fBetaLepCut          = fBetaLepCut;
	HParticleTrackSorter::kUseFakeRejection    = kUseFakeRejection;
	HParticleTrackSorter::kUseMETAQA           = kUseMETAQA;
	HParticleTrackSorter::fMETAQACut           = fMETAQACut;
	HParticleTrackSorter::fTOFMAXCut           = fTOFMAXCut;
	HParticleTrackSorter::fDeltaSegCut         = fDeltaSegCut;
};