18 #if !defined(__CINT__) || defined(__MAKECINT__) 21 #include <TDirectoryFile.h> 39 #include <TGraphErrors.h> 55 TString titleData=
"LHC11cd_EMC7",
66 TString daLeg =
"pp@7 TeV, LHC11c+d EMC7";
67 if(titleData ==
"LHC11cd_EMC7")
68 daLeg =
"pp@7 TeV, LHC11c+d EMC7";
69 if(titleData ==
"LHC11cd_INT7")
70 daLeg =
"pp@7 TeV, LHC11c+d INT7";
73 TString mcLeg =
"MC: #gammaJ+JJ({p^{EMCal}_{T,#gamma}>3.5 GeV/#it{c})";
74 if ( titleMC.Contains(
"Low") )
77 Double_t energyBins[] = { 8, 10, 12, 14, 16, 18, 20 };
78 for(
Int_t ie = 0; ie < nEBins+1; ie++) binE[ie] = energyBins[ie];
80 mcLeg =
"#gammaJ+JJ(p^{EMCal}_{T,#gamma}>3.5 GeV/#it{c})";
83 if ( titleMC.Contains(
"High") )
86 Double_t energyBins[] = { 18, 20, 25, 30, 35, 40, 50 };
87 for(
Int_t ie = 0; ie < nEBins+1; ie++) binE[ie] = energyBins[ie];
88 mcDir =
"JJ_Dec_High_GJ";
89 mcLeg =
"#gammaJ+JJ(p^{EMCal}_{T,#gamma}>7 GeV/#it{c})}";
92 if ( titleMC.Contains(
"MB") )
96 for(
Int_t ie = 0; ie < nEBins+1; ie++) binE[ie] = energyBins[ie];
103 Bool_t bIndFracCase = kFALSE;
104 if ( bIndFracCase ) nprodtmp = 7;
110 if ( debug ) printf(
"N prod %d\n",nProd);
115 prod [nproditer++] = Form(
"data/module/TCardChannel3/%s",titleData.Data());
116 prodLeg[nproditer-1] =
"Data";
117 prod [nproditer++] = Form(
"simu/module//pp_7TeV_%s/TCardChannel_Mimic0_Scaled2_v2",mcDir.Data());
118 prodLeg[nproditer-1] =
"MC, default";
119 prod [nproditer++] = Form(
"simu/module//pp_7TeV_%s/TCardChannel_Mimic10c_EcellCut_Scaled2_v2",mcDir.Data());
120 prodLeg[nproditer-1] =
"MC, mimic";
154 printf(
"INDUCED Fraction cases!!!! %d\n",nproditer);
155 if ( titleMC.Contains(
"Low") )
157 prod [nproditer++] =
"data/LHC11cd_EMC7";
158 prodLeg[nproditer-1] =
"Data";
159 prod [nproditer++] =
"simu/mimic/pp_7TeV_JJ_Dec_GJ/default/sumw2on";
160 prodLeg[nproditer-1] =
"MC, Default";
169 prod [nproditer++] =
"simu/mimic/pp_7TeV_JJ_Dec_GJ/TCard_C_70/sumw2on";
170 prodLeg[nproditer-1] =
"MC, mimic #mu_{1}=0.5%";
179 prod [nproditer++] =
"simu/mimic/pp_7TeV_JJ_Dec_GJ/TCard_C_68/sumw2on";
180 prodLeg[nproditer-1] =
"MC, mimic #mu_{1}=1.0%";
195 prod [nproditer++] =
"simu/mimic/pp_7TeV_JJ_Dec_GJ/TCard_C_86/sumw2on";
196 prodLeg[nproditer-1] =
"MC, mimic #mu_{1}=1.25%";
207 prod [nproditer++] =
"simu/mimic/pp_7TeV_JJ_Dec_GJ/TCard_C_65/sumw2on";
208 prodLeg[nproditer-1] =
"MC, mimic #mu_{1}=1.4%";
211 prod [nproditer++] =
"simu/mimic/pp_7TeV_JJ_Dec_GJ/TCard_C_64/sumw2on";
212 prodLeg[nproditer-1] =
"MC, mimic #mu_{1}=1.5%";
222 printf(
"iprod %d,\n \t file: %s,\n \t legend: %s\n",
223 iprod,prod[iprod].
Data(),prodLeg[iprod].
Data());
229 "SMM02NoCut",
"SMM02",
"SMM20LowM02NoCut",
"SMM20LowM02",
"SMM20HighM02NoCut",
"SMM20HighM02",
"SMNCell",
230 "SMEMaxEClusterRat",
"SMNCellModuleMax",
"SMNCellModuleOut",
231 "SMECellModuleMax",
"SMECellModuleOut",
"SMNCellModuleMaxOutRat",
"SMECellModuleMaxRat",
232 "SMECellModuleMaxOutRat",
"SMECellModuleMaxTot",
"SMECellModuleMaxTotRat",
"SMECellModuleMaxTotRatClus",
233 "SMNCellModuleOutModDiff",
"SMNCellModuleOutModSame",
"SMECellModuleOutModDiff",
"SMECellModuleOutModSame",
234 "TCardChannelNCellModMax",
"TCardChannelNCell",
"TCardChannelM02",
"TCardChannelM02NoCut" 240 Bool_t plotRatio = kFALSE;
241 Bool_t saveHisto = kFALSE;
249 if ( clusterization ==
"" ) clusterization =
"_V1_Ecell100_Eseed500";
250 opt =
"_IndFracCases";
253 for(
Int_t ihisto = 0; ihisto < 26; ihisto++)
255 if(histoName[ihisto].Contains(
"TCard")) lastP = 15;
256 printf(
"histogram %d %s\n",ihisto,histoName[ihisto].
Data());
260 histoName[ihisto], nProd, prod, prodLeg,
261 clusterization,tm,pid,
266 plotRatio, saveHisto,
Compare a cluster parameter (shower shape) between different productions.
void MakeDataMCComparisonPerSMClusterEbin(TString titleMC="JJDecLow", TString titleData="LHC11cd_EMC7", Bool_t debug=kFALSE)
Bool_t Data(TH1F *h, Double_t *rangefit, Bool_t writefit, Double_t &sgn, Double_t &errsgn, Double_t &bkg, Double_t &errbkg, Double_t &sgnf, Double_t &errsgnf, Double_t &sigmafit, Int_t &status)
void CompareTH3DataAndMCProd(TString &outputFileName, TString histoName="SMM02NoCut", const Int_t nProd=3, TString *prod=0x0, TString *prodLeg=0x0, TString clusterization="", TString tm="_TMDep", TString pid="_Neutral", TString titleMC="JJDecLow", TString mcLeg="MC: #gammaJ+JJ({p^{EMCal}_{T,#gamma}>3.5 GeV/#it{c})", TString titleData="LHC11cd_EMC7", TString daLeg="pp@7 TeV, LHC11c+d EMC7", Int_t firstP=0, Int_t lastP=9, TArrayD binE=0, Int_t firstMC=1, Bool_t plotRatio=kFALSE, Bool_t saveHisto=kFALSE, TString opt="", Bool_t debug=kFALSE)
TString prod[]
productions to be compared, directory name