EMCal/EMCALPositionScanCodes/MakeFile.C

0001 void MakeFile(){
0002
0003   gSystem->Load("libg4eval.so");
0004   gSystem->Load("libqa_modules.so");
0005   gSystem->Load("libPrototype3.so");
0006  gStyle->SetOptFit(0);
0007    gStyle->SetOptStat(0);
0008
0009
0010 //TFile *fin = new TFile("2nd_positionscan.lst_EMCalCalib.root");
0011
0012
0013    TFile *fin = new TFile("3rd_positionscan.lst_EMCalCalib.root");
0014
0015 TTree *t = (TTree *)fin->Get("T");
0016
0017   t->SetAlias("C2_Inner_e", "1*abs(TOWER_RAW_C2[2].energy)");
0018   t->SetAlias("C2_Outer_e", "1*abs(TOWER_RAW_C2[3].energy)");
0019 t->SetAlias("Average_column", "Sum$(TOWER_CALIB_CEMC.get_column() * TOWER_CALIB_CEMC.get_energy())/Sum$(TOWER_CALIB_CEMC.get_energy())");
0020   t->SetAlias("Average_HODO_HORIZONTAL", "Sum$(TOWER_CALIB_HODO_HORIZONTAL.towerid * (abs(TOWER_CALIB_HODO_HORIZONTAL.energy)>30) * abs(TOWER_CALIB_HODO_HORIZONTAL.energy))/Sum$((abs(TOWER_CALIB_HODO_HORIZONTAL.energy)>30) * abs(TOWER_CALIB_HODO_HORIZONTAL.energy))");
0021   t->SetAlias("Valid_HODO_HORIZONTAL", "Sum$(abs(TOWER_CALIB_HODO_HORIZONTAL.energy)>30) > 0");
0022   t->SetAlias("No_Triger_VETO", "Sum$(abs(TOWER_RAW_TRIGGER_VETO.energy)>15)==0");
0023  t->SetAlias("Valid_HODO_VERTICAL", "Sum$(abs(TOWER_CALIB_HODO_VERTICAL.energy)>30) > 0");
0024   t->SetAlias("C2_Sum_e", "C2_Inner_e + C2_Outer_e");
0025   t->SetAlias("Average_HODO_VERTICAL","Sum$(TOWER_CALIB_HODO_VERTICAL.towerid * (abs(TOWER_CALIB_HODO_VERTICAL.energy)>30) * abs(TOWER_CALIB_HODO_VERTICAL.energy))/Sum$((abs(TOWER_CALIB_HODO_VERTICAL.energy)>30) * abs(TOWER_CALIB_HODO_VERTICAL.energy))");
0026 //return;
0027   t->SetAlias("Energy_Sum_CEMC", "1*Sum$(TOWER_CALIB_CEMC.get_energy())");
0028
0029 TH3F *Energyhis = new TH3F("Energyhis","",170,170,340,185,60,245,200,0,10);
0030
0031  t->Draw("clus_5x5_prod.sum_E:info.beam_2CV_mm + 5*info.hodo_v:info.beam_2CH_mm - 5*info.hodo_h>>Energyhis","good_e");
0032
0033
0034
0035 //TFile *fout =  new TFile("His2.root","RECREATE");
0036 TFile *fout =  new TFile("His3.root","RECREATE");
0037
0038
0039  Energyhis->Write();
0040
0041     double mean = Energyhis->GetMean();
0042
0043     cout << "Mean is =  "<< mean << endl;
0044
0045
0046 }