autocutx.C File Reference

#include <algorithm>
#include "TFile.h"
#include "TTree.h"
#include "TLeaf.h"
#include "TString.h"
#include "TH1F.h"
#include "TCanvas.h"
#include "TROOT.h"
#include "TEventList.h"
#include "TDirectory.h"
#include <iostream>
#include "TLine.h"
#include "TLatex.h"
#include "TStyle.h"
#include "TObjArray.h"
#include "TPRegexp.h"
#include "TRegexp.h"
#include <map>
#include <utility>

Go to the source code of this file.

Macros
#define	MAX   1000
#define	BINS   500
#define	NCAN   3
#define	MAXVARS   8
#define	MAXCUTS   12

Typedefs
typedef std::vector< pair < double, int > >	ValueMap
typedef std::map< int, int >	CountMap

Functions
int	gettype (TTree *t, TString varname)
bool	mycompare (int i, int j)
int	init (TTree *t)
int	uid (int lev, int lrun, int lmode)
int	countEvents (TTree *t, TEventList &el)
void	makeMaps (TTree *t, TString varname, TEventList &els, TEventList &elb, ValueMap &sig, ValueMap &bg, int id)
double	bestEffEvt (TTree *t, TString varname, TEventList &els, TEventList &elb, double &bestcut, double supr, int id)
double	bestSuppressionEvt (TTree *t, TString varname, TEventList &els, TEventList &elb, double &bestcut, double eff, int id)
int	findcut (TTree *t, TEventList &els, TEventList &elb, double supr, double &bestqa)
int	countVars (TString s)
int	countCuts (TString s)
int	autocutx (TString fname, TString precut="", double supr=0.95, double target=0.0001, double mineff=0.1, double minreleff=0.0, int evmult=10000, double norm=1.0, int n0s=-1)

Variables
CountMap	evcnt
CountMap	evcntrec [10]
std::map< TString, TString >	mctvar
TString	vars [MAX]
TString	cuts [MAX]
double	cut [MAX]
double	qual [MAX]
int	idx [MAX]
double	N0_sig
double	N0_bg
double	Nsigev
double	Nbgev
bool	dstarmode
Float_t	fbranch [MAX]
Int_t	ibranch [MAX]
Bool_t	bbranch [MAX]
Int_t	ev
Int_t	run
Int_t	mode
Int_t	rec
Int_t	nbranch

Macro Definition Documentation

#define BINS   500

Definition at line 23 of file autocutx.C.

Referenced by bestCombiEvt(), and cutfinderx().

#define MAX   1000

Definition at line 22 of file autocutx.C.

#define MAXCUTS   12

Definition at line 26 of file autocutx.C.

#define MAXVARS   8

Definition at line 25 of file autocutx.C.

#define NCAN   3

Definition at line 24 of file autocutx.C.

Typedef Documentation

typedef std::map<int, int> CountMap

Definition at line 29 of file autocutx.C.

typedef std::vector<pair<double, int> > ValueMap

Definition at line 28 of file autocutx.C.

Function Documentation

int autocutx	(	TString	fname,
		TString	precut = "",
		double	supr = 0.95,
		double	target = 0.0001,
		double	mineff = 0.1,
		double	minreleff = 0.0,
		int	evmult = 10000,
		double	norm = 1.0,
		int	n0s = -1
	)

Definition at line 401 of file autocutx.C.

References b, c1, cnt, countCuts(), countEvents(), cuts, dstarmode, elb, els, evcntrec, f, findcut(), i, idx, init, mycompare(), N0_bg, N0_sig, Nbgev, nbranch, nsig, Nsigev, printf(), qual, s, t, TString, and vars.

