full_core_ntp.C File Reference

#include <iostream>
#include <fstream>
#include <vector>
#include <map>
#include "TFile.h"
#include "TTree.h"
#include "TChain.h"
#include "TRandom3.h"
#include "TDatabasePDG.h"
#include "TParticlePDG.h"
#include "TCanvas.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TLegend.h"
#include "TParticle.h"
#include "TClonesArray.h"
#include "TF1.h"
#include "TStyle.h"
#include "TNtuple.h"
#include "TLatex.h"
#include "TGaxis.h"
#include "TEventList.h"

Go to the source code of this file.

Functions
void	setStyle ()
double	getPos (TH1 *h)
void	configHisto (TH1 *h)
void	init_modes ()
void	init_cuts_1000 (TString energy)
void	init_cuts_95 (TString energy)
int	cnt_events (TChain *t, TString cut, bool globcnt=false)
void	fillHistos (TChain *t, TString var, TString cut, TString all, TString sel, TString sig)
double	DbMass (TString name)
int	full_core_ntp (TString fname, int nev=50, bool tagonly=false)

Variables
TRandom3	fRand
TLorentzVector	fIni
const int	fMAX = 100
int	nbins = 300
double	low = 0.0
double	high = 6.
std::map< int, int >	evcnt
std::map< int, int >	treeidx
std::map< TString, TString >	cuts
std::map< TString, TString >	modes
std::map< TString, TString >	mctvar
TString	ntpname [10] = {"nphi","nlam","n2e","njpsi1","njpsi2","netac2","nd01","ndpm1","nds1","nlamc"}
TString	plotname [10] = {"phim","lamm","eem","jpsim","jpsim","etacm","d0m","dpmm","dsm","lamcm"}

Function Documentation

int cnt_events	(	TChain *	t,
		TString	cut,
		bool	globcnt = false
	)

Definition at line 285 of file full_core_ntp.C.

References ev, evcnt, and treeidx.

Referenced by full_core_ntp().

{
        treeidx.clear();
        Float_t ev;
        t->Draw(">>el",cut);
        
        TEventList *el=(TEventList*)gDirectory->Get("el");

        t->SetBranchStatus("*",0); 
        t->SetBranchStatus("ev",1);
        t->SetBranchAddress("ev",&ev);
        
        std::map<int, int> sigcnt;
        for (int j=0;j<el->GetN();++j) 
        {
                t->GetEntry(el->GetEntry(j));
                int tidx = t->GetTreeNumber();
                sigcnt[(int)ev+tidx*1000]+=1;
                if (globcnt) evcnt[(int)ev+tidx*1000]+=1;
        }
        t->SetBranchStatus("*",1); 
        return sigcnt.size();
}

void configHisto

(

TH1 *

h

)

Definition at line 76 of file full_core_ntp.C.

Referenced by full_core_ntp().

{
        h->SetFillStyle(3001);
        h->SetFillColor(602);
        h->SetLineStyle(3);
        h->SetLineColor(602);
/*      TGaxis *ax = (TGaxis*)h->GetYaxis();
        ax->SetMaxDigits(2);*/
}

double DbMass

(

TString

name

)

Definition at line 316 of file full_core_ntp.C.

Referenced by full_core_ntp().

{
        TDatabasePDG *fPdg = TDatabasePDG::Instance();
        if (fPdg->GetParticle(name)) return fPdg->GetParticle(name)->Mass();
}

void fillHistos	(	TChain *	t,
		TString	var,
		TString	cut,
		TString	all,
		TString	sel,
		TString	sig
	)

Definition at line 309 of file full_core_ntp.C.

References mctvar.

Referenced by full_core_ntp().

{
        t->Project(all,var,cut);
        t->Project(sel,var,cut+" && tag");
        t->Project(sig,var,cut+" && "+mctvar[t->GetName()]+"mct");
        
}

int full_core_ntp	(	TString	fname,
		int	nev = 50,
		bool	tagonly = false
	)

Definition at line 322 of file full_core_ntp.C.

References c2, cnt_events(), configHisto(), cuts, DbMass(), energy, evcnt, fillHistos(), i, init_cuts_1000(), init_modes(), mode, nbins, ntpname, offset(), out, setStyle(), t, and TString.

