DCHgeometryconstants_bmn.h

Go to the documentation of this file.

#ifndef DCHGEOMETRYCONSTANTS_BMN_H
#define DCHGEOMETRYCONSTANTS_BMN_H
 
#include "TMath.h"
 
// Beam entering position
// static Double_t beamPos[3] = {0.,0.,-210.};
/*static Double_t beamPos[3][numRunPeriods] = {
                                             {0.,0.,-210.};
                                             {0.,0.,-210.};
                                             {0.,0.,0.};
                                             {0.,0.,0.};
                                             {0.,0.,0.};
                                             {0.,0.,-210.};
                                             {0.,0.,0.};
                                            }*/
// Basic detector parameters
const UShort_t numChambers = 2;                          // number of drift chambers
const UShort_t numLayers = 8;                            // number of anode wire planes in the drift chamber
const UShort_t numLayers_half = numLayers / 2;           // half number of anode wire planes in the drift chamber
const UShort_t twoNumLayers = numLayers * numChambers;   // number of anode wire planes in both the drift chambers
const UShort_t ngaps = numLayers / 2 - 1;                // number of gaps in the drift chamber
// const UShort_t dimzarray=numChambers*numLayers; // number of anode wire planes in the drift chambers
const UShort_t numRunPeriods = 7;
 
// const UShort_t runType=1; // run type: 1,2,3
// Detector's position
// const Double_t deflectAngle = 0.3/5.7 ;
// const Double_t deflectAngle = 0. ;
// const Double_t angleXZ = 0.029845130209; //1.9 grad = 1.71 degree
// const Double_t angleXZ = 0.04712389; //2.7 degree
// const Double_t angleXZ = -0.04712389; //2.7 degree
// const Double_t angleXZ = 0.087266463; //5 degrees
const Double_t angleXZ = 0.0;   // 0. grad = 0. degree
// static Double_t DCH1_Xpos[numRunPeriods] = {0.,0.,0.,0.,0.,0.,0.};
// static Double_t DCH1_Ypos[numRunPeriods] = {0.,0.,0.,0.,0.,0.,0.};
// static Double_t DCH1_Xpos[numRunPeriods] = {0.,0.,0.,0.,0.,-39.381,0.};
// static Double_t DCH1_Ypos[numRunPeriods] = {0.,0.,0.,0.,0.,1.2917,0.};
const Double_t DCH1_Zpos[numRunPeriods] = {610., 610., 662.05, 0., 0., 515.5, 0.};   // cm
// const Double_t DCH1_Zpos[3] = {550.0,637.85,637.85}; //cm
// static Double_t DCH2_Xpos[numRunPeriods] = {0.,0.,0.,0.,0.,0.,0.};
// static Double_t DCH2_Ypos[numRunPeriods] = {0.,0.,0.,0.,0.,0.,0.};
// static Double_t DCH2_Xpos[numRunPeriods] = {0.,0.,0.,0.,0.,-41.299,0.};
// static Double_t DCH2_Ypos[numRunPeriods] = {0.,0.,0.,0.,0.,2.2445,0.};
const Double_t DCH2_Zpos[numRunPeriods] = {720., 720., 772.05, 0., 0., 715., 0.};   // cm
// Double_t DCH2_Zpos[3] = {650.0,748.35,748.35}; //cm
// const Double_t DCH2_Zpos[3] = {650.0,673.,748.35}; //cm
 
// const Double_t DCH_ZlayerPos_local[numLayers] = {-7.8, -6.6, -3.0, -1.8, 1.8, 3.0, 6.6, 7.8}; // position of gas
// volumes (layers) within DCH
const Double_t DCH_ZlayerPos_local[numLayers] = {7.8,  6.6,  3.0, 1.8, -1.8,
                                                 -3.0, -6.6, -7.8};   // position of gas volumes (layers) within DCH
// const Double_t DCH_ZlayerPos_local[numLayers] = {3.0, 1.8, 7.8, 6.6, -6.6, -7.8, -1.8, -3.0}; // position of gas
// volumes (layers) within DCH const Double_t DCH_ZlayerPos_local[numLayers] = {9.3, 8.1, 3.5, 2.3, -2.3, -3.5, -8.1,
// -9.3}; // position of gas volumes (layers) within DCH const Double_t DCH_ZlayerPos_local[numLayers] =
// {10.2, 9., 3.8, 2.6, -2.6, -3.8, -9., -10.2}; // position of gas volumes (layers) within DCH const Double_t
// DCH_ZlayerPos_local[numLayers] = {-9.3, -8.1, -3.5, -2.3, 2.3, 3.5, 8.1, 9.3}; // position of gas volumes (layers)
// within DCH
 
// const Double_t DCH1_ZlayerPos_global[numLayers] =
// {DCH_ZlayerPos_local[0]+DCH1_Zpos,DCH_ZlayerPos_local[1]+DCH1_Zpos,DCH_ZlayerPos_local[2]+DCH1_Zpos,DCH_ZlayerPos_local[3]+DCH1_Zpos,DCH_ZlayerPos_local[4]+DCH1_Zpos,DCH_ZlayerPos_local[5]+DCH1_Zpos,DCH_ZlayerPos_local[6]+DCH1_Zpos,DCH_ZlayerPos_local[7]+DCH1_Zpos};
// // position of gas volumes (layers) of DCH1 in laboratory CS const Double_t DCH2_ZlayerPos_global[numLayers] =
// {DCH_ZlayerPos_local[0]+DCH2_Zpos,DCH_ZlayerPos_local[1]+DCH2_Zpos,DCH_ZlayerPos_local[2]+DCH2_Zpos,DCH_ZlayerPos_local[3]+DCH2_Zpos,DCH_ZlayerPos_local[4]+DCH2_Zpos,DCH_ZlayerPos_local[5]+DCH2_Zpos,DCH_ZlayerPos_local[6]+DCH2_Zpos,DCH_ZlayerPos_local[7]+DCH2_Zpos};
// // position of gas volumes (layers) of DCH2 in laboratory CS
 