{
//      gStyle->SetTitleX(0.2);
//      gStyle->SetTitleY(0.993);
//      gStyle->SetTitleH(0.07);
//      gStyle->SetPadTopMargin(0.08);
//      gStyle->SetTitleBorderSize(0);
        gStyle->SetOptStat(0);
        
        TCanvas *c1=new TCanvas("c1","c1",10,10,1800,550);
        c1->Divide(7,2);
        
        TString tagcut = "tag";
        
        TRegexp ren("M[0-9][0-9][0-9]"); 
        TRegexp regntp("n[0-9]+");
        TRegexp regS("[0-9]+S");
        TRegexp regB("[0-9]+B");
        
        TString en = fname(ren);  en = en(1,5); 
        TString s  = fname(regS); s  = s(0,s.Length()-1);
        TString b  = fname(regB); b  = b(0,b.Length()-1);
        
        TString ntp   = fname(regntp);
        TString smode = ntp(1,ntp.Length());
        int dmode     = smode.Atoi()/10;
        dstarmode     = (dmode==11 || dmode==13 || dmode==15);
        
        if (dmode>=10 && dmode<20) norm = 0.5;
        
        if (dmode/10 == 6) norm = 0.2;
        
        if (n0s<0) n0s = s.Atoi();
        int n0b = b.Atoi();
        
        N0_sig = n0s*evmult;
        N0_bg = n0b*evmult;
        
        cout <<"tree '"<<ntp.Data()<<"' with N_sig = "<<N0_sig<<", N_bg = "<<N0_bg<<endl;

        TFile *f=new TFile(fname,"READ");
        TTree *t=(TTree*)f->Get(ntp);
        
        init(t);
        
        TEventList elsall("elsall"), elball("elball");
        TEventList els("els"), elb("elb");
        
        if (precut=="") precut = tagcut;
        else precut = tagcut+"&&"+precut;
        
        TString sigcut = "xmct";
        TString bgcut  = "!xmct";
        
        t->Draw(">>elsall",tagcut+"&&"+sigcut);
        t->Draw(">>elball",tagcut+"&&"+bgcut);
        float nsig = countEvents(t,elsall);
        float nbg = countEvents(t,elball);
        
        double beff=1., seff=1., rseff=1.;
        
        TRegexp rnum("[\\-]*[0-9]+\\.[0-9]+$");
        int cnt=0;
        
        while (cnt<150 && beff>target && seff>=mineff && rseff>=minreleff) 
        {               
                c1->cd(cnt%7+1); 
                t->SetLineColor(1);     t->Draw("xm",precut);
                t->SetLineColor(2);     t->Draw("xm",precut+"&&"+sigcut,"same");
                t->SetLineColor(4);     t->Draw("xm",precut+"&&"+bgcut,"same");
                
                c1->Update();
                                
                t->Draw(">>els",sigcut+"&&"+precut);
                t->Draw(">>elb",bgcut+"&&"+precut);

                Nbgev  = countEvents(t, elb);
                Nsigev = countEvents(t, els);
                
                beff  = Nbgev/N0_bg;
                seff  = Nsigev/N0_sig;
                rseff = Nsigev/nsig;
                
                printf("S=%6d  B=%6d (%.2f): ",(int)Nsigev, (int)Nbgev, supr);
                
                // stop if target reached or signal efficiency too small
                if (beff<=target || seff<mineff || rseff<minreleff) continue;
                
                double qa=0.;
                
                int bestid = findcut(t, els, elb, supr, qa);
                
                TString bestcut = cuts[bestid];
                
                std::vector<int> myidx (idx, idx+nbranch);      
                std::sort(myidx.begin(), myidx.end(), mycompare);       
                
                if (precut.Length()>3 && precut.Index(vars[bestid]+">")<0 && precut.Index(vars[bestid]+"<")<0)
                {
                        int i = 0, ncuts = countCuts(precut);
                        
                        while (i<nbranch && precut.Index(vars[idx[i]]+">")<0 && precut.Index(vars[idx[i]]+"<")<0) ++i;
                        
                        if (qual[idx[i]]/qa>0.8 ) bestcut = cuts[idx[i]];
                }               
                
                TString thenum = bestcut(rnum);
                double cutval = thenum.Atof();
                TString thevar = bestcut(0,bestcut.Index(thenum));
                TString thebarevar = thevar;
                thebarevar.ReplaceAll(">","");
                thebarevar.ReplaceAll("<","");
                thebarevar.ReplaceAll("=","");
                                
                TRegexp rcut(thevar+"[\\-]*[0-9]+\\.[0-9]+");
                TString theoldcut = precut(rcut);
                
                c1->cd(cnt%7+8);
                
                t->SetLineColor(1);     t->Draw(thebarevar,precut);
                t->SetLineColor(2);     t->Draw(thebarevar,precut+"&&"+sigcut,"same");
                t->SetLineColor(4);     t->Draw(thebarevar,precut+"&&"+bgcut,"same");
                c1->Update();
                TLine l; l.SetLineColor(6); l.SetLineWidth(2);
                l.DrawLine(cutval,0,cutval,gPad->GetUymax()*0.9);
                
                c1->Update();
                
                if (theoldcut!="") 
                        precut.ReplaceAll(theoldcut,bestcut);
                else
                        precut += "&&"+bestcut;         
                        
                cout <<" : "<<precut<<endl;
                cnt++;
        }
        
        cout <<"\nCUT : "<<en<<smode<<" : "<<precut.Data()<<endl;
        cout <<"SIG EVT: "<<Nsigev<<" ev    BG: "<<Nbgev<<" ev"<<endl;

        printf("SIG EFF : %6.1f%%    BG EFF : %7.3f%%\n",Nsigev/N0_sig*100.,Nbgev/N0_bg*100.);
        printf("SIG REL : %6.1f%%    BG REL : %7.3f%%\n\n",Nsigev/nsig*100.,Nbgev/nbg*100.);

        cout<<"Recoil : "; for (int k=0;k<10;++k) printf("    %02d ",k);cout <<endl;
        cout<<"Eff    : "; for (int k=0;k<10;++k) printf("%6.1f%%",(double) evcntrec[k].size()/evmult*norm*100.);cout <<endl<<endl;

        t->SetEventList(0);
  return 0;
}

