BmnPreparationTOF700.C

Go to the documentation of this file.

// file: full path to raw-file
// nEvents: if 0 then decode all events
 
void BmnPreparationTOF700(TString file = "mpd_run_trigCode_8760.data", Long_t nEvents = 0)
{
 
    gSystem->ExpandPathName(file);
    gSystem->ExpandPathName(outfile);
 
    Int_t iVerbose = 1;   ///<- Verbosity level: 0 - Progress Bar; 1 - short info on passed events
    UInt_t period = 7;
 
    TStopwatch timer;
    timer.Start();
 
    BmnRawDataDecoder* decoder = new BmnRawDataDecoder(file, outfile, nEvents, period);
    // use kSRCSETUP for Short-Range Correlation program and kBMNSETUP otherwise
    BmnSetup stp = (decoder->GetRunId() >= 2041 && decoder->GetRunId() <= 3588) ? kSRCSETUP : kBMNSETUP;
    decoder->SetBmnSetup(stp);
    decoder->SetVerbose(iVerbose);
    decoder->SetAdcDecoMode(period < 6 ? kBMNADCSM : kBMNADCMK);
 
    TString PeriodSetupExt = Form("%d%s.txt", period, ((stp == kBMNSETUP) ? "" : "_SRC"));
    decoder->SetTrigPlaceMapping(TString("Trig_PlaceMap_Run") + PeriodSetupExt);
    decoder->SetTrigChannelMapping(TString("Trig_map_Run") + PeriodSetupExt);
    decoder->SetTOF700ReferenceRun(-1);
    decoder->SetTof700Geom(TString("TOF700_geometry_run") + PeriodSetupExt);
    if (decoder->GetRunId() >= 4278 && decoder->GetPeriodId() == 7)
        decoder->SetTof700Mapping(TString("TOF700_map_period_") + Form("%d_from_run_4278.txt", period));
    else
        decoder->SetTof700Mapping(TString("TOF700_map_period_") + Form("%d.txt", period));
 
    decoder->PreparationTOF700Init();   // Decode data into detector-digits using current mappings.
    BmnTof2Raw2DigitNew* tof700m = decoder->GetTof700Mapper();
 
    if (period <= 7) {
        if (stp == kSRCSETUP) {
            tof700m->SetT0shift(-15000.);
        } else {
            tof700m->SetT0shift(+147000.);
        }
    } else if (period == 8) {
        tof700m->SetT0shift(+10000.);
    } else {
        printf
        {"wrong perion number - %d !\n", period);
        }
 
        tof700m->Book();
        decoder->PreparationTOF700();   // Analyze TOF700 data for slewing time limits.
        // draw quality control histogram - comment if not neccessary
        tof700m->WritePreparationHists();
        tof700m->drawprep();
        delete decoder;
    }
 
    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();
    }