change Gauss-Job.py LbDelphes.py and add xyReader.C

Sim/LbDelphes/options/Gauss-Job.py

@@ -19,12 +19,12 @@ from Configurables import LHCbApp
 LHCbApp().DDDBtag   = "dddb-20170721-3"
 LHCbApp().CondDBtag = "sim-20170721-2-vc-md100"
 
-#eventType = '56000400'#photon particle gun
+eventType = '56000400'#photon particle gun
 #eventType = '30000000'#min bias
 #eventType = '13102201'#Bs-> phi gamma
 #eventType = '562000001'#pi0 particle gun
 #eventType = '12163401'#B+-> D0 K, D0 -> K pi pi0 (CF)
-eventType = '13102201'#Bs->phi gamma, phi->kk
+#eventType = '13102201'#Bs->phi gamma, phi->kk
 #eventType="13112001"#Bs->MuMu
 
 #LHCbApp().OutputLevel=DEBUG

Sim/LbDelphes/options/LbDelphes.py

@@ -34,7 +34,10 @@ def delphesForGauss():
     #PGPrimaryVertex().OutputLevel=DEBUG
     PGPrimaryVertex().OutputVerticesName = '/Event/Rec/Vertex/Primary'
     PGPrimaryVertex().InputVerticesName = '/Event/MCFast/MCVertices'
-    DelphesAlg().LHCbDelphesCardLocation = '$LBDELPHESROOT/options/'+c.name+'.tcl'
+    # DelphesAlg().LHCbDelphesCardLocation = '$LBDELPHESROOT/options/'+c.name+'.tcl'
+    DelphesAlg().LHCbDelphesCardLocation = '$LBDELPHESROOT/options/cards/delphes_card_LHCb_EndVelo_Rich2_withEcal.tcl'
+    
+    # delphes_card_LHCb_EndVelo_Rich2_withEcal.tcl
     DelphesAlg().DelphesTrackLocation = 'FinalTrackMerger/tracks'
     #DelphesAlg().OutputLevel=DEBUG
     

Sim/LbDelphes/options/xyReader.C

+#include "TFile.h"
+#include "TTree.h"
+#include "TStyle.h"
+#include "TH2.h"
+#include "TGaxis.h"
+
+void xyReader(const char* filename){
+  TFile* file1 = TFile::Open(filename);
+  TDirectory* dir = file1->GetDirectory("DelphesTuple");
+  TTree* rec_tree = (TTree*)dir->Get("Rec");
+  TTree* gen_tree = (TTree*)dir->Get("Gen");
+  
+  rec_tree->SetMarkerStyle(2);
+  rec_tree->SetMarkerColor(kRed);
+  gen_tree->SetMarkerStyle(2);
+  gen_tree->SetMarkerColor(kBlue);
+  rec_tree->Draw("x:y>>hist1(80,-3880,3880,80,-3880,3880)","","");
+  gen_tree->Draw("(x+thetaX*(12650.5-z)):(y+thetaY*(12650.5-z))>>hist2(80,-3880,3880,80,-3880,3880)","","same");
+  
+  TLegend *legend = new TLegend(0.8,0.8,0.95,0.95);
+  legend->AddEntry("hist2", "Gen", "le");
+  legend->AddEntry("hist1", "Rec", "P");
+  legend->Draw();
+
+}