flyopunta/src/src_dev62/devrich.cpp

/*
    <one line to give the program's name and a brief idea of what it does.>
    Copyright (C) <year>  <name of author>

    This program is free software; you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation; either version 2 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License along
    with this program; if not, write to the Free Software Foundation, Inc.,
    51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.

*/
#include "parm.h"
#include "part_db.h"
#include "devrich.h"

using namespace std;

DevRich::DevRich()
        :DevCrn()
{
    devtype=TypDevRich;
    devclass="Ring Imaging Cherenkov";
Gout<<"\n\n < "<<nome <<" > Dev Id "<<idev<<" Device type <"<<devclass<<"> typ "<<devtype; 

}

DevRich::~DevRich()
{
}

int DevRich::SimulaDev()
{
    if(mhit<1)return 0;
    double Thres_pi=15.;
    double Thres_mu=12.;
    qvet Eimp;
    int loop=0;
    const char *Part_nome;
    while (loop<mhit)
    {
        Eimp=M_Hits[loop].Pdev;
        Part_nome= M_Hits[loop].pnome;

        if ( Part_nome==Pip|| Part_nome==Pim)
        {
            if (M_Hits[loop].Pdev.mom <Thres_pi)
                M_Hits[loop].idead=-99999;
            /*uccido la particella*/
        } if (M_Hits[loop].Pdev.mom <Thres_pi)
            M_Hits[loop].idead=-99999;

        else  if ( Part_nome==Mum|| Part_nome==Mup)
        {
            if (M_Hits[loop].Pdev.mom <Thres_mu)
                M_Hits[loop].idead=-99999;
        }






//Gout<<"\n Nome "<<Part_nome<<" e "<<Eimp;
//Gout<<"\n Nome "<<Part_nome<<" e "<<Eimp;
        loop++;
    }
//==============  La parte che segue è fissa!
// elimino le paricelle morte, compattizzo i dati e  ridefinisco nhit e m

    //Gout<<"\n Prima  del coppattamento..."<<nome<<" "<< nhit<<" " <<mhit<< "  "<<M_Hits[nhit].pnome;

    loop =0;
    int newhit=mhit;
    while (loop<newhit)
    {
        if (M_Hits[loop].idead!=-99999)loop++;
        else if (newhit==1)newhit=0;
        else {
            newhit--;
            for (int i=loop;i<newhit;i++)
            {
                M_Hits[i]=M_Hits[i+1];
            }
        }
    }
    mhit=newhit;
    nhit=mhit-1;



    return 0;
}