Back to home page

Project CMSSW displayed by LXR

 
 

    


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"