Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 12:48:14

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 // Choose depth
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 //    if( k == 2 && j == 1 && i == 16 ) cout<<plmeang_nn[i][j][idep][k]<<endl;
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 //    cout<<" mean "<<plmean1[i][idep][k]<<" "<<minmean1[i][idep][k]<<endl;
0096 //    cout<<" mean_nn "<<plmean1_nn[i][idep][k]<<" "<<minmean1_nn[i][idep][k]<<endl;
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 // ieta=16 is in HB and HE
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 //   ieta=16 is in HB and HE 
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 // Two dimensional histogram 
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 //      if(i==1) cout<<" First "<<plmeang_nn[i][j][idep][k]<<" "<<plnoise[i][j][idep][k]<<" "<<plmeang_nn[i][j][idep][k] - plnoise[i][j][idep][k]<<" "<<plmean1[1][idep][k]<<
0251 //          " "<<plmean1_nn[1][idep][k]<<endl;
0252     
0253       
0254 //      if( i == 16) cout<<j<<" "<<k<<" "<<plmeang_nn[i][j][idep][k]-plnoise[i][j][idep][k]<<" "<<plmeang_nn[i][j][idep][k]<<" "<<plnoise[i][j][idep][k]<<endl;     
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     //       if( i == 16 && k == 1 ) cout<<" "<<plmeang_nn[i][j][idep][k]-plnoise[i][j][idep][k]<<" "<<plmean1_nn[i][idep][k]<<" "<<tt0<<endl;
0259        
0260        if( i == 1  ) {
0261 //        cout<<" 0J= "<<j<<endl;
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 //        cout<<" 1J= "<<j<<endl;
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 //        cout<<" 1nJ= "<<j<<endl;
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 //        cout<<" 2J= "<<j<<endl;
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 //        cout<<" 3J= "<<j<<endl;
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 // Print histograms    
0440     }
0441     }
0442 //    h12mlnetacoefpl1->Draw();
0443     fclose(Out1);
0444     fclose(Out2);
0445     
0446 //    h12mlnetacoefdist->Draw();
0447 //    h12mlnetacoefdist->Fit("gaus","","",0.8,1.1);
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 }