newmen/versatile-diamond

#include "des_methyl_from_dimer.h"

const char DesMethylFromDimer::__name[] = "desorption methyl from dimer";

double DesMethylFromDimer::RATE()
{
    static double value = getRate("DES_METHYL_FROM_DIMER", Env::cH());
    return value;
}

void DesMethylFromDimer::find(MethylOnDimerCMiu *target)
{
    create<DesMethylFromDimer>(target);
}

void DesMethylFromDimer::doIt()
{
    assert(target()->type() == MethylOnDimerCMiu::ID);

    Atom *atoms[2] = { target()->atom(1), target()->atom(0) };
    Atom *a = atoms[0], *b = atoms[1];

    assert(a->is(23));
    assert(b->is(25));

    Handbook::amorph().erase(b);

    a->unbondFrom(b);

    a->changeType(21);

    Handbook::scavenger().markAtom(b);

    Finder::findAll(atoms, 2);
}