drawMvdHistos.C

Go to the documentation of this file.

double CalcEventsPerYear(Int_t simEvents, Int_t interactionRate = 2*1E7, Double_t dutyCycle = 0.5)
{
        Int_t secondsPerYear = 60*60*24*365;
        return (double)interactionRate/simEvents * secondsPerYear * dutyCycle;
}


int drawMvdHistos(TString fileName, TString selector, Int_t columns = 2, Int_t mult = 1)
{
        gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
        gROOT->LoadMacro("$VMCWORKDIR/macro/mvd/TransformHisto.C");
        gStyle->SetPalette(1);
//      gSystem->Load("libEve");
//      gROOT->LoadMacro("DrawEveHisto.C");

        TFile f(fileName);
        TGeoManager* fGeoMan = new TGeoManager("geoMan","geoMan");
        fGeoMan->Import("Mvd_GridUrqmdSim.root");
        PndGeoHandling gH(gGeoManager);

        TString id = gH->GetVolumeID(selector);

        //TEveManager::Create();

        TList* l = f.GetListOfKeys();
        TIter iter(l);
        TIter iter2(l);
        TObject* ob;
        int NHistos = 0;
        double max = 0, oldMax = 0;

        TString path;
        while(ob = iter()){
                if (TString(ob->GetName()).Contains(id)){
                                NHistos++;
                                TH2D* histo = (TH2D*)(f.Get(ob->GetName()));
                                max = histo->GetMaximum();
                                if (max > oldMax)
                                        oldMax = max;
                }
        }
        max = oldMax*1.1;

        TCanvas *c1 = new TCanvas();
        c1->Divide(columns, NHistos/columns);
        TH1D h1("h1",selector, 1000,0,max*CalcEventsPerYear(mult, 1E6));

        int i = 0;
        while(ob = iter2()){
                if (TString(ob->GetName()).Contains(id)){
                        TH2D* histo = (TH2D*)(f.Get(ob->GetName()));
                        TF1* f1 = new TF1("f1","1",-100,100);
                        histo->Multiply(f1, CalcEventsPerYear(mult, 1E6));
                        c1->cd(++i);
                        histo->DrawCopy("colz");
                        TH1D htemp(TransformHisto(histo,0,max*CalcEventsPerYear(mult,1E6)));
                        h1.Add(&htemp);
                }
        }
        TCanvas *c2 = new TCanvas();
        c2->cd(1);
        h1.DrawCopy();
  return 0;

}