MakePlots Class Reference

#include <MakePlots.h>

List of all members.

Public Member Functions
	MakePlots ()
	~MakePlots ()
void	FillHistos (int set)
void	SetupHistos ()
void	Run ()
void	SavePlots (int norm=0)
void	SaveRatioPlots ()
void	SetRatioPair (int set1, int set2, const char *name)
void	AddSet (const char *filename, const char *typesname, int ismc, int ismrcc=0)
void	SetOutDir (const char *o)
void	SetOutFile (const char *o)
void	SetPots (int chain, double p)
Public Attributes
std::vector< double >	pots
Private Member Functions
void	SortOutStats (TPad *pad, float fracPadx=0.3, float fracPady=0.5, int numPerCol=4, float startx=0.98, float starty=0.98)
void	SetBranches (TChain *c)
void	FillHistos (int type, int cut)
Private Attributes
std::vector< std::vector< std::vector< TH1F * > > >	histos
std::vector< std::pair< int, int > >	ratiopair
std::vector< std::string >	rationame
char	names [nHistos][200]
std::vector< int >	ismc
std::vector< int >	ismrcc
char	cutnames [nCut][200]
std::vector< TChain * >	chains
std::vector< std::string >	typesname
int	nchains
char	outdir [200]
char	outfile [200]
Double_t	particles_primary_long_e
Double_t	particles_longest_s_particle_s
Double_t	particles_elec_vise
Double_t	particles_primary_phi
Double_t	particles_prim_par_b
Double_t	particles_mol_rad_r
Double_t	particles_longest_z
Double_t	particles_emfrac
Int_t	particles_ntot
Double_t	particles_total_long_e
Double_t	particles_weighted_phi
Double_t	particles_largest_particle_e
Double_t	particles_totvise
Double_t	particles_prim_vise
Double_t	event_visenergy
Double_t	event_max_z
Double_t	event_min_z
Int_t	event_nstrips
Int_t	event_inFiducial
Int_t	event_contained
Double_t	event_pidA
Double_t	event_pidB
Double_t	event_pidC
Double_t	event_pidD
Int_t	mctrue_type
Double_t	mctrue_oscprob
Double_t	mctrue_totbeamweight
Int_t	mctrue_iresonance
Double_t	trainweight
Int_t	type
Double_t	largest_frac
Double_t	reco_frac
Int_t	isnc
Int_t	isnue
Int_t	iscc
Int_t	istau
Int_t	isbeamve
Int_t	isndis
Int_t	isdis
Double_t	particles_rms_r
Double_t	particles_prim_par_a
Double_t	particles_prim_par_e0
Double_t	particles_prim_par_chisq
Double_t	particles_largest_particle_peakdiff
Double_t	particles_largest_particle_cmp_chisq
Int_t	particles_largest_particlecmp_ndf
Double_t	particles_longest_particle_type
Double_t	mrccinfo_particle_s
Double_t	mrccinfo_sum_e

---

Constructor & Destructor Documentation

MakePlots::MakePlots (   )

Definition at line 14 of file MakePlots.cxx. References nchains, outdir, and outfile. { nchains=0; sprintf(outdir,"plots"); sprintf(outfile,"plots.root"); }

MakePlots::~MakePlots (   )  [inline]

Definition at line 19 of file MakePlots.h. {};

---

Member Function Documentation

void MakePlots::AddSet (  const char *  filename, const char *  typesname, int  ismc, int  ismrcc = 0 )

Definition at line 21 of file MakePlots.cxx. References chains, histos, ismc, ismrcc, nchains, pots, and typesname. { nchains=chains.size(); chains.push_back(new TChain("PA")); int added = chains[nchains]->Add(filename); if(!added) { printf("error adding chain %s '%s' mc %d\n",filename,tn, mc); delete chains[nchains]; return; } printf("adding chain %s '%s' mc %d\n",filename,tn, mc); char tmpchar[200]; sprintf(tmpchar,"%s",tn); typesname.push_back(tmpchar); ismc.push_back(mc); ismrcc.push_back(mrcc); TChain * POTchain=new TChain("pottree"); POTchain->Add(filename); int ent=POTchain->GetEntries(); pots.push_back(0); for(int i=0;i<ent;i++) { POTchain->GetEntry(i); pots[nchains]+=POTchain->GetLeaf("pot")->GetValue(); } printf("%s has %f POTs \n",tn,pots[nchains]); nchains++; std::vector<std::vector<TH1F *> > a; histos.push_back(a); }

void MakePlots::FillHistos (  int  type, int  cut )  [private]

