FAMOS PU
summary at 15.09.2005 can be found at fpu.ppt
FAMOS_1_3_2_pre2
V.Zhukov 12.09.2005
version tested with FAMOS_1_3_2_pre2 where the structure of the
reader was changed, in
addition it operates with scramv1 which also impose some modifications
in BuildFiles.
Version is available at cern /afs/cern.ch/user/j/joukov/public/FAMOS_1_3_2_pre2
or tar file of src/ famos132pre2.tar.gz
Modifications are analogous to previous version:
locally checkout COBRA_8_5_1
( cvs -d
:kserver:cmscvs.cern.ch:/cvs_server/repositories/COBRA co -r
COBRA_8_5_1 xxx):
Utilities/CHBook4
CHBookFile
add folders (SIG and PU)
and public cd() to navigate in different folder
CHBookCWNtuple
add offset to allow the ntuples with the same Id
GeneratorInterface/HBookInterface
HepEventCmkinNtupleReader add files->cd() to
select the ntuple file for signal and PU.
Famos_1_3_2_pre2:
FamosGeneric/FamosManager
FamosPUGenerator
analogous to previous version, reads the PU ntuples.new: the ntuple Id
can be defined, multifiles.
FamosEvtMgr
modification for the PU generator, now the PU ntuples is read
first
FamosGeneric/FamosSimEven
FBaseSimEvent
addPU() method which merges sig and PU HEP events
FamosGeneric/FamosManager/test
CmkinEventReader
add folder name, some extra prints in debug mode
testFamosWithPU.cpp,
testCmkinReader.cpp, famosrc test files(no
RecCollection in these examples)
The orcarc interface remains the same (in addition the ntuplId can be
specified for PU ntuples), see below
Note that in order to run you have to cleanup your previous libs
and bins by scramv1 b clean, then
scramv1 b all new libs and binaries. Mofify the famosrc file for you
local file locations (CMKIN and PUCollection)
and run testFamosWithPU -c famosrc. The root out.root will be created.
The official minbias ntuples can be accessed from cern via rfcp
(rfdir) at
lxcmsa:/data1/cmsprod/generation/mu05b_MBforPU/mu05b_MBforPU_20200000*.ntpl.gz
(each 1000 events)
FAMOS_1_2_0
V.Zhukov 15.06.2005
Modifications:
- CHBook4 checkout locally at Famos and modified.
This allowed to get rid of
FamosHepEventCmkinNtupleReader
FamosCHBookFile
- the PollCollection can contain list of files
PUGenerator:PUCollection
=
/portal/ekpcms1/scratch/zhukov/gendata/minbias1.ntpl
/portal/ekpcms1/scratch/zhukov/gendata/minbias2.ntpl
/portal/ekpcms1/scratch/zhukov/gendata/minbias3.ntpl
...
modified version can be found at /afs/cern.ch/user/j/joukov/public/
FAMOS_1_2_0
Zhukov 02.06.2005
1. Status
First simple version of the minimbias event pile up is implimented in FAMOS_1_2_0
Can be found at /afs/cern.ch/user/j/joukov/public/
FAMOS_1_2_0
or as a tar file famos.tar.gz
plans:
version with simpler structure
is under way
to do:
-should be impimented in the head version.
- detector PU ? who and how to do that?
Main problem:
COBRA hbook interface has to be slightly modified or the
local FAMOS hbook should be used.
2. Implimentation
Modifications and new
classes:
FamosGeneric/
FamosEventReader/
FamosPUGenerator
(reads PU events and aad them to existing sim event)
FamosHepEventCmkinNtupleReader (just extension of the
HepEventCmkinNtupleReader)
FamosCHBookFile (add cd() to folder)
CmkinEventReader
FamosManager/
FamosEventMgr (acreate PU generator inside)
test/
testFamosWithPU ( test )
FamosSimEvent/
FBaseSimEvent ( add PU events to
the sim event)
Note that because hbook wrapper in COBRA is almost all private ,
most of the basic hbook classes
are just copied to
FamosEventReader/
CHObject
CHBook
CHBookCWNtuple
or reimplimented (as FamosCHBookFile)
. To use the includes the scram has
to be reconfigured
otherwise only src is used.
All this can be avoided if we will make changes directly in COBRA
and
synchronize with the FAMOS or rename and reimpliment our
hbook version.
What actually is needed from hbook :
-larger pawc size (defined in CHObject)
-public cd() function in CHFile
- or/and public access to the hbook offset in hropen in
CHBookCWNtuple
(Another rpossibility different ntpl ID for minbias. This can simplify
two hbook files simulataneous reading)
3. Running
the settings config is similar to the PUGenerator in ORCA
in the orcarc :
# name of the minibias ntpl, the id =101 is fixed, if ntpl doesnt exist
# no PU event are added, if error during redout ->exit
PUGenerator:PUCollection
=
/portal/ekpcms1/scratch/zhukov/gendata/minbias1.ntpl
/portal/ekpcms1/scratch/zhukov/gendata/minbias2.ntpl
# the pileup event is selected randomly form the whole PU file
PUGenerator:RandomShuffle = false
# dont reuse the event. If there are no event left the
file is starting from the beginning or
# in case of random access the reuse is disabled when all events are
used
PUGenerator:NoReuse
= false
# poisson average
PUGenerator:AverageEvents = 3.5
#minimum and maximum bunch works as a scale to the PU events(at that
moment).
# i Npuevents=(MaxBunch-MinBunch)*Poisson(average)
PUGenerator:MinBunch
= 0
PUGenerator:MaxBunch
= 1
# defines slice of the minibias file to be used
PUGenerator:FirstEvent
= 1
PUGenerator:LastEvent
= 999999
# extra prints
PUGenerator:debug
= true
# PU ntuple ID (for FAMOS_1_3_2)
PUGenerator:ntplId
= 101
The example how to run and access information can be
found at FamosManager/test/testFamosWithPU.cpp
4. Examples
signal: SUSY LM9 point with trileptonic final state;
pp->ch20+chi+-->3l+chi0+nu
minibias: CMKIN minbias.txt
two PU considered low luminosity Npu=3.5 and high luminosity 17.5.
Reconstruction was done in FAMOS_1_2_0 with simple test testFamosWithPU.cpp
Number of PU events in each signal event
Number of reconstructed jets
Missing ET
Pt of all reconstructed muons (no selection)
Pt of reconstructed electrons
number of reconstructed tracks