double bestEffEvt	(	TTree *	t,
		TString	varname,
		TEventList &	els,
		TEventList &	elb,
		double &	bestcut,
		double	supr,
		int	id
	)

Definition at line 189 of file autocutx.C.

References cuts, fabs(), gettype(), i, makeMaps(), Nbgev, Nsigev, qual, and TString.

Referenced by cutfinder_toyev(), cutfinderx(), and findcut().

{       
        CountMap bgcnt, sigcnt, sigcnt2;        
        ValueMap sigvals, bgvals;
        
        int Nbgval = elb.GetN();
        int dtype = gettype(t,varname);
        
        makeMaps(t, varname, els, elb, sigvals, bgvals, id);
        
        
        // sort signals by first element = variable value
        sort(bgvals.begin(), bgvals.end());

        int i=0;
        bgcnt.clear();
        // find cut for current efficiency requirement left cut
        while ( (bgcnt.size()/(double)Nbgev)<(1-supr) && i<Nbgval ) bgcnt[bgvals[i++].second]+=1;
        double leftcut = bgvals[i].first;
        
        bgcnt.clear();
        i=bgvals.size()-1;
        // find cut for current efficiency requirement right cut
        while ( (bgcnt.size()/(double)Nbgev)<(1-supr) && i>=0) bgcnt[bgvals[i--].second]+=1; 
        double rightcut = bgvals[i].first;
        
        sigcnt.clear();
        sigcnt2.clear();
        
        for (i=0;i<(int)sigvals.size();++i)
        {
                if (sigvals[i].first<=leftcut)  sigcnt[sigvals[i].second]+=1;
                if (sigvals[i].first>=rightcut) sigcnt2[sigvals[i].second]+=1;
        }
        
        t->SetBranchStatus("*",1);
                
        double lefteff  = sigcnt.size()/Nsigev;
        double righteff = sigcnt2.size()/Nsigev;

        //cout <<varname<<" "<<flush;
        
        TString prec="%.4f";
        
        bestcut = rightcut;
        if (lefteff>righteff) 
        {
                bestcut = leftcut;
                
                if (bestcut<1e-3 || fabs(1.-bestcut)<1e-3) prec="%.6f";
                if (dtype==0) cuts[id] = TString::Format(TString("%s<="+prec).Data(),varname.Data(),leftcut);
                else cuts[id] = TString::Format("%s<=%.1f",varname.Data(),leftcut);
                
                qual[id] = lefteff;

                return lefteff;
        }
        
        if (bestcut<1e-3 || fabs(1.-bestcut)<1e-3) prec="%.6f";

        if (dtype==0) cuts[id] = TString::Format(TString("%s>="+prec).Data(),varname.Data(),rightcut);
        else cuts[id] = TString::Format("%s>=%.1f",varname.Data(),rightcut);

        qual[id] = righteff;

        return righteff;
}

double bestSuppressionEvt	(	TTree *	t,
		TString	varname,
		TEventList &	els,
		TEventList &	elb,
		double &	bestcut,
		double	eff,
		int	id
	)

Definition at line 259 of file autocutx.C.

References cuts, fabs(), gettype(), i, makeMaps(), Nbgev, Nsigev, qual, and TString.

Referenced by cutfinder_fullev(), cutfinder_toyev(), cutfinderx(), and findcut().