{
        int i,j;
        nev*=1000;
        
        setStyle();     
        init_modes();
        
        TString mode(fname(fname.Last('/')+4,3));

        TString energy(fname(fname.Last('/')+2,2));
        double sqrts=energy.Atof();
        
        if (tagonly)
                init_cuts_1000("");
        else
                init_cuts_1000(energy);
        
        cout <<energy<<endl;
        
        
        //TFile *f=new TFile(fname);
        TChain *t[10];
        
        for (i=0;i<10;++i)
        {
                t[i]=new TChain(ntpname[i]);
                t[i]->Add(fname);
/*              t[i]=(TTree*)f->Get(ntpname[i]);
                if (t==0) cout <<"not found :"<<ntpname[i].Data()<<endl;*/
        }
        
        TCanvas *c2=new TCanvas("c2","c2",20,20,1800,800);
        c2->Divide(5,2,0.005,0.01);
        
        
        double PhiMass = DbMass("phi");
        
        double D0Mass  = DbMass("D0");
        double DpmMass = DbMass("D+");
        double DsMass  = DbMass("D_s+");
        
        double JpsiMass  = DbMass("J/psi");
        double EtacMass  = DbMass("eta_c");
        double Chic0Mass = DbMass("chi_0c");
        
        double LamMass  = DbMass("Lambda0");
        double LamcMass = DbMass("Lambda_c+");
        
        // *** resonances
        double fPhiQaMin = PhiMass - 0.1;               
        double fPhiQaMax = PhiMass + 0.1;               
        
        // *** open charm
        double fD0QaMin = D0Mass - 0.3;         
        double fD0QaMax = D0Mass + 0.3;         
        
        double fDpmQaMin = DpmMass - 0.3;               
        double fDpmQaMax = DpmMass + 0.3;               
        
        double fDsQaMin = DsMass - 0.25;                
        double fDsQaMax = DsMass + 0.25;                
        
        // *** baryons
        double fLamcQaMin = LamcMass - 0.2;             
        double fLamcQaMax = LamcMass + 0.2;             

        double fLamQaMin = LamMass - 0.15;              
        double fLamQaMax = LamMass + 0.15;              
        
        // *** charmonia
        double fEtacQaMin = EtacMass - 0.65;            
        double fEtacQaMax = EtacMass + 0.65;            
        
        double fJpsiQaMin = JpsiMass - 1.00;            
        double fJpsiQaMax = JpsiMass + 0.75;            
        
        
        // *** electro-magnetic 
        double fEcm = sqrts/10.;
        double f2eQaMin = fEcm - 1.5;           
        double f2eQaMax = fEcm + 1.0;           
        
        
        double low2=0.5;
        // ****** Jpsi -> ee
        TH1F *h_jpsi1=new TH1F("h_jpsi1","J/#psi #rightarrow e^{+}e^{-}",nbins,fJpsiQaMin,fJpsiQaMax);
        TH1F *h_jpsi1sel=new TH1F("h_jpsi1sel","J/psi",nbins,fJpsiQaMin,fJpsiQaMax);
        TH1F *h_jpsi1sig=new TH1F("h_jpsi1sig","J/psi",nbins,fJpsiQaMin,fJpsiQaMax);
        h_jpsi1sig->SetLineColor(2);
        configHisto(h_jpsi1sel);
        
        // ****** Jpsi -> mu mu
        TH1F *h_jpsi2=new TH1F("h_jpsi2","J/#psi #rightarrow #mu^{+}#mu^{-}",nbins,fJpsiQaMin,fJpsiQaMax);
        TH1F *h_jpsi2sel=new TH1F("h_jpsi2sel","J/psi",nbins,fJpsiQaMin,fJpsiQaMax);
        TH1F *h_jpsi2sig=new TH1F("h_jpsi2sig","J/psi",nbins,fJpsiQaMin,fJpsiQaMax);
        h_jpsi2sig->SetLineColor(2);
        configHisto(h_jpsi2sel);
        
        // ****** eta_c -> Ks K pi
        TH1F *h_etac=new TH1F("h_etac","#eta_{c} #rightarrow K_{S} K^{+} #pi^{-}",nbins,fEtacQaMin,fEtacQaMax);
        TH1F *h_etacsel=new TH1F("h_etacsel","etac",nbins,fEtacQaMin,fEtacQaMax);
        TH1F *h_etacsig=new TH1F("h_etacsig","etac",nbins,fEtacQaMin,fEtacQaMax);
        h_etacsig->SetLineColor(2);
        configHisto(h_etacsel);
        
        // ****** D0 -> K- pi+
        TH1F *h_d0=new TH1F("h_d0","D^{0} #rightarrow K^{-}#pi^{+}",nbins,fD0QaMin,fD0QaMax);
        TH1F *h_d0sel=new TH1F("h_d0sel","D0",nbins,fD0QaMin,fD0QaMax);
        TH1F *h_d0sig=new TH1F("h_d0sig","D0",nbins,fD0QaMin,fD0QaMax);;
        h_d0sig->SetLineColor(2);
        configHisto(h_d0sel);
                
        // ****** D+ -> K- pi+ pi+
        TH1F *h_dpm=new TH1F("h_dpm","D^{+} #rightarrow K^{-}#pi^{+}#pi^{+}",nbins,fDpmQaMin,fDpmQaMax);
        TH1F *h_dpmsel=new TH1F("h_dpmsel","D#pm",nbins,fDpmQaMin,fDpmQaMax);
        TH1F *h_dpmsig=new TH1F("h_dpmsig","D#pm",nbins,fDpmQaMin,fDpmQaMax);
        h_dpmsig->SetLineColor(2);
        configHisto(h_dpmsel);
                
        // ****** Ds+ -> K+ K- pi+
        TH1F *h_ds=new TH1F("h_ds","D_{s}^{+} #rightarrow K^{+}K^{-}#pi^{+}",nbins,fDsQaMin,fDsQaMax);
        TH1F *h_dssel=new TH1F("h_dssel","Ds",nbins,fDsQaMin,fDsQaMax);
        TH1F *h_dssig=new TH1F("h_dssig","Ds",nbins,fDsQaMin,fDsQaMax);
        h_dssig->SetLineColor(2);
        configHisto(h_dssel);
        
        // ****** phi -> K+ K-
        TH1F *h_phi=new TH1F("h_phi","#phi #rightarrow K^{+}K^{-}",nbins,fPhiQaMin,fPhiQaMax);
        TH1F *h_phisel=new TH1F("h_phisel","Phi",nbins,fPhiQaMin,fPhiQaMax);
        TH1F *h_phisig=new TH1F("h_phisig","Phi",nbins,fPhiQaMin,fPhiQaMax);
        h_phisig->SetLineColor(2);
        configHisto(h_phisel);
        
        // ****** Lambda_c -> p K- pi+
        TH1F *h_lamc=new TH1F("h_lamc","#Lambda_{c} #rightarrow pK^{-}#pi^{+}",nbins,fLamcQaMin,fLamcQaMax);
        TH1F *h_lamcsel=new TH1F("h_lamcsel","Lam_c",nbins,fLamcQaMin,fLamcQaMax);
        TH1F *h_lamcsig=new TH1F("h_lamcsig","Lam_c",nbins,fLamcQaMin,fLamcQaMax);
        h_lamcsig->SetLineColor(2);
        configHisto(h_lamcsel);
        
        // ****** Lambda -> p pi-
        TH1F *h_lam=new TH1F("h_lam","#Lambda #rightarrow p#pi^{-}",nbins,fLamQaMin,fLamQaMax);
        TH1F *h_lamsel=new TH1F("h_lamsel","Lam",nbins,fLamQaMin,fLamQaMax);
        TH1F *h_lamsig=new TH1F("h_lamsig","Lam",nbins,fLamQaMin,fLamQaMax);
        h_lamsig->SetLineColor(2);
        configHisto(h_lamsel);
        
        // ****** pp -> e+ e-
        TH1F *h_ee=new TH1F("h_ee","p#bar{p} #rightarrow e^{+}e^{-}",nbins,f2eQaMin,f2eQaMax);
        TH1F *h_eesel=new TH1F("h_eesel","ee",nbins,f2eQaMin,f2eQaMax);
        TH1F *h_eesig=new TH1F("h_eesig","ee",nbins,f2eQaMin,f2eQaMax);
        h_eesig->SetLineColor(2);
        configHisto(h_eesel);
        
        h_jpsi1->SetXTitle("m(e^{+}e^{-}) [GeV/c^{2}]");
        h_jpsi2->SetXTitle("m(#mu^{+}#mu^{-}) [GeV/c^{2}]");
        h_etac->SetXTitle("m(K_{S} K^{+}#pi^{-}) [GeV/c^{2}]");

        h_d0->SetXTitle("m(K^{-}#pi^{+}) [GeV/c^{2}]");

        h_dpm->SetXTitle("m(K^{-}#pi^{+}#pi^{+}) [GeV/c^{2}]");

        h_ds->SetXTitle("m(K^{+}K^{-}#pi^{+}) [GeV/c^{2}]");

        h_phi->SetXTitle("m(K^{+}K^{-}) [GeV/c^{2}]");
        h_lamc->SetXTitle("m(pK^{-}#pi^{+}) [GeV/c^{2}]");
        h_lam->SetXTitle("m(p#pi^{-}) [GeV/c^{2}]");
        h_ee->SetXTitle("m(e^{+}e^{-}) [GeV/c^{2}]");
        
        // *** make histos
        fillHistos(t[0], "phim", cuts["nphi"], "h_phi", "h_phisel", "h_phisig");
        fillHistos(t[1], "lamm", cuts["nlam"], "h_lam", "h_lamsel", "h_lamsig");
        fillHistos(t[2], "eem",  cuts["n2e"], "h_ee",  "h_eesel",  "h_eesig");
        
        if (sqrts>24)
        {
                fillHistos(t[3], "jpsim", cuts["njpsi1"], "h_jpsi1", "h_jpsi1sel", "h_jpsi1sig");
                fillHistos(t[4], "jpsim", cuts["njpsi2"], "h_jpsi2", "h_jpsi2sel", "h_jpsi2sig");
                fillHistos(t[5], "etacm",  cuts["netac2"],  "h_etac",  "h_etacsel",  "h_etacsig");
                fillHistos(t[6], "d0m",    cuts["nd01"],    "h_d0",    "h_d0sel",    "h_d0sig");
                fillHistos(t[7], "dpmm",     cuts["ndpm1"],   "h_dpm",   "h_dpmsel",   "h_dpmsig");
        }
        
        if (sqrts>38)
                fillHistos(t[8], "dsm", cuts["nds1"], "h_ds", "h_dssel", "h_dssig");
        
        if (sqrts>45)
                fillHistos(t[9], "lamcm", cuts["nlamc"], "h_lamc", "h_lamcsel", "h_lamcsig");
        
        // event counts
        
        double channelcount[10]={0,0,0,0,0,0,0,0,0,0};
        double eff[10]={0,0,0,0,0,0,0,0,0,0};
        
        evcnt.clear();
        for (i=0;i<3;++i)
        {
                channelcount[i]=cnt_events(t[i], cuts[ntpname[i]]+"&&tag", true);
                cout <<i<< "  "<<channelcount[i]<<endl;
        }
        
        if (sqrts>24)
        for (i=3;i<8;++i)
        {
                channelcount[i]=cnt_events(t[i], cuts[ntpname[i]]+"&&tag", true);
                cout <<i<< "  "<<channelcount[i]<<endl;
        }
        if (sqrts>38)
        for (i=8;i<9;++i)
        {
                channelcount[i]=cnt_events(t[i], cuts[ntpname[i]]+"&&tag", true);
                cout <<i<< "  "<<channelcount[i]<<endl;
        }
        if (sqrts>45)
        for (i=9;i<10;++i)
        {
                channelcount[i]=cnt_events(t[i], cuts[ntpname[i]]+"&&tag", true);
                cout <<i<< "  "<<channelcount[i]<<endl;
        }
        
        double evCnt = evcnt.size();
        for (i=0;i<10;++i) eff[i]=channelcount[i]/nev*100.;
        
        cout <<"all:"<<evCnt<<endl;
        
        // ** draw the plots
        TLatex latex;
        latex.SetTextSize(0.08);
        char tmp[200];
        double offset = 0;//0.1;
        double posx = 0.24;
        double posy = 0.79;
        
//      // *** plot all stuff and write out info
//      c1->cd(1); h_pi0->Draw(); h_pi0sel->Draw("same"); h_pi0sig->Draw("same");
//      c1->cd(2); h_ks->Draw(); h_kssel->Draw("same"); h_kssig->Draw("same");

        c2->cd(2); h_phi->Draw(); h_phisel->Draw("same");   h_phisig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[0]);
        latex.DrawLatexNDC(posx,posy,tmp);
        
        c2->cd(1); h_ee->Draw(); h_eesel->Draw("same"); h_eesig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[2]);
        latex.DrawLatexNDC(posx,posy,tmp);
        
        c2->cd(9); h_lam->Draw(); h_lamsel->Draw("same"); h_lamsig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[1]);
        latex.DrawLatexNDC(posx,posy,tmp);
        
        c2->cd(3); h_etac->Draw(); h_etacsel->Draw("same"); h_etacsig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[5]);
        latex.DrawLatexNDC(posx,posy,tmp);
        
        c2->cd(4); h_jpsi1->Draw(); h_jpsi1sel->Draw("same"); h_jpsi1sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[3]);
        latex.DrawLatexNDC(posx,posy,tmp);
        
        c2->cd(5); h_jpsi2->Draw(); h_jpsi2sel->Draw("same"); h_jpsi2sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[4]);
        latex.DrawLatexNDC(posx,posy,tmp);
        
        c2->cd(7); h_dpm->Draw();  h_dpmsel->Draw("same");  h_dpmsig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[7]);
        latex.DrawLatexNDC(posx,posy,tmp);
        
        c2->cd(6); h_d0->Draw();   h_d0sel->Draw("same");   h_d0sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[6]);
        latex.DrawLatexNDC(posx,posy,tmp);

        c2->cd(8); h_ds->Draw();   h_dssel->Draw("same");   h_dssig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[8]);
        latex.DrawLatexNDC(posx,posy,tmp);
        
        c2->cd(10); h_lamc->Draw(); h_lamcsel->Draw("same"); h_lamcsig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",eff[9]);
        latex.DrawLatexNDC(posx,posy,tmp);
        
        // total efficiency in first plot (phi)
        c2->cd(5);
        posy = h_jpsi2->GetMaximum()*1.12;
        if (posy==0) posy=1.12;
        sprintf(tmp,"(#epsilon_{t} = %4.1f%%)",double(evCnt)/double(nev)*100.);
        latex.SetTextColor(kMagenta+2);
        latex.DrawLatex(3.0,posy,tmp);
        
        c2->cd();
                
        //TString plot(fname(fname.Last('/')+1,fname.Last('.')-fname.Last('/')-1));
