softrig/qa_softtrig.C

Go to the documentation of this file.

bool checkfile(TString fn)
{
        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 qa_softtrig(TString outpre="M9999", int mode, double pmom, int from=1, int to=1, int nEvents=0)
{
        gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
        gSystem->Load("libSofTrig");

        TString OutFile   = outpre+"_st_qa.root";
        TString inParFile = TString::Format("%s_%d_par.root",outpre.Data(),from);

        bool finefile[1000];
        for (int i=0;i<1000;++i) finefile[i]=false;

        FairRunAna *fRun= new FairRunAna();

        bool firstfile=true;

        // add pid files
        for (int i=from;i<=to;++i)
        {
                TString fname = outpre+"_"; fname+=i; fname+="_pid.root";
                //TString fname2 = outpre+"_"; fname2+=i; fname2+="_sim.root";

                if (checkfile(fname) )//&& checkfile(fname2))
                {
                        finefile[i]=true;
                        if (firstfile) fRun->SetInputFile(fname);
                        else fRun->AddFile(fname);
                        firstfile=false;
                }

        }
        FairLogger::GetLogger()->SetLogToFile(kFALSE);

        // *** initialization
/*      FairRuntimeDb* rtdb = fRun->GetRuntimeDb();

        FairParRootFileIo* parIO = new FairParRootFileIo();
        parIO->open(inParFile);
        rtdb->setFirstInput(parIO);
        rtdb->setOutput(parIO);  */

        fRun->SetOutputFile(OutFile);

        //---------------------Create and Set the Field(s)----------
        //PndMultiField *fField= new PndMultiField("FULL");
        //fRun->SetField(fField);

        RhoCalculationTools::ForceConstantBz(20.0);


        // *** HERE OUR TASK GOES!
        PndSoftTriggerTask *stTask = new PndSoftTriggerTask(pmom, mode);

        stTask->SetPi0SignalParams(0.135, 0.01);
        stTask->SetEtaSignalParams(0.547, 0.03);
        stTask->SetKs0SignalParams(0.493, 0.05);

        stTask->SetGammaMinE(0.15);
        stTask->SetTrackMinP(0.15);

        TString algo = "PidAlgoEmcBayes;PidAlgoDrc;PidAlgoDisc;PidAlgoStt;PidAlgoMdtHardCuts";

        stTask->SetAlgoElectron(algo);
        stTask->SetAlgoMuon(algo);
        stTask->SetAlgoPion(algo);
        stTask->SetAlgoKaon(algo);
        stTask->SetAlgoProton(algo);

        fRun->AddTask(stTask);

        // *** and run analysis
        fRun->Init();
        fRun->Run(0,nEvents);

        //gObjectTable->Print();
}