prepare anasen monte carlo simulation

.vscode/c_cpp_properties.json

@@ -11,6 +11,19 @@
             "cStandard": "c17",
             "cppStandard": "gnu++17",
             "intelliSenseMode": "linux-gcc-x64"
+        },
+        {
+            "name": "Ryan",
+            "includePath": [
+                "${workspaceFolder}/**",
+                "/home/ryan/Downloads/root_build/**",
+                "/home/ryan/Downloads/root_build/include"
+            ],
+            "defines": [],
+            "compilerPath": "/usr/bin/gcc",
+            "cStandard": "c17",
+            "cppStandard": "gnu++17",
+            "intelliSenseMode": "linux-gcc-x64"
         }
     ],
     "version": 4

.vscode/settings.json

@@ -77,6 +77,7 @@
     "typeinfo": "cpp",
     "variant": "cpp",
     "bitset": "cpp",
-    "forward_list": "cpp"
+    "forward_list": "cpp",
+    "anasenMS.C": "cpp"
   }
 }

ClassAnasen.h

@@ -433,9 +433,11 @@ inline  void ANASEN::DrawTrack(TVector3 pos, TVector3 direction){
   startPos->SetLineColor(kBlack);
   worldBox->AddNode(startPos, 3, new TGeoCombiTrans( pos.X(), pos.Y(), pos.Z(), new TGeoRotation("rotA", 0, 0, 0.)));
 
-  TGeoVolume * hit = geom->MakeSphere("hitpos", 0, 0, 3);
-  hit->SetLineColor(kRed);
-  worldBox->AddNode(hit, 2, new TGeoCombiTrans( sx3.hitPos.X(), sx3.hitPos.Y(), sx3.hitPos.Z(), new TGeoRotation("rotA", 0, 0, 0.)));
+  if( sx3.id >= 0 ){
+    TGeoVolume * hit = geom->MakeSphere("hitpos", 0, 0, 3);
+    hit->SetLineColor(kRed);
+    worldBox->AddNode(hit, 2, new TGeoCombiTrans( sx3.hitPos.X(), sx3.hitPos.Y(), sx3.hitPos.Z(), new TGeoRotation("rotA", 0, 0, 0.)));
+  }
   
   geom->CloseGeometry();
   geom->SetVisLevel(4);

anasenMS.C

@@ -0,0 +1,23 @@
+#include "ClassAnasen.h"
+
+#include "TRandom.h"
+#include "TH1.h"
+#include "TH2.h"
+#include "TStyle.h"
+#include "TCanvas.h"
+
+
+//======== Gerneate light particle based on reaction
+// find out the CalTrack and the real track
+// find out the Q-value uncertaintly
+
+void anasenMS(){
+
+
+  //Reaction
+
+
+
+
+
+}

script.C

@@ -257,7 +257,7 @@ void script(TString fileName = "", int maxEvent = -1){
 
   haha->CalTrack(sx3.hitPos, wireID.first, wireID.second, true);
 
-  haha->DrawDeducedTrack(sx3.hitPos, wireID.first, wireID.second);
+  //haha->DrawDeducedTrack(sx3.hitPos, wireID.first, wireID.second);
 
   // haha->DrawTrack(pos, dir);