New recoil detector for RGH (2 sectors, 3 layers/sector of urwell, x/…

common-tools/clas-detector/src/main/java/org/jlab/detector/base/DetectorType.java

@@ -31,6 +31,7 @@ public enum DetectorType {
     URWELL    (23, "URWELL"),
     AHDC      (24, "AHDC"),
     ATOF      (25, "ATOF"),
+    RECOIL    (26, "RECOIL"),
     TARGET    (100, "TARGET"),
     MAGNETS   (101, "MAGNETS"),
     ECIN      (110, "ECIN"),

common-tools/clas-jcsg/src/main/java/org/jlab/detector/geant4/v2/recoil/RecoilConstants.java

@@ -0,0 +1,108 @@
+package org.jlab.detector.geant4.v2.recoil;
+
+
+import org.jlab.detector.calib.utils.DatabaseConstantProvider; 
+import org.jlab.geom.prim.Point3D;
+
+
+public class RecoilConstants {
+
+    private final static String CCDBPATH = "/geometry/recoil/";
+
+    public final static int NMAXREGIONS = 3;    //max number of regions 
+    public final static int NREGIONS    = 3;    //number of regions 
+    public final static int NSECTORS    = 2;    //number of sectors
+    public final static int NLAYERS     = 2;    //number of layers
+    public final static int NCHAMBERS   = 1;    //number of chambers in a sector
+
+    public final static double HORIZONTHAL_OPENING_ANGLE = 34.;
+    public final static double VERTICAL_OPENING_ANGLE = 50.;
+    public final static double RADIUS[] = {33.5,55.5,79.5};
+    public final static double WIDTH[] = new double[NMAXREGIONS];
+    public final static double HEIGHT[] = new double[NMAXREGIONS];
+
+    public final static double THTILT = 0;            // theta tilt (deg)
+    /*    public final static double XENLARGEMENT = 0.5; // cm
+    public final static double YENLARGEMENT = 1.;  // cm
+    public final static double ZENLARGEMENT = 0.1; // cm
+
+    // Sector geometrical parameters
+    public final static double THOPEN = 34.;           // opening angle between endplate planes (deg)
+    public final static double THTILT = 0;            // theta tilt (deg)
+    public final static double THMIN  = 4.694;         // polar angle to the base of first chamber (deg)
+    public final static double SECTORHEIGHT = 146.21;  //height of each sector (cm)
+    public final static double DX0CHAMBER0  = 5.197;   // halfbase of chamber 1  (cm)*/
+
+    // Chamber volumes  and materials (units are cm)
+    public final static double[] CHAMBERVOLUMESTHICKNESS = {0.0025, 0.0005,0.3,                                // window
+                                                            0.0025, 0.0005,0.4,                                // cathode
+                                                            0.0005, 0.005, 0.0005,                             // uRWell + DlC
+                                                            0.0005, 0.005, 0.0005,                             // Capacitive sharing layer1
+                                                            0.0005, 0.005, 0.0005,                             // Capacitive sharing layer2
+                                                            0.005,  0.0005,0.005, 0.005,  0.0005,0.005, 0.005, // Readout
+                                                            0.0127, 0.3, 0.0125};                              // support
+    public final static String[] CHAMBERVOLUMESNAME = {"window_kapton", "window_Al", "window_gas",
+           "cathode_kapton", "cathode_Al", "cathode_gas",
+           "muRwell_Cu", "muRwell_kapton", "muRwell_dlc", 
+           "capa_sharing_layer1_glue","capa_sharing_layer1_Cr","capa_sharing_layer1_kapton",
+           "capa_sharing_layer2_glue","capa_sharing_layer2_Cr","capa_sharing_layer2_kapton",
+           "readout1_glue", "readout1_Cu", "readout1_kapton", "readout2_glue", "readout2_Cu", "readout2_kapton", "readout3_glue",
+           "support_skin1_g10", "support_honeycomb_nomex", "support_skin2_g10"};
+
+    // URWELL position in the CLAS12 frame 
+    /*    public final static double TGT2DC0    = 228.078; // cm            
+   // public final static double URWELL2DC0 = 2;       // cm
+    public final static double URWELL2DC0[] = new double[NMAXREGIONS];
+    public final static double DIST2TGT[] = new double[NMAXREGIONS];
+    public final static double W2TGT[] = new double[NMAXREGIONS];; 
+    public final static double YMIN[] = new double[NMAXREGIONS];
+    public final static double ZMIN[] = new double[NMAXREGIONS];*/
+
+    public final static double PITCH = 0.1 ;       // cm
+    public final static double STEREOANGLE = 90;   // deg
+
+     /*
+     * @return String a path to a directory in CCDB of the format {@code "/geometry/detector/"}
+     */
+    public static String getCcdbPath()
+    {
+            return CCDBPATH;
+    }
+
+     /**
+     * Loads the the necessary tables for the URWELL geometry for a given DatabaseConstantProvider.
+     * 
+     * @return DatabaseConstantProvider the same thing
+     */
+    public static DatabaseConstantProvider connect( DatabaseConstantProvider cp )
+    {
+          //  cp.loadTable( CCDBPATH +"RWELL");
+
+            load(cp  );
+            return cp;
+    }
+
+    /**
+     * Reads all the necessary constants from CCDB into static variables.
+     * Please use a DatabaseConstantProvider to access CCDB and load the following tables:
+     * @param cp a ConstantProvider that has loaded the necessary tables
+     */
+
+    public static synchronized void load( DatabaseConstantProvider cp )
+    {
+            // read constants from svt table
+//            NREGIONS = cp.getInteger( CCDBPATH+"svt/nRegions", 0 );
+
+             for (int i=0; i<NMAXREGIONS; i++){
+
+/*                URWELL2DC0[i] =  -2.+i*1.3;
+		  DIST2TGT[i]   = (TGT2DC0+URWELL2DC0[i]);
+		  W2TGT[i] = DIST2TGT[i]/Math.cos(Math.toRadians(THTILT-THMIN));
+		  YMIN[i]= W2TGT[i]*Math.sin(Math.toRadians(THMIN)); // distance from the base chamber1 and beamline
+		  ZMIN[i] = W2TGT[i]*Math.cos(Math.toRadians(THMIN));  
+		 */
+		 WIDTH[i]=2.*RADIUS[i]*Math.sin(Math.toRadians(HORIZONTHAL_OPENING_ANGLE)/2);
+		 HEIGHT[i]=RADIUS[i]*Math.tan(Math.toRadians(VERTICAL_OPENING_ANGLE)/2);
+	     }
+    }
+}