{       
        CountMap bgcnt, bgcnt2, sigcnt; 
        ValueMap sigvals, bgvals;
        
        int Nsigval = els.GetN();
        int dtype = gettype(t,varname);
        
        makeMaps(t, varname, els, elb, sigvals, bgvals, id);
        
        // sort signals by first element = variable value
        sort(sigvals.begin(), sigvals.end());

        int i=0;
        sigcnt.clear();
        // find cut for current efficiency requirement left cut
        while ( (sigcnt.size()/(double)Nsigev)<eff && i<Nsigval ) sigcnt[sigvals[i++].second]+=1;
        double leftcut = sigvals[i].first;
        
        sigcnt.clear();
        i=sigvals.size()-1;
        // find cut for current efficiency requirement right cut
        while ( (sigcnt.size()/(double)Nsigev)<eff && i>=0) sigcnt[sigvals[i--].second]+=1; 
        double rightcut = sigvals[i].first;
        
        bgcnt.clear();
        bgcnt2.clear();
        
        for (i=0;i<(int)bgvals.size();++i)
        {
                if (bgvals[i].first<=leftcut)  bgcnt[bgvals[i].second];
                if (bgvals[i].first>=rightcut) bgcnt2[bgvals[i].second];
        }
        
        t->SetBranchStatus("*",1);
        
        
        double leftsupr = (Nbgev-bgcnt.size())/Nbgev;
        double rightsupr = (Nbgev-bgcnt2.size())/Nbgev;

        //cout <<varname<<" "<<flush;
        
        if (leftcut!=leftcut || rightcut!=rightcut) return 0;
        
        TString prec="%.4f";
        
        bestcut = rightcut;
        if (leftsupr>rightsupr) 
        {
                bestcut = leftcut;
                
                if (bestcut<1e-3 || fabs(1.-bestcut)<1e-3) prec="%.6f";
                if (dtype==0) cuts[id] = TString::Format(TString("%s<="+prec).Data(),varname.Data(),leftcut);
                else cuts[id] = TString::Format("%s<=%.1f",varname.Data(),leftcut);
                
                qual[id] = leftsupr;

                return leftsupr;
        }
        
        if (bestcut<1e-3 || fabs(1.-bestcut)<1e-3) prec="%.6f";

        if (dtype==0) cuts[id] = TString::Format(TString("%s>="+prec).Data(),varname.Data(),rightcut);
        else cuts[id] = TString::Format("%s>=%.1f",varname.Data(),rightcut);

        qual[id] = rightsupr;

        return rightsupr;
}

int countCuts

(

TString

s

)

Definition at line 393 of file autocutx.C.

Referenced by autocutx().

{
        TObjArray *tok = s.Tokenize("&&");
        return tok->GetEntries();       
}

int countEvents	(	TTree *	t,
		TEventList &	el
	)

Definition at line 129 of file autocutx.C.

References ev, evcnt, evcntrec, i, mode, rec, run, and uid().

Referenced by autocutx(), cutfinder(), cutfinder_fullev(), cutfinder_toyev(), cutfinderx(), cutqa(), drawHistos(), showvarexample(), simubg(), and TMVAApply().

{
        t->SetEventList(0);
        t->SetBranchStatus("*",0);
        t->SetBranchStatus("ev",1);
        t->SetBranchStatus("run",1);
        t->SetBranchStatus("mode",1);
        t->SetBranchStatus("recmode",1);
                
        evcnt.clear();
        for (int j=0;j<10;++j) evcntrec[j].clear();
        
        for (int i=0;i<el.GetN();++i)
        {
                t->GetEntry(el.GetEntry(i));
                evcnt[uid(ev,run,mode)]+=1;
                if (rec<10) evcntrec[rec][uid(ev,run,mode)]+=1;
        }
        t->SetBranchStatus("*",1);

        return evcnt.size();
}

int countVars

(

TString

s

)

Definition at line 368 of file autocutx.C.

References cnt, i, and TString.

{
        int cnt = 0;
        TString svar="";
        TRegexp rvar("-[a-zA-Z_][a-zA-Z0-9_]+-");
        
        TObjArray *tok = s.Tokenize("&&");
        int N = tok->GetEntries();      
        for (int i=0;i<N;++i) 
        {
                if (i<100) 
                {
                        TString toks = ((TObjString*)tok->At(i))->String();
                        toks.ReplaceAll("\t","");
                        toks = toks.Strip(TString::kBoth);
                        if (toks.Index("<")>0) toks = toks(0,toks.Index("<"));
                        if (toks.Index(">")>0) toks = toks(0,toks.Index(">"));
                        svar += "-"+toks+"-";
                }
        }
        while (svar(rvar)!="") {svar.ReplaceAll(svar(rvar),""); cnt++;}
        return cnt;
}

int findcut	(	TTree *	t,
		TEventList &	els,
		TEventList &	elb,
		double	supr,
		double &	bestqa
	)

Definition at line 331 of file autocutx.C.

References bestEffEvt(), bestSuppressionEvt(), cut, cuts, fabs(), i, idx, nbranch, TString, and vars.

Referenced by autocutx().

