| 1 | //**************************************************************************************************************
|
|---|
| 2 | //
|
|---|
| 3 | // Sample A ON = Runs 17206 - 17265
|
|---|
| 4 | // Sample A OFF = Runs 17269 - 17283
|
|---|
| 5 | //
|
|---|
| 6 | // Sample B ON = Runs 17286 - 17456
|
|---|
| 7 | // Sample B OFF = Runs 17459 - 17543
|
|---|
| 8 | //
|
|---|
| 9 | //**************************************************************************************************************
|
|---|
| 10 |
|
|---|
| 11 | //**************************************************************************************************************
|
|---|
| 12 | // general INPUT
|
|---|
| 13 | //**************************************************************************************************************
|
|---|
| 14 |
|
|---|
| 15 | // pi
|
|---|
| 16 | const double PI=3.14159265358979333; //atan(1)*4.0
|
|---|
| 17 |
|
|---|
| 18 | // conversion from deg to mm and vice versa
|
|---|
| 19 | const Double_t mm2deg = 180./17000./PI;
|
|---|
| 20 | const Double_t deg2mm = 1./mm2deg;
|
|---|
| 21 |
|
|---|
| 22 | // conversion from rad to deg and vice versa
|
|---|
| 23 | const Double_t rad2deg = 180./PI;
|
|---|
| 24 | const Double_t deg2rad = 1./rad2deg;
|
|---|
| 25 |
|
|---|
| 26 | // conversion from millisec to min
|
|---|
| 27 | const Double_t msec2min=1./1000./60.;
|
|---|
| 28 |
|
|---|
| 29 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 30 | // main path
|
|---|
| 31 |
|
|---|
| 32 | // HUB
|
|---|
| 33 | TString inpath = "/users/eeh/hengsteb/MAGIC/Mkn421Feb2004/SlidingWindow/";
|
|---|
| 34 | TString path = "/users/eeh/hengsteb/MAGIC/Mkn421Feb2004/";
|
|---|
| 35 |
|
|---|
| 36 | // laptop
|
|---|
| 37 | //TString inpath = "/home/hengsteb/MAGIC/Mkn421Feb2004/";
|
|---|
| 38 | //TString path = "/home/hengsteb/MAGIC/Mkn421Feb2004/";
|
|---|
| 39 |
|
|---|
| 40 | //**************************************************************************************************************
|
|---|
| 41 | // INPUT/OUTPUT for MagicHillas.C, ScaleHillas.C, SrcCorrect.C, FSrcPlot.C
|
|---|
| 42 | //**************************************************************************************************************
|
|---|
| 43 |
|
|---|
| 44 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 45 | // Mkn421 ON A
|
|---|
| 46 | TString dataOnA="SLIDINGRun17206-17265CalibRun17193_D_Mkn421On_E";
|
|---|
| 47 |
|
|---|
| 48 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 49 | // Mkn421 OFF A
|
|---|
| 50 | TString dataOffA="SLIDINGRun17269-17283CalibRun17268_D_OffMkn421_E";
|
|---|
| 51 |
|
|---|
| 52 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 53 | // Mkn421 ON B
|
|---|
| 54 | TString dataOnB="SLIDINGRun17286-17456CalibRun17285_D_Mkn421On_E";
|
|---|
| 55 |
|
|---|
| 56 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 57 | // Mkn421 OFF B
|
|---|
| 58 | TString dataOffB="SLIDINGRun17459-17543CalibRun17458_D_OffMkn421_E";
|
|---|
| 59 |
|
|---|
| 60 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 61 | // select the sample!
|
|---|
| 62 | TString data=dataOnB;
|
|---|
| 63 |
|
|---|
| 64 | TString fileCalibData = inpath + data + ".root";
|
|---|
| 65 | TString fileHillas = path + data + "_Hillas.root";
|
|---|
| 66 | TString fileHScaled = path + data + "_HScaled.root";
|
|---|
| 67 | TString fileHScalHist = path + data + "_HScalHist.root";
|
|---|
| 68 |
|
|---|
| 69 | TString fileFSrcPlot = path + data + "_FSrcPlot.root";
|
|---|
| 70 | TString fileSrcPos = path + data + "_SrcPos.dat";
|
|---|
| 71 |
|
|---|
| 72 | TString fileFSrcPlotC = path + data + "_FSrcPlotC.root"; //corrected with known src pos
|
|---|
| 73 | TString fileHScaledC = path + data + "_HScaledC.root"; //corrected with known src pos
|
|---|
| 74 |
|
|---|
| 75 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 76 | // dead pixels
|
|---|
| 77 |
|
|---|
| 78 | const Int_t NDEAD=54;
|
|---|
| 79 | const Short_t DEADPIX[NDEAD]=
|
|---|
| 80 | { 0, 25, 26, 34, 38, 46, 54, 124, 125, 155, 156, 157, 158, 162, 170, 172, 201, 202, 203, 204, 205,
|
|---|
| 81 | 206, 207, 208, 224, 388, 395 ,397, 440, 441, 442, 443, 444, 445, 446, 483, 484, 485, 486, 487, 488, 489,
|
|---|
| 82 | 490, 532, 533, 534, 535, 536, 537, 538, 539, 540, 543, 559 };
|
|---|
| 83 | // 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
|
|---|
| 84 |
|
|---|
| 85 | TArrayS blindpix(NDEAD,DEADPIX);
|
|---|
| 86 |
|
|---|
| 87 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 88 | // default cuts for false source plots
|
|---|
| 89 | /*
|
|---|
| 90 | //Double_t SizeLo = 1000.; // #ph
|
|---|
| 91 | //Double_t SizeLo = 1111.; // #ph
|
|---|
| 92 | Double_t SizeLo = 2000.; // #ph
|
|---|
| 93 | Double_t DistLo = 0.6; // degree
|
|---|
| 94 | Double_t DistUp = 1.4; // degree
|
|---|
| 95 |
|
|---|
| 96 | // static cuts
|
|---|
| 97 | //Double_t LengthLo = 0.14*deg2mm; // degree
|
|---|
| 98 | //Double_t LengthUp = 0.26*deg2mm; // degree
|
|---|
| 99 | //Double_t WidthLo = 0.05*deg2mm; // degree
|
|---|
| 100 | //Double_t WidthUp = 0.14*deg2mm; // degree
|
|---|
| 101 |
|
|---|
| 102 | // scaled Hillas
|
|---|
| 103 |
|
|---|
| 104 | Double_t LengthLo = 0.2525;
|
|---|
| 105 | Double_t LengthUp = 0.9312;
|
|---|
| 106 | Double_t WidthLo = 0.2875;
|
|---|
| 107 | Double_t WidthUp = 0.738;
|
|---|
| 108 | //*/
|
|---|
| 109 |
|
|---|
| 110 | //**************************************************************************************************************
|
|---|
| 111 | // INPUT/OUTPUT for AlphaPlot.C, CutOptim.C
|
|---|
| 112 | //**************************************************************************************************************
|
|---|
| 113 |
|
|---|
| 114 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 115 | // Mkn421 A + B
|
|---|
| 116 |
|
|---|
| 117 | // AlphaPlot + LightCurve, CutOptim
|
|---|
| 118 | TString fileOn1 = path + dataOnA + "_HScaledC.root";
|
|---|
| 119 | TString fileOn2 = path + dataOnB + "_HScaledC.root";
|
|---|
| 120 |
|
|---|
| 121 | TString fileOff1 = path + dataOffA + "_HScaled.root";
|
|---|
| 122 | TString fileOff2 = path + dataOffB + "_HScaled.root";
|
|---|
| 123 |
|
|---|
| 124 | TString fileRateCFac = path + "Mkn421A+B_CFac.dat";
|
|---|
| 125 | TString fileLCurve = path + "Mkn421A+B_LCurve.dat";
|
|---|
| 126 | TString fileAPlots = path + "Mkn421A+B_APlot.root";
|
|---|
| 127 |
|
|---|
| 128 | // matrices needed for cut optimization, optimize on sample A
|
|---|
| 129 | TString fileOptOn = path + dataOnA + "_HScaledC.root";
|
|---|
| 130 | TString fileOptOff = path + dataOffA + "_HScaled.root";
|
|---|
| 131 |
|
|---|
| 132 | TString fileMatOn = path + dataOnA + "_Matrix.root";
|
|---|
| 133 | TString fileMatOff = path + dataOffA + "_Matrix.root";
|
|---|
| 134 |
|
|---|
| 135 | //--------------------------------------------------------------------------------------------------------------
|
|---|
| 136 | // optimized cuts for 1-dim aplot
|
|---|
| 137 |
|
|---|
| 138 | //Double_t SizeLo = 1111.; // #ph
|
|---|
| 139 | Double_t SizeLo = 2000.; // #ph
|
|---|
| 140 | Double_t DistLo = 0.55;
|
|---|
| 141 | Double_t DistUp = 1.27;
|
|---|
| 142 | //Double_t DistLo = 1.;
|
|---|
| 143 | //Double_t DistUp = 1.5;
|
|---|
| 144 | //Double_t DistLo = 0.5;
|
|---|
| 145 | //Double_t DistUp = 1.22;
|
|---|
| 146 |
|
|---|
| 147 | // optimized on half of sample A for scaled Hillas with dist-cut 0.6 - 1.25, size > 1300
|
|---|
| 148 | Double_t LengthLo = 0.318;
|
|---|
| 149 | Double_t LengthUp = 0.809;
|
|---|
| 150 | Double_t WidthLo = 0.3;
|
|---|
| 151 | Double_t WidthUp = 0.705;
|
|---|
| 152 | //*/
|
|---|