//fname(pos1+1,pos2-1));
    TString plot(fname(fname.Last('/')+1,6));
        if (!tagonly) plot+="_opt_supr";
        //plot="figfull/"+plot+".pdf";
        
        //cout <<plot<<endl;
        
        c2->SaveAs("figfull_v2/full_"+plot+".gif");
        c2->SaveAs("figfull_v2/full_"+plot+".pdf");
        
        std::ofstream ofs;
        
        ofs.open("full_supr1000.txt", std::ofstream::out | std::ofstream::app);
        TString result=TString::Format("%s & %3.1f & %6.2f & %6.2f & %6.2f & %6.2f & %6.2f & %6.2f & %6.2f & %6.2f & %6.2f & %6.2f & %6.2f \\\\\n",
                                                                   mode.Data(), sqrts/10., eff[2], eff[0], eff[5], eff[3], eff[4], eff[6], eff[7], eff[8], eff[1], eff[9], evCnt/nev*100);
        //cout <<result;                                                        
        ofs <<result;
        ofs.close();
  return 0;
}

double getPos

(

TH1 *

h

)

Definition at line 71 of file full_core_ntp.C.

{
        return (h->GetBinCenter(h->GetNbinsX())-h->GetBinCenter(1))*0.05+h->GetBinCenter(1);
}

