RecoEgamma/test/bkgValidation.csh

0001 #!/bin/csh
0002
0003 #This script can be used to generate a web page to compare histograms from
0004 #two input root files produced using the EDAnalyzers in RecoEgamma/Examples,
0005 #by running one of:
0006 #
0007 #
0008 #
0009 #  "Validation/RecoEgamma/test/PhotonValidator_cfg.py
0010 #
0011 # The default list of histograms (configurable) is based on version VXX-XX-XX
0012 # of Validation/RecoEgamma
0013 #
0014 #Two files are created by this script: validation.C and validation.html.
0015 #validation.C should be run inside root to greate a set of gif images
0016 #which can then be viewed in a web browser using validation.html.
0017
0018 #=============BEGIN CONFIGURATION=================
0019 setenv TYPE Photons
0020 setenv RUNTYPE Central
0021 setenv STARTUP True
0022 setenv CMSSWver1 4_2_0
0023 setenv CMSSWver2 4_2_2
0024 setenv OLDRELEASE 4_2_0
0025 setenv NEWRELEASE 4_2_2
0026 setenv OLDPRERELEASE
0027 setenv NEWPRERELEASE
0028
0029
0030 if ( $STARTUP == True) then
0031 setenv OLDGLOBALTAG START42_V9-v1
0032 setenv NEWGLOBALTAG START42_V11-v1
0033 else
0034 setenv OLDGLOBALTAG MC_42_V9-v1
0035 setenv NEWGLOBALTAG MC_42_V11-v1
0036 endif
0037
0038
0039
0040 #setenv OLDRELEASE ${OLDRELEASE}_${OLDPRERELEASE}
0041 setenv OLDRELEASE ${OLDRELEASE}
0042 #setenv NEWRELEASE ${NEWRELEASE}_${NEWPRERELEASE}
0043 setenv NEWRELEASE ${NEWRELEASE}
0044
0045
0046
0047 #setenv WorkDir1   /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver1}/src/Validation/RecoEgamma/test
0048 #setenv WorkDir2   /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver2}_${NEWPRERELEASE}/src/Validation/RecoEgamma/test
0049
0050
0051 #setenv WorkDir1   /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver1}_${OLDPRERELEASE}/src/Validation/RecoEgamma/test
0052 #setenv WorkDir2   /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver2}_${NEWPRERELEASE}/src/Validation/RecoEgamma/test
0053
0054
0055 #setenv WorkDir1   /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver1}_${OLDPRERELEASE}/src/Validation/RecoEgamma/test
0056 #setenv WorkDir2   /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver2}/src/Validation/RecoEgamma/test
0057
0058
0059 setenv WorkDir1   /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver1}/src/Validation/RecoEgamma/test
0060 setenv WorkDir2   /afs/cern.ch/user/n/nancy/scratch0/CMSSW/test/CMSSW_${CMSSWver2}/src/Validation/RecoEgamma/test
0061
0062 #Name of sample (affects output directory name and htmldescription only)
0063 setenv SAMPLE QCD_Pt_80_120STARTUP
0064 #setenv SAMPLE QCD_Pt_20_30STARTUP
0065
0066
0067 if ( $RUNTYPE == Central ) then
0068 setenv HISTOPATHNAME_Efficiencies DQMData/Run\ 1/EgammaV/Run\ summary/PhotonValidator/Efficiencies
0069 setenv HISTOPATHNAME_Photons DQMData/Run\ 1/EgammaV/Run\ summary/PhotonValidator/Photons
0070 setenv HISTOPATHNAME_Conversions DQMData/Run\ 1/EgammaV/Run\ summary/PhotonValidator/ConversionInfo
0071 setenv HISTOPATHNAME_Background DQMData/Run\ 1/EgammaV/Run\ summary/PhotonValidator/Background
0072 endif
0073
0074 if ( $RUNTYPE == Local ) then
0075 setenv HISTOPATHNAME_Efficiencies DQMData/EgammaV/PhotonValidator/Efficiencies
0076 setenv HISTOPATHNAME_Photons DQMData/EgammaV/PhotonValidator/Photons
0077 setenv HISTOPATHNAME_Conversions DQMData/EgammaV/PhotonValidator/ConversionInfo
0078 setenv HISTOPATHNAME_Background DQMData/EgammaV/PhotonValidator/Background
0079 endif
0080 #==============END BASIC CONFIGURATION==================
0081
0082
0083 #Input root trees for the two cases to be compared
0084
0085
0086 if ($SAMPLE == QCD_Pt_20_30STARTUP) then
0087
0088 setenv OLDFILE ${WorkDir1}/PhotonValidationRelVal${OLDRELEASE}_QCD_Pt_20_30.root
0089 setenv NEWFILE ${WorkDir2}/PhotonValidationRelVal${NEWRELEASE}_QCD_Pt_20_30.root
0090
0091 else if ($SAMPLE == QCD_Pt_80_120STARTUP) then
0092 if ( $RUNTYPE == Local ) then
0093 setenv OLDFILE ${WorkDir1}/PhotonValidationRelVal${OLDRELEASE}_QCD_Pt_80_120.root
0094 setenv NEWFILE ${WorkDir2}/PhotonValidationRelVal${NEWRELEASE}_QCD_Pt_80_120.root
0095 else if ( $RUNTYPE == Central ) then
0096 setenv OLDFILE ${WorkDir1}/DQM_V0003_R000000001__RelValQCD_Pt_80_120__CMSSW_${OLDRELEASE}-${OLDGLOBALTAG}__DQM.root
0097 setenv NEWFILE ${WorkDir2}/DQM_V0001_R000000001__RelValQCD_Pt_80_120__CMSSW_${NEWRELEASE}-${NEWGLOBALTAG}__DQM.root
0098 endif
0099
0100
0101
0102 #Location of output.  The default will put your output in:
0103 #http://cmsdoc.cern.ch/Physics/egamma/www/validation/
0104
0105 setenv CURRENTDIR $PWD
0106 setenv OUTPATH /afs/cern.ch/cms/Physics/egamma/www/validation
0107 cd $OUTPATH
0108 if (! -d $NEWRELEASE) then
0109   mkdir $NEWRELEASE
0110 endif
0111 setenv OUTPATH $OUTPATH/$NEWRELEASE
0112 cd $OUTPATH
0113
0114 if (! -d ${TYPE}) then
0115   mkdir ${TYPE}
0116 endif
0117 setenv OUTPATH $OUTPATH/${TYPE}
0118 cd  $OUTPATH
0119
0120 if (! -d vs${OLDRELEASE}) then
0121   mkdir vs${OLDRELEASE}
0122 endif
0123 setenv OUTPATH $OUTPATH/vs${OLDRELEASE}
0124
0125
0126 setenv OUTDIR $OUTPATH/${SAMPLE}
0127 if (! -d $OUTDIR) then
0128   cd $OUTPATH
0129   mkdir $OUTDIR
0130   cd $OUTDIR
0131   mkdir gifs
0132 endif
0133 cd $OUTDIR
0134
0135
0136 #The list of histograms to be compared for each TYPE can be configured below:
0137
0138
0139 if ( $TYPE == Photons ) then
0140
0141 cat > efficiencyForBkg <<EOF
0142   bkgEffVsEta
0143   bkgEffVsPhi
0144   bkgEffVsEt
0145   deadChVsEtaBkg
0146   deadChVsPhiBkg
0147   deadChVsEtBkg
0148 EOF
0149
0150
0151
0152 cat > scaledhistosForBkg <<EOF
0153
0154   nOfPhotons
0155   scBkgEta
0156   scBkgPhi
0157   scBkgEAll
0158   scBkgEtAll
0159   phoBkgEta
0160   phoBkgPhi
0161   phoBkgEAll
0162   phoBkgEtAll
0163   phoBkgDEta
0164   phoBkgDPhi
0165   isoTrkSolidConeDR04BkgBarrel
0166   isoTrkSolidConeDR04BkgEndcap
0167   nTrkSolidConeDR04BkgBarrel
0168   nTrkSolidConeDR04BkgEndcap
0169
0170   convEtaBkg
0171   convPhiBkg
0172   PoverEtracksBkgAll
0173   PoverEtracksBkgBarrel
0174   PoverEtracksBkgEndcap
0175   mvaOutBkgAll
0176   mvaOutBkgBarrel
0177   mvaOutBkgEndcap
0178
0179   hDPhiTracksAtVtxBkgAll
0180   hDCotTracksBkgAll
0181
0182 EOF
0183
0184
0185 cat > scaledhistosForBkgLogScale <<EOF
0186   hOverEBkgAll
0187   hOverEBkgBarrel
0188   hOverEBkgEndcap
0189   hcalTowerSumEtConeDR04BkgBarrel
0190   hcalTowerSumEtConeDR04BkgEndcap
0191   EoverPtracksBkgAll
0192   EoverPtracksBkgBarrel
0193   EoverPtracksBkgEndcap
0194   r9BkgBarrel
0195   r9BkgEndcap
0196   r1BkgBarrel
0197   r1BkgEndcap
0198   r2BkgBarrel
0199   r2BkgEndcap
0200   sigmaIetaIetaBkgBarrel
0201   sigmaIetaIetaBkgEndcap
0202   ecalRecHitSumEtConeDR04BkgBarrel
0203   ecalRecHitSumEtConeDR04BkgEndcap
0204
0205
0206 EOF
0207
0208
0209
0210 cat > unscaledhistosForBkg <<EOF
0211
0212   pR1VsEtaBkgAll
0213   pR2VsEtaBkgAll
0214   pR1VsEtBkgAll
0215   pR2VsEtBkgAll
0216   pSigmaIetaIetaVsEtaBkgAll
0217   pSigmaIetaIetaVsEtBkgAll
0218   pHOverEVsEtaBkgAll
0219   pHOverEVsEtBkgAll
0220   pEcalRecHitSumEtConeDR04VsEtBkgBarrel
0221   pEcalRecHitSumEtConeDR04VsEtBkgEndcap
0222   pEcalRecHitSumEtConeDR04VsEtaBkgAll
0223   pHcalTowerSumEtConeDR04VsEtBkgBarrel
0224   pHcalTowerSumEtConeDR04VsEtBkgEndcap
0225   pHcalTowerSumEtConeDR04VsEtaBkgAll
0226   pIsoTrkSolidConeDR04VsEtBkgBarrel
0227   pIsoTrkSolidConeDR04VsEtBkgEndcap
0228   pIsoTrkSolidConeDR04VsEtaBkgAll
0229   pnTrkSolidConeDR04VsEtBkgBarrel
0230   pnTrkSolidConeDR04VsEtBkgEndcap
0231   p_nTrkSolidConeDR04VsEtaBkgAll
0232
0233
0234
0235 EOF
0236
0237 cat > 2DhistosForBkg <<EOF
0238   R9VsEtaBkgAll
0239   R9VsEtBkgAll
0240   hOverEVsEtaBkgAll
0241   hOverEVsEtBkgAll
0242
0243
0244 EOF
0245
0246
0247 endif
0248
0249 #=================END CONFIGURATION=====================
0250
0251 if (-e validation.C) rm validation.C
0252 touch validation.C
0253 cat > begin.C <<EOF
0254 {
0255 TFile *file_old = TFile::Open("$OLDFILE");
0256 TFile *file_new = TFile::Open("$NEWFILE");
0257
0258 EOF
0259 cat begin.C >>& validation.C
0260 rm begin.C
0261
0262 setenv N 1
0263
0264
0265 foreach i (`cat efficiencyForBkg`)
0266   cat > temp$N.C <<EOF
0267
0268 TCanvas *c$i = new TCanvas("c$i");
0269 c$i->SetFillColor(10);
0270 //file_old->cd("DQMData/EgammaV/PhotonValidator/Efficiencies");
0271 file_old->cd("$HISTOPATHNAME_Efficiencies");
0272 $i->SetStats(0);
0273 if ( $i==deadChVsEtaBkg ||  $i==deadChVsPhiBkg ||  $i==deadChVsEtBkg ) {
0274 $i->GetYaxis()->SetRangeUser(0.,0.2);
0275
0276 }else if ( $i==bkgEffVsEta ||  $i==bkgEffVsPhi ) {
0277 $i->GetYaxis()->SetRangeUser(0.,0.4);
0278 }else if (  $i==bkgEffVsEt ) {
0279 $i->GetYaxis()->SetRangeUser(0.,1.);
0280 } else {
0281 $i->GetYaxis()->SetRangeUser(0.,1.1);
0282 }
0283 $i->SetLineColor(kPink+8);
0284 $i->SetMarkerColor(kPink+8);
0285 $i->SetMarkerStyle(20);
0286 $i->SetMarkerSize(1);
0287 $i->SetLineWidth(1);
0288 $i->Draw("e1");
0289
0290 //file_new->cd("DQMData/EgammaV/PhotonValidator/Efficiencies");
0291 file_new->cd("$HISTOPATHNAME_Efficiencies");
0292 $i->SetStats(0);
0293 $i->SetMinimum(0.);
0294 $i->SetMaximum(1.1);
0295 $i->SetLineColor(kBlack);
0296 $i->SetMarkerColor(kBlack);
0297 $i->SetMarkerStyle(20);
0298 $i->SetMarkerSize(1);
0299 $i->SetLineWidth(1);
0300 $i->Draw("e1same");
0301 c$i->SaveAs("gifs/$i.gif");
0302
0303 EOF
0304   setenv N `expr $N + 1`
0305 end
0306
0307
0308
0309
0310
0311 foreach i (`cat scaledhistosForBkg`)
0312   cat > temp$N.C <<EOF
0313 TCanvas *c$i = new TCanvas("c$i");
0314 c$i->SetFillColor(10);
0315 //file_new->cd("DQMData/EgammaV/PhotonValidator/Background");
0316 file_new->cd("$HISTOPATHNAME_Background");
0317 Double_t mnew=$i->GetMaximum();
0318 Double_t nnew=$i->GetEntries();
0319 //file_old->cd("DQMData/EgammaV/PhotonValidator/Background");
0320 file_old->cd("$HISTOPATHNAME_Background");
0321 Double_t mold=$i->GetMaximum();
0322 Double_t nold=$i->GetEntries();
0323 $i->SetStats(0);
0324 $i->SetMinimum(0.);
0325 //if ( mnew > mold)
0326 // $i->SetMaximum(mnew+mnew*0.2);
0327 //else
0328 //$i->SetMaximum(mold+mold*0.2);
0329 //$i->SetMaximum(mold+mold*0.2);
0330 $i->SetLineColor(kPink+8);
0331 $i->SetFillColor(kPink+8);
0332 //$i->SetLineWidth(3);
0333 $i->Draw();
0334 //file_new->cd("DQMData/EgammaV/PhotonValidator/Background");
0335 file_new->cd("$HISTOPATHNAME_Background");
0336 Double_t nnew=$i->GetEntries();
0337 $i->SetStats(0);
0338 $i->SetLineColor(kBlack);
0339 $i->SetMarkerColor(kBlack);
0340 $i->SetMarkerStyle(20);
0341 $i->SetMarkerSize(1);
0342 //$i->SetLineWidth(1);
0343 $i->Scale(nold/nnew);
0344 $i->Draw("e1same");
0345 c$i->SaveAs("gifs/$i.gif");
0346
0347 EOF
0348   setenv N `expr $N + 1`
0349 end
0350
0351
0352
0353 foreach i (`cat unscaledhistosForBkg`)
0354   cat > temp$N.C <<EOF
0355 TCanvas *c$i = new TCanvas("c$i");
0356 c$i->SetFillColor(10);
0357 //file_old->cd("DQMData/EgammaV/PhotonValidator/Background");
0358 file_old->cd("$HISTOPATHNAME_Background");
0359 $i->SetStats(0);
0360 if ( $i==pEcalRecHitSumEtConeDR04VsEtaBkgAll ||  $i==pHcalTowerSumEtConeDR04VsEtaBkgAll  ) {
0361 $i->GetYaxis()->SetRangeUser(0.,25.);
0362 } else if ( $i==pEcalRecHitSumEtConeDR04VsEtBkgBarrel || $i==pHcalTowerSumEtConeDR04VsEtBkgBarrel )
0363 { $i->GetYaxis()->SetRangeUser(0.,30.);
0364 } else if (  $i==p_nTrkSolidConeDR04VsEtaBkgAll || $i==pnTrkSolidConeDR04VsEtBkgBarrel ||   $i==pnTrkSolidConeDR04VsEtBkgEndcap )
0365 { $i->GetYaxis()->SetRangeUser(0.,20.);
0366 } else if (   $i==pIsoTrkSolidConeDR04VsEtaBkgAll ||  $i==pIsoTrkSolidConeDR04VsEtBkgBarrel || $i==pIsoTrkSolidConeDR04VsEtBkgEndcap)
0367 { $i->GetYaxis()->SetRangeUser(0.,100.);
0368 } else if ( $i==pEcalRecHitSumEtConeDR04VsEtBkgEndcap || $i==pHcalTowerSumEtConeDR04VsEtBkgEndcap )
0369 {$i->GetYaxis()->SetRangeUser(0.,30.);
0370 } else if ( $i==pSigmaIetaIetaVsEtaBkgAll || $i==pSigmaIetaIetaVsEtBkgAll ||  $i==pHOverEVsEtaBkgAll ||  $i==pHOverEVsEtBkgAll )
0371 { $i->GetYaxis()->SetRangeUser(0.,0.1);
0372 } else  {
0373 $i->SetMinimum(0.);
0374 $i->SetMaximum(1.);
0375 }
0376 $i->SetLineColor(kPink+8);
0377 $i->SetMarkerColor(kPink+8);
0378 $i->SetMarkerStyle(20);
0379 $i->SetMarkerSize(1);
0380 $i->SetLineWidth(1);
0381 $i->Draw();
0382 //file_new->cd("DQMData/EgammaV/PhotonValidator/Background");
0383 file_new->cd("$HISTOPATHNAME_Background");
0384 $i->SetStats(0);
0385 $i->SetLineColor(kBlack);
0386 $i->SetMarkerColor(kBlack);
0387 $i->SetMarkerStyle(20);
0388 $i->SetMarkerSize(1);
0389 $i->SetLineWidth(1);
0390 $i->Draw("e1same");
0391 c$i->SaveAs("gifs/$i.gif");
0392
0393 EOF
0394   setenv N `expr $N + 1`
0395 end
0396
0397
0398
0399 foreach i (`cat scaledhistosForBkgLogScale`)
0400   cat > temp$N.C <<EOF
0401 TCanvas *c$i = new TCanvas("c$i");
0402 c$i->SetFillColor(10);
0403 c$i->SetLogy(1);
0404 //file_new->cd("DQMData/EgammaV/PhotonValidator/Background");
0405 file_new->cd("$HISTOPATHNAME_Background");
0406 Double_t nnew=$i->GetEntries();
0407 //file_old->cd("DQMData/EgammaV/PhotonValidator/Background");
0408 file_old->cd("$HISTOPATHNAME_Background");
0409 if ( $i==hcalTowerSumEtConeDR04BkgBarrel ||  $i==hcalTowerSumEtConeDR04BkgEndcap  ) {
0410 $i->GetXaxis()->SetRangeUser(0.,10.);
0411 } else if ( $i==hOverEBkgBarrel || $i==hOverEBkgEndcap ) {
0412 $i->GetXaxis()->SetRangeUser(0.,1.);
0413 }
0414
0415 Double_t nold=$i->GetEntries();
0416 $i->SetStats(0);
0417 $i->SetLineColor(kPink+8);
0418 $i->SetFillColor(kPink+8);
0419 $i->Draw();
0420 //file_new->cd("DQMData/EgammaV/PhotonValidator/Background");
0421 file_new->cd("$HISTOPATHNAME_Background");
0422 Double_t nnew=$i->GetEntries();
0423 $i->SetStats(0);
0424 $i->SetLineColor(kBlack);
0425 $i->SetMarkerColor(kBlack);
0426 $i->SetMarkerStyle(20);
0427 $i->SetMarkerSize(1);
0428 $i->Scale(nold/nnew);
0429 $i->Draw("e1same");
0430 c$i->SaveAs("gifs/$i.gif");
0431
0432 EOF
0433   setenv N `expr $N + 1`
0434 end
0435
0436
0437
0438
0439
0440 foreach i (`cat 2DhistosForBkg`)
0441   cat > temp$N.C <<EOF
0442 TCanvas *c$i = new TCanvas("c$i");
0443 c$i->SetFillColor(10);
0444 //file_old->cd("DQMData/EgammaV/PhotonValidator/Background");
0445 file_old->cd("$HISTOPATHNAME_Background");
0446 $i->SetStats(0);
0447 $i->SetMinimum(0.);
0448 $i->SetMarkerColor(kPink+8);
0449 $i->SetMarkerStyle(2);
0450 $i->SetMarkerSize(0.2);
0451 $i->Draw();
0452 //file_new->cd("DQMData/EgammaV/PhotonValidator/Background");
0453 file_new->cd("$HISTOPATHNAME_Background");
0454 $i->SetStats(0);
0455 $i->SetMarkerColor(kBlack);
0456 $i->SetMarkerStyle(2);
0457 $i->SetMarkerSize(0.2);
0458 $i->Draw("same");
0459 c$i->SaveAs("gifs/$i.gif");
0460
0461 EOF
0462   setenv N `expr $N + 1`
0463 end
0464
0465
0466
0467
0468 setenv NTOT `expr $N - 1`
0469 setenv N 1
0470 while ( $N <= $NTOT )
0471   cat temp$N.C >>& validation.C
0472   rm temp$N.C
0473   setenv N `expr $N + 1`
0474 end
0475
0476 cat > end.C <<EOF
0477 }
0478 EOF
0479 cat end.C >>& validation.C
0480 rm end.C
0481
0482
0483 if ( $TYPE == PixelMatchGsfElectron ) then
0484   setenv ANALYZER PixelMatchGsfElectronAnalyzer
0485   setenv CFG read_gsfElectrons
0486 else if ( $TYPE == Photons ) then
0487   setenv ANALYZER PhotonValidator
0488   setenv CFG PhotonValidator_cfg
0489 endif
0490
0491 if (-e validation.html) rm validation.html
0492 if (-e bkgValidationPlotsTemplate.html) rm bkgValidationPlotsTemplate.html
0493 cp ${CURRENTDIR}/bkgValidationPlotsTemplate.html bkgValidationPlotsTemplate.html
0494 touch validation.html
0495 cat > begin.html <<EOF
0496 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2 Final//EN">
0497 <html>
0498 <head>
0499 <meta http-equiv="content-type" content="text/html; charset=UTF-8" />
0500 <title>$NEWRELEASE vs $OLDRELEASE $TYPE validation</title>
0501 </head>
0502
0503 <h1>$NEWRELEASE vs $OLDRELEASE $TYPE validation
0504 <br>
0505  $SAMPLE
0506 </h1>
0507
0508 <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>,
0509 using <a href="http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/Validation/RecoEgamma/test/$CFG.py">Validation/RecoEgamma/test/$CFG.py</a>
0510 <p>The script used to make the plots is <a href="validation.C">here</a>.
0511 <br>
0512 In all plots below, $OLDRELEASE is in purple , $NEWRELEASE in black.
0513 <br>
0514 Click on the plots to see them enlarged.
0515 <br>
0516 Responsible: N. Marinelli
0517 <br>
0518 <br>
0519
0520
0521
0522 EOF
0523 cat begin.html >>& validation.html
0524 rm begin.html
0525 cat  bkgValidationPlotsTemplate.html >>& validation.html
0526 rm   bkgValidationPlotsTemplate.html
0527
0528 rm efficiencyForBkg
0529 rm scaledhistosForBkg
0530 rm scaledhistosForBkgLogScale
0531 rm unscaledhistosForBkg
0532 rm 2DhistosForBkg
0533
0534 #echo "Now paste the following into your terminal window:"
0535 #echo ""
0536 echo "cd $OUTDIR"
0537 #echo " root -b"
0538 #echo ".x validation.C"
0539 #echo ".q"
0540 #echo "cd $CURRENTDIR"
0541 #echo ""
0542
0543
0544 root -b -l -q validation.C
0545 cd $CURRENTDIR
0546 echo "Then you can view your valdation plots here:"
0547 echo "http://cmsdoc.cern.ch/Physics/egamma/www/$OUTPATH/validation.html"