LMU/musrsim
LMU/musrsim
5
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Added Muonium TOF

Browse Source
This commit is contained in:
prokscha 2020-01-25 15:39:17 +01:00
parent 7ab43a5410
commit 0e0a1d64c6
1 changed files with 8 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
run/macros/td_mcp2_tof.C
View File

@ -17,7 +17,7 @@ void td_mcp2_tof(const char* fname)
t1 = (TTree*) gDirectory->FindObjectAny("t1"); t1 = (TTree*) gDirectory->FindObjectAny("t1");
// muon time-of-flight to TD // muon time-of-flight to TD
t1->Draw("eventID:1000*save_time","save_detID==890&&save_particleID==-13"); //0, 2 t1->Draw("eventID:1000*save_time","save_detID==890&&((save_particleID==-13)||(save_particleID==-1313))"); //0, 2
nEventsTD = t1->GetSelectedRows(); nEventsTD = t1->GetSelectedRows();
ptofTD = t1->GetV2(); ptofTD = t1->GetV2();
pEventIDTD = t1->GetV1(); pEventIDTD = t1->GetV1();
@ -27,11 +27,12 @@ void td_mcp2_tof(const char* fname)
tofTD[i] = ptofTD[i]; tofTD[i] = ptofTD[i];
} }
// muon time-of-flight to MCP2, mu+ // muon time-of-flight to MCP2, mu+ or Mu
t1->Draw("eventID:1000*save_time","save_detID==902&&save_particleID==-13"); //0, 2 t1->Draw("eventID:1000*save_time","save_detID==902&&((save_particleID==-13)||(save_particleID==-1313))"); //0, 2
nEventsMCP2 = t1->GetSelectedRows(); nEventsMCP2 = t1->GetSelectedRows();
ptofMCP2 = t1->GetV2(); ptofMCP2 = t1->GetV2();
pEventIDMCP2 = t1->GetV1(); pEventIDMCP2 = t1->GetV1();
// cout << "nEventsMCP2 = " << nEventsMCP2 << endl;
for (i = 0; i<nEventsMCP2; i++){ for (i = 0; i<nEventsMCP2; i++){
EventIDMCP2[i] = pEventIDMCP2[i]; EventIDMCP2[i] = pEventIDMCP2[i];
@ -39,11 +40,11 @@ void td_mcp2_tof(const char* fname)
} }
// TOF TD-MCP2 // TOF TD-MCP2
htof = new TH1D("Time-of-flight TD-MCP2", "Time-of-flight TD-MCP2", 1000, 0.25, 500.25); htof = new TH1D("tof_TDMCP2", "Time-of-flight TD-MCP2", 1000, 0.25, 500.25);
k = 0; k = 0;
for (i = 0; i<nEventsMCP2; i++){ for (i = 0; i<nEventsMCP2; i++){
for (j = k; j<nEventsTD; j++){ for (j = k; j<nEventsTD; j++){
// cout << i << ", " << j << ", " << EventIDMCP2[i] << ", " << EventIDTD[j] << ", " << tofMCP2[i] << ", " << tofTD[j] << endl; // cout << i << ", " << j << ", " << EventIDMCP2[i] << ", " << EventIDTD[j] << ", " << tofMCP2[i] << ", " << tofTD[j] << endl;
if (EventIDMCP2[i] == EventIDTD[j]){ if (EventIDMCP2[i] == EventIDTD[j]){
htof->Fill(tofMCP2[i] - tofTD[j]); htof->Fill(tofMCP2[i] - tofTD[j]);
k = j + 1; k = j + 1;