Add LUT maker for ITS1,2,3 (#26)

* Add Proto-LUT maker for ITS1,2,3

* Add carbon foam to the ITS3

Author: njacazio

src/lutWrite.its1.cc

@@ -0,0 +1,29 @@
+/// @author Nicolo' Jacazio
+/// @email: [email protected]
+/// @since  30/10/2020
+/// @file   lutWrite.its1.cc
+
+#include "lutWrite.cc"
+
+void fatInit_its1(float field = 0.5, float rmin = 100.) {
+
+  fat.SetBField(field);
+  fat.SetdNdEtaCent(400.);
+  fat.MakeAliceCurrent();
+  fat.SetAtLeastHits(4);
+  fat.SetAtLeastCorr(4);
+  fat.SetAtLeastFake(0);
+  //
+  fat.SetMinRadTrack(rmin);
+  //
+  fat.PrintLayout();
+}
+
+void lutWrite_its1(const char *filename = "lutCovm.dat", int pdg = 211,
+                   float field = 0.2, float rmin = 20.) {
+
+  // init FAT
+  fatInit_its1(field, rmin);
+  // write
+  lutWrite(filename, pdg, field, rmin);
+}

src/lutWrite.its2.cc

@@ -0,0 +1,47 @@
+/// @author Nicolo' Jacazio
+/// @email: [email protected]
+/// @since  30/10/2020
+/// @file   lutWrite.its2.cc
+
+#include "lutWrite.cc"
+
+void fatInit_its2(float field = 0.5, float rmin = 100.) {
+  fat.SetBField(field);
+  fat.SetdNdEtaCent(400.);
+
+  fat.AddLayer((char *)"bpipe", 2.0, 0.0022); // beam pipe
+  fat.AddLayer((char *)"vertex", 0, 0); // dummy vertex for matrix calculation
+
+  // new ideal Pixel properties?
+  Double_t x0 = 0.0050;
+  Double_t resRPhi = 0.0006;
+  Double_t resZ = 0.0006;
+  Double_t eff = 0.98;
+  Double_t xrhoOB = 1.1646e-01; // 500 mum Si
+
+  fat.AddLayer((char *)"ddd1", 2.3, x0, xrhoOB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"ddd2", 3.1, x0, xrhoOB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"ddd3", 3.9, x0, xrhoOB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"ddd4", 19.4, x0, xrhoOB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"ddd5", 24.7, x0, xrhoOB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"ddd6", 35.3, x0, xrhoOB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"ddd7", 40.5, x0, xrhoOB, resRPhi, resZ, eff);
+
+  fat.AddTPC(0.1, 0.1); // TPC
+  fat.SetAtLeastHits(4);
+  fat.SetAtLeastCorr(4);
+  fat.SetAtLeastFake(0);
+  //
+  fat.SetMinRadTrack(rmin);
+  //
+  fat.PrintLayout();
+}
+
+void lutWrite_its2(const char *filename = "lutCovm.dat", int pdg = 211,
+                   float field = 0.2, float rmin = 20.) {
+
+  // init FAT
+  fatInit_its2(field, rmin);
+  // write
+  lutWrite(filename, pdg, field, rmin);
+}

src/lutWrite.its3.cc

@@ -0,0 +1,52 @@
+/// @author Nicolo' Jacazio
+/// @email: [email protected]
+/// @since  30/10/2020
+/// @file   lutWrite.its3.cc
+
+#include "lutWrite.cc"
+
+void fatInit_its3(float field = 0.5, float rmin = 100.) {
+  fat.SetBField(field);
+  fat.SetdNdEtaCent(400.);
+
+  fat.AddLayer((char *)"bpipe", 1.6, 0.0022); // beam pipe
+  fat.AddLayer((char *)"vertex", 0, 0); // dummy vertex for matrix calculation
+
+  // new ideal Pixel properties?
+  Double_t x0IB = 0.0005;       // X/X0 of the inner barrel (first three layers)
+  Double_t x0OB = 0.0035;       // X/X0 of the outer barrel
+  Double_t resRPhi = 0.0006;    // Resolution in Rphi
+  Double_t resZ = 0.0006;       // Resolution in Z
+  Double_t eff = 0.98;          // Efficiency
+  Double_t xrhoIB = 1.1646e-02; // Surface density for 50 mum thick Si
+  Double_t xrhoOB = 1.1646e-01; // Surface density for 500 mum thick Si
+
+  fat.AddLayer((char *)"ddd1", 1.8, x0IB, xrhoIB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"foam1", 1.8 + 9.370 * x0IB, 0.0008); // Foam spacer
+  fat.AddLayer((char *)"ddd2", 2.4, x0IB, xrhoIB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"foam2", 2.4 + 9.370 * x0IB, 0.0008); // Foam spacer
+  fat.AddLayer((char *)"ddd3", 3.0, x0IB, xrhoIB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"foam3", 3.0 + 9.370 * x0IB, 0.0008); // Foam spacer
+  fat.AddLayer((char *)"ddd4", 19.4, x0OB, xrhoOB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"ddd5", 24.7, x0OB, xrhoOB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"ddd6", 35.3, x0OB, xrhoOB, resRPhi, resZ, eff);
+  fat.AddLayer((char *)"ddd7", 40.5, x0OB, xrhoOB, resRPhi, resZ, eff);
+
+  fat.AddTPC(0.1, 0.1); // TPC
+  fat.SetAtLeastHits(4);
+  fat.SetAtLeastCorr(4);
+  fat.SetAtLeastFake(0);
+  //
+  fat.SetMinRadTrack(rmin);
+  //
+  fat.PrintLayout();
+}
+
+void lutWrite_its3(const char *filename = "lutCovm.dat", int pdg = 211,
+                   float field = 0.2, float rmin = 20.) {
+
+  // init FAT
+  fatInit_its3(field, rmin);
+  // write
+  lutWrite(filename, pdg, field, rmin);
+}