EcalTestBeam/test/write_crystal_cfg.pl

0001 # build the cff file for all the crystals corresponding to
0002 # their nominal maximum reference point as stored in the table
0003 # F. Cossutti - 23-Aug-2006 15:19
0004
0005
0006 use diagnostics;
0007 use strict;
0008
0009 my $file = "BarrelSM1CrystalCenterElectron120GeV.dat";
0010
0011 if (! open(INPUT,"<$file") ) {
0012     print STDERR "Can't open input file $file: $!\n";
0013     exit 1;
0014    }
0015
0016     while(<INPUT>){
0017         chomp(my $line = $_);
0018
0019         print $line, "\n";
0020
0021         my $crystal = my $eta = my $phi = ();
0022
0023         ($crystal,$eta,$phi) = split(' ', $line, 3);
0024
0025         my $filename = "crystal".$crystal."_cff.py";
0026
0027         open(OUTFILE, ">$filename");
0028         print OUTFILE "import FWCore.ParameterSet.Config as cms\n";
0029         print OUTFILE "                                        \n";
0030         print OUTFILE "common_beam_direction_parameters = cms.PSet(\n";
0031         print OUTFILE "    MinEta = cms.double(".$eta,"),\n";
0032         print OUTFILE "    MaxEta = cms.double(".$eta,"),\n";
0033         print OUTFILE "    MinPhi = cms.double(".$phi,"),\n";
0034         print OUTFILE "    MaxPhi = cms.double(".$phi,"),\n";
0035         print OUTFILE "    Psi    = cms.double(999.9),\n";
0036         print OUTFILE "    BeamMeanX = cms.double(0.0),\n";
0037         print OUTFILE "    BeamMeanY = cms.double(0.0),\n";
0038         print OUTFILE "    BeamPosition = cms.double(-26733.5)\n";
0039         print OUTFILE ")\n";
0040
0041     }
0042     close(INPUT);
0043
0044 exit 0;