Warning, /Validation/RecoEgamma/test/fastSimValidation.csh is written in an unsupported language. File is not indexed.
0001 ##!/bin/csh -x
0002 #!/bin/csh
0003
0004 #This script can be used to generate a web page to compare histograms from
0005 #two input root files produced using the EDAnalyzers in RecoEgamma/Examples,
0006 #by running one of:
0007 #
0008 #
0009 #
0010 # "Validation/RecoEgamma/test/PhotonValidator_cfg.py
0011 #
0012 # The default list of histograms (configurable) is based on version VXX-XX-XX
0013 # of Validation/RecoEgamma
0014 #
0015 #Two files are created by this script: validation.C and validation.html.
0016 #validation.C should be run inside root to greate a set of gif images
0017 #which can then be viewed in a web browser using validation.html.
0018
0019 #=============BEGIN CONFIGURATION=================
0020 setenv TYPE GEDPhotons
0021 setenv ANALYZERNAME pfPhotonValidator
0022 #setenv TYPE Photons
0023 #setenv ANALYZERNAME PhotonValidator
0024 setenv CMSSWver1 7_0_0
0025 setenv RELEASE 7_0_0
0026 setenv PRERELEASE
0027 setenv UPGRADE True
0028 setenv LHCENERGY 13
0029
0030 setenv FULLGLOBALTAG POSTLS170_V3-v2
0031 setenv FASTGLOBALTAG POSTLS170_V3_FastSim-v2
0032
0033 #setenv RELEASE ${RELEASE}_${PRERELEASE}
0034 setenv RELEASE ${RELEASE}
0035
0036
0037 setenv WorkDir1 /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver1}/src/Validation/RecoEgamma/test
0038 #setenv WorkDir1 /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver1}_${PRERELEASE}/src/Validation/RecoEgamma/test
0039
0040
0041 #Name of sample (affects output directory name and htmldescription only)
0042
0043
0044 setenv SAMPLE H130GGgluonfusion
0045
0046 #TYPE must be one ofPixelMatchGsfElectron, Photon
0047
0048 #==============END BASIC CONFIGURATION==================
0049
0050
0051 #Input root trees for the two cases to be compared
0052
0053
0054 if ($SAMPLE == SingleGammaPt10IDEAL) then
0055
0056 setenv FULLSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_SingleGammaPt10.root
0057 setenv FASTSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_SingleGammaPt10_FastSim.root
0058
0059 else if ($SAMPLE == SingleGammaPt35IDEAL) then
0060
0061 setenv FULLSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_SingleGammaPt35.root
0062 setenv FASTSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_SingleGammaPt35_FastSim.root
0063
0064
0065 else if ($SAMPLE == SingleGammaFlatPt10_100IDEAL) then
0066
0067 setenv FULLSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_SingleGammaFlatPt10To100.root
0068 setenv FASTSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_SingleGammaFlatPt10To100_FastSim.root
0069
0070
0071 else if ($SAMPLE == H130GGgluonfusion) then
0072
0073 setenv HISTOPATHNAME_Efficiencies DQMData/Run\ 1/EgammaV/Run\ summary/${ANALYZERNAME}/Efficiencies
0074 setenv HISTOPATHNAME_Photons DQMData/Run\ 1/EgammaV/Run\ summary/${ANALYZERNAME}/Photons
0075 setenv HISTOPATHNAME_Conversions DQMData/Run\ 1/EgammaV/Run\ summary/${ANALYZERNAME}/ConversionInfo
0076
0077 if ( $UPGRADE == True ) then
0078 setenv FULLSIM ${WorkDir1}/DQM_V0001_R000000001__RelValH130GGgluonfusion_${LHCENERGY}__CMSSW_${RELEASE}-${FULLGLOBALTAG}__DQM.root
0079 setenv FASTSIM ${WorkDir1}/DQM_V0001_R000000001__RelValH130GGgluonfusion_${LHCENERGY}__CMSSW_${RELEASE}-${FASTGLOBALTAG}__DQM.root
0080 endif
0081
0082
0083 else if ($SAMPLE == PhotonJets_Pt_10) then
0084
0085 setenv FULLSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_PhotonJets_Pt_10.root
0086 setenv FASTSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_PhotonJets_Pt_10_FastSim.root
0087
0088
0089 else if ($SAMPLE == GammaJets_Pt_80_120STARTUP) then
0090
0091 setenv FULLSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_GammaJets_Pt_80_120.root
0092 setenv FASTSIM ${WorkDir1}/PhotonValidationRelVal${RELEASE}_GammaJets_Pt_80_120.root
0093
0094
0095 else if ($SAMPLE == QCD_Pt_80_120STARTUP) then
0096
0097 endif
0098
0099
0100
0101 #Location of output. The default will put your output in:
0102 #http://cmsdoc.cern.ch/Physics/egamma/www/validation/
0103
0104 setenv CURRENTDIR $PWD
0105 setenv OUTPATH /afs/cern.ch/cms/Physics/egamma/www/validation/Photons
0106 cd $OUTPATH
0107 if (! -d $RELEASE) then
0108 mkdir $RELEASE
0109 endif
0110 setenv OUTPATH $OUTPATH/$RELEASE
0111 cd $OUTPATH
0112
0113 if (! -d ${TYPE}) then
0114 mkdir ${TYPE}
0115 endif
0116 setenv OUTPATH $OUTPATH/${TYPE}
0117 cd $OUTPATH
0118
0119 if (! -d vsFastSim) then
0120 mkdir vsFastSim
0121 endif
0122 setenv OUTPATH $OUTPATH/vsFastSim
0123
0124 if ( $UPGRADE == True) then
0125 setenv OUTDIR $OUTPATH/${SAMPLE}_${LHCENERGY}TeV
0126 else
0127 setenv OUTDIR $OUTPATH/${SAMPLE}
0128 endif
0129
0130 if (! -d $OUTDIR) then
0131 cd $OUTPATH
0132 mkdir $OUTDIR
0133 cd $OUTDIR
0134 mkdir gifs
0135 endif
0136 cd $OUTDIR
0137
0138
0139 #The list of histograms to be compared for each TYPE can be configured below:
0140
0141
0142 if ( $TYPE == Photons || $TYPE == GEDPhotons ) then
0143
0144
0145 #cat > efficiencyForPhotons <<EOF
0146 # recoEffVsEta
0147 # recoEffVsPhi
0148 # recoEffVsEt
0149 # deadChVsEta
0150 # deadChVsPhi
0151 # deadChVsEt
0152 #EOF
0153
0154
0155 # gamgamMassAll
0156 # gamgamMassBarrel
0157 # gamgamMassEndcap
0158 # gamgamMassNoConvAll
0159 # gamgamMassNoConvBarrel
0160 # gamgamMassNoConvEndcap
0161 # gamgamMassConvAll
0162 # gamgamMassConvBarrel
0163 # gamgamMassConvEndcap
0164
0165
0166 cat > scaledhistosForPhotons <<EOF
0167
0168 scEta
0169 scPhi
0170 scEAll
0171 scEtAll
0172 phoEta
0173 phoPhi
0174 phoDEta
0175 phoDPhi
0176 phoEAll
0177 phoEtAll
0178 eResAll
0179 eResBarrel
0180 eResEndcap
0181 eResunconvAll
0182 eResunconvBarrel
0183 eResunconvEndcap
0184 eResconvAll
0185 eResconvBarrel
0186 eResconvEndcap
0187 psEEndcap
0188 isoTrkSolidConeDR04All
0189 isoTrkSolidConeDR04Barrel
0190 isoTrkSolidConeDR04Endcap
0191 nTrkSolidConeDR04All
0192 nTrkSolidConeDR04Barrel
0193 nTrkSolidConeDR04Endcap
0194 r9Barrel
0195 r9Endcap
0196 r1Barrel
0197 r1Endcap
0198 r2Barrel
0199 r2Endcap
0200 sigmaIetaIetaBarrel
0201 sigmaIetaIetaEndcap
0202 hOverEAll
0203 hOverEBarrel
0204 hOverEEndcap
0205 newhOverEAll
0206 newhOverEBarrel
0207 newhOverEEndcap
0208 hcalTowerSumEtConeDR04Barrel
0209 hcalTowerSumEtConeDR04Endcap
0210 hcalTowerBcSumEtConeDR04Barrel
0211 hcalTowerBcSumEtConeDR04Endcap
0212 ecalRecHitSumEtConeDR04Barrel
0213 ecalRecHitSumEtConeDR04Endcap
0214
0215
0216
0217
0218 EOF
0219
0220 cat > scaledhistosForPhotonsLogScale <<EOF
0221 hOverEAll
0222 hOverEBarrel
0223 hOverEEndcap
0224 newhOverEAll
0225 newhOverEBarrel
0226 newhOverEEndcap
0227 hcalTowerSumEtConeDR04Barrel
0228 hcalTowerSumEtConeDR04Endcap
0229 hcalTowerBcSumEtConeDR04Barrel
0230 hcalTowerBcSumEtConeDR04Endcap
0231 ecalRecHitSumEtConeDR04Barrel
0232 ecalRecHitSumEtConeDR04Endcap
0233 r9Barrel
0234 r9Endcap
0235 r1Barrel
0236 r1Endcap
0237 r2Barrel
0238 r2Endcap
0239 sigmaIetaIetaAll
0240 sigmaIetaIetaBarrel
0241 sigmaIetaIetaEndcap
0242
0243
0244
0245 EOF
0246
0247
0248 cat > unscaledhistosForPhotons <<EOF
0249 pEResVsR9Barrel
0250 pEResVsR9Endcap
0251 scpEResVsR9Barrel
0252 scpEResVsR9Endcap
0253 pEResVsEtAll
0254 pEResVsEtBarrel
0255 pEResVsEtEndcap
0256 pEResVsEtaAll
0257 pEResVsEtaUnconv
0258 pEResVsEtaConv
0259 pEcalRecHitSumEtConeDR04VsEtaAll
0260 pEcalRecHitSumEtConeDR04VsEtBarrel
0261 pEcalRecHitSumEtConeDR04VsEtEndcap
0262 pHcalTowerSumEtConeDR04VsEtaAll
0263 pHcalTowerSumEtConeDR04VsEtBarrel
0264 pHcalTowerSumEtConeDR04VsEtEndcap
0265 pHcalTowerBcSumEtConeDR04VsEtaAll
0266 pHcalTowerBcSumEtConeDR04VsEtBarrel
0267 pHcalTowerBcSumEtConeDR04VsEtEndcap
0268 pHoverEVsEtaAll
0269 pHoverEVsEtAll
0270 pnewHoverEVsEtaAll
0271 pnewHoverEVsEtAll
0272
0273
0274 EOF
0275
0276
0277
0278
0279
0280
0281 endif
0282
0283
0284
0285
0286
0287
0288 #=================END CONFIGURATION=====================
0289
0290 if (-e validation.C) rm validation.C
0291 touch validation.C
0292 cat > begin.C <<EOF
0293 {
0294 TFile *file_old = TFile::Open("$FASTSIM");
0295 TFile *file_new = TFile::Open("$FULLSIM");
0296
0297 EOF
0298 cat begin.C >>& validation.C
0299 rm begin.C
0300
0301 setenv N 1
0302
0303
0304
0305
0306
0307 foreach i (`cat efficiencyForPhotons`)
0308 cat > temp$N.C <<EOF
0309
0310 TCanvas *c$i = new TCanvas("c$i");
0311 c$i->SetFillColor(10);
0312 file_old->cd("$HISTOPATHNAME_Efficiencies");
0313 $i->SetStats(0);
0314 if ( $i==deadChVsEta || $i==deadChVsPhi || $i==deadChVsEt ) {
0315 $i->SetMinimum(0.);
0316 $i->SetMaximum(0.2);
0317 } else if ( $i==recoEffVsEt ) {
0318 $i->GetXaxis()->SetRangeUser(0.,200.);
0319 } else {
0320 $i->SetMinimum(0.);
0321 $i->SetMaximum(1.1);
0322 }
0323 $i->SetLineColor(kPink+8);
0324 $i->SetMarkerColor(kPink+8);
0325 $i->SetMarkerStyle(20);
0326 $i->SetMarkerSize(0.8);
0327 $i->SetLineWidth(1);
0328 $i->Draw();
0329 file_new->cd("$HISTOPATHNAME_Efficiencies");
0330 $i->SetStats(0);
0331 $i->SetMinimum(0.);
0332 $i->SetMaximum(1.1);
0333 $i->SetLineColor(kBlack);
0334 $i->SetMarkerColor(kBlack);
0335 $i->SetMarkerStyle(20);
0336 $i->SetMarkerSize(0.8);
0337 $i->SetLineWidth(1);
0338 $i->Draw("same");
0339 c$i->SaveAs("gifs/$i.gif");
0340
0341 EOF
0342 setenv N `expr $N + 1`
0343 end
0344
0345
0346 foreach i (`cat scaledhistosForPhotons`)
0347 cat > temp$N.C <<EOF
0348 TCanvas *c$i = new TCanvas("c$i");
0349 c$i->SetFillColor(10);
0350 c$i->Divide(1,2);
0351 c$i->cd(1);
0352 //file_new->cd("DQMData/EgammaV/PhotonValidator/Photons");
0353 file_new->cd("$HISTOPATHNAME_Photons");
0354 int nBins = $i->GetNbinsX();
0355 float xMin=$i->GetBinLowEdge(1);
0356 float xMax=$i->GetBinLowEdge(nBins)+$i->GetBinWidth(nBins);
0357 Double_t mnew=$i->GetMaximum();
0358 Double_t nnew=$i->GetEntries();
0359 //file_old->cd("DQMData/EgammaV/PhotonValidator/Photons");
0360 file_old->cd("$HISTOPATHNAME_Photons");
0361
0362 TH1F* hold=new TH1F("hold"," ",nBins,xMin,xMax);
0363 hold=$i;
0364 Double_t mold=$i->GetMaximum();
0365 Double_t nold=$i->GetEntries();
0366 if ( $i==scEAll || $i==phoEAll ) {
0367 $i->GetYaxis()->SetRangeUser(0.,6000.);
0368 } else if ( $i==eResAll || $i==eResunconvAll ) {
0369 $i->SetMaximum(12000.);
0370 } else if ( $i==eResEndcap ) {
0371 $i->SetMaximum(2500.);
0372 } else if ( $i==eResBarrel || $i==eResunconvBarrel ) {
0373 $i->SetMaximum(10000.);
0374 } else if ( $i==eResunconvEndcap ) {
0375 $i->SetMaximum(1500.);
0376 }
0377 $i->SetStats(0);
0378 $i->SetMinimum(0.);
0379 //if ( mnew > mold)
0380 // $i->SetMaximum(mnew+mnew*0.2);
0381 //else
0382 //$i->SetMaximum(mold+mold*0.2);
0383 //$i->SetMaximum(mold+mold*0.2);
0384 $i->SetLineColor(kPink+8);
0385 $i->SetFillColor(kPink+8);
0386 //$i->SetLineWidth(3);
0387 $i->Draw();
0388 //file_new->cd("DQMData/EgammaV/PhotonValidator/Photons");
0389 file_new->cd("$HISTOPATHNAME_Photons");
0390 Double_t nnew=$i->GetEntries();
0391 $i->SetStats(0);
0392 $i->SetLineColor(kBlack);
0393 $i->SetMarkerColor(kBlack);
0394 $i->SetMarkerStyle(20);
0395 $i->SetMarkerSize(0.8);
0396 //$i->SetLineWidth(1);
0397 $i->Scale(nold/nnew);
0398 TH1F* hnew=new TH1F("hnew"," ",nBins,xMin,xMax);
0399 hnew=$i;
0400 $i->Draw("e1same");
0401 c$i->cd(2);
0402 TH1F* ratio=new TH1F("ratio"," ",nBins,xMin,xMax);
0403 ratio->Divide(hnew,hold);
0404 for ( int i=1; i<=ratio->GetNbinsX(); i++ ) {
0405 float num=hnew->GetBinContent(i);
0406 float den=hold->GetBinContent(i);
0407 float dNum=hnew->GetBinError(i);
0408 float dDen=hold->GetBinError(i);
0409 float erro=0;
0410 if ( num!=0 && den!=0) {
0411 erro= ((1./den)*(1./den)*dNum*dNum) + ((num*num)/(den*den*den*den) * (dDen*dDen));
0412 erro=sqrt(erro);
0413 }
0414 ratio->SetBinError(i, erro);
0415 }
0416 ratio->SetStats(0);
0417 ratio->SetLineColor(1);
0418 ratio->SetLineWidth(2);
0419 if ( $i==psEEndcap) {
0420 ratio->GetYaxis()->SetRangeUser(0.,5.);
0421 }
0422 ratio->Draw("e");
0423 TLine *l = new TLine(xMin,1.,xMax,1.);
0424 l->Draw();
0425 c$i->SaveAs("gifs/$i.gif");
0426 //TString gifName=TString("gifs/$i")+"_ratio.gif";
0427 //c$i->SaveAs(gifName);
0428 EOF
0429 setenv N `expr $N + 1`
0430 end
0431
0432
0433
0434
0435
0436 foreach i (`cat scaledhistosForPhotonsLogScale`)
0437 cat > temp$N.C <<EOF
0438 TCanvas *cc$i = new TCanvas("cc$i");
0439 cc$i->cd();
0440 cc$i->SetFillColor(10);
0441 cc$i->SetLogy();
0442 file_new->cd("$HISTOPATHNAME_Photons");
0443 Double_t nnew=$i->GetEntries();
0444 file_old->cd("$HISTOPATHNAME_Photons");
0445 if ( $i==hcalTowerSumEtConeDR04Barrel || $i==hcalTowerSumEtConeDR04Endcap ) {
0446 $i->GetXaxis()->SetRangeUser(0.,10.);
0447 }
0448 Double_t nold=$i->GetEntries();
0449 $i->SetStats(0);
0450 $i->SetMinimum(1);
0451 $i->SetLineColor(kPink+8);
0452 $i->SetFillColor(kPink+8);
0453 $i->Draw();
0454 file_new->cd("$HISTOPATHNAME_Photons");
0455 Double_t nnew=$i->GetEntries();
0456 $i->SetStats(0);
0457 $i->SetLineColor(kBlack);
0458 $i->SetMarkerColor(kBlack);
0459 $i->SetMarkerStyle(20);
0460 $i->SetMarkerSize(0.8);
0461 $i->Draw("e1same");
0462 cc$i->SaveAs("gifs/log$i.gif");
0463
0464 EOF
0465 setenv N `expr $N + 1`
0466 end
0467
0468
0469
0470
0471
0472 foreach i (`cat unscaledhistosForPhotons`)
0473 cat > temp$N.C <<EOF
0474 TCanvas *c$i = new TCanvas("c$i");
0475 c$i->SetFillColor(10);
0476 file_old->cd("$HISTOPATHNAME_Photons");
0477 $i->SetStats(0);
0478 if ( $i==pEcalRecHitSumEtConeDR04VsEtaAll ) {
0479 $i->GetYaxis()->SetRangeUser(0.,5.);
0480 } else if ( $i==pEcalRecHitSumEtConeDR04VsEtBarrel )
0481 { $i->GetYaxis()->SetRangeUser(0.,20.);
0482 } else if ( $i==pEcalRecHitSumEtConeDR04VsEtEndcap )
0483 { $i->GetYaxis()->SetRangeUser(0.,20.);
0484 } else if ( $i==pHcalTowerSumEtConeDR04VsEtaAll )
0485 { $i->GetYaxis()->SetRangeUser(0.,0.3);
0486 } else if ( $i==pHcalTowerBcSumEtConeDR04VsEtaAll )
0487 { $i->GetYaxis()->SetRangeUser(0.,1.);
0488 } else if ( $i==pHcalTowerSumEtConeDR04VsEtBarrel || $i==pHcalTowerBcSumEtConeDR04VsEtBarrel)
0489 { $i->GetYaxis()->SetRangeUser(0.,5.);
0490 } else if ( $i==pHcalTowerSumEtConeDR04VsEtEndcap || $i==pHcalTowerBcSumEtConeDR04VsEtEndcap )
0491 { $i->GetYaxis()->SetRangeUser(0.,5.);
0492 } else if ( $i==pHoverEVsEtaAll || $i==pnewHoverEVsEtaAll )
0493 { $i->GetYaxis()->SetRangeUser(-0.05,0.05);
0494 } else if ( $i==pHoverEVsEtAll || $i==pnewHoverEVsEtAll )
0495 { $i->GetYaxis()->SetRangeUser(-0.05,0.05);
0496 } else {
0497 $i->SetMinimum(0.8);
0498 $i->SetMaximum(1.1);
0499 }
0500 $i->SetLineColor(kPink+8);
0501 $i->SetMarkerColor(kPink+8);
0502 $i->SetMarkerStyle(20);
0503 $i->SetMarkerSize(0.8);
0504 $i->SetLineWidth(1);
0505 $i->Draw();
0506 file_new->cd("$HISTOPATHNAME_Photons");
0507 $i->SetStats(0);
0508 $i->SetLineColor(kBlack);
0509 $i->SetMarkerColor(kBlack);
0510 $i->SetMarkerStyle(20);
0511 $i->SetMarkerSize(0.8);
0512 $i->SetLineWidth(1);
0513 $i->Draw("e1same");
0514 c$i->SaveAs("gifs/$i.gif");
0515
0516 EOF
0517 setenv N `expr $N + 1`
0518 end
0519
0520
0521
0522
0523
0524
0525 setenv NTOT `expr $N - 1`
0526 setenv N 1
0527 while ( $N <= $NTOT )
0528 cat temp$N.C >>& validation.C
0529 rm temp$N.C
0530 setenv N `expr $N + 1`
0531 end
0532
0533 cat > end.C <<EOF
0534 }
0535 EOF
0536 cat end.C >>& validation.C
0537 rm end.C
0538
0539
0540 if ( $TYPE == PixelMatchGsfElectron ) then
0541 setenv ANALYZER PixelMatchGsfElectronAnalyzer
0542 setenv CFG read_gsfElectrons
0543 else if ( $TYPE == Photons || $TYPE == GEDPhotons ) then
0544 setenv ANALYZER ${ANALYZERNAME}
0545 setenv CFG PhotonValidator_cfg
0546 endif
0547
0548 if (-e validation.html) rm validation.html
0549 if (-e validationPlotsTemplate.html) rm validationPlotsTemplate.html
0550 cp ${CURRENTDIR}/validationPlotsTemplate.html validationPlotsTemplate.html
0551 touch validation.html
0552
0553
0554 cat > begin.html <<EOF
0555 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2 Final//EN">
0556 <html>
0557 <head>
0558 <meta http-equiv="content-type" content="text/html; charset=UTF-8" />
0559 <title>$RELEASE : Full Sim vs Fast Sim $TYPE validation</title>
0560 </head>
0561
0562 <h1>$RELEASE : Full Sim vs Fast Sim $TYPE validation
0563 <br>
0564 $SAMPLE
0565 </h1>
0566
0567 <p>The following plots were made using <a href="http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/Validation/RecoEgamma/src/$ANALYZER.cc">Validation/RecoEgamma/src/$ANALYZER</a>,
0568 using <a href="http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/Validation/RecoEgamma/test/$CFG.py">Validation/RecoEgamma/test/$CFG.py</a>
0569 <p>The script used to make the plots is <a href="validation.C">here</a>.
0570 <br>
0571 In all plots below, FastSim is in purple, Full Sim in black.<br>
0572 WARNING: no plots are not yet available for ecal seeded conversions in FastSim.
0573 <br>
0574 Click on the plots to see them enlarged.
0575 <br>
0576 Responsible: N. Marinelli
0577 <br>
0578 <br>
0579
0580 EOF
0581
0582
0583 cat begin.html >>& validation.html
0584 rm begin.html
0585 cat validationPlotsTemplate.html >>& validation.html
0586 rm validationPlotsTemplate.html
0587
0588
0589 rm scaledhistosForPhotons
0590 rm unscaledhistosForPhotons
0591 rm efficiencyForPhotons
0592 rm scaledhistosForPhotonsLogScale
0593
0594 #echo "Now paste the following into your terminal window:"
0595 #echo ""
0596 echo "cd $OUTDIR"
0597 #echo " root -b"
0598 #echo ".x validation.C"
0599 #echo ".q"
0600 #echo "cd $CURRENTDIR"
0601 #echo ""
0602
0603
0604 root -b -l -q validation.C
0605 cd $CURRENTDIR
0606 echo "Then you can view your valdation plots here:"
0607 echo "http://cmsdoc.cern.ch/Physics/egamma/www/$OUTPATH/validation.html"