void init_cuts_1000

(

TString

energy

)

Definition at line 113 of file full_core_ntp.C.

References cuts.

Referenced by full_core_ntp().

{
        if (energy=="24")
        {
                cuts["nphi"]    = "eslnpidk>1&&phipcm>0.5&&phipcm<0.75";
                cuts["nlam"]    = "lampcm>0.38&&lampcm<0.48&&eslnpidp>0&&lamp>0.6";
                cuts["n2e"]     = "eslnpide>1";

        }
        else if (energy=="38")
        {
                cuts["nphi"]    = "eslnpidk>1&&phipcm>1.4";
                cuts["nlam"]    = "eslnpidp>0&&lampcm>1.25&&esthr>0.93&&lampt>0.75";
                cuts["n2e"]     = "eslnpide>1";
                cuts["njpsi1"]  = "eslnpide>1&&jpsipcm<0.9";
                cuts["njpsi2"]  = "eslnpidmu>1&&jpsipcm<0.9";
                cuts["netac2"]  = "eslnpidk>0&&eslnpidpi>2";
                cuts["nd01"]    = "eslnpidk>0&&d0pcm<0.4&&d0d0pt>0.4";
                cuts["ndpm1"]   = "eslnpidk>0&&dpmpcm<0.3";
                
        }
        else if (energy=="45")
        {
                cuts["nphi"]    = "eslnpidk>1&&phipcm>1.55&&esthr>0.85";
                cuts["nlam"]    = "esl1npidp>0&&esthr>0.93&&lampt>1";
                cuts["n2e"]     = "eslnpide>1&&mmiss>-1";
                cuts["njpsi1"]  = "eslnpide>1&&jpsipcm<1.8&&detemcmax>2.5";
                cuts["njpsi2"]  = "eslnpidmu>1&&jpsipcm<1.5";
                cuts["netac2"]  = "eslnpidk>0&&eslnpidpi>2&&etacd1pt>0.35&&etacd0pt>0.35";
                cuts["nd01"]    = "eslnpidk>0&&d0pt>1&&mmiss>1.5&&d0d0pidk>0.25";
                cuts["ndpm1"]   = "eslnpidk>0&&dpmpcm>1.1&&dpmpcm<1.5&&dpmpt>0.9&&dpmd0pidk>0.3";
                cuts["nds1"]    = "eslnpidk>1&&dspcm>1&&dspcm<1.17&&dspt>0.5";
                
        }
        else if (energy=="55")
        {
                cuts["nphi"]    = "eslnpidk>1&&phipcm>1.5&&esthr>0.89";
                cuts["nlam"]    = "eslnpidp>0&&esthr>0.97&&lampcm>2.1&&lampt>1";
                cuts["n2e"]     = "eslnpide>1&&mmiss>-1";
                cuts["njpsi1"]  = "eslnpide>1&&jpsipcm<2.2&&detemcmax>4&&jpsid1tht>0.07";
                cuts["njpsi2"]  = "eslnpidmu>1&&jpsipcm<2";
                cuts["netac2"]  = "eslnpidk>0&&eslnpidpi>2&&esptmax>1.1&&esnpart<15";
                cuts["nd01"]    = "eslnpidk>0&&d0pcm>1.9&&mmiss>1.5&&d0pt>1.6&&d0pcm<2.15";
                cuts["ndpm1"]   = "eslnpidk>0&&dpmpcm>1.95&&dpmpcm<2.15&&dpmpt>0.8&&dpmd0pt>0.4";
                cuts["nds1"]    = "eslnpidk>1&&dspcm>1.85&&dspcm<2.1&&dspt>1.3";
                cuts["nlamc"]   = "eslnpidk>0&&eslnpidp>0&&lamcpt>.9&&lamcd0pidp>0.21&&lamcd0thtcm<2";
                
        }
        else
        {
                cuts["nphi"]    = "1";
                cuts["nlam"]    = "1";
                cuts["n2e"]     = "1";
                cuts["njpsi1"]  = "1";
                cuts["njpsi2"]  = "1";
                cuts["netac2"]  = "1";
                cuts["nd01"]    = "1";
                cuts["ndpm1"]   = "1";
                cuts["nds1"]    = "1";
                cuts["nlamc"]   = "1";
        }
                // factor 1000
/*              cuts["nphi"]    = "fw1>0.0129387&&phipt>0.2";
                cuts["nlam"]    = "fw3>0.2&&fw4>0.3&&sumpc>0.9";
                cuts["nee"]     = "1";*/

                // factor 1000
/*              cuts["nphi"]    = "phipcm>1.31489&&phipt>0.698733";
                cuts["nlam"]    = "fw2>0.688083&&lampt>0.990166";
                cuts["nee"]     = "1";
                cuts["njpsi1"]  = "1";
                cuts["njpsi2"]  = "1";
                cuts["netac"]   = "ptmax>0.84&&sumptc>2.9";
                cuts["nd0"]     = "d0pcm<0.35&&ptmax>0.8";
                cuts["ndpm"]    = "dpcm<0.34&&ptmax>0.67&&fw1>0";*/
                
                // first attempt
/*              cuts["netac"]   = "ptmax>0.7&&sumptc>2.7";
                cuts["nd0"]     = "d0pcm<0.33&&ptmax>0.75";
                cuts["ndpm"]    = "dpcm<0.31&&ptmax>0.6&&pmin>0.1";*/
                // factor 1000
//              cuts["nphi"]    = "phipcm>1.9&&phipt>1";
//              cuts["nlam"]    = "lampcm>1.8&&thr>0.95&&lampt>1";
//              cuts["nee"]     = "1";
//              cuts["njpsi1"]  = "sumpc>3.5";
//              cuts["njpsi2"]  = "sumpc>3.5";
//              cuts["netac"]   = "ptmax>1&&sumpc>4";
//              cuts["nd0"]     = "abs(d0pcm-1.24)<0.15&&ptmax>1.05&&sumptc>3.1";
//              cuts["ndpm"]    = "abs(dpcm-1.24)<0.12&&dpt>1&&sumptc>3";
//              cuts["nds"]     = "abs(dspcm-1.1)<0.09&&ptmax>0.93";

                // first attempt                                cuts["nphi"]    = "phipcm>1.9&&phipt>1";
/*              cuts["nd0"]     = "abs(d0pcm-1.24)<0.15&&ptmax>1.2";
                cuts["ndpm"]    = "abs(dpcm-1.24)<0.15&&dpt>1&&pmin>0.15";
                cuts["nds"]     = "abs(dspcm-1.1)<0.1&&ptmax>0.8";*/
                // factor 1000
/*              cuts["nphi"]    = "phipcm>2.4&&phipt>0.5";
                cuts["nlam"]    = "lampcm>0.23&&thr>0.98&&lampt>1.2";
                cuts["nee"]     = "1";
                cuts["njpsi1"]  = "sumpc>4";
                cuts["njpsi2"]  = "sumpc>4";
                cuts["netac"]   = "ptmax>1.1&&sumpc>4.8";
                cuts["nd0"]     = "ptmax>1&&d0pcm>1.8&&d0pt>1.5";
                cuts["ndpm"]    = "ptmax>1&&dpcm>1.8&&dpt>1.5";
                cuts["nds"]     = "dspcm>1.8&&ptmax>0.8&&dspt>1.2";
                cuts["nlamc"]   = "abs(lamcpcm-1.53)<0.12&&lamcpt>1.3&&fw1>0";*/
                // first attempt
/*              cuts["nlamc"]   = "abs(lamcpcm-1.53)<0.15&&lamcpt>1.3";*/
}

