[Halld-offline] pythia event filter

[Aram Teymurazyan]

Hi Will,


here is brute force example of selecting p eta events:

==============================================================
  vector<const DMCThrown*> throwns;
  eventLoop->Get(throwns);

int N_FINAL_STATES =0;

for(unsigned int i=0; i<throwns.size(); i++)
     {
       if(throwns[i]->mech==1)
         {
           N_FINAL_STATES++;
           if(throwns[i]->pdgtype==2212)flag_proton  = true;
           if(throwns[i]->pdgtype==221)flag_eta = true;
         }
     }

   if(N_FINAL_STATES!=2 && flag_eta!=true && flag_proton!=true) return 
NOERROR;
==============================================================

this should be part of:

jerror_t MyProcessor::evnt(JEventLoop *eventLoop, int eventnumber);

-Aram

On Mon, 22 Mar 2010, Will Levine wrote:

> Hi,
> I'm looking for a way to filter out all events of a specific final
> states from bggen/pythia output. Curtis said there had been some
> discussion of this, so I was hoping someone could help me out.
>
> Thanks,
> Will
> _______________________________________________
> Halld-offline mailing list
> Halld-offline at jlab.org
> https://mailman.jlab.org/mailman/listinfo/halld-offline
>