Preparation_ZDC.C

Go to the documentation of this file.

 
using namespace std;
 
void Preparation_ZDC(char* fname = "../raw/bmn_run0871.root")
{
    gROOT->LoadMacro("$VMCWORKDIR/macro/run/bmnloadlibs.C");
    bmnloadlibs();
    /////////////////////////////////////////////////////////////////////////////////////
    int RUN, n_notrig = 0, n_phys = 0, n_led = 0, n_bad = 0;
    UInt_t TRIGWORD = 0;
    sscanf(&fname[strlen(fname) - 8], "%d", &RUN);
 
    const char* mapping;
    mapping = "ZDC_map_period_5.txt";
 
    BmnZDCRaw2Digit ZDC(mapping, fname, "zdc_muon_calibration.txt");
    ZDC.print();
 
    cout << "Process RUN:  " << RUN << endl;
 
    /////////////////////////////////////////////////////////////////////////////////////
    TFile* _f_in = new TFile(fname, "READ");
    TTree* _t_in = (TTree*)_f_in->Get("BMN_RAW");
    TClonesArray* zdc_raw = new TClonesArray("BmnADCDigit");
    _t_in->SetBranchAddress("bmn_zdc", &zdc_raw);
    _t_in->SetBranchAddress("bmn_trigword", &TRIGWORD);
    /////////////////////////////////////////////////////////////////////////////////////
 
    for (int ev = 0; ev < _t_in->GetEntries(); ev++) {
 
        if ((ev % 1000) == 0)
            printf("Preparation ZDC, event %d, phys %d, led %d\n", ev, n_phys, n_led);
 
        TRIGWORD = 0;
        zdc_raw->Clear();
 
        _t_in->GetEntry(ev);
 
        if (TRIGWORD == 0) {
            n_notrig++;
        } else if ((TRIGWORD & 0x8) == 0) {
            n_phys++;
            ZDC.fillAmplitudes(zdc_raw);
            ZDC.fillSampleProfilesAll(zdc_raw, 0., 0., 48.);
            // ZDC.fillSampleProfiles(zdc_raw, 0., 0., 48., 3);
        } else if ((TRIGWORD & 0x8) == 0x8) {
            n_led++;
            // ZDC.fillAmplitudes(zdc_raw);
            // ZDC.fillSampleProfilesAll(zdc_raw, 0., 0., 48.);
 
            // ZDC.fillSampleProfiles(zdc_raw, 0., 0., 48., 3);
            // if ((ev % 1000) == 0 && ZDC.getLogChan(0) >= 0) printf("     Log channel 0 (Id = 0x%0x, Ch = %d)
            // amplitude is %f\n", ZDC.getLogId(0), ZDC.getLogChan(0), ZDC.getLogAmp(0)); if ((ev % 1000) == 0 &&
            // ZDC.getLogChan(1) >= 0) printf("     Log channel 1 (Id = 0x%0x, Ch = %d) amplitude is %f\n",
            // ZDC.getLogId(1), ZDC.getLogChan(1), ZDC.getLogAmp(1));
        } else {
            n_bad++;
        }
    }
    /////////////////////////////////////////////////////////////////////////////////////
    ZDC.drawprof();
    ZDC.drawtest();
    /////////////////////////////////////////////////////////////////////////////////////
 
    _f_in->Close();
}
 
void select_hist()
{
    // example of macro called when a mouse event occurs in a pad.
    //  Example:
    //  Root > TFile f("hsimple.root");
    //  Root > hpxpy.Draw();
    //  Root > c1.AddExec("exsel","select_hist()");
    //  When moving the mouse in the canvas, a second canvas shows the
    //  projection along X of the bin corresponding to the Y position
    //  of the mouse. The resulting histogram is fitted with a gaussian.
    //  A "dynamic" line shows the current bin position in Y.
    //  This more elaborated example can be used as a starting point
    //  to develop more powerful interactive applications exploiting CINT
    //  as a development engine.
    // Author: Rene Brun
 
    int event = gPad->GetEvent();
    if (event != 11)
        return;
 
    TObject* select = gPad->GetSelected();
    if (!select)
        return;
    if (!(select->InheritsFrom("TH1")))
        return;
 
    // create or set the new canvas c2
    TCanvas* csingle = 0;
    if ((csingle = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("csingle")))
        ;
    else
        csingle = new TCanvas("csingle");
    csingle->cd();
 
    // draw slice corresponding to mouse position
    if (!(select->InheritsFrom("TH1F"))) {
        TH1F* h = (TH1F*)select;
        h->Draw();
    } else if (!(select->InheritsFrom("TH2F"))) {
        TH2F* h = (TH2F*)select;
        h->Draw();
    } else if (!(select->InheritsFrom("TProfile"))) {
        TProfile* h = (TProfile*)select;
        h->Draw();
    } else if (!(select->InheritsFrom("TH1D"))) {
        TH1D* h = (TH1D*)select;
        h->Draw();
    } else if (!(select->InheritsFrom("TH2D"))) {
        TH2D* h = (TH2D*)select;
        h->Draw();
    } else if (!(select->InheritsFrom("TProfile2D"))) {
        TProfile2D* h = (TProfile2D*)select;
        h->Draw();
    } else
        return;
    csingle->Update();
}