// Detector's construct parameters
const Double_t angleLayer[numLayers] = {0.,  0.,  90.,  90.,
                                        45., 45., -45., -45.};   // angles of inclined wires, degrees // for exp.
// const Double_t angleLayer[numLayers] = {-45., -45., 45., 45., 90., 90., 0., 0.}; // angles of inclined wires, degrees
// // for exp. const Double_t angleLayer[numLayers] = {0., 0., 90., 90., -45., -45., 45., 45.}; // angles of inclined
// wires, degrees // for sim. const Double_t angleStereoRad=angleStereo*TMath::DegToRad();
const Double_t angleLayerRad[numLayers] = {
    angleLayer[0] * TMath::DegToRad(), angleLayer[1] * TMath::DegToRad(), angleLayer[2] * TMath::DegToRad(),
    angleLayer[3] * TMath::DegToRad(), angleLayer[4] * TMath::DegToRad(), angleLayer[5] * TMath::DegToRad(),
    angleLayer[6] * TMath::DegToRad(), angleLayer[7] * TMath::DegToRad()};   // angles  of inclined wires, radians
const Double_t dirCosLayerX[numLayers] = {
    cos(angleLayerRad[0]), cos(angleLayerRad[1]), cos(angleLayerRad[2]), cos(angleLayerRad[3]), cos(angleLayerRad[4]),
    cos(angleLayerRad[5]), cos(angleLayerRad[6]), cos(angleLayerRad[7])};   // x direction  cosines of inclined wires
const Double_t dirCosLayerY[numLayers] = {
    sin(angleLayerRad[0]), sin(angleLayerRad[1]), sin(angleLayerRad[2]), sin(angleLayerRad[3]), sin(angleLayerRad[4]),
    sin(angleLayerRad[5]), sin(angleLayerRad[6]), sin(angleLayerRad[7])};   // y direction  cosines of inclined wires
 
const Double_t ZLength_DCH = 20.0;
const Double_t ZLength_DCHActiveVolume = 0.6;
const Double_t ZhalfLength_DCHActiveVolume = ZLength_DCHActiveVolume / 2.;
const Double_t ZhalfLength_DCHActiveVolume_W = ZhalfLength_DCHActiveVolume + 0.1;   // wider
// const Double_t DCH1_ZlayerPos_global_End[numLayers] =
// {DCH1_ZlayerPos_global[0]+ZhalfLength_DCHActiveVolume,DCH1_ZlayerPos_global[1]+ZhalfLength_DCHActiveVolume,DCH1_ZlayerPos_global[2]+ZhalfLength_DCHActiveVolume,DCH1_ZlayerPos_global[3]+ZhalfLength_DCHActiveVolume,DCH1_ZlayerPos_global[4]+ZhalfLength_DCHActiveVolume,DCH1_ZlayerPos_global[5]+ZhalfLength_DCHActiveVolume,DCH1_ZlayerPos_global[6]+ZhalfLength_DCHActiveVolume,DCH1_ZlayerPos_global[7]+ZhalfLength_DCHActiveVolume};
// // ends of gas volumes (layers) of DCH1 in laboratory CS const Double_t DCH2_ZlayerPos_global_End[numLayers] =
// {DCH2_ZlayerPos_global[0]+ZhalfLength_DCHActiveVolume,DCH2_ZlayerPos_global[1]+ZhalfLength_DCHActiveVolume,DCH2_ZlayerPos_global[2]+ZhalfLength_DCHActiveVolume,DCH2_ZlayerPos_global[3]+ZhalfLength_DCHActiveVolume,DCH2_ZlayerPos_global[4]+ZhalfLength_DCHActiveVolume,DCH2_ZlayerPos_global[5]+ZhalfLength_DCHActiveVolume,DCH2_ZlayerPos_global[6]+ZhalfLength_DCHActiveVolume,DCH2_ZlayerPos_global[7]+ZhalfLength_DCHActiveVolume};
// // ends of gas volumes (layers) of DCH2 in laboratory CS
const Double_t HoleSize_DCH = 12.0;
 
const Double_t SideLengthOfOctagon = 120.0;
const Double_t InnerRadiusOfOctagon = (2.414 * SideLengthOfOctagon) / 2.0;
 
const Double_t MaxRadiusOfActiveVolume = 120.0;
const Double_t MinRadiusOfActiveVolume = 12.0;
 
// const UShort_t numWiresPerLayer = 256;
const UInt_t numWiresPerLayer = 240;
const UInt_t halfNumWiresPerLayer = numWiresPerLayer / 2;
// const Double_t stepXY=2.*MaxRadiusOfActiveVolume/Double_t(numWiresPerLayer);
const Double_t stepXY = 1.;
const Double_t stepXYhalf = stepXY / 2.;
const Int_t detIdDch1 = 13;
const Int_t detIdDch2 = 17;
// const Double_t diagonal=sqrt(stepXYhalf*stepXYhalf+ZhalfLength_DCHActiveVolume*ZhalfLength_DCHActiveVolume);
 
#endif