void init_cuts_95

(

TString

energy

)

Definition at line 222 of file full_core_ntp.C.

References cuts.

{
        if (energy=="24")
        {
                cuts["nphi"]    = "eslnpidk>1&&phipcm>0.45&&phipcm<0.78";
                cuts["nlam"]    = "eslnpidp>0&&lampcm<0.55&&lampcm>0.25&&mmiss<1.2";
                cuts["n2e"]     = "eslnpide>1&&mmiss>-1";

        }
        else if (energy=="38")
        {
                cuts["nphi"]    = "eslnpidk>1&&phipcm>1.2";
                cuts["nlam"]    = "esthr>0.85&&lampcm>1.15&&lampt>0.2&&eslnpidp>0";
                cuts["n2e"]     = "eslnpide>1&&mmiss>-1";
                cuts["njpsi1"]  = "eslnpide>1&&jpsipcm<1";
                cuts["njpsi2"]  = "eslnpidmu>1&&jpsipcm<1";
                cuts["netac2"]  = "etacpcm<1.25&&essumptc>1.75";
                cuts["nd01"]    = "eslnpidk>0&&d0pcm<0.4&&d0d0pt>0.25";
                cuts["ndpm1"]   = "eslnpidk>0";
                
        }
        else if (energy=="45")
        {
                cuts["nphi"]    = "eslnpidk>1&&phipcm>1.55";
                cuts["nlam"]    = "esthr>0.87&&eslnpidp>0&&lampcm>1.5";
                cuts["n2e"]     = "eslnpide>1&&mmiss>-1";
                cuts["njpsi1"]  = "eslnpide>1&&jpsipcm<1.8&&detemcmax>2.5";
                cuts["njpsi2"]  = "eslnpidmu>1&&jpsipcm<1.5";
                cuts["netac2"]  = "etacpcm<1.25&&etacd1p>0.25";
                cuts["nd01"]    = "eslnpidk>0&&d0pcm<1.6&&d0pcm>0.8&&d0d1p<6";
                cuts["ndpm1"]   = "eslnpidk>0";
                cuts["nds1"]    = "eslnpidk>0&&dspcm>0.4&&dsp<8&&dsthtcm>0.5";
                
        }
        else if (energy=="55")
        {
                cuts["nphi"]    = "eslnpidk>1&&phipcm>1.75&&phipt>0.5";
                cuts["nlam"]    = "esthr>0.9&&lampcm>1&&mmiss<2.4";
                cuts["n2e"]     = "eslnpide>0&&mmiss>-1";
                cuts["njpsi1"]  = "eslnpide>1&&jpsipcm<1.9&&detemcmax>3.5";
                cuts["njpsi2"]  = "eslnpidmu>1&&jpsipcm<2.5";
                cuts["netac2"]  = "etacpcm<2&&etacd1p>0.4&&etacecm>2.85&&eslnpide<1";
                cuts["nd01"]    = "eslnpidk>0&&d0pcm<2.3&&d0pcm>0.95";
                cuts["ndpm1"]   = "eslnpidk>0";
                cuts["nds1"]    = "eslnpidk>0&&dspcm>0.5&&dsthtcm>0.5";
                cuts["nlamc"]   = "mmiss>2";
                
        }
        else
        {
                cuts["nphi"]    = "1";
                cuts["nlam"]    = "1";
                cuts["n2e"]     = "1";
                cuts["njpsi1"]  = "1";
                cuts["njpsi2"]  = "1";
                cuts["netac2"]  = "1";
                cuts["nd01"]    = "1";
                cuts["ndpm1"]   = "1";
                cuts["nds1"]    = "1";
                cuts["nlamc"]   = "1";
        }
}

