merge_Energy_Statistic_And_SytematicSources_Gamma.C File Reference

#include <iomanip.h>

Include dependency graph for merge_Energy_Statistic_And_SytematicSources_Gamma.C:

Functions
void	merge_Energy_Statistic_And_SytematicSources_Gamma ()

Function Documentation

merge_Energy_Statistic_And_SytematicSources_Gamma()

void merge_Energy_Statistic_And_SytematicSources_Gamma

(

)

ATTENTION: THE READION INPUT FILES MAY NOOOOOT CONTAIN NUMBERS WITHOUD DIGITS FOR THE VALUES: for example 2 -> has to be 2.0 !!!

---

---

---

---

The this output goes into: libShower_Energy_STATandSYS_Gamma.txt

{
    
    
    
//  Bool_t LINEARERRORADDION=kTRUE;
    Bool_t LINEARERRORADDION=kFALSE;
 
    Int_t plate;
    Int_t dummy;
    Int_t plateArray[13]={10,12,14,16,18,20,23,26,29,32,35,40,45};
    Float_t dummyf;
    Int_t plate_i;
    Float_t sigmasysoutput;
    Float_t sigmasysinput[12][13];
  Float_t sigmasysoutput_tot[12][13];
    
    for (int i=0;i<13;i++) {
        for (int j=0;j<12;j++) { 
            sigmasysoutput_tot[j][i]=0;
        }
    }
    
// libShower_Energy_Statistics_Electron.txt
// libShower_Energy_Systematics_Electron.txt
 
// libShower_Energy_Statistics_Gamma.txt
// libShower_Energy_Systematics_Gamma.txt
 
 
  const char* name = "libShower_Energy_Statistics_Gamma.txt";
    pFile = fopen (name,"r");
    for (int i=0;i<13;i++) {
    plate_i=i;
    cout << plate_i << "---------------------------------" << endl;
    fscanf (pFile,"%d",&dummy);
    for (int j=0;j<12;j++) { 
        fscanf (pFile," %f ", &dummyf); cout << dummyf<< "  ";  sigmasysinput[j][plate_i]=dummyf; 
        if (!LINEARERRORADDION) sigmasysoutput_tot[j][plate_i]+=dummyf*dummyf;
        if (LINEARERRORADDION) sigmasysoutput_tot[j][plate_i]+=dummyf;
    }
    cout << endl;
    }
    fclose (pFile);
    
 
  const char* name = "libShower_Energy_Systematics_Gamma.txt";
    pFile = fopen (name,"r");
    for (int i=0;i<13;i++) {
    plate_i=i;
    cout << plate_i << "---------------------------------" << endl;
    fscanf (pFile,"%d",&dummy);
    for (int j=0;j<12;j++) { 
        fscanf (pFile," %f ", &dummyf); cout << dummyf<< "  "; sigmasysinput[j][plate_i]=dummyf; 
        if (!LINEARERRORADDION) sigmasysoutput_tot[j][plate_i]+=dummyf*dummyf;
        if (LINEARERRORADDION) sigmasysoutput_tot[j][plate_i]+=dummyf;
    }
    cout << endl;
    }
    fclose (pFile);
 
    
 
    cout << endl << endl;
    
    for (int i=0;i<13;i++) {
        cout << plateArray[i] << "  ";
        for (int j=0;j<12;j++) { 
            if (!LINEARERRORADDION) dummyf=TMath::Sqrt(sigmasysoutput_tot[j][i]);
            if (LINEARERRORADDION) dummyf=sigmasysoutput_tot[j][i];
            sigmasysoutput_tot[j][i]=dummyf;
            printf("%2.2f  ",sigmasysoutput_tot[j][i]);
        }
        cout << endl;
    }
    
}