drawMvdHistos.C File Reference

Go to the source code of this file.

Functions
double	CalcEventsPerYear (Int_t simEvents, Int_t interactionRate=2 *1E7, Double_t dutyCycle=0.5)
int	drawMvdHistos (TString fileName, TString selector, Int_t columns=2, Int_t mult=1)

Function Documentation

double CalcEventsPerYear	(	Int_t	simEvents,
		Int_t	interactionRate = 2*1E7,
		Double_t	dutyCycle = 0.5
	)

Definition at line 1 of file drawMvdHistos.C.

{
        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
	)

Definition at line 8 of file drawMvdHistos.C.

References c1, c2, CalcEventsPerYear(), f, f1, PndGeoHandling::GetVolumeID(), gGeoManager, h1, htemp, i, max(), mult, TransformHisto(), and TString.

{
        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;

}