//Create by Christine Nattrass, Rebecca Scott, Irakli Martashvili //University of Tennessee at Knoxville //This macro takes an input file created by AliAnalysisHadEtMonteCarlo and creates an AliAnalysisHadEtCorrections for the determination of the corrected Et // #include "TFile.h" // #include <iostream> // #include "AliAnalysisHadEtCorrections.h" // #include <iostream> // #include <TROOT.h> // #include <TSystem.h> // #include "TStopwatch.h" //Corrections added in by hand to deal with the inadequacies of PYTHIA and HIJING TH1D *pp276TPCBkgd(); TH1D *pp276ITSBkgd(); TH1D *PbPb276TPCBkgd(); TH1D *PbPb276ITSBkgd(); Double_t xAxis1[112] = {0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1, 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.2, 0.22, 0.24, 0.26, 0.28, 0.3, 0.32, 0.34, 0.36, 0.38, 0.4, 0.42, 0.44, 0.46, 0.48, 0.5, 0.52, 0.54, 0.56, 0.58, 0.6, 0.62, 0.64, 0.66, 0.68, 0.7, 0.72, 0.74, 0.76, 0.78, 0.8, 0.82, 0.84, 0.86, 0.88, 0.9, 0.92, 0.94, 0.96, 0.98, 1, 1.05, 1.1, 1.15, 1.2, 1.25, 1.3, 1.35, 1.4, 1.45, 1.5, 1.55, 1.6, 1.65, 1.7, 1.75, 1.8, 1.85, 1.9, 1.95, 2, 2.2, 2.4, 2.6, 2.8, 3, 3.2, 3.4, 3.6, 3.8, 4, 4.2, 4.4, 4.6, 4.8, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10, 12, 14, 16, 18, 20, 25}; Float_t CorrNeutral(float ptcut, char *prodname, char *shortprodname, bool ispp = true, bool forSim = true, bool TPC, bool hadronic = false, float etacut = 0.7, int dataset); TH1D *GetHistoCorrNeutral(float cut, char *name, bool ispp, bool forSim, int mycase, bool eta, int color, int marker, bool hadronic = false); Float_t CorrPtCut(float ptcut, char *prodname = "Enter Production Name", char *shortprodname = "EnterProductionName", bool ispp = true, bool forSim = true, int mycase = 0); TH1D *GetHistoCorrPtCut(float ptcut = 0.15, char *name, bool ispp = true, bool forSim = true, int mycase); TH1D *GetHistoCorrNotID(float etacut,char *name, bool TPC, bool eta, bool ispp = true, bool forSim = true); TH1D *CorrNotID(float etacut,char *name, char *prodname, char *shortprodname, bool TPC, bool ispp = true, bool forSim = true); Float_t CorrNotIDConst(float ptcut, float etacut,char *name, char *prodname, char *shortprodname, bool TPC, bool ispp, bool forSim); TH1D *GetHistoNoID(float etacut, char *name, bool eta, bool TPC, bool ispp, bool forSim); TH1D *CorrNoID(float etacut,char *name, char *prodname, char *shortprodname, bool ispp, bool forSim); Float_t CorrNoIDConst(float etacut, float ptcut,char *name, char *prodname, char *shortprodname, bool ispp, bool forSim); TH1D* bayneseffdiv(TH1D* numerator, TH1D* denominator,Char_t* name); TH1D *GetHistoEfficiency(float cut, char *name, int mycase, int color, int marker,bool TPC, bool ITS, int cb = -1, int cblast = -1); void CorrEfficiencyPlots(bool TPC, char *prodname, char *shortprodname); TH1D *GetHistoCorrBkgd(float etacut,char *name, bool TPC,bool ispp,bool forSim); void CorrBkgdPlots(char *prodname, char *shortprodname, bool TPC,bool ispp,bool forSim); //Some variables that we'll use multiple times. We'll declare them here since they don't seem to delete right in the functions char prefix[100]; char histoname[100]; char epsname[100]; char pngname[100]; TFile *file = NULL;//initiated in main function const char *mynameTPC = "TPC"; const char *mynameITS = "ITS"; const char *mynameTPCITS = "TPCITS"; const char *detectorEMCAL = "EMCAL"; const char *detectorPHOS = "PHOS"; const char *reweightedNo = ""; const char *reweightedYes = "Reweighted"; //=========================================================================================== void GetCorrections(char *prodname = "Enter Production Name", char *shortprodname = "EnterProductionName", bool ispp = true, bool forSim = true, bool TPC = true, char *infilename="Et.ESD.new.sim.merged.root", int dataset = 2009){ TStopwatch timer; timer.Start(); gSystem->Load("libTree.so"); gSystem->Load("libGeom.so"); gSystem->Load("libVMC.so"); gSystem->Load("libXMLIO.so"); gSystem->Load("libSTEERBase.so"); gSystem->Load("libESD.so"); gSystem->Load("libAOD.so"); gSystem->Load("libANALYSIS"); gSystem->Load("libANALYSISalice"); gSystem->AddIncludePath("-I$ALICE_ROOT/include"); gROOT->ProcessLine(".L AliAnalysisEtCuts.cxx+g"); gROOT->ProcessLine(".L AliAnalysisHadEtCorrections.cxx+g"); file = new TFile(infilename); char outfilename[200]; char *sim = "ForData"; if(forSim) sim = "ForSimulations"; char *system = "PbPb"; if(ispp){ system = "pp"; if(dataset==2013) system = "pPb"; } sprintf(outfilename,"rootFiles/corrections/corrections.%s.%s.%s.root",shortprodname,system,sim); TFile *outfile = new TFile(outfilename,"RECREATE"); AliAnalysisHadEtCorrections *hadCorrectionEMCAL = new AliAnalysisHadEtCorrections(); hadCorrectionEMCAL->SetName("hadCorrectionEMCAL"); float etacut = 0.7; hadCorrectionEMCAL->SetEtaCut(etacut); hadCorrectionEMCAL->IsData(!forSim); hadCorrectionEMCAL->IsEMCal(kTRUE); hadCorrectionEMCAL->SetProduction(shortprodname); hadCorrectionEMCAL->SetProductionDescription(prodname); hadCorrectionEMCAL->SetDataSet(dataset); //float etacut = hadCorrectionEMCAL->GetEtaCut(); //cout<<"eta cut is "<<etacut<<endl; hadCorrectionEMCAL->SetAcceptanceCorrectionFull(1.0); cout<<"Warning: Acceptance corrections will have to be updated to include real acceptance maps of the EMCAL"<<endl; hadCorrectionEMCAL->SetAcceptanceCorrectionPHOS(360.0/60.0); hadCorrectionEMCAL->SetAcceptanceCorrectionEMCAL(360.0/60.0); float ptcut = 0.1; float neutralCorr = CorrNeutral(ptcut,prodname,shortprodname,ispp,forSim,TPC,false,etacut,dataset); hadCorrectionEMCAL->SetNeutralCorrection(neutralCorr); //Using error from data, see analysis note for details if(ispp){ hadCorrectionEMCAL->SetNeutralCorrectionLowBound(neutralCorr*(1.0-0.014)); hadCorrectionEMCAL->SetNeutralCorrectionHighBound(neutralCorr*(1.0+0.014)); } else{ hadCorrectionEMCAL->SetNeutralCorrectionLowBound(neutralCorr*(1.0-0.049)); hadCorrectionEMCAL->SetNeutralCorrectionHighBound(neutralCorr*(1.0+0.049)); } float hadronicCorr = CorrNeutral(ptcut,prodname,shortprodname,ispp,forSim,TPC,true,etacut,dataset); hadCorrectionEMCAL->SetNotHadronicCorrection(hadronicCorr); if(ispp){ hadCorrectionEMCAL->SetNotHadronicCorrectionLowBound(hadronicCorr*(1.0-0.008)); hadCorrectionEMCAL->SetNotHadronicCorrectionHighBound(hadronicCorr*(1.0+0.008)); } else{ hadCorrectionEMCAL->SetNotHadronicCorrectionLowBound(hadronicCorr*(1.0-0.023)); hadCorrectionEMCAL->SetNotHadronicCorrectionHighBound(hadronicCorr*(1.0+0.023)); } float ptcutITS = CorrPtCut(0.1,prodname,shortprodname,ispp,forSim); hadCorrectionEMCAL->SetpTCutCorrectionITS(ptcutITS); float ptcutTPC = CorrPtCut(0.15,prodname,shortprodname,ispp,forSim); hadCorrectionEMCAL->SetpTCutCorrectionTPC(ptcutTPC); float ptcutITSLow = CorrPtCut(0.1,prodname,shortprodname,ispp,forSim,-1); float ptcutTPCLow = CorrPtCut(0.15,prodname,shortprodname,ispp,forSim,-1); hadCorrectionEMCAL->SetpTCutCorrectionITSLowBound(ptcutITSLow); hadCorrectionEMCAL->SetpTCutCorrectionTPCLowBound(ptcutTPCLow); float ptcutITSHigh = CorrPtCut(0.1,prodname,shortprodname,ispp,forSim,1); float ptcutTPCHigh = CorrPtCut(0.15,prodname,shortprodname,ispp,forSim,1); hadCorrectionEMCAL->SetpTCutCorrectionITSHighBound(ptcutITSHigh); hadCorrectionEMCAL->SetpTCutCorrectionTPCHighBound(ptcutTPCHigh); TH1D *NotIDTPC = CorrNotID(etacut,"CorrNotIDEMCALTPC",prodname,shortprodname,true,ispp,forSim); TH1D *NotIDITS = CorrNotID(etacut,"CorrNotIDEMCALITS",prodname,shortprodname,false,ispp,forSim); hadCorrectionEMCAL->SetNotIDCorrectionTPC(NotIDTPC); hadCorrectionEMCAL->SetNotIDCorrectionITS(NotIDITS); Float_t NotIDConstTPC = CorrNotIDConst(0.15,etacut,"CorrNotIDEMCALTPC2",prodname,shortprodname,true,ispp,forSim); Float_t NotIDConstITS = CorrNotIDConst(0.10,etacut,"CorrNotIDEMCALTPC2",prodname,shortprodname,true,ispp,forSim); hadCorrectionEMCAL->SetNotIDConstCorrectionTPC(1.0/NotIDConstTPC); hadCorrectionEMCAL->SetNotIDConstCorrectionITS(1.0/NotIDConstITS); if(ispp){ hadCorrectionEMCAL->SetNotIDConstCorrectionTPCLowBound(1./NotIDConstTPC*(1.0-0.010)); hadCorrectionEMCAL->SetNotIDConstCorrectionITSLowBound(1./NotIDConstITS*(1.0-0.010)); hadCorrectionEMCAL->SetNotIDConstCorrectionTPCHighBound(1./NotIDConstTPC*(1.0+0.010)); hadCorrectionEMCAL->SetNotIDConstCorrectionITSHighBound(1./NotIDConstITS*(1.0+0.010)); } else{ hadCorrectionEMCAL->SetNotIDConstCorrectionTPCLowBound(1./NotIDConstTPC*(1.0-0.022)); hadCorrectionEMCAL->SetNotIDConstCorrectionITSLowBound(1./NotIDConstITS*(1.0-0.022)); hadCorrectionEMCAL->SetNotIDConstCorrectionTPCHighBound(1./NotIDConstTPC*(1.0+0.022)); hadCorrectionEMCAL->SetNotIDConstCorrectionITSHighBound(1./NotIDConstITS*(1.0+0.022)); } TH1D *NoID = CorrNoID(etacut,"CorrNoIDEMCAL",prodname,shortprodname,ispp,forSim); hadCorrectionEMCAL->SetNotIDCorrectionNoPID(NoID); Float_t NoIDTPC = CorrNoIDConst(etacut,0.15,"CorrNoIDEMCAL2",prodname,shortprodname,ispp,forSim); Float_t NoIDITS = CorrNoIDConst(etacut,0.1,"CorrNoIDEMCAL2",prodname,shortprodname,ispp,forSim); hadCorrectionEMCAL->SetNotIDConstCorrectionTPCNoID(1./NoIDTPC); hadCorrectionEMCAL->SetNotIDConstCorrectionITSNoID(1./NoIDITS); hadCorrectionEMCAL->SetNotIDConstCorrectionTPCNoIDLowBound(1./NoIDTPC*.98); hadCorrectionEMCAL->SetNotIDConstCorrectionITSNoIDLowBound(1./NoIDITS*.98); hadCorrectionEMCAL->SetNotIDConstCorrectionTPCNoIDHighBound(1./NoIDTPC*1.02); hadCorrectionEMCAL->SetNotIDConstCorrectionITSNoIDHighBound(1./NoIDITS*1.02); //Here we're going to do a bit of a cheat. We want the efficiency for ITS standalone tracks + TPC+ITS tracks. This is returned by the function for the efficiency function if I ask for ITS only efficiency. Had I known how this worked I probably would have written the code differently... but... //anyhow I left a switch for changing it back. bool useITSStandalone = false; TH1D *efficiencyPionTPC = GetHistoEfficiency(etacut,"hEfficiencyPionTPC",1,1,20,useITSStandalone,true); hadCorrectionEMCAL->SetEfficiencyPionTPC(efficiencyPionTPC); TH1D *efficiencyKaonTPC = GetHistoEfficiency(etacut,"hEfficiencyKaonTPC",2,1,20,useITSStandalone,true); hadCorrectionEMCAL->SetEfficiencyKaonTPC(efficiencyKaonTPC); TH1D *efficiencyProtonTPC = GetHistoEfficiency(etacut,"hEfficiencyProtonTPC",3,1,20,useITSStandalone,true); hadCorrectionEMCAL->SetEfficiencyProtonTPC(efficiencyProtonTPC); TH1D *efficiencyHadronTPC = GetHistoEfficiency(etacut,"hEfficiencyHadronTPC",0,1,20,useITSStandalone,true); hadCorrectionEMCAL->SetEfficiencyHadronTPC(efficiencyHadronTPC); TH1D *efficiencyPionITS = GetHistoEfficiency(etacut,"hEfficiencyPionITS",1,1,20,false,true); hadCorrectionEMCAL->SetEfficiencyPionITS(efficiencyPionITS); TH1D *efficiencyKaonITS = GetHistoEfficiency(etacut,"hEfficiencyKaonITS",2,1,20,false,true); hadCorrectionEMCAL->SetEfficiencyKaonITS(efficiencyKaonITS); TH1D *efficiencyProtonITS = GetHistoEfficiency(etacut,"hEfficiencyProtonITS",3,1,20,false,true); hadCorrectionEMCAL->SetEfficiencyProtonITS(efficiencyProtonITS); TH1D *efficiencyHadronITS = GetHistoEfficiency(etacut,"hEfficiencyHadronITS",0,1,20,false,true); hadCorrectionEMCAL->SetEfficiencyHadronITS(efficiencyHadronITS); if(!ispp){ TH1D *efficiencyPionTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB0",1,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB5",1,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB10",1,1,20,useITSStandalone,true,10,15); for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB10->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB0",2,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB5",2,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB10",2,1,20,useITSStandalone,true,10,15); for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB10->Clone(Form("Test%i",i)),i);//Kaon TH1D *efficiencyProtonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB0",3,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB5",3,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB10",3,1,20,useITSStandalone,true,10,15); for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB10->Clone(Form("Test%i",i)),i);//Proton TH1D *efficiencyHadronTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB0",0,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB5",0,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB10",0,1,20,useITSStandalone,true,10,15); for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB10->Clone(Form("Test%i",i)),i);//Hadron TH1D *efficiencyPionITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB0",1,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB5",1,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB10",1,1,20,false,true,10,15); for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB10->Clone(Form("Test%i",i)),i);//Pion TH1D *efficiencyKaonITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB0",2,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB5",2,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB10",2,1,20,false,true,10,15); for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB10->Clone(Form("Test%i",i)),i);//Kaon TH1D *efficiencyProtonITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB0",3,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB5",3,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB10",3,1,20,false,true,10,15); for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB10->Clone(Form("Test%i",i)),i);//Proton TH1D *efficiencyHadronITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB0",0,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB5",0,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB10",0,1,20,false,true,10,15); for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB10->Clone(Form("Test%i",i)),i);//Hadron }//EMCAL hadCorrectionEMCAL->SetEfficiencyErrorLowBound(0.95); hadCorrectionEMCAL->SetEfficiencyErrorHighBound(1.05); //CorrEfficiencyPlots(true,prodname,shortprodname); //CorrEfficiencyPlots(false,prodname,shortprodname,infilename); //hadCorrectionEMCAL->GetEfficiencyHadronTPC()->Draw(); TH1D *backgroundTPC; TH1D *backgroundITS; if((dataset==20111 || dataset==20100|| dataset==2011) && !forSim){//2.76 TeV p+p or Pb+Pb if(dataset==20111){ cout<<"Fixing 2.76 TeV p+p background to be average of 900 GeV and 7 TeV scaling"<<endl; backgroundTPC = pp276TPCBkgd(); backgroundTPC->SetName("hBackgroundTPC"); backgroundITS = pp276ITSBkgd(); backgroundITS->SetName("hBackgroundITS"); } else{//PbPb cout<<"Fixing 2.76 TeV Pb+Pb background to be average of 900 GeV and 7 TeV scaling with baryon enhancement"<<endl; backgroundTPC = pp276TPCBkgd(); backgroundTPC->SetName("hBackgroundTPC"); //ITS background is currently a placeholder backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); } } else{ backgroundTPC = GetHistoCorrBkgd(etacut,"hBackgroundTPC",true,ispp,forSim); backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); } float bkgdpcterror = 0.0; switch(dataset){ case 2009: bkgdpcterror = 0.37; break; case 20111: bkgdpcterror = 0.38; break; case 2010: case 2012: case 2013: bkgdpcterror = 0.13; break; case 20100: case 2011: bkgdpcterror = 0.76; break; } hadCorrectionEMCAL->SetBackgroundCorrectionTPC(backgroundTPC); hadCorrectionEMCAL->SetBackgroundCorrectionITS(backgroundITS); hadCorrectionEMCAL->SetBackgroundErrorLowBound(1.0-bkgdpcterror/100.0); hadCorrectionEMCAL->SetBackgroundErrorHighBound(1.0+bkgdpcterror/100.0); //CorrBkgdPlots(prodname,shortprodname,true,ispp,forSim); //CorrBkgdPlots(prodname,shortprodname,false,ispp,forSim); hadCorrectionEMCAL->Report(); outfile->cd(); hadCorrectionEMCAL->Write(); outfile->Write(); delete hadCorrectionEMCAL; AliAnalysisHadEtCorrections *hadCorrectionPHOS = new AliAnalysisHadEtCorrections(); hadCorrectionPHOS->SetName("hadCorrectionPHOS"); float etacut = 0.12; hadCorrectionPHOS->SetEtaCut(etacut); hadCorrectionPHOS->IsData(!forSim); hadCorrectionPHOS->IsEMCal(kTRUE); hadCorrectionPHOS->SetProduction(shortprodname); hadCorrectionPHOS->SetProductionDescription(prodname); hadCorrectionPHOS->SetDataSet(dataset); //float etacut = hadCorrectionPHOS->GetEtaCut(); //cout<<"eta cut is "<<etacut<<endl; hadCorrectionPHOS->SetAcceptanceCorrectionFull(1.0); cout<<"Warning: Acceptance corrections will have to be updated to include real acceptance maps of the PHOS"<<endl; hadCorrectionPHOS->SetAcceptanceCorrectionPHOS(360.0/60.0); hadCorrectionPHOS->SetAcceptanceCorrectionEMCAL(360.0/60.0); float ptcut = 0.1; float neutralCorr = CorrNeutral(ptcut,prodname,shortprodname,ispp,forSim,TPC,false,etacut,dataset); hadCorrectionPHOS->SetNeutralCorrection(neutralCorr); //Using error from data, see analysis note for details if(ispp){ hadCorrectionPHOS->SetNeutralCorrectionLowBound(neutralCorr*(1.0-.014)); hadCorrectionPHOS->SetNeutralCorrectionHighBound(neutralCorr*(1.0+.014)); } else{ hadCorrectionPHOS->SetNeutralCorrectionLowBound(neutralCorr*(1.0-0.049)); hadCorrectionPHOS->SetNeutralCorrectionHighBound(neutralCorr*(1.0+0.049)); } float hadronicCorr = CorrNeutral(ptcut,prodname,shortprodname,ispp,forSim,TPC,true,etacut,dataset); hadCorrectionPHOS->SetNotHadronicCorrection(hadronicCorr); if(ispp){ hadCorrectionPHOS->SetNotHadronicCorrectionLowBound(neutralCorr*(1.0-0.008)); hadCorrectionPHOS->SetNotHadronicCorrectionHighBound(neutralCorr*(1.0+0.008)); } else{ hadCorrectionPHOS->SetNotHadronicCorrectionLowBound(neutralCorr*(1.0-0.023)); hadCorrectionPHOS->SetNotHadronicCorrectionHighBound(neutralCorr*(1.0+0.023)); } float ptcutITS = CorrPtCut(0.1,prodname,shortprodname,ispp,forSim); hadCorrectionPHOS->SetpTCutCorrectionITS(ptcutITS); float ptcutTPC = CorrPtCut(0.15,prodname,shortprodname,ispp,forSim); hadCorrectionPHOS->SetpTCutCorrectionTPC(ptcutTPC); float ptcutITSLow = CorrPtCut(0.1,prodname,shortprodname,ispp,forSim,-1); float ptcutTPCLow = CorrPtCut(0.15,prodname,shortprodname,ispp,forSim,-1); hadCorrectionPHOS->SetpTCutCorrectionITSLowBound(ptcutITSLow); hadCorrectionPHOS->SetpTCutCorrectionTPCLowBound(ptcutTPCLow); float ptcutITSHigh = CorrPtCut(0.1,prodname,shortprodname,ispp,forSim,1); float ptcutTPCHigh = CorrPtCut(0.15,prodname,shortprodname,ispp,forSim,1); hadCorrectionPHOS->SetpTCutCorrectionITSHighBound(ptcutITSHigh); hadCorrectionPHOS->SetpTCutCorrectionTPCHighBound(ptcutTPCHigh); TH1D *NotIDTPC = CorrNotID(etacut,"CorrNotIDPHOSTPC",prodname,shortprodname,true,ispp,forSim); TH1D *NotIDITS = CorrNotID(etacut,"CorrNotIDPHOSITS",prodname,shortprodname,false,ispp,forSim); hadCorrectionPHOS->SetNotIDCorrectionTPC(NotIDTPC); hadCorrectionPHOS->SetNotIDCorrectionITS(NotIDITS); Float_t NotIDConstTPC = CorrNotIDConst(0.15,etacut,"CorrNotIDPHOSTPC2",prodname,shortprodname,true,ispp,forSim); Float_t NotIDConstITS = CorrNotIDConst(0.10,etacut,"CorrNotIDPHOSTPC2",prodname,shortprodname,true,ispp,forSim); hadCorrectionPHOS->SetNotIDConstCorrectionTPC(1./NotIDConstTPC); hadCorrectionPHOS->SetNotIDConstCorrectionITS(1./NotIDConstITS); if(ispp){ hadCorrectionPHOS->SetNotIDConstCorrectionTPCLowBound(1./NotIDConstTPC*(1.0-0.010)); hadCorrectionPHOS->SetNotIDConstCorrectionITSLowBound(1./NotIDConstITS*(1.0-0.010)); hadCorrectionPHOS->SetNotIDConstCorrectionTPCHighBound(1./NotIDConstTPC*(1.0+0.010)); hadCorrectionPHOS->SetNotIDConstCorrectionITSHighBound(1./NotIDConstITS*(1.0+0.010)); } else{ hadCorrectionPHOS->SetNotIDConstCorrectionTPCLowBound(1./NotIDConstTPC*(1.0-0.022)); hadCorrectionPHOS->SetNotIDConstCorrectionITSLowBound(1./NotIDConstITS*(1.0-0.022)); hadCorrectionPHOS->SetNotIDConstCorrectionTPCHighBound(1./NotIDConstTPC*(1.0+0.022)); hadCorrectionPHOS->SetNotIDConstCorrectionITSHighBound(1./NotIDConstITS*(1.0+0.022)); } TH1D *NoID = CorrNoID(etacut,"CorrNoIDPHOS",prodname,shortprodname,ispp,forSim); hadCorrectionPHOS->SetNotIDCorrectionNoPID(NoID); Float_t NoIDTPC = CorrNoIDConst(etacut,0.15,"CorrNoIDPHOS2",prodname,shortprodname,ispp,forSim); Float_t NoIDITS = CorrNoIDConst(etacut,0.1,"CorrNoIDPHOS2",prodname,shortprodname,ispp,forSim); hadCorrectionPHOS->SetNotIDConstCorrectionTPCNoID(1./NoIDTPC); hadCorrectionPHOS->SetNotIDConstCorrectionITSNoID(1./NoIDITS); hadCorrectionPHOS->SetNotIDConstCorrectionTPCNoIDLowBound(1./NoIDTPC*.98); hadCorrectionPHOS->SetNotIDConstCorrectionITSNoIDLowBound(1./NoIDITS*.98); hadCorrectionPHOS->SetNotIDConstCorrectionTPCNoIDHighBound(1./NoIDTPC*1.02); hadCorrectionPHOS->SetNotIDConstCorrectionITSNoIDHighBound(1./NoIDITS*1.02); TH1D *efficiencyPionTPC = GetHistoEfficiency(etacut,"hEfficiencyPionTPC",1,1,20,useITSStandalone,true); TH1D *efficiencyKaonTPC = GetHistoEfficiency(etacut,"hEfficiencyKaonTPC",2,1,20,useITSStandalone,true); TH1D *efficiencyProtonTPC = GetHistoEfficiency(etacut,"hEfficiencyProtonTPC",3,1,20,useITSStandalone,true); TH1D *efficiencyHadronTPC = GetHistoEfficiency(etacut,"hEfficiencyHadronTPC",0,1,20,useITSStandalone,true); TH1D *efficiencyPionITS = GetHistoEfficiency(etacut,"hEfficiencyPionITS",1,1,20,false,true); TH1D *efficiencyKaonITS = GetHistoEfficiency(etacut,"hEfficiencyKaonITS",2,1,20,false,true); TH1D *efficiencyProtonITS = GetHistoEfficiency(etacut,"hEfficiencyProtonITS",3,1,20,false,true); TH1D *efficiencyHadronITS = GetHistoEfficiency(etacut,"hEfficiencyHadronITS",0,1,20,false,true); //CorrEfficiencyPlots(true,prodname,shortprodname); //CorrEfficiencyPlots(false,prodname,shortprodname); hadCorrectionPHOS->SetEfficiencyPionTPC(efficiencyPionTPC); hadCorrectionPHOS->SetEfficiencyKaonTPC(efficiencyKaonTPC); hadCorrectionPHOS->SetEfficiencyProtonTPC(efficiencyProtonTPC); hadCorrectionPHOS->SetEfficiencyHadronTPC(efficiencyHadronTPC); hadCorrectionPHOS->SetEfficiencyPionITS(efficiencyPionITS); hadCorrectionPHOS->SetEfficiencyKaonITS(efficiencyKaonITS); hadCorrectionPHOS->SetEfficiencyProtonITS(efficiencyProtonITS); hadCorrectionPHOS->SetEfficiencyHadronITS(efficiencyHadronITS); if(!ispp){ TH1D *efficiencyPionTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB0",1,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB5",1,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB10",1,1,20,useITSStandalone,true,10,15); for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB10->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB0",2,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB5",2,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB10",2,1,20,useITSStandalone,true,10,15); for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB10->Clone(Form("Test%i",i)),i);//Kaon TH1D *efficiencyProtonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB0",3,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB5",3,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB10",3,1,20,useITSStandalone,true,10,15); for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB10->Clone(Form("Test%i",i)),i);//Proton TH1D *efficiencyHadronTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB0",0,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB5",0,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB10",0,1,20,useITSStandalone,true,10,15); for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB10->Clone(Form("Test%i",i)),i);//Hadron TH1D *efficiencyPionITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB0",1,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB5",1,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB10",1,1,20,false,true,10,15); for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB10->Clone(Form("Test%i",i)),i);//Pion TH1D *efficiencyKaonITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB0",2,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB5",2,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB10",2,1,20,false,true,10,15); for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB10->Clone(Form("Test%i",i)),i);//Kaon TH1D *efficiencyProtonITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB0",3,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB5",3,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB10",3,1,20,false,true,10,15); for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB10->Clone(Form("Test%i",i)),i);//Proton TH1D *efficiencyHadronITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB0",0,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB5",0,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB10",0,1,20,false,true,10,15); for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB10->Clone(Form("Test%i",i)),i);//Hadron }//EMCAL hadCorrectionPHOS->SetEfficiencyErrorLowBound(0.95); hadCorrectionPHOS->SetEfficiencyErrorHighBound(1.05); TH1D *backgroundTPC; TH1D *backgroundITS; if((dataset==20111 || dataset==20100|| dataset==2011) && !forSim){//2.76 TeV p+p or Pb+Pb if(dataset==20111){ cout<<"Fixing 2.76 TeV p+p background to be average of 900 GeV and 7 TeV scaling"<<endl; backgroundTPC = pp276TPCBkgd(); backgroundTPC->SetName("hBackgroundTPCPHOS"); backgroundITS = pp276ITSBkgd(); backgroundITS->SetName("hBackgroundITSPHOS"); } else{//PbPb cout<<"Fixing 2.76 TeV Pb+Pb background to be average of 900 GeV and 7 TeV scaling with baryon enhancement"<<endl; backgroundTPC = pp276TPCBkgd(); backgroundTPC->SetName("hBackgroundTPCPHOS"); //ITS background is currently a placeholder backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); } } else{ backgroundTPC = GetHistoCorrBkgd(etacut,"hBackgroundTPC",true,ispp,forSim); backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); } hadCorrectionPHOS->SetBackgroundCorrectionTPC(backgroundTPC); hadCorrectionPHOS->SetBackgroundCorrectionITS(backgroundITS); hadCorrectionPHOS->SetBackgroundErrorLowBound(1.0-0.001); hadCorrectionPHOS->SetBackgroundErrorHighBound(1.0+0.001); //CorrBkgdPlots(prodname,shortprodname,true,ispp,forSim); //CorrBkgdPlots(prodname,shortprodname,false,ispp,forSim); hadCorrectionPHOS->Report(); //Write the output outfile->cd(); hadCorrectionPHOS->Write(); outfile->Write(); outfile->Close(); timer.Stop(); timer.Print(); } //==================================CorrNeutral============================================== Float_t CorrNeutral(float ptcut, char *prodname, char *shortprodname, bool ispp, bool forSim, bool TPC, bool hadronic, float etacut, int dataset){ if(!forSim){//for data we have evaluated the neutral correction from ALICE data if(hadronic){//for tot et from had et return 1/0.5665;//0.5665 0.0145 <-- Correct - what should be applied } else{//for had et only if(dataset==2009){ return 1.0/0.7571; } else{ return 1.0/0.751; } } } gStyle->SetOptTitle(0); gStyle->SetOptStat(0); gStyle->SetOptFit(0); TCanvas *c = new TCanvas("c","c",400,400); c->SetTopMargin(0.0); c->SetRightMargin(0.0); c->SetBorderSize(0); c->SetFillColor(0); c->SetFillColor(0); c->SetBorderMode(0); c->SetFrameFillColor(0); c->SetFrameBorderMode(0); TPad *ptpad = c->cd(1); ptpad->SetTopMargin(0.04); ptpad->SetRightMargin(0.04); ptpad->SetLeftMargin(0.149288); ptpad->SetBorderSize(0); ptpad->SetFillColor(0); ptpad->SetFillColor(0); ptpad->SetBorderMode(0); ptpad->SetFrameFillColor(0); ptpad->SetFrameBorderMode(0); int phosmarker = 20; sprintf(prefix,"%s%2.1f",shortprodname,ptcut); sprintf(histoname,"%stotal",histoname); int colortotal = 1; int casetotal = 4; if(hadronic) casetotal = 8; TH1D *total = GetHistoCorrNeutral(ptcut,histoname,ispp,forSim,casetotal,false,colortotal,phosmarker,hadronic); int colorallneutral = 2; TH1D *allneutral = GetHistoCorrNeutral(ptcut,"allneutral",ispp,forSim,3,false,colorallneutral,phosmarker,hadronic); int colorchargedsecondary = TColor::kViolet-3; TH1D *chargedsecondary = GetHistoCorrNeutral(ptcut,"chargedsecondary",ispp,forSim,2,false,colorchargedsecondary,phosmarker,hadronic); int colorneutralUndet = 4; TH1D *neutralUndet = GetHistoCorrNeutral(ptcut,"neutralUndet",ispp,forSim,1,false,colorneutralUndet,phosmarker,hadronic); int colorv0 = TColor::kGreen+2; TH1D *v0 = GetHistoCorrNeutral(ptcut,"v0",ispp,forSim,0,false,colorv0,phosmarker,hadronic); int colorem = TColor::kCyan; TH1D *em = GetHistoCorrNeutral(ptcut,"em",ispp,forSim,9,false,colorem,phosmarker,hadronic); TF1 *func = new TF1("func","[0]",-.7,.7); func->SetParameter(0,0.2); total->Fit(func,"","",-etacut,etacut); //total->SetAxisRange(0.0,4); total->GetXaxis()->SetLabelSize(0.05); total->GetYaxis()->SetLabelSize(0.045); total->GetXaxis()->SetTitleSize(0.05); total->GetYaxis()->SetTitleSize(0.06); if(hadronic){ total->SetMaximum(0.6); total->SetMinimum(0.0); } else{ total->SetMaximum(0.3); total->SetMinimum(0.0); } total->Draw(); allneutral->Draw("same"); chargedsecondary->Draw("same"); neutralUndet->Draw("same"); v0->Draw("same"); if(hadronic) em->Draw("same"); TLatex *tex = new TLatex(0.161478,1.0835,prodname); tex->SetTextSize(0.0537634); tex->Draw(); TLegend *leg2 = new TLegend(0.518321,0.746873,0.774812,0.955343); leg2->AddEntry(total,"Total"); leg2->AddEntry(allneutral,"#Lambda,#bar{#Lambda},K^{0}_{S},K^{0}_{L},n,#bar{n}"); leg2->AddEntry(neutralUndet,"K^{0}_{L},n,#bar{n}"); if(hadronic) leg2->AddEntry(em,"e^{#pm},#gamma,#eta,#pi^{0},#omega"); leg2->SetFillStyle(0); leg2->SetFillColor(0); leg2->SetBorderSize(0); leg2->SetTextSize(0.0548607); leg2->Draw(); if(hadronic){ sprintf(epsname,"pics/%s/fhadronic.eps",shortprodname); sprintf(pngname,"pics/%s/fhadronic.png",shortprodname); } else{ sprintf(epsname,"pics/%s/fneutral.eps",shortprodname); sprintf(pngname,"pics/%s/fneutral.png",shortprodname); } c->SaveAs(epsname); c->SaveAs(pngname); delete total; delete allneutral; delete chargedsecondary; delete neutralUndet; delete v0; delete em; delete c; float corr = func->GetParameter(0); delete func; delete tex; delete leg2; return 1.0/(1.0-corr); } TH1D *GetHistoCorrNeutral(float cut, char *name, bool ispp, bool forSim, int mycase, bool eta, int color, int marker, bool hadronic){ file->cd(); char *reweightname = ""; if(!forSim) reweightname = "Reweighted"; TH2F *numeratorParent; switch(mycase){ case 0: numeratorParent= (TH2F*)((TH2F*) out2->FindObject(Form("EtSimulatedLambda%s",reweightname)))->Clone("v0"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedAntiLambda%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedK0S%s",reweightname))); break; case 1: numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtSimulatedK0L%s",reweightname)))->Clone("Knnbar"); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedNeutron")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiNeutron")); break; case 2: numeratorParent = (TH2F*)((TH2F*) out2->FindObject("EtSimulatedOmega"))->Clone("ch2ndary"); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiOmega")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedXi")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedXi0")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi0")); break; case 3: numeratorParent= (TH2F*)((TH2F*) out2->FindObject(Form("EtSimulatedLambda%s",reweightname)))->Clone("allneutral"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedAntiLambda%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedK0S%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedK0L%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedNeutron")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiNeutron")); break; case 4: numeratorParent= (TH2F*)((TH2F*) out2->FindObject(Form("EtSimulatedLambda%s",reweightname)))->Clone("allneutral"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedAntiLambda%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedK0S%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedK0L%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedNeutron")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiNeutron")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedOmega")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiOmega")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedXi")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedXi0")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi0")); break; case 5: numeratorParent= (TH2F*)((TH2F*) out2->FindObject("EtSimulatedXi"))->Clone("allxi"); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedXi0")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi0")); break; case 6: numeratorParent = (TH2F*)((TH2F*) out2->FindObject("EtSimulatedOmega"))->Clone("allomega"); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiOmega")); break; case 7: numeratorParent = (TH2F*)((TH2F*) out2->FindObject("EtSimulatedSigma"))->Clone("allsigma"); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiSigma")); break; case 8: numeratorParent= (TH2F*)((TH2F*) out2->FindObject(Form("EtSimulatedLambda%s",reweightname)))->Clone("allneutral"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedAntiLambda%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedK0S%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtSimulatedK0L%s",reweightname))); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedNeutron")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiNeutron")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedOmega")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiOmega")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedXi")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedXi0")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi0")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedGamma")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedEta")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedPi0")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedOmega0")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedEPlus")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedEMinus")); break; case 9: numeratorParent= (TH2F*)((TH2F*) out2->FindObject("EtSimulatedGamma"))->Clone("allem"); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedEta")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedPi0")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedOmega0")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedEPlus")); numeratorParent->Add((TH2F*) out2->FindObject("EtSimulatedEMinus")); break; } TH2F *allhad; //allhad=(TH2F*) ((TH2F*) out2->FindObject("EtSimulatedAllHadron"))->Clone("id"); allhad=(TH2F*) ((TH2F*) out2->FindObject("EtSimulatedPiPlus"))->Clone("id"); allhad->Add((TH2F*) out2->FindObject("EtSimulatedPiMinus")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedKMinus")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedKPlus")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedProton")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedAntiProton")); allhad->Add((TH2F*) out2->FindObject(Form("EtSimulatedLambda%s",reweightname))); allhad->Add((TH2F*) out2->FindObject(Form("EtSimulatedAntiLambda%s",reweightname))); allhad->Add((TH2F*) out2->FindObject(Form("EtSimulatedK0S%s",reweightname))); allhad->Add((TH2F*) out2->FindObject(Form("EtSimulatedK0L%s",reweightname))); allhad->Add((TH2F*) out2->FindObject("EtSimulatedNeutron")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedAntiNeutron")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedEPlus")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedEMinus")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedOmega")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedAntiOmega")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedXi")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedSigma")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedAntiSigma")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedXi0")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedAntiXi0")); if(hadronic){//if we are getting the correction for the hadronic only case... allhad->Add((TH2F*) out2->FindObject("EtSimulatedGamma")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedEta")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedPi0")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedOmega0")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedEPlus")); allhad->Add((TH2F*) out2->FindObject("EtSimulatedEMinus")); } TH1D *denominator; TH1D *numerator; if(eta){ int lowbin = numeratorParent->GetYaxis()->FindBin(-cut+.001);//make sure we don't accv0entally get the wrong bin int highbin = numeratorParent->GetYaxis()->FindBin(cut-.001); //cout<<"Projecting from "<<numeratorParent->GetYaxis()->GetBinLowEdge(lowbin)<<" to "<<numeratorParent->GetYaxis()->GetBinLowEdge(highbin+1)<<endl; denominator = allhad->ProjectionX("name",lowbin,highbin); numerator = numeratorParent->ProjectionX("numerator",lowbin,highbin); } else{ int lowbin = allhad->GetXaxis()->FindBin(cut);//make sure we don't accidentally get the wrong bin int highbin = allhad->GetXaxis()->GetNbins(); //cout<<"Projecting from "<<allhad->GetXaxis()->GetBinLowEdge(lowbin)<<" to "<<allhad->GetXaxis()->GetBinLowEdge(highbin+1)<<endl; numerator = numeratorParent->ProjectionY("name",lowbin,highbin); denominator = allhad->ProjectionY("denominator",lowbin,highbin); } numerator->Divide(denominator); //numerator->Rebin(2); //numerator->Scale(0.5); numerator->SetYTitle("E_{T}^{had,sample}/E_{T}^{had,total}"); numerator->GetYaxis()->SetTitleOffset(1.2); numerator->SetMarkerColor(color); numerator->SetLineColor(color); numerator->SetMarkerStyle(marker); delete denominator; delete numeratorParent; delete allhad; //file->Close(); numerator->SetName(name); return numerator; } //===============================CorrPtCut========================================= TH1D *GetHistoCorrPtCut(float ptcut, char *name, bool ispp, bool forSim, int mycase){ file->cd(); TH2F *allhad = ((TH2F*) out2->FindObject("EtSimulatedAllHadron"))->Clone("allhad"); TH2F *ptlow = ((TH2F*) out2->FindObject("EtSimulatedChargedHadronAssumingNoPt"))->Clone("ptlow"); TH2F *pthigh; if(ptcut>0.14){//TPC cut off (TH2F*)pthigh =(TH2F*) ((TH2F*) out2->FindObject("EtSimulatedChargedHadronAssumingPtTPCCut"))->Clone("pthigh"); } else{ (TH2F*)pthigh =(TH2F*) ((TH2F*) out2->FindObject("EtSimulatedChargedHadronAssumingPtITSCut"))->Clone("pthigh"); } int lowbin = allhad->GetXaxis()->FindBin(0.0);//make sure we don't accidentally get the wrong bin int highbin = allhad->GetXaxis()->FindBin(ptcut); int nbins = allhad->GetXaxis()->GetNbins(); //cout<<"Projecting from "<<allhad->GetXaxis()->GetBinLowEdge(lowbin)<<" to "<<allhad->GetXaxis()->GetBinLowEdge(highbin+1)<<endl; //cout<<"Projecting from "<<allhad->GetXaxis()->GetBinLowEdge(lowbin)<<" to "<<allhad->GetXaxis()->GetBinLowEdge(nbins)<<endl; //allhad->Sumw2(); TH1D *numerator; TH1D *denominator; switch(mycase){ case -1: numerator = ptlow->ProjectionY("nameLow",lowbin,highbin); denominator = allhad->ProjectionY("denominatorLow",highbin,nbins); denominator->Add(numerator); break; case 1: numerator = pthigh->ProjectionY("nameHigh",lowbin,highbin); denominator = allhad->ProjectionY("denominatorHigh",highbin,nbins); denominator->Add(numerator); break; default: numerator = allhad->ProjectionY("name",lowbin,highbin); denominator = allhad->ProjectionY("denominator",lowbin,nbins); } numerator->Divide(denominator); numerator->SetYTitle("E_{T}^{had, p_{T}<cut-off}/E_{T}^{had, all p_{T}}"); numerator->GetYaxis()->SetTitleOffset(1.); numerator->GetYaxis()->SetTitleSize(0.08); numerator->GetYaxis()->SetLabelSize(0.05); numerator->GetXaxis()->SetTitleSize(0.08); numerator->GetXaxis()->SetLabelSize(0.05); numerator->GetXaxis()->SetTitleOffset(.6); //numerator->Rebin(2); //numerator->Scale(0.5); //numerator->Draw("e"); delete allhad; delete denominator; delete ptlow; delete pthigh; numerator->SetName(name); return numerator; } Float_t CorrPtCut(float ptcut, char *prodname, char *shortprodname, bool ispp, bool forSim, int mycase){ gStyle->SetOptTitle(0); gStyle->SetOptStat(0); gStyle->SetOptFit(0); TCanvas *c = new TCanvas("c","c",500,400); c->SetTopMargin(0.04); c->SetRightMargin(0.04); c->SetLeftMargin(0.181452); c->SetBottomMargin(0.134409); c->SetBorderSize(0); c->SetFillColor(0); c->SetFillColor(0); c->SetBorderMode(0); c->SetFrameFillColor(0); c->SetFrameBorderMode(0); TH1D *High = GetHistoCorrPtCut(0.15-.001,"High",ispp,forSim,mycase); TH1D *Low = GetHistoCorrPtCut(0.1-.001,"Low",ispp,forSim,mycase); TH1D *Lowest = GetHistoCorrPtCut(0.05-.001,"Lowest",ispp,forSim,mycase); TF1 *func = new TF1("func","[0]",-.7,.7); func->SetParameter(0,0.2); if(ptcut<.125){//its cuts Low->Fit(func); } else{//tpc cuts High->Fit(func); } High->SetMaximum(0.04); High->SetMinimum(0.0); High->SetMarkerColor(2); Low->SetMarkerColor(4); High->SetLineColor(2); Low->SetLineColor(4); High->SetMinimum(0.0); High->SetMarkerStyle(20); Low->SetMarkerStyle(21); Lowest->SetMarkerStyle(22); High->Draw(); Low->Draw("same"); Lowest->Draw("same"); TLatex *tex = new TLatex(-0.723444,0.0373593,prodname); tex->SetTextSize(0.0537634); tex->Draw(); TLegend *leg = new TLegend(0.217742,0.66129,0.477823,0.873656); leg->AddEntry(High,"p_{T} cut-off = 0.15 GeV/c"); leg->AddEntry(Low,"p_{T} cut-off = 0.1 GeV/c"); leg->AddEntry(Lowest,"p_{T} cut-off = 0.05 GeV/c"); leg->SetFillStyle(0); leg->SetFillColor(0); leg->SetBorderSize(0); leg->SetTextSize(0.0537634); leg->Draw(); if(ptcut<.125){//its cuts c->SaveAs(Form("pics/%s/fptcutITS.eps",shortprodname)); c->SaveAs(Form("pics/%s/fptcutITS.png",shortprodname)); } else{//tpc cuts c->SaveAs(Form("pics/%s/fptcutTPC.eps",shortprodname)); c->SaveAs(Form("pics/%s/fptcutTPC.png",shortprodname)); } float corr = func->GetParameter(0); //cout<<"Pt cut correction: "<<1.0/(1.0-corr)<<endl; delete High; delete Low; delete Lowest; delete func; delete c; delete tex; delete leg; return 1.0/(1.0-corr); } //==================================CorrNotID================================================= TH1D *GetHistoCorrNotID(float etacut,char *name, bool TPC, bool eta, bool ispp, bool forSim){ file->cd(); char *myname = mynameITS; if(TPC) myname = mynameTPCITS; TH2F *notid = ((TH2F*) out2->FindObject(Form("EtReconstructed%sUnidentifiedAssumingPion",myname)))->Clone("notid"); TH2F *nNotid = ((TH2F*) out2->FindObject(Form("EtNReconstructed%sUnidentified",myname)))->Clone("nNotid"); if(!eta){ //cout<<"Correction determined for all charged hadrons"<<endl; notid->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiPlus",myname))); notid->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiMinus",myname))); notid->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedKPlus",myname))); notid->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedKMinus",myname))); notid->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedProton",myname))); notid->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedAntiProton",myname))); nNotid->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%sPiPlus",myname))); nNotid->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%sPiMinus",myname))); nNotid->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%sKPlus",myname))); nNotid->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%sKMinus",myname))); nNotid->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%sProton",myname))); nNotid->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%sAntiProton",myname))); } TH2F *id = ((TH2F*) out2->FindObject(Form("EtReconstructed%sUnidentified",myname)))->Clone("id"); if(!eta){ id->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiPlus",myname))); id->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiMinus",myname))); id->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedKPlus",myname))); id->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedKMinus",myname))); id->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedProton",myname))); id->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedAntiProton",myname))); } TH1D *nNotidProj; TH1D *denominator; TH1D *numerator; if(eta){ int lowbin = notid->GetYaxis()->FindBin(-etacut+.001);//make sure we don't accv0entally get the wrong bin int highbin = notid->GetYaxis()->FindBin(etacut-.001); //cout<<"Projecting from "<<notid->GetYaxis()->GetBinLowEdge(lowbin)<<" to "<<notid->GetYaxis()->GetBinLowEdge(highbin+1)<<endl; denominator = id->ProjectionX("name",lowbin,highbin); numerator = notid->ProjectionX("numerator",lowbin,highbin); nNotidProj = nNotid->ProjectionX("nNotidProj",lowbin,highbin); } else{ //cout<<"Getting eta dependence"<<endl; int lowbin = id->GetXaxis()->FindBin(etacut);//make sure we don't accidentally get the wrong bin int highbin; if(etacut<0.15){//then we actually have ITS standalone tracks and we only want this to run from 0.1 to 0.15 because this will be used only for ITS standalone tracks highbin = id->GetXaxis()->FindBin(0.15); } else{ highbin = id->GetXaxis()->GetNbins(); } //cout<<"Projecting from "<<id->GetXaxis()->GetBinLowEdge(lowbin)<<" to "<<id->GetXaxis()->GetBinLowEdge(highbin+1)<<endl; numerator = notid->ProjectionY("name",lowbin,highbin); denominator = id->ProjectionY("denominator",lowbin,highbin); nNotidProj = nNotid->ProjectionY("nNotidProj",lowbin,highbin); } TH1D *result = numerator; if(!denominator){ cerr<<"Uh-oh! Can't find denominator!!"; return numerator; } else{result->Divide(denominator);} if(result->GetNbinsX() != nNotidProj->GetNbinsX()){ cerr<<"Uh-oh! Can't rescale errors! "<<result->GetNbinsX()<<"!="<<nNotidProj->GetNbinsX()<<endl; return result; } if(!nNotidProj){ cerr<<"Uh-oh! Can't find histogram!!"; return numerator; } //fixing the errors for(int i=1;i<=result->GetNbinsX();i++){ Float_t value = result->GetBinContent(i); Float_t valueerr = result->GetBinError(i); Float_t n = nNotidProj->GetBinContent(i); Float_t err; if(n<=0) err = 0.0; else{err= value/TMath::Power(n,0.5);} result->SetBinError(i,err); //cout<<"Was "<<valueerr<<", setting to "<<err<<endl; } result->SetYTitle("Ratio of E_{T}^{assuming pion}/E_{T}^{real}"); result->GetYaxis()->SetTitleOffset(1.2); delete denominator; delete nNotidProj; delete notid; delete nNotid; delete id; result->SetName(name); return result; } TH1D *CorrNotID(float etacut,char *name, char *prodname, char *shortprodname, bool TPC, bool ispp, bool forSim){ gStyle->SetOptTitle(0); gStyle->SetOptStat(0); gStyle->SetOptFit(0); TCanvas *c = new TCanvas("c","c",500,400); c->SetTopMargin(0.04); c->SetRightMargin(0.04); c->SetBorderSize(0); c->SetFillColor(0); c->SetFillColor(0); c->SetBorderMode(0); c->SetFrameFillColor(0); c->SetFrameBorderMode(0); TH1D *PHOS = GetHistoCorrNotID(etacut,name,TPC,true,ispp,forSim); PHOS->SetMarkerColor(2); PHOS->SetLineColor(2); PHOS->SetAxisRange(0.0,4); if(TPC){ PHOS->SetMaximum(1.1); PHOS->SetMinimum(0.85); } else{ PHOS->SetAxisRange(0.0,0.5); } PHOS->SetMarkerStyle(20); PHOS->Draw(); TLatex *tex = new TLatex(0.161478,1.0835,prodname); tex->SetTextSize(0.0537634); tex->Draw(); char *detector = detectorEMCAL; if(etacut<0.2) detector = detectorPHOS; if(TPC){ sprintf(epsname,"pics/%s/fnotidTPC%s.eps",shortprodname,detector); sprintf(pngname,"pics/%s/fnotidTPC%s.png",shortprodname,detector); } else{ sprintf(epsname,"pics/%s/fnotidITS%s.eps",shortprodname,detector); sprintf(pngname,"pics/%s/fnotidITS%s.png",shortprodname,detector); } c->SaveAs(epsname); c->SaveAs(pngname); delete c; delete tex; PHOS->SetName(name); return PHOS; } Float_t CorrNotIDConst(float ptcut, float etacut,char *name, char *prodname, char *shortprodname, bool TPC, bool ispp, bool forSim){ if(!forSim){ if(ispp){ return 0.996; } else{ return 0.992; } } gStyle->SetOptTitle(0); gStyle->SetOptStat(0); gStyle->SetOptFit(0); TCanvas *c = new TCanvas("c","c",500,400); c->SetTopMargin(0.04); c->SetRightMargin(0.04); c->SetBorderSize(0); c->SetFillColor(0); c->SetFillColor(0); c->SetBorderMode(0); c->SetFrameFillColor(0); c->SetFrameBorderMode(0); TH1D *PHOS = GetHistoCorrNotID(ptcut,name,TPC,false,ispp,forSim); PHOS->SetMarkerColor(2); PHOS->SetLineColor(2); PHOS->SetMaximum(1.01); PHOS->SetMinimum(0.98); TF1 *func = new TF1("func","[0]",-etacut,etacut); PHOS->Fit(func,"","",-etacut,etacut); PHOS->SetMarkerStyle(20); PHOS->Draw(); TLatex *tex = new TLatex(0.161478,1.0835,prodname); tex->SetTextSize(0.0537634); tex->Draw(); char *detector = detectorEMCAL; if(etacut<0.2) detector = detectorPHOS; if(TPC){ sprintf(epsname,"pics/%s/fnotidConstTPC%s.eps",shortprodname,detector); sprintf(pngname,"pics/%s/fnotidConstTPC%s.png",shortprodname,detector); } else{ sprintf(epsname,"pics/%s/fnotidConstITS%s.eps",shortprodname,detector); sprintf(pngname,"pics/%s/fnotidConstITS%s.png",shortprodname,detector); } c->SaveAs(epsname); c->SaveAs(pngname); delete c; delete PHOS; float value = func->GetParameter(0); delete func; delete tex; return value; } //==================================CorrNoID================================================= TH1D *GetHistoNoID(float etacut, char *name, bool eta, bool TPC, bool ispp, bool forSim){ file->cd(); char *myname = mynameITS; if(TPC) myname = mynameTPCITS; TH2F *notid = ((TH2F*) out2->FindObject(Form("EtReconstructed%sChargedHadronAssumingPion",myname)))->Clone("notid"); TH2F *nNotid = ((TH2F*) out2->FindObject(Form("EtNReconstructed%sChargedHadron",myname)))->Clone("nNotid"); TH2F *id = ((TH2F*) out2->FindObject(Form("EtReconstructed%sChargedHadron",myname)))->Clone("id"); int lowbin = id->GetYaxis()->FindBin(-etacut+.001);//make sure we don't accidentally get the wrong bin int highbin = id->GetYaxis()->FindBin(etacut-.001); TH1D *nNotidProj; TH1D *denominator; TH1D *numerator; if(eta){ int lowbin = notid->GetYaxis()->FindBin(-etacut+.001);//make sure we don't accv0entally get the wrong bin int highbin = notid->GetYaxis()->FindBin(etacut-.001); //cout<<"Projecting from "<<notid->GetYaxis()->GetBinLowEdge(lowbin)<<" to "<<notid->GetYaxis()->GetBinLowEdge(highbin+1)<<endl; denominator = id->ProjectionX("name",lowbin,highbin); numerator = notid->ProjectionX("numerator",lowbin,highbin); nNotidProj = nNotid->ProjectionX("nNotidProj",lowbin,highbin); } else{ //cout<<"Getting eta dependence"<<endl; int lowbin = id->GetXaxis()->FindBin(etacut);//make sure we don't accidentally get the wrong bin int highbin; if(etacut<0.15){//then we actually have ITS standalone tracks and we only want this to run from 0.1 to 0.15 because this will be used only for ITS standalone tracks highbin = id->GetXaxis()->FindBin(0.15); } else{ highbin = id->GetXaxis()->GetNbins(); } //cout<<"Projecting from "<<id->GetXaxis()->GetBinLowEdge(lowbin)<<" to "<<id->GetXaxis()->GetBinLowEdge(highbin+1)<<endl; numerator = notid->ProjectionY("name",lowbin,highbin); denominator = id->ProjectionY("denominator",lowbin,highbin); nNotidProj = nNotid->ProjectionY("nNotidProj",lowbin,highbin); } if(denominator) numerator->Divide(denominator); if(numerator->GetNbinsX() != nNotidProj->GetNbinsX()){ cerr<<"Uh-oh! Can't rescale errors! "<<numerator->GetNbinsX()<<"!="<<nNotidProj->GetNbinsX()<<endl; return numerator; } //fixing the errors for(int i=1;i<=numerator->GetNbinsX();i++){ Float_t value = numerator->GetBinContent(i); Float_t valueerr = numerator->GetBinError(i); Float_t n = nNotidProj->GetBinContent(i); Float_t err=0.; if(n>0.0){ err = value/TMath::Power(n,0.5); } numerator->SetBinError(i,err);; } numerator->SetYTitle("Ratio of E_{T}^{assuming pion}/E_{T}^{real}"); numerator->GetYaxis()->SetTitleOffset(1.2); delete denominator; delete nNotidProj; delete notid; delete nNotid; delete id; numerator->SetName(name); return numerator; } TH1D *CorrNoID(float etacut,char *name, char *prodname, char *shortprodname, bool ispp, bool forSim){ gStyle->SetOptTitle(0); gStyle->SetOptStat(0); gStyle->SetOptFit(0); TCanvas *c = new TCanvas("c","c",500,400); c->SetTopMargin(0.04); c->SetRightMargin(0.04); c->SetBorderSize(0); c->SetFillColor(0); c->SetFillColor(0); c->SetBorderMode(0); c->SetFrameFillColor(0); c->SetFrameBorderMode(0); TH1D *PHOS = GetHistoNoID(etacut,name,true,true,ispp,forSim); PHOS->SetMarkerColor(2); PHOS->SetLineColor(2); PHOS->SetAxisRange(0.0,4); PHOS->SetMaximum(1.1); PHOS->SetMinimum(0.85); PHOS->SetMarkerStyle(20);; PHOS->Draw(); TLatex *tex = new TLatex(0.161478,1.0835,prodname); tex->SetTextSize(0.0537634); tex->Draw(); char *detector = detectorEMCAL; if(etacut<0.2) detector = detectorPHOS; sprintf(epsname,"pics/%s/fnoid%s.eps",shortprodname,detector); sprintf(pngname,"pics/%s/fnoid%s.png",shortprodname,detector); c->SaveAs(epsname); c->SaveAs(pngname); delete c; delete tex; PHOS->SetName(name); return PHOS; } Float_t CorrNoIDConst(float etacut, float ptcut,char *name, char *prodname, char *shortprodname, bool ispp, bool forSim){ gStyle->SetOptTitle(0); gStyle->SetOptStat(0); gStyle->SetOptFit(0); TCanvas *c = new TCanvas("c","c",500,400); c->SetTopMargin(0.04); c->SetRightMargin(0.04); c->SetBorderSize(0); c->SetFillColor(0); c->SetFillColor(0); c->SetBorderMode(0); c->SetFrameFillColor(0); c->SetFrameBorderMode(0); bool TPC = true; if(ptcut<.15) TPC = false; TH1D *PHOS = GetHistoNoID(ptcut,name,false,TPC,ispp,forSim); TF1 *func = new TF1("func","[0]",-etacut,etacut); PHOS->Fit(func,"","",-etacut,etacut); PHOS->SetMarkerColor(2); PHOS->SetLineColor(2); PHOS->SetMaximum(1.1); PHOS->SetMinimum(0.85); PHOS->SetMarkerStyle(20);; PHOS->Draw(); TLatex *tex = new TLatex(0.161478,1.0835,prodname); tex->SetTextSize(0.0537634); tex->Draw(); char *detector = detectorEMCAL; if(etacut<0.2) detector = detectorPHOS; if(TPC){ sprintf(epsname,"pics/%s/fnoid%sTPC.eps",shortprodname,detector); sprintf(pngname,"pics/%s/fnoid%sTPC.png",shortprodname,detector); } else{ sprintf(epsname,"pics/%s/fnoid%sITS.eps",shortprodname,detector); sprintf(pngname,"pics/%s/fnoid%sITS.png",shortprodname,detector); } c->SaveAs(epsname); c->SaveAs(pngname); delete c; delete PHOS; float value = func->GetParameter(0); delete func; delete tex; return value; } //==================================Efficiency================================================= TH1D* bayneseffdiv(TH1D* numerator, TH1D* denominator,Char_t* name) { if(!numerator){ cerr<<"Error: numerator does not exist!"<<endl; return NULL; } if(!denominator){ cerr<<"Error: denominator does not exist!"<<endl; return NULL; } TH1D* result = (TH1D*)numerator->Clone(name); Int_t nbins = numerator->GetNbinsX(); for (Int_t ibin=0; ibin<= nbins+1; ++ibin) { Double_t numeratorVal = numerator->GetBinContent(ibin); Double_t denominatorVal = denominator->GetBinContent(ibin); // Check if the errors are right or the thing is scaled Double_t numeratorValErr = numerator->GetBinError(ibin); if (!(numeratorValErr==0. || numeratorVal ==0.) ) { Double_t rescale = numeratorValErr*numeratorValErr/numeratorVal; numeratorVal /= rescale; } Double_t denominatorValErr = denominator->GetBinError(ibin); if (!(denominatorValErr==0. || denominatorVal==0. )) { Double_t rescale = denominatorValErr*denominatorValErr/denominatorVal; denominatorVal /= rescale; } Double_t quotient = 0.; if (denominatorVal!=0.) { quotient = numeratorVal/denominatorVal; } Double_t quotientErr=0; if(numeratorVal>0.0 && denominatorVal>0.0 && (denominatorVal+2.0)>0.0 && (denominatorVal+3.0) >0.0){ double argument = (numeratorVal+1.0)/(denominatorVal+2.0)* ((numeratorVal+2.0)/(denominatorVal+3.0)-(numeratorVal+1.0)/(denominatorVal+2.0)); double test = TMath::Power(TMath::Abs(argument),0.5); quotientErr = TMath::Sqrt( TMath::Abs( (numeratorVal+1.0)/(denominatorVal+2.0)* ((numeratorVal+2.0)/(denominatorVal+3.0)-(numeratorVal+1.0)/(denominatorVal+2.0)))); } result->SetBinContent(ibin,quotient); result->SetBinError(ibin,quotientErr); //cout<<"Setting bin "<<ibin<<" to "<<quotient<<" "<<numeratorVal<<"/"<<denominatorVal<<endl; } return result; } TH1D *GetHistoEfficiency(float cut, char *name, int mycase, int color, int marker,bool TPC,bool ITS, int cb, int cblast){ bool eta = true; file->cd(); char *myname = mynameITS; if(TPC && !ITS){ myname = mynameTPCITS;cout<<"I was here I should not be called"<<endl;} if(TPC&&ITS) myname = mynameTPCITS; char *mynamealt = mynameTPCITS; TH2F *numeratorParent; switch(mycase){ case 0: if(cblast != -1){//add more centrality bins for(int i=cb;i<=cblast;i++){ if(i==cb) numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiPlus",i)))->Clone("RecoHadron"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiPlus",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiMinus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KMinus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KPlus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"Proton",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"AntiProton",i))); if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiPlus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiMinus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KMinus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KPlus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"Proton",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"AntiProton",i))); } } } else{ numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"PiPlus")))->Clone("RecoHadron"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"PiMinus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"KMinus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"KPlus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"Proton"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"AntiProton"))); if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiPlus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiMinus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KMinus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KPlus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"Proton"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"AntiProton"))); } } break; case 1://pion if(cblast != -1){//add more centrality bins for(int i=cb;i<=cblast;i++){ if(i==cb) numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiPlus",i)))->Clone("RecoPion"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiPlus",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiMinus",i))); if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiMinus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiPlus",i))); } } } else{ numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"PiPlus")))->Clone("RecoPion"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"PiMinus"))); if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiPlus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiMinus"))); } } break; case 2://kaon if(cblast != -1){//add more centrality bins for(int i=cb;i<=cblast;i++){ if(i==cb) numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KPlus",i)))->Clone("RecoKaon"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KPlus",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KMinus",i))); if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KPlus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KMinus",i))); } } } else{ // cout<<"I am kaoning here !"<<endl; numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"KPlus")))->Clone("RecoKaon"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"KMinus"))); if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KPlus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KMinus"))); } //cout<<"Numerator "<<numeratorParent->GetEntries()<<endl; } break; case 3://proton if(cblast != -1){//add more centrality bins for(int i=cb;i<=cblast;i++){ if(i==cb) numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"Proton",i)))->Clone("RecoProton"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"Proton",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"AntiProton",i))); if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"Proton",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"AntiProton",i))); } } } else{ //cout<<"I am protoning here !"<<endl; numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"Proton")))->Clone("RecoProton"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"AntiProton"))); if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"Proton"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"AntiProton"))); } //cout<<"Numerator "<<numeratorParent->GetEntries()<<endl; } break; case 4://electron numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",myname,"EPlus",cbname)))->Clone("RecoElectron"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",myname,"EMinus",cbname))); if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",mynamealt,"EPlus",cbname))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",mynamealt,"EMinus",cbname))); } break; } TH2F *denominatorParent; switch(mycase){ case 0: if(cblast != -1){//add more centrality bins denominatorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNSimulatedChargedHadronCB%i",cb)))->Clone("RecoHadron"); for(int i=cb+1;i<=cblast;i++){ denominatorParent->Add((TH2F*) out2->FindObject(Form("EtNSimulatedChargedHadronCB%i",i))); } } else{ denominatorParent = (TH2F*)((TH2F*) out2->FindObject("EtNSimulatedChargedHadron"))->Clone("RecoHadron"); } break; case 1://pion if(cblast != -1){//add more centrality bins denominatorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNSimulatedPiPlusCB%i",cb)))->Clone("RecoPion"); denominatorParent->Add((TH2F*) out2->FindObject(Form("EtNSimulatedPiMinusCB%i",cb))); for(int i=cb+1;i<=cblast;i++){ denominatorParent->Add((TH2F*) out2->FindObject(Form("EtNSimulatedPiPlusCB%i",i))); denominatorParent->Add((TH2F*) out2->FindObject(Form("EtNSimulatedPiMinusCB%i",i))); } } else{ denominatorParent = (TH2F*)((TH2F*) out2->FindObject("EtNSimulatedPiPlus"))->Clone("RecoPion"); denominatorParent->Add((TH2F*) out2->FindObject("EtNSimulatedPiMinus")); } break; case 2://kaon if(cblast != -1){//add more centrality bins denominatorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNSimulatedKPlusCB%i",cb)))->Clone("RecoKaon"); denominatorParent->Add((TH2F*) out2->FindObject(Form("EtNSimulatedKMinusCB%i",cb))); for(int i=cb+1;i<=cblast;i++){ denominatorParent->Add((TH2F*) out2->FindObject(Form("EtNSimulatedKPlusCB%i",i))); denominatorParent->Add((TH2F*) out2->FindObject(Form("EtNSimulatedKMinusCB%i",i))); } } else{ denominatorParent = (TH2F*)((TH2F*) out2->FindObject("EtNSimulatedKPlus"))->Clone("RecoKaon"); denominatorParent->Add((TH2F*) out2->FindObject("EtNSimulatedKMinus")); } break; case 3://proton if(cblast != -1){//add more centrality bins for(int i=cb;i<=cblast;i++){ if(cb==i)denominatorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNSimulatedProtonCB%i",i)))->Clone("RecoProton"); else{denominatorParent->Add((TH2F*) out2->FindObject(Form("EtNSimulatedProtonCB%i",i)));} denominatorParent->Add((TH2F*) out2->FindObject(Form("EtNSimulatedAntiProtonCB%i",i))); } } else{ denominatorParent = (TH2F*)((TH2F*) out2->FindObject("EtNSimulatedProton"))->Clone("RecoProton"); denominatorParent->Add((TH2F*) out2->FindObject("EtNSimulatedAntiProton")); } break; case 4://electron denominatorParent = (TH2F*)((TH2F*) out2->FindObject("EtNSimulatedEPlus"))->Clone("RecoElectron"); denominatorParent->Add((TH2F*) out2->FindObject("EtNSimulatedEMinus")); break; } TH1D *denominator; TH1D *numerator; if(eta){ int lowbin = numeratorParent->GetYaxis()->FindBin(-cut+.001);//make sure we don't accv0entally get the wrong bin int highbin = numeratorParent->GetYaxis()->FindBin(cut-.001); //cout<<"Projecting from "<<numeratorParent->GetYaxis()->GetBinLowEdge(lowbin)<<" to "<<numeratorParent->GetYaxis()->GetBinLowEdge(highbin+1)<<endl; denominator = denominatorParent->ProjectionX(Form("garbage%s",name),lowbin,highbin); numerator = numeratorParent->ProjectionX(name,lowbin,highbin); } else{ int lowbin = denominatorParent->GetXaxis()->FindBin(cut);//make sure we don't accidentally get the wrong bin int highbin = denominatorParent->GetXaxis()->GetNbins(); //cout<<"Projecting from "<<denominatorParent->GetXaxis()->GetBinLowEdge(lowbin)<<" to "<<denominatorParent->GetXaxis()->GetBinLowEdge(highbin+1)<<endl; numerator = numeratorParent->ProjectionY(name,lowbin,highbin); denominator = denominatorParent->ProjectionY(Form("denominator%s",name),lowbin,highbin); } delete numeratorParent; delete denominatorParent; //numerator->Divide(denominator); TH1D *result = bayneseffdiv((TH1D*) numerator,(TH1D*)denominator,name); //result->Rebin(2); //result->Scale(0.5); result->SetYTitle("Efficiency"); result->GetYaxis()->SetTitleOffset(0.8); result->GetXaxis()->SetTitleOffset(0.8); result->GetYaxis()->SetLabelSize(0.05); result->GetXaxis()->SetLabelSize(0.05); result->GetYaxis()->SetTitleSize(0.05); result->GetXaxis()->SetTitleSize(0.05); result->SetMarkerColor(color); result->SetLineColor(color); result->SetMarkerStyle(marker); //result->SetName(name); //result->Draw("e"); delete denominator; delete numerator; delete numeratorParent; delete denominatorParent; result->SetName(name); return result; } void CorrEfficiencyPlots(bool TPC, char *prodname, char *shortprodname){ bool ITS = true; gStyle->SetOptTitle(0); gStyle->SetOptStat(0); gStyle->SetOptFit(0); TCanvas *c = new TCanvas("c","c",600,400); c->SetTopMargin(0.02); c->SetRightMargin(0.02); c->SetBorderSize(0); c->SetFillColor(0); c->SetFillColor(0); c->SetBorderMode(0); c->SetFrameFillColor(0); c->SetFrameBorderMode(0); //c->SetLogx(); int colortotal = 1; int colorpi = 2; int colork = 3; int colorp = 4; int phosmarker = 20; int emcalmarker = 24; float ptcut1 = 0.05; float ptcut2 = 0.1; TH1D *PHOStotal = GetHistoEfficiency(0.12,"PHOStotal",0,colortotal,phosmarker,TPC,ITS); TH1D *PHOSpi = GetHistoEfficiency(0.12,"PHOSpi",1,colorpi,phosmarker,TPC,ITS); TH1D *PHOSp = GetHistoEfficiency(0.12,"PHOSp",2,colork,phosmarker,TPC,ITS); TH1D *PHOSk = GetHistoEfficiency(0.12,"PHOSk",3,colorp,phosmarker,TPC,ITS); if(!TPC){PHOStotal->GetXaxis()->SetRange(PHOStotal->GetXaxis()->FindBin(0.05),PHOStotal->GetXaxis()->FindBin(1.0));} else{PHOStotal->GetXaxis()->SetRange(PHOStotal->GetXaxis()->FindBin(0.15),PHOStotal->GetXaxis()->FindBin(3.0));} PHOStotal->SetMinimum(0.0); PHOStotal->SetMaximum(1.0); PHOStotal->Draw(); PHOSpi->Draw("same"); PHOSp->Draw("same"); PHOSk->Draw("same"); TH1D *EMCALtotal = GetHistoEfficiency(0.7,"EMCALtotal",0,colortotal,emcalmarker,TPC,ITS); TH1D *EMCALpi = GetHistoEfficiency(0.7,"EMCALpi",1,colorpi,emcalmarker,TPC,ITS); TH1D *EMCALp = GetHistoEfficiency(0.7,"EMCALp",2,colork,emcalmarker,TPC,ITS); TH1D *EMCALk = GetHistoEfficiency(0.7,"EMCALk",3,colorp,emcalmarker,TPC,ITS); EMCALtotal->Draw("same"); EMCALpi->Draw("same"); EMCALp->Draw("same"); EMCALk->Draw("same"); TLegend *leg = new TLegend(0.22651,0.247312,0.370805,0.438172); leg->AddEntry(PHOStotal,"#pi,K,p"); leg->AddEntry(PHOSpi,"#pi^{#pm}"); leg->AddEntry(PHOSk,"K^{#pm}"); leg->AddEntry(PHOSp,"p,#bar{p}"); leg->SetFillStyle(0); leg->SetFillColor(0); leg->SetBorderSize(0); leg->SetTextSize(0.06); leg->Draw(); TLine *line = new TLine(0.2,0.0,0.2,1.0); line->Draw(); line->SetLineWidth(3.0); //line->SetLineColor(TColor::kYellow); line->SetLineStyle(2); TLatex *tex = new TLatex(0.497269,0.0513196,prodname); tex->SetTextSize(0.0537634); tex->Draw(); TLatex *tex3 = new TLatex(1.16186,0.28348,"Closed symbols |#eta|<0.12 (PHOS)"); tex3->SetTextSize(0.0537634); tex3->Draw(); TLatex *tex4 = new TLatex(1.16186,0.213221,"Open symbols |#eta|<0.70 (EMCal)"); tex4->SetTextSize(0.0537634); tex4->Draw(); TLatex *tex2 = new TLatex(0.241937,0.448436,"Likely TPC cut-off 200 MeV/c"); tex2->SetTextSize(0.0537634); tex2->Draw(); if(TPC){ if(ITS){ sprintf(epsname,"pics/%s/CorrEfficiencyITSTPC.eps",shortprodname); sprintf(pngname,"pics/%s/CorrEfficiencyITSTPC.png",shortprodname); } else{ sprintf(epsname,"pics/%s/CorrEfficiencyTPC.eps",shortprodname); sprintf(pngname,"pics/%s/CorrEfficiencyTPC.png",shortprodname); } } else{ sprintf(epsname,"pics/%s/CorrEfficiencyITS.eps",shortprodname); sprintf(pngname,"pics/%s/CorrEfficiencyITS.png",shortprodname); } delete PHOStotal; delete PHOSpi; delete PHOSp; delete PHOSk; delete EMCALtotal; delete EMCALpi; delete EMCALp; delete EMCALk; delete leg; c->SaveAs(epsname); c->SaveAs(pngname); delete c; delete line; delete tex; delete tex3; } //==================================CorrBkgd================================================= TH1D *GetHistoCorrBkgd(float etacut,char *name, bool TPC,bool ispp,bool forSim){ file->cd(); char *reweightname = reweightedNo; if(!forSim) reweightname = reweightedYes; char *myname = mynameITS; if(TPC) myname = mynameTPCITS; TH2F *signal = ((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiPlus",myname)))->Clone("signal"); signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiMinus",myname))); signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedKMinus",myname))); signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedKPlus",myname))); signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedProton",myname))); signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedAntiProton",myname))); signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sUnidentifiedAssumingPion",myname))); //Et of all unidentified hadrons (plus hadrons identified as pions) calculated assuming their true mass TH2F *bkgd = ((TH2F*) out2->FindObject(Form("EtReconstructed%sMisidentifiedElectrons",myname)))->Clone("bkgd"); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sLambdaDaughters%s",myname,reweightname))); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiLambdaDaughters%s",myname,reweightname))); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sK0SDaughters%s",myname,reweightname))); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sXiDaughters",myname))); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiXiDaughters",myname))); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sOmegaDaughters",myname))); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiOmegaDaughters",myname))); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sConversionElectrons",myname)) ); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryMuons",myname)) ); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryPions",myname)) ); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryProtons",myname)) ); int lowbin = bkgd->GetYaxis()->FindBin(-etacut+.001);//make sure we don't accidentally get the wrong bin int highbin = bkgd->GetYaxis()->FindBin(etacut-.001); //cout<<"Projecting from "<<bkgd->GetYaxis()->GetBinLowEdge(lowbin)<<" to "<<bkgd->GetYaxis()->GetBinLowEdge(highbin+1)<<endl; TH1D *denominator = signal->ProjectionX(name,lowbin,highbin); TH1D *numerator = bkgd->ProjectionX("numerator",lowbin,highbin); denominator->Add(numerator); numerator->Divide(denominator); numerator->SetYTitle("Ratio of E_{T}^{background}/E_{T}^{real}"); numerator->GetYaxis()->SetTitleOffset(1.2); delete signal; delete bkgd; delete denominator; numerator->SetName(name); return numerator; } void CorrBkgdPlots(char *prodname, char *shortprodname, bool TPC,bool ispp,bool forSim){ gStyle->SetOptTitle(0); gStyle->SetOptStat(0); gStyle->SetOptFit(0); TCanvas *c = new TCanvas("c","c",500,400); c->SetTopMargin(0.04); c->SetRightMargin(0.04); c->SetBorderSize(0); c->SetFillColor(0); c->SetFillColor(0); c->SetBorderMode(0); c->SetFrameFillColor(0); c->SetFrameBorderMode(0); TH1D *PHOS = GetHistoCorrBkgd(0.12,"PHOS2",TPC,ispp,forSim); TH1D *EMCAL = GetHistoCorrBkgd(0.7,"EMCAL2",TPC,ispp,forSim); PHOS->SetMarkerColor(2); EMCAL->SetMarkerColor(4); PHOS->SetLineColor(2); EMCAL->SetLineColor(4); //EMCAL->SetLineWidth(2); //PHOS->SetAxisRange(0.0,4); PHOS->SetMaximum(0.2); PHOS->SetMinimum(0.0); PHOS->SetMarkerStyle(20);; EMCAL->SetMarkerStyle(21); // TF1 *funcEMCAL = new TF1("funcEMCAL","[0]+0.0*x",0.05,4); // funcEMCAL->SetParameter(0,0.95); // funcEMCAL->SetParLimits(0,0.9,1.1); //EMCAL->Fit(funcEMCAL);//,"","",0.05,3.0); // TF1 *funcPHOS = new TF1("funcPHOS","[0]+0.0*x",0.05,4); // funcPHOS->SetParameter(0,1.0); //PHOS->Fit(funcPHOS); if(TPC) PHOS->GetXaxis()->SetRange(PHOS->GetXaxis()->FindBin(0.0),PHOS->GetXaxis()->FindBin(4.)); else{ PHOS->GetXaxis()->SetRange(PHOS->GetXaxis()->FindBin(0.0),PHOS->GetXaxis()->FindBin(1.));} PHOS->Draw(); EMCAL->Draw("same"); TLatex *tex = new TLatex(0.161478,1.0835,prodname); tex->SetTextSize(0.0537634); tex->Draw(); TLegend *leg = new TLegend(0.145161,0.604839,0.40121,0.860215); leg->AddEntry(PHOS,"|#eta|<0.12"); leg->AddEntry(EMCAL,"|#eta|<0.70"); leg->SetFillStyle(0); leg->SetFillColor(0); leg->SetBorderSize(0); leg->Draw(); if(TPC){ c->SaveAs(Form("pics/%s/bkgdTPC.eps",shortprodname)); c->SaveAs(Form("pics/%s/bkgdTPC.png",shortprodname)); } else{ c->SaveAs(Form("pics/%s/bkgdITS.eps",shortprodname)); c->SaveAs(Form("pics/%s/bkgdITS.png",shortprodname)); } delete c; delete PHOS; delete EMCAL; delete tex; delete leg; } //TH1D *ppBkgd() TH1D *pp276TPCBkgd(){ //=========Macro generated from canvas: c4/c4 //========= (Thu May 5 23:07:14 2011) by ROOT version5.28/00c TH1D *Allpt7TeVScaling = new TH1D("Allpt7TeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); Allpt7TeVScaling->SetBinContent(13,0.1673787); Allpt7TeVScaling->SetBinContent(14,0.1548408); Allpt7TeVScaling->SetBinContent(15,0.1328529); Allpt7TeVScaling->SetBinContent(16,0.1220669); Allpt7TeVScaling->SetBinContent(17,0.1109226); Allpt7TeVScaling->SetBinContent(18,0.0981996); Allpt7TeVScaling->SetBinContent(19,0.09152354); Allpt7TeVScaling->SetBinContent(20,0.07937801); Allpt7TeVScaling->SetBinContent(21,0.07163593); Allpt7TeVScaling->SetBinContent(22,0.0562236); Allpt7TeVScaling->SetBinContent(23,0.04398756); Allpt7TeVScaling->SetBinContent(24,0.039459); Allpt7TeVScaling->SetBinContent(25,0.03618901); Allpt7TeVScaling->SetBinContent(26,0.0322147); Allpt7TeVScaling->SetBinContent(27,0.03294796); Allpt7TeVScaling->SetBinContent(28,0.03114461); Allpt7TeVScaling->SetBinContent(29,0.03039579); Allpt7TeVScaling->SetBinContent(30,0.03261539); Allpt7TeVScaling->SetBinContent(31,0.03239832); Allpt7TeVScaling->SetBinContent(32,0.03055859); Allpt7TeVScaling->SetBinContent(33,0.03181271); Allpt7TeVScaling->SetBinContent(34,0.02882665); Allpt7TeVScaling->SetBinContent(35,0.03102149); Allpt7TeVScaling->SetBinContent(36,0.0272696); Allpt7TeVScaling->SetBinContent(37,0.02783727); Allpt7TeVScaling->SetBinContent(38,0.02334413); Allpt7TeVScaling->SetBinContent(39,0.02357657); Allpt7TeVScaling->SetBinContent(40,0.02630686); Allpt7TeVScaling->SetBinContent(41,0.02372341); Allpt7TeVScaling->SetBinContent(42,0.02279457); Allpt7TeVScaling->SetBinContent(43,0.01999143); Allpt7TeVScaling->SetBinContent(44,0.02019579); Allpt7TeVScaling->SetBinContent(45,0.02029719); Allpt7TeVScaling->SetBinContent(46,0.0166231); Allpt7TeVScaling->SetBinContent(47,0.01824842); Allpt7TeVScaling->SetBinContent(48,0.0188995); Allpt7TeVScaling->SetBinContent(49,0.01836067); Allpt7TeVScaling->SetBinContent(50,0.02113978); Allpt7TeVScaling->SetBinContent(51,0.02236854); Allpt7TeVScaling->SetBinContent(52,0.0204874); Allpt7TeVScaling->SetBinContent(53,0.02236788); Allpt7TeVScaling->SetBinContent(54,0.01862457); Allpt7TeVScaling->SetBinContent(55,0.0215448); Allpt7TeVScaling->SetBinContent(56,0.01757402); Allpt7TeVScaling->SetBinContent(57,0.01899477); Allpt7TeVScaling->SetBinContent(58,0.01766989); Allpt7TeVScaling->SetBinContent(59,0.01984435); Allpt7TeVScaling->SetBinContent(60,0.01842092); Allpt7TeVScaling->SetBinContent(61,0.01668053); Allpt7TeVScaling->SetBinContent(62,0.01708603); Allpt7TeVScaling->SetBinContent(63,0.01557306); Allpt7TeVScaling->SetBinContent(64,0.01878942); Allpt7TeVScaling->SetBinContent(65,0.01458043); Allpt7TeVScaling->SetBinContent(66,0.0151346); Allpt7TeVScaling->SetBinContent(67,0.01565407); Allpt7TeVScaling->SetBinContent(68,0.01327014); Allpt7TeVScaling->SetBinContent(69,0.01284378); Allpt7TeVScaling->SetBinContent(70,0.01046102); Allpt7TeVScaling->SetBinContent(71,0.009963117); Allpt7TeVScaling->SetBinContent(72,0.01359609); Allpt7TeVScaling->SetBinContent(73,0.01088263); Allpt7TeVScaling->SetBinContent(74,0.01317689); Allpt7TeVScaling->SetBinContent(75,0.01433201); Allpt7TeVScaling->SetBinContent(76,0.01215939); Allpt7TeVScaling->SetBinContent(77,0.01380766); Allpt7TeVScaling->SetBinContent(78,0.01499505); Allpt7TeVScaling->SetBinContent(79,0.015013); Allpt7TeVScaling->SetBinContent(80,0.01140629); Allpt7TeVScaling->SetBinContent(81,0.009461979); Allpt7TeVScaling->SetBinContent(82,0.007775343); Allpt7TeVScaling->SetBinContent(83,0.01092198); Allpt7TeVScaling->SetBinContent(84,0.009924767); Allpt7TeVScaling->SetBinContent(85,0.006399066); Allpt7TeVScaling->SetBinContent(86,0.008370865); Allpt7TeVScaling->SetBinContent(87,0.01115192); Allpt7TeVScaling->SetBinContent(88,0.003781796); Allpt7TeVScaling->SetBinContent(89,0.01596916); Allpt7TeVScaling->SetBinContent(90,0.005959723); Allpt7TeVScaling->SetBinContent(91,0.006885685); Allpt7TeVScaling->SetBinContent(93,0.006011974); Allpt7TeVScaling->SetBinContent(94,0.008862645); Allpt7TeVScaling->SetBinContent(95,0.009963223); Allpt7TeVScaling->SetBinContent(96,0.01457121); Allpt7TeVScaling->SetBinContent(97,0.007500135); Allpt7TeVScaling->SetBinContent(98,0.02636933); Allpt7TeVScaling->SetBinContent(103,0.06353767); Allpt7TeVScaling->SetBinError(13,0.02859277); Allpt7TeVScaling->SetBinError(14,0.006796474); Allpt7TeVScaling->SetBinError(15,0.004234415); Allpt7TeVScaling->SetBinError(16,0.0036414); Allpt7TeVScaling->SetBinError(17,0.003333835); Allpt7TeVScaling->SetBinError(18,0.003009034); Allpt7TeVScaling->SetBinError(19,0.002839933); Allpt7TeVScaling->SetBinError(20,0.002588262); Allpt7TeVScaling->SetBinError(21,0.001703954); Allpt7TeVScaling->SetBinError(22,0.001495734); Allpt7TeVScaling->SetBinError(23,0.001297475); Allpt7TeVScaling->SetBinError(24,0.001213434); Allpt7TeVScaling->SetBinError(25,0.001162398); Allpt7TeVScaling->SetBinError(26,0.001137158); Allpt7TeVScaling->SetBinError(27,0.001208358); Allpt7TeVScaling->SetBinError(28,0.001183675); Allpt7TeVScaling->SetBinError(29,0.001263691); Allpt7TeVScaling->SetBinError(30,0.001320041); Allpt7TeVScaling->SetBinError(31,0.001362399); Allpt7TeVScaling->SetBinError(32,0.001338259); Allpt7TeVScaling->SetBinError(33,0.001369813); Allpt7TeVScaling->SetBinError(34,0.001358652); Allpt7TeVScaling->SetBinError(35,0.00142864); Allpt7TeVScaling->SetBinError(36,0.001378454); Allpt7TeVScaling->SetBinError(37,0.001384729); Allpt7TeVScaling->SetBinError(38,0.001257903); Allpt7TeVScaling->SetBinError(39,0.001295459); Allpt7TeVScaling->SetBinError(40,0.001546992); Allpt7TeVScaling->SetBinError(41,0.001496938); Allpt7TeVScaling->SetBinError(42,0.001486813); Allpt7TeVScaling->SetBinError(43,0.00140119); Allpt7TeVScaling->SetBinError(44,0.001478331); Allpt7TeVScaling->SetBinError(45,0.001566594); Allpt7TeVScaling->SetBinError(46,0.001407568); Allpt7TeVScaling->SetBinError(47,0.001509742); Allpt7TeVScaling->SetBinError(48,0.001583129); Allpt7TeVScaling->SetBinError(49,0.001668789); Allpt7TeVScaling->SetBinError(50,0.001915902); Allpt7TeVScaling->SetBinError(51,0.001968585); Allpt7TeVScaling->SetBinError(52,0.00195473); Allpt7TeVScaling->SetBinError(53,0.00207896); Allpt7TeVScaling->SetBinError(54,0.001937327); Allpt7TeVScaling->SetBinError(55,0.002108783); Allpt7TeVScaling->SetBinError(56,0.001816997); Allpt7TeVScaling->SetBinError(57,0.001892795); Allpt7TeVScaling->SetBinError(58,0.001803616); Allpt7TeVScaling->SetBinError(59,0.001994761); Allpt7TeVScaling->SetBinError(60,0.001964909); Allpt7TeVScaling->SetBinError(61,0.001255034); Allpt7TeVScaling->SetBinError(62,0.001312232); Allpt7TeVScaling->SetBinError(63,0.001358296); Allpt7TeVScaling->SetBinError(64,0.001569736); Allpt7TeVScaling->SetBinError(65,0.001469072); Allpt7TeVScaling->SetBinError(66,0.001588803); Allpt7TeVScaling->SetBinError(67,0.001683862); Allpt7TeVScaling->SetBinError(68,0.00164699); Allpt7TeVScaling->SetBinError(69,0.001760645); Allpt7TeVScaling->SetBinError(70,0.00160423); Allpt7TeVScaling->SetBinError(71,0.001693146); Allpt7TeVScaling->SetBinError(72,0.002042705); Allpt7TeVScaling->SetBinError(73,0.001965856); Allpt7TeVScaling->SetBinError(74,0.002245841); Allpt7TeVScaling->SetBinError(75,0.002520926); Allpt7TeVScaling->SetBinError(76,0.002402524); Allpt7TeVScaling->SetBinError(77,0.002627656); Allpt7TeVScaling->SetBinError(78,0.002963027); Allpt7TeVScaling->SetBinError(79,0.003088673); Allpt7TeVScaling->SetBinError(80,0.002958055); Allpt7TeVScaling->SetBinError(81,0.001434729); Allpt7TeVScaling->SetBinError(82,0.001531537); Allpt7TeVScaling->SetBinError(83,0.002200666); Allpt7TeVScaling->SetBinError(84,0.00242244); Allpt7TeVScaling->SetBinError(85,0.002270338); Allpt7TeVScaling->SetBinError(86,0.002972849); Allpt7TeVScaling->SetBinError(87,0.003966084); Allpt7TeVScaling->SetBinError(88,0.002679263); Allpt7TeVScaling->SetBinError(89,0.006083804); Allpt7TeVScaling->SetBinError(90,0.004227019); Allpt7TeVScaling->SetBinError(91,0.004886052); Allpt7TeVScaling->SetBinError(93,0.006030514); Allpt7TeVScaling->SetBinError(94,0.008892003); Allpt7TeVScaling->SetBinError(95,0.01001297); Allpt7TeVScaling->SetBinError(96,0.008480241); Allpt7TeVScaling->SetBinError(97,0.007527505); Allpt7TeVScaling->SetBinError(98,0.01889477); Allpt7TeVScaling->SetBinError(103,0.065502); Allpt7TeVScaling->SetMinimum(0); Allpt7TeVScaling->SetMaximum(0.04); Allpt7TeVScaling->SetEntries(1049.145); Allpt7TeVScaling->SetStats(0); Allpt7TeVScaling->SetMarkerStyle(20); Allpt7TeVScaling->GetXaxis()->SetTitle("p_{T}"); Allpt7TeVScaling->GetXaxis()->SetRange(1,91); Allpt7TeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); Allpt7TeVScaling->GetYaxis()->SetTitleOffset(1.2); TH1D *Allpt900GeVScaling = new TH1D("Allpt900GeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); Allpt900GeVScaling->SetBinContent(13,0.1950247); Allpt900GeVScaling->SetBinContent(14,0.147561); Allpt900GeVScaling->SetBinContent(15,0.136378); Allpt900GeVScaling->SetBinContent(16,0.1233744); Allpt900GeVScaling->SetBinContent(17,0.1126807); Allpt900GeVScaling->SetBinContent(18,0.09909555); Allpt900GeVScaling->SetBinContent(19,0.09191239); Allpt900GeVScaling->SetBinContent(20,0.08438053); Allpt900GeVScaling->SetBinContent(21,0.07119319); Allpt900GeVScaling->SetBinContent(22,0.05773434); Allpt900GeVScaling->SetBinContent(23,0.04292845); Allpt900GeVScaling->SetBinContent(24,0.0396127); Allpt900GeVScaling->SetBinContent(25,0.03395077); Allpt900GeVScaling->SetBinContent(26,0.03305251); Allpt900GeVScaling->SetBinContent(27,0.03297093); Allpt900GeVScaling->SetBinContent(28,0.031215); Allpt900GeVScaling->SetBinContent(29,0.03178372); Allpt900GeVScaling->SetBinContent(30,0.03236985); Allpt900GeVScaling->SetBinContent(31,0.03271276); Allpt900GeVScaling->SetBinContent(32,0.02993124); Allpt900GeVScaling->SetBinContent(33,0.03004746); Allpt900GeVScaling->SetBinContent(34,0.03077143); Allpt900GeVScaling->SetBinContent(35,0.03211214); Allpt900GeVScaling->SetBinContent(36,0.02807512); Allpt900GeVScaling->SetBinContent(37,0.02651954); Allpt900GeVScaling->SetBinContent(38,0.02503885); Allpt900GeVScaling->SetBinContent(39,0.02325777); Allpt900GeVScaling->SetBinContent(40,0.02773712); Allpt900GeVScaling->SetBinContent(41,0.02491155); Allpt900GeVScaling->SetBinContent(42,0.02389028); Allpt900GeVScaling->SetBinContent(43,0.02085532); Allpt900GeVScaling->SetBinContent(44,0.01959721); Allpt900GeVScaling->SetBinContent(45,0.01937255); Allpt900GeVScaling->SetBinContent(46,0.01745646); Allpt900GeVScaling->SetBinContent(47,0.01890853); Allpt900GeVScaling->SetBinContent(48,0.02053692); Allpt900GeVScaling->SetBinContent(49,0.01813791); Allpt900GeVScaling->SetBinContent(50,0.02187495); Allpt900GeVScaling->SetBinContent(51,0.02409687); Allpt900GeVScaling->SetBinContent(52,0.0206565); Allpt900GeVScaling->SetBinContent(53,0.02269165); Allpt900GeVScaling->SetBinContent(54,0.01850795); Allpt900GeVScaling->SetBinContent(55,0.02115427); Allpt900GeVScaling->SetBinContent(56,0.01543302); Allpt900GeVScaling->SetBinContent(57,0.01642095); Allpt900GeVScaling->SetBinContent(58,0.01542331); Allpt900GeVScaling->SetBinContent(59,0.01931254); Allpt900GeVScaling->SetBinContent(60,0.01650654); Allpt900GeVScaling->SetBinContent(61,0.01827512); Allpt900GeVScaling->SetBinContent(62,0.01700125); Allpt900GeVScaling->SetBinContent(63,0.01544321); Allpt900GeVScaling->SetBinContent(64,0.01830588); Allpt900GeVScaling->SetBinContent(65,0.01584369); Allpt900GeVScaling->SetBinContent(66,0.0145421); Allpt900GeVScaling->SetBinContent(67,0.01382206); Allpt900GeVScaling->SetBinContent(68,0.01358825); Allpt900GeVScaling->SetBinContent(69,0.01294888); Allpt900GeVScaling->SetBinContent(70,0.009346655); Allpt900GeVScaling->SetBinContent(71,0.01235781); Allpt900GeVScaling->SetBinContent(72,0.0127376); Allpt900GeVScaling->SetBinContent(73,0.008968166); Allpt900GeVScaling->SetBinContent(74,0.01467765); Allpt900GeVScaling->SetBinContent(75,0.01432309); Allpt900GeVScaling->SetBinContent(76,0.01204856); Allpt900GeVScaling->SetBinContent(77,0.01436071); Allpt900GeVScaling->SetBinContent(78,0.01451222); Allpt900GeVScaling->SetBinContent(79,0.01378967); Allpt900GeVScaling->SetBinContent(80,0.01040148); Allpt900GeVScaling->SetBinContent(81,0.01041213); Allpt900GeVScaling->SetBinContent(82,0.009971816); Allpt900GeVScaling->SetBinContent(83,0.01417668); Allpt900GeVScaling->SetBinContent(84,0.008713292); Allpt900GeVScaling->SetBinContent(85,0.005425386); Allpt900GeVScaling->SetBinContent(86,0.005303683); Allpt900GeVScaling->SetBinContent(87,0.006759432); Allpt900GeVScaling->SetBinContent(88,0.004445477); Allpt900GeVScaling->SetBinContent(89,0.01810702); Allpt900GeVScaling->SetBinContent(90,0.0108909); Allpt900GeVScaling->SetBinContent(91,0.01186746); Allpt900GeVScaling->SetBinContent(92,0.006128581); Allpt900GeVScaling->SetBinContent(93,0.007452721); Allpt900GeVScaling->SetBinContent(94,0.0259471); Allpt900GeVScaling->SetBinContent(95,0.0114321); Allpt900GeVScaling->SetBinContent(98,0.03264295); Allpt900GeVScaling->SetBinContent(103,0.07829265); Allpt900GeVScaling->SetBinError(13,0.03504783); Allpt900GeVScaling->SetBinError(14,0.007336667); Allpt900GeVScaling->SetBinError(15,0.004765958); Allpt900GeVScaling->SetBinError(16,0.004028143); Allpt900GeVScaling->SetBinError(17,0.003760166); Allpt900GeVScaling->SetBinError(18,0.003357858); Allpt900GeVScaling->SetBinError(19,0.003156502); Allpt900GeVScaling->SetBinError(20,0.002970822); Allpt900GeVScaling->SetBinError(21,0.001886555); Allpt900GeVScaling->SetBinError(22,0.001676045); Allpt900GeVScaling->SetBinError(23,0.001411794); Allpt900GeVScaling->SetBinError(24,0.001346425); Allpt900GeVScaling->SetBinError(25,0.001254258); Allpt900GeVScaling->SetBinError(26,0.001299384); Allpt900GeVScaling->SetBinError(27,0.001328066); Allpt900GeVScaling->SetBinError(28,0.001330088); Allpt900GeVScaling->SetBinError(29,0.001446767); Allpt900GeVScaling->SetBinError(30,0.00145402); Allpt900GeVScaling->SetBinError(31,0.001463994); Allpt900GeVScaling->SetBinError(32,0.001449209); Allpt900GeVScaling->SetBinError(33,0.001414847); Allpt900GeVScaling->SetBinError(34,0.001598429); Allpt900GeVScaling->SetBinError(35,0.001590854); Allpt900GeVScaling->SetBinError(36,0.001526798); Allpt900GeVScaling->SetBinError(37,0.001486126); Allpt900GeVScaling->SetBinError(38,0.001457142); Allpt900GeVScaling->SetBinError(39,0.001474238); Allpt900GeVScaling->SetBinError(40,0.001757902); Allpt900GeVScaling->SetBinError(41,0.001712601); Allpt900GeVScaling->SetBinError(42,0.001671385); Allpt900GeVScaling->SetBinError(43,0.001570556); Allpt900GeVScaling->SetBinError(44,0.001627725); Allpt900GeVScaling->SetBinError(45,0.001637262); Allpt900GeVScaling->SetBinError(46,0.001641227); Allpt900GeVScaling->SetBinError(47,0.001716345); Allpt900GeVScaling->SetBinError(48,0.001907778); Allpt900GeVScaling->SetBinError(49,0.001842589); Allpt900GeVScaling->SetBinError(50,0.00215445); Allpt900GeVScaling->SetBinError(51,0.002310602); Allpt900GeVScaling->SetBinError(52,0.002171046); Allpt900GeVScaling->SetBinError(53,0.002386521); Allpt900GeVScaling->SetBinError(54,0.002149237); Allpt900GeVScaling->SetBinError(55,0.002233898); Allpt900GeVScaling->SetBinError(56,0.001856697); Allpt900GeVScaling->SetBinError(57,0.001891742); Allpt900GeVScaling->SetBinError(58,0.001840334); Allpt900GeVScaling->SetBinError(59,0.002168741); Allpt900GeVScaling->SetBinError(60,0.002044863); Allpt900GeVScaling->SetBinError(61,0.001431972); Allpt900GeVScaling->SetBinError(62,0.001452706); Allpt900GeVScaling->SetBinError(63,0.001477499); Allpt900GeVScaling->SetBinError(64,0.001722091); Allpt900GeVScaling->SetBinError(65,0.001674803); Allpt900GeVScaling->SetBinError(66,0.001737491); Allpt900GeVScaling->SetBinError(67,0.001783912); Allpt900GeVScaling->SetBinError(68,0.00186443); Allpt900GeVScaling->SetBinError(69,0.001965873); Allpt900GeVScaling->SetBinError(70,0.001687347); Allpt900GeVScaling->SetBinError(71,0.002072538); Allpt900GeVScaling->SetBinError(72,0.002233273); Allpt900GeVScaling->SetBinError(73,0.001965905); Allpt900GeVScaling->SetBinError(74,0.002669571); Allpt900GeVScaling->SetBinError(75,0.00278613); Allpt900GeVScaling->SetBinError(76,0.002661412); Allpt900GeVScaling->SetBinError(77,0.002952605); Allpt900GeVScaling->SetBinError(78,0.003268721); Allpt900GeVScaling->SetBinError(79,0.003273694); Allpt900GeVScaling->SetBinError(80,0.002900898); Allpt900GeVScaling->SetBinError(81,0.001677933); Allpt900GeVScaling->SetBinError(82,0.001933408); Allpt900GeVScaling->SetBinError(83,0.002759328); Allpt900GeVScaling->SetBinError(84,0.002526604); Allpt900GeVScaling->SetBinError(85,0.002297818); Allpt900GeVScaling->SetBinError(86,0.002659213); Allpt900GeVScaling->SetBinError(87,0.003392141); Allpt900GeVScaling->SetBinError(88,0.003150511); Allpt900GeVScaling->SetBinError(89,0.007499872); Allpt900GeVScaling->SetBinError(90,0.006322856); Allpt900GeVScaling->SetBinError(91,0.00689283); Allpt900GeVScaling->SetBinError(92,0.006146806); Allpt900GeVScaling->SetBinError(93,0.007481131); Allpt900GeVScaling->SetBinError(94,0.01853139); Allpt900GeVScaling->SetBinError(95,0.01149767); Allpt900GeVScaling->SetBinError(98,0.02345991); Allpt900GeVScaling->SetBinError(103,0.08126059); Allpt900GeVScaling->SetMinimum(0); Allpt900GeVScaling->SetMaximum(0.04); Allpt900GeVScaling->SetEntries(724.3671); Allpt900GeVScaling->SetStats(0); Allpt900GeVScaling->SetMarkerStyle(20); Allpt900GeVScaling->GetXaxis()->SetTitle("p_{T}"); Allpt900GeVScaling->GetXaxis()->SetRange(1,91); Allpt900GeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); Allpt900GeVScaling->GetYaxis()->SetTitleOffset(1.2); Allpt7TeVScaling->Add(Allpt900GeVScaling); Allpt7TeVScaling->Scale(0.5); delete Allpt900GeVScaling; return Allpt7TeVScaling; } TH1D *PbPb276TPCBkgd(){ TH1D *Allpt7TeVScaling = new TH1D("Allpt7TeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); Allpt7TeVScaling->SetBinContent(13,0.1745755); Allpt7TeVScaling->SetBinContent(14,0.1624186); Allpt7TeVScaling->SetBinContent(15,0.1431048); Allpt7TeVScaling->SetBinContent(16,0.1331614); Allpt7TeVScaling->SetBinContent(17,0.121262); Allpt7TeVScaling->SetBinContent(18,0.1082923); Allpt7TeVScaling->SetBinContent(19,0.1011451); Allpt7TeVScaling->SetBinContent(20,0.08875923); Allpt7TeVScaling->SetBinContent(21,0.08100933); Allpt7TeVScaling->SetBinContent(22,0.06527292); Allpt7TeVScaling->SetBinContent(23,0.05220065); Allpt7TeVScaling->SetBinContent(24,0.04752904); Allpt7TeVScaling->SetBinContent(25,0.04421568); Allpt7TeVScaling->SetBinContent(26,0.03879551); Allpt7TeVScaling->SetBinContent(27,0.03944679); Allpt7TeVScaling->SetBinContent(28,0.03707215); Allpt7TeVScaling->SetBinContent(29,0.03663599); Allpt7TeVScaling->SetBinContent(30,0.03937437); Allpt7TeVScaling->SetBinContent(31,0.03888919); Allpt7TeVScaling->SetBinContent(32,0.03646502); Allpt7TeVScaling->SetBinContent(33,0.03784992); Allpt7TeVScaling->SetBinContent(34,0.03360079); Allpt7TeVScaling->SetBinContent(35,0.03761566); Allpt7TeVScaling->SetBinContent(36,0.03276653); Allpt7TeVScaling->SetBinContent(37,0.03401566); Allpt7TeVScaling->SetBinContent(38,0.02803893); Allpt7TeVScaling->SetBinContent(39,0.0290018); Allpt7TeVScaling->SetBinContent(40,0.03323201); Allpt7TeVScaling->SetBinContent(41,0.03006077); Allpt7TeVScaling->SetBinContent(42,0.02867674); Allpt7TeVScaling->SetBinContent(43,0.02471421); Allpt7TeVScaling->SetBinContent(44,0.02610375); Allpt7TeVScaling->SetBinContent(45,0.02601621); Allpt7TeVScaling->SetBinContent(46,0.02184781); Allpt7TeVScaling->SetBinContent(47,0.02479646); Allpt7TeVScaling->SetBinContent(48,0.02519956); Allpt7TeVScaling->SetBinContent(49,0.02376095); Allpt7TeVScaling->SetBinContent(50,0.02836189); Allpt7TeVScaling->SetBinContent(51,0.0339595); Allpt7TeVScaling->SetBinContent(52,0.03057126); Allpt7TeVScaling->SetBinContent(53,0.03238621); Allpt7TeVScaling->SetBinContent(54,0.0264272); Allpt7TeVScaling->SetBinContent(55,0.03203146); Allpt7TeVScaling->SetBinContent(56,0.02554187); Allpt7TeVScaling->SetBinContent(57,0.02901947); Allpt7TeVScaling->SetBinContent(58,0.02584016); Allpt7TeVScaling->SetBinContent(59,0.03055042); Allpt7TeVScaling->SetBinContent(60,0.02658193); Allpt7TeVScaling->SetBinContent(61,0.02553988); Allpt7TeVScaling->SetBinContent(62,0.02787508); Allpt7TeVScaling->SetBinContent(63,0.02836785); Allpt7TeVScaling->SetBinContent(64,0.03449766); Allpt7TeVScaling->SetBinContent(65,0.02710397); Allpt7TeVScaling->SetBinContent(66,0.02741966); Allpt7TeVScaling->SetBinContent(67,0.03184318); Allpt7TeVScaling->SetBinContent(68,0.02993716); Allpt7TeVScaling->SetBinContent(69,0.03113033); Allpt7TeVScaling->SetBinContent(70,0.01946293); Allpt7TeVScaling->SetBinContent(71,0.02890383); Allpt7TeVScaling->SetBinContent(72,0.02602078); Allpt7TeVScaling->SetBinContent(73,0.02063028); Allpt7TeVScaling->SetBinContent(74,0.03186679); Allpt7TeVScaling->SetBinContent(75,0.0399062); Allpt7TeVScaling->SetBinContent(76,0.02051914); Allpt7TeVScaling->SetBinContent(77,0.04454688); Allpt7TeVScaling->SetBinContent(78,0.04798449); Allpt7TeVScaling->SetBinContent(79,0.0311417); Allpt7TeVScaling->SetBinContent(80,0.03774521); Allpt7TeVScaling->SetBinContent(81,0.0219128); Allpt7TeVScaling->SetBinContent(82,0.02323492); Allpt7TeVScaling->SetBinContent(83,0.02412182); Allpt7TeVScaling->SetBinContent(84,0.03636555); Allpt7TeVScaling->SetBinContent(85,0.007488934); Allpt7TeVScaling->SetBinContent(86,0.01676086); Allpt7TeVScaling->SetBinContent(87,0.02427021); Allpt7TeVScaling->SetBinContent(88,0.02482882); Allpt7TeVScaling->SetBinContent(89,0.03680926); Allpt7TeVScaling->SetBinContent(90,0.01600685); Allpt7TeVScaling->SetBinContent(91,0.006885685); Allpt7TeVScaling->SetBinContent(93,0.0195411); Allpt7TeVScaling->SetBinContent(94,0.008862645); Allpt7TeVScaling->SetBinContent(95,0.009963223); Allpt7TeVScaling->SetBinContent(96,0.01399091); Allpt7TeVScaling->SetBinContent(97,0.007500135); Allpt7TeVScaling->SetBinContent(98,0.01922651); Allpt7TeVScaling->SetBinContent(103,0.03511665); Allpt7TeVScaling->SetBinError(13,0.03042345); Allpt7TeVScaling->SetBinError(14,0.007201474); Allpt7TeVScaling->SetBinError(15,0.004623859); Allpt7TeVScaling->SetBinError(16,0.004039713); Allpt7TeVScaling->SetBinError(17,0.003694238); Allpt7TeVScaling->SetBinError(18,0.003377785); Allpt7TeVScaling->SetBinError(19,0.00319956); Allpt7TeVScaling->SetBinError(20,0.00294949); Allpt7TeVScaling->SetBinError(21,0.001965652); Allpt7TeVScaling->SetBinError(22,0.001771706); Allpt7TeVScaling->SetBinError(23,0.0015726); Allpt7TeVScaling->SetBinError(24,0.001495611); Allpt7TeVScaling->SetBinError(25,0.00145454); Allpt7TeVScaling->SetBinError(26,0.001398941); Allpt7TeVScaling->SetBinError(27,0.001495129); Allpt7TeVScaling->SetBinError(28,0.001474393); Allpt7TeVScaling->SetBinError(29,0.001608602); Allpt7TeVScaling->SetBinError(30,0.001684714); Allpt7TeVScaling->SetBinError(31,0.001755251); Allpt7TeVScaling->SetBinError(32,0.001658064); Allpt7TeVScaling->SetBinError(33,0.001734642); Allpt7TeVScaling->SetBinError(34,0.001624169); Allpt7TeVScaling->SetBinError(35,0.001893747); Allpt7TeVScaling->SetBinError(36,0.001748243); Allpt7TeVScaling->SetBinError(37,0.001824646); Allpt7TeVScaling->SetBinError(38,0.001587514); Allpt7TeVScaling->SetBinError(39,0.001697203); Allpt7TeVScaling->SetBinError(40,0.002112545); Allpt7TeVScaling->SetBinError(41,0.002060819); Allpt7TeVScaling->SetBinError(42,0.002026571); Allpt7TeVScaling->SetBinError(43,0.001821759); Allpt7TeVScaling->SetBinError(44,0.002051793); Allpt7TeVScaling->SetBinError(45,0.002193957); Allpt7TeVScaling->SetBinError(46,0.001991362); Allpt7TeVScaling->SetBinError(47,0.002242771); Allpt7TeVScaling->SetBinError(48,0.002305595); Allpt7TeVScaling->SetBinError(49,0.002335076); Allpt7TeVScaling->SetBinError(50,0.002798601); Allpt7TeVScaling->SetBinError(51,0.003536788); Allpt7TeVScaling->SetBinError(52,0.003440952); Allpt7TeVScaling->SetBinError(53,0.003495971); Allpt7TeVScaling->SetBinError(54,0.003067046); Allpt7TeVScaling->SetBinError(55,0.003553133); Allpt7TeVScaling->SetBinError(56,0.003084681); Allpt7TeVScaling->SetBinError(57,0.003207651); Allpt7TeVScaling->SetBinError(58,0.002958898); Allpt7TeVScaling->SetBinError(59,0.003486212); Allpt7TeVScaling->SetBinError(60,0.003160784); Allpt7TeVScaling->SetBinError(61,0.00221753); Allpt7TeVScaling->SetBinError(62,0.002573063); Allpt7TeVScaling->SetBinError(63,0.003014171); Allpt7TeVScaling->SetBinError(64,0.003535776); Allpt7TeVScaling->SetBinError(65,0.003472028); Allpt7TeVScaling->SetBinError(66,0.003649103); Allpt7TeVScaling->SetBinError(67,0.004547596); Allpt7TeVScaling->SetBinError(68,0.004987629); Allpt7TeVScaling->SetBinError(69,0.005726305); Allpt7TeVScaling->SetBinError(70,0.003978767); Allpt7TeVScaling->SetBinError(71,0.00668463); Allpt7TeVScaling->SetBinError(72,0.005573138); Allpt7TeVScaling->SetBinError(73,0.005310436); Allpt7TeVScaling->SetBinError(74,0.007698246); Allpt7TeVScaling->SetBinError(75,0.01002912); Allpt7TeVScaling->SetBinError(76,0.006029285); Allpt7TeVScaling->SetBinError(77,0.0117225); Allpt7TeVScaling->SetBinError(78,0.01333887); Allpt7TeVScaling->SetBinError(79,0.009949805); Allpt7TeVScaling->SetBinError(80,0.01335105); Allpt7TeVScaling->SetBinError(81,0.004966615); Allpt7TeVScaling->SetBinError(82,0.006669427); Allpt7TeVScaling->SetBinError(83,0.007591192); Allpt7TeVScaling->SetBinError(84,0.01217267); Allpt7TeVScaling->SetBinError(85,0.002708538); Allpt7TeVScaling->SetBinError(86,0.008143581); Allpt7TeVScaling->SetBinError(87,0.01149673); Allpt7TeVScaling->SetBinError(88,0.01759547); Allpt7TeVScaling->SetBinError(89,0.01766708); Allpt7TeVScaling->SetBinError(90,0.01339557); Allpt7TeVScaling->SetBinError(91,0.004886052); Allpt7TeVScaling->SetBinError(93,0.01960311); Allpt7TeVScaling->SetBinError(94,0.008892003); Allpt7TeVScaling->SetBinError(95,0.01001297); Allpt7TeVScaling->SetBinError(96,0.008178781); Allpt7TeVScaling->SetBinError(97,0.007527505); Allpt7TeVScaling->SetBinError(98,0.01463862); Allpt7TeVScaling->SetBinError(103,0.03621621); Allpt7TeVScaling->SetMinimum(0); Allpt7TeVScaling->SetMaximum(0.04); Allpt7TeVScaling->SetEntries(1949.423); Allpt7TeVScaling->SetStats(0); Allpt7TeVScaling->SetMarkerStyle(20); Allpt7TeVScaling->GetXaxis()->SetTitle("p_{T}"); Allpt7TeVScaling->GetXaxis()->SetRange(1,91); Allpt7TeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); Allpt7TeVScaling->GetYaxis()->SetTitleOffset(1.2); TH1D *Allpt900GeVScaling = new TH1D("Allpt900GeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); Allpt900GeVScaling->SetBinContent(13,0.2098361); Allpt900GeVScaling->SetBinContent(14,0.1557776); Allpt900GeVScaling->SetBinContent(15,0.1480364); Allpt900GeVScaling->SetBinContent(16,0.1368096); Allpt900GeVScaling->SetBinContent(17,0.1244875); Allpt900GeVScaling->SetBinContent(18,0.1119436); Allpt900GeVScaling->SetBinContent(19,0.104192); Allpt900GeVScaling->SetBinContent(20,0.0938359); Allpt900GeVScaling->SetBinContent(21,0.08228509); Allpt900GeVScaling->SetBinContent(22,0.06783117); Allpt900GeVScaling->SetBinContent(23,0.05079209); Allpt900GeVScaling->SetBinContent(24,0.04795673); Allpt900GeVScaling->SetBinContent(25,0.04181242); Allpt900GeVScaling->SetBinContent(26,0.03994469); Allpt900GeVScaling->SetBinContent(27,0.03981588); Allpt900GeVScaling->SetBinContent(28,0.0377358); Allpt900GeVScaling->SetBinContent(29,0.0399667); Allpt900GeVScaling->SetBinContent(30,0.04095473); Allpt900GeVScaling->SetBinContent(31,0.04058456); Allpt900GeVScaling->SetBinContent(32,0.03641218); Allpt900GeVScaling->SetBinContent(33,0.03599966); Allpt900GeVScaling->SetBinContent(34,0.03853838); Allpt900GeVScaling->SetBinContent(35,0.04105409); Allpt900GeVScaling->SetBinContent(36,0.03461075); Allpt900GeVScaling->SetBinContent(37,0.03478729); Allpt900GeVScaling->SetBinContent(38,0.03150168); Allpt900GeVScaling->SetBinContent(39,0.03121329); Allpt900GeVScaling->SetBinContent(40,0.0370806); Allpt900GeVScaling->SetBinContent(41,0.03322445); Allpt900GeVScaling->SetBinContent(42,0.03200198); Allpt900GeVScaling->SetBinContent(43,0.02753292); Allpt900GeVScaling->SetBinContent(44,0.02776155); Allpt900GeVScaling->SetBinContent(45,0.02671865); Allpt900GeVScaling->SetBinContent(46,0.02612601); Allpt900GeVScaling->SetBinContent(47,0.02779079); Allpt900GeVScaling->SetBinContent(48,0.02762603); Allpt900GeVScaling->SetBinContent(49,0.0288977); Allpt900GeVScaling->SetBinContent(50,0.03181282); Allpt900GeVScaling->SetBinContent(51,0.04194707); Allpt900GeVScaling->SetBinContent(52,0.03287613); Allpt900GeVScaling->SetBinContent(53,0.04363778); Allpt900GeVScaling->SetBinContent(54,0.03003709); Allpt900GeVScaling->SetBinContent(55,0.03286254); Allpt900GeVScaling->SetBinContent(56,0.02585759); Allpt900GeVScaling->SetBinContent(57,0.03390592); Allpt900GeVScaling->SetBinContent(58,0.03212414); Allpt900GeVScaling->SetBinContent(59,0.03559123); Allpt900GeVScaling->SetBinContent(60,0.02464371); Allpt900GeVScaling->SetBinContent(61,0.03680997); Allpt900GeVScaling->SetBinContent(62,0.03569024); Allpt900GeVScaling->SetBinContent(63,0.03662885); Allpt900GeVScaling->SetBinContent(64,0.04690057); Allpt900GeVScaling->SetBinContent(65,0.04036927); Allpt900GeVScaling->SetBinContent(66,0.03237555); Allpt900GeVScaling->SetBinContent(67,0.04285642); Allpt900GeVScaling->SetBinContent(68,0.04778185); Allpt900GeVScaling->SetBinContent(69,0.0411135); Allpt900GeVScaling->SetBinContent(70,0.02465524); Allpt900GeVScaling->SetBinContent(71,0.03686565); Allpt900GeVScaling->SetBinContent(72,0.03694553); Allpt900GeVScaling->SetBinContent(73,0.02206743); Allpt900GeVScaling->SetBinContent(74,0.04854334); Allpt900GeVScaling->SetBinContent(75,0.04190287); Allpt900GeVScaling->SetBinContent(76,0.03255752); Allpt900GeVScaling->SetBinContent(77,0.0715847); Allpt900GeVScaling->SetBinContent(78,0.07329647); Allpt900GeVScaling->SetBinContent(79,0.02727783); Allpt900GeVScaling->SetBinContent(80,0.041527); Allpt900GeVScaling->SetBinContent(81,0.0400001); Allpt900GeVScaling->SetBinContent(82,0.0246785); Allpt900GeVScaling->SetBinContent(83,0.03546198); Allpt900GeVScaling->SetBinContent(84,0.0244789); Allpt900GeVScaling->SetBinContent(85,0.01692187); Allpt900GeVScaling->SetBinContent(86,0.01669025); Allpt900GeVScaling->SetBinContent(87,0.01506698); Allpt900GeVScaling->SetBinContent(88,0.02896953); Allpt900GeVScaling->SetBinContent(89,0.04019321); Allpt900GeVScaling->SetBinContent(90,0.0108909); Allpt900GeVScaling->SetBinContent(91,0.02059646); Allpt900GeVScaling->SetBinContent(92,0.01781294); Allpt900GeVScaling->SetBinContent(93,0.01930956); Allpt900GeVScaling->SetBinContent(94,0.0259471); Allpt900GeVScaling->SetBinContent(95,0.0114321); Allpt900GeVScaling->SetBinContent(98,0.04961994); Allpt900GeVScaling->SetBinContent(103,0.1436233); Allpt900GeVScaling->SetBinError(13,0.04135227); Allpt900GeVScaling->SetBinError(14,0.007955716); Allpt900GeVScaling->SetBinError(15,0.005353158); Allpt900GeVScaling->SetBinError(16,0.00465994); Allpt900GeVScaling->SetBinError(17,0.004317818); Allpt900GeVScaling->SetBinError(18,0.003992759); Allpt900GeVScaling->SetBinError(19,0.003794229); Allpt900GeVScaling->SetBinError(20,0.003428406); Allpt900GeVScaling->SetBinError(21,0.002309662); Allpt900GeVScaling->SetBinError(22,0.002082593); Allpt900GeVScaling->SetBinError(23,0.001768867); Allpt900GeVScaling->SetBinError(24,0.001734781); Allpt900GeVScaling->SetBinError(25,0.001668325); Allpt900GeVScaling->SetBinError(26,0.001736728); Allpt900GeVScaling->SetBinError(27,0.001793642); Allpt900GeVScaling->SetBinError(28,0.001824914); Allpt900GeVScaling->SetBinError(29,0.002216558); Allpt900GeVScaling->SetBinError(30,0.00224095); Allpt900GeVScaling->SetBinError(31,0.002244602); Allpt900GeVScaling->SetBinError(32,0.002073202); Allpt900GeVScaling->SetBinError(33,0.001913981); Allpt900GeVScaling->SetBinError(34,0.002456659); Allpt900GeVScaling->SetBinError(35,0.002596584); Allpt900GeVScaling->SetBinError(36,0.002232065); Allpt900GeVScaling->SetBinError(37,0.00251725); Allpt900GeVScaling->SetBinError(38,0.002196274); Allpt900GeVScaling->SetBinError(39,0.002513529); Allpt900GeVScaling->SetBinError(40,0.002973807); Allpt900GeVScaling->SetBinError(41,0.002807704); Allpt900GeVScaling->SetBinError(42,0.002773487); Allpt900GeVScaling->SetBinError(43,0.002527062); Allpt900GeVScaling->SetBinError(44,0.002921618); Allpt900GeVScaling->SetBinError(45,0.002907984); Allpt900GeVScaling->SetBinError(46,0.003115439); Allpt900GeVScaling->SetBinError(47,0.003224686); Allpt900GeVScaling->SetBinError(48,0.00312276); Allpt900GeVScaling->SetBinError(49,0.00393121); Allpt900GeVScaling->SetBinError(50,0.003948797); Allpt900GeVScaling->SetBinError(51,0.005565144); Allpt900GeVScaling->SetBinError(52,0.004615559); Allpt900GeVScaling->SetBinError(53,0.006431509); Allpt900GeVScaling->SetBinError(54,0.004787345); Allpt900GeVScaling->SetBinError(55,0.004506561); Allpt900GeVScaling->SetBinError(56,0.004277351); Allpt900GeVScaling->SetBinError(57,0.005084143); Allpt900GeVScaling->SetBinError(58,0.005085463); Allpt900GeVScaling->SetBinError(59,0.005312616); Allpt900GeVScaling->SetBinError(60,0.003899113); Allpt900GeVScaling->SetBinError(61,0.003977416); Allpt900GeVScaling->SetBinError(62,0.004292332); Allpt900GeVScaling->SetBinError(63,0.005072534); Allpt900GeVScaling->SetBinError(64,0.006271891); Allpt900GeVScaling->SetBinError(65,0.006463988); Allpt900GeVScaling->SetBinError(66,0.005869634); Allpt900GeVScaling->SetBinError(67,0.008405251); Allpt900GeVScaling->SetBinError(68,0.009847454); Allpt900GeVScaling->SetBinError(69,0.009749354); Allpt900GeVScaling->SetBinError(70,0.007179523); Allpt900GeVScaling->SetBinError(71,0.01004867); Allpt900GeVScaling->SetBinError(72,0.01065873); Allpt900GeVScaling->SetBinError(73,0.007849093); Allpt900GeVScaling->SetBinError(74,0.01402708); Allpt900GeVScaling->SetBinError(75,0.01330162); Allpt900GeVScaling->SetBinError(76,0.01245402); Allpt900GeVScaling->SetBinError(77,0.02130282); Allpt900GeVScaling->SetBinError(78,0.02296044); Allpt900GeVScaling->SetBinError(79,0.01175655); Allpt900GeVScaling->SetBinError(80,0.0187519); Allpt900GeVScaling->SetBinError(81,0.009932965); Allpt900GeVScaling->SetBinError(82,0.007798875); Allpt900GeVScaling->SetBinError(83,0.01126926); Allpt900GeVScaling->SetBinError(84,0.01123177); Allpt900GeVScaling->SetBinError(85,0.01013651); Allpt900GeVScaling->SetBinError(86,0.01187595); Allpt900GeVScaling->SetBinError(87,0.009213867); Allpt900GeVScaling->SetBinError(88,0.02054123); Allpt900GeVScaling->SetBinError(89,0.01976358); Allpt900GeVScaling->SetBinError(90,0.006322856); Allpt900GeVScaling->SetBinError(91,0.01400211); Allpt900GeVScaling->SetBinError(92,0.01786717); Allpt900GeVScaling->SetBinError(93,0.01938449); Allpt900GeVScaling->SetBinError(94,0.01853139); Allpt900GeVScaling->SetBinError(95,0.01149767); Allpt900GeVScaling->SetBinError(98,0.03786105); Allpt900GeVScaling->SetBinError(103,0.1493976); Allpt900GeVScaling->SetMinimum(0); Allpt900GeVScaling->SetMaximum(0.04); Allpt900GeVScaling->SetEntries(482.6447); Allpt900GeVScaling->SetStats(0); Allpt900GeVScaling->SetMarkerStyle(20); Allpt900GeVScaling->GetXaxis()->SetTitle("p_{T}"); Allpt900GeVScaling->GetXaxis()->SetRange(1,91); Allpt900GeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); Allpt900GeVScaling->GetYaxis()->SetTitleOffset(1.2); Allpt7TeVScaling->Add(Allpt900GeVScaling); Allpt7TeVScaling->Scale(0.5); delete Allpt900GeVScaling; return Allpt7TeVScaling; } TH1D *pp276ITSBkgd(){ TH1D *Allpt7TeVScaling = new TH1D("Allpt7TeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); Allpt7TeVScaling->SetBinContent(11,0.2557976); Allpt7TeVScaling->SetBinContent(12,0.1130502); Allpt7TeVScaling->SetBinContent(13,0.07504393); Allpt7TeVScaling->SetBinContent(14,0.05774832); Allpt7TeVScaling->SetBinContent(15,0.06389806); Allpt7TeVScaling->SetBinContent(16,0.0468956); Allpt7TeVScaling->SetBinContent(17,0.04355167); Allpt7TeVScaling->SetBinContent(18,0.04035013); Allpt7TeVScaling->SetBinContent(19,0.04359472); Allpt7TeVScaling->SetBinContent(20,0.03392261); Allpt7TeVScaling->SetBinContent(21,0.0327676); Allpt7TeVScaling->SetBinContent(22,0.02916985); Allpt7TeVScaling->SetBinContent(23,0.02603138); Allpt7TeVScaling->SetBinContent(24,0.02298032); Allpt7TeVScaling->SetBinContent(25,0.02366939); Allpt7TeVScaling->SetBinContent(26,0.01941947); Allpt7TeVScaling->SetBinContent(27,0.02077003); Allpt7TeVScaling->SetBinContent(28,0.01985847); Allpt7TeVScaling->SetBinContent(29,0.0177459); Allpt7TeVScaling->SetBinContent(30,0.01889314); Allpt7TeVScaling->SetBinContent(31,0.01783637); Allpt7TeVScaling->SetBinContent(32,0.01742894); Allpt7TeVScaling->SetBinContent(33,0.01400672); Allpt7TeVScaling->SetBinContent(34,0.01509014); Allpt7TeVScaling->SetBinContent(35,0.01473409); Allpt7TeVScaling->SetBinContent(36,0.01337399); Allpt7TeVScaling->SetBinContent(37,0.01458302); Allpt7TeVScaling->SetBinContent(38,0.01103552); Allpt7TeVScaling->SetBinContent(39,0.01440217); Allpt7TeVScaling->SetBinContent(40,0.0147644); Allpt7TeVScaling->SetBinContent(41,0.01459072); Allpt7TeVScaling->SetBinContent(42,0.01318416); Allpt7TeVScaling->SetBinContent(43,0.01364365); Allpt7TeVScaling->SetBinContent(44,0.01186023); Allpt7TeVScaling->SetBinContent(45,0.01218222); Allpt7TeVScaling->SetBinContent(46,0.01015807); Allpt7TeVScaling->SetBinContent(47,0.01158454); Allpt7TeVScaling->SetBinContent(48,0.009338749); Allpt7TeVScaling->SetBinContent(49,0.01320016); Allpt7TeVScaling->SetBinContent(50,0.01029357); Allpt7TeVScaling->SetBinContent(51,0.01238829); Allpt7TeVScaling->SetBinContent(52,0.008714745); Allpt7TeVScaling->SetBinContent(53,0.01300408); Allpt7TeVScaling->SetBinContent(54,0.009445714); Allpt7TeVScaling->SetBinContent(55,0.01301231); Allpt7TeVScaling->SetBinContent(56,0.01010809); Allpt7TeVScaling->SetBinContent(57,0.01181843); Allpt7TeVScaling->SetBinContent(58,0.01037317); Allpt7TeVScaling->SetBinContent(59,0.01172911); Allpt7TeVScaling->SetBinContent(60,0.007576451); Allpt7TeVScaling->SetBinContent(61,0.011185); Allpt7TeVScaling->SetBinContent(62,0.01009882); Allpt7TeVScaling->SetBinContent(63,0.01042078); Allpt7TeVScaling->SetBinContent(64,0.01079128); Allpt7TeVScaling->SetBinContent(65,0.008657009); Allpt7TeVScaling->SetBinContent(66,0.009641373); Allpt7TeVScaling->SetBinContent(67,0.01201938); Allpt7TeVScaling->SetBinContent(68,0.01031088); Allpt7TeVScaling->SetBinContent(69,0.009905995); Allpt7TeVScaling->SetBinContent(70,0.009569216); Allpt7TeVScaling->SetBinContent(71,0.009473246); Allpt7TeVScaling->SetBinContent(72,0.00825013); Allpt7TeVScaling->SetBinContent(73,0.008405926); Allpt7TeVScaling->SetBinContent(74,0.006479467); Allpt7TeVScaling->SetBinContent(75,0.01135131); Allpt7TeVScaling->SetBinContent(76,0.007318432); Allpt7TeVScaling->SetBinContent(77,0.01145386); Allpt7TeVScaling->SetBinContent(78,0.01425465); Allpt7TeVScaling->SetBinContent(79,0.01205772); Allpt7TeVScaling->SetBinContent(80,0.009454997); Allpt7TeVScaling->SetBinContent(81,0.006863416); Allpt7TeVScaling->SetBinContent(82,0.006285407); Allpt7TeVScaling->SetBinContent(83,0.004974776); Allpt7TeVScaling->SetBinContent(84,0.00748478); Allpt7TeVScaling->SetBinContent(85,0.002485323); Allpt7TeVScaling->SetBinContent(86,0.009973802); Allpt7TeVScaling->SetBinContent(87,0.00543998); Allpt7TeVScaling->SetBinContent(89,0.01703395); Allpt7TeVScaling->SetBinContent(90,0.004394557); Allpt7TeVScaling->SetBinContent(96,0.01833441); Allpt7TeVScaling->SetBinContent(98,0.02964621); Allpt7TeVScaling->SetBinContent(101,0.0574538); Allpt7TeVScaling->SetBinError(11,0.07407797); Allpt7TeVScaling->SetBinError(12,0.01275005); Allpt7TeVScaling->SetBinError(13,0.00608927); Allpt7TeVScaling->SetBinError(14,0.00428747); Allpt7TeVScaling->SetBinError(15,0.004071853); Allpt7TeVScaling->SetBinError(16,0.003225169); Allpt7TeVScaling->SetBinError(17,0.002968269); Allpt7TeVScaling->SetBinError(18,0.002761875); Allpt7TeVScaling->SetBinError(19,0.002800096); Allpt7TeVScaling->SetBinError(20,0.002386005); Allpt7TeVScaling->SetBinError(21,0.001589813); Allpt7TeVScaling->SetBinError(22,0.001475329); Allpt7TeVScaling->SetBinError(23,0.001364666); Allpt7TeVScaling->SetBinError(24,0.00126007); Allpt7TeVScaling->SetBinError(25,0.001265789); Allpt7TeVScaling->SetBinError(26,0.001156208); Allpt7TeVScaling->SetBinError(27,0.001182566); Allpt7TeVScaling->SetBinError(28,0.001171098); Allpt7TeVScaling->SetBinError(29,0.001122203); Allpt7TeVScaling->SetBinError(30,0.001173633); Allpt7TeVScaling->SetBinError(31,0.001154546); Allpt7TeVScaling->SetBinError(32,0.001162445); Allpt7TeVScaling->SetBinError(33,0.001067574); Allpt7TeVScaling->SetBinError(34,0.001116243); Allpt7TeVScaling->SetBinError(35,0.001144851); Allpt7TeVScaling->SetBinError(36,0.001110251); Allpt7TeVScaling->SetBinError(37,0.001192711); Allpt7TeVScaling->SetBinError(38,0.001062774); Allpt7TeVScaling->SetBinError(39,0.001255752); Allpt7TeVScaling->SetBinError(40,0.001319021); Allpt7TeVScaling->SetBinError(41,0.001323284); Allpt7TeVScaling->SetBinError(42,0.001303463); Allpt7TeVScaling->SetBinError(43,0.001388458); Allpt7TeVScaling->SetBinError(44,0.001312344); Allpt7TeVScaling->SetBinError(45,0.001383885); Allpt7TeVScaling->SetBinError(46,0.001289152); Allpt7TeVScaling->SetBinError(47,0.001430287); Allpt7TeVScaling->SetBinError(48,0.001302716); Allpt7TeVScaling->SetBinError(49,0.001578988); Allpt7TeVScaling->SetBinError(50,0.0014581); Allpt7TeVScaling->SetBinError(51,0.00165325); Allpt7TeVScaling->SetBinError(52,0.001404793); Allpt7TeVScaling->SetBinError(53,0.001757303); Allpt7TeVScaling->SetBinError(54,0.001562987); Allpt7TeVScaling->SetBinError(55,0.001911932); Allpt7TeVScaling->SetBinError(56,0.001695077); Allpt7TeVScaling->SetBinError(57,0.001932989); Allpt7TeVScaling->SetBinError(58,0.001844963); Allpt7TeVScaling->SetBinError(59,0.002000117); Allpt7TeVScaling->SetBinError(60,0.001689391); Allpt7TeVScaling->SetBinError(61,0.001386339); Allpt7TeVScaling->SetBinError(62,0.001349763); Allpt7TeVScaling->SetBinError(63,0.001454327); Allpt7TeVScaling->SetBinError(64,0.001584791); Allpt7TeVScaling->SetBinError(65,0.00151686); Allpt7TeVScaling->SetBinError(66,0.001668507); Allpt7TeVScaling->SetBinError(67,0.001963962); Allpt7TeVScaling->SetBinError(68,0.001970941); Allpt7TeVScaling->SetBinError(69,0.001992913); Allpt7TeVScaling->SetBinError(70,0.002006771); Allpt7TeVScaling->SetBinError(71,0.002187254); Allpt7TeVScaling->SetBinError(72,0.002140473); Allpt7TeVScaling->SetBinError(73,0.002256925); Allpt7TeVScaling->SetBinError(74,0.002056219); Allpt7TeVScaling->SetBinError(75,0.002854892); Allpt7TeVScaling->SetBinError(76,0.002450159); Allpt7TeVScaling->SetBinError(77,0.003197387); Allpt7TeVScaling->SetBinError(78,0.003840947); Allpt7TeVScaling->SetBinError(79,0.003660464); Allpt7TeVScaling->SetBinError(80,0.003363956); Allpt7TeVScaling->SetBinError(81,0.001624793); Allpt7TeVScaling->SetBinError(82,0.00181065); Allpt7TeVScaling->SetBinError(83,0.001888966); Allpt7TeVScaling->SetBinError(84,0.002671173); Allpt7TeVScaling->SetBinError(85,0.001759843); Allpt7TeVScaling->SetBinError(86,0.004094911); Allpt7TeVScaling->SetBinError(87,0.003871198); Allpt7TeVScaling->SetBinError(89,0.008630108); Allpt7TeVScaling->SetBinError(90,0.004405367); Allpt7TeVScaling->SetBinError(96,0.01308755); Allpt7TeVScaling->SetBinError(98,0.02127892); Allpt7TeVScaling->SetBinError(101,0.05904744); Allpt7TeVScaling->SetMinimum(0); Allpt7TeVScaling->SetMaximum(0.04); Allpt7TeVScaling->SetEntries(300.6659); Allpt7TeVScaling->SetStats(0); Allpt7TeVScaling->SetMarkerStyle(20); Allpt7TeVScaling->GetXaxis()->SetTitle("p_{T}"); Allpt7TeVScaling->GetXaxis()->SetRange(1,61); Allpt7TeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); Allpt7TeVScaling->GetYaxis()->SetTitleOffset(1.2); TH1D *Allpt900GeVScaling = new TH1D("Allpt900GeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); Allpt900GeVScaling->SetBinContent(9,1); Allpt900GeVScaling->SetBinContent(11,0.2834258); Allpt900GeVScaling->SetBinContent(12,0.1188776); Allpt900GeVScaling->SetBinContent(13,0.08085647); Allpt900GeVScaling->SetBinContent(14,0.05677858); Allpt900GeVScaling->SetBinContent(15,0.0616757); Allpt900GeVScaling->SetBinContent(16,0.04844882); Allpt900GeVScaling->SetBinContent(17,0.03999142); Allpt900GeVScaling->SetBinContent(18,0.03890007); Allpt900GeVScaling->SetBinContent(19,0.0406008); Allpt900GeVScaling->SetBinContent(20,0.02900627); Allpt900GeVScaling->SetBinContent(21,0.03088961); Allpt900GeVScaling->SetBinContent(22,0.03033992); Allpt900GeVScaling->SetBinContent(23,0.02521969); Allpt900GeVScaling->SetBinContent(24,0.02143008); Allpt900GeVScaling->SetBinContent(25,0.02266187); Allpt900GeVScaling->SetBinContent(26,0.01981764); Allpt900GeVScaling->SetBinContent(27,0.02127366); Allpt900GeVScaling->SetBinContent(28,0.0192753); Allpt900GeVScaling->SetBinContent(29,0.01854744); Allpt900GeVScaling->SetBinContent(30,0.01903471); Allpt900GeVScaling->SetBinContent(31,0.01819287); Allpt900GeVScaling->SetBinContent(32,0.01715784); Allpt900GeVScaling->SetBinContent(33,0.01488607); Allpt900GeVScaling->SetBinContent(34,0.0156348); Allpt900GeVScaling->SetBinContent(35,0.01699235); Allpt900GeVScaling->SetBinContent(36,0.01398954); Allpt900GeVScaling->SetBinContent(37,0.01435865); Allpt900GeVScaling->SetBinContent(38,0.01236315); Allpt900GeVScaling->SetBinContent(39,0.01282288); Allpt900GeVScaling->SetBinContent(40,0.01668125); Allpt900GeVScaling->SetBinContent(41,0.01372736); Allpt900GeVScaling->SetBinContent(42,0.01433745); Allpt900GeVScaling->SetBinContent(43,0.01477496); Allpt900GeVScaling->SetBinContent(44,0.0116776); Allpt900GeVScaling->SetBinContent(45,0.01245205); Allpt900GeVScaling->SetBinContent(46,0.009475905); Allpt900GeVScaling->SetBinContent(47,0.01117031); Allpt900GeVScaling->SetBinContent(48,0.009558273); Allpt900GeVScaling->SetBinContent(49,0.01404047); Allpt900GeVScaling->SetBinContent(50,0.01037464); Allpt900GeVScaling->SetBinContent(51,0.01195013); Allpt900GeVScaling->SetBinContent(52,0.01155682); Allpt900GeVScaling->SetBinContent(53,0.0137127); Allpt900GeVScaling->SetBinContent(54,0.009255665); Allpt900GeVScaling->SetBinContent(55,0.01133598); Allpt900GeVScaling->SetBinContent(56,0.009080946); Allpt900GeVScaling->SetBinContent(57,0.01124535); Allpt900GeVScaling->SetBinContent(58,0.008398657); Allpt900GeVScaling->SetBinContent(59,0.008324869); Allpt900GeVScaling->SetBinContent(60,0.007103489); Allpt900GeVScaling->SetBinContent(61,0.01142644); Allpt900GeVScaling->SetBinContent(62,0.0110906); Allpt900GeVScaling->SetBinContent(63,0.01076479); Allpt900GeVScaling->SetBinContent(64,0.01168746); Allpt900GeVScaling->SetBinContent(65,0.008513672); Allpt900GeVScaling->SetBinContent(66,0.00828279); Allpt900GeVScaling->SetBinContent(67,0.01140648); Allpt900GeVScaling->SetBinContent(68,0.01111323); Allpt900GeVScaling->SetBinContent(69,0.007425678); Allpt900GeVScaling->SetBinContent(70,0.009847098); Allpt900GeVScaling->SetBinContent(71,0.009834419); Allpt900GeVScaling->SetBinContent(72,0.008930991); Allpt900GeVScaling->SetBinContent(73,0.00669819); Allpt900GeVScaling->SetBinContent(74,0.00662618); Allpt900GeVScaling->SetBinContent(75,0.01370167); Allpt900GeVScaling->SetBinContent(76,0.006451024); Allpt900GeVScaling->SetBinContent(77,0.01307605); Allpt900GeVScaling->SetBinContent(78,0.01579264); Allpt900GeVScaling->SetBinContent(79,0.01255205); Allpt900GeVScaling->SetBinContent(80,0.01010045); Allpt900GeVScaling->SetBinContent(81,0.01082155); Allpt900GeVScaling->SetBinContent(82,0.006495255); Allpt900GeVScaling->SetBinContent(83,0.005512525); Allpt900GeVScaling->SetBinContent(84,0.003527945); Allpt900GeVScaling->SetBinContent(85,0.001806946); Allpt900GeVScaling->SetBinContent(86,0.004225982); Allpt900GeVScaling->SetBinContent(89,0.01106803); Allpt900GeVScaling->SetBinContent(90,0.005044269); Allpt900GeVScaling->SetBinContent(98,0.04456838); Allpt900GeVScaling->SetBinContent(101,0.07960857); Allpt900GeVScaling->SetBinError(9,1.414214); Allpt900GeVScaling->SetBinError(11,0.08936516); Allpt900GeVScaling->SetBinError(12,0.01474769); Allpt900GeVScaling->SetBinError(13,0.006947232); Allpt900GeVScaling->SetBinError(14,0.00470276); Allpt900GeVScaling->SetBinError(15,0.004394222); Allpt900GeVScaling->SetBinError(16,0.003598075); Allpt900GeVScaling->SetBinError(17,0.00312581); Allpt900GeVScaling->SetBinError(18,0.003002622); Allpt900GeVScaling->SetBinError(19,0.002976063); Allpt900GeVScaling->SetBinError(20,0.002438669); Allpt900GeVScaling->SetBinError(21,0.00172144); Allpt900GeVScaling->SetBinError(22,0.001661658); Allpt900GeVScaling->SetBinError(23,0.001477901); Allpt900GeVScaling->SetBinError(24,0.001329128); Allpt900GeVScaling->SetBinError(25,0.00137131); Allpt900GeVScaling->SetBinError(26,0.001299862); Allpt900GeVScaling->SetBinError(27,0.001332285); Allpt900GeVScaling->SetBinError(28,0.001274491); Allpt900GeVScaling->SetBinError(29,0.001266741); Allpt900GeVScaling->SetBinError(30,0.001305369); Allpt900GeVScaling->SetBinError(31,0.001295113); Allpt900GeVScaling->SetBinError(32,0.001278961); Allpt900GeVScaling->SetBinError(33,0.00122243); Allpt900GeVScaling->SetBinError(34,0.001256858); Allpt900GeVScaling->SetBinError(35,0.001357681); Allpt900GeVScaling->SetBinError(36,0.001267233); Allpt900GeVScaling->SetBinError(37,0.00131839); Allpt900GeVScaling->SetBinError(38,0.001256835); Allpt900GeVScaling->SetBinError(39,0.001306764); Allpt900GeVScaling->SetBinError(40,0.001544152); Allpt900GeVScaling->SetBinError(41,0.001438378); Allpt900GeVScaling->SetBinError(42,0.001507702); Allpt900GeVScaling->SetBinError(43,0.001615741); Allpt900GeVScaling->SetBinError(44,0.00143817); Allpt900GeVScaling->SetBinError(45,0.001536765); Allpt900GeVScaling->SetBinError(46,0.001377493); Allpt900GeVScaling->SetBinError(47,0.001536932); Allpt900GeVScaling->SetBinError(48,0.001466076); Allpt900GeVScaling->SetBinError(49,0.001814185); Allpt900GeVScaling->SetBinError(50,0.001620601); Allpt900GeVScaling->SetBinError(51,0.001814658); Allpt900GeVScaling->SetBinError(52,0.001822325); Allpt900GeVScaling->SetBinError(53,0.001976947); Allpt900GeVScaling->SetBinError(54,0.001729338); Allpt900GeVScaling->SetBinError(55,0.001988634); Allpt900GeVScaling->SetBinError(56,0.001792157); Allpt900GeVScaling->SetBinError(57,0.00206644); Allpt900GeVScaling->SetBinError(58,0.001842037); Allpt900GeVScaling->SetBinError(59,0.001826512); Allpt900GeVScaling->SetBinError(60,0.001732292); Allpt900GeVScaling->SetBinError(61,0.001512696); Allpt900GeVScaling->SetBinError(62,0.001543542); Allpt900GeVScaling->SetBinError(63,0.001633358); Allpt900GeVScaling->SetBinError(64,0.001816012); Allpt900GeVScaling->SetBinError(65,0.001680175); Allpt900GeVScaling->SetBinError(66,0.001692314); Allpt900GeVScaling->SetBinError(67,0.00213108); Allpt900GeVScaling->SetBinError(68,0.002295308); Allpt900GeVScaling->SetBinError(69,0.001926907); Allpt900GeVScaling->SetBinError(70,0.002272706); Allpt900GeVScaling->SetBinError(71,0.002475109); Allpt900GeVScaling->SetBinError(72,0.002399123); Allpt900GeVScaling->SetBinError(73,0.002241209); Allpt900GeVScaling->SetBinError(74,0.002365641); Allpt900GeVScaling->SetBinError(75,0.003450507); Allpt900GeVScaling->SetBinError(76,0.002503392); Allpt900GeVScaling->SetBinError(77,0.003801781); Allpt900GeVScaling->SetBinError(78,0.004417338); Allpt900GeVScaling->SetBinError(79,0.00421245); Allpt900GeVScaling->SetBinError(80,0.003844397); Allpt900GeVScaling->SetBinError(81,0.002274688); Allpt900GeVScaling->SetBinError(82,0.002043975); Allpt900GeVScaling->SetBinError(83,0.002262563); Allpt900GeVScaling->SetBinError(84,0.00204239); Allpt900GeVScaling->SetBinError(85,0.00180851); Allpt900GeVScaling->SetBinError(86,0.002994568); Allpt900GeVScaling->SetBinError(89,0.007939965); Allpt900GeVScaling->SetBinError(90,0.005058548); Allpt900GeVScaling->SetBinError(98,0.03222111); Allpt900GeVScaling->SetBinError(101,0.08264771); Allpt900GeVScaling->SetMinimum(0); Allpt900GeVScaling->SetMaximum(0.04); Allpt900GeVScaling->SetEntries(3.837105); Allpt900GeVScaling->SetStats(0); Allpt900GeVScaling->SetMarkerStyle(20); Allpt900GeVScaling->GetXaxis()->SetTitle("p_{T}"); Allpt900GeVScaling->GetXaxis()->SetRange(1,61); Allpt900GeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); Allpt900GeVScaling->GetYaxis()->SetTitleOffset(1.2); Allpt7TeVScaling->Add(Allpt900GeVScaling); Allpt7TeVScaling->Scale(0.5); delete Allpt900GeVScaling; return Allpt7TeVScaling; }
GetCorrections.C:1
GetCorrections.C:2
GetCorrections.C:3
GetCorrections.C:4
GetCorrections.C:5
GetCorrections.C:6
GetCorrections.C:7
GetCorrections.C:8
GetCorrections.C:9
GetCorrections.C:10
GetCorrections.C:11
GetCorrections.C:12
GetCorrections.C:13
GetCorrections.C:14
GetCorrections.C:15
GetCorrections.C:16
GetCorrections.C:17
GetCorrections.C:18
GetCorrections.C:19
GetCorrections.C:20
GetCorrections.C:21
GetCorrections.C:22
GetCorrections.C:23
GetCorrections.C:24
GetCorrections.C:25
GetCorrections.C:26
GetCorrections.C:27
GetCorrections.C:28
GetCorrections.C:29
GetCorrections.C:30
GetCorrections.C:31
GetCorrections.C:32
GetCorrections.C:33
GetCorrections.C:34
GetCorrections.C:35
GetCorrections.C:36
GetCorrections.C:37
GetCorrections.C:38
GetCorrections.C:39
GetCorrections.C:40
GetCorrections.C:41
GetCorrections.C:42
GetCorrections.C:43
GetCorrections.C:44
GetCorrections.C:45
GetCorrections.C:46
GetCorrections.C:47
GetCorrections.C:48
GetCorrections.C:49
GetCorrections.C:50
GetCorrections.C:51
GetCorrections.C:52
GetCorrections.C:53
GetCorrections.C:54
GetCorrections.C:55
GetCorrections.C:56
GetCorrections.C:57
GetCorrections.C:58
GetCorrections.C:59
GetCorrections.C:60
GetCorrections.C:61
GetCorrections.C:62
GetCorrections.C:63
GetCorrections.C:64
GetCorrections.C:65
GetCorrections.C:66
GetCorrections.C:67
GetCorrections.C:68
GetCorrections.C:69
GetCorrections.C:70
GetCorrections.C:71
GetCorrections.C:72
GetCorrections.C:73
GetCorrections.C:74
GetCorrections.C:75
GetCorrections.C:76
GetCorrections.C:77
GetCorrections.C:78
GetCorrections.C:79
GetCorrections.C:80
GetCorrections.C:81
GetCorrections.C:82
GetCorrections.C:83
GetCorrections.C:84
GetCorrections.C:85
GetCorrections.C:86
GetCorrections.C:87
GetCorrections.C:88
GetCorrections.C:89
GetCorrections.C:90
GetCorrections.C:91
GetCorrections.C:92
GetCorrections.C:93
GetCorrections.C:94
GetCorrections.C:95
GetCorrections.C:96
GetCorrections.C:97
GetCorrections.C:98
GetCorrections.C:99
GetCorrections.C:100
GetCorrections.C:101
GetCorrections.C:102
GetCorrections.C:103
GetCorrections.C:104
GetCorrections.C:105
GetCorrections.C:106
GetCorrections.C:107
GetCorrections.C:108
GetCorrections.C:109
GetCorrections.C:110
GetCorrections.C:111
GetCorrections.C:112
GetCorrections.C:113
GetCorrections.C:114
GetCorrections.C:115
GetCorrections.C:116
GetCorrections.C:117
GetCorrections.C:118
GetCorrections.C:119
GetCorrections.C:120
GetCorrections.C:121
GetCorrections.C:122
GetCorrections.C:123
GetCorrections.C:124
GetCorrections.C:125
GetCorrections.C:126
GetCorrections.C:127
GetCorrections.C:128
GetCorrections.C:129
GetCorrections.C:130
GetCorrections.C:131
GetCorrections.C:132
GetCorrections.C:133
GetCorrections.C:134
GetCorrections.C:135
GetCorrections.C:136
GetCorrections.C:137
GetCorrections.C:138
GetCorrections.C:139
GetCorrections.C:140
GetCorrections.C:141
GetCorrections.C:142
GetCorrections.C:143
GetCorrections.C:144
GetCorrections.C:145
GetCorrections.C:146
GetCorrections.C:147
GetCorrections.C:148
GetCorrections.C:149
GetCorrections.C:150
GetCorrections.C:151
GetCorrections.C:152
GetCorrections.C:153
GetCorrections.C:154
GetCorrections.C:155
GetCorrections.C:156
GetCorrections.C:157
GetCorrections.C:158
GetCorrections.C:159
GetCorrections.C:160
GetCorrections.C:161
GetCorrections.C:162
GetCorrections.C:163
GetCorrections.C:164
GetCorrections.C:165
GetCorrections.C:166
GetCorrections.C:167
GetCorrections.C:168
GetCorrections.C:169
GetCorrections.C:170
GetCorrections.C:171
GetCorrections.C:172
GetCorrections.C:173
GetCorrections.C:174
GetCorrections.C:175
GetCorrections.C:176
GetCorrections.C:177
GetCorrections.C:178
GetCorrections.C:179
GetCorrections.C:180
GetCorrections.C:181
GetCorrections.C:182
GetCorrections.C:183
GetCorrections.C:184
GetCorrections.C:185
GetCorrections.C:186
GetCorrections.C:187
GetCorrections.C:188
GetCorrections.C:189
GetCorrections.C:190
GetCorrections.C:191
GetCorrections.C:192
GetCorrections.C:193
GetCorrections.C:194
GetCorrections.C:195
GetCorrections.C:196
GetCorrections.C:197
GetCorrections.C:198
GetCorrections.C:199
GetCorrections.C:200
GetCorrections.C:201
GetCorrections.C:202
GetCorrections.C:203
GetCorrections.C:204
GetCorrections.C:205
GetCorrections.C:206
GetCorrections.C:207
GetCorrections.C:208
GetCorrections.C:209
GetCorrections.C:210
GetCorrections.C:211
GetCorrections.C:212
GetCorrections.C:213
GetCorrections.C:214
GetCorrections.C:215
GetCorrections.C:216
GetCorrections.C:217
GetCorrections.C:218
GetCorrections.C:219
GetCorrections.C:220
GetCorrections.C:221
GetCorrections.C:222
GetCorrections.C:223
GetCorrections.C:224
GetCorrections.C:225
GetCorrections.C:226
GetCorrections.C:227
GetCorrections.C:228
GetCorrections.C:229
GetCorrections.C:230
GetCorrections.C:231
GetCorrections.C:232
GetCorrections.C:233
GetCorrections.C:234
GetCorrections.C:235
GetCorrections.C:236
GetCorrections.C:237
GetCorrections.C:238
GetCorrections.C:239
GetCorrections.C:240
GetCorrections.C:241
GetCorrections.C:242
GetCorrections.C:243
GetCorrections.C:244
GetCorrections.C:245
GetCorrections.C:246
GetCorrections.C:247
GetCorrections.C:248
GetCorrections.C:249
GetCorrections.C:250
GetCorrections.C:251
GetCorrections.C:252
GetCorrections.C:253
GetCorrections.C:254
GetCorrections.C:255
GetCorrections.C:256
GetCorrections.C:257
GetCorrections.C:258
GetCorrections.C:259
GetCorrections.C:260
GetCorrections.C:261
GetCorrections.C:262
GetCorrections.C:263
GetCorrections.C:264
GetCorrections.C:265
GetCorrections.C:266
GetCorrections.C:267
GetCorrections.C:268
GetCorrections.C:269
GetCorrections.C:270
GetCorrections.C:271
GetCorrections.C:272
GetCorrections.C:273
GetCorrections.C:274
GetCorrections.C:275
GetCorrections.C:276
GetCorrections.C:277
GetCorrections.C:278
GetCorrections.C:279
GetCorrections.C:280
GetCorrections.C:281
GetCorrections.C:282
GetCorrections.C:283
GetCorrections.C:284
GetCorrections.C:285
GetCorrections.C:286
GetCorrections.C:287
GetCorrections.C:288
GetCorrections.C:289
GetCorrections.C:290
GetCorrections.C:291
GetCorrections.C:292
GetCorrections.C:293
GetCorrections.C:294
GetCorrections.C:295
GetCorrections.C:296
GetCorrections.C:297
GetCorrections.C:298
GetCorrections.C:299
GetCorrections.C:300
GetCorrections.C:301
GetCorrections.C:302
GetCorrections.C:303
GetCorrections.C:304
GetCorrections.C:305
GetCorrections.C:306
GetCorrections.C:307
GetCorrections.C:308
GetCorrections.C:309
GetCorrections.C:310
GetCorrections.C:311
GetCorrections.C:312
GetCorrections.C:313
GetCorrections.C:314
GetCorrections.C:315
GetCorrections.C:316
GetCorrections.C:317
GetCorrections.C:318
GetCorrections.C:319
GetCorrections.C:320
GetCorrections.C:321
GetCorrections.C:322
GetCorrections.C:323
GetCorrections.C:324
GetCorrections.C:325
GetCorrections.C:326
GetCorrections.C:327
GetCorrections.C:328
GetCorrections.C:329
GetCorrections.C:330
GetCorrections.C:331
GetCorrections.C:332
GetCorrections.C:333
GetCorrections.C:334
GetCorrections.C:335
GetCorrections.C:336
GetCorrections.C:337
GetCorrections.C:338
GetCorrections.C:339
GetCorrections.C:340
GetCorrections.C:341
GetCorrections.C:342
GetCorrections.C:343
GetCorrections.C:344
GetCorrections.C:345
GetCorrections.C:346
GetCorrections.C:347
GetCorrections.C:348
GetCorrections.C:349
GetCorrections.C:350
GetCorrections.C:351
GetCorrections.C:352
GetCorrections.C:353
GetCorrections.C:354
GetCorrections.C:355
GetCorrections.C:356
GetCorrections.C:357
GetCorrections.C:358
GetCorrections.C:359
GetCorrections.C:360
GetCorrections.C:361
GetCorrections.C:362
GetCorrections.C:363
GetCorrections.C:364
GetCorrections.C:365
GetCorrections.C:366
GetCorrections.C:367
GetCorrections.C:368
GetCorrections.C:369
GetCorrections.C:370
GetCorrections.C:371
GetCorrections.C:372
GetCorrections.C:373
GetCorrections.C:374
GetCorrections.C:375
GetCorrections.C:376
GetCorrections.C:377
GetCorrections.C:378
GetCorrections.C:379
GetCorrections.C:380
GetCorrections.C:381
GetCorrections.C:382
GetCorrections.C:383
GetCorrections.C:384
GetCorrections.C:385
GetCorrections.C:386
GetCorrections.C:387
GetCorrections.C:388
GetCorrections.C:389
GetCorrections.C:390
GetCorrections.C:391
GetCorrections.C:392
GetCorrections.C:393
GetCorrections.C:394
GetCorrections.C:395
GetCorrections.C:396
GetCorrections.C:397
GetCorrections.C:398
GetCorrections.C:399
GetCorrections.C:400
GetCorrections.C:401
GetCorrections.C:402
GetCorrections.C:403
GetCorrections.C:404
GetCorrections.C:405
GetCorrections.C:406
GetCorrections.C:407
GetCorrections.C:408
GetCorrections.C:409
GetCorrections.C:410
GetCorrections.C:411
GetCorrections.C:412
GetCorrections.C:413
GetCorrections.C:414
GetCorrections.C:415
GetCorrections.C:416
GetCorrections.C:417
GetCorrections.C:418
GetCorrections.C:419
GetCorrections.C:420
GetCorrections.C:421
GetCorrections.C:422
GetCorrections.C:423
GetCorrections.C:424
GetCorrections.C:425
GetCorrections.C:426
GetCorrections.C:427
GetCorrections.C:428
GetCorrections.C:429
GetCorrections.C:430
GetCorrections.C:431
GetCorrections.C:432
GetCorrections.C:433
GetCorrections.C:434
GetCorrections.C:435
GetCorrections.C:436
GetCorrections.C:437
GetCorrections.C:438
GetCorrections.C:439
GetCorrections.C:440
GetCorrections.C:441
GetCorrections.C:442
GetCorrections.C:443
GetCorrections.C:444
GetCorrections.C:445
GetCorrections.C:446
GetCorrections.C:447
GetCorrections.C:448
GetCorrections.C:449
GetCorrections.C:450
GetCorrections.C:451
GetCorrections.C:452
GetCorrections.C:453
GetCorrections.C:454
GetCorrections.C:455
GetCorrections.C:456
GetCorrections.C:457
GetCorrections.C:458
GetCorrections.C:459
GetCorrections.C:460
GetCorrections.C:461
GetCorrections.C:462
GetCorrections.C:463
GetCorrections.C:464
GetCorrections.C:465
GetCorrections.C:466
GetCorrections.C:467
GetCorrections.C:468
GetCorrections.C:469
GetCorrections.C:470
GetCorrections.C:471
GetCorrections.C:472
GetCorrections.C:473
GetCorrections.C:474
GetCorrections.C:475
GetCorrections.C:476
GetCorrections.C:477
GetCorrections.C:478
GetCorrections.C:479
GetCorrections.C:480
GetCorrections.C:481
GetCorrections.C:482
GetCorrections.C:483
GetCorrections.C:484
GetCorrections.C:485
GetCorrections.C:486
GetCorrections.C:487
GetCorrections.C:488
GetCorrections.C:489
GetCorrections.C:490
GetCorrections.C:491
GetCorrections.C:492
GetCorrections.C:493
GetCorrections.C:494
GetCorrections.C:495
GetCorrections.C:496
GetCorrections.C:497
GetCorrections.C:498
GetCorrections.C:499
GetCorrections.C:500
GetCorrections.C:501
GetCorrections.C:502
GetCorrections.C:503
GetCorrections.C:504
GetCorrections.C:505
GetCorrections.C:506
GetCorrections.C:507
GetCorrections.C:508
GetCorrections.C:509
GetCorrections.C:510
GetCorrections.C:511
GetCorrections.C:512
GetCorrections.C:513
GetCorrections.C:514
GetCorrections.C:515
GetCorrections.C:516
GetCorrections.C:517
GetCorrections.C:518
GetCorrections.C:519
GetCorrections.C:520
GetCorrections.C:521
GetCorrections.C:522
GetCorrections.C:523
GetCorrections.C:524
GetCorrections.C:525
GetCorrections.C:526
GetCorrections.C:527
GetCorrections.C:528
GetCorrections.C:529
GetCorrections.C:530
GetCorrections.C:531
GetCorrections.C:532
GetCorrections.C:533
GetCorrections.C:534
GetCorrections.C:535
GetCorrections.C:536
GetCorrections.C:537
GetCorrections.C:538
GetCorrections.C:539
GetCorrections.C:540
GetCorrections.C:541
GetCorrections.C:542
GetCorrections.C:543
GetCorrections.C:544
GetCorrections.C:545
GetCorrections.C:546
GetCorrections.C:547
GetCorrections.C:548
GetCorrections.C:549
GetCorrections.C:550
GetCorrections.C:551
GetCorrections.C:552
GetCorrections.C:553
GetCorrections.C:554
GetCorrections.C:555
GetCorrections.C:556
GetCorrections.C:557
GetCorrections.C:558
GetCorrections.C:559
GetCorrections.C:560
GetCorrections.C:561
GetCorrections.C:562
GetCorrections.C:563
GetCorrections.C:564
GetCorrections.C:565
GetCorrections.C:566
GetCorrections.C:567
GetCorrections.C:568
GetCorrections.C:569
GetCorrections.C:570
GetCorrections.C:571
GetCorrections.C:572
GetCorrections.C:573
GetCorrections.C:574
GetCorrections.C:575
GetCorrections.C:576
GetCorrections.C:577
GetCorrections.C:578
GetCorrections.C:579
GetCorrections.C:580
GetCorrections.C:581
GetCorrections.C:582
GetCorrections.C:583
GetCorrections.C:584
GetCorrections.C:585
GetCorrections.C:586
GetCorrections.C:587
GetCorrections.C:588
GetCorrections.C:589
GetCorrections.C:590
GetCorrections.C:591
GetCorrections.C:592
GetCorrections.C:593
GetCorrections.C:594
GetCorrections.C:595
GetCorrections.C:596
GetCorrections.C:597
GetCorrections.C:598
GetCorrections.C:599
GetCorrections.C:600
GetCorrections.C:601
GetCorrections.C:602
GetCorrections.C:603
GetCorrections.C:604
GetCorrections.C:605
GetCorrections.C:606
GetCorrections.C:607
GetCorrections.C:608
GetCorrections.C:609
GetCorrections.C:610
GetCorrections.C:611
GetCorrections.C:612
GetCorrections.C:613
GetCorrections.C:614
GetCorrections.C:615
GetCorrections.C:616
GetCorrections.C:617
GetCorrections.C:618
GetCorrections.C:619
GetCorrections.C:620
GetCorrections.C:621
GetCorrections.C:622
GetCorrections.C:623
GetCorrections.C:624
GetCorrections.C:625
GetCorrections.C:626
GetCorrections.C:627
GetCorrections.C:628
GetCorrections.C:629
GetCorrections.C:630
GetCorrections.C:631
GetCorrections.C:632
GetCorrections.C:633
GetCorrections.C:634
GetCorrections.C:635
GetCorrections.C:636
GetCorrections.C:637
GetCorrections.C:638
GetCorrections.C:639
GetCorrections.C:640
GetCorrections.C:641
GetCorrections.C:642
GetCorrections.C:643
GetCorrections.C:644
GetCorrections.C:645
GetCorrections.C:646
GetCorrections.C:647
GetCorrections.C:648
GetCorrections.C:649
GetCorrections.C:650
GetCorrections.C:651
GetCorrections.C:652
GetCorrections.C:653
GetCorrections.C:654
GetCorrections.C:655
GetCorrections.C:656
GetCorrections.C:657
GetCorrections.C:658
GetCorrections.C:659
GetCorrections.C:660
GetCorrections.C:661
GetCorrections.C:662
GetCorrections.C:663
GetCorrections.C:664
GetCorrections.C:665
GetCorrections.C:666
GetCorrections.C:667
GetCorrections.C:668
GetCorrections.C:669
GetCorrections.C:670
GetCorrections.C:671
GetCorrections.C:672
GetCorrections.C:673
GetCorrections.C:674
GetCorrections.C:675
GetCorrections.C:676
GetCorrections.C:677
GetCorrections.C:678
GetCorrections.C:679
GetCorrections.C:680
GetCorrections.C:681
GetCorrections.C:682
GetCorrections.C:683
GetCorrections.C:684
GetCorrections.C:685
GetCorrections.C:686
GetCorrections.C:687
GetCorrections.C:688
GetCorrections.C:689
GetCorrections.C:690
GetCorrections.C:691
GetCorrections.C:692
GetCorrections.C:693
GetCorrections.C:694
GetCorrections.C:695
GetCorrections.C:696
GetCorrections.C:697
GetCorrections.C:698
GetCorrections.C:699
GetCorrections.C:700
GetCorrections.C:701
GetCorrections.C:702
GetCorrections.C:703
GetCorrections.C:704
GetCorrections.C:705
GetCorrections.C:706
GetCorrections.C:707
GetCorrections.C:708
GetCorrections.C:709
GetCorrections.C:710
GetCorrections.C:711
GetCorrections.C:712
GetCorrections.C:713
GetCorrections.C:714
GetCorrections.C:715
GetCorrections.C:716
GetCorrections.C:717
GetCorrections.C:718
GetCorrections.C:719
GetCorrections.C:720
GetCorrections.C:721
GetCorrections.C:722
GetCorrections.C:723
GetCorrections.C:724
GetCorrections.C:725
GetCorrections.C:726
GetCorrections.C:727
GetCorrections.C:728
GetCorrections.C:729
GetCorrections.C:730
GetCorrections.C:731
GetCorrections.C:732
GetCorrections.C:733
GetCorrections.C:734
GetCorrections.C:735
GetCorrections.C:736
GetCorrections.C:737
GetCorrections.C:738
GetCorrections.C:739
GetCorrections.C:740
GetCorrections.C:741
GetCorrections.C:742
GetCorrections.C:743
GetCorrections.C:744
GetCorrections.C:745
GetCorrections.C:746
GetCorrections.C:747
GetCorrections.C:748
GetCorrections.C:749
GetCorrections.C:750
GetCorrections.C:751
GetCorrections.C:752
GetCorrections.C:753
GetCorrections.C:754
GetCorrections.C:755
GetCorrections.C:756
GetCorrections.C:757
GetCorrections.C:758
GetCorrections.C:759
GetCorrections.C:760
GetCorrections.C:761
GetCorrections.C:762
GetCorrections.C:763
GetCorrections.C:764
GetCorrections.C:765
GetCorrections.C:766
GetCorrections.C:767
GetCorrections.C:768
GetCorrections.C:769
GetCorrections.C:770
GetCorrections.C:771
GetCorrections.C:772
GetCorrections.C:773
GetCorrections.C:774
GetCorrections.C:775
GetCorrections.C:776
GetCorrections.C:777
GetCorrections.C:778
GetCorrections.C:779
GetCorrections.C:780
GetCorrections.C:781
GetCorrections.C:782
GetCorrections.C:783
GetCorrections.C:784
GetCorrections.C:785
GetCorrections.C:786
GetCorrections.C:787
GetCorrections.C:788
GetCorrections.C:789
GetCorrections.C:790
GetCorrections.C:791
GetCorrections.C:792
GetCorrections.C:793
GetCorrections.C:794
GetCorrections.C:795
GetCorrections.C:796
GetCorrections.C:797
GetCorrections.C:798
GetCorrections.C:799
GetCorrections.C:800
GetCorrections.C:801
GetCorrections.C:802
GetCorrections.C:803
GetCorrections.C:804
GetCorrections.C:805
GetCorrections.C:806
GetCorrections.C:807
GetCorrections.C:808
GetCorrections.C:809
GetCorrections.C:810
GetCorrections.C:811
GetCorrections.C:812
GetCorrections.C:813
GetCorrections.C:814
GetCorrections.C:815
GetCorrections.C:816
GetCorrections.C:817
GetCorrections.C:818
GetCorrections.C:819
GetCorrections.C:820
GetCorrections.C:821
GetCorrections.C:822
GetCorrections.C:823
GetCorrections.C:824
GetCorrections.C:825
GetCorrections.C:826
GetCorrections.C:827
GetCorrections.C:828
GetCorrections.C:829
GetCorrections.C:830
GetCorrections.C:831
GetCorrections.C:832
GetCorrections.C:833
GetCorrections.C:834
GetCorrections.C:835
GetCorrections.C:836
GetCorrections.C:837
GetCorrections.C:838
GetCorrections.C:839
GetCorrections.C:840
GetCorrections.C:841
GetCorrections.C:842
GetCorrections.C:843
GetCorrections.C:844
GetCorrections.C:845
GetCorrections.C:846
GetCorrections.C:847
GetCorrections.C:848
GetCorrections.C:849
GetCorrections.C:850
GetCorrections.C:851
GetCorrections.C:852
GetCorrections.C:853
GetCorrections.C:854
GetCorrections.C:855
GetCorrections.C:856
GetCorrections.C:857
GetCorrections.C:858
GetCorrections.C:859
GetCorrections.C:860
GetCorrections.C:861
GetCorrections.C:862
GetCorrections.C:863
GetCorrections.C:864
GetCorrections.C:865
GetCorrections.C:866
GetCorrections.C:867
GetCorrections.C:868
GetCorrections.C:869
GetCorrections.C:870
GetCorrections.C:871
GetCorrections.C:872
GetCorrections.C:873
GetCorrections.C:874
GetCorrections.C:875
GetCorrections.C:876
GetCorrections.C:877
GetCorrections.C:878
GetCorrections.C:879
GetCorrections.C:880
GetCorrections.C:881
GetCorrections.C:882
GetCorrections.C:883
GetCorrections.C:884
GetCorrections.C:885
GetCorrections.C:886
GetCorrections.C:887
GetCorrections.C:888
GetCorrections.C:889
GetCorrections.C:890
GetCorrections.C:891
GetCorrections.C:892
GetCorrections.C:893
GetCorrections.C:894
GetCorrections.C:895
GetCorrections.C:896
GetCorrections.C:897
GetCorrections.C:898
GetCorrections.C:899
GetCorrections.C:900
GetCorrections.C:901
GetCorrections.C:902
GetCorrections.C:903
GetCorrections.C:904
GetCorrections.C:905
GetCorrections.C:906
GetCorrections.C:907
GetCorrections.C:908
GetCorrections.C:909
GetCorrections.C:910
GetCorrections.C:911
GetCorrections.C:912
GetCorrections.C:913
GetCorrections.C:914
GetCorrections.C:915
GetCorrections.C:916
GetCorrections.C:917
GetCorrections.C:918
GetCorrections.C:919
GetCorrections.C:920
GetCorrections.C:921
GetCorrections.C:922
GetCorrections.C:923
GetCorrections.C:924
GetCorrections.C:925
GetCorrections.C:926
GetCorrections.C:927
GetCorrections.C:928
GetCorrections.C:929
GetCorrections.C:930
GetCorrections.C:931
GetCorrections.C:932
GetCorrections.C:933
GetCorrections.C:934
GetCorrections.C:935
GetCorrections.C:936
GetCorrections.C:937
GetCorrections.C:938
GetCorrections.C:939
GetCorrections.C:940
GetCorrections.C:941
GetCorrections.C:942
GetCorrections.C:943
GetCorrections.C:944
GetCorrections.C:945
GetCorrections.C:946
GetCorrections.C:947
GetCorrections.C:948
GetCorrections.C:949
GetCorrections.C:950
GetCorrections.C:951
GetCorrections.C:952
GetCorrections.C:953
GetCorrections.C:954
GetCorrections.C:955
GetCorrections.C:956
GetCorrections.C:957
GetCorrections.C:958
GetCorrections.C:959
GetCorrections.C:960
GetCorrections.C:961
GetCorrections.C:962
GetCorrections.C:963
GetCorrections.C:964
GetCorrections.C:965
GetCorrections.C:966
GetCorrections.C:967
GetCorrections.C:968
GetCorrections.C:969
GetCorrections.C:970
GetCorrections.C:971
GetCorrections.C:972
GetCorrections.C:973
GetCorrections.C:974
GetCorrections.C:975
GetCorrections.C:976
GetCorrections.C:977
GetCorrections.C:978
GetCorrections.C:979
GetCorrections.C:980
GetCorrections.C:981
GetCorrections.C:982
GetCorrections.C:983
GetCorrections.C:984
GetCorrections.C:985
GetCorrections.C:986
GetCorrections.C:987
GetCorrections.C:988
GetCorrections.C:989
GetCorrections.C:990
GetCorrections.C:991
GetCorrections.C:992
GetCorrections.C:993
GetCorrections.C:994
GetCorrections.C:995
GetCorrections.C:996
GetCorrections.C:997
GetCorrections.C:998
GetCorrections.C:999
GetCorrections.C:1000
GetCorrections.C:1001
GetCorrections.C:1002
GetCorrections.C:1003
GetCorrections.C:1004
GetCorrections.C:1005
GetCorrections.C:1006
GetCorrections.C:1007
GetCorrections.C:1008
GetCorrections.C:1009
GetCorrections.C:1010
GetCorrections.C:1011
GetCorrections.C:1012
GetCorrections.C:1013
GetCorrections.C:1014
GetCorrections.C:1015
GetCorrections.C:1016
GetCorrections.C:1017
GetCorrections.C:1018
GetCorrections.C:1019
GetCorrections.C:1020
GetCorrections.C:1021
GetCorrections.C:1022
GetCorrections.C:1023
GetCorrections.C:1024
GetCorrections.C:1025
GetCorrections.C:1026
GetCorrections.C:1027
GetCorrections.C:1028
GetCorrections.C:1029
GetCorrections.C:1030
GetCorrections.C:1031
GetCorrections.C:1032
GetCorrections.C:1033
GetCorrections.C:1034
GetCorrections.C:1035
GetCorrections.C:1036
GetCorrections.C:1037
GetCorrections.C:1038
GetCorrections.C:1039
GetCorrections.C:1040
GetCorrections.C:1041
GetCorrections.C:1042
GetCorrections.C:1043
GetCorrections.C:1044
GetCorrections.C:1045
GetCorrections.C:1046
GetCorrections.C:1047
GetCorrections.C:1048
GetCorrections.C:1049
GetCorrections.C:1050
GetCorrections.C:1051
GetCorrections.C:1052
GetCorrections.C:1053
GetCorrections.C:1054
GetCorrections.C:1055
GetCorrections.C:1056
GetCorrections.C:1057
GetCorrections.C:1058
GetCorrections.C:1059
GetCorrections.C:1060
GetCorrections.C:1061
GetCorrections.C:1062
GetCorrections.C:1063
GetCorrections.C:1064
GetCorrections.C:1065
GetCorrections.C:1066
GetCorrections.C:1067
GetCorrections.C:1068
GetCorrections.C:1069
GetCorrections.C:1070
GetCorrections.C:1071
GetCorrections.C:1072
GetCorrections.C:1073
GetCorrections.C:1074
GetCorrections.C:1075
GetCorrections.C:1076
GetCorrections.C:1077
GetCorrections.C:1078
GetCorrections.C:1079
GetCorrections.C:1080
GetCorrections.C:1081
GetCorrections.C:1082
GetCorrections.C:1083
GetCorrections.C:1084
GetCorrections.C:1085
GetCorrections.C:1086
GetCorrections.C:1087
GetCorrections.C:1088
GetCorrections.C:1089
GetCorrections.C:1090
GetCorrections.C:1091
GetCorrections.C:1092
GetCorrections.C:1093
GetCorrections.C:1094
GetCorrections.C:1095
GetCorrections.C:1096
GetCorrections.C:1097
GetCorrections.C:1098
GetCorrections.C:1099
GetCorrections.C:1100
GetCorrections.C:1101
GetCorrections.C:1102
GetCorrections.C:1103
GetCorrections.C:1104
GetCorrections.C:1105
GetCorrections.C:1106
GetCorrections.C:1107
GetCorrections.C:1108
GetCorrections.C:1109
GetCorrections.C:1110
GetCorrections.C:1111
GetCorrections.C:1112
GetCorrections.C:1113
GetCorrections.C:1114
GetCorrections.C:1115
GetCorrections.C:1116
GetCorrections.C:1117
GetCorrections.C:1118
GetCorrections.C:1119
GetCorrections.C:1120
GetCorrections.C:1121
GetCorrections.C:1122
GetCorrections.C:1123
GetCorrections.C:1124
GetCorrections.C:1125
GetCorrections.C:1126
GetCorrections.C:1127
GetCorrections.C:1128
GetCorrections.C:1129
GetCorrections.C:1130
GetCorrections.C:1131
GetCorrections.C:1132
GetCorrections.C:1133
GetCorrections.C:1134
GetCorrections.C:1135
GetCorrections.C:1136
GetCorrections.C:1137
GetCorrections.C:1138
GetCorrections.C:1139
GetCorrections.C:1140
GetCorrections.C:1141
GetCorrections.C:1142
GetCorrections.C:1143
GetCorrections.C:1144
GetCorrections.C:1145
GetCorrections.C:1146
GetCorrections.C:1147
GetCorrections.C:1148
GetCorrections.C:1149
GetCorrections.C:1150
GetCorrections.C:1151
GetCorrections.C:1152
GetCorrections.C:1153
GetCorrections.C:1154
GetCorrections.C:1155
GetCorrections.C:1156
GetCorrections.C:1157
GetCorrections.C:1158
GetCorrections.C:1159
GetCorrections.C:1160
GetCorrections.C:1161
GetCorrections.C:1162
GetCorrections.C:1163
GetCorrections.C:1164
GetCorrections.C:1165
GetCorrections.C:1166
GetCorrections.C:1167
GetCorrections.C:1168
GetCorrections.C:1169
GetCorrections.C:1170
GetCorrections.C:1171
GetCorrections.C:1172
GetCorrections.C:1173
GetCorrections.C:1174
GetCorrections.C:1175
GetCorrections.C:1176
GetCorrections.C:1177
GetCorrections.C:1178
GetCorrections.C:1179
GetCorrections.C:1180
GetCorrections.C:1181
GetCorrections.C:1182
GetCorrections.C:1183
GetCorrections.C:1184
GetCorrections.C:1185
GetCorrections.C:1186
GetCorrections.C:1187
GetCorrections.C:1188
GetCorrections.C:1189
GetCorrections.C:1190
GetCorrections.C:1191
GetCorrections.C:1192
GetCorrections.C:1193
GetCorrections.C:1194
GetCorrections.C:1195
GetCorrections.C:1196
GetCorrections.C:1197
GetCorrections.C:1198
GetCorrections.C:1199
GetCorrections.C:1200
GetCorrections.C:1201
GetCorrections.C:1202
GetCorrections.C:1203
GetCorrections.C:1204
GetCorrections.C:1205
GetCorrections.C:1206
GetCorrections.C:1207
GetCorrections.C:1208
GetCorrections.C:1209
GetCorrections.C:1210
GetCorrections.C:1211
GetCorrections.C:1212
GetCorrections.C:1213
GetCorrections.C:1214
GetCorrections.C:1215
GetCorrections.C:1216
GetCorrections.C:1217
GetCorrections.C:1218
GetCorrections.C:1219
GetCorrections.C:1220
GetCorrections.C:1221
GetCorrections.C:1222
GetCorrections.C:1223
GetCorrections.C:1224
GetCorrections.C:1225
GetCorrections.C:1226
GetCorrections.C:1227
GetCorrections.C:1228
GetCorrections.C:1229
GetCorrections.C:1230
GetCorrections.C:1231
GetCorrections.C:1232
GetCorrections.C:1233
GetCorrections.C:1234
GetCorrections.C:1235
GetCorrections.C:1236
GetCorrections.C:1237
GetCorrections.C:1238
GetCorrections.C:1239
GetCorrections.C:1240
GetCorrections.C:1241
GetCorrections.C:1242
GetCorrections.C:1243
GetCorrections.C:1244
GetCorrections.C:1245
GetCorrections.C:1246
GetCorrections.C:1247
GetCorrections.C:1248
GetCorrections.C:1249
GetCorrections.C:1250
GetCorrections.C:1251
GetCorrections.C:1252
GetCorrections.C:1253
GetCorrections.C:1254
GetCorrections.C:1255
GetCorrections.C:1256
GetCorrections.C:1257
GetCorrections.C:1258
GetCorrections.C:1259
GetCorrections.C:1260
GetCorrections.C:1261
GetCorrections.C:1262
GetCorrections.C:1263
GetCorrections.C:1264
GetCorrections.C:1265
GetCorrections.C:1266
GetCorrections.C:1267
GetCorrections.C:1268
GetCorrections.C:1269
GetCorrections.C:1270
GetCorrections.C:1271
GetCorrections.C:1272
GetCorrections.C:1273
GetCorrections.C:1274
GetCorrections.C:1275
GetCorrections.C:1276
GetCorrections.C:1277
GetCorrections.C:1278
GetCorrections.C:1279
GetCorrections.C:1280
GetCorrections.C:1281
GetCorrections.C:1282
GetCorrections.C:1283
GetCorrections.C:1284
GetCorrections.C:1285
GetCorrections.C:1286
GetCorrections.C:1287
GetCorrections.C:1288
GetCorrections.C:1289
GetCorrections.C:1290
GetCorrections.C:1291
GetCorrections.C:1292
GetCorrections.C:1293
GetCorrections.C:1294
GetCorrections.C:1295
GetCorrections.C:1296
GetCorrections.C:1297
GetCorrections.C:1298
GetCorrections.C:1299
GetCorrections.C:1300
GetCorrections.C:1301
GetCorrections.C:1302
GetCorrections.C:1303
GetCorrections.C:1304
GetCorrections.C:1305
GetCorrections.C:1306
GetCorrections.C:1307
GetCorrections.C:1308
GetCorrections.C:1309
GetCorrections.C:1310
GetCorrections.C:1311
GetCorrections.C:1312
GetCorrections.C:1313
GetCorrections.C:1314
GetCorrections.C:1315
GetCorrections.C:1316
GetCorrections.C:1317
GetCorrections.C:1318
GetCorrections.C:1319
GetCorrections.C:1320
GetCorrections.C:1321
GetCorrections.C:1322
GetCorrections.C:1323
GetCorrections.C:1324
GetCorrections.C:1325
GetCorrections.C:1326
GetCorrections.C:1327
GetCorrections.C:1328
GetCorrections.C:1329
GetCorrections.C:1330
GetCorrections.C:1331
GetCorrections.C:1332
GetCorrections.C:1333
GetCorrections.C:1334
GetCorrections.C:1335
GetCorrections.C:1336
GetCorrections.C:1337
GetCorrections.C:1338
GetCorrections.C:1339
GetCorrections.C:1340
GetCorrections.C:1341
GetCorrections.C:1342
GetCorrections.C:1343
GetCorrections.C:1344
GetCorrections.C:1345
GetCorrections.C:1346
GetCorrections.C:1347
GetCorrections.C:1348
GetCorrections.C:1349
GetCorrections.C:1350
GetCorrections.C:1351
GetCorrections.C:1352
GetCorrections.C:1353
GetCorrections.C:1354
GetCorrections.C:1355
GetCorrections.C:1356
GetCorrections.C:1357
GetCorrections.C:1358
GetCorrections.C:1359
GetCorrections.C:1360
GetCorrections.C:1361
GetCorrections.C:1362
GetCorrections.C:1363
GetCorrections.C:1364
GetCorrections.C:1365
GetCorrections.C:1366
GetCorrections.C:1367
GetCorrections.C:1368
GetCorrections.C:1369
GetCorrections.C:1370
GetCorrections.C:1371
GetCorrections.C:1372
GetCorrections.C:1373
GetCorrections.C:1374
GetCorrections.C:1375
GetCorrections.C:1376
GetCorrections.C:1377
GetCorrections.C:1378
GetCorrections.C:1379
GetCorrections.C:1380
GetCorrections.C:1381
GetCorrections.C:1382
GetCorrections.C:1383
GetCorrections.C:1384
GetCorrections.C:1385
GetCorrections.C:1386
GetCorrections.C:1387
GetCorrections.C:1388
GetCorrections.C:1389
GetCorrections.C:1390
GetCorrections.C:1391
GetCorrections.C:1392
GetCorrections.C:1393
GetCorrections.C:1394
GetCorrections.C:1395
GetCorrections.C:1396
GetCorrections.C:1397
GetCorrections.C:1398
GetCorrections.C:1399
GetCorrections.C:1400
GetCorrections.C:1401
GetCorrections.C:1402
GetCorrections.C:1403
GetCorrections.C:1404
GetCorrections.C:1405
GetCorrections.C:1406
GetCorrections.C:1407
GetCorrections.C:1408
GetCorrections.C:1409
GetCorrections.C:1410
GetCorrections.C:1411
GetCorrections.C:1412
GetCorrections.C:1413
GetCorrections.C:1414
GetCorrections.C:1415
GetCorrections.C:1416
GetCorrections.C:1417
GetCorrections.C:1418
GetCorrections.C:1419
GetCorrections.C:1420
GetCorrections.C:1421
GetCorrections.C:1422
GetCorrections.C:1423
GetCorrections.C:1424
GetCorrections.C:1425
GetCorrections.C:1426
GetCorrections.C:1427
GetCorrections.C:1428
GetCorrections.C:1429
GetCorrections.C:1430
GetCorrections.C:1431
GetCorrections.C:1432
GetCorrections.C:1433
GetCorrections.C:1434
GetCorrections.C:1435
GetCorrections.C:1436
GetCorrections.C:1437
GetCorrections.C:1438
GetCorrections.C:1439
GetCorrections.C:1440
GetCorrections.C:1441
GetCorrections.C:1442
GetCorrections.C:1443
GetCorrections.C:1444
GetCorrections.C:1445
GetCorrections.C:1446
GetCorrections.C:1447
GetCorrections.C:1448
GetCorrections.C:1449
GetCorrections.C:1450
GetCorrections.C:1451
GetCorrections.C:1452
GetCorrections.C:1453
GetCorrections.C:1454
GetCorrections.C:1455
GetCorrections.C:1456
GetCorrections.C:1457
GetCorrections.C:1458
GetCorrections.C:1459
GetCorrections.C:1460
GetCorrections.C:1461
GetCorrections.C:1462
GetCorrections.C:1463
GetCorrections.C:1464
GetCorrections.C:1465
GetCorrections.C:1466
GetCorrections.C:1467
GetCorrections.C:1468
GetCorrections.C:1469
GetCorrections.C:1470
GetCorrections.C:1471
GetCorrections.C:1472
GetCorrections.C:1473
GetCorrections.C:1474
GetCorrections.C:1475
GetCorrections.C:1476
GetCorrections.C:1477
GetCorrections.C:1478
GetCorrections.C:1479
GetCorrections.C:1480
GetCorrections.C:1481
GetCorrections.C:1482
GetCorrections.C:1483
GetCorrections.C:1484
GetCorrections.C:1485
GetCorrections.C:1486
GetCorrections.C:1487
GetCorrections.C:1488
GetCorrections.C:1489
GetCorrections.C:1490
GetCorrections.C:1491
GetCorrections.C:1492
GetCorrections.C:1493
GetCorrections.C:1494
GetCorrections.C:1495
GetCorrections.C:1496
GetCorrections.C:1497
GetCorrections.C:1498
GetCorrections.C:1499
GetCorrections.C:1500
GetCorrections.C:1501
GetCorrections.C:1502
GetCorrections.C:1503
GetCorrections.C:1504
GetCorrections.C:1505
GetCorrections.C:1506
GetCorrections.C:1507
GetCorrections.C:1508
GetCorrections.C:1509
GetCorrections.C:1510
GetCorrections.C:1511
GetCorrections.C:1512
GetCorrections.C:1513
GetCorrections.C:1514
GetCorrections.C:1515
GetCorrections.C:1516
GetCorrections.C:1517
GetCorrections.C:1518
GetCorrections.C:1519
GetCorrections.C:1520
GetCorrections.C:1521
GetCorrections.C:1522
GetCorrections.C:1523
GetCorrections.C:1524
GetCorrections.C:1525
GetCorrections.C:1526
GetCorrections.C:1527
GetCorrections.C:1528
GetCorrections.C:1529
GetCorrections.C:1530
GetCorrections.C:1531
GetCorrections.C:1532
GetCorrections.C:1533
GetCorrections.C:1534
GetCorrections.C:1535
GetCorrections.C:1536
GetCorrections.C:1537
GetCorrections.C:1538
GetCorrections.C:1539
GetCorrections.C:1540
GetCorrections.C:1541
GetCorrections.C:1542
GetCorrections.C:1543
GetCorrections.C:1544
GetCorrections.C:1545
GetCorrections.C:1546
GetCorrections.C:1547
GetCorrections.C:1548
GetCorrections.C:1549
GetCorrections.C:1550
GetCorrections.C:1551
GetCorrections.C:1552
GetCorrections.C:1553
GetCorrections.C:1554
GetCorrections.C:1555
GetCorrections.C:1556
GetCorrections.C:1557
GetCorrections.C:1558
GetCorrections.C:1559
GetCorrections.C:1560
GetCorrections.C:1561
GetCorrections.C:1562
GetCorrections.C:1563
GetCorrections.C:1564
GetCorrections.C:1565
GetCorrections.C:1566
GetCorrections.C:1567
GetCorrections.C:1568
GetCorrections.C:1569
GetCorrections.C:1570
GetCorrections.C:1571
GetCorrections.C:1572
GetCorrections.C:1573
GetCorrections.C:1574
GetCorrections.C:1575
GetCorrections.C:1576
GetCorrections.C:1577
GetCorrections.C:1578
GetCorrections.C:1579
GetCorrections.C:1580
GetCorrections.C:1581
GetCorrections.C:1582
GetCorrections.C:1583
GetCorrections.C:1584
GetCorrections.C:1585
GetCorrections.C:1586
GetCorrections.C:1587
GetCorrections.C:1588
GetCorrections.C:1589
GetCorrections.C:1590
GetCorrections.C:1591
GetCorrections.C:1592
GetCorrections.C:1593
GetCorrections.C:1594
GetCorrections.C:1595
GetCorrections.C:1596
GetCorrections.C:1597
GetCorrections.C:1598
GetCorrections.C:1599
GetCorrections.C:1600
GetCorrections.C:1601
GetCorrections.C:1602
GetCorrections.C:1603
GetCorrections.C:1604
GetCorrections.C:1605
GetCorrections.C:1606
GetCorrections.C:1607
GetCorrections.C:1608
GetCorrections.C:1609
GetCorrections.C:1610
GetCorrections.C:1611
GetCorrections.C:1612
GetCorrections.C:1613
GetCorrections.C:1614
GetCorrections.C:1615
GetCorrections.C:1616
GetCorrections.C:1617
GetCorrections.C:1618
GetCorrections.C:1619
GetCorrections.C:1620
GetCorrections.C:1621
GetCorrections.C:1622
GetCorrections.C:1623
GetCorrections.C:1624
GetCorrections.C:1625
GetCorrections.C:1626
GetCorrections.C:1627
GetCorrections.C:1628
GetCorrections.C:1629
GetCorrections.C:1630
GetCorrections.C:1631
GetCorrections.C:1632
GetCorrections.C:1633
GetCorrections.C:1634
GetCorrections.C:1635
GetCorrections.C:1636
GetCorrections.C:1637
GetCorrections.C:1638
GetCorrections.C:1639
GetCorrections.C:1640
GetCorrections.C:1641
GetCorrections.C:1642
GetCorrections.C:1643
GetCorrections.C:1644
GetCorrections.C:1645
GetCorrections.C:1646
GetCorrections.C:1647
GetCorrections.C:1648
GetCorrections.C:1649
GetCorrections.C:1650
GetCorrections.C:1651
GetCorrections.C:1652
GetCorrections.C:1653
GetCorrections.C:1654
GetCorrections.C:1655
GetCorrections.C:1656
GetCorrections.C:1657
GetCorrections.C:1658
GetCorrections.C:1659
GetCorrections.C:1660
GetCorrections.C:1661
GetCorrections.C:1662
GetCorrections.C:1663
GetCorrections.C:1664
GetCorrections.C:1665
GetCorrections.C:1666
GetCorrections.C:1667
GetCorrections.C:1668
GetCorrections.C:1669
GetCorrections.C:1670
GetCorrections.C:1671
GetCorrections.C:1672
GetCorrections.C:1673
GetCorrections.C:1674
GetCorrections.C:1675
GetCorrections.C:1676
GetCorrections.C:1677
GetCorrections.C:1678
GetCorrections.C:1679
GetCorrections.C:1680
GetCorrections.C:1681
GetCorrections.C:1682
GetCorrections.C:1683
GetCorrections.C:1684
GetCorrections.C:1685
GetCorrections.C:1686
GetCorrections.C:1687
GetCorrections.C:1688
GetCorrections.C:1689
GetCorrections.C:1690
GetCorrections.C:1691
GetCorrections.C:1692
GetCorrections.C:1693
GetCorrections.C:1694
GetCorrections.C:1695
GetCorrections.C:1696
GetCorrections.C:1697
GetCorrections.C:1698
GetCorrections.C:1699
GetCorrections.C:1700
GetCorrections.C:1701
GetCorrections.C:1702
GetCorrections.C:1703
GetCorrections.C:1704
GetCorrections.C:1705
GetCorrections.C:1706
GetCorrections.C:1707
GetCorrections.C:1708
GetCorrections.C:1709
GetCorrections.C:1710
GetCorrections.C:1711
GetCorrections.C:1712
GetCorrections.C:1713
GetCorrections.C:1714
GetCorrections.C:1715
GetCorrections.C:1716
GetCorrections.C:1717
GetCorrections.C:1718
GetCorrections.C:1719
GetCorrections.C:1720
GetCorrections.C:1721
GetCorrections.C:1722
GetCorrections.C:1723
GetCorrections.C:1724
GetCorrections.C:1725
GetCorrections.C:1726
GetCorrections.C:1727
GetCorrections.C:1728
GetCorrections.C:1729
GetCorrections.C:1730
GetCorrections.C:1731
GetCorrections.C:1732
GetCorrections.C:1733
GetCorrections.C:1734
GetCorrections.C:1735
GetCorrections.C:1736
GetCorrections.C:1737
GetCorrections.C:1738
GetCorrections.C:1739
GetCorrections.C:1740
GetCorrections.C:1741
GetCorrections.C:1742
GetCorrections.C:1743
GetCorrections.C:1744
GetCorrections.C:1745
GetCorrections.C:1746
GetCorrections.C:1747
GetCorrections.C:1748
GetCorrections.C:1749
GetCorrections.C:1750
GetCorrections.C:1751
GetCorrections.C:1752
GetCorrections.C:1753
GetCorrections.C:1754
GetCorrections.C:1755
GetCorrections.C:1756
GetCorrections.C:1757
GetCorrections.C:1758
GetCorrections.C:1759
GetCorrections.C:1760
GetCorrections.C:1761
GetCorrections.C:1762
GetCorrections.C:1763
GetCorrections.C:1764
GetCorrections.C:1765
GetCorrections.C:1766
GetCorrections.C:1767
GetCorrections.C:1768
GetCorrections.C:1769
GetCorrections.C:1770
GetCorrections.C:1771
GetCorrections.C:1772
GetCorrections.C:1773
GetCorrections.C:1774
GetCorrections.C:1775
GetCorrections.C:1776
GetCorrections.C:1777
GetCorrections.C:1778
GetCorrections.C:1779
GetCorrections.C:1780
GetCorrections.C:1781
GetCorrections.C:1782
GetCorrections.C:1783
GetCorrections.C:1784
GetCorrections.C:1785
GetCorrections.C:1786
GetCorrections.C:1787
GetCorrections.C:1788
GetCorrections.C:1789
GetCorrections.C:1790
GetCorrections.C:1791
GetCorrections.C:1792
GetCorrections.C:1793
GetCorrections.C:1794
GetCorrections.C:1795
GetCorrections.C:1796
GetCorrections.C:1797
GetCorrections.C:1798
GetCorrections.C:1799
GetCorrections.C:1800
GetCorrections.C:1801
GetCorrections.C:1802
GetCorrections.C:1803
GetCorrections.C:1804
GetCorrections.C:1805
GetCorrections.C:1806
GetCorrections.C:1807
GetCorrections.C:1808
GetCorrections.C:1809
GetCorrections.C:1810
GetCorrections.C:1811
GetCorrections.C:1812
GetCorrections.C:1813
GetCorrections.C:1814
GetCorrections.C:1815
GetCorrections.C:1816
GetCorrections.C:1817
GetCorrections.C:1818
GetCorrections.C:1819
GetCorrections.C:1820
GetCorrections.C:1821
GetCorrections.C:1822
GetCorrections.C:1823
GetCorrections.C:1824
GetCorrections.C:1825
GetCorrections.C:1826
GetCorrections.C:1827
GetCorrections.C:1828
GetCorrections.C:1829
GetCorrections.C:1830
GetCorrections.C:1831
GetCorrections.C:1832
GetCorrections.C:1833
GetCorrections.C:1834
GetCorrections.C:1835
GetCorrections.C:1836
GetCorrections.C:1837
GetCorrections.C:1838
GetCorrections.C:1839
GetCorrections.C:1840
GetCorrections.C:1841
GetCorrections.C:1842
GetCorrections.C:1843
GetCorrections.C:1844
GetCorrections.C:1845
GetCorrections.C:1846
GetCorrections.C:1847
GetCorrections.C:1848
GetCorrections.C:1849
GetCorrections.C:1850
GetCorrections.C:1851
GetCorrections.C:1852
GetCorrections.C:1853
GetCorrections.C:1854
GetCorrections.C:1855
GetCorrections.C:1856
GetCorrections.C:1857
GetCorrections.C:1858
GetCorrections.C:1859
GetCorrections.C:1860
GetCorrections.C:1861
GetCorrections.C:1862
GetCorrections.C:1863
GetCorrections.C:1864
GetCorrections.C:1865
GetCorrections.C:1866
GetCorrections.C:1867
GetCorrections.C:1868
GetCorrections.C:1869
GetCorrections.C:1870
GetCorrections.C:1871
GetCorrections.C:1872
GetCorrections.C:1873
GetCorrections.C:1874
GetCorrections.C:1875
GetCorrections.C:1876
GetCorrections.C:1877
GetCorrections.C:1878
GetCorrections.C:1879
GetCorrections.C:1880
GetCorrections.C:1881
GetCorrections.C:1882
GetCorrections.C:1883
GetCorrections.C:1884
GetCorrections.C:1885
GetCorrections.C:1886
GetCorrections.C:1887
GetCorrections.C:1888
GetCorrections.C:1889
GetCorrections.C:1890
GetCorrections.C:1891
GetCorrections.C:1892
GetCorrections.C:1893
GetCorrections.C:1894
GetCorrections.C:1895
GetCorrections.C:1896
GetCorrections.C:1897
GetCorrections.C:1898
GetCorrections.C:1899
GetCorrections.C:1900
GetCorrections.C:1901
GetCorrections.C:1902
GetCorrections.C:1903
GetCorrections.C:1904
GetCorrections.C:1905
GetCorrections.C:1906
GetCorrections.C:1907
GetCorrections.C:1908
GetCorrections.C:1909
GetCorrections.C:1910
GetCorrections.C:1911
GetCorrections.C:1912
GetCorrections.C:1913
GetCorrections.C:1914
GetCorrections.C:1915
GetCorrections.C:1916
GetCorrections.C:1917
GetCorrections.C:1918
GetCorrections.C:1919
GetCorrections.C:1920
GetCorrections.C:1921
GetCorrections.C:1922
GetCorrections.C:1923
GetCorrections.C:1924
GetCorrections.C:1925
GetCorrections.C:1926
GetCorrections.C:1927
GetCorrections.C:1928
GetCorrections.C:1929
GetCorrections.C:1930
GetCorrections.C:1931
GetCorrections.C:1932
GetCorrections.C:1933
GetCorrections.C:1934
GetCorrections.C:1935
GetCorrections.C:1936
GetCorrections.C:1937
GetCorrections.C:1938
GetCorrections.C:1939
GetCorrections.C:1940
GetCorrections.C:1941
GetCorrections.C:1942
GetCorrections.C:1943
GetCorrections.C:1944
GetCorrections.C:1945
GetCorrections.C:1946
GetCorrections.C:1947
GetCorrections.C:1948
GetCorrections.C:1949
GetCorrections.C:1950
GetCorrections.C:1951
GetCorrections.C:1952
GetCorrections.C:1953
GetCorrections.C:1954
GetCorrections.C:1955
GetCorrections.C:1956
GetCorrections.C:1957
GetCorrections.C:1958
GetCorrections.C:1959
GetCorrections.C:1960
GetCorrections.C:1961
GetCorrections.C:1962
GetCorrections.C:1963
GetCorrections.C:1964
GetCorrections.C:1965
GetCorrections.C:1966
GetCorrections.C:1967
GetCorrections.C:1968
GetCorrections.C:1969
GetCorrections.C:1970
GetCorrections.C:1971
GetCorrections.C:1972
GetCorrections.C:1973
GetCorrections.C:1974
GetCorrections.C:1975
GetCorrections.C:1976
GetCorrections.C:1977
GetCorrections.C:1978
GetCorrections.C:1979
GetCorrections.C:1980
GetCorrections.C:1981
GetCorrections.C:1982
GetCorrections.C:1983
GetCorrections.C:1984
GetCorrections.C:1985
GetCorrections.C:1986
GetCorrections.C:1987
GetCorrections.C:1988
GetCorrections.C:1989
GetCorrections.C:1990
GetCorrections.C:1991
GetCorrections.C:1992
GetCorrections.C:1993
GetCorrections.C:1994
GetCorrections.C:1995
GetCorrections.C:1996
GetCorrections.C:1997
GetCorrections.C:1998
GetCorrections.C:1999
GetCorrections.C:2000
GetCorrections.C:2001
GetCorrections.C:2002
GetCorrections.C:2003
GetCorrections.C:2004
GetCorrections.C:2005
GetCorrections.C:2006
GetCorrections.C:2007
GetCorrections.C:2008
GetCorrections.C:2009
GetCorrections.C:2010
GetCorrections.C:2011
GetCorrections.C:2012
GetCorrections.C:2013
GetCorrections.C:2014
GetCorrections.C:2015
GetCorrections.C:2016
GetCorrections.C:2017
GetCorrections.C:2018
GetCorrections.C:2019
GetCorrections.C:2020
GetCorrections.C:2021
GetCorrections.C:2022
GetCorrections.C:2023
GetCorrections.C:2024
GetCorrections.C:2025
GetCorrections.C:2026
GetCorrections.C:2027
GetCorrections.C:2028
GetCorrections.C:2029
GetCorrections.C:2030
GetCorrections.C:2031
GetCorrections.C:2032
GetCorrections.C:2033
GetCorrections.C:2034
GetCorrections.C:2035
GetCorrections.C:2036
GetCorrections.C:2037
GetCorrections.C:2038
GetCorrections.C:2039
GetCorrections.C:2040
GetCorrections.C:2041
GetCorrections.C:2042
GetCorrections.C:2043
GetCorrections.C:2044
GetCorrections.C:2045
GetCorrections.C:2046
GetCorrections.C:2047
GetCorrections.C:2048
GetCorrections.C:2049
GetCorrections.C:2050
GetCorrections.C:2051
GetCorrections.C:2052
GetCorrections.C:2053
GetCorrections.C:2054
GetCorrections.C:2055
GetCorrections.C:2056
GetCorrections.C:2057
GetCorrections.C:2058
GetCorrections.C:2059
GetCorrections.C:2060
GetCorrections.C:2061
GetCorrections.C:2062
GetCorrections.C:2063
GetCorrections.C:2064
GetCorrections.C:2065
GetCorrections.C:2066
GetCorrections.C:2067
GetCorrections.C:2068
GetCorrections.C:2069
GetCorrections.C:2070
GetCorrections.C:2071
GetCorrections.C:2072
GetCorrections.C:2073
GetCorrections.C:2074
GetCorrections.C:2075
GetCorrections.C:2076
GetCorrections.C:2077
GetCorrections.C:2078
GetCorrections.C:2079
GetCorrections.C:2080
GetCorrections.C:2081
GetCorrections.C:2082
GetCorrections.C:2083
GetCorrections.C:2084
GetCorrections.C:2085
GetCorrections.C:2086
GetCorrections.C:2087
GetCorrections.C:2088
GetCorrections.C:2089
GetCorrections.C:2090
GetCorrections.C:2091
GetCorrections.C:2092
GetCorrections.C:2093
GetCorrections.C:2094
GetCorrections.C:2095
GetCorrections.C:2096
GetCorrections.C:2097
GetCorrections.C:2098
GetCorrections.C:2099
GetCorrections.C:2100
GetCorrections.C:2101
GetCorrections.C:2102
GetCorrections.C:2103
GetCorrections.C:2104
GetCorrections.C:2105
GetCorrections.C:2106
GetCorrections.C:2107
GetCorrections.C:2108
GetCorrections.C:2109
GetCorrections.C:2110
GetCorrections.C:2111
GetCorrections.C:2112
GetCorrections.C:2113
GetCorrections.C:2114
GetCorrections.C:2115
GetCorrections.C:2116
GetCorrections.C:2117
GetCorrections.C:2118
GetCorrections.C:2119
GetCorrections.C:2120
GetCorrections.C:2121
GetCorrections.C:2122
GetCorrections.C:2123
GetCorrections.C:2124
GetCorrections.C:2125
GetCorrections.C:2126
GetCorrections.C:2127
GetCorrections.C:2128
GetCorrections.C:2129
GetCorrections.C:2130
GetCorrections.C:2131
GetCorrections.C:2132
GetCorrections.C:2133
GetCorrections.C:2134
GetCorrections.C:2135
GetCorrections.C:2136
GetCorrections.C:2137
GetCorrections.C:2138
GetCorrections.C:2139
GetCorrections.C:2140
GetCorrections.C:2141
GetCorrections.C:2142
GetCorrections.C:2143
GetCorrections.C:2144
GetCorrections.C:2145
GetCorrections.C:2146
GetCorrections.C:2147
GetCorrections.C:2148
GetCorrections.C:2149
GetCorrections.C:2150
GetCorrections.C:2151
GetCorrections.C:2152
GetCorrections.C:2153
GetCorrections.C:2154
GetCorrections.C:2155
GetCorrections.C:2156
GetCorrections.C:2157
GetCorrections.C:2158
GetCorrections.C:2159
GetCorrections.C:2160
GetCorrections.C:2161
GetCorrections.C:2162
GetCorrections.C:2163
GetCorrections.C:2164
GetCorrections.C:2165
GetCorrections.C:2166
GetCorrections.C:2167
GetCorrections.C:2168
GetCorrections.C:2169
GetCorrections.C:2170
GetCorrections.C:2171
GetCorrections.C:2172
GetCorrections.C:2173
GetCorrections.C:2174
GetCorrections.C:2175
GetCorrections.C:2176
GetCorrections.C:2177
GetCorrections.C:2178
GetCorrections.C:2179
GetCorrections.C:2180
GetCorrections.C:2181
GetCorrections.C:2182
GetCorrections.C:2183
GetCorrections.C:2184
GetCorrections.C:2185
GetCorrections.C:2186
GetCorrections.C:2187
GetCorrections.C:2188
GetCorrections.C:2189
GetCorrections.C:2190
GetCorrections.C:2191
GetCorrections.C:2192
GetCorrections.C:2193
GetCorrections.C:2194
GetCorrections.C:2195
GetCorrections.C:2196
GetCorrections.C:2197
GetCorrections.C:2198
GetCorrections.C:2199
GetCorrections.C:2200
GetCorrections.C:2201
GetCorrections.C:2202
GetCorrections.C:2203
GetCorrections.C:2204
GetCorrections.C:2205
GetCorrections.C:2206
GetCorrections.C:2207
GetCorrections.C:2208
GetCorrections.C:2209
GetCorrections.C:2210
GetCorrections.C:2211
GetCorrections.C:2212
GetCorrections.C:2213
GetCorrections.C:2214
GetCorrections.C:2215
GetCorrections.C:2216
GetCorrections.C:2217
GetCorrections.C:2218
GetCorrections.C:2219
GetCorrections.C:2220
GetCorrections.C:2221
GetCorrections.C:2222
GetCorrections.C:2223
GetCorrections.C:2224
GetCorrections.C:2225
GetCorrections.C:2226
GetCorrections.C:2227
GetCorrections.C:2228
GetCorrections.C:2229
GetCorrections.C:2230
GetCorrections.C:2231
GetCorrections.C:2232
GetCorrections.C:2233
GetCorrections.C:2234
GetCorrections.C:2235
GetCorrections.C:2236
GetCorrections.C:2237
GetCorrections.C:2238
GetCorrections.C:2239
GetCorrections.C:2240
GetCorrections.C:2241
GetCorrections.C:2242
GetCorrections.C:2243
GetCorrections.C:2244
GetCorrections.C:2245
GetCorrections.C:2246
GetCorrections.C:2247
GetCorrections.C:2248
GetCorrections.C:2249
GetCorrections.C:2250
GetCorrections.C:2251
GetCorrections.C:2252
GetCorrections.C:2253
GetCorrections.C:2254
GetCorrections.C:2255
GetCorrections.C:2256
GetCorrections.C:2257
GetCorrections.C:2258
GetCorrections.C:2259
GetCorrections.C:2260
GetCorrections.C:2261
GetCorrections.C:2262
GetCorrections.C:2263
GetCorrections.C:2264
GetCorrections.C:2265
GetCorrections.C:2266
GetCorrections.C:2267
GetCorrections.C:2268
GetCorrections.C:2269
GetCorrections.C:2270
GetCorrections.C:2271
GetCorrections.C:2272
GetCorrections.C:2273
GetCorrections.C:2274
GetCorrections.C:2275
GetCorrections.C:2276
GetCorrections.C:2277
GetCorrections.C:2278
GetCorrections.C:2279
GetCorrections.C:2280
GetCorrections.C:2281
GetCorrections.C:2282
GetCorrections.C:2283
GetCorrections.C:2284
GetCorrections.C:2285
GetCorrections.C:2286
GetCorrections.C:2287
GetCorrections.C:2288
GetCorrections.C:2289
GetCorrections.C:2290
GetCorrections.C:2291
GetCorrections.C:2292
GetCorrections.C:2293
GetCorrections.C:2294
GetCorrections.C:2295
GetCorrections.C:2296
GetCorrections.C:2297
GetCorrections.C:2298
GetCorrections.C:2299
GetCorrections.C:2300
GetCorrections.C:2301
GetCorrections.C:2302
GetCorrections.C:2303
GetCorrections.C:2304
GetCorrections.C:2305
GetCorrections.C:2306
GetCorrections.C:2307
GetCorrections.C:2308
GetCorrections.C:2309
GetCorrections.C:2310
GetCorrections.C:2311
GetCorrections.C:2312
GetCorrections.C:2313
GetCorrections.C:2314
GetCorrections.C:2315
GetCorrections.C:2316
GetCorrections.C:2317
GetCorrections.C:2318
GetCorrections.C:2319
GetCorrections.C:2320
GetCorrections.C:2321
GetCorrections.C:2322
GetCorrections.C:2323
GetCorrections.C:2324
GetCorrections.C:2325
GetCorrections.C:2326
GetCorrections.C:2327
GetCorrections.C:2328
GetCorrections.C:2329
GetCorrections.C:2330
GetCorrections.C:2331
GetCorrections.C:2332
GetCorrections.C:2333
GetCorrections.C:2334
GetCorrections.C:2335
GetCorrections.C:2336
GetCorrections.C:2337
GetCorrections.C:2338
GetCorrections.C:2339
GetCorrections.C:2340
GetCorrections.C:2341
GetCorrections.C:2342
GetCorrections.C:2343
GetCorrections.C:2344
GetCorrections.C:2345
GetCorrections.C:2346
GetCorrections.C:2347
GetCorrections.C:2348
GetCorrections.C:2349
GetCorrections.C:2350
GetCorrections.C:2351
GetCorrections.C:2352
GetCorrections.C:2353
GetCorrections.C:2354
GetCorrections.C:2355
GetCorrections.C:2356
GetCorrections.C:2357
GetCorrections.C:2358
GetCorrections.C:2359
GetCorrections.C:2360
GetCorrections.C:2361
GetCorrections.C:2362
GetCorrections.C:2363
GetCorrections.C:2364
GetCorrections.C:2365
GetCorrections.C:2366
GetCorrections.C:2367
GetCorrections.C:2368
GetCorrections.C:2369
GetCorrections.C:2370
GetCorrections.C:2371
GetCorrections.C:2372
GetCorrections.C:2373
GetCorrections.C:2374
GetCorrections.C:2375
GetCorrections.C:2376
GetCorrections.C:2377
GetCorrections.C:2378
GetCorrections.C:2379
GetCorrections.C:2380
GetCorrections.C:2381
GetCorrections.C:2382
GetCorrections.C:2383
GetCorrections.C:2384
GetCorrections.C:2385
GetCorrections.C:2386
GetCorrections.C:2387
GetCorrections.C:2388
GetCorrections.C:2389
GetCorrections.C:2390
GetCorrections.C:2391
GetCorrections.C:2392
GetCorrections.C:2393
GetCorrections.C:2394
GetCorrections.C:2395
GetCorrections.C:2396
GetCorrections.C:2397
GetCorrections.C:2398
GetCorrections.C:2399
GetCorrections.C:2400
GetCorrections.C:2401
GetCorrections.C:2402
GetCorrections.C:2403
GetCorrections.C:2404
GetCorrections.C:2405
GetCorrections.C:2406
GetCorrections.C:2407
GetCorrections.C:2408
GetCorrections.C:2409
GetCorrections.C:2410
GetCorrections.C:2411
GetCorrections.C:2412
GetCorrections.C:2413
GetCorrections.C:2414
GetCorrections.C:2415
GetCorrections.C:2416
GetCorrections.C:2417
GetCorrections.C:2418
GetCorrections.C:2419
GetCorrections.C:2420
GetCorrections.C:2421
GetCorrections.C:2422
GetCorrections.C:2423
GetCorrections.C:2424
GetCorrections.C:2425
GetCorrections.C:2426
GetCorrections.C:2427
GetCorrections.C:2428
GetCorrections.C:2429
GetCorrections.C:2430
GetCorrections.C:2431
GetCorrections.C:2432
GetCorrections.C:2433
GetCorrections.C:2434
GetCorrections.C:2435
GetCorrections.C:2436
GetCorrections.C:2437
GetCorrections.C:2438
GetCorrections.C:2439
GetCorrections.C:2440
GetCorrections.C:2441
GetCorrections.C:2442
GetCorrections.C:2443
GetCorrections.C:2444
GetCorrections.C:2445
GetCorrections.C:2446
GetCorrections.C:2447
GetCorrections.C:2448
GetCorrections.C:2449
GetCorrections.C:2450
GetCorrections.C:2451
GetCorrections.C:2452
GetCorrections.C:2453
GetCorrections.C:2454
GetCorrections.C:2455
GetCorrections.C:2456
GetCorrections.C:2457
GetCorrections.C:2458
GetCorrections.C:2459
GetCorrections.C:2460
GetCorrections.C:2461
GetCorrections.C:2462
GetCorrections.C:2463
GetCorrections.C:2464
GetCorrections.C:2465
GetCorrections.C:2466
GetCorrections.C:2467
GetCorrections.C:2468
GetCorrections.C:2469
GetCorrections.C:2470
GetCorrections.C:2471
GetCorrections.C:2472
GetCorrections.C:2473
GetCorrections.C:2474
GetCorrections.C:2475
GetCorrections.C:2476
GetCorrections.C:2477
GetCorrections.C:2478
GetCorrections.C:2479
GetCorrections.C:2480
GetCorrections.C:2481
GetCorrections.C:2482
GetCorrections.C:2483
GetCorrections.C:2484
GetCorrections.C:2485
GetCorrections.C:2486
GetCorrections.C:2487
GetCorrections.C:2488
GetCorrections.C:2489
GetCorrections.C:2490
GetCorrections.C:2491
GetCorrections.C:2492
GetCorrections.C:2493
GetCorrections.C:2494
GetCorrections.C:2495
GetCorrections.C:2496
GetCorrections.C:2497
GetCorrections.C:2498
GetCorrections.C:2499
GetCorrections.C:2500
GetCorrections.C:2501
GetCorrections.C:2502
GetCorrections.C:2503
GetCorrections.C:2504
GetCorrections.C:2505
GetCorrections.C:2506
GetCorrections.C:2507
GetCorrections.C:2508
GetCorrections.C:2509
GetCorrections.C:2510
GetCorrections.C:2511
GetCorrections.C:2512
GetCorrections.C:2513
GetCorrections.C:2514
GetCorrections.C:2515
GetCorrections.C:2516
GetCorrections.C:2517
GetCorrections.C:2518
GetCorrections.C:2519
GetCorrections.C:2520
GetCorrections.C:2521
GetCorrections.C:2522
GetCorrections.C:2523
GetCorrections.C:2524
GetCorrections.C:2525
GetCorrections.C:2526
GetCorrections.C:2527
GetCorrections.C:2528
GetCorrections.C:2529
GetCorrections.C:2530
GetCorrections.C:2531
GetCorrections.C:2532
GetCorrections.C:2533
GetCorrections.C:2534
GetCorrections.C:2535
GetCorrections.C:2536
GetCorrections.C:2537
GetCorrections.C:2538
GetCorrections.C:2539
GetCorrections.C:2540
GetCorrections.C:2541
GetCorrections.C:2542
GetCorrections.C:2543
GetCorrections.C:2544
GetCorrections.C:2545
GetCorrections.C:2546
GetCorrections.C:2547
GetCorrections.C:2548
GetCorrections.C:2549
GetCorrections.C:2550
GetCorrections.C:2551
GetCorrections.C:2552
GetCorrections.C:2553
GetCorrections.C:2554
GetCorrections.C:2555
GetCorrections.C:2556
GetCorrections.C:2557
GetCorrections.C:2558
GetCorrections.C:2559
GetCorrections.C:2560
GetCorrections.C:2561
GetCorrections.C:2562
GetCorrections.C:2563
GetCorrections.C:2564
GetCorrections.C:2565
GetCorrections.C:2566
GetCorrections.C:2567
GetCorrections.C:2568
GetCorrections.C:2569
GetCorrections.C:2570
GetCorrections.C:2571
GetCorrections.C:2572
GetCorrections.C:2573
GetCorrections.C:2574
GetCorrections.C:2575
GetCorrections.C:2576
GetCorrections.C:2577
GetCorrections.C:2578
GetCorrections.C:2579
GetCorrections.C:2580
GetCorrections.C:2581
GetCorrections.C:2582
GetCorrections.C:2583
GetCorrections.C:2584
GetCorrections.C:2585
GetCorrections.C:2586
GetCorrections.C:2587
GetCorrections.C:2588
GetCorrections.C:2589
GetCorrections.C:2590
GetCorrections.C:2591
GetCorrections.C:2592
GetCorrections.C:2593
GetCorrections.C:2594
GetCorrections.C:2595
GetCorrections.C:2596
GetCorrections.C:2597
GetCorrections.C:2598
GetCorrections.C:2599
GetCorrections.C:2600
GetCorrections.C:2601
GetCorrections.C:2602
GetCorrections.C:2603
GetCorrections.C:2604
GetCorrections.C:2605
GetCorrections.C:2606
GetCorrections.C:2607
GetCorrections.C:2608
GetCorrections.C:2609
GetCorrections.C:2610
GetCorrections.C:2611
GetCorrections.C:2612
GetCorrections.C:2613
GetCorrections.C:2614
GetCorrections.C:2615
GetCorrections.C:2616
GetCorrections.C:2617
GetCorrections.C:2618
GetCorrections.C:2619
GetCorrections.C:2620
GetCorrections.C:2621
GetCorrections.C:2622
GetCorrections.C:2623
GetCorrections.C:2624
GetCorrections.C:2625
GetCorrections.C:2626
GetCorrections.C:2627
GetCorrections.C:2628
GetCorrections.C:2629
GetCorrections.C:2630
GetCorrections.C:2631
GetCorrections.C:2632
GetCorrections.C:2633
GetCorrections.C:2634
GetCorrections.C:2635
GetCorrections.C:2636
GetCorrections.C:2637
GetCorrections.C:2638
GetCorrections.C:2639
GetCorrections.C:2640
GetCorrections.C:2641
GetCorrections.C:2642
GetCorrections.C:2643
GetCorrections.C:2644
GetCorrections.C:2645
GetCorrections.C:2646
GetCorrections.C:2647
GetCorrections.C:2648
GetCorrections.C:2649
GetCorrections.C:2650
GetCorrections.C:2651
GetCorrections.C:2652
GetCorrections.C:2653
GetCorrections.C:2654
GetCorrections.C:2655
GetCorrections.C:2656
GetCorrections.C:2657
GetCorrections.C:2658
GetCorrections.C:2659
GetCorrections.C:2660
GetCorrections.C:2661
GetCorrections.C:2662
GetCorrections.C:2663
GetCorrections.C:2664
GetCorrections.C:2665
GetCorrections.C:2666
GetCorrections.C:2667
GetCorrections.C:2668
GetCorrections.C:2669
GetCorrections.C:2670
GetCorrections.C:2671
GetCorrections.C:2672
GetCorrections.C:2673
GetCorrections.C:2674
GetCorrections.C:2675
GetCorrections.C:2676
GetCorrections.C:2677
GetCorrections.C:2678
GetCorrections.C:2679
GetCorrections.C:2680
GetCorrections.C:2681
GetCorrections.C:2682
GetCorrections.C:2683
GetCorrections.C:2684
GetCorrections.C:2685
GetCorrections.C:2686
GetCorrections.C:2687
GetCorrections.C:2688
GetCorrections.C:2689
GetCorrections.C:2690
GetCorrections.C:2691
GetCorrections.C:2692
GetCorrections.C:2693
GetCorrections.C:2694
GetCorrections.C:2695
GetCorrections.C:2696
GetCorrections.C:2697
GetCorrections.C:2698
GetCorrections.C:2699
GetCorrections.C:2700
GetCorrections.C:2701
GetCorrections.C:2702
GetCorrections.C:2703
GetCorrections.C:2704
GetCorrections.C:2705
GetCorrections.C:2706
GetCorrections.C:2707
GetCorrections.C:2708
GetCorrections.C:2709
GetCorrections.C:2710
GetCorrections.C:2711
GetCorrections.C:2712
GetCorrections.C:2713
GetCorrections.C:2714
GetCorrections.C:2715
GetCorrections.C:2716
GetCorrections.C:2717
GetCorrections.C:2718
GetCorrections.C:2719
GetCorrections.C:2720
GetCorrections.C:2721
GetCorrections.C:2722
GetCorrections.C:2723
GetCorrections.C:2724
GetCorrections.C:2725
GetCorrections.C:2726
GetCorrections.C:2727
GetCorrections.C:2728
GetCorrections.C:2729
GetCorrections.C:2730
GetCorrections.C:2731
GetCorrections.C:2732
GetCorrections.C:2733
GetCorrections.C:2734
GetCorrections.C:2735
GetCorrections.C:2736
GetCorrections.C:2737
GetCorrections.C:2738
GetCorrections.C:2739
GetCorrections.C:2740
GetCorrections.C:2741
GetCorrections.C:2742
GetCorrections.C:2743
GetCorrections.C:2744
GetCorrections.C:2745
GetCorrections.C:2746
GetCorrections.C:2747
GetCorrections.C:2748
GetCorrections.C:2749
GetCorrections.C:2750
GetCorrections.C:2751
GetCorrections.C:2752
GetCorrections.C:2753
GetCorrections.C:2754
GetCorrections.C:2755
GetCorrections.C:2756
GetCorrections.C:2757
GetCorrections.C:2758
GetCorrections.C:2759
GetCorrections.C:2760
GetCorrections.C:2761
GetCorrections.C:2762
GetCorrections.C:2763
GetCorrections.C:2764
GetCorrections.C:2765
GetCorrections.C:2766
GetCorrections.C:2767
GetCorrections.C:2768
GetCorrections.C:2769
GetCorrections.C:2770
GetCorrections.C:2771
GetCorrections.C:2772
GetCorrections.C:2773
GetCorrections.C:2774
GetCorrections.C:2775
GetCorrections.C:2776
GetCorrections.C:2777
GetCorrections.C:2778
GetCorrections.C:2779
GetCorrections.C:2780
GetCorrections.C:2781
GetCorrections.C:2782
GetCorrections.C:2783
GetCorrections.C:2784
GetCorrections.C:2785
GetCorrections.C:2786
GetCorrections.C:2787
GetCorrections.C:2788
GetCorrections.C:2789
GetCorrections.C:2790
GetCorrections.C:2791
GetCorrections.C:2792
GetCorrections.C:2793
GetCorrections.C:2794
GetCorrections.C:2795
GetCorrections.C:2796
GetCorrections.C:2797
GetCorrections.C:2798
GetCorrections.C:2799
GetCorrections.C:2800
GetCorrections.C:2801
GetCorrections.C:2802
GetCorrections.C:2803
GetCorrections.C:2804
GetCorrections.C:2805
GetCorrections.C:2806
GetCorrections.C:2807
GetCorrections.C:2808
GetCorrections.C:2809
GetCorrections.C:2810
GetCorrections.C:2811
GetCorrections.C:2812
GetCorrections.C:2813
GetCorrections.C:2814
GetCorrections.C:2815
GetCorrections.C:2816
GetCorrections.C:2817
GetCorrections.C:2818
GetCorrections.C:2819
GetCorrections.C:2820
GetCorrections.C:2821
GetCorrections.C:2822
GetCorrections.C:2823
GetCorrections.C:2824
GetCorrections.C:2825
GetCorrections.C:2826
GetCorrections.C:2827
GetCorrections.C:2828
GetCorrections.C:2829
GetCorrections.C:2830
GetCorrections.C:2831
GetCorrections.C:2832
GetCorrections.C:2833
GetCorrections.C:2834
GetCorrections.C:2835
GetCorrections.C:2836
GetCorrections.C:2837
GetCorrections.C:2838
GetCorrections.C:2839
GetCorrections.C:2840
GetCorrections.C:2841
GetCorrections.C:2842
GetCorrections.C:2843
GetCorrections.C:2844
GetCorrections.C:2845
GetCorrections.C:2846
GetCorrections.C:2847
GetCorrections.C:2848
GetCorrections.C:2849
GetCorrections.C:2850
GetCorrections.C:2851
GetCorrections.C:2852
GetCorrections.C:2853
GetCorrections.C:2854
GetCorrections.C:2855
GetCorrections.C:2856
GetCorrections.C:2857
GetCorrections.C:2858
GetCorrections.C:2859
GetCorrections.C:2860
GetCorrections.C:2861
GetCorrections.C:2862
GetCorrections.C:2863
GetCorrections.C:2864
GetCorrections.C:2865
GetCorrections.C:2866
GetCorrections.C:2867
GetCorrections.C:2868
GetCorrections.C:2869
GetCorrections.C:2870
GetCorrections.C:2871
GetCorrections.C:2872
GetCorrections.C:2873
GetCorrections.C:2874
GetCorrections.C:2875
GetCorrections.C:2876
GetCorrections.C:2877
GetCorrections.C:2878
GetCorrections.C:2879
GetCorrections.C:2880
GetCorrections.C:2881
GetCorrections.C:2882
GetCorrections.C:2883
GetCorrections.C:2884
GetCorrections.C:2885
GetCorrections.C:2886
GetCorrections.C:2887
GetCorrections.C:2888
GetCorrections.C:2889
GetCorrections.C:2890
GetCorrections.C:2891
GetCorrections.C:2892