{
        int i;
        
        bestqa = -999.;
        
        TString bestcut="";
        
        int bestid = -1;
        
        //cout <<"000/000"<<flush;

        for(i=0; i<nbranch; ++i)
        {
                idx[i] = i;
                
                if (i%(nbranch/10)==0){ cout <<"#"<<flush;}
//              printf("\b\b\b\b\b\b\b%03d/%03d",i+1,nbranch);
//              fflush(stdout);
                
                double qa=0;
                
                if (supr<0)     qa = bestEffEvt(t, vars[i], els, elb, cut[i], -supr, i);
                else        qa = bestSuppressionEvt(t, vars[i], els, elb, cut[i], supr, i);     
                
                if (qa>bestqa && fabs(cut[i])>1e-5 && fabs(cut[i]-1.0)>1e-5) 
                {
                        bestqa  = qa;
                        bestcut = cuts[i];
                        bestid = i;
                }
        }
        
        return bestid;
}

int gettype	(	TTree *	t,
		TString	varname
	)

Definition at line 51 of file autocutx.C.

References TString.

Referenced by bestCombiEvt(), bestEffEvt(), bestSuppressionEvt(), init(), makeMaps(), TMVATrainer(), TMVATrainer_608(), and TMVATraining().

{
        TString leaftype = t->GetLeaf(varname)->GetTypeName();
        
        if (leaftype=="Float_t")         return 0; 
        else if (leaftype=="Int_t")      return 1;  
        else if (leaftype=="Bool_t") return 2;
        
        return -1;
}

int init

(

TTree *

t

)

Definition at line 71 of file autocutx.C.

References bbranch, dstarmode, ev, fbranch, gettype(), i, ibranch, mode, nbranch, ok, rec, run, TString, v, and vars.

{
        t->SetBranchAddress("ev",&ev);
        t->SetBranchAddress("run",&run);
        t->SetBranchAddress("mode",&mode);
        t->SetBranchAddress("recmode",&rec);
        
        TObjArray* branches = t->GetListOfBranches();   
        
        nbranch = 0;
        
        for(int i=0; i<=branches->GetLast(); ++i)
        {
                TBranch* branch = (TBranch*)branches->UncheckedAt(i);
                TString v=branch->GetName(); 

                if ( v=="ev" || v=="mode" || v=="run" || v=="nsig" ) continue;
                if ( v.Contains("pdg") || v.Contains("beam") || v.Contains("mct") ) continue;
                if ( v.BeginsWith("t") && v!="thr") continue;
                if ( v.EndsWith("vx") || v.EndsWith("vy") || v.EndsWith("vz") || v.EndsWith("pocmag")) continue;

                bool ok=false;
                
                if ( v.BeginsWith("es") || v=="mmiss" || v.EndsWith("d0m") ) ok=true;
                if ( v.EndsWith("p") || v.EndsWith("tht") || v.EndsWith("pcm") || v.EndsWith("thtcm") || v.EndsWith("pt") ) ok=true;
                if ( ok || v.EndsWith("ang") || v.Contains("poc") ) ok=true;
                if ( ok || v.Contains("pid") || v.Contains("min") || v.Contains("max") || v.Contains("sum") || v.Contains("fw") ) ok=true;
                if ( ok || v.EndsWith("sph") || v.EndsWith("apl") || v.EndsWith("pla") || v.EndsWith("thr") || v.EndsWith("cir") ) ok=true;
                
                if ( v=="xmdif" && dstarmode ) ok = true;
                if ( v=="xd0m" && dstarmode ) ok = false;
                
                if (!ok) continue;
                
                Int_t dtype=gettype(t, v);
                switch (dtype)
                {
                        case 0: t->SetBranchAddress(v, &(fbranch[nbranch])); break;
                        case 1: t->SetBranchAddress(v, &(ibranch[nbranch])); break;
                        case 2: t->SetBranchAddress(v, &(bbranch[nbranch])); break;
                }
                
                vars[nbranch] = v;
                nbranch++;
        }
        
        return nbranch;
}

void makeMaps	(	TTree *	t,
		TString	varname,
		TEventList &	els,
		TEventList &	elb,
		ValueMap &	sig,
		ValueMap &	bg,
		int	id
	)

Definition at line 153 of file autocutx.C.

References bbranch, ev, fbranch, gettype(), i, ibranch, mode, run, and uid().

Referenced by bestCombiEvt(), bestEffEvt(), and bestSuppressionEvt().