void init_modes

(

)

Definition at line 86 of file full_core_ntp.C.

References mctvar, and modes.

Referenced by full_core_ntp().

{
        modes["nphi"]   = "000";        
        modes["nlam"]   = "400";
        modes["n2e"]    = "300";
        modes["njpsi1"] = "200";
        modes["njpsi2"] = "201";
        modes["netac2"] = "221";        
        modes["nd01"]   = "100";
        modes["ndpm1"]  = "120";        
        modes["nds1"]   = "140";        
        modes["nlamc"]  = "420";

        modes["DPM"]    = "900";
        
        mctvar["nphi"]  = "phi";        
        mctvar["nlam"]  = "lam";
        mctvar["n2e"]   = "ee";
        mctvar["njpsi1"]= "jpsi";
        mctvar["njpsi2"]= "jpsi";
        mctvar["netac2"]= "etac";       
        mctvar["nd01"]  = "d0";
        mctvar["ndpm1"] = "dpm";        
        mctvar["nds1"]  = "ds"; 
        mctvar["nlamc"] = "lamc";
}

void setStyle

(

)

Definition at line 47 of file full_core_ntp.C.

Referenced by full_core_ntp().

{
        gStyle->SetPadTopMargin(0.13);
        gStyle->SetPadBottomMargin(0.18);
        gStyle->SetPadLeftMargin(0.2);
        gStyle->SetPadRightMargin(0.016);

        gStyle->SetLabelSize(0.076,"X");
        gStyle->SetLabelSize(0.07,"Y");

        gStyle->SetStatY(0.87);
        gStyle->SetStatX(0.986);
        gStyle->SetStatW(0.4);
        gStyle->SetStatH(0.12);
        gStyle->SetOptStat("e"); // only entries
        
        gStyle->SetTitleH(0.1);
        gStyle->SetTitleX(0.08);
        gStyle->SetTitleOffset(1.1,"x");
        gStyle->SetTitleXSize(0.076);
        
        gStyle->SetNdivisions(505);
}

