From 0e0a1d64c6ab622913f6e7a7372f302a5628315d Mon Sep 17 00:00:00 2001
From: prokscha <thomas.prokscha@psi.ch>
Date: Sat, 25 Jan 2020 15:39:17 +0100
Subject: [PATCH] Added Muonium TOF

---
 run/macros/td_mcp2_tof.C | 15 ++++++++-------
 1 file changed, 8 insertions(+), 7 deletions(-)

diff --git a/run/macros/td_mcp2_tof.C b/run/macros/td_mcp2_tof.C
index d97b1e4..2c33aa9 100644
--- a/run/macros/td_mcp2_tof.C
+++ b/run/macros/td_mcp2_tof.C
@@ -17,21 +17,22 @@ void td_mcp2_tof(const char* fname)
   t1 = (TTree*) gDirectory->FindObjectAny("t1");
 
   // 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();
   ptofTD     = t1->GetV2();
   pEventIDTD = t1->GetV1();
-
+  
   for (i = 0; i<nEventsTD; i++){
    EventIDTD[i] = pEventIDTD[i];
    tofTD[i]     = ptofTD[i];
   }
-
-  // muon time-of-flight to MCP2, mu+
-  t1->Draw("eventID:1000*save_time","save_detID==902&&save_particleID==-13"); //0, 2
+  
+  // muon time-of-flight to MCP2, mu+ or Mu
+  t1->Draw("eventID:1000*save_time","save_detID==902&&((save_particleID==-13)||(save_particleID==-1313))"); //0, 2
   nEventsMCP2  = t1->GetSelectedRows();
   ptofMCP2     = t1->GetV2();
   pEventIDMCP2 = t1->GetV1();
+//   cout << "nEventsMCP2 = " << nEventsMCP2 << endl;
 
   for (i = 0; i<nEventsMCP2; i++){
    EventIDMCP2[i] = pEventIDMCP2[i];
@@ -39,11 +40,11 @@ void td_mcp2_tof(const char* fname)
   }
 
   // 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;
   for (i = 0; i<nEventsMCP2; i++){
    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]){
        htof->Fill(tofMCP2[i] - tofTD[j]);
        k = j + 1;