class HParticleGeantEvent: public TObject

_HADES_CLASS_DESCRIPTION



 HParticleGeantEvent

 Simulation object keeping HGeant pointers and decays to ease
 the geant use.



 USAGE:


  HParticleGeantEvent geantevent; // create the object out side event loop



  for (Int_t i=0; i < entries; i++) {
    Int_t nbytes =  loop.nextEvent(i);             // get next event. categories will be cleared before
    if(nbytes <= 0) { cout<<nbytes<<endl; break; } // last event reached


    geantevent.nextEvent(); // fill the geant event structure

    cout<<"event "<<i<<"------------------------------------------"<<endl;
    cout<<"impact : "<<geantevent.getImpactParam()<<", beam E: "<< geantevent.getBeamEnergy()<< ", event plane: "<<geantevent.getEventPlane()<<endl;
    cout<<"mult charged prim 4pi "<<geantevent.getMultiplicity(-1,0,0,0)             // id,generation,charge,type,detbits
        <<" acc "<< geantevent.getMultiplicity(-1,0,0,1)
        <<" reco true iMDC+oMDC+META "<<geantevent.getMultiplicity(-1,0,0,2,kIsInInnerMDC|kIsInOuterMDC|kIsInMETA)
        <<" reco true all "<<geantevent.getMultiplicity(-1,0,0,2,0)
        <<endl;

    //----------------------------------------------------
    // example 1:
    for(UInt_t j = 0; j < geantevent.getDecays().size(); j++ ){
        HParticleGeantDecay* decay = geantevent.getDecays()[j];
        if(decay->getDistFromVertex()<100 ) decay->print(); // print all decays max 10cm arround 0,0,0
    }

    //----------------------------------------------------
    // example 2 :
    vector<HParticleGeantDecay*> decays;      // a decay keeps pointer to mother and daughters and the mother decay
    geantevent.isDecay(decays,52,-1);         // all omgegas from PLUTO (external : generation = -1) filter params :  motherid, generation, med, dist
    for(UInt_t j = 0; j < decays.size(); j++) {
       decays[j] ->print();
    }
    //----------------------------------------------------
    // example 3 :
    geantevent.isDecay(decays,18,0);          // all primary lambda decays (urqmd+PLUTO(is not decayed outside GEANT))
    vector<Int_t > pids;
    pids.push_back(14);  // proton
    pids.push_back(9);   // pion
    for(UInt_t j = 0; j < decays.size(); j++) {
       HParticleGeantDecay* lambda = decays[j];
       if(!lambda->isExternalSource()) continue;  // keep only PLUTO
       if(lambda->isDaughter(pids)) // both particles found

       vector<HGeantKine*>& daughters = lambda->getDaughters();
       for(UInt_t l = 0; l < daughters.size(); j++){
           HParticleGeant* p = geantevent.getParticle(daughters[l]);
           if(p->isInAcceptance()){
             Int_t nReco p->getNRecoCand(); // number of HParticleCandSim keeping this geant track
             for(Int_t k=0; k < nReco; k++){
                 if(!p->getRecoCand(k)->isGhostTrack() && p->getRecoCand(k)->isIndDetectors(kIsInInnerMDC||kIsInOuterMDC|kIsInMETA)) { // true reco

                 }
             }
           }
        }

    }

    //----------------------------------------------------
    // example 4 :
    vector<HParticleGeant*> particles;
    geantevent.isParticle(particles,2,-2,0); // all primary positrons  possible filter params : id,motherid,generation,med,dist
    for(UInt_t j = 0; j < particles.size(); j++) {
        particles[j]->print();
        HParticleGeantDecay* d = particles[j]->getMotherDecay(); // will even work for PLUTO decayed particles
        if(d) d->print();
    }
    //----------------------------------------------------
    // example 5 :
    HParticleCandSim* cand1=0;
    for(Int_t j = 0; j < particleCat->getEntries(); j ++){
        cand1 = HCategoryManager::getObject(cand1,particleCat,j);
        if(!cand1->isFlagBit(Particle::kIsLepton)) continue;

        Int_t tr = cand1->getGeantTrack();
        HParticleGeant* p      = geantevent.getParticle(tr);
        HParticleGeantDecay* d = p->getMotherDecay();    // mother decay
        Int_t nDaughter        = d ->getNDaughters();
        HGeantKine* sister = 0;
        for(Int_t k=0; k < nDaughter; k++){              // find sister track
            sister = d->getDaughter(k);
            if(sister->getTrack() != p->getParticle()->getTrack()) break;
        }
        HParticleGeant* pSister = geantevent.getParticle(sister);

        Int_t nReco p->getNRecoCand(); // number of HParticleCandSim keeping this geant track
        for(Int_t k=0; k < nReco; k++){
           p->getRecoCand(k)->print(); // all candidates sharing the same geant track
        }
        nReco pSister->getNRecoCand(); // number of HParticleCandSim keeping this geant track
        for(Int_t k=0; k < nReco; k++){
           p->getRecoCand(k)->print(); // all candidates sharing the same geant track
        }
    }
  } // end eventloop