BmnPreparationECAL.C

Go to the documentation of this file.

 
void BmnPreparationECAL(char* fname = "bmn_run4926_raw.root")
{
    gROOT->LoadMacro("$VMCWORKDIR/macro/run/bmnloadlibs.C");
    bmnloadlibs();
    /////////////////////////////////////////////////////////////////////////////////////
    int RUN;
    sscanf(&fname[strlen(fname) - 13], "%d", &RUN);
 
    char mapping[128] = {"ECAL_map_period_7.txt"};
 
    BmnECALRaw2Digit* ECAL = new BmnECALRaw2Digit(mapping, fname);
    ECAL->print();
 
    cout << "Process RUN:  " << RUN << endl;
 
    /////////////////////////////////////////////////////////////////////////////////////
    TFile* _f_in = new TFile(fname, "READ");
    TTree* _t_in = (TTree*)_f_in->Get("BMN_RAW");
    TClonesArray* ecal_raw = new TClonesArray("BmnADCDigit");
    _t_in->SetBranchAddress("ADC", &ecal_raw);
    /////////////////////////////////////////////////////////////////////////////////////
 
    for (int ev = 0; ev < _t_in->GetEntries(); ev++) {
 
        if ((ev % 5000) == 0)
            printf("Preparation ECAL, event %d\n", ev);
 
        ecal_raw->Clear();
 
        _t_in->GetEntry(ev);
 
        ECAL->fillAmplitudes(ecal_raw);
        ECAL->fillSampleProfilesAll(ecal_raw, 0., 0., 48.);
    }
    /////////////////////////////////////////////////////////////////////////////////////
    ECAL->drawecal(1);     // draw ECAL modules with module number and calibration coefficient
    ECAL->drawprof(229);   // average waveform for module 229
    /////////////////////////////////////////////////////////////////////////////////////
 
    _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();
}