{
        int i;
        t->SetBranchStatus("*",0);
        t->SetBranchStatus("ev",1);
        t->SetBranchStatus("mode",1);
        t->SetBranchStatus("run",1);
        t->SetBranchStatus(varname,1);
        
        Float_t var;
        Int_t dtype=gettype(t, varname);
        
        sig.clear();
        bg.clear();

        // prepare the pairs of variable value, eventnumber for signal and count signals
        for (i=0;i<els.GetN();++i)
        {
                t->GetEntry(els.GetEntry(i));
                switch (dtype) {case 2: var = (Float_t)bbranch[id]; break; case 1: var = (Float_t) ibranch[id]; break; default: var=fbranch[id]; }
                
                sig.push_back(std::make_pair(var, uid(ev,run,mode)  ));
        }
                
        // prepare the pairs of variable value, eventnumber for background and count background
        for (i=0;i<elb.GetN();++i)
        {
                t->GetEntry(elb.GetEntry(i));
                switch (dtype) {case 2: var = (Float_t)bbranch[id]; break; case 1: var = (Float_t) ibranch[id]; break; default: var=fbranch[id]; }
                
                bg.push_back(std::make_pair(var, uid(ev,run,mode) ));
        }
}

bool mycompare	(	int	i,
		int	j
	)

Definition at line 64 of file autocutx.C.

References i, and qual.

Referenced by autocutx(), cutfinder(), cutfinder_fullev(), cutfinder_toyev(), and cutfinderx().

{
  return qual[i]>qual[j];
}

int uid	(	int	lev,
		int	lrun,
		int	lmode
	)

Definition at line 122 of file autocutx.C.

Referenced by PndDrcReco::CanvasSave(), RhoFitterBase::Chi2Contribution(), countevents(), countEvents(), PndSimpleAnalysis::FillGenericLists(), PndAnaTask::FillList(), makeMaps(), RhoCandidate::SetUid(), and TMVAApply().

{
        return lev+10000*lrun+(((lmode/100)%10)*20+lmode%10)*100000;
}

Variable Documentation

Bool_t bbranch[MAX]

Definition at line 46 of file autocutx.C.

Referenced by init(), makeMaps(), TMVAApply(), and TMVATester().

double cut[MAX]

Definition at line 36 of file autocutx.C.

Referenced by analyse_etac1(), analyse_etac1_slc(), analyse_J(), analyse_J_slc(), analyse_phi_base(), PndTrkCleanup::BadTrack_ParStt(), checkphsp2_2(), Config(), CreateChannelHistos(), cutfinder(), cutfinder_fullev(), cutfinder_toyev(), cutfinderx(), cutqa(), FCDrawVariable(), findcut(), findcuts(), PndFileNameCreator::GetCustomFileName(), histos_ca(), PndSimpleCombinerTask::InitParms(), QA_histos(), readConfig(), RhoPdtLoader::ReadPDGTable(), PndSoftTriggerTask::SelectPid(), PndSoftTriggerTask::SelectPidProb(), PndSoftTriggerTask::SetDstMDiffCut(), PndSoftTriggerTask::SetInitialPidCut(), PndSciTAnaIdeal::SetLHKaonCut(), PndSdsStripCorrelator::Setup(), PndMvdRiemannVertexFinderTask::SetVertexCut(), PndRiemannTrack::SetVertexCut(), simubg(), PndTrkCleanup::SttSkewCleanup(), suppression(), and TMVAApply().

TString cuts[MAX]

Definition at line 35 of file autocutx.C.

Referenced by PndFilteredPrimaryGenerator::AddFilter(), autocutx(), bestCombiEvt(), bestEffEvt(), bestSuppressionEvt(), Config(), cutfinderx(), findcut(), full_core_ntp(), get_tuples(), init_cuts(), init_cuts_1000(), init_cuts_95(), init_cuts_old(), PndSoftTriggerTask::ReadConfiguration(), readmodes(), PndFtfDirect::Setup(), simu_tag(), simubg(), and toy_core_ntp().

bool dstarmode

Definition at line 42 of file autocutx.C.

Referenced by autocutx(), cutfinderx(), and init().

Int_t ev

Definition at line 47 of file autocutx.C.

CountMap evcnt

Definition at line 31 of file autocutx.C.

Referenced by cnt_events(), countevents(), countEvents(), full_core_ntp(), simu_tag(), and toy_core_ntp().

CountMap evcntrec[10]

Definition at line 31 of file autocutx.C.

Referenced by autocutx(), countEvents(), and cutfinderx().

Float_t fbranch[MAX]

Definition at line 44 of file autocutx.C.

Referenced by init(), makeMaps(), TMVAApply(), and TMVATester().

Int_t ibranch[MAX]

Definition at line 45 of file autocutx.C.

Referenced by init(), makeMaps(), TMVAApply(), and TMVATester().

int idx[MAX]

