ls -1 /nfs/slac/g/lcd/ilc_data4/snowmass/ILC350/flavorTraining/sidloi3/slcio/reco/*.slcio > ilc350flav.lst
source subdstlst-tup-from-rec-job ilc350flav.lst
WHERE:
subdstlst-tup-from-rec-job:
cat $1 | xargs -n 1 bsub -q long -R rhel50 "echo sleep 0.1;source subflavtup-from-rec-job.sh"
subflavtup-from-rec-job.sh: export mydir=`echo $1 | awk -F '/' '{print substr($0,1,length($0)-length($NF))}'`
export myfil=`echo $1 | awk -F '/' '{print $NF}'`
echo "mydir = "$mydir
echo "myfil = "$myfil
mkdir /u/ey/homer/sidhome/lcfi/steering/tupparts/$myfil-dir
cd /u/ey/homer/sidhome/lcfi/steering/tupparts/$myfil-dir
rm /u/ey/homer/sidhome/lcfi/steering/tupparts/$myfil-dir/*.{slcio,log}
ln -s $1 input-rec.slcio
Marlin /u/ey/homer/sidhome/lcfi/steering/vertexing-basic.xml > /u/ey/homer/sidhome/lcfi/steering/tupparts/$myfil-dir/pretup-vtx.log
Marlin /u/ey/homer/sidhome/lcfi/steering/makentuple-all-batch.xml > /u/ey/homer/sidhome/lcfi/steering/tupparts/$myfil-dir/tup.log
rm input-dst.slcio
ln -s output-ntpl.root $myfil-ntpl.root
The XML's accessed above are shown below:
vertexing-basic.xml: <marlin> <execute> <processor name="VertexFinder"/> <processor name="MyLCIOOutputProcessor"/> </execute> <global> <parameter name="LCIOInputFiles"> input-rec.slcio </parameter> <parameter name="GearXMLFile">/u/ey/homer/sidhome/lcfi/steering/clic_sid.gear</parameter> <parameter name="MaxRecordNumber" value="-1" /> <parameter name="SkipNEvents" value="0" /> <parameter name="SupressCheck" value="false" /> <parameter name="Verbosity" options="DEBUG0-4,MESSAGE0-4,WARNING0-4,ERROR0-4,SILENT">WARNING</parameter> </global> <processor name="VertexFinder" type="LcfiplusProcessor"> <!-- run primary and secondary vertex finders --> <parameter name="Algorithms" type="stringVec"> PrimaryVertexFinder BuildUpVertex </parameter> <parameter name="ReadSubdetectorEnergies" type="int" value="0"/> <!-- true for ILD --> <parameter name="UpdateVertexRPDaughters" type="int" value="0"/> <!-- false for non-updative PandoraPFOs --> <parameter name="PrintEventNumber" type="int" value="10"/> <!-- 0 for not printing event number, n for printing every n events --> <!-- specify input collection names --> <parameter name="PFOCollection" type="string" value="PandoraPFOCollection" /> <parameter name="PrimaryVertexCollectionName" type="string" value="PrimaryVertex" /> <parameter name="BuildUpVertexCollectionName" type="string" value="BuildUpVertex" /> <parameter name="BuildUpVertex.V0VertexCollectionName" type="string" value="BuildUpVertex_V0" /> <parameter name="MagneticField" type="float" value="5"/> <parameter name="BeamSizeX" type="float" value="335E-6"/> <parameter name="BeamSizeY" type="float" value="2.7E-6"/> <parameter name="BeamSizeZ" type="float" value="0.225"/> <!-- parameters for primary vertex finder --> <parameter name="PrimaryVertexFinder.TrackMaxD0" type="double" value="20." /> <parameter name="PrimaryVertexFinder.TrackMaxZ0" type="double" value="20." /> <parameter name="PrimaryVertexFinder.TrackMaxInnermostHitRadius" type="double" value="20." /> <parameter name="PrimaryVertexFinder.TrackMinVtxFtdHits" type="int" value="3" /> <parameter name="PrimaryVertexFinder.Chi2Threshold" type="double" value="25." /> <parameter name="PrimaryVertexFinder.UseBeamConstraint" type="int" value="1" /> <!-- parameters for secondary vertex finder --> <parameter name="BuildUpVertex.TrackMaxD0" type="double" value="10." /> <parameter name="BuildUpVertex.TrackMaxZ0" type="double" value="20." /> <parameter name="BuildUpVertex.TrackMinPt" type="double" value="0.1" /> <parameter name="BuildUpVertex.TrackMaxD0Err" type="double" value="0.1" /> <parameter name="BuildUpVertex.TrackMaxZ0Err" type="double" value="0.1" /> <parameter name="BuildUpVertex.TrackMinTpcHits" type="int" value="4" /> <parameter name="BuildUpVertex.TrackMinFtdHits" type="int" value="3" /> <parameter name="BuildUpVertex.TrackMinVxdHits" type="int" value="3" /> <parameter name="BuildUpVertex.TrackMinVxdFtdHits" type="int" value="0" /> <parameter name="BuildUpVertex.PrimaryChi2Threshold" type="double" value="25." /> <parameter name="BuildUpVertex.SecondaryChi2Threshold" type="double" value="9." /> <parameter name="BuildUpVertex.MassThreshold" type="double" value="10." /> <parameter name="BuildUpVertex.MinDistFromIP" type="double" value="0.3" /> <parameter name="BuildUpVertex.MaxChi2ForDistOrder" type="double" value="1.0" /> <parameter name="BuildUpVertex.AssocIPTracks" type="int" value="1" /> <parameter name="BuildUpVertex.AssocIPTracksMinDist" type="double" value="0." /> <parameter name="BuildUpVertex.AssocIPTracksChi2RatioSecToPri" type="double" value="2.0" /> <parameter name="BuildUpVertex.UseV0Selection" type="int" value="1" /> </processor> <processor name="MyLCIOOutputProcessor" type="LCIOOutputProcessor"> <parameter name="LCIOOutputFile" type="string"> input-dst.slcio </parameter> <parameter name="LCIOWriteMode" type="string" value="WRITE_NEW"/> </processor> </marlin>
makentuple-all-batch.xml: <marlin> <execute> <processor name="JetClustering"/> <processor name="MakeNtuple"/> </execute> <global> <parameter name="LCIOInputFiles"> input-dst.slcio </parameter> <parameter name="GearXMLFile">/u/ey/homer/sidhome/lcfi/steering/clic_sid.gear</parameter> <parameter name="MaxRecordNumber" value="-1" /> <parameter name="SkipNEvents" value="0" /> <parameter name="SupressCheck" value="false" /> <parameter name="Verbosity" options="DEBUG0-4,MESSAGE0-4,WARNING0-4,ERROR0-4,SILENT">WARNING</parameter> </global> <processor name="JetClustering" type="LcfiplusProcessor"> <!-- run primary and secondary vertex finders --> <parameter name="Algorithms" type="stringVec"> JetClustering JetVertexRefiner</parameter> <!-- general parameters --> <parameter name="PFOCollection" type="string" value="PandoraPFOCollection" /> <!-- input PFO collection --> <parameter name="UseMCP" type="int" value="0" /> <!-- MC info not used --> <parameter name="MCPCollection" type="string" value="" /> <!-- not used --> <parameter name="MCPFORelation" type="string" value="" /> <!-- not used --> <parameter name="ReadSubdetectorEnergies" type="int" value="0"/> <!-- true for ILD --> <parameter name="UpdateVertexRPDaughters" type="int" value="0"/> <!-- false for non-updative PandoraPFOs --> <parameter name="MagneticField" type="float" value="5"/> <!-- jet clustering parameters --> <parameter name="JetClustering.InputVertexCollectionName" type="string" value="BuildUpVertex" /> <!-- vertex collections to be used in JC --> <parameter name="JetClustering.OutputJetCollectionName" type="stringVec" value="VertexJets" /> <!-- output collection name, may be multiple --> <parameter name="JetClustering.NJetsRequested" type="intVec" value="2" /> <!-- Multiple NJets can be specified --> <parameter name="JetClustering.YCut" type="doubleVec" value="0." /> <!-- specify 0 if not used --> <parameter name="JetClustering.UseMuonID" type="int" value="1" /> <!-- jet-muon ID for jet clustering --> <parameter name="JetClustering.VertexSelectionMinimumDistance" type="double" value="0.3" /> <!-- in mm --> <parameter name="JetClustering.VertexSelectionMaximumDistance" type="double" value="30." /> <!-- in mm --> <parameter name="JetClustering.VertexSelectionK0MassWidth" type="double" value="0.02" /> <!-- in GeV --> <parameter name="JetClustering.YAddedForJetVertexVertex" type="double" value="100"/> <!-- add penalty for combining vertices --> <parameter name="JetClustering.YAddedForJetLeptonVertex" type="double" value="100"/> <!-- add penalty for combining lepton and vertex --> <parameter name="JetClustering.YAddedForJetLeptonLepton" type="double" value="100"/> <!-- add penalty for combining leptons --> <!-- vertex refiner parameters --> <parameter name="JetVertexRefiner.InputJetCollectionName" type="string" value="VertexJets" /> <parameter name="JetVertexRefiner.OutputJetCollectionName" type="string" value="RefinedJets" /> <parameter name="JetVertexRefiner.PrimaryVertexCollectionName" type="string" value="PrimaryVertex" /> <parameter name="JetVertexRefiner.InputVertexCollectionName" type="string" value="BuildUpVertex" /> <parameter name="JetVertexRefiner.V0VertexCollectionName" type="string" value="BuildUpVertex_V0" /> <parameter name="JetVertexRefiner.OutputVertexCollectionName" type="string" value="RefinedVertex" /> <parameter name="JetVertexRefiner.MinPosSingle" type="double" value="0.3" /> <parameter name="JetVertexRefiner.MaxPosSingle" type="double" value="30." /> <parameter name="JetVertexRefiner.MinEnergySingle" type="double" value="1." /> <parameter name="JetVertexRefiner.MaxAngleSingle" type="double" value="0.5" /> <parameter name="JetVertexRefiner.MaxSeparationPerPosSingle" type="double" value="0.1" /> <parameter name="JetVertexRefiner.mind0sigSingle" type="double" value="5." /> <parameter name="JetVertexRefiner.minz0sigSingle" type="double" value="5." /> <parameter name="JetVertexRefiner.OneVertexProbThreshold" type="double" value="0.001" /> <parameter name="JetVertexRefiner.MaxCharmFlightLengthPerJetEnergy" type="double" value="0.1" /> </processor> <processor name="MakeNtuple" type="LcfiplusProcessor"> <parameter name="Algorithms" type="stringVec">FlavorTag MakeNtuple</parameter> <!-- general parameters --> <parameter name="PFOCollection" type="string" value="PandoraPFOCollection" /> <!-- input PFO collection --> <parameter name="UseMCP" type="int" value="0" /> <!-- MC info not used --> <parameter name="MCPCollection" type="string" value="" /> <!-- not used --> <parameter name="MCPFORelation" type="string" value="" /> <!-- not used --> <parameter name="ReadSubdetectorEnergies" type="int" value="0"/> <!-- true for ILD --> <parameter name="UpdateVertexRPDaughters" type="int" value="0"/> <!-- false for non-updative PandoraPFOs --> <parameter name="PrimaryVertexCollectionName" type="string" value="PrimaryVertex" /> <parameter name="FlavorTag.JetCollectionName" type="string" value="RefinedJets" /> <parameter name="MakeNtuple.AuxiliaryInfo" type="int" value="-1" /> <parameter name="FlavorTag.D0ProbFileName" type="string" value="/u/ey/homer/sidhome/lcfi/steering/d0prob_zpole.root"/> <parameter name="FlavorTag.Z0ProbFileName" type="string" value="/u/ey/homer/sidhome/lcfi/steering/z0prob_zpole.root"/> <parameter name="MakeNtuple.OutputRootFileName" type="string" value="output-ntpl.root" /> </processor> </marlin>