Definition at line 656 of file MakePlots.cxx. References event_nstrips, event_visenergy, histos, ismc, mctrue_totbeamweight, particles_elec_vise, particles_emfrac, particles_largest_particle_cmp_chisq, particles_largest_particle_e, particles_largest_particle_peakdiff, particles_largest_particlecmp_ndf, particles_longest_s_particle_s, particles_longest_z, particles_mol_rad_r, particles_ntot, particles_prim_par_a, particles_prim_par_b, particles_prim_par_chisq, particles_prim_par_e0, particles_prim_vise, particles_primary_long_e, particles_primary_phi, particles_rms_r, particles_total_long_e, particles_totvise, particles_weighted_phi, and pots. { if(!pots[type])return; double weight = 1; if(ismc[type]==1)weight*=mctrue_totbeamweight*mctrue_oscprob; histos[type][cut][0]->Fill(particles_longest_z,weight); histos[type][cut][1]->Fill(particles_longest_s_particle_s,weight); histos[type][cut][2]->Fill(particles_primary_long_e,weight); histos[type][cut][3]->Fill(particles_elec_vise/25.,weight); histos[type][cut][4]->Fill(particles_primary_phi,weight); if(particles_prim_par_e0)histos[type][cut][5]->Fill(particles_prim_par_b,weight); histos[type][cut][6]->Fill(particles_mol_rad_r,weight); histos[type][cut][7]->Fill(particles_emfrac,weight); histos[type][cut][8]->Fill(particles_ntot,weight); histos[type][cut][9]->Fill(particles_total_long_e,weight); histos[type][cut][10]->Fill(particles_weighted_phi,weight); histos[type][cut][11]->Fill((particles_largest_particle_e*0.03379042+0.5866994)/particles_totvise,weight); histos[type][cut][12]->Fill(particles_prim_vise/25./event_visenergy,weight); histos[type][cut][13]->Fill(particles_totvise,weight); histos[type][cut][14]->Fill(event_visenergy,weight); histos[type][cut][15]->Fill(particles_rms_r,weight); if(particles_prim_par_e0)histos[type][cut][16]->Fill(particles_prim_par_a,weight); if(particles_prim_par_e0)histos[type][cut][17]->Fill(particles_prim_par_e0/25./60.,weight); histos[type][cut][18]->Fill(particles_prim_par_chisq,weight); histos[type][cut][19]->Fill(particles_largest_particle_peakdiff/25.,weight); histos[type][cut][20]->Fill(particles_largest_particle_cmp_chisq/particles_largest_particlecmp_ndf,weight); histos[type][cut][21]->Fill(event_nstrips,weight); }

void MakePlots::FillHistos (  int  set  )

