epic-lfhcal-tbana/NewStructure/runCAENCalibration_2024.sh

0001 #! /bin/bash
0002
0003 PlotBaseDir=..
0004
0005 function MuonCalib()
0006 {
0007         runNrFile='../configs/TB2024/DataTakingDB_202409_CAEN.csv'
0008         echo "=================================================================================="
0009         echo "option $1"
0010         echo "run Nr Pedestal: $2"
0011         echo "run Nr Muon: $3"
0012         echo "dataRawDir: $4"
0013         echo "dataOutDir: $5"
0014         echo "OutNameRun:" $6
0015         if [ $1 == "transfer" ]; then
0016                 echo "badchannelMap:" $7
0017         else
0018     echo "layer skip:" $8
0019         fi
0020         echo "=================================================================================="
0021         if [ $1 == "transfer" ]; then
0022                 ./DataPrep -d 1 -e -f -P $5/PedestalCalib_$2.root -i $4/raw_$3.root -o $4/rawPed_$3.root -O $PlotBaseDir/CAEN_PlotsCalibTransfer_2024/$6 -r $runNrFile
0023         elif [ $1 == "default" ]; then
0024                 time ./DataPrep -f -d 1  -s -i $4/rawPed_$3.root -o $5/rawPedAndMuon_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuon_2024/$6 -r $runNrFile -l $8
0025         elif [ $1 == "imp1st" ]; then
0026                 time ./DataPrep -f -d 1  -S -i $5/rawPedAndMuon_$3.root -o $5/rawPedAndMuonImp_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImproved_2024/$6 -r $runNrFile -l $8
0027         elif [ $1 == "imp2nd" ]; then
0028                 ./DataPrep -f -d 1  -S -i $5/rawPedAndMuonImp_$3.root -o $5/rawPedAndMuonImp2nd_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImproved_2024/$6_2ndIte -r $runNrFile -l $8
0029         elif [ $1 == "imp3rd" ]; then
0030                 ./DataPrep -f -d 1  -S -i $5/rawPedAndMuonImp2nd_$3.root -o $5/rawPedAndMuonImp3rd_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImproved_2024/$6_3rdIte -r $runNrFile -l $8
0031         elif [ $1 == "transferWithBC" ]; then
0032                 time ./DataPrep -d 1 -a -e -f -P $5/PedestalCalib_$2.root -i $4/raw_$3.root -o $4/rawPedWBC_$3.root -B $7 -O $PlotBaseDir/CAEN_PlotsCalibTransferBC_2024/$6 -r $runNrFile -l $8
0033         elif [ $1 == "defaultWithBC" ]; then
0034                 time ./DataPrep -e -a -f -d 1  -s -i $4/rawPedWBC_$3.root -o $5/rawPedAndMuonWBC_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonBC_2024/$6 -r $runNrFile -l $8
0035         elif [ $1 == "impWBC1st" ]; then
0036                 time ./DataPrep -f -d 1 -a -S -i $5/rawPedAndMuonWBC_$3.root -o $5/rawPedAndMuonWBCImp_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImprovedBC_2024/$6 -r $runNrFile -l $8
0037         elif [ $1 == "impWBC2nd" ]; then
0038                 time ./DataPrep -f -d 1 -a -S -i $5/rawPedAndMuonWBCImp_$3.root -o $5/rawPedAndMuonWBCImp2nd_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImprovedBC_2024/$6_2ndIte -r $runNrFile -l $8
0039         elif [ $1 == "impWBC3rd" ]; then
0040                 time ./DataPrep -f -d 1 -a -S -i $5/rawPedAndMuonWBCImp2nd_$3.root -o $5/rawPedAndMuonWBCImp3rd_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImprovedBC_2024/$6_3rdIte -r $runNrFile -l $8
0041         elif [ $1 == "impWBC4th" ]; then
0042                 time ./DataPrep -f -e -d 1 -a -S -i $5/rawPedAndMuonWBCImp3rd_$3.root -o $5/rawPedAndMuonWBCImp4th_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImprovedBC_2024/$6_4thIte -r $runNrFile -l $8
0043         elif [ $1 == "impWBC5th" ]; then
0044                 time ./DataPrep -f -e -d 1 -a -S -i $5/rawPedAndMuonWBCImp4th_$3.root -o $5/rawPedAndMuonWBCImp5th_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImprovedBC_2024/$6_5thIte -r $runNrFile -l $8
0045         elif [ $1 == "impWBC6th" ]; then
0046                 time ./DataPrep -f -d 1 -a -S -i $5/rawPedAndMuonWBCImp5th_$3.root -o $5/rawPedAndMuonWBCImp6th_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImprovedBC_2024/$6_6thIte -r $runNrFile -l $8
0047         elif [ $1 == "noise" ]; then
0048                 ./DataPrep -f -d 1  -n -i $5/rawPedAndMuon_$3.root -o $5/rawPedAndMuonNoise_$2.root -O $PlotBaseDir/CAEN_PlotsCalibNoiseRe_2024/$6 -r $runNrFile
0049         elif [ $1 == "transferAlt" ]; then
0050                 ./DataPrep -d 1  -f -P $5/rawPedAndMuonNoise_$3.root -i $4/raw_$3.root -o $4/rawPedImp_$3.root -r $runNrFile
0051         elif [ $1 == "defaultImpPed" ]; then
0052                 ./DataPrep -f -d 1  -s -i $4/rawPedImp_$3.root -o $5/rawPedImpAndMuon_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonPedImp_2024/$6 -r $runNrFile
0053         elif [ $1 == "saveNewPed" ]; then
0054                 ./DataPrep -f -d 1 -N -i $5/rawPedAndMuonNoise_$3.root -o $4/raw_pedonly_$3.root
0055         elif [ $1 == "saveNewMuon" ]; then
0056                 time ./DataPrep -f -d 1 -M -i $5/rawPedAndMuonWBC_$3.root -o $4/raw_muononly_$3.root
0057         elif [ $1 == "imp1st_red" ]; then
0058                 time ./DataPrep -f -d 1  -S -i $4/raw_muononly_$3.root -o $5/rawPedAndMuonImp1st_red_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImproved_2024/$6_Red1st -r $runNrFile -l $8
0059         elif [ $1 == "imp2nd_red" ]; then
0060                 time ./DataPrep -f -d 1  -S -i $5/rawPedAndMuonImp1st_red_$3.root -o $5/rawPedAndMuonImp2nd_red_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImproved_2024/$6_Red2nd -r $runNrFile -l $8
0061         elif [ $1 == "imp3rd_red" ]; then
0062                 time ./DataPrep -f -d 1  -S -i $5/rawPedAndMuonImp2nd_red_$3.root -o $5/rawPedAndMuonImp3rd_red_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImproved_2024/$6_Red3rd -r $runNrFile -l $8
0063         elif [ $1 == "imp4th_red" ]; then
0064                 time ./DataPrep -f -d 1  -S -i $5/rawPedAndMuonImp3rd_red_$3.root -o $5/rawPedAndMuonImp4th_red_$3.root -O $PlotBaseDir/CAEN_PlotsCalibMuonImproved_2024/$6_Red4th -r $runNrFile -l $8
0065         fi
0066 }
0067
0068 # running example:
0069 # bash runCalibration_2024.sh fbockExt2 muoncalibA1 improvedWBC4th
0070
0071
0072 dataDirRaw=""
0073 dataDirOut=""
0074 if [ $1 = "fbock" ]; then
0075         dataDirRaw=/media/fbock/Samsung_T5/LFHCAL_TB/202408_PST09/CAENdata/MuonRuns
0076         dataDirRawE=/media/fbock/Samsung_T5/LFHCAL_TB/202408_PST09/CAENdata/ElectronRuns
0077         dataDirRawH=/media/fbock/Samsung_T5/LFHCAL_TB/202408_PST09/CAENdata/HadronRuns
0078         dataDirOut=/home/fbock/EIC/Analysis/LFHCalTB2024/CAENdata/MuonRuns
0079         dataDirOutE=/home/fbock/EIC/Analysis/LFHCalTB2024/CAENdata/ElectronRuns
0080         dataDirOutH=/home/fbock/EIC/Analysis/LFHCalTB2024/CAENdata/HadronRuns
0081         PlotBaseDir=..
0082 elif [ $1 = "fbockExt" ]; then
0083         dataDirRaw=/media/fbock/T7/LFHCalTBData/202408_PST09/CAENData
0084         dataDirRawE=/media/fbock/T7/LFHCalTBData/202408_PST09/CAENData
0085         dataDirRawH=/media/fbock/T7/LFHCalTBData/202408_PST09/CAENData
0086         dataDirOut=/media/fbock/T7/LFHCalTBData/202408_PST09/CAENData
0087         dataDirOutE=/media/fbock/T7/LFHCalTBData/202408_PST09/CAENData
0088         dataDirOutH=/media/fbock/T7/LFHCalTBData/202408_PST09/CAENData
0089         PlotBaseDir=/media/fbock/T7/LFHCalTBData/202408_PST09/ReanlysisCAEN
0090 elif [ $1 = "eglimos" ]; then
0091         dataDirRaw=/home/ewa/EIC/test_beam2024/fullScanC
0092         dataDirRawE=/home/ewa/EIC/test_beam2024/fullScanC
0093         dataDirRawH=/home/ewa/EIC/test_beam2024/fullScanC
0094         dataDirOut=/home/ewa/EIC/test_beam2024/fullScanC/Output/MuonRuns
0095         dataDirOutE=/home/ewa/EIC/test_beam2024/fullScanC/Output/ElectronRuns
0096         dataDirOutH=/home/ewa/EIC/test_beam2024/fullScanC/Output/HadronRuns
0097         PlotBaseDir=..
0098 elif [ $1 = "rjh78" ]; then
0099         dataDirRaw=/Users/ryanhamilton/Documents/Research/data.nosync/202408_PST09/CAENData/outfiles/MuonRuns
0100         dataDirRawE=/Users/ryanhamilton/Documents/Research/data.nosync/202408_PST09/CAENData/outfiles/ElectronRuns
0101         dataDirRawH=/Users/ryanhamilton/Documents/Research/data.nosync/202408_PST09/CAENData/outfiles/HadronRuns
0102         dataDirOut=/Users/ryanhamilton/Documents/Research/data.nosync/202408_PST09/CAENData/outfiles/MuonRuns
0103         dataDirOutE=/Users/ryanhamilton/Documents/Research/data.nosync/202408_PST09/CAENData/outfiles/ElectronRuns
0104         dataDirOutH=/Users/ryanhamilton/Documents/Research/data.nosync/202408_PST09/CAENData/outfiles/HadronRuns
0105         PlotBaseDir=/Users/ryanhamilton/Documents/Research/data.nosync/202408_PST09/CAENData/outfiles/plots
0106
0107 elif [ $1 = "egpott" ]; then
0108         dataDirRaw=/Users/egpott/rhig/lfhcal/data/outfiles/MuonRuns
0109         dataDirRawE=/Users/egpott/rhig/lfhcal/data/outfiles/ElectronRuns
0110         dataDirRawH=/Users/egpott/rhig/lfhcal/data/outfiles/HadronRuns
0111         dataDirOut=/Users/egpott/rhig/lfhcal/data/outfiles/MuonRuns
0112         dataDirOutE=/Users/egpott/rhig/lfhcal/data/outfiles/ElectronRuns
0113         dataDirOutH=/Users/egpott/rhig/lfhcal/data/outfiles/ElectronRuns
0114         PlotBaseDir=/Users/egpott/rhig/lfhcal/data/outfiles/plots
0115
0116 else
0117         echo "Please select a known user name, otherwise I don't know where the data is"
0118         exit
0119 fi
0120
0121 skipLayer=0
0122 if [ $# -eq 4 ]; then
0123   skipLayer=$4
0124 fi
0125
0126 runList=../configs/TB2024/DataTakingDB_202409_CAEN.csv
0127
0128 # *-- Enable lines for the runs of interest
0129 # pedestal runs
0130 # pedestalRuns='271 277 303 306 308 311 315 332 369 377 404 420 454 465 476 492 505 521 528 552 553 ' # all pedestal runs
0131 #pedestalRuns='303 306 308 311 315 420 553 332 369 377 404 465 476 492 505 521' # all pedestal runs
0132 # pedestalRuns='271 277 454 528 552' # pedestal runs 45V
0133 # pedestalRuns='377 404'
0134 if [ $2 = "pedestal" ]; then
0135
0136   if [ $3 = "FullSetA" ]; then
0137     pedestalRuns='277 271'
0138   elif [ $3 = "FullSetB" ]; then
0139     pedestalRuns='332 369'
0140   fi
0141
0142         for runNr in $pedestalRuns; do
0143                 ./DataPrep -d 1 -p -i $dataDirRaw/raw_$runNr.root -f -o $dataDirOut/PedestalCalib_$runNr.root -O $PlotBaseDir/CAEN_PlotsPedestal_2024/Run$runNr -r $runList
0144         done;
0145 fi
0146
0147 if [ $2 == "mergemuons" ]; then
0148
0149
0150         hadd -f $dataDirRaw/raw_muonScanA1_45V.root $dataDirRaw/raw_244.root $dataDirRaw/raw_250.root
0151 #       hadd -f $dataDirRaw/raw_muonScanA2_45V.root $dataDirRaw/raw_283.root $dataDirRaw/raw_282.root
0152
0153         # hadd -f $dataDirRaw/raw_muonScanD1_45V.root $dataDirRaw/raw_412.root $dataDirRaw/raw_417.root
0154 #       hadd -f $dataDirRaw/raw_muonScanD2_45V.root $dataDirRaw/raw_460.root $dataDirRaw/raw_456.root $dataDirRaw/raw_457.root
0155 #       hadd -f $dataDirRaw/raw_muonScanH1_45V.root $dataDirRaw/raw_526.root $dataDirRaw/raw_527.root
0156 #       hadd -f $dataDirRaw/raw_muonScanH2_45V.root $dataDirRaw/raw_554.root $dataDirRaw/raw_559.root
0157 #       hadd -f $dataDirRaw/raw_muonScanC1_43_5V.root $dataDirRaw/raw_376.root $dataDirRaw/raw_375.root
0158 #       hadd -f $dataDirRaw/raw_muonScanC2_43_5V.root $dataDirRaw/raw_405.root $dataDirRaw/raw_410.root $dataDirRaw/raw_408.root
0159 #       hadd -f $dataDirRaw/raw_muonScanE1_40V.root $dataDirRaw/raw_464.root $dataDirRaw/raw_463.root
0160 #       hadd -f $dataDirRaw/raw_muonScanE2_40V.root $dataDirRaw/raw_481.root $dataDirRaw/raw_478.root
0161 #       hadd -f $dataDirRaw/raw_muonScanF1_41V.root $dataDirRaw/raw_486.root $dataDirRaw/raw_489.root
0162 #       hadd -f $dataDirRaw/raw_muonScanF2_41V.root $dataDirRaw/raw_507.root $dataDirRaw/raw_506.root
0163 #       hadd -f $dataDirRaw/raw_muonScanG_46V.root $dataDirRaw/raw_508.root $dataDirRaw/raw_510.root $dataDirRaw/raw_511.root $dataDirRaw/raw_525.root
0164 fi
0165
0166
0167 badChannelMap='../configs/TB2024/badChannelMap_TBSetup_CAEN_202408.txt'
0168 # # muon runs different scans in groups with separate pedestal
0169 #50.3 events
0170 muonHVScan_44V='305'
0171 pedHVScan_44V='303'
0172 if [ $2 == "muoncalibHV" ] || [ $2 == "muoncalibHV44" ]; then
0173         echo "running muon calib for 44V runs"
0174         MuonCalib $3 $pedHVScan_44V $muonHVScan_44V $dataDirRaw $dataDirOut muonHVScan_44V $badChannelMap $skipLayer
0175 fi
0176
0177 #50.1K events
0178 muonHVScan_43V='307'
0179 pedHVScan_43V='306'
0180 if [ $2 == "muoncalibHV" ] || [ $2 == "muoncalibHV43" ]; then
0181         echo "running muon calib for 43V runs"
0182         MuonCalib $3 $pedHVScan_43V $muonHVScan_43V $dataDirRaw $dataDirOut muonHVScan_43V $badChannelMap $skipLayer
0183 fi
0184
0185 #50.6K events
0186 muonHVScan_42V='309'
0187 pedHVScan_42V='308'
0188 if [ $2 == "muoncalibHV" ] || [ $2 == "muoncalibHV42" ]; then
0189         echo "running muon calib for 42V runs"
0190         MuonCalib $3 $pedHVScan_42V $muonHVScan_42V $dataDirRaw $dataDirOut muonHVScan_42V $badChannelMap $skipLayer
0191 fi
0192
0193 #51K events
0194 muonHVScan_41V='312'
0195 pedHVScan_41V='311'
0196 if [ $2 == "muoncalibHV" ] || [ $2 == "muoncalibHV41" ]; then
0197         echo "running muon calib for 41V runs"
0198         MuonCalib $3 $pedHVScan_41V $muonHVScan_41V $dataDirRaw $dataDirOut muonHVScan_41V $badChannelMap $skipLayer
0199 fi
0200
0201 #50.8K events
0202 muonHVScan_40V='316'
0203 pedHVScan_40V='315'
0204 if [ $2 == "muoncalibHV" ] || [ $2 == "muoncalibHV40" ]; then
0205         echo "running muon calib for 40V runs"
0206         MuonCalib $3 $pedHVScan_40V $muonHVScan_40V $dataDirRaw $dataDirOut muonHVScan_40V $badChannelMap $skipLayer
0207 fi
0208
0209 muonScanA_45V='244 250 282 283'
0210 pedScanA_45V='271 277'
0211 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibA1" ] || [ $2 == "muoncalib45V" ]; then
0212         # 192K events
0213         echo "running muon calib for 45V runs, campaing A1"
0214         MuonCalib $3 271 muonScanA1_45V $dataDirRaw $dataDirOut muonScanA1_45V $badChannelMap $skipLayer
0215 fi
0216 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibA2" ] || [ $2 == "muoncalib45V" ]; then
0217         # 201.6K events
0218         echo "running muon calib for 45V runs, campaing A2"
0219         MuonCalib $3 277 muonScanA2_45V $dataDirRaw $dataDirOut muonScanA2_45V $badChannelMap $skipLayer
0220 fi
0221
0222 #102.2K events
0223 muonScanD1_45V='412 417'                # these runs are messed up don't analyze
0224 pedScanD1_45V='420'
0225 # if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibD1" ]; then
0226 #       echo "running muon calib for 45V runs, campaing D1"
0227 # #     MuonCalib $3 $pedScanD1_45V muonScanD1_45V $dataDirRaw $dataDirOut muonScanD1_45V  $badChannelMap $skipLayer
0228 # #     MuonCalib $3 $pedScanD1_45V 412 $dataDirRaw $dataDirOut 412
0229 # #     MuonCalib $3 $pedScanD1_45V 417 $dataDirRaw $dataDirOut 417
0230 # fi
0231
0232 # 29.3K events
0233 muonScanD2_45V='460 456 457'
0234 pedScanD2_45V='454'
0235 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibD2" ] || [ $2 == "muoncalib45V" ]; then
0236         echo "running muon calib for 45V runs, campaing D2"
0237         MuonCalib $3 $pedScanD2_45V muonScanD2_45V $dataDirRaw $dataDirOut muonScanD2_45V $badChannelMap $skipLayer
0238 fi
0239
0240 # 50.6K events 2nd column underrespresented
0241 muonScanH1_45V='526 527'
0242 pedScanH1_45V='528'
0243 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibH1" ] || [ $2 == "muoncalib45V" ]; then
0244         echo "running muon calib for 45V runs, campaing H1"
0245         MuonCalib $3 $pedScanH1_45V muonScanH1_45V $dataDirRaw $dataDirOut muonScanH1_45V $badChannelMap $skipLayer
0246 fi
0247
0248 #33.5K events 1st-2nd column only
0249 muonScanH2_45V='554 559'
0250 # pedScanH2_45V='552 553'
0251 pedScanH2_45V='552'
0252 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibH2" ] || [ $2 == "muoncalib45V" ]; then
0253         echo "running muon calib for 45V runs, campaing H2"
0254         MuonCalib $3 $pedScanH2_45V muonScanH2_45V $dataDirRaw $dataDirOut muonScanH2_45V $badChannelMap $skipLayer
0255 fi
0256
0257 # 202.6K events
0258 muonScanB1_42V='331 322'
0259 pedScanB1_42V='332'
0260 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibB1" ] || [ $2 == "muoncalib42V" ]; then
0261         echo "running muon calib for 42V runs, campaing B1"
0262         MuonCalib $3 $pedScanB1_42V muonScanB1_42V $dataDirRaw $dataDirOut muonScanB1_42V $badChannelMap $skipLayer
0263 fi
0264
0265 # 214.8k events
0266 muonScanB2_42V='370 371 374'
0267 pedScanB2_42V='369'
0268 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibB2" ] || [ $2 == "muoncalib42V" ]; then
0269         echo "running muon calib for 42V runs, campaing B2"
0270         MuonCalib $3 $pedScanB2_42V muonScanB2_42V $dataDirRaw $dataDirOut muonScanB2_42V $badChannelMap $skipLayer
0271 fi
0272
0273 # 240.2K events
0274 muonScanC1_43_5V='376 375'
0275 pedScanC1_43_5V='377'
0276 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibC1" ] || [ $2 == "muoncalibAllC" ]; then
0277         echo "running muon calib for 43.5V runs, campaing C1"
0278         MuonCalib $3 $pedScanC1_43_5V muonScanC1_43_5V $dataDirRaw $dataDirOut muonScanC1_43_5V $badChannelMap $skipLayer
0279 fi
0280 muonScanC2_43_5V='405 410 408'
0281 pedScanC2_43_5V='404'
0282 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibC2" ] || [ $2 == "muoncalibAllC" ]; then
0283         echo "running muon calib for 43.5V runs, campaing C2"
0284         MuonCalib $3 $pedScanC2_43_5V muonScanC2_43_5V $dataDirRaw $dataDirOut muonScanC2_43_5V $badChannelMap $skipLayer
0285 fi
0286
0287 # 40.7K events
0288 muonScanE1_40V='463 464'
0289 pedScanE1_40V='465'
0290 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibE1" ] || [ $2 == "muoncalib40V" ] ; then
0291         echo "running muon calib for 40V runs, campaing E1"
0292         MuonCalib $3 $pedScanE1_40V muonScanE1_40V $dataDirRaw $dataDirOut muonScanE1_40V $badChannelMap $skipLayer
0293 fi
0294
0295 muonScanE2_40V='481 478'
0296 pedScanE2_40V='476'
0297 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibE2" ]  || [ $2 == "muoncalib40V" ]; then
0298         MuonCalib $3 $pedScanE2_40V muonScanE2_40V $dataDirRaw $dataDirOut muonScanE2_40V $badChannelMap $skipLayer
0299 fi
0300
0301 # 22.6K events
0302 muonScanF1_41V='486 489'
0303 pedScanF1_41V='492'
0304 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibF1" ] || [ $2 == "muoncalib41V" ]; then
0305         echo "running muon calib for 41V runs, campaing F1"
0306         MuonCalib $3 $pedScanF1_41V muonScanF1_41V $dataDirRaw $dataDirOut muonScanF1_41V $badChannelMap $skipLayer
0307 fi
0308
0309 # 42.8K events
0310 muonScanF2_41V='507 506'
0311 pedScanF2_41V='505'
0312 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibF2" ] || [ $2 == "muoncalib41V" ]; then
0313         echo "running muon calib for 41V runs, campaing F2"
0314         MuonCalib $3 $pedScanF2_41V muonScanF2_41V $dataDirRaw $dataDirOut muonScanF2_41V $badChannelMap $skipLayer
0315 fi
0316
0317 # 101.7K events
0318 muonScanG_46V='508 510 511 525'
0319 pedScanG_46V='521'
0320 if [ $2 == "muoncalibAll" ] || [ $2 == "muoncalibG" ] || [ $2 == "muoncalib46V" ]; then
0321         echo "running muon calib for 46V runs, campaing G"
0322         MuonCalib $3 $pedScanG_46V muonScanG_46V $dataDirRaw $dataDirOut muonScanG_46V $badChannelMap $skipLayer
0323 fi
0324
0325
0326
0327 if [ $2 == "reducemuons" ]; then
0328         runs='261 264 265 269 270 272 274 275 '
0329         for runNr in $runs; do
0330                 time ./DataPrep -f -d 1 -M -i $dataDirOutH/calibrated_Run_$runNr.root -o $dataDirOutH/calibrated_localmuononly_Run_$runNr.root
0331         done
0332         runs='251 252 254 257 258 '
0333         for runNr in $runs; do
0334                 time ./DataPrep -f -d 1 -M -i $dataDirOutE/calibrated_Run_$runNr.root -o $dataDirOutE/calibrated_localmuononly_Run_$runNr.root
0335         done
0336         time ./DataPrep -f -d 1 -M -i $dataDirOut/calibratedMuon_muonScanA1_45V.root -o $dataDirOut/calibratedMuon_localmuononly_muonScanA1_45V.root
0337         time ./DataPrep -f -d 1 -M -i $dataDirOut/calibratedMuon_muonScanA2_45V.root -o $dataDirOut/calibratedMuon_localmuononly_muonScanA2_45V.root
0338 fi