ensdfdemo.C

Go to the documentation of this file.

using namespace Gw;

void oneLS(const Char_t *filename, const Char_t *nucleus, Int_t nbcas = 100)
{
        // load the AGS file
        BaseGEM random;
        if ( random.Import(filename,nucleus) == 0 ) {
                cout << " Level Scheme not correctly loaded " << endl; return;
        }

        // some spectra to be filled
        TH1F *K = new TH1F("K","K",100,0,100);
        TH1F *h = new TH1F("ProjTot","ProjTot",8192,0,4096);

        // start the Monte-Carlo
        Cascade cas; GammaLink *gam;
        for (Int_t i = 0; i < nbcas; i++ ) {

                random.DoCascade(cas);

                K->Fill(cas.GetSize()) ;
                for (Int_t j = 0; j < cas.GetSize(); j++ ) {
                        gam = (GammaLink *)cas.At(j); h->Fill(gam->GetEnergy().Get());
                } // j
        } //i
        h->Draw();
}

void testfeeding(const Char_t *filename, const Char_t *nucleus, Int_t nbcas = 100)
{
        // load the level scheme from an external file
        GEM random;
        if ( random.Import(filename,nucleus) == 0 ) {
                cout << " Level Scheme not correctly loaded " << endl; return;
        }
        // Adds a feeding to the GEM process
        random.SetFeeding("TestFeeding"); random.SetExI(80000,60); // entry point 50 MeV, spin 60
        
        // some spectra to be filled
        TH1F *K = new TH1F("K","K",100,0,100);
        TH1F *h = new TH1F("ProjTot","ProjTot",8192,0,4096);

        // start the Monte-Carlo
        Cascade cas; GammaLink *gam;
        for (Int_t i = 0; i < nbcas; i++ ) {

                random.DoCascade(cas);

                K->Fill(cas.GetSize()) ;
                for (Int_t j = 0; j < cas.GetSize(); j++ ) {
                        gam = (GammaLink *)cas.At(j); h->Fill(gam->GetEnergy().Get());
                } // j
        } //i
        h->Draw();
}