File indexing completed on 2024-04-06 11:58:53
0001 {
0002 gROOT->Reset();
0003 gROOT->SetStyle("Plain");
0004
0005 gStyle->SetOptStat(1111);
0006 gStyle->SetOptFit(111);
0007 Float_t plmeang_nn[120][120][5][5],plnoise[120][120][5][5];
0008 Float_t minmeang_nn[120][120][5][5],minnoise[120][120][5][5];
0009
0010 FILE *Out1 = fopen("coef_with_noise.txt", "w+");
0011 FILE *Out2 = fopen("coef_without_noise_12mln.txt", "w+");
0012
0013 cout<<" Read "<<endl;
0014
0015 std::ifstream in20( "var_noise_12mln.txt" );
0016 std::string line;
0017 while( std::getline( in20, line)){
0018 istringstream linestream(line);
0019 Float_t var,err;
0020 int subd,eta,phi,dep;
0021 linestream>>subd>>eta>>phi>>dep>>var;
0022
0023
0024 if( eta > 0 )
0025 {
0026 plnoise[eta][phi][dep][subd] = var;
0027 } else
0028 {
0029 minnoise[abs(eta)][phi][dep][subd] = var;
0030 }
0031
0032 if( eta == 16) cout<<"subd "<<subd<<" "<<eta<<" "<<phi<<" "<<dep<<" "<<plnoise[eta][phi][dep][subd]<<endl;
0033
0034
0035 }
0036 cout<<" End of noise read "<<endl;
0037
0038 std::ifstream in21( "var.txt" );
0039
0040 while( std::getline( in21, line)){
0041 istringstream linestream(line);
0042 Float_t var,err;
0043 int subd,eta,phi,dep;
0044 linestream>>subd>>eta>>phi>>dep>>var>>err;
0045 if( eta > 0 )
0046 {
0047 plmeang_nn[eta][phi][dep][subd] = var;
0048 } else
0049 {
0050 minmeang_nn[abs(eta)][phi][dep][subd] = var;
0051 }
0052
0053 if( eta == 16) cout<<"subd "<<subd<<" "<<eta<<" "<<phi<<" "<<dep<<" "<<plmeang_nn[eta][phi][dep][subd]<<endl;
0054
0055 }
0056
0057
0058
0059 Int_t idep = 1;
0060
0061 Float_t plmean1[30][5][5];
0062 Float_t minmean1[30][5][5];
0063 Float_t plmean1_nn[30][5][5];
0064 Float_t minmean1_nn[30][5][5];
0065
0066 for(Int_t k=1; k<3; k++)
0067 {
0068
0069 for(Int_t i=1; i<30; i++) {
0070
0071 Int_t nch1 = 72.;
0072 if( i > 20 ) nch1 = 36.;
0073
0074 plmean1[i][idep][k] = 0.;
0075 minmean1[i][idep][k] = 0.;
0076 plmean1_nn[i][idep][k] = 0.;
0077 minmean1_nn[i][idep][k] = 0.;
0078
0079 for(Int_t j=1; j<73; j++){
0080
0081
0082
0083 plmean1[i][idep][k] = plmean1[i][idep][k] + plmeang_nn[i][j][idep][k];
0084 minmean1[i][idep][k] = minmean1[i][idep][k] + minmeang_nn[i][j][idep][k];
0085 plmean1_nn[i][idep][k] = plmean1_nn[i][idep][k] + plmeang_nn[i][j][idep][k] - plnoise[i][j][idep][k];
0086 minmean1_nn[i][idep][k] = minmean1_nn[i][idep][k] + minmeang_nn[i][j][idep][k] - minnoise[i][j][idep][k];
0087
0088 }
0089
0090 plmean1[i][idep][k] = plmean1[i][idep][k]/nch1;
0091 minmean1[i][idep][k] = minmean1[i][idep][k]/nch1;
0092
0093 plmean1_nn[i][idep][k] = plmean1_nn[i][idep][k]/nch1;
0094 minmean1_nn[i][idep][k] = minmean1_nn[i][idep][k]/nch1;
0095
0096
0097 }
0098 }
0099
0100 TH1F *h12mlnetacoefpl1 = new TH1F("h12mlnetacoefpl1", "h12mlnetacoefpl1", 100, 0., 2.);
0101 TH1F *h12mlnetacoefpl1_noise = new TH1F("h12mlnetacoefpl1_noise", "h12mlnetacoefpl1_noise", 100, 0., 2.);
0102 TH1F *h12mlnetacoefpl2 = new TH1F("h12mlnetacoefpl2", "h12mlnetacoefpl2", 100, 0., 2.);
0103 TH1F *h12mlnetacoefpl2_noise = new TH1F("h12mlnetacoefpl2_noise", "h12mlnetacoefpl2_noise", 100, 0., 2.);
0104 TH1F *h12mlnetacoefpl3 = new TH1F("h12mlnetacoefpl3", "h12mlnetacoefpl3", 100, 0., 2.);
0105 TH1F *h12mlnetacoefpl3_noise = new TH1F("h12mlnetacoefpl3_noise", "h12mlnetacoefpl3_noise", 100, 0., 2.);
0106 TH1F *h12mlnetacoefpl4 = new TH1F("h12mlnetacoefpl4", "h12mlnetacoefpl4", 100, 0., 2.);
0107 TH1F *h12mlnetacoefpl4_noise = new TH1F("h12mlnetacoefpl4_noise", "h12mlnetacoefpl4_noise", 100, 0., 2.);
0108 TH1F *h12mlnetacoefpl5 = new TH1F("h12mlnetacoefpl5", "h12mlnetacoefpl5", 100, 0., 2.);
0109 TH1F *h12mlnetacoefpl5_noise = new TH1F("h12mlnetacoefpl5_noise", "h12mlnetacoefpl5_noise", 100, 0., 2.);
0110 TH1F *h12mlnetacoefpl6 = new TH1F("h12mlnetacoefpl6", "h12mlnetacoefpl6", 100, 0., 2.);
0111 TH1F *h12mlnetacoefpl6_noise = new TH1F("h12mlnetacoefpl6_noise", "h12mlnetacoefpl6_noise", 100, 0., 2.);
0112 TH1F *h12mlnetacoefpl7 = new TH1F("h12mlnetacoefpl7", "h12mlnetacoefpl7", 100, 0., 2.);
0113 TH1F *h12mlnetacoefpl7_noise = new TH1F("h12mlnetacoefpl7_noise", "h12mlnetacoefpl7_noise", 100, 0., 2.);
0114 TH1F *h12mlnetacoefpl8 = new TH1F("h12mlnetacoefpl8", "h12mlnetacoefpl8", 100, 0., 2.);
0115 TH1F *h12mlnetacoefpl8_noise = new TH1F("h12mlnetacoefpl8_noise", "h12mlnetacoefpl8_noise", 100, 0., 2.);
0116 TH1F *h12mlnetacoefpl9 = new TH1F("h12mlnetacoefpl9", "h12mlnetacoefpl9", 100, 0., 2.);
0117 TH1F *h12mlnetacoefpl9_noise = new TH1F("h12mlnetacoefpl9_noise", "h12mlnetacoefpl9_noise", 100, 0., 2.);
0118 TH1F *h12mlnetacoefpl10 = new TH1F("h12mlnetacoefpl10", "h12mlnetacoefpl10", 100, 0., 2.);
0119 TH1F *h12mlnetacoefpl10_noise = new TH1F("h12mlnetacoefpl10_noise", "h12mlnetacoefpl10_noise", 100, 0., 2.);
0120 TH1F *h12mlnetacoefpl11 = new TH1F("h12mlnetacoefpl11", "h12mlnetacoefpl11", 100, 0., 2.);
0121 TH1F *h12mlnetacoefpl11_noise = new TH1F("h12mlnetacoefpl11_noise", "h12mlnetacoefpl11_noise", 100, 0., 2.);
0122 TH1F *h12mlnetacoefpl12 = new TH1F("h12mlnetacoefpl12", "h12mlnetacoefpl12", 100, 0., 2.);
0123 TH1F *h12mlnetacoefpl12_noise = new TH1F("h12mlnetacoefpl12_noise", "h12mlnetacoefpl12_noise", 100, 0., 2.);
0124 TH1F *h12mlnetacoefpl13 = new TH1F("h12mlnetacoefpl13", "h12mlnetacoefpl13", 100, 0., 2.);
0125 TH1F *h12mlnetacoefpl13_noise = new TH1F("h12mlnetacoefpl13_noise", "h12mlnetacoefpl13_noise", 100, 0., 2.);
0126 TH1F *h12mlnetacoefpl14 = new TH1F("h12mlnetacoefpl14", "h12mlnetacoefpl14", 100, 0., 2.);
0127 TH1F *h12mlnetacoefpl14_noise = new TH1F("h12mlnetacoefpl14_noise", "h12mlnetacoefpl14_noise", 100, 0., 2.);
0128 TH1F *h12mlnetacoefpl15 = new TH1F("h12mlnetacoefpl15", "h12mlnetacoefpl15", 100, 0., 2.);
0129 TH1F *h12mlnetacoefpl15_noise = new TH1F("h12mlnetacoefpl15_noise", "h12mlnetacoefpl15_noise", 100, 0., 2.);
0130
0131 TH1F *h12mlnetacoefpl16_HB = new TH1F("h12mlnetacoefpl16_HB", "h12mlnetacoefpl16_HB", 100, 0., 2.);
0132 TH1F *h12mlnetacoefpl16_HB_noise = new TH1F("h12mlnetacoefpl16_HB_noise", "h12mlnetacoefpl16_HB_noise", 100, 0., 2.);
0133 TH1F *h12mlnetacoefpl16_HE = new TH1F("h12mlnetacoefpl16_HE", "h12mlnetacoefpl16_HE", 100, 0., 2.);
0134 TH1F *h12mlnetacoefpl16_HE_noise = new TH1F("h12mlnetacoefpl16_HE_noise", "h12mlnetacoefpl16_HE_noise", 100, 0., 2.);
0135
0136 TH1F *h12mlnetacoefpl17 = new TH1F("h12mlnetacoefpl17", "h12mlnetacoefpl17", 100, 0., 2.);
0137 TH1F *h12mlnetacoefpl17_noise = new TH1F("h12mlnetacoefpl17_noise", "h12mlnetacoefpl17_noise", 100, 0., 2.);
0138 TH1F *h12mlnetacoefpl18 = new TH1F("h12mlnetacoefpl18", "h12mlnetacoefpl18", 100, 0., 2.);
0139 TH1F *h12mlnetacoefpl18_noise = new TH1F("h12mlnetacoefpl18_noise", "h12mlnetacoefpl18_noise", 100, 0., 2.);
0140 TH1F *h12mlnetacoefpl19 = new TH1F("h12mlnetacoefpl19", "h12mlnetacoefpl19", 100, 0., 2.);
0141 TH1F *h12mlnetacoefpl19_noise = new TH1F("h12mlnetacoefpl19_noise", "h12mlnetacoefpl19_noise", 100, 0., 2.);
0142 TH1F *h12mlnetacoefpl20 = new TH1F("h12mlnetacoefpl20", "h12mlnetacoefpl20", 100, 0., 2.);
0143 TH1F *h12mlnetacoefpl20_noise = new TH1F("h12mlnetacoefpl20_noise", "h12mlnetacoefpl20_noise", 100, 0., 2.);
0144 TH1F *h12mlnetacoefpl21 = new TH1F("h12mlnetacoefpl21", "h12mlnetacoefpl21", 100, 0., 2.);
0145 TH1F *h12mlnetacoefpl21_noise = new TH1F("h12mlnetacoefpl21_noise", "h12mlnetacoefpl21_noise", 100, 0., 2.);
0146 TH1F *h12mlnetacoefpl22 = new TH1F("h12mlnetacoefpl22", "h12mlnetacoefpl22", 100, 0., 2.);
0147 TH1F *h12mlnetacoefpl22_noise = new TH1F("h12mlnetacoefpl22_noise", "h12mlnetacoefpl22_noise", 100, 0., 2.);
0148 TH1F *h12mlnetacoefpl23 = new TH1F("h12mlnetacoefpl23", "h12mlnetacoefpl23", 100, 0., 2.);
0149 TH1F *h12mlnetacoefpl23_noise = new TH1F("h12mlnetacoefpl23_noise", "h12mlnetacoefpl23_noise", 100, 0., 2.);
0150 TH1F *h12mlnetacoefpl24 = new TH1F("h12mlnetacoefpl24", "h12mlnetacoefpl24", 100, 0., 2.);
0151 TH1F *h12mlnetacoefpl24_noise = new TH1F("h12mlnetacoefpl24_noise", "h12mlnetacoefpl24_noise", 100, 0., 2.);
0152 TH1F *h12mlnetacoefpl25 = new TH1F("h12mlnetacoefpl25", "h12mlnetacoefpl25", 100, 0., 2.);
0153 TH1F *h12mlnetacoefpl25_noise = new TH1F("h12mlnetacoefpl25_noise", "h12mlnetacoefpl25_noise", 100, 0., 2.);
0154 TH1F *h12mlnetacoefpl26 = new TH1F("h12mlnetacoefpl26", "h12mlnetacoefpl26", 100, 0., 2.);
0155 TH1F *h12mlnetacoefpl26_noise = new TH1F("h12mlnetacoefpl26_noise", "h12mlnetacoefpl26_noise", 100, 0., 2.);
0156 TH1F *h12mlnetacoefpl27 = new TH1F("h12mlnetacoefpl27", "h12mlnetacoefpl27", 100, 0., 2.);
0157 TH1F *h12mlnetacoefpl27_noise = new TH1F("h12mlnetacoefpl27_noise", "h12mlnetacoefpl27_noise", 100, 0., 2.);
0158 TH1F *h12mlnetacoefpl28 = new TH1F("h12mlnetacoefpl28", "h12mlnetacoefpl28", 100, 0., 2.);
0159 TH1F *h12mlnetacoefpl28_noise = new TH1F("h12mlnetacoefpl28_noise", "h12mlnetacoefpl28_noise", 100, 0., 2.);
0160 TH1F *h12mlnetacoefpl29 = new TH1F("h12mlnetacoefpl29", "h12mlnetacoefpl29", 100, 0., 2.);
0161 TH1F *h12mlnetacoefpl29_noise = new TH1F("h12mlnetacoefpl29_noise", "h12mlnetacoefpl29_noise", 100, 0., 2.);
0162
0163 TH1F *h12mlnetacoefmin1 = new TH1F("h12mlnetacoefmin1", "h12mlnetacoefmin1", 100, 0., 2.);
0164 TH1F *h12mlnetacoefmin1_noise = new TH1F("h12mlnetacoefmin1_noise", "h12mlnetacoefmin1_noise", 100, 0., 2.);
0165 TH1F *h12mlnetacoefmin2 = new TH1F("h12mlnetacoefmin2", "h12mlnetacoefmin2", 100, 0., 2.);
0166 TH1F *h12mlnetacoefmin2_noise = new TH1F("h12mlnetacoefmin2_noise", "h12mlnetacoefmin2_noise", 100, 0., 2.);
0167 TH1F *h12mlnetacoefmin3 = new TH1F("h12mlnetacoefmin3", "h12mlnetacoefmin3", 100, 0., 2.);
0168 TH1F *h12mlnetacoefmin3_noise = new TH1F("h12mlnetacoefmin3_noise", "h12mlnetacoefmin3_noise", 100, 0., 2.);
0169 TH1F *h12mlnetacoefmin4 = new TH1F("h12mlnetacoefmin4", "h12mlnetacoefmin4", 100, 0., 2.);
0170 TH1F *h12mlnetacoefmin4_noise = new TH1F("h12mlnetacoefmin4_noise", "h12mlnetacoefmin4_noise", 100, 0., 2.);
0171 TH1F *h12mlnetacoefmin5 = new TH1F("h12mlnetacoefmin5", "h12mlnetacoefmin5", 100, 0., 2.);
0172 TH1F *h12mlnetacoefmin5_noise = new TH1F("h12mlnetacoefmin5_noise", "h12mlnetacoefmin5_noise", 100, 0., 2.);
0173 TH1F *h12mlnetacoefmin6 = new TH1F("h12mlnetacoefmin6", "h12mlnetacoefmin6", 100, 0., 2.);
0174 TH1F *h12mlnetacoefmin6_noise = new TH1F("h12mlnetacoefmin6_noise", "h12mlnetacoefmin6_noise", 100, 0., 2.);
0175 TH1F *h12mlnetacoefmin7 = new TH1F("h12mlnetacoefmin7", "h12mlnetacoefmin7", 100, 0., 2.);
0176 TH1F *h12mlnetacoefmin7_noise = new TH1F("h12mlnetacoefmin7_noise", "h12mlnetacoefmin7_noise", 100, 0., 2.);
0177 TH1F *h12mlnetacoefmin8 = new TH1F("h12mlnetacoefmin8", "h12mlnetacoefmin8", 100, 0., 2.);
0178 TH1F *h12mlnetacoefmin8_noise = new TH1F("h12mlnetacoefmin8_noise", "h12mlnetacoefmin8_noise", 100, 0., 2.);
0179 TH1F *h12mlnetacoefmin9 = new TH1F("h12mlnetacoefmin9", "h12mlnetacoefmin9", 100, 0., 2.);
0180 TH1F *h12mlnetacoefmin9_noise = new TH1F("h12mlnetacoefmin9_noise", "h12mlnetacoefmin9_noise", 100, 0., 2.);
0181 TH1F *h12mlnetacoefmin10 = new TH1F("h12mlnetacoefmin10", "h12mlnetacoefmin10", 100, 0., 2.);
0182 TH1F *h12mlnetacoefmin10_noise = new TH1F("h12mlnetacoefmin10_noise", "h12mlnetacoefmin10_noise", 100, 0., 2.);
0183 TH1F *h12mlnetacoefmin11 = new TH1F("h12mlnetacoefmin11", "h12mlnetacoefmin11", 100, 0., 2.);
0184 TH1F *h12mlnetacoefmin11_noise = new TH1F("h12mlnetacoefmin11_noise", "h12mlnetacoefmin11_noise", 100, 0., 2.);
0185 TH1F *h12mlnetacoefmin12 = new TH1F("h12mlnetacoefmin12", "h12mlnetacoefmin12", 100, 0., 2.);
0186 TH1F *h12mlnetacoefmin12_noise = new TH1F("h12mlnetacoefmin12_noise", "h12mlnetacoefmin12_noise", 100, 0., 2.);
0187 TH1F *h12mlnetacoefmin13 = new TH1F("h12mlnetacoefmin13", "h12mlnetacoefmin13", 100, 0., 2.);
0188 TH1F *h12mlnetacoefmin13_noise = new TH1F("h12mlnetacoefmin13_noise", "h12mlnetacoefmin13_noise", 100, 0., 2.);
0189 TH1F *h12mlnetacoefmin14 = new TH1F("h12mlnetacoefmin14", "h12mlnetacoefmin14", 100, 0., 2.);
0190 TH1F *h12mlnetacoefmin14_noise = new TH1F("h12mlnetacoefmin14_noise", "h12mlnetacoefmin14_noise", 100, 0., 2.);
0191 TH1F *h12mlnetacoefmin15 = new TH1F("h12mlnetacoefmin15", "h12mlnetacoefmin15", 100, 0., 2.);
0192 TH1F *h12mlnetacoefmin15_noise = new TH1F("h12mlnetacoefmin15_noise", "h12mlnetacoefmin15_noise", 100, 0., 2.);
0193
0194 TH1F *h12mlnetacoefmin16_HB = new TH1F("h12mlnetacoefmin16_HB", "h12mlnetacoefmin16_HB", 100, 0., 2.);
0195 TH1F *h12mlnetacoefmin16_HB_noise = new TH1F("h12mlnetacoefmin16_HB_noise", "h12mlnetacoefmin16_HB_noise", 100, 0., 2.);
0196 TH1F *h12mlnetacoefmin16_HE = new TH1F("h12mlnetacoefmin16_HE", "h12mlnetacoefmin16_HE", 100, 0., 2.);
0197 TH1F *h12mlnetacoefmin16_HE_noise = new TH1F("h12mlnetacoefmin16_HE_noise", "h12mlnetacoefmin16_HE_noise", 100, 0., 2.);
0198
0199 TH1F *h12mlnetacoefmin17 = new TH1F("h12mlnetacoefmin17", "h12mlnetacoefmin17", 100, 0., 2.);
0200 TH1F *h12mlnetacoefmin17_noise = new TH1F("h12mlnetacoefmin17_noise", "h12mlnetacoefmin17_noise", 100, 0., 2.);
0201 TH1F *h12mlnetacoefmin18 = new TH1F("h12mlnetacoefmin18", "h12mlnetacoefmin18", 100, 0., 2.);
0202 TH1F *h12mlnetacoefmin18_noise = new TH1F("h12mlnetacoefmin18_noise", "h12mlnetacoefmin18_noise", 100, 0., 2.);
0203 TH1F *h12mlnetacoefmin19 = new TH1F("h12mlnetacoefmin19", "h12mlnetacoefmin19", 100, 0., 2.);
0204 TH1F *h12mlnetacoefmin19_noise = new TH1F("h12mlnetacoefmin19_noise", "h12mlnetacoefmin19_noise", 100, 0., 2.);
0205 TH1F *h12mlnetacoefmin20 = new TH1F("h12mlnetacoefmin20", "h12mlnetacoefmin20", 100, 0., 2.);
0206 TH1F *h12mlnetacoefmin20_noise = new TH1F("h12mlnetacoefmin20_noise", "h12mlnetacoefmin20_noise", 100, 0., 2.);
0207 TH1F *h12mlnetacoefmin21 = new TH1F("h12mlnetacoefmin21", "h12mlnetacoefmin21", 100, 0., 2.);
0208 TH1F *h12mlnetacoefmin21_noise = new TH1F("h12mlnetacoefmin21_noise", "h12mlnetacoefmin21_noise", 100, 0., 2.);
0209 TH1F *h12mlnetacoefmin22 = new TH1F("h12mlnetacoefmin22", "h12mlnetacoefmin22", 100, 0., 2.);
0210 TH1F *h12mlnetacoefmin22_noise = new TH1F("h12mlnetacoefmin22_noise", "h12mlnetacoefmin22_noise", 100, 0., 2.);
0211 TH1F *h12mlnetacoefmin23 = new TH1F("h12mlnetacoefmin23", "h12mlnetacoefmin23", 100, 0., 2.);
0212 TH1F *h12mlnetacoefmin23_noise = new TH1F("h12mlnetacoefmin23_noise", "h12mlnetacoefmin23_noise", 100, 0., 2.);
0213 TH1F *h12mlnetacoefmin24 = new TH1F("h12mlnetacoefmin24", "h12mlnetacoefmin24", 100, 0., 2.);
0214 TH1F *h12mlnetacoefmin24_noise = new TH1F("h12mlnetacoefmin24_noise", "h12mlnetacoefmin24_noise", 100, 0., 2.);
0215 TH1F *h12mlnetacoefmin25 = new TH1F("h12mlnetacoefmin25", "h12mlnetacoefmin25", 100, 0., 2.);
0216 TH1F *h12mlnetacoefmin25_noise = new TH1F("h12mlnetacoefmin25_noise", "h12mlnetacoefmin25_noise", 100, 0., 2.);
0217 TH1F *h12mlnetacoefmin26 = new TH1F("h12mlnetacoefmin26", "h12mlnetacoefmin26", 100, 0., 2.);
0218 TH1F *h12mlnetacoefmin26_noise = new TH1F("h12mlnetacoefmin26_noise", "h12mlnetacoefmin26_noise", 100, 0., 2.);
0219 TH1F *h12mlnetacoefmin27 = new TH1F("h12mlnetacoefmin27", "h12mlnetacoefmin27", 100, 0., 2.);
0220 TH1F *h12mlnetacoefmin27_noise = new TH1F("h12mlnetacoefmin27_noise", "h12mlnetacoefmin27_noise", 100, 0., 2.);
0221 TH1F *h12mlnetacoefmin28 = new TH1F("h12mlnetacoefmin28", "h12mlnetacoefmin28", 100, 0., 2.);
0222 TH1F *h12mlnetacoefmin28_noise = new TH1F("h12mlnetacoefmin28_noise", "h12mlnetacoefmin28_noise", 100, 0., 2.);
0223 TH1F *h12mlnetacoefmin29 = new TH1F("h12mlnetacoefmin29", "h12mlnetacoefmin29", 100, 0., 2.);
0224 TH1F *h12mlnetacoefmin29_noise = new TH1F("h12mlnetacoefmin29_noise", "h12mlnetacoefmin29_noise", 100, 0., 2.);
0225
0226
0227
0228
0229 TH2F *h12mlnetacoefpl2D1 = new TH2F("h12mlnetacoefpl2D1", "h12mlnetacoefpl2D1", 72, 0.5, 72.5, 100, 0., 2.);
0230 TH2F *h12mlnetacoefpl2D1_noise = new TH2F("h12mlnetacoefpl2D1_noise", "h12mlnetacoefpl2D1_noise",72, 0.5, 72.5, 100, 0., 2.);
0231 TH2F *h12mlnetacoefmin2D1 = new TH2F("h12mlnetacoefmin2D1", "h12mlnetacoefmin2D1",72, 0.5, 72.5, 100, 0., 2.);
0232 TH2F *h12mlnetacoefmin2D1_noise = new TH2F("h12mlnetacoefmin2D1_noise", "h12mlnetacoefmin2D1_noise",72, 0.5, 72.5, 100, 0., 2.);
0233
0234
0235 TH2F *h12mlnetacoefpl2D16 = new TH2F("h12mlnetacoefpl2D16", "h12mlnetacoefpl2D16", 72, 0.5, 72.5, 100, 0., 2.);
0236 TH2F *h12mlnetacoefpl2D16_noise = new TH2F("h12mlnetacoefpl2D16_noise", "h12mlnetacoefpl2D16_noise",72, 0.5, 72.5, 100, 0., 2.);
0237 TH2F *h12mlnetacoefmin2D16 = new TH2F("h12mlnetacoefmin2D16", "h12mlnetacoefmin2D16",72, 0.5, 72.5, 100, 0., 2.);
0238 TH2F *h12mlnetacoefmin2D16_noise = new TH2F("h12mlnetacoefmin2D16_noise", "h12mlnetacoefmin2D16_noise",72, 0.5, 72.5, 100, 0., 2.);
0239
0240
0241
0242
0243
0244 for(Int_t k=1; k<3; k++)
0245 {
0246 for (Int_t i=1; i<30; i++)
0247 {
0248 for (Int_t j=1; j<73; j++)
0249 {
0250
0251
0252
0253
0254
0255 if (plmeang_nn[i][j][idep][k]-plnoise[i][j][idep][k] > 0) {
0256 Float_t tt0 = sqrt(plmean1_nn[i][idep][k]/(plmeang_nn[i][j][idep][k] - plnoise[i][j][idep][k]));
0257
0258
0259
0260 if( i == 1 ) {
0261
0262 h12mlnetacoefpl1->Fill(tt0);
0263 h12mlnetacoefpl2D1->Fill((float)j, tt0);
0264 }
0265
0266 if( i == 2 ) h12mlnetacoefpl2->Fill(tt0);
0267 if( i == 3 ) h12mlnetacoefpl3->Fill(tt0);
0268 if( i == 4 ) h12mlnetacoefpl4->Fill(tt0);
0269 if( i == 5 ) h12mlnetacoefpl5->Fill(tt0);
0270 if( i == 6 ) h12mlnetacoefpl6->Fill(tt0);
0271 if( i == 7 ) h12mlnetacoefpl7->Fill(tt0);
0272 if( i == 8 ) h12mlnetacoefpl8->Fill(tt0);
0273 if( i == 9 ) h12mlnetacoefpl9->Fill(tt0);
0274 if( i == 10 ) h12mlnetacoefpl10->Fill(tt0);
0275 if( i == 11 ) {
0276 cout<<" J= "<<j<<" Mean "<<plmean1_nn[i][idep][k]<<" Value "<<
0277 plmeang_nn[i][j][idep][k] - plnoise[i][j][idep][k]<<" tt0 "<<tt0<<endl;
0278 h12mlnetacoefpl11->Fill(tt0);
0279 }
0280 if( i == 12 ) h12mlnetacoefpl12->Fill(tt0);
0281 if( i == 13 ) h12mlnetacoefpl13->Fill(tt0);
0282 if( i == 14 ) h12mlnetacoefpl14->Fill(tt0);
0283 if( i == 15 ) h12mlnetacoefpl15->Fill(tt0);
0284 if( i == 16 ) { if (k == 1) h12mlnetacoefpl16_HB->Fill(tt0); if( k == 2 ) h12mlnetacoefpl16_HE->Fill(tt0);}
0285 if( i == 16 ) h12mlnetacoefpl2D16->Fill((float)j, tt0);
0286 if( i == 17 ) h12mlnetacoefpl17->Fill(tt0);
0287 if( i == 18 ) h12mlnetacoefpl18->Fill(tt0);
0288 if( i == 19 ) h12mlnetacoefpl19->Fill(tt0);
0289 if( i == 20 ) h12mlnetacoefpl20->Fill(tt0);
0290 if( i == 21 ) h12mlnetacoefpl21->Fill(tt0);
0291 if( i == 22 ) h12mlnetacoefpl22->Fill(tt0);
0292 if( i == 23 ) h12mlnetacoefpl23->Fill(tt0);
0293 if( i == 24 ) h12mlnetacoefpl24->Fill(tt0);
0294 if( i == 25 ) h12mlnetacoefpl25->Fill(tt0);
0295 if( i == 26 ) h12mlnetacoefpl26->Fill(tt0);
0296 if( i == 27 ) h12mlnetacoefpl27->Fill(tt0);
0297 if( i == 28 ) h12mlnetacoefpl28->Fill(tt0);
0298 if( i == 29 ) h12mlnetacoefpl29->Fill(tt0);
0299
0300 fprintf(Out2,"%s %d %d %d %.5f\n","HB",i,j,idep,tt0);
0301
0302 }
0303
0304 if (plmeang_nn[i][j][idep][k] > 0) {
0305 if( i == 1 ) {
0306
0307 }
0308
0309 Float_t tt0 = sqrt(plmean1[i][idep][k]/plmeang_nn[i][j][idep][k]);
0310 if( i == 1 ) h12mlnetacoefpl1_noise->Fill(tt0);
0311 if( i == 2 ) h12mlnetacoefpl2_noise->Fill(tt0);
0312 if( i == 3 ) h12mlnetacoefpl3_noise->Fill(tt0);
0313 if( i == 4 ) h12mlnetacoefpl4_noise->Fill(tt0);
0314 if( i == 5 ) h12mlnetacoefpl5_noise->Fill(tt0);
0315 if( i == 6 ) h12mlnetacoefpl6_noise->Fill(tt0);
0316 if( i == 7 ) h12mlnetacoefpl7_noise->Fill(tt0);
0317 if( i == 8 ) h12mlnetacoefpl8_noise->Fill(tt0);
0318 if( i == 9 ) h12mlnetacoefpl9_noise->Fill(tt0);
0319 if( i == 10 ) h12mlnetacoefpl10_noise->Fill(tt0);
0320 if( i == 11 ) h12mlnetacoefpl11_noise->Fill(tt0);
0321 if( i == 12 ) h12mlnetacoefpl12_noise->Fill(tt0);
0322 if( i == 13 ) h12mlnetacoefpl13_noise->Fill(tt0);
0323 if( i == 14 ) h12mlnetacoefpl14_noise->Fill(tt0);
0324 if( i == 15 ) h12mlnetacoefpl15_noise->Fill(tt0);
0325 if( i == 16 ) {if (k == 1) h12mlnetacoefpl16_HB_noise->Fill(tt0); if( k == 2 ) h12mlnetacoefpl16_HE_noise->Fill(tt0);}
0326 if( i == 17 ) h12mlnetacoefpl17_noise->Fill(tt0);
0327 if( i == 18 ) h12mlnetacoefpl18_noise->Fill(tt0);
0328 if( i == 19 ) h12mlnetacoefpl19_noise->Fill(tt0);
0329 if( i == 20 ) h12mlnetacoefpl20_noise->Fill(tt0);
0330 if( i == 21 ) h12mlnetacoefpl21_noise->Fill(tt0);
0331 if( i == 22 ) h12mlnetacoefpl22_noise->Fill(tt0);
0332 if( i == 23 ) h12mlnetacoefpl23_noise->Fill(tt0);
0333 if( i == 24 ) h12mlnetacoefpl24_noise->Fill(tt0);
0334 if( i == 25 ) h12mlnetacoefpl25_noise->Fill(tt0);
0335 if( i == 26 ) h12mlnetacoefpl26_noise->Fill(tt0);
0336 if( i == 27 ) h12mlnetacoefpl27_noise->Fill(tt0);
0337 if( i == 28 ) h12mlnetacoefpl28_noise->Fill(tt0);
0338 if( i == 29 ) h12mlnetacoefpl29_noise->Fill(tt0);
0339
0340 if( i == 1 ) h12mlnetacoefpl2D1_noise->Fill((float)j, tt0);
0341 if( i == 16 ) h12mlnetacoefpl2D16_noise->Fill((float)j, tt0);
0342
0343 fprintf(Out1,"%s %d %d %d %.5f\n","HB",i,j,idep,tt0);
0344 if( i == 1 ) {
0345
0346 }
0347
0348 }
0349
0350 Int_t ii = -1*i;
0351
0352 if (minmeang_nn[i][j][idep][k]-minnoise[i][j][idep][k] > 0) {
0353 if( i == 1 ) {
0354
0355 }
0356
0357 Float_t tt0 = sqrt(minmean1_nn[i][idep][k]/(minmeang_nn[i][j][idep][k] - minnoise[i][j][idep][k]));
0358
0359 if( i == 1 ) h12mlnetacoefmin1->Fill(tt0);
0360 if( i == 2 ) h12mlnetacoefmin2->Fill(tt0);
0361 if( i == 3 ) h12mlnetacoefmin3->Fill(tt0);
0362 if( i == 4 ) h12mlnetacoefmin4->Fill(tt0);
0363 if( i == 5 ) h12mlnetacoefmin5->Fill(tt0);
0364 if( i == 6 ) h12mlnetacoefmin6->Fill(tt0);
0365 if( i == 7 ) h12mlnetacoefmin7->Fill(tt0);
0366 if( i == 8 ) h12mlnetacoefmin8->Fill(tt0);
0367 if( i == 9 ) h12mlnetacoefmin9->Fill(tt0);
0368 if( i == 10 ) h12mlnetacoefmin10->Fill(tt0);
0369 if( i == 11 ) h12mlnetacoefmin11->Fill(tt0);
0370 if( i == 12 ) h12mlnetacoefmin12->Fill(tt0);
0371 if( i == 13 ) h12mlnetacoefmin13->Fill(tt0);
0372 if( i == 14 ) h12mlnetacoefmin14->Fill(tt0);
0373 if( i == 15 ) h12mlnetacoefmin15->Fill(tt0);
0374 if( i == 16 ) { if (k == 1) h12mlnetacoefmin16_HB->Fill(tt0); if( k == 2 ) h12mlnetacoefmin16_HE->Fill(tt0);}
0375 if( i == 17 ) h12mlnetacoefmin17->Fill(tt0);
0376 if( i == 18 ) h12mlnetacoefmin18->Fill(tt0);
0377 if( i == 19 ) h12mlnetacoefmin19->Fill(tt0);
0378 if( i == 20 ) h12mlnetacoefmin20->Fill(tt0);
0379 if( i == 21 ) h12mlnetacoefmin21->Fill(tt0);
0380 if( i == 22 ) h12mlnetacoefmin22->Fill(tt0);
0381 if( i == 23 ) h12mlnetacoefmin23->Fill(tt0);
0382 if( i == 24 ) h12mlnetacoefmin24->Fill(tt0);
0383 if( i == 25 ) h12mlnetacoefmin25->Fill(tt0);
0384 if( i == 26 ) h12mlnetacoefmin26->Fill(tt0);
0385 if( i == 27 ) h12mlnetacoefmin27->Fill(tt0);
0386 if( i == 28 ) h12mlnetacoefmin28->Fill(tt0);
0387 if( i == 29 ) h12mlnetacoefmin29->Fill(tt0);
0388
0389 fprintf(Out2,"%s %d %d %d %.5f\n","HB",ii,j,idep,tt0);
0390 if( i == 1 ) h12mlnetacoefmin2D1->Fill((float)j, tt0);
0391 if( i == 16 ) h12mlnetacoefmin2D16->Fill((float)j, tt0);
0392
0393 }
0394
0395 if (minmeang_nn[i][j][idep][k] > 0) {
0396 if( i == 1 ) {
0397
0398 }
0399
0400 Float_t tt0 = sqrt(minmean1[i][idep][k]/minmeang_nn[i][j][idep][k]);
0401 if( i == 1 ) h12mlnetacoefmin1_noise->Fill(tt0);
0402 if( i == 2 ) h12mlnetacoefmin2_noise->Fill(tt0);
0403 if( i == 3 ) h12mlnetacoefmin3_noise->Fill(tt0);
0404 if( i == 4 ) h12mlnetacoefmin4_noise->Fill(tt0);
0405 if( i == 5 ) h12mlnetacoefmin5_noise->Fill(tt0);
0406 if( i == 6 ) h12mlnetacoefmin6_noise->Fill(tt0);
0407 if( i == 7 ) h12mlnetacoefmin7_noise->Fill(tt0);
0408 if( i == 8 ) h12mlnetacoefmin8_noise->Fill(tt0);
0409 if( i == 9 ) h12mlnetacoefmin9_noise->Fill(tt0);
0410 if( i == 10 ) h12mlnetacoefmin10_noise->Fill(tt0);
0411 if( i == 11 ) h12mlnetacoefmin11_noise->Fill(tt0);
0412 if( i == 12 ) h12mlnetacoefmin12_noise->Fill(tt0);
0413 if( i == 13 ) h12mlnetacoefmin13_noise->Fill(tt0);
0414 if( i == 14 ) h12mlnetacoefmin14_noise->Fill(tt0);
0415 if( i == 15 ) h12mlnetacoefmin15_noise->Fill(tt0);
0416 if( i == 16 ) { if (k == 1) h12mlnetacoefmin16_HB_noise->Fill(tt0); if( k == 2 ) h12mlnetacoefmin16_HE_noise->Fill(tt0);}
0417 if( i == 17 ) h12mlnetacoefmin17_noise->Fill(tt0);
0418 if( i == 18 ) h12mlnetacoefmin18_noise->Fill(tt0);
0419 if( i == 19 ) h12mlnetacoefmin19_noise->Fill(tt0);
0420 if( i == 20 ) h12mlnetacoefmin20_noise->Fill(tt0);
0421 if( i == 21 ) h12mlnetacoefmin21_noise->Fill(tt0);
0422 if( i == 22 ) h12mlnetacoefmin22_noise->Fill(tt0);
0423 if( i == 23 ) h12mlnetacoefmin23_noise->Fill(tt0);
0424 if( i == 24 ) h12mlnetacoefmin24_noise->Fill(tt0);
0425 if( i == 25 ) h12mlnetacoefmin25_noise->Fill(tt0);
0426 if( i == 26 ) h12mlnetacoefmin26_noise->Fill(tt0);
0427 if( i == 27 ) h12mlnetacoefmin27_noise->Fill(tt0);
0428 if( i == 28 ) h12mlnetacoefmin28_noise->Fill(tt0);
0429 if( i == 29 ) h12mlnetacoefmin29_noise->Fill(tt0);
0430
0431 fprintf(Out1,"%s %d %d %d %.5f\n","HB",ii,j,idep,tt0);
0432 if( i == 1 ) h12mlnetacoefmin2D1_noise->Fill((float)j, tt0);
0433 if( i == 16 ) h12mlnetacoefmin2D16_noise->Fill((float)j, tt0);
0434
0435 }
0436
0437
0438 }
0439
0440 }
0441 }
0442
0443 fclose(Out1);
0444 fclose(Out2);
0445
0446
0447
0448
0449
0450 TFile efile("coefficients_12mln.root","recreate");
0451
0452 h12mlnetacoefpl1_noise->Write();
0453 h12mlnetacoefpl2_noise->Write();
0454 h12mlnetacoefpl3_noise->Write();
0455 h12mlnetacoefpl4_noise->Write();
0456 h12mlnetacoefpl5_noise->Write();
0457 h12mlnetacoefpl6_noise->Write();
0458 h12mlnetacoefpl7_noise->Write();
0459 h12mlnetacoefpl8_noise->Write();
0460 h12mlnetacoefpl9_noise->Write();
0461 h12mlnetacoefpl10_noise->Write();
0462 h12mlnetacoefpl11_noise->Write();
0463 h12mlnetacoefpl12_noise->Write();
0464 h12mlnetacoefpl13_noise->Write();
0465 h12mlnetacoefpl14_noise->Write();
0466 h12mlnetacoefpl15_noise->Write();
0467 h12mlnetacoefpl16_HB_noise->Write();
0468 h12mlnetacoefpl16_HE_noise->Write();
0469 h12mlnetacoefpl17_noise->Write();
0470 h12mlnetacoefpl18_noise->Write();
0471 h12mlnetacoefpl19_noise->Write();
0472 h12mlnetacoefpl20_noise->Write();
0473 h12mlnetacoefpl21_noise->Write();
0474 h12mlnetacoefpl22_noise->Write();
0475 h12mlnetacoefpl23_noise->Write();
0476 h12mlnetacoefpl24_noise->Write();
0477 h12mlnetacoefpl25_noise->Write();
0478 h12mlnetacoefpl26_noise->Write();
0479 h12mlnetacoefpl27_noise->Write();
0480 h12mlnetacoefpl28_noise->Write();
0481 h12mlnetacoefpl29_noise->Write();
0482
0483 h12mlnetacoefpl1->Write();
0484 h12mlnetacoefpl2->Write();
0485 h12mlnetacoefpl3->Write();
0486 h12mlnetacoefpl4->Write();
0487 h12mlnetacoefpl5->Write();
0488 h12mlnetacoefpl6->Write();
0489 h12mlnetacoefpl7->Write();
0490 h12mlnetacoefpl8->Write();
0491 h12mlnetacoefpl9->Write();
0492 h12mlnetacoefpl10->Write();
0493 h12mlnetacoefpl11->Write();
0494 h12mlnetacoefpl12->Write();
0495 h12mlnetacoefpl13->Write();
0496 h12mlnetacoefpl14->Write();
0497 h12mlnetacoefpl15->Write();
0498 h12mlnetacoefpl16_HB->Write();
0499 h12mlnetacoefpl16_HE->Write();
0500 h12mlnetacoefpl17->Write();
0501 h12mlnetacoefpl18->Write();
0502 h12mlnetacoefpl19->Write();
0503 h12mlnetacoefpl20->Write();
0504 h12mlnetacoefpl21->Write();
0505 h12mlnetacoefpl22->Write();
0506 h12mlnetacoefpl23->Write();
0507 h12mlnetacoefpl24->Write();
0508 h12mlnetacoefpl25->Write();
0509 h12mlnetacoefpl26->Write();
0510 h12mlnetacoefpl27->Write();
0511 h12mlnetacoefpl28->Write();
0512 h12mlnetacoefpl29->Write();
0513
0514 h12mlnetacoefpl1_noise->Clear();
0515 h12mlnetacoefpl2_noise->Clear();
0516 h12mlnetacoefpl3_noise->Clear();
0517 h12mlnetacoefpl4_noise->Clear();
0518 h12mlnetacoefpl5_noise->Clear();
0519 h12mlnetacoefpl6_noise->Clear();
0520 h12mlnetacoefpl7_noise->Clear();
0521 h12mlnetacoefpl8_noise->Clear();
0522 h12mlnetacoefpl9_noise->Clear();
0523 h12mlnetacoefpl10_noise->Clear();
0524 h12mlnetacoefpl11_noise->Clear();
0525 h12mlnetacoefpl12_noise->Clear();
0526 h12mlnetacoefpl13_noise->Clear();
0527 h12mlnetacoefpl14_noise->Clear();
0528 h12mlnetacoefpl15_noise->Clear();
0529 h12mlnetacoefpl16_HB_noise->Clear();
0530 h12mlnetacoefpl16_HE_noise->Clear();
0531 h12mlnetacoefpl17_noise->Clear();
0532 h12mlnetacoefpl18_noise->Clear();
0533 h12mlnetacoefpl19_noise->Clear();
0534 h12mlnetacoefpl20_noise->Clear();
0535 h12mlnetacoefpl21_noise->Clear();
0536 h12mlnetacoefpl22_noise->Clear();
0537 h12mlnetacoefpl23_noise->Clear();
0538 h12mlnetacoefpl24_noise->Clear();
0539 h12mlnetacoefpl25_noise->Clear();
0540 h12mlnetacoefpl26_noise->Clear();
0541 h12mlnetacoefpl27_noise->Clear();
0542 h12mlnetacoefpl28_noise->Clear();
0543 h12mlnetacoefpl29_noise->Clear();
0544
0545 h12mlnetacoefpl1->Clear();
0546 h12mlnetacoefpl2->Clear();
0547 h12mlnetacoefpl3->Clear();
0548 h12mlnetacoefpl4->Clear();
0549 h12mlnetacoefpl5->Clear();
0550 h12mlnetacoefpl6->Clear();
0551 h12mlnetacoefpl7->Clear();
0552 h12mlnetacoefpl8->Clear();
0553 h12mlnetacoefpl9->Clear();
0554 h12mlnetacoefpl10->Clear();
0555 h12mlnetacoefpl11->Clear();
0556 h12mlnetacoefpl12->Clear();
0557 h12mlnetacoefpl13->Clear();
0558 h12mlnetacoefpl14->Clear();
0559 h12mlnetacoefpl15->Clear();
0560 h12mlnetacoefpl16_HB->Clear();
0561 h12mlnetacoefpl16_HE->Clear();
0562 h12mlnetacoefpl17->Clear();
0563 h12mlnetacoefpl18->Clear();
0564 h12mlnetacoefpl19->Clear();
0565 h12mlnetacoefpl20->Clear();
0566 h12mlnetacoefpl21->Clear();
0567 h12mlnetacoefpl22->Clear();
0568 h12mlnetacoefpl23->Clear();
0569 h12mlnetacoefpl24->Clear();
0570 h12mlnetacoefpl25->Clear();
0571 h12mlnetacoefpl26->Clear();
0572 h12mlnetacoefpl27->Clear();
0573 h12mlnetacoefpl28->Clear();
0574 h12mlnetacoefpl29->Clear();
0575
0576 h12mlnetacoefmin1_noise->Write();
0577 h12mlnetacoefmin2_noise->Write();
0578 h12mlnetacoefmin3_noise->Write();
0579 h12mlnetacoefmin4_noise->Write();
0580 h12mlnetacoefmin5_noise->Write();
0581 h12mlnetacoefmin6_noise->Write();
0582 h12mlnetacoefmin7_noise->Write();
0583 h12mlnetacoefmin8_noise->Write();
0584 h12mlnetacoefmin9_noise->Write();
0585 h12mlnetacoefmin10_noise->Write();
0586 h12mlnetacoefmin11_noise->Write();
0587 h12mlnetacoefmin12_noise->Write();
0588 h12mlnetacoefmin13_noise->Write();
0589 h12mlnetacoefmin14_noise->Write();
0590 h12mlnetacoefmin15_noise->Write();
0591 h12mlnetacoefmin16_HB_noise->Write();
0592 h12mlnetacoefmin16_HE_noise->Write();
0593 h12mlnetacoefmin17_noise->Write();
0594 h12mlnetacoefmin18_noise->Write();
0595 h12mlnetacoefmin19_noise->Write();
0596 h12mlnetacoefmin20_noise->Write();
0597 h12mlnetacoefmin21_noise->Write();
0598 h12mlnetacoefmin22_noise->Write();
0599 h12mlnetacoefmin23_noise->Write();
0600 h12mlnetacoefmin24_noise->Write();
0601 h12mlnetacoefmin25_noise->Write();
0602 h12mlnetacoefmin26_noise->Write();
0603 h12mlnetacoefmin27_noise->Write();
0604 h12mlnetacoefmin28_noise->Write();
0605 h12mlnetacoefmin29_noise->Write();
0606
0607 h12mlnetacoefmin1->Write();
0608 h12mlnetacoefmin2->Write();
0609 h12mlnetacoefmin3->Write();
0610 h12mlnetacoefmin4->Write();
0611 h12mlnetacoefmin5->Write();
0612 h12mlnetacoefmin6->Write();
0613 h12mlnetacoefmin7->Write();
0614 h12mlnetacoefmin8->Write();
0615 h12mlnetacoefmin9->Write();
0616 h12mlnetacoefmin10->Write();
0617 h12mlnetacoefmin11->Write();
0618 h12mlnetacoefmin12->Write();
0619 h12mlnetacoefmin13->Write();
0620 h12mlnetacoefmin14->Write();
0621 h12mlnetacoefmin15->Write();
0622 h12mlnetacoefmin16_HB->Write();
0623 h12mlnetacoefmin16_HE->Write();
0624 h12mlnetacoefmin17->Write();
0625 h12mlnetacoefmin18->Write();
0626 h12mlnetacoefmin19->Write();
0627 h12mlnetacoefmin20->Write();
0628 h12mlnetacoefmin21->Write();
0629 h12mlnetacoefmin22->Write();
0630 h12mlnetacoefmin23->Write();
0631 h12mlnetacoefmin24->Write();
0632 h12mlnetacoefmin25->Write();
0633 h12mlnetacoefmin26->Write();
0634 h12mlnetacoefmin27->Write();
0635 h12mlnetacoefmin28->Write();
0636 h12mlnetacoefmin29->Write();
0637
0638 h12mlnetacoefmin1_noise->Clear();
0639 h12mlnetacoefmin2_noise->Clear();
0640 h12mlnetacoefmin3_noise->Clear();
0641 h12mlnetacoefmin4_noise->Clear();
0642 h12mlnetacoefmin5_noise->Clear();
0643 h12mlnetacoefmin6_noise->Clear();
0644 h12mlnetacoefmin7_noise->Clear();
0645 h12mlnetacoefmin8_noise->Clear();
0646 h12mlnetacoefmin9_noise->Clear();
0647 h12mlnetacoefmin10_noise->Clear();
0648 h12mlnetacoefmin11_noise->Clear();
0649 h12mlnetacoefmin12_noise->Clear();
0650 h12mlnetacoefmin13_noise->Clear();
0651 h12mlnetacoefmin14_noise->Clear();
0652 h12mlnetacoefmin15_noise->Clear();
0653 h12mlnetacoefmin16_HB_noise->Clear();
0654 h12mlnetacoefmin16_HE_noise->Clear();
0655 h12mlnetacoefmin17_noise->Clear();
0656 h12mlnetacoefmin18_noise->Clear();
0657 h12mlnetacoefmin19_noise->Clear();
0658 h12mlnetacoefmin20_noise->Clear();
0659 h12mlnetacoefmin21_noise->Clear();
0660 h12mlnetacoefmin22_noise->Clear();
0661 h12mlnetacoefmin23_noise->Clear();
0662 h12mlnetacoefmin24_noise->Clear();
0663 h12mlnetacoefmin25_noise->Clear();
0664 h12mlnetacoefmin26_noise->Clear();
0665 h12mlnetacoefmin27_noise->Clear();
0666 h12mlnetacoefmin28_noise->Clear();
0667 h12mlnetacoefmin29_noise->Clear();
0668
0669 h12mlnetacoefmin1->Clear();
0670 h12mlnetacoefmin2->Clear();
0671 h12mlnetacoefmin3->Clear();
0672 h12mlnetacoefmin4->Clear();
0673 h12mlnetacoefmin5->Clear();
0674 h12mlnetacoefmin6->Clear();
0675 h12mlnetacoefmin7->Clear();
0676 h12mlnetacoefmin8->Clear();
0677 h12mlnetacoefmin9->Clear();
0678 h12mlnetacoefmin10->Clear();
0679 h12mlnetacoefmin11->Clear();
0680 h12mlnetacoefmin12->Clear();
0681 h12mlnetacoefmin13->Clear();
0682 h12mlnetacoefmin14->Clear();
0683 h12mlnetacoefmin15->Clear();
0684 h12mlnetacoefmin16_HB->Clear();
0685 h12mlnetacoefmin16_HB->Clear();
0686 h12mlnetacoefmin17->Clear();
0687 h12mlnetacoefmin18->Clear();
0688 h12mlnetacoefmin19->Clear();
0689 h12mlnetacoefmin20->Clear();
0690 h12mlnetacoefmin21->Clear();
0691 h12mlnetacoefmin22->Clear();
0692 h12mlnetacoefmin23->Clear();
0693 h12mlnetacoefmin24->Clear();
0694 h12mlnetacoefmin25->Clear();
0695 h12mlnetacoefmin26->Clear();
0696 h12mlnetacoefmin27->Clear();
0697 h12mlnetacoefmin28->Clear();
0698 h12mlnetacoefmin29->Clear();
0699
0700 h12mlnetacoefpl2D1->Write();
0701 h12mlnetacoefpl2D1_noise->Write();
0702 h12mlnetacoefmin2D1->Write();
0703 h12mlnetacoefmin2D1_noise->Write();
0704 h12mlnetacoefpl2D1->Clear();
0705 h12mlnetacoefpl2D1_noise->Clear();
0706 h12mlnetacoefmin2D1->Clear();
0707 h12mlnetacoefmin2D1_noise->Clear();
0708
0709 h12mlnetacoefpl2D16->Write();
0710 h12mlnetacoefpl2D16_noise->Write();
0711 h12mlnetacoefmin2D16->Write();
0712 h12mlnetacoefmin2D16_noise->Write();
0713 h12mlnetacoefpl2D16->Clear();
0714 h12mlnetacoefpl2D16_noise->Clear();
0715 h12mlnetacoefmin2D16->Clear();
0716 h12mlnetacoefmin2D16_noise->Clear();
0717
0718 }