1 #if !defined(__CLING__) 11 #include "~/Acculinna/ER/beamtime/data/ERBeamTimeEventHeader.h" 23 TString inFile =
"he8_10_0010.reco.root";
24 TString inPath =
"~/Acculinna/ER/showBeam_VC/reco/";
26 TFile *fr =
new TFile(inPath + inFile);
27 TTree *tr = (TTree*)fr->Get(
"reco");
29 TFile *frIn =
new TFile(
"~/Acculinna/ER/showBeam_VC/analysed/he8_10_0010.Digi.root");
30 TTree *trIn = (TTree*)frIn->Get(
"er");
40 TCanvas *c1 =
new TCanvas(
"c1",
"diagnostics", 1000, 1000);
44 tr->Draw(
"trigger",
"");
49 tr->Draw(
"BeamDetMWPCDigiY1.fWireNb:BeamDetMWPCDigiX1.fWireNb",
"BeamDetMWPCDigiY1.fWireNb",
"col");
62 TCanvas *c2 =
new TCanvas(
"c2",
"MWPC projections", 1500, 1000);
66 tr->Draw(
"trigger",
"");
67 tr->Draw(
"BeamDetMWPCDigiY1.fWireNb:BeamDetMWPCDigiX1.fWireNb",
"trigger==2 && BeamDetMWPCDigiY1.fWireNb>0",
"col");
72 tr->Draw(
"BeamDetMWPCDigiY1.fWireNb:BeamDetMWPCDigiX1.fWireNb",
"trigger==3 && BeamDetMWPCDigiY1.fWireNb>0",
"col");
77 tr->Draw(
"BeamDetMWPCDigiY2.fWireNb:BeamDetMWPCDigiX2.fWireNb",
"trigger==2 && BeamDetMWPCDigiY2.fWireNb>0",
"col");
81 tr->Draw(
"BeamDetMWPCDigiY2.fWireNb:BeamDetMWPCDigiX2.fWireNb",
"trigger==3 && BeamDetMWPCDigiY2.fWireNb>0",
"col");
85 tr->Draw(
"fYt:fXt",
"trigger==2 && fYt>-40. && fXt>-40.");
86 tr->SetMarkerColor(kRed);
87 tr->Draw(
"fYt:fXt",
"trigger==2 && ( (fYt)^2 + (fXt)^2 ) < 8.^2",
"same");
91 tr->SetMarkerColor(kBlack);
92 tr->Draw(
"fYt:fXt",
"trigger==3 && fYt>-40. && fXt>-40.");
93 tr->SetMarkerColor(kRed);
94 tr->Draw(
"fYt:fXt",
"trigger==3 && ( (fYt-2)^2 + (fXt+3)^2 ) < 9.^2",
"same");