Definition at line 587 of file MakePlots.cxx. References chains, event_contained, event_inFiducial, event_pidB, event_visenergy, histos, ismrcc, mrccinfo_particle_s, mrccinfo_sum_e, particles_largest_particle_cmp_chisq, particles_longest_particle_type, particles_longest_s_particle_s, particles_ntot, particles_totvise, pots, and SetBranches(). Referenced by Run(). { TChain *c=chains[set]; printf("make plots set %d\n",set); SetBranches(c); if(!pots[set])return; int ent=c->GetEntries(); printf("%d entries\n",ent); for(int i=0;i<ent;i++) { // if(i>1000)break;//testing if(i%10000==0)printf("entry %d\n",i); c->GetEntry(i); //is this an mrcc file? if so do the mrcc cut if(ismrcc[set]) { if(mrccinfo_particle_s<3)continue;//no long muon track event_visenergy-= (mrccinfo_sum_e*0.03799819+0.4803261); //adjust! } FillHistos(set,0); if(event_inFiducial!=1 || event_contained!=1)continue; FillHistos(set,1); //presel if(particles_ntot<1)continue; if(particles_longest_s_particle_s>1.2 ||particles_longest_s_particle_s<0.2)continue; //broken in mrcc files if(event_max_z-event_min_z>0.9 || event_max_z-event_min_z < 0.2)continue; //need to use particle total energy rather than strip energy for mrcc... if(particles_totvise<0.2 || particles_totvise>8)continue; //if(event_nstrips<10 || event_nstrips>40)continue; //don't want a long muon.. if(particles_longest_particle_type>12)continue; if(particles_largest_particle_cmp_chisq>30)continue; FillHistos(set,2); if(event_pidB<0.5)continue; // if(event_pidA<0.5)continue; FillHistos(set,3); } for(int j=0;j<nCut;j++) for(int k=0;k<nHistos;k++) { if(pots[set])histos[set][j][k]->Scale(3.25e8/pots[set]); } }

void MakePlots::Run (   )

Definition at line 56 of file MakePlots.cxx. References cutnames, FillHistos(), outdir, and SetupHistos(). { // AddSet("/minos/data/dogwood0/near/data/PO*.root","near data",0); // AddSet("/minos/data/dogwood0/near/mc/PO*.root","near mc",1); /* TChain * near_mc=new TChain("PA"); TChain * near_data=new TChain("PA"); TChain * near_mc_POT=new TChain("pottree"); TChain * near_data_POT=new TChain("pottree"); / *near_data->Add("/minos/data/dogwood0/near/data/PO*.root"); near_data_POT->Add("/minos/data/dogwood0/near/data/PO*.root"); near_mc->Add("/minos/data/dogwood0/near/mc/PO*.root"); near_mc_POT->Add("/minos/data/dogwood0/near/mc/PO*.root"); */ /* near_data->Add("/data/minos/ManagedPIDCVS/PO-N00011449_0000.spill.sntp.dogwood0.0-Trim-NN.root"); near_data_POT->Add("/data/minos/ManagedPIDCVS/PO-N00011449_0000.spill.sntp.dogwood0.0-Trim-NN.root"); near_mc->Add("/data/minos/ManagedPIDCVS/PO-n13037040_0000_L010185N_D04.sntp.dogwood0-Trim-NN.root"); near_mc_POT->Add("/data/minos/ManagedPIDCVS/PO-n13037040_0000_L010185N_D04.sntp.dogwood0-Trim-NN.root"); for(int i=0;i<nType;i++)pots[i]=0; int ent; ent=near_data_POT->GetEntries(); for(int i=0;i<ent;i++) { near_data_POT->GetEntry(i); pots[0]+=near_data_POT->GetLeaf("pot")->GetValue(); } printf("%f near data POTs\n",pots[0]); ismc[0]=0; ent=near_mc_POT->GetEntries(); for(int i=0;i<ent;i++) { near_mc_POT->GetEntry(i); pots[1]+=near_mc_POT->GetLeaf("pot")->GetValue(); } printf("%f near mc POTs\n",pots[1]); ismc[1]=1; */ if(!nchains) { printf("You must specify some input sets...\n"); return; } printf("outdir: %s\n",outdir); //define cut names sprintf(cutnames[0],"No Cuts"); sprintf(cutnames[1],"Fiducial"); sprintf(cutnames[2],"Preselection"); sprintf(cutnames[3],"PID"); SetupHistos(); for(int i=0;i<nchains;i++) FillHistos(i); }

void MakePlots::SavePlots (  int  norm = 0  )

Definition at line 466 of file MakePlots.cxx. References cutnames, histos, max, min, names, outdir, outfile, SortOutStats(), and typesname. { TFile *of = TFile::Open(outfile,"RECREATE"); of->cd(); for(int k=0;k<nHistos;k++) for(int j=0;j<nCut;j++) for(int i=0;i<(int)histos.size();i++) { if(norm && histos[i][j][k]->Integral())histos[i][j][k]->Scale(1./histos[i][j][k]->Integral()); histos[i][j][k]->Write(); } of->Close(); TCanvas c1; TPad p1("p1","",0,0.4,1,1); TPad p2("p2","",0,0,1,0.4); c1.Draw(); p1.Draw(); p2.Draw(); for(int k=0;k<nHistos;k++) for(int j=0;j<nCut;j++) { p1.cd(); double maxb=0; for(int i=0;i<(int)histos.size();i++) { char tmpchar[200]; sprintf(tmpchar,"%s %s",names[k],cutnames[j]); histos[i][j][k]->SetTitle(tmpchar); double maxt=histos[i][j][k]->GetMaximum(); histos[i][j][k]->SetXTitle(names[k]); if(maxt>maxb)maxb=maxt; } for(int i=0;i<(int)histos.size();i++) { histos[i][j][k]->GetYaxis()->SetRangeUser(0,maxb*1.2); histos[i][j][k]->SetName(typesname[i].c_str()); if(i==0) { histos[i][j][k]->Draw("E"); continue; } histos[i][j][k]->SetLineColor(1+i); histos[i][j][k]->Draw("SAMESE"); } p2.cd(); double extt=0; std::vector<TH1F *> tmp; for(int i=1;i<(int)histos.size();i++) { tmp.push_back((TH1F*)histos[0][j][k]->Clone("tmp")); tmp[i-1]->Add(histos[i][j][k],-1); tmp[i-1]->Divide(histos[i][j][k]); tmp[i-1]->SetLineColor(i+1); int minbin=tmp[i-1]->GetMinimumBin(); int maxbin=tmp[i-1]->GetMaximumBin(); double min=tmp[i-1]->GetBinContent(minbin); double max=tmp[i-1]->GetBinContent(maxbin); double ext = TMath::Abs(min); ext = TMath::Abs(max) < ext?ext:TMath::Abs(max); if(ext>extt)extt=ext; } extt=extt>0.2?0.2:extt; extt*=1.25; for(int i=0;i<(int)histos.size()-1;i++) { tmp[i]->GetYaxis()->SetRangeUser(-extt,extt); tmp[i]->SetStats(0); tmp[i]->SetTitle(""); if(i==0) { tmp[i]->Draw(); continue; } tmp[i]->Draw("same"); } TLine l(tmp[0]->GetXaxis()->GetXmin(),0,tmp[0]->GetXaxis()->GetXmax(),0); l.Draw("same"); SortOutStats(&c1); char tmpchar[200]; printf("%d %s/%s_%d.gif\n",k,outdir,names[k],j); sprintf(tmpchar,"%s/%s_%d.gif",outdir,names[k],j); c1.SaveAs(tmpchar); for(int i=0;i<(int)histos.size()-1;i++)delete tmp[i]; } }

void MakePlots::SaveRatioPlots (   )

Definition at line 341 of file MakePlots.cxx. References Nav::GetName(), histos, max, min, names, outdir, outfile, rationame, ratiopair, and SortOutStats(). { if(ratiopair.size()<1)return; const int nHist=ratiopair.size(); std::vector<TH1F*> histos[nCut][nHistos]; for(int j=0;j<nCut;j++) for(int k=0;k<nHistos;k++) for(int i=0;i<nHist;i++) { histos[j][k].push_back((TH1F*)0); printf("making %s with %s and %s %d %d\n",rationame[i].c_str(),this->histos[ratiopair[i].first][j][k]->GetName(),this->histos[ratiopair[i].second][j][k]->GetName(),j,k); histos[j][k][i]=(TH1F*)this->histos[ratiopair[i].first][j][k]->Clone(rationame[i].c_str()); histos[j][k][i]->SetDirectory(0); histos[j][k][i]->Divide((TH1F*)this->histos[ratiopair[i].second][j][k]); printf("made %s max %f\n",histos[j][k][i]->GetName(),histos[j][k][i]->GetMaximum()); } TFile *of = TFile::Open(outfile,"RECREATE"); of->cd(); for(int k=0;k<nHistos;k++) for(int j=0;j<nCut;j++) for(int i=0;i<nHist;i++) histos[j][k][i]->Write(); of->Close(); TCanvas c1; TPad p1("p1","",0,0.4,1,1); TPad p2("p2","",0,0,1,0.4); c1.Draw(); p1.Draw(); p2.Draw(); for(int k=0;k<nHistos;k++) for(int j=0;j<nCut;j++) { p1.cd(); double maxb=0; for(int i=0;i<nHist;i++) { //char tmpchar[200]; double maxt=histos[j][k][i]->GetMaximum(); if(maxt>maxb)maxb=maxt; } for(int i=0;i<nHist;i++) { // histos[j][k][i]->GetYaxis()->SetRangeUser(0.5,1.5);//0,maxb*1.2); if(i==0) { histos[j][k][i]->Draw("E"); continue; } histos[j][k][i]->SetLineColor(1+i); histos[j][k][i]->Draw("SAMESE"); } p2.cd(); double extt=0; std::vector<TH1F*> tmp; if(nHist>1) { for(int i=1;i<nHist;i++) { tmp.push_back((TH1F*)0); tmp[i-1]= (TH1F*)histos[j][k][0]->Clone("tmp"); tmp[i-1]->Add(histos[j][k][i],-1); tmp[i-1]->Divide(histos[j][k][i]); tmp[i-1]->SetLineColor(i+1); int minbin=tmp[i-1]->GetMinimumBin(); int maxbin=tmp[i-1]->GetMaximumBin(); double min=tmp[i-1]->GetBinContent(minbin); double max=tmp[i-1]->GetBinContent(maxbin); double ext = TMath::Abs(min); ext = TMath::Abs(max) < ext?ext:TMath::Abs(max); if(ext>extt)extt=ext; } extt=extt>0.2?0.2:extt; extt*=1.25; for(int i=0;i<nHist-1;i++) { tmp[i]->GetYaxis()->SetRangeUser(-extt,extt); tmp[i]->SetStats(0); tmp[i]->SetTitle(""); if(i==0) { tmp[i]->Draw(); continue; } tmp[i]->Draw("same"); } TLine l(tmp[0]->GetXaxis()->GetXmin(),0,tmp[0]->GetXaxis()->GetXmax(),0); l.Draw("same"); } SortOutStats(&c1); char tmpchar[200]; printf("%d %s/%s_%d.gif\n",k,outdir,names[k],j); sprintf(tmpchar,"%s/%s_%d.gif",outdir,names[k],j); c1.SaveAs(tmpchar); for(int i=0;i<(int)tmp.size()-1;i++){delete tmp[i];tmp[i]=0;} } }

void MakePlots::SetBranches (  TChain *  c  )  [private]

Definition at line 253 of file MakePlots.cxx. References event_contained, event_inFiducial, event_max_z, event_min_z, event_nstrips, event_pidA, event_pidB, event_pidC, event_pidD, event_visenergy, mctrue_iresonance, mctrue_oscprob, mctrue_totbeamweight, mctrue_type, mrccinfo_particle_s, mrccinfo_sum_e, particles_elec_vise, particles_emfrac, particles_largest_particle_cmp_chisq, particles_largest_particle_e, particles_largest_particle_peakdiff, particles_largest_particlecmp_ndf, particles_longest_particle_type, particles_longest_s_particle_s, particles_longest_z, particles_mol_rad_r, particles_ntot, particles_prim_par_a, particles_prim_par_b, particles_prim_par_chisq, particles_prim_par_e0, particles_prim_vise, particles_primary_long_e, particles_primary_phi, particles_rms_r, particles_total_long_e, particles_totvise, and particles_weighted_phi. Referenced by FillHistos(). { c->SetMakeClass(1); c->SetBranchStatus("*",0); c->SetBranchStatus("particles.primary_long_e",1); c->SetBranchStatus("particles.longest_s_particle_s",1); c->SetBranchStatus("particles.elec_vise",1); c->SetBranchStatus("particles.primary_phi",1); c->SetBranchStatus("particles.prim_par_b",1); c->SetBranchStatus("particles.mol_rad_r",1); c->SetBranchStatus("particles.longest_z",1); c->SetBranchStatus("particles.emfrac",1); c->SetBranchStatus("particles.ntot",1); c->SetBranchStatus("particles.total_long_e",1); c->SetBranchStatus("particles.weighted_phi",1); c->SetBranchStatus("particles.largest_particle_e",1); c->SetBranchStatus("particles.totvise",1); c->SetBranchStatus("particles.prim_vise",1); c->SetBranchStatus("event.visenergy",1); c->SetBranchStatus("event.nstrips",1); c->SetBranchStatus("event.inFiducial",1); c->SetBranchStatus("event.max_z",1); c->SetBranchStatus("event.min_z",1); c->SetBranchStatus("mctrue.type",1); c->SetBranchStatus("mctrue.oscprob",1); c->SetBranchStatus("mctrue.totbeamweight",1); c->SetBranchStatus("mctrue.iresonance",1); c->SetBranchStatus("event.contained",1); c->SetBranchStatus("event.pidA",1); c->SetBranchStatus("event.pidB",1); c->SetBranchStatus("event.pidC",1); c->SetBranchStatus("event.pidD",1); c->SetBranchStatus("mrccinfo.particle_s",1); c->SetBranchStatus("mrccinfo.sum_e",1); c->SetBranchStatus("particles.rms_r",1); c->SetBranchStatus("particles.prim_par_a",1); c->SetBranchStatus("particles.prim_par_e0",1); c->SetBranchStatus("particles.prim_par_chisq",1); c->SetBranchStatus("particles.largest_particle_peakdiff",1); c->SetBranchStatus("particles.largest_particle_cmp_chisq",1); c->SetBranchStatus("particles.largest_particle_cmp_ndf",1); c->SetBranchStatus("particles.longest_particle_type",1); c->SetBranchAddress("particles.primary_long_e",&particles_primary_long_e); c->SetBranchAddress("particles.longest_s_particle_s",&particles_longest_s_particle_s); c->SetBranchAddress("particles.elec_vise",&particles_elec_vise); c->SetBranchAddress("particles.primary_phi",&particles_primary_phi); c->SetBranchAddress("particles.prim_par_b",&particles_prim_par_b); c->SetBranchAddress("particles.mol_rad_r",&particles_mol_rad_r); c->SetBranchAddress("particles.longest_z",&particles_longest_z); c->SetBranchAddress("particles.emfrac",&particles_emfrac); c->SetBranchAddress("particles.ntot",&particles_ntot); c->SetBranchAddress("particles.total_long_e",&particles_total_long_e); c->SetBranchAddress("particles.weighted_phi",&particles_weighted_phi); c->SetBranchAddress("particles.largest_particle_e",&particles_largest_particle_e); c->SetBranchAddress("particles.totvise",&particles_totvise); c->SetBranchAddress("particles.prim_vise",&particles_prim_vise); c->SetBranchAddress("event.visenergy",&event_visenergy); c->SetBranchAddress("event.nstrips",&event_nstrips); c->SetBranchAddress("event.inFiducial",&event_inFiducial); c->SetBranchAddress("event.max_z",&event_max_z); c->SetBranchAddress("event.min_z",&event_min_z); c->SetBranchAddress("mctrue.type",&mctrue_type); c->SetBranchAddress("mctrue.oscprob",&mctrue_oscprob); c->SetBranchAddress("mctrue.totbeamweight",&mctrue_totbeamweight); c->SetBranchAddress("mctrue.iresonance",&mctrue_iresonance); c->SetBranchAddress("event.contained",&event_contained); c->SetBranchAddress("event.pidA",&event_pidA); c->SetBranchAddress("event.pidB",&event_pidB); c->SetBranchAddress("event.pidC",&event_pidC); c->SetBranchAddress("event.pidD",&event_pidD); c->SetBranchAddress("mrccinfo.particle_s",&mrccinfo_particle_s); c->SetBranchAddress("mrccinfo.sum_e",&mrccinfo_sum_e); c->SetBranchAddress("particles.rms_r",&particles_rms_r); c->SetBranchAddress("particles.prim_par_a",&particles_prim_par_a); c->SetBranchAddress("particles.prim_par_e0",&particles_prim_par_e0); c->SetBranchAddress("particles.prim_par_chisq",&particles_prim_par_chisq); c->SetBranchAddress("particles.largest_particle_peakdiff",&particles_largest_particle_peakdiff); c->SetBranchAddress("particles.largest_particle_cmp_chisq",&particles_largest_particle_cmp_chisq); c->SetBranchAddress("particles.largest_particle_cmp_ndf",&particles_largest_particlecmp_ndf); c->SetBranchAddress("particles.longest_particle_type",&particles_longest_particle_type); }

void MakePlots::SetOutDir (  const char *  o  )

Definition at line 131 of file MakePlots.cxx. References outdir. { sprintf(outdir,"%s",o); }

void MakePlots::SetOutFile (  const char *  o  )

Definition at line 138 of file MakePlots.cxx. References outfile. { sprintf(outfile,"%s",o); }

void MakePlots::SetPots (  int  chain, double  p )  [inline]

Definition at line 46 of file MakePlots.h. {if(chain>-1 && chain<nchains)pots[chain]=p;};

void MakePlots::SetRatioPair (  int  set1, int  set2, const char *  name )  [inline]

Definition at line 35 of file MakePlots.h. { ratiopair.push_back(std::pair<int,int>(set1,set2)); rationame.push_back(name); };

void MakePlots::SetupHistos (   )

Definition at line 145 of file MakePlots.cxx. References histos, and names. Referenced by Run(). { char tmp[200]; for(int i=0;i<(int)histos.size();i++) for(int j=0;j<nCut;j++) { std::vector<TH1F*> a; histos[i].push_back(a); sprintf(tmp,"particles_longest_z_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,10,0,2)); sprintf(names[0],"particles_longest_z"); sprintf(tmp,"particles_longest_s_particle_s_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,2)); sprintf(names[1],"particles_longest_s_particle_s"); sprintf(tmp,"particles_primary_long_e_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,10,0,2)); sprintf(names[2],"particles_primary_long_e"); sprintf(tmp,"particles_elec_vise_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,10)); sprintf(names[3],"particles_elec_vise"); sprintf(tmp,"particles_primary_phi_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,1.6)); sprintf(names[4],"particles_primary_phi"); sprintf(tmp,"particles_prim_par_b_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,1.5)); sprintf(names[5],"particles_prim_par_b"); sprintf(tmp,"particles_mol_rad_r_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,0.2)); sprintf(names[6],"particles_mol_rad_r"); sprintf(tmp,"particles_emfrac_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,1.1)); sprintf(names[7],"particles_emfrac"); sprintf(tmp,"particles_ntot_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,20)); sprintf(names[8],"particles_ntot"); sprintf(tmp,"particles_total_long_e_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,200)); sprintf(names[9],"particles_total_long_e"); sprintf(tmp,"particles_weighted_phi_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,1.6)); sprintf(names[10],"particles_weighted_phi"); sprintf(tmp,"largest_frac_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,50,0,2)); sprintf(names[11],"largest_frac"); sprintf(tmp,"reco_frac_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,50,0,2)); sprintf(names[12],"reco_frac"); sprintf(tmp,"particles_totvise_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,20)); sprintf(names[13],"particles_totvise"); sprintf(tmp,"event_visenergy_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,20)); sprintf(names[14],"event_visenergy"); sprintf(tmp,"particles_rms_r_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,4)); sprintf(names[15],"particles_rms_r"); sprintf(tmp,"particles_prim_par_a_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,10)); sprintf(names[16],"particles_prim_par_a"); sprintf(tmp,"particles_prim_par_e0_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,10)); sprintf(names[17],"particles_prim_par_e0"); sprintf(tmp,"particles_prim_par_chisq_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0.1,500)); sprintf(names[18],"particles_prim_par_chisq"); sprintf(tmp,"particles_largest_particle_peakdiff_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,-5,5)); sprintf(names[19],"particles_largest_particle_peakdiff"); sprintf(tmp,"largest_cmp_chisqndf_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,2)); sprintf(names[20],"largest_cmp_chisqndf"); sprintf(tmp,"event_nstrips_%d_%d",i,j); histos[i][j].push_back(new TH1F(tmp,tmp,20,0,80)); sprintf(names[21],"event_nstrips"); } for(int j=0;j<nCut;j++) for(int i=0;i<(int)histos.size();i++) for(int k=0;k<nHistos;k++) histos[i][j][k]->Sumw2(); }

void MakePlots::SortOutStats (  TPad *  pad, float  fracPadx = 0.3, float  fracPady = 0.5, int  numPerCol = 4, float  startx = 0.98, float  starty = 0.98 )  [private]

Definition at line 691 of file MakePlots.cxx. Referenced by SavePlots(), and SaveRatioPlots(). { pad->Update(); pad->Modified(); TList *theList = pad->GetListOfPrimitives(); TIterator *iter = theList->MakeIterator(); TObject *ob; int numStats = 0; while((ob = iter->Next())){ if(ob->InheritsFrom("TPad")) { TPad *p = (TPad*) ob; SortOutStats(p,fracPadx,fracPady,numPerCol,startx,starty); } else if(ob->InheritsFrom("TH1")||ob->InheritsFrom("TH2") ||ob->InheritsFrom("TGraph")||ob->InheritsFrom("TProfile")){ if(ob->FindObject("stats")) numStats+=1; } else if(ob->InheritsFrom("TLegend")){ numStats+=1; } } if(numStats!=0){ float xwidth = fracPadx/float(int((numStats-1)/numPerCol)+1); float ywidth = 0.2; if(numStats>=numPerCol) ywidth = fracPady/numPerCol; else ywidth = fracPady/float(numStats%numPerCol); float xgap = 0.01; float ygap = 0.01; float x2 = startx; float y2 = starty; float x1 = x2-xwidth; float y1 = y2-ywidth; int cnt = 0; iter->Reset(); while((ob = iter->Next())){ int col = 1; if(ob->InheritsFrom("TH1")){ TH1 *h1 = (TH1*) ob; col = h1->GetLineColor(); } else if(ob->InheritsFrom("TProfile")){ TProfile *p1 = (TProfile*) ob; col = p1->GetLineColor(); } else if(ob->InheritsFrom("TGraph")){ TGraph *g1 = (TGraph*) ob; col = g1->GetLineColor(); } TPaveStats *stat = 0; if((stat = (TPaveStats*) ob->FindObject("stats"))){ stat->SetTextColor(col); stat->SetX1NDC(x1-int(cnt/numPerCol)*(xwidth+xgap)); stat->SetX2NDC(x2-int(cnt/numPerCol)*(xwidth+xgap)); stat->SetY1NDC(y1-float(cnt%numPerCol)*(ywidth+ygap)); stat->SetY2NDC(y2-float(cnt%numPerCol)*(ywidth+ygap)); cnt+=1; } else if(ob->InheritsFrom("TLegend")){ TLegend *leg = (TLegend*) ob; leg->SetX1NDC(x1-int(cnt/numPerCol)*(xwidth+xgap)); leg->SetX2NDC(x2-int(cnt/numPerCol)*(xwidth+xgap)); leg->SetY1NDC(y1-float(cnt%numPerCol)*(ywidth+ygap)); leg->SetY2NDC(y2-float(cnt%numPerCol)*(ywidth+ygap)); cnt+=1; } } } }

---

Member Data Documentation

std::vector<TChain *> MakePlots::chains [private]

Definition at line 67 of file MakePlots.h. Referenced by AddSet(), and FillHistos().

char MakePlots::cutnames[nCut][200] [private]

Definition at line 65 of file MakePlots.h. Referenced by Run(), and SavePlots().

Int_t MakePlots::event_contained [private]

Definition at line 95 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Int_t MakePlots::event_inFiducial [private]

Definition at line 94 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::event_max_z [private]

Definition at line 91 of file MakePlots.h. Referenced by SetBranches().

Double_t MakePlots::event_min_z [private]

Definition at line 92 of file MakePlots.h. Referenced by SetBranches().

Int_t MakePlots::event_nstrips [private]

Definition at line 93 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::event_pidA [private]

Definition at line 97 of file MakePlots.h. Referenced by SetBranches().

Double_t MakePlots::event_pidB [private]

Definition at line 98 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::event_pidC [private]

Definition at line 99 of file MakePlots.h. Referenced by SetBranches().

Double_t MakePlots::event_pidD [private]

Definition at line 100 of file MakePlots.h. Referenced by SetBranches().

Double_t MakePlots::event_visenergy [private]

Definition at line 90 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

std::vector<std::vector<std::vector<TH1F *> > > MakePlots::histos [private]

Definition at line 56 of file MakePlots.h. Referenced by AddSet(), FillHistos(), SavePlots(), SaveRatioPlots(), and SetupHistos().

Int_t MakePlots::isbeamve [private]

Definition at line 117 of file MakePlots.h.

Int_t MakePlots::iscc [private]

Definition at line 115 of file MakePlots.h.

Int_t MakePlots::isdis [private]

Definition at line 119 of file MakePlots.h.

std::vector<int> MakePlots::ismc [private]

Definition at line 62 of file MakePlots.h. Referenced by AddSet(), and FillHistos().

std::vector<int> MakePlots::ismrcc [private]

Definition at line 63 of file MakePlots.h. Referenced by AddSet(), and FillHistos().

Int_t MakePlots::isnc [private]

Definition at line 113 of file MakePlots.h.

Int_t MakePlots::isndis [private]

Definition at line 118 of file MakePlots.h.

Int_t MakePlots::isnue [private]

Definition at line 114 of file MakePlots.h.

Int_t MakePlots::istau [private]

Definition at line 116 of file MakePlots.h.

Double_t MakePlots::largest_frac [private]

Definition at line 110 of file MakePlots.h.

Int_t MakePlots::mctrue_iresonance [private]

Definition at line 105 of file MakePlots.h. Referenced by SetBranches().

Double_t MakePlots::mctrue_oscprob [private]

Definition at line 103 of file MakePlots.h. Referenced by SetBranches().

Double_t MakePlots::mctrue_totbeamweight [private]

Definition at line 104 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Int_t MakePlots::mctrue_type [private]

Definition at line 102 of file MakePlots.h. Referenced by SetBranches().

Double_t MakePlots::mrccinfo_particle_s [private]

Definition at line 130 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::mrccinfo_sum_e [private]

Definition at line 131 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

char MakePlots::names[nHistos][200] [private]

Definition at line 61 of file MakePlots.h. Referenced by SavePlots(), SaveRatioPlots(), and SetupHistos().

int MakePlots::nchains [private]

Definition at line 69 of file MakePlots.h. Referenced by AddSet(), and MakePlots().

char MakePlots::outdir[200] [private]

Definition at line 71 of file MakePlots.h. Referenced by MakePlots(), Run(), SavePlots(), SaveRatioPlots(), and SetOutDir().

char MakePlots::outfile[200] [private]

Definition at line 72 of file MakePlots.h. Referenced by MakePlots(), SavePlots(), SaveRatioPlots(), and SetOutFile().

Double_t MakePlots::particles_elec_vise [private]

Definition at line 76 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_emfrac [private]

Definition at line 81 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_largest_particle_cmp_chisq [private]

Definition at line 126 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_largest_particle_e [private]

Definition at line 86 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_largest_particle_peakdiff [private]

Definition at line 125 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Int_t MakePlots::particles_largest_particlecmp_ndf [private]

Definition at line 127 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_longest_particle_type [private]

Definition at line 129 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_longest_s_particle_s [private]

Definition at line 75 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_longest_z [private]

Definition at line 80 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_mol_rad_r [private]

Definition at line 79 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Int_t MakePlots::particles_ntot [private]

Definition at line 82 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_prim_par_a [private]

Definition at line 122 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_prim_par_b [private]

Definition at line 78 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_prim_par_chisq [private]

Definition at line 124 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_prim_par_e0 [private]

Definition at line 123 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_prim_vise [private]

Definition at line 89 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_primary_long_e [private]

Definition at line 74 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_primary_phi [private]

Definition at line 77 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_rms_r [private]

Definition at line 121 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_total_long_e [private]

Definition at line 83 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_totvise [private]

Definition at line 87 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

Double_t MakePlots::particles_weighted_phi [private]

Definition at line 84 of file MakePlots.h. Referenced by FillHistos(), and SetBranches().

std::vector<double> MakePlots::pots

Definition at line 24 of file MakePlots.h. Referenced by AddSet(), and FillHistos().

std::vector<std::string> MakePlots::rationame [private]

Definition at line 59 of file MakePlots.h. Referenced by SaveRatioPlots().

std::vector<std::pair<int,int> > MakePlots::ratiopair [private]

Definition at line 58 of file MakePlots.h. Referenced by SaveRatioPlots().

Double_t MakePlots::reco_frac [private]

Definition at line 111 of file MakePlots.h.

Double_t MakePlots::trainweight [private]

Definition at line 107 of file MakePlots.h.

Int_t MakePlots::type [private]

Definition at line 108 of file MakePlots.h.

std::vector<std::string> MakePlots::typesname [private]

Definition at line 68 of file MakePlots.h. Referenced by AddSet(), and SavePlots().

---

The documentation for this class was generated from the following files:

• MakePlots.h
• MakePlots.cxx

---