QA/auxi.C File Reference

#include "TGeoManager.h"

Go to the source code of this file.

Functions
void	CloseGeoManager ()
bool	checkfile (TString fn)
FairRunAna *	initrun (TString prefix, TString outfile, int min=-1, int max=-1)
void	plotmyhistos (std::vector< TH1 * > h, int maxy=700, double asp=1.1)
void	plotmyhistos (int maxy=700, double asp=1.1)
int	writemyhistos (int maxy=800, double asp=1.1)
int	fillM (RhoCandList &l, TH1 *h)
int	fillP (RhoCandList &l, TH1 *h)

Function Documentation

bool checkfile

(

TString

fn

)

Definition at line 20 of file QA/auxi.C.

References t, and x0.

{
        bool fileok=true;
        TFile fff(fn);
        if (fff.IsZombie()) fileok=false;
        TTree *t=(TTree*)fff.Get("pndsim");
        if (t==0x0) fileok=false;

        if (!fileok) cout <<"Skipping broken file '"<<fn<<"'"<<endl;
        return fileok;
}

void CloseGeoManager

(

)

Definition at line 11 of file QA/auxi.C.

References gGeoManager.

Referenced by ana_complete(), ana_day1(), ana_multi(), anaideal_complete(), benchmark(), digi(), digi_complete(), digi_complete_tb(), digi_LMD(), pid_complete(), pid_multi(), QA_LMD(), QAmacro_drc_1(), QAmacro_drc_2(), QAmacro_drc_3(), QAmacro_drc_4(), QAmacro_emc_2(), QAmacro_emc_3(), QAmacro_emc_4(), QAmacro_fastsim_2(), QAmacro_gem_2(), QAmacro_hyp_2(), QAmacro_mvd_2(), QAmacro_mvd_ana(), QAmacro_stt_2(), QAmacro_stt_3(), reco(), reco_complete(), reco_LMD(), reco_multi(), recoideal_complete(), recoqa_ca(), run_digi(), run_reco(), and run_sim().

{
  if (gROOT->GetVersionInt() >= 60602) {
    gGeoManager->GetListOfVolumes()->Delete();
    gGeoManager->GetListOfShapes()->Delete();
    delete gGeoManager;
  }
}

int fillM	(	RhoCandList &	l,
		TH1 *	h
	)

Definition at line 139 of file QA/auxi.C.

References RhoCandList::GetLength(), and i.

Referenced by tut_ana(), tut_ana_comb(), and tut_ana_pid().

{
        for (int i=0;i<l.GetLength();++i) h->Fill(l[i]->M());

        return l.GetLength();
}

int fillP	(	RhoCandList &	l,
		TH1 *	h
	)

Definition at line 146 of file QA/auxi.C.

References RhoCandList::GetLength(), and i.

{
        for (int i=0;i<l.GetLength();++i) h->Fill(l[i]->P());

        return l.GetLength();
}

FairRunAna* initrun	(	TString	prefix,
		TString	outfile,
		int	min = -1,
		int	max = -1
	)

Definition at line 32 of file QA/auxi.C.

References fRun, rtdb, and TString.

Referenced by tut_ana(), tut_ana_comb(), tut_ana_fit(), tut_ana_mclist(), tut_ana_mcmatch(), tut_ana_ntp(), tut_ana_ntp_noqa(), tut_ana_ntp_qa(), tut_ana_pid(), and tut_ana_task().

{
        // *** the files coming from the simulation
        TString inPidFile  = prefix+"_pid.root";    // this file contains the PndPidCandidates and McTruth
        TString inParFile  = prefix+"_par.root";

        // *** PID table with selection thresholds; can be modified by the user
        TString pidParFile = TString(gSystem->Getenv("VMCWORKDIR"))+"/macro/params/all.par";

        // *** initialization
        FairLogger::GetLogger()->SetLogToFile(kFALSE);
        FairRunAna* fRun = new FairRunAna();
        FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
        fRun->SetSource(new FairFileSource(inPidFile));

        // *** setup parameter database
        FairParRootFileIo* parIO = new FairParRootFileIo();
        parIO->open(inParFile);
        FairParAsciiFileIo* parIOPid = new FairParAsciiFileIo();
        parIOPid->open(pidParFile.Data(),"in");

        rtdb->setFirstInput(parIO);
        rtdb->setSecondInput(parIOPid);
        rtdb->setOutput(parIO);

        fRun->SetOutputFile(outfile);

        return fRun;
}

void plotmyhistos	(	std::vector< TH1 * >	h,
		int	maxy = 700,
		double	asp = 1.1
	)

Definition at line 62 of file QA/auxi.C.

References c, i, and sqrt().

Referenced by plotmyhistos(), tut_ana(), tut_ana_comb(), tut_ana_fast(), tut_ana_fit(), tut_ana_mcmatch(), and tut_ana_pid().

{
        int N = h.size();
        if (N<=0) return;

        int nx=sqrt(N);
        int ny=nx;
        while(nx*ny<N) nx++;

        int dim = maxy/ny;
        if (dim>400) dim=400;

        TCanvas *c=new TCanvas("c","c",10,10,dim*nx*asp, dim*ny);

        c->Divide(nx,ny, 0.0005,0.0005);

        for (int i=0;i<N;++i)
        {
                c->cd(i+1);
                h[i]->Draw();
        }
}

void plotmyhistos	(	int	maxy = 700,
		double	asp = 1.1
	)

Definition at line 85 of file QA/auxi.C.

References h, hist, i, next, plotmyhistos(), and TString.

{
    std::vector<TH1*> h;

    //somehow either one or the other is filled, depending on
    // whether a TFile was opened, or we write to one
    int Nkey = gDirectory->GetListOfKeys()->GetSize();
    int Nlst = gDirectory->GetList()->GetSize();

    if( Nkey==0 && Nlst==0) return;

    if (Nkey>0)
    {
        TIter next(gDirectory->GetListOfKeys());
        TKey *key;
        while ((key = (TKey*)next()))
        {
            TClass *cl = gROOT->GetClass(key->GetClassName());
            if (!cl->InheritsFrom("TH1")) continue;
            TH1 *hist = (TH1*)key->ReadObj();
            h.push_back(hist);
        }
    }
    else
    {
        TList *hl=gDirectory->GetList();
        for (int i=0;i<hl->GetSize();++i)
        {
            TString cn = hl->At(i)->ClassName();
            if (cn.BeginsWith("TH1") || cn.BeginsWith("TH2")) h.push_back((TH1*)hl->At(i));
        }
    }

    plotmyhistos(h, maxy, asp);
}

int writemyhistos	(	int	maxy = 800,
		double	asp = 1.1
	)

Definition at line 121 of file QA/auxi.C.

References h, i, TString, and Write().

Referenced by tut_ana(), tut_ana_comb(), tut_ana_fast(), tut_ana_fit(), tut_ana_mcmatch(), and tut_ana_pid().

{
        TList *hl=gDirectory->GetList();
        std::vector<TH1*> h;
        for (int i=0;i<hl->GetSize();++i)
        {
                TString cn = hl->At(i)->ClassName();
                if (cn.BeginsWith("TH1") || cn.BeginsWith("TH2")) h.push_back((TH1*)hl->At(i));
        }

        int N = h.size();
        if (N==0) return 0;

        for (int i=0;i<N;++i) h[i]->Write();

        return N;
}