Variable Documentation

std::map<TString, TString> cuts

Definition at line 39 of file full_core_ntp.C.

Referenced by full_core_ntp(), init_cuts_1000(), and init_cuts_95().

std::map<int, int> evcnt

Definition at line 37 of file full_core_ntp.C.

Referenced by cnt_events(), and full_core_ntp().

TLorentzVector fIni

Definition at line 29 of file full_core_ntp.C.

Referenced by prod_fsim(), prod_sim(), QAmacro_mvd_1(), quickfsimana(), run_sim(), PndMasterRunSim::Setup(), sim(), sim_complete(), sim_complete_pbarA(), sim_complete_runs(), sim_evtgen(), sim_rich(), sim_scit(), sim_test(), simfast_opt(), tb_sim_complete(), tb_sim_scit(), and tut_fastsim().

const int fMAX = 100

Definition at line 32 of file full_core_ntp.C.

TRandom3 fRand

Definition at line 28 of file full_core_ntp.C.

double high = 6.

Definition at line 35 of file full_core_ntp.C.

Referenced by plotTrackCands().

double low = 0.0

Definition at line 34 of file full_core_ntp.C.

std::map<TString, TString> mctvar

Definition at line 41 of file full_core_ntp.C.

Referenced by fillHistos(), and init_modes().

std::map<TString, TString> modes

Definition at line 40 of file full_core_ntp.C.

Referenced by init_modes().

int nbins = 300

Definition at line 33 of file full_core_ntp.C.

Referenced by anaclust(), PndHypGeSpectrumAnalyser::DrawCalibratedSpectrum(), fillUserHisto(), full_core_ntp(), PndDrcRecoLookupMap::ProcessPhotonHit(), and TransformHisto().

TString ntpname[10] = {"nphi","nlam","n2e","njpsi1","njpsi2","netac2","nd01","ndpm1","nds1","nlamc"}

Definition at line 43 of file full_core_ntp.C.

Referenced by full_core_ntp().

TString plotname[10] = {"phim","lamm","eem","jpsim","jpsim","etacm","d0m","dpmm","dsm","lamcm"}

Definition at line 44 of file full_core_ntp.C.

Referenced by toy_core().

std::map<int, int> treeidx

Definition at line 38 of file full_core_ntp.C.

Referenced by cnt_events().