Changeset 137 in JEM-EUSO

Timestamp:

May 7, 2013, 2:42:24 PM (11 years ago)

Author:

moretto

Message:

new reco root files and modification of the lecture macro to read those new files

Location:

esaf_lal/branches/camille/camille_work/comparison_new_old_code

Files:

• lecture.C (modified) (4 diffs)
• simu.mac.reco.CfgDump.cfg (added)
• simu.mac.reco.log (added)
• simu.mac.reco.root (added)
• simu.old_linux.reco.CfgDump.cfg (added)
• simu.old_linux.reco.log (added)
• simu.old_linux.reco.root (added)

esaf_lal/branches/camille/camille_work/comparison_new_old_code/lecture.C

@@ -24 +24 @@
     
     
-    TCanvas *c1 = new TCanvas("c1", "Comparison new code vs old => Red: new on mac / Black: old on linux");
+    TCanvas *c1 = new TCanvas("c1", "Comparison new code vs old (Simulation) => Red: new on mac / Black: old on linux");
     c1->Divide(2,2);
     
@@ -74 +74 @@
     
     
-    TCanvas *c2 = new TCanvas("c2", "Comparison new code vs old => Red: new on mac / Black: new on linux");
+    TCanvas *c2 = new TCanvas("c2", "Comparison new code vs old (Simulation) => Red: new on mac / Black: new on linux");
     c2->Divide(2,2);
     
@@ -124 +124 @@
     
     
-    TCanvas *c3 = new TCanvas("c3", "Comparison new code vs old => Red: new on linux / Black: old on linux");
+    TCanvas *c3 = new TCanvas("c3", "Comparison new code vs old (Simulation) => Red: new on linux / Black: old on linux");
     c3->Divide(2,2);
     
@@ -150 +150 @@
     h23->Draw();
     h24->Draw("Same");
+    
+    //---------------------------------------
+    
+    TFile *f7 = new TFile("/Users/moretto/Documents/thesis/ESAF/Code/esaf_lal/esaf/camille_work/comparison_new_old_code/simu.mac.reco.root");
+    TTree *recotree7 = (TTree*)f7->Get("recotree");
+    
+    TFile *f8 = new TFile("/Users/moretto/Documents/thesis/ESAF/Code/esaf_lal/esaf/camille_work/comparison_new_old_code/simu.old_linux.reco.root");
+    TTree *recotree8 = (TTree*)f8->Get("recotree");
+    
+    //option 1:
+    
+    recotree7->Project("h25", "fTrackDirection2.fTheta");
+    recotree8->Project("h26", "fTrackDirection2.fTheta");
+    
+    recotree7->Project("h27", "fTrackDirection2.fPhi");
+    recotree8->Project("h28", "fTrackDirection2.fPhi");
+    
+    recotree7->Project("h29", "fTrackDirection2.fErrorTheta");
+    recotree8->Project("h30", "fTrackDirection2.fErrorTheta");
+    
+    recotree7->Project("h31", "fTrackDirection2.fErrorPhi");
+    recotree8->Project("h32", "fTrackDirection2.fErrorPhi");
+    
+    
+    TCanvas *c4 = new TCanvas("c4", "Comparison new code vs old (Reconstruction) => Red: new on mac / Black: old on linux");
+    c4->Divide(2,2);
+    
+    c4->cd(1);
+    h25->SetLineColor(kRed);
+    h26->SetLineColor(kBlack);
+    h25->Draw();
+    h26->Draw("Same");
+    
+    c4->cd(2);
+    h27->SetLineColor(kRed);
+    h28->SetLineColor(kBlack);
+    h27->Draw();
+    h28->Draw("Same");
+    
+    c4->cd(3);
+    h29->SetLineColor(kRed);
+    h30->SetLineColor(kBlack);
+    h29->Draw();
+    h30->Draw("Same");
+    
+    c4->cd(4);
+    h31->SetLineColor(kRed);
+    h32->SetLineColor(kBlack);
+    h31->Draw();
+    h32->Draw("Same");
+