FEDRA emulsion software from the OPERA Collaboration
check_shower.C File Reference

Functions

cd (1)
 
cd (2)
 
cd (3)
 
cd (4)
 
cd (5)
 
cd (6)
 
cout<< endl;cout<< "Plot 1 : Number of BT per shower"<< endl;cout<< "Plot 2 : Neural Network output with 15 plates (1:e-; 0:pion)"<< endl;cout<< "Plot 3 : NN output versus Nb of BT (showers with 90% efficiency & 1% pion contamination)"<< endl;cout<< "Plot 4 : Shower (X-Y) with 90% eff & 1% pion contamination & highest Nb of BT"<< endl;cout<< "Plot 5 : Shower (Y-Z) with 90% eff & 1% pion contamination & highest Nb of BT"<< endl;cout<< "Plot 6 : Shower (X-Y-Z) with 90% eff & 1% pion contamination & highest Nb of BT"<< endl;TCanvas *c=new TCanvas("check_shower","check_shower", 1024, 900);c-> Divide (2, 3)
 
t Draw ("output15:sizeb>>NN15vsNbBT", cut3 &&cut2 &&cut)
 
t Draw ("output15>>NN15", cut2)
 
t Draw ("sizeb>>NbBT", cut2)
 
t Draw ("xb:yb:zb>>xyz", cut &&cut2 &&cut3 &&cut4)
 
t Draw ("xb:yb>>xy", cut &&cut2 &&cut3 &&cut4)
 
t Draw ("yb:zb>>yz", cut &&cut2 &&cut3 &&cut4)
 
NN15vsNbBT Draw ()
 
 for (int i=0;i< nentries;i++)
 
 if (max<=10) return
 
t SetBranchAddress ("sizeb", &sizeb)
 
NN15vsNbBT SetMarkerStyle (20)
 
NN15 SetTitle ("NN ouput with 15 plates")
 
NN15vsNbBT SetTitle ("NN output (15 plates) VS Nb of BT (90% Eff & 1% Cont)")
 
NbBT SetTitle ("Number of BT per shower")
 

Variables

TCut cut = "eProb90==1&&eProb1==1"
 
TCut cut2 = "output15>=0"
 
TCut cut3 = "sizeb>10"
 
TCut cut4 =s
 
int max =0
 
int nentries = t->GetEntries()
 
 s = max
 
int sizeb
 
TTree * t = (TTree*)f.Get("treebranch")
 

Function Documentation

◆ cd() [1/6]

canv cd ( )

◆ cd() [2/6]

canv cd ( )

◆ cd() [3/6]

canv cd ( )

◆ cd() [4/6]

c cd ( )

◆ cd() [5/6]

c cd ( )

◆ cd() [6/6]

c cd ( )

◆ Divide()

cout<< endl; cout<< "Plot 1 : Number of BT per shower"<< endl; cout<< "Plot 2 : Neural Network output with 15 plates (1:e-; 0:pion)"<< endl; cout<< "Plot 3 : NN output versus Nb of BT (showers with 90% efficiency & 1% pion contamination)"<< endl; cout<< "Plot 4 : Shower (X-Y) with 90% eff & 1% pion contamination & highest Nb of BT"<< endl; cout<< "Plot 5 : Shower (Y-Z) with 90% eff & 1% pion contamination & highest Nb of BT"<< endl; cout<< "Plot 6 : Shower (X-Y-Z) with 90% eff & 1% pion contamination & highest Nb of BT"<< endl; TCanvas *c=new TCanvas("check_shower","check_shower", 1024, 900); c-> Divide ( ,
 
)

◆ Draw() [1/7]

t Draw ( "output15:sizeb>>NN15vsNbBT"  ,
cut3 &&cut2 &&  cut 
)

◆ Draw() [2/7]

t Draw ( "output15>>NN15"  ,
cut2   
)

◆ Draw() [3/7]

t Draw ( "sizeb>>NbBT"  ,
cut2   
)

◆ Draw() [4/7]

t Draw ( "xb:yb:zb>>xyz"  ,
cut &&cut2 &&cut3 &&  cut4 
)

◆ Draw() [5/7]

t Draw ( "xb:yb>>xy ,
cut &&cut2 &&cut3 &&  cut4 
)

◆ Draw() [6/7]

t Draw ( "yb:zb>>yz"  ,
cut &&cut2 &&cut3 &&  cut4 
)

◆ Draw() [7/7]

pt Draw ( )

◆ for()

for ( )

43 {
44 t->GetEntry(i);
45
46 if (sizeb>max)
47 max=sizeb;
48 }
TTree * t
Definition: check_shower.C:4
int sizeb
Definition: check_shower.C:38
int max
Definition: check_shower.C:41

◆ if()

if ( max<=  10)

◆ SetBranchAddress()

t SetBranchAddress ( "sizeb"  ,
sizeb 
)

◆ SetMarkerStyle()

xyz SetMarkerStyle ( 20  )

◆ SetTitle() [1/3]

NN15 SetTitle ( "NN ouput with 15 plates )

◆ SetTitle() [2/3]

NN15vsNbBT SetTitle ( "NN output (15 plates) VS Nb of BT (90% Eff & 1% Cont)"  )

◆ SetTitle() [3/3]

NbBT SetTitle ( "Number of BT per shower"  )

Variable Documentation

◆ cut

TCut cut = "eProb90==1&&eProb1==1"

◆ cut2

TCut cut2 = "output15>=0"

◆ cut3

TCut cut3 = "sizeb>10"

◆ cut4

cut4 =s

◆ max

int max =0

◆ nentries

int nentries = t->GetEntries()

◆ s

TString s = max

◆ sizeb

int sizeb

◆ t

TTree* t = (TTree*)f.Get("treebranch")