Definition at line 38 of file autocutx.C.

Referenced by PndFastSim::acceptFilters(), ana_dsdsj2(), ana_dsdsj2_EvtLoop(), ana_dsdsj_full(), autocutx(), bestEffEvt(), bestSeparation(), bestSepEvt(), bestSuppressionEvt(), cutfinder(), cutfinder_fullev(), cutfinder_toyev(), cutfinderx(), SimpleCand::Dau(), PndSmpCand::Dau(), PndSoftTriggerTask::DoCombinatorics(), PndMdtHitProducer::Exec(), PndEmcHitProducer::Exec(), PndEmcAnalysis::Exec(), findcut(), PndEmcMakeClusterOnline::FinishCluster(), PndEmcDigiCalibrator::GetIdxByEnergy(), PndEmcDigiCalibrator::GetIdxByEnergyForBump(), PndSimpleCombiner::InitDecayInfo(), PndSttMvdGemTracking::OrderGemHits(), PndFsmMdtPid::respond(), PndFsmEmcPid::respond(), Photospp::PH_HEPEVT_Interface::set(), PndPidProbability::SetIndex(), PndFastSim::SetInvMassFilter(), PndPidCandidate::SetMcIndex(), SimpleCand::SetMotherIdx(), PndPidBremCorrected4Mom::SetPidCandIdx(), PndPidCandidate::SetTrackBranch(), PndPidCandidate::SetTrackIndex(), and PndSimpleAnalysis::SetupAnalysis().

std::map<TString, TString> mctvar

Definition at line 32 of file autocutx.C.

Referenced by fillHistos(), and init_modes().

Int_t mode

Definition at line 47 of file autocutx.C.

Referenced by PndSoftTriggerTask::AcceptCandidate(), compress(), countevents(), countEvents(), crosstag(), cutfinderx(), cutqa(), evaltrig(), full_core_ntp(), getpars(), init(), inspectmode(), inspectvar(), makeMaps(), PndParticleQATask::PndParticleQATask(), PndProdAnaTask::PndProdAnaTask(), prod_ana(), prod_fsim(), prod_myana(), prod_sof(), quickfsimana(), PndSoftTriggerTask::ReadTriggerLines(), run0AllSysDPMbkg(), run_sim_stt_dpm(), run_sim_sttcombi_dpm(), run_sim_tpc_dpm(), run_sim_tpccombi_dpm(), runLumi0DPMbkg(), runLumi0DPMbkgSte(), runLumiPixel0SimDPM(), runLumiPixel0SimDPMDirect(), runLumiPixel0SimDPMDirect_all(), runLumiPixel0SimFTF(), runMC_dpm(), PndPidCorrelator::SetDrcMode(), PndPidCorrelator::SetDskMode(), PndPidCorrelator::SetEmcMode(), PndPidCorrelator::SetFtofMode(), PndPidCorrelator::SetFtsMode(), PndPidCorrelator::SetMdtMode(), PndSoftTriggerLine::SetModeCode(), PndPidCorrelator::SetMvdMode(), PndOnlineFilterInfo::SetNTag(), PndKFParticleFinderPID::SetPIDMode(), PndMdtDigiProducer::SetStripMode(), PndPidCorrelator::SetSttMode(), PndSoftTriggerTask::SetTagNSigMode(), PndPidCorrelator::SetTofMode(), PndSdsStripCorrelator::Setup(), sim_LMD(), simfast_dpm(), simfast_dpm_cmp(), softtrigger_kin5(), PndSoftTriggerTask::TagMode(), TMVAApply(), toy_core(), toy_core_ntp(), tut_fastsim(), and UserDecayConfig().

double N0_bg

Definition at line 40 of file autocutx.C.

Referenced by autocutx(), cutfinder(), cutfinder_fullev(), cutfinder_toyev(), and cutfinderx().

double N0_sig

Definition at line 40 of file autocutx.C.

Referenced by autocutx(), cutfinder(), cutfinder_fullev(), cutfinder_toyev(), and cutfinderx().

double Nbgev

Definition at line 41 of file autocutx.C.

Referenced by autocutx(), bestCombiEvt(), bestEffEvt(), bestSeparation(), bestSepEvt(), bestSuppressionEvt(), cutfinder_fullev(), cutfinder_toyev(), and cutfinderx().

Int_t nbranch

Definition at line 47 of file autocutx.C.

Referenced by autocutx(), cutfinderx(), findcut(), and init().

double Nsigev

Definition at line 41 of file autocutx.C.

