HcalCalibAlgos/macros/print_coef_endcap_plus_8.9mln.C

0001 {
0002 gROOT->Reset();
0003 gROOT->SetStyle("Plain");
0004
0005 gStyle->SetOptStat(1111);
0006 gStyle->SetOptFit(111);
0007
0008
0009 TH1F  *h1etacoefmin16_3a = new TH1F("h1etacoefmin16_3a", "h1etacoefmin16_3a", 100, 0., 2.);
0010 TH1F  *h1etacoefmin17a = new TH1F("h1etacoefmin17a", "h1etacoefmin17a", 100, 0.6, 1.4);
0011 TH1F  *h1etacoefmin18a = new TH1F("h1etacoefmin18a", "h1etacoefmin18a", 100, 0.6, 1.4);
0012 TH1F  *h1etacoefmin19a = new TH1F("h1etacoefmin19a", "h1etacoefmin19a", 100, 0.6, 1.4);
0013 TH1F  *h1etacoefmin20a = new TH1F("h1etacoefmin20a", "h1etacoefmin20a", 100, 0.6, 1.4);
0014 TH1F  *h1etacoefmin21a = new TH1F("h1etacoefmin21a", "h1etacoefmin21a", 100, 0.6, 1.4);
0015 TH1F  *h1etacoefmin22a = new TH1F("h1etacoefmin22a", "h1etacoefmin22a", 100, 0.6, 1.4);
0016 TH1F  *h1etacoefmin23a = new TH1F("h1etacoefmin23a", "h1etacoefmin23a", 100, 0.6, 1.4);
0017 TH1F  *h1etacoefmin24a = new TH1F("h1etacoefmin24a", "h1etacoefmin24a", 100, 0.6, 1.4);
0018 TH1F  *h1etacoefmin25a = new TH1F("h1etacoefmin25a", "h1etacoefmin25a", 100, 0.7, 1.3);
0019 TH1F  *h1etacoefmin26a = new TH1F("h1etacoefmin26a", "h1etacoefmin26a", 100, 0.7, 1.3);
0020 TH1F  *h1etacoefmin27a = new TH1F("h1etacoefmin27a", "h1etacoefmin27a", 100, 0.7, 1.3);
0021 TH1F  *h1etacoefmin28a = new TH1F("h1etacoefmin28a", "h1etacoefmin28a", 100, 0.7, 1.3);
0022 TH1F  *h1etacoefmin29a = new TH1F("h1etacoefmin29a", "h1etacoefmin29a", 100, 0.7, 1.3);
0023
0024 TH1F  *h1etacoefmin16_3b = new TH1F("h1etacoefmin16_3b", "h1etacoefmin16_3b", 100, 0., 2.);
0025 TH1F  *h1etacoefmin17b = new TH1F("h1etacoefmin17b", "h1etacoefmin17b", 100, 0.6, 1.4);
0026 TH1F  *h1etacoefmin18b = new TH1F("h1etacoefmin18b", "h1etacoefmin18b", 100, 0.6, 1.4);
0027 TH1F  *h1etacoefmin19b = new TH1F("h1etacoefmin19b", "h1etacoefmin19b", 100, 0.6, 1.4);
0028 TH1F  *h1etacoefmin20b = new TH1F("h1etacoefmin20b", "h1etacoefmin20b", 100, 0.6, 1.4);
0029 TH1F  *h1etacoefmin21b = new TH1F("h1etacoefmin21b", "h1etacoefmin21b", 100, 0.6, 1.4);
0030 TH1F  *h1etacoefmin22b = new TH1F("h1etacoefmin22b", "h1etacoefmin22b", 100, 0.6, 1.4);
0031 TH1F  *h1etacoefmin23b = new TH1F("h1etacoefmin23b", "h1etacoefmin23b", 100, 0.6, 1.4);
0032 TH1F  *h1etacoefmin24b = new TH1F("h1etacoefmin24b", "h1etacoefmin24b", 100, 0.6, 1.4);
0033 TH1F  *h1etacoefmin25b = new TH1F("h1etacoefmin25b", "h1etacoefmin25b", 100, 0.7, 1.3);
0034 TH1F  *h1etacoefmin26b = new TH1F("h1etacoefmin26b", "h1etacoefmin26b", 100, 0.7, 1.3);
0035 TH1F  *h1etacoefmin27b = new TH1F("h1etacoefmin27b", "h1etacoefmin27b", 100, 0.7, 1.3);
0036 TH1F  *h1etacoefmin28b = new TH1F("h1etacoefmin28b", "h1etacoefmin28b", 100, 0.7, 1.3);
0037 TH1F  *h1etacoefmin29b = new TH1F("h1etacoefmin29b", "h1etacoefmin29b", 100, 0.7, 1.3);
0038
0039 TH1F  *h1etacoefmin16_3 = new TH1F("h1etacoefmin16_3", "h1etacoefmin16_3", 100, 0., 2.);
0040 TH1F  *h1etacoefmin17 = new TH1F("h1etacoefmin17", "h1etacoefmin17", 100, 0.6, 1.4);
0041 TH1F  *h1etacoefmin18 = new TH1F("h1etacoefmin18", "h1etacoefmin18", 100, 0.6, 1.4);
0042 TH1F  *h1etacoefmin19 = new TH1F("h1etacoefmin19", "h1etacoefmin19", 100, 0.6, 1.4);
0043 TH1F  *h1etacoefmin20 = new TH1F("h1etacoefmin20", "h1etacoefmin20", 100, 0.6, 1.4);
0044 TH1F  *h1etacoefmin21 = new TH1F("h1etacoefmin21", "h1etacoefmin21", 100, 0.6, 1.4);
0045 TH1F  *h1etacoefmin22 = new TH1F("h1etacoefmin22", "h1etacoefmin22", 100, 0.6, 1.4);
0046 TH1F  *h1etacoefmin23 = new TH1F("h1etacoefmin23", "h1etacoefmin23", 100, 0.6, 1.4);
0047 TH1F  *h1etacoefmin24 = new TH1F("h1etacoefmin24", "h1etacoefmin24", 100, 0.6, 1.4);
0048 TH1F  *h1etacoefmin25 = new TH1F("h1etacoefmin25", "h1etacoefmin25", 100, 0.7, 1.3);
0049 TH1F  *h1etacoefmin26 = new TH1F("h1etacoefmin26", "h1etacoefmin26", 100, 0.7, 1.3);
0050 TH1F  *h1etacoefmin27 = new TH1F("h1etacoefmin27", "h1etacoefmin27", 100, 0.7, 1.3);
0051 TH1F  *h1etacoefmin28 = new TH1F("h1etacoefmin28", "h1etacoefmin28", 100, 0.7, 1.3);
0052 TH1F  *h1etacoefmin29 = new TH1F("h1etacoefmin29", "h1etacoefmin29", 100, 0.7, 1.3);
0053 // Two-dim
0054
0055
0056 TH2F  *h2etacoefmin16_3a = new TH2F("h2etacoefmin16_3a", "h2etacoefmin16_3a",72, 0.5, 72.5, 100, 0., 2.);
0057 TH2F  *h2etacoefmin17a = new TH2F("h2etacoefmin17a", "h2etacoefmin17a",72, 0.5, 72.5, 100, 0.6, 1.4);
0058 TH2F  *h2etacoefmin18a = new TH2F("h2etacoefmin18a", "h2etacoefmin18a",72, 0.5, 72.5, 100, 0.6, 1.4);
0059 TH2F  *h2etacoefmin19a = new TH2F("h2etacoefmin19a", "h2etacoefmin19a",72, 0.5, 72.5, 100, 0.6, 1.4);
0060 TH2F  *h2etacoefmin20a = new TH2F("h2etacoefmin20a", "h2etacoefmin20a",72, 0.5, 72.5,100, 0.6, 1.4);
0061 TH2F  *h2etacoefmin21a = new TH2F("h2etacoefmin21a", "h2etacoefmin21a",72, 0.5, 72.5,100, 0.6, 1.4);
0062 TH2F  *h2etacoefmin22a = new TH2F("h2etacoefmin22a", "h2etacoefmin22a",72, 0.5, 72.5,100, 0.6, 1.4);
0063 TH2F  *h2etacoefmin23a = new TH2F("h2etacoefmin23a", "h2etacoefmin23a",72, 0.5, 72.5,100, 0.6, 1.4);
0064 TH2F  *h2etacoefmin24a = new TH2F("h2etacoefmin24a", "h2etacoefmin24a",72, 0.5, 72.5, 100, 0.6, 1.4);
0065 TH2F  *h2etacoefmin25a = new TH2F("h2etacoefmin25a", "h2etacoefmin25a",72, 0.5, 72.5, 100, 0.7, 1.3);
0066 TH2F  *h2etacoefmin26a = new TH2F("h2etacoefmin26a", "h2etacoefmin26a",72, 0.5, 72.5, 100, 0.7, 1.3);
0067 TH2F  *h2etacoefmin27a = new TH2F("h2etacoefmin27a", "h2etacoefmin27a",72, 0.5, 72.5, 100, 0.7, 1.3);
0068 TH2F  *h2etacoefmin28a = new TH2F("h2etacoefmin28a", "h2etacoefmin28a",72, 0.5, 72.5, 100, 0.7, 1.3);
0069 TH2F  *h2etacoefmin29a = new TH2F("h2etacoefmin29a", "h2etacoefmin29a", 72, 0.5, 72.5,100, 0.7, 1.3);
0070
0071 TH2F  *h2etacoefmin16_3b = new TH2F("h2etacoefmin16_3b", "h2etacoefmin16_3b",72, 0.5, 72.5, 100, 0., 2.);
0072 TH2F  *h2etacoefmin17b = new TH2F("h2etacoefmin17b", "h2etacoefmin17b",72, 0.5, 72.5, 100, 0.6, 1.4);
0073 TH2F  *h2etacoefmin18b = new TH2F("h2etacoefmin18b", "h2etacoefmin18b",72, 0.5, 72.5, 100, 0.6, 1.4);
0074 TH2F  *h2etacoefmin19b = new TH2F("h2etacoefmin19b", "h2etacoefmin19b",72, 0.5, 72.5,100, 0.6, 1.4);
0075 TH2F  *h2etacoefmin20b = new TH2F("h2etacoefmin20b", "h2etacoefmin20b",72, 0.5, 72.5, 100, 0.6, 1.4);
0076 TH2F  *h2etacoefmin21b = new TH2F("h2etacoefmin21b", "h2etacoefmin21b",72, 0.5, 72.5, 100, 0.6, 1.4);
0077 TH2F  *h2etacoefmin22b = new TH2F("h2etacoefmin22b", "h2etacoefmin22b",72, 0.5, 72.5, 100, 0.6, 1.4);
0078 TH2F  *h2etacoefmin23b = new TH2F("h2etacoefmin23b", "h2etacoefmin23b",72, 0.5, 72.5, 100, 0.6, 1.4);
0079 TH2F  *h2etacoefmin24b = new TH2F("h2etacoefmin24b", "h2etacoefmin24b",72, 0.5, 72.5, 100, 0.6, 1.4);
0080 TH2F  *h2etacoefmin25b = new TH2F("h2etacoefmin25b", "h2etacoefmin25b",72, 0.5, 72.5, 100, 0.7, 1.3);
0081 TH2F  *h2etacoefmin26b = new TH2F("h2etacoefmin26b", "h2etacoefmin26b",72, 0.5, 72.5, 100, 0.7, 1.3);
0082 TH2F  *h2etacoefmin27b = new TH2F("h2etacoefmin27b", "h2etacoefmin27b",72, 0.5, 72.5, 100, 0.7, 1.3);
0083 TH2F  *h2etacoefmin28b = new TH2F("h2etacoefmin28b", "h2etacoefmin28b",72, 0.5, 72.5, 100, 0.7, 1.3);
0084 TH2F  *h2etacoefmin29b = new TH2F("h2etacoefmin29b", "h2etacoefmin29b",72, 0.5, 72.5, 100, 0.7, 1.3);
0085
0086 TH2F  *h2etacoefmin16_3 = new TH2F("h2etacoefmin16_3", "h2etacoefmin16_3",72, 0.5, 72.5, 100, 0., 2.);
0087 TH2F  *h2etacoefmin17 = new TH2F("h2etacoefmin17", "h2etacoefmin17",72, 0.5, 72.5, 100, 0.6, 1.4);
0088 TH2F  *h2etacoefmin18 = new TH2F("h2etacoefmin18", "h2etacoefmin18",72, 0.5, 72.5, 100, 0.6, 1.4);
0089 TH2F  *h2etacoefmin19 = new TH2F("h2etacoefmin19", "h2etacoefmin19",72, 0.5, 72.5, 100, 0.6, 1.4);
0090 TH2F  *h2etacoefmin20 = new TH2F("h2etacoefmin20", "h2etacoefmin20",72, 0.5, 72.5, 100, 0.6, 1.4);
0091 TH2F  *h2etacoefmin21 = new TH2F("h2etacoefmin21", "h2etacoefmin21",72, 0.5, 72.5, 100, 0.6, 1.4);
0092 TH2F  *h2etacoefmin22 = new TH2F("h2etacoefmin22", "h2etacoefmin22",72, 0.5, 72.5, 100, 0.6, 1.4);
0093 TH2F  *h2etacoefmin23 = new TH2F("h2etacoefmin23", "h2etacoefmin23",72, 0.5, 72.5,100, 0.6, 1.4);
0094 TH2F  *h2etacoefmin24 = new TH2F("h2etacoefmin24", "h2etacoefmin24",72, 0.5, 72.5, 100, 0.6, 1.4);
0095 TH2F  *h2etacoefmin25 = new TH2F("h2etacoefmin25", "h2etacoefmin25",72, 0.5, 72.5, 100, 0.7, 1.3);
0096 TH2F  *h2etacoefmin26 = new TH2F("h2etacoefmin26", "h2etacoefmin26",72, 0.5, 72.5, 100, 0.7, 1.3);
0097 TH2F  *h2etacoefmin27 = new TH2F("h2etacoefmin27", "h2etacoefmin27",72, 0.5, 72.5, 100, 0.7, 1.3);
0098 TH2F  *h2etacoefmin28 = new TH2F("h2etacoefmin28", "h2etacoefmin28",72, 0.5, 72.5, 100, 0.7, 1.3);
0099 TH2F  *h2etacoefmin29 = new TH2F("h2etacoefmin29", "h2etacoefmin29",72, 0.5, 72.5, 100, 0.7, 1.3);
0100
0101
0102 cout<<" Book histos "<<endl;
0103
0104 std::string line;
0105 std::ifstream in20( "coefficients_8.9mln.txt" );
0106
0107 Int_t i11 = 0;
0108
0109 Int_t maxc[36] = {1,3,5,7,9,11,13,15,17,19,21,23,25,27,29,31,33,35,37,39,41,43,45,47,49,51,53,55,57,59,61,63,65,67,69,71};
0110 Int_t maxc1[18] = {1,5,9,13,17,21,25,29,33,37,41,45,49,53,57,61,65,69};
0111
0112 while( std::getline( in20, line)){
0113 int subd,eta,phi,dep;
0114 Float_t var,err;
0115 istringstream linestream(line);
0116 linestream>>subd>>dep>>eta>>phi>>var>>err;
0117   if( subd == 2 && eta > 0 && abs(eta)<30) {
0118
0119 cout<<"eta "<<subd<<" "<<eta<<" phi "<<phi<<endl;
0120      if( phi == maxc[i11]) {
0121
0122        if(dep == 3 && eta == 16) {
0123        cout<<" 1 Phi saved "<< phi<<endl;
0124        h1etacoefmin16_3a->Fill(var);
0125        h2etacoefmin16_3a->Fill(phi,var);
0126        i11++;
0127        if(i11 == 36) i11=0;
0128        }
0129
0130        if(dep == 1 && eta == 17) {
0131        cout<<" 2 Phi saved "<< phi<<endl;
0132        h1etacoefmin17a->Fill(var);
0133        h2etacoefmin17a->Fill(phi,var);
0134        i11++;
0135        if(i11 == 36) i11=0;
0136
0137        }
0138        if(dep == 1 && eta == 18) {
0139        cout<<" 3 Phi saved "<< phi<<endl;
0140        h1etacoefmin18a->Fill(var);
0141        h2etacoefmin18a->Fill(phi,var);
0142        i11++;
0143        if(i11 == 36) i11=0;
0144
0145        }
0146        if(dep == 1 && eta == 19) {
0147        cout<<" 4 Phi saved "<< phi<<endl;
0148        h1etacoefmin19a->Fill(var);
0149        h2etacoefmin19a->Fill(phi,var);
0150        i11++;
0151        if(i11 == 36) i11=0;
0152
0153        }
0154        if(dep == 1 && eta == 20) {
0155        cout<<" 5 Phi saved "<< phi<<endl;
0156        h1etacoefmin20a->Fill(var);
0157        h2etacoefmin20a->Fill(phi,var);
0158        i11++;
0159        if(i11 == 36) i11=0;
0160        }
0161        }else{
0162
0163        if(dep == 3 && eta == 16) {
0164        h1etacoefmin16_3b->Fill(var);
0165        h2etacoefmin16_3b->Fill(phi,var);
0166        }
0167        if(dep == 1 && eta == 17) {
0168        h1etacoefmin17b->Fill(var);
0169        h2etacoefmin17b->Fill(phi,var);
0170        //if( var > 0.95 ) h1etacoefmin2bn->Fill(var);
0171        }
0172        if(dep == 1 && eta == 18) {
0173        h1etacoefmin18b->Fill(var);
0174        h2etacoefmin18b->Fill(phi,var);
0175        //if( var > 0.9 ) h1etacoefmin3bn->Fill(var);
0176        }
0177        if(dep == 1 && eta == 19) {
0178        h1etacoefmin19b->Fill(var);
0179        h2etacoefmin19b->Fill(phi,var);
0180        //if( var > 0.9 ) h1etacoefmin4bn->Fill(var);
0181        }
0182        if(dep == 1 && eta == 20) {
0183        h1etacoefmin20b->Fill(var);
0184        h2etacoefmin20b->Fill(phi,var);
0185        }
0186
0187        }
0188
0189       if(abs(eta)>20 && phi == maxc1[i11]) {
0190
0191        if(dep == 1 && eta == 21) {
0192        cout<<" 6 Phi saved "<< phi<<endl;
0193        h1etacoefmin21a->Fill(var);
0194        h2etacoefmin21a->Fill(phi,var);
0195        i11++;
0196        if(i11 == 18) i11=0;
0197
0198        }
0199        if(dep == 1 && eta == 22) {
0200        cout<<" 7 Phi saved "<< phi<<endl;
0201        h1etacoefmin22a->Fill(var);
0202        h2etacoefmin22a->Fill(phi,var);
0203        i11++;
0204        if(i11 == 18) i11=0;
0205
0206        }
0207        if(dep == 1 && eta == 23) {
0208        cout<<" 8 Phi saved "<< phi<<endl;
0209        h1etacoefmin23a->Fill(var);
0210        h2etacoefmin23a->Fill(phi,var);
0211        i11++;
0212        if(i11 == 18) i11=0;
0213
0214        }
0215        if(dep == 1 && eta == 24) {
0216        cout<<" 9 Phi saved "<< phi<<endl;
0217        h1etacoefmin24a->Fill(var);
0218        h2etacoefmin24a->Fill(phi,var);
0219        i11++;
0220        if(i11 == 18) i11=0;
0221
0222        }
0223        if(dep == 1 && eta == 25) {
0224        cout<<" 10 Phi saved "<< phi<<endl;
0225        h1etacoefmin25a->Fill(var);
0226        h2etacoefmin25a->Fill(phi,var);
0227        i11++;
0228        if(i11 == 18) i11=0;
0229
0230        }
0231        if(dep == 1 && eta == 26) {
0232        cout<<" 11 Phi saved "<< phi<<endl;
0233        h1etacoefmin26a->Fill(var);
0234        h2etacoefmin26a->Fill(phi,var);
0235        i11++;
0236        if(i11 == 18) i11=0;
0237
0238        }
0239        if(dep == 1 && eta == 27) {
0240        cout<<" 12 Phi saved "<< phi<<endl;
0241        h1etacoefmin27a->Fill(var);
0242        h2etacoefmin27a->Fill(phi,var);
0243        i11++;
0244        if(i11 == 18) i11=0;
0245
0246        }
0247        if(dep == 1 && eta == 28) {
0248        cout<<" 13 Phi saved "<< phi<<endl;
0249        h1etacoefmin28a->Fill(var);
0250        h2etacoefmin28a->Fill(phi,var);
0251        i11++;
0252        if(i11 == 18) i11=0;
0253
0254        }
0255        if(dep == 1 && eta == 29) {
0256        cout<<" 14 Phi saved "<< phi<<endl;
0257        h1etacoefmin29a->Fill(var);
0258        h2etacoefmin29a->Fill(phi,var);
0259        i11++;
0260        if(i11 == 18) i11=0;
0261        }
0262      } else{
0263         if(dep == 1 && eta == 21) {
0264        h1etacoefmin21b->Fill(var);
0265        h2etacoefmin21b->Fill(phi,var);
0266        }
0267        if(dep == 1 && eta == 22) {
0268        h1etacoefmin22b->Fill(var);
0269        h2etacoefmin22b->Fill(phi,var);
0270        }
0271        if(dep == 1 && eta == 23) {
0272        h1etacoefmin23b->Fill(var);
0273        h2etacoefmin23b->Fill(phi,var);
0274        }
0275        if(dep == 1 && eta == 24) {
0276        h1etacoefmin24b->Fill(var);
0277        h2etacoefmin24b->Fill(phi,var);
0278        }
0279        if(dep == 1 && eta == 25) {
0280        h1etacoefmin25b->Fill(var);
0281        h2etacoefmin25b->Fill(phi,var);
0282        }
0283        if(dep == 1 && eta == 26) {
0284        h1etacoefmin26b->Fill(var);
0285        h2etacoefmin26b->Fill(phi,var);
0286        }
0287        if(dep == 1 && eta == 27) {
0288        h1etacoefmin27b->Fill(var);
0289        h2etacoefmin27b->Fill(phi,var);
0290        }
0291        if(dep == 1 && eta == 28) {
0292        h1etacoefmin28b->Fill(var);
0293        h2etacoefmin28b->Fill(phi,var);
0294        }
0295        if(dep == 1 && eta == 29) {
0296        h1etacoefmin29b->Fill(var);
0297        h2etacoefmin29b->Fill(phi,var);
0298        }
0299
0300      }
0301
0302
0303     if(dep == 3 && eta == 16) {h2etacoefmin16_3->Fill(phi,var);}
0304     if(dep == 1 && eta == 17) {h2etacoefmin17->Fill(phi,var);}
0305     if(dep == 1 && eta == 18) {h2etacoefmin18->Fill(phi,var);}
0306     if(dep == 1 && eta == 19) {h2etacoefmin19->Fill(phi,var);}
0307     if(dep == 1 && eta == 20) {h2etacoefmin20->Fill(phi,var);}
0308     if(dep == 1 && eta == 21) {h2etacoefmin21->Fill(phi,var);}
0309     if(dep == 1 && eta == 22) {h2etacoefmin22->Fill(phi,var);}
0310     if(dep == 1 && eta == 23) {h2etacoefmin23->Fill(phi,var);}
0311     if(dep == 1 && eta == 24) {h2etacoefmin24->Fill(phi,var);}
0312     if(dep == 1 && eta == 25) {h2etacoefmin25->Fill(phi,var);}
0313     if(dep == 1 && eta == 26) {h2etacoefmin26->Fill(phi,var);}
0314     if(dep == 1 && eta == 27) {h2etacoefmin27->Fill(phi,var);}
0315     if(dep == 1 && eta == 28) {h2etacoefmin28->Fill(phi,var);}
0316     if(dep == 1 && eta == 29) {h2etacoefmin29->Fill(phi,var);}
0317
0318 /*
0319     if( phi == 70 || var < 0.95 ) continue;
0320 */
0321
0322     if(dep == 3 && eta == 16) {h1etacoefmin16_3->Fill(var);}
0323     if(dep == 1 && eta == 17) {h1etacoefmin17->Fill(var);}
0324     if(dep == 1 && eta == 18) {h1etacoefmin18->Fill(var);}
0325     if(dep == 1 && eta == 19) {h1etacoefmin19->Fill(var);}
0326     if(dep == 1 && eta == 20) {h1etacoefmin20->Fill(var);}
0327     if(dep == 1 && eta == 21) {h1etacoefmin21->Fill(var);}
0328     if(dep == 1 && eta == 22) {h1etacoefmin22->Fill(var);}
0329     if(dep == 1 && eta == 23) {h1etacoefmin23->Fill(var);}
0330     if(dep == 1 && eta == 24) {h1etacoefmin24->Fill(var);}
0331     if(dep == 1 && eta == 25) {h1etacoefmin25->Fill(var);}
0332     if(dep == 1 && eta == 26) {h1etacoefmin26->Fill(var);}
0333     if(dep == 1 && eta == 27) {h1etacoefmin27->Fill(var);}
0334     if(dep == 1 && eta == 28) {h1etacoefmin28->Fill(var);}
0335     if(dep == 1 && eta == 29) {h1etacoefmin29->Fill(var);}
0336
0337
0338   } // subd = 2
0339   if( subd > 2 ) break;
0340 }
0341
0342
0343 TFile efile("coefficients_219_val_endcap_plus_8.9mln.root","recreate");
0344 /*
0345 h1etacoefmin1a->Write();
0346 h1etacoefmin1a->Write();
0347 h1etacoefmin2a->Write();
0348 h1etacoefmin3a->Write();
0349 h1etacoefmin4a->Write();
0350 h1etacoefmin5a->Write();
0351 h1etacoefmin6a->Write();
0352 h1etacoefmin7a->Write();
0353 h1etacoefmin8a->Write();
0354 h1etacoefmin9a->Write();
0355 h1etacoefmin10a->Write();
0356 h1etacoefmin11a->Write();
0357 h1etacoefmin12a->Write();
0358 h1etacoefmin13a->Write();
0359 h1etacoefmin14a->Write();
0360 h1etacoefmin15a->Write();
0361 h1etacoefmin16_1a->Write();
0362 h1etacoefmin16_2a->Write();
0363
0364 h1etacoefmin1b->Write();
0365 h1etacoefmin2b->Write();h1etacoefmin2bn->Write();
0366 h1etacoefmin3b->Write();h1etacoefmin3bn->Write();
0367 h1etacoefmin4b->Write();h1etacoefmin4bn->Write();
0368 h1etacoefmin5b->Write();
0369 h1etacoefmin6b->Write();
0370 h1etacoefmin7b->Write();
0371 h1etacoefmin8b->Write();
0372 h1etacoefmin9b->Write();
0373 h1etacoefmin10b->Write();
0374 h1etacoefmin11b->Write();
0375 h1etacoefmin12b->Write();
0376 h1etacoefmin13b->Write();
0377 h1etacoefmin14b->Write();
0378 h1etacoefmin15b->Write();
0379 h1etacoefmin16_1b->Write();
0380 h1etacoefmin16_2b->Write();
0381
0382 h2etacoefmin1a->Write();
0383 h2etacoefmin1a->Write();
0384 h2etacoefmin2a->Write();
0385 h2etacoefmin3a->Write();
0386 h2etacoefmin4a->Write();
0387 h2etacoefmin5a->Write();
0388 h2etacoefmin6a->Write();
0389 h2etacoefmin7a->Write();
0390 h2etacoefmin8a->Write();
0391 h2etacoefmin9a->Write();
0392 h2etacoefmin10a->Write();
0393 h2etacoefmin11a->Write();
0394 h2etacoefmin12a->Write();
0395 h2etacoefmin13a->Write();
0396 h2etacoefmin14a->Write();
0397 h2etacoefmin15a->Write();
0398 h2etacoefmin16_1a->Write();
0399 h2etacoefmin16_2a->Write();
0400
0401 h2etacoefmin1b->Write();
0402 h2etacoefmin2b->Write();
0403 h2etacoefmin3b->Write();
0404 h2etacoefmin4b->Write();
0405 h2etacoefmin5b->Write();
0406 h2etacoefmin6b->Write();
0407 h2etacoefmin7b->Write();
0408 h2etacoefmin8b->Write();
0409 h2etacoefmin9b->Write();
0410 h2etacoefmin10b->Write();
0411 h2etacoefmin11b->Write();
0412 h2etacoefmin12b->Write();
0413 h2etacoefmin13b->Write();
0414 h2etacoefmin14b->Write();
0415 h2etacoefmin15b->Write();
0416 h2etacoefmin16_1b->Write();
0417 h2etacoefmin16_2b->Write();
0418 */
0419
0420 h2etacoefmin16_3->Write();
0421 h2etacoefmin17->Write();
0422 h2etacoefmin18->Write();
0423 h2etacoefmin19->Write();
0424 h2etacoefmin20->Write();
0425 h2etacoefmin21->Write();
0426 h2etacoefmin22->Write();
0427 h2etacoefmin23->Write();
0428 h2etacoefmin24->Write();
0429 h2etacoefmin25->Write();
0430 h2etacoefmin26->Write();
0431 h2etacoefmin27->Write();
0432 h2etacoefmin28->Write();
0433 h2etacoefmin29->Write();
0434
0435 h2etacoefmin16_3a->Write();
0436 h2etacoefmin17a->Write();
0437 h2etacoefmin18a->Write();
0438 h2etacoefmin19a->Write();
0439 h2etacoefmin20a->Write();
0440 h2etacoefmin21a->Write();
0441 h2etacoefmin22a->Write();
0442 h2etacoefmin23a->Write();
0443 h2etacoefmin24a->Write();
0444 h2etacoefmin25a->Write();
0445 h2etacoefmin26a->Write();
0446 h2etacoefmin27a->Write();
0447 h2etacoefmin28a->Write();
0448 h2etacoefmin29a->Write();
0449
0450 h2etacoefmin16_3b->Write();
0451 h2etacoefmin17b->Write();
0452 h2etacoefmin18b->Write();
0453 h2etacoefmin19b->Write();
0454 h2etacoefmin20b->Write();
0455 h2etacoefmin21b->Write();
0456 h2etacoefmin22b->Write();
0457 h2etacoefmin23b->Write();
0458 h2etacoefmin24b->Write();
0459 h2etacoefmin25b->Write();
0460 h2etacoefmin26b->Write();
0461 h2etacoefmin27b->Write();
0462 h2etacoefmin28b->Write();
0463 h2etacoefmin29b->Write();
0464
0465
0466 h1etacoefmin16_3->Write();
0467 h1etacoefmin17->Write();
0468 h1etacoefmin18->Write();
0469 h1etacoefmin19->Write();
0470 h1etacoefmin20->Write();
0471 h1etacoefmin21->Write();
0472 h1etacoefmin22->Write();
0473 h1etacoefmin23->Write();
0474 h1etacoefmin24->Write();
0475 h1etacoefmin25->Write();
0476 h1etacoefmin26->Write();
0477 h1etacoefmin27->Write();
0478 h1etacoefmin28->Write();
0479 h1etacoefmin29->Write();
0480
0481 h1etacoefmin16_3a->Write();
0482 h1etacoefmin17a->Write();
0483 h1etacoefmin18a->Write();
0484 h1etacoefmin19a->Write();
0485 h1etacoefmin20a->Write();
0486 h1etacoefmin21a->Write();
0487 h1etacoefmin22a->Write();
0488 h1etacoefmin23a->Write();
0489 h1etacoefmin24a->Write();
0490 h1etacoefmin25a->Write();
0491 h1etacoefmin26a->Write();
0492 h1etacoefmin27a->Write();
0493 h1etacoefmin28a->Write();
0494 h1etacoefmin29a->Write();
0495
0496 h1etacoefmin16_3b->Write();
0497 h1etacoefmin17b->Write();
0498 h1etacoefmin18b->Write();
0499 h1etacoefmin19b->Write();
0500 h1etacoefmin20b->Write();
0501 h1etacoefmin21b->Write();
0502 h1etacoefmin22b->Write();
0503 h1etacoefmin23b->Write();
0504 h1etacoefmin24b->Write();
0505 h1etacoefmin25b->Write();
0506 h1etacoefmin26b->Write();
0507 h1etacoefmin27b->Write();
0508 h1etacoefmin28b->Write();
0509 h1etacoefmin29b->Write();
0510
0511 }