Referenced by autocutx(), bestCombiEvt(), bestEffEvt(), bestSeparation(), bestSepEvt(), bestSuppressionEvt(), cutfinder_fullev(), cutfinder_toyev(), and cutfinderx().

double qual[MAX]

Definition at line 37 of file autocutx.C.

Referenced by autocutx(), bestEffEvt(), bestSuppressionEvt(), and mycompare().

Int_t rec

Definition at line 47 of file autocutx.C.

Referenced by countevents(), countEvents(), init(), PndCANPlets::PndCANPlets(), and DecayTreeFitter::RecoPhoton::updCache().

Int_t run

Definition at line 47 of file autocutx.C.

Referenced by bestSeparation(), bestSuppressionEvt(), countevents(), countEvents(), init(), makeMaps(), Next(), PndGemMonitor::PndGemMonitor(), PndGeoHandling::PndGeoHandling(), prod_ana(), prod_fsim(), prod_sof(), qa_softtrig(), quickana(), quickfsimana(), PndEmcReadProto192Data::SetParContainers(), PndEmcReadProtoData::SetParContainers(), PndPidFtofAssociatorTask::SetParContainers(), PndMdtTrkFinder::SetParContainers(), PndMdtTrkProducer::SetParContainers(), PndMcListConverter::SetParContainers(), PndEmcFWEndcapTimebasedWaveforms::SetParContainers(), PndEmcMakeDigi::SetParContainers(), PndEmcMakeRecoHit::SetParContainers(), PndEmc2DLocMaxFinder::SetParContainers(), PndJPsiExample::SetParContainers(), PndTutAnaTask::SetParContainers(), PndEventBuilderManager::SetParContainers(), PndEmcHitsToWaveform::SetParContainers(), PndTutAnaTaskD0::SetParContainers(), PndEmcFWEndcapDigi::SetParContainers(), PndTutThaiTask::SetParContainers(), PndEmcMakeCluster::SetParContainers(), PndGemMatchHits::SetParContainers(), PndGeoHandling::SetParContainers(), PndGemDigiAna::SetParContainers(), PndScrutAnaTask::SetParContainers(), PndGemIdealHitProducer::SetParContainers(), PndGemEmpty::SetParContainers(), PndMdtPointsToWaveform::SetParContainers(), PndEventBuilderOnMCEvents::SetParContainers(), PndEmcExpClusterSplitter::SetParContainers(), PndGemFindHitsAna::SetParContainers(), PndEmcCorrBump::SetParContainers(), PndEmcMakeCorr::SetParContainers(), PndEmcMakeClusterOnline::SetParContainers(), PndEmcTmpWaveformToDigi::SetParContainers(), PndEmcPhiBumpSplitter::SetParContainers(), PndAnalysisTaskExample::SetParContainers(), PndEmcAnalysis::SetParContainers(), PndGemEventBuilderOnTracks::SetParContainers(), PndAnalysisTask::SetParContainers(), PndTripleAnaTask::SetParContainers(), PndGemFindClusters::SetParContainers(), PndMyAnalysisTask::SetParContainers(), PndEmcHitProducer::SetParContainers(), PndEmcWaveformToDigi::SetParContainers(), PndEmcWaveformToCalibratedDigi::SetParContainers(), PndHypFullIdealAna::SetParContainers(), PndGemDigitize::SetParContainers(), PndHypFullAna::SetParContainers(), PndHypSimpleAna::SetParContainers(), PndEmcMultiWaveformToCalibratedDigi::SetParContainers(), PndDrcRecoLookupMapS::SetParContainers(), PndGemFindHits::SetParContainers(), PndGemFindClustersTB::SetParContainers(), PndGemFindHitsQA::SetParContainers(), PndGemMagneticFieldVsTrackParameters::SetParContainers(), PndGemTrackFinderOnHits::SetParContainers(), PndMvdGemTrackFinderOnHits::SetParContainers(), PndGemTrackFinderQA::SetParContainers(), PndGemTrackFinderOnHitsTB::SetParContainers(), PndDrcRecoLookupMap::SetParContainers(), PndSimpleAnalysis::SetParContainers(), PndFastSim::SetParContainers(), PndLLbarAnaTask::SetParContainers(), PndBarrelTrackFinder::SetParContainers(), PndPidCorrelator::SetParContainers(), PndSoftTriggerTask::SetParContainers(), and TMVAApply().

TString vars[MAX]

Definition at line 34 of file autocutx.C.

Referenced by autocutx(), cutfinder(), cutfinder_fullev(), cutfinder_toyev(), cutfinderx(), dalitz(), findcut(), init(), makefml(), plotTMVADemoData(), pz(), TMVATester(), TMVATrainer(), and TMVATrainer_608().