1 # Input file
   2 
   3 #InputCollections = ["/afs/ifh.de/group/atlas/scratch/topdata/MC08/AODs/mc08.105200.T1_McAtNlo_Jimmy.merge.AOD.e357_s462_r635_t53_tid064824/AOD.064824._00203.pool.root.1"]
   4 #InputCollections = ["/afs/ifh.de/group/atlas/scratch/topdata/MC08/AODs/mc09_valid.105200.T1_McAtNlo_Jimmy.recon.AOD.e357_s462_s520_r757_tid080079/AOD.080079._000925.pool.root.1"]
   5 #InputCollections = ["/afs/ifh.de/group/atlas/scratch/topdata/MC08/AODs/mc08.107682.AlpgenJimmyWenuNp2_pt20.merge.AOD.e368_s462_r635_t53_tid068238/AOD.068238._00005.pool.root.1"]
   6 #InputCollections = ["/afs/ifh.de/group/atlas/scratch/topdata/cosmics/data09_cos.00121251.physics_CosmicMuons.merge.AOD.r733_p37_tid073935/AOD.073935._000001.pool.root.1"]
   7 
   8 InputCollections = ["AOD.pool.root"]
   9 
  10 
  11 # Output file
  12 
  13 #OutputFileName = "/usr1/scratch/userscratch/clange/SingleTopD3PD.root"
  14 #OutputFileName = "/usr1/scratch/userscratch/clange/SingleTopCosmicsD3PD.root"
  15 OutputFileName = "SingleTopD3PD.root"
  16 
  17 
  18 # Event numbers
  19 
  20 SkipEvents = 0
  21 EvtMax     = -1
  22 
  23 
  24 # Data type flags : RealData/Simulation and FullReco/FastSim
  25 
  26 Data     = False
  27 Cosmics = False
  28 FullReco = True
  29 
  30 
  31 # Execution options :
  32 
  33 # Event selections
  34 
  35 ExecuteTriggerDecisionFilter = False  #If set to True only events with passed chains including e20i or mu20 are selected
  36 ExecuteEventFilter           = False 
  37 
  38 # Reconstructed Object selections
  39 
  40 ExecuteOverlapRemoval   = False #have to use flag to remove overlaps
  41 ExecuteTrackSelection   = False
  42 ExecuteClusterSelection = False 
  43 ExecuteVertexSelection  = False   
  44 
  45 # Reconstructed Object associations
  46 
  47 ExecuteRecoAssociation = False
  48 
  49 # New Jet reconstructions
  50 
  51 ExecuteNewJetAlgorithm = True
  52 ExecuteNewJetSelection = True
  53 
  54 # New and existing Jet tagging
  55 
  56 ExecuteNewJetTagging        = True
  57 ExecuteJetReTagging         = False
  58 ExecuteReTaggedJetSelection = False
  59 
  60 # Trigger Object selections
  61 
  62 ExecuteTriggerSelection  = True #Object selection (without Decision filtering)
  63 ExecuteDecisionSelection = False #Object selection (with Decision filtering)
  64 
  65 # Reconstructed Object and Trigger Object associations
  66 
  67 ExecuteRecoToTriggerAssociation = False 
  68 ExecuteTriggerToRecoAssociation = False
  69 
  70 # Raw Object selections
  71 
  72 ExecuteCellSelection  = False
  73 ExecuteTowerSelection = False
  74 
  75 # Truth Particle/Object selections
  76 
  77 ExecuteTruthSelection      = True    
  78 ExecuteTruthOverlapRemoval = False   
  79 RemoveStatusThree          = True
  80 
  81 # Reconstructed Object and Truth Object/Particle associations
  82 
  83 ExecuteRecoToTruthAssociation = False   
  84 ExecuteTruthToRecoAssociation = False 
  85 
  86 # User Data calculations
  87 
  88 ExecuteEventVariables = False
  89 
  90 # Tree dumps
  91 
  92 ExecuteInfoTreeMaker       = True
  93 ExecuteTriggerTreeMaker    = True
  94 ExecuteJetTaggingTreeMaker = False
  95 ExecutePdfTreeMaker        = True
  96 SimulationVersion          = "mc08"
  97 
  98 # If ExecuteInfoTreeMaker = False, the EventInfo container is dumped in the RecoTree
  99 # If ExecuteInfoTreeMaker = True, the EventInfo container is dumped in a separate InfoTree
 100 # The Trigger containers are dumped in the TriggerTree only if the ExecuteSelection and ExecuteTriggerTreeMaker flags are set to True
 101 # The TrackParticle, CaloCluster and Vertex containers are dumped in the RecoTree only if the ExecuteSelection flags are set to True
 102 # The Cell and Tower containers are dumped in the RawTree only if the ExecuteSelection flags are set to True
 103 # If ExecutePdfTreeMaker = True, the Pdf info is dumped according to the SimulationVersion flag ("mc08" or "mc09")
 104 
 105 # Reconstructed Object container dumps
 106 
 107 DumpNewJets      = True
 108 DumpReTaggedJets = False
 109 
 110 # if DumpNewJets = True, the new jet containers are dumped in dedicated branches
 111 # if DumpReTaggedJets = True, the newly tagged jet containers are dumped in dedicated branches
 112 # To use one of these new containers as the current one, assign its name to the variable JetContainerName instead of the default
 113 # "Cone4H1TowerJets" (name with extension "AOD" if tagged). Note that in this case, the container will be dumped twice if the
 114 # dump flags are set to True
 115 
 116 # Truth Particle container dumps
 117 
 118 DumpLightQuarkTruth = False
 119 DumpPhotonTruth     = True
 120 
 121 
 122 # Input container names
 123     
 124 if Data:
 125   EventInfoName = "ByteStreamEventInfo"
 126   
 127 else :
 128   EventInfoName = "McEventInfo"
 129 
 130 if FullReco:
 131 
 132   MissingEtName              = "MET_RefFinal"
 133   PhotonContainerName        = "PhotonAODCollection"
 134   ElectronContainerName      = "ElectronAODCollection"
 135   MuonContainerName          = "StacoMuonCollection"
 136   TauJetContainerName        = "TauRecContainer"
 137   JetContainerName           = "Cone4H1TowerJets"
 138   TrackParticleContainerName = "TrackParticleCandidate"
 139   CaloClusterContainerName   = "egClusterCollection"
 140   VertexContainerName        = "VxPrimaryCandidate"
 141   EgammaDetailContainerName  = "egDetailAOD"
 142   TauJetDetailContainerName  = "TauRecDetailsContainer"
 143   
 144 else:
 145 
 146   MissingEtName              = "AtlfastMissingEt"
 147   PhotonContainerName        = "AtlfastPhotonCollection"
 148   ElectronContainerName      = "AtlfastElectronCollection"
 149   MuonContainerName          = "AtlfastMuonCollection"
 150   TauJetContainerName        = "AtlfastTauJet1p3pContainer"
 151   JetContainerName           = "AtlfastJetContainer"
 152   TrackParticleContainerName = "AtlfastTrackParticles"
 153   CaloClusterContainerName   = "egClusterCollection"
 154   VertexContainerName        = "VxPrimaryCandidate"
 155   EgammaDetailContainerName  = "egDetailAOD"
 156   TauJetDetailContainerName  = "TauRecDetailsContainer"
 157 
 158 if Cosmics:
 159   JetContainerName           = "AntiKt4H1TowerAODJets"
 160 
 161 TruthMissingEtName              = "MET_Truth"
 162 TruthParticleContainerName      = "SpclMC"
 163 TruthJetContainerName           = "Cone4TruthJets"
 164 TruthTrackParticleContainerName = "TrackParticleTruthCollection"
 165 McEventContainerName            = "GEN_AOD"
 166 
 167 CaloCellContainerName  = "AllCalo" 
 168 CaloTowerContainerName = "CombinedTower"
 169 
 170 
 171 # Reconstructed Object selection parameters
 172 
 173 if FullReco:
 174   #ElectronSelection_Author = "Electron"
 175   #ElectronSelection_EmMask = "Medium"
 176   ElectronSelection_Author = "All"
 177   ElectronSelection_EmMask = "None"
 178 else:
 179   ElectronSelection_Author = "All"
 180   ElectronSelection_EmMask = "None"
 181 
 182 ElectronSelection_ApplyEtaCrackCut = False
 183 ElectronSelection_ApplyEtConeCut   = False
 184 ElectronSelection_ConeSize         = 0.2 
 185 ElectronSelection_EtConeOption     = "OnlyEt"
 186 ElectronSelection_EtConeCut        = 6.0*GeV
 187 ElectronSelection_PtMin            = 5.0*GeV
 188 ElectronSelection_EtaMax           = 2.5
 189 ElectronSelection_EtaCrackMin      = 1.37
 190 ElectronSelection_EtaCrackMax      = 1.52
 191 
 192 MuonSelection_ApplyIsHighPtCut       = False
 193 MuonSelection_ApplyBestMatchCut      = False
 194 MuonSelection_ApplyIsCombinedMuonCut = False
 195 MuonSelection_Author                 = ["All"]
 196 MuonSelection_ApplyEtConeCut         = False
 197 MuonSelection_ConeSize               = 0.2
 198 MuonSelection_EtConeCut              = 6.0*GeV
 199 MuonSelection_PtMin                  = 5.0*GeV
 200 MuonSelection_EtaMax                 = 2.7
 201 
 202 if FullReco:
 203   #PhotonSelection_Author = "Photon"
 204   PhotonSelection_Author = "All"
 205 else:
 206   PhotonSelection_Author = "All"
 207 
 208 PhotonSelection_EmMask           = "None"
 209 PhotonSelection_ApplyEtaCrackCut = False
 210 PhotonSelection_ApplyEtConeCut   = False
 211 PhotonSelection_ConeSize         = 0.2
 212 PhotonSelection_EtConeCut        = 6.0*GeV
 213 PhotonSelection_PtMin            = 5.0*GeV
 214 PhotonSelection_EtaMax           = 2.5
 215 PhotonSelection_EtaCrackMin      = 1.37
 216 PhotonSelection_EtaCrackMax      = 1.52
 217  
 218 TauJetSelection_Author           = "All"
 219 TauJetSelection_TauFlag          = "Loose"
 220 TauJetSelection_TrackNumberMin   = 0
 221 TauJetSelection_TrackNumberMax   = 100
 222 TauJetSelection_DiscriminantName = []
 223 TauJetSelection_DiscriminantType = []
 224 TauJetSelection_DiscriminantCut  = []
 225 TauJetSelection_PtMin            = 10.0*GeV
 226 TauJetSelection_EtaMax           = 2.5
 227 
 228 JetSelection_PtMin  = 15.0*GeV
 229 JetSelection_EtaMax = 5.0
 230 
 231 TaggedJetSelection_TaggingAlgorithm = "Baseline" 
 232 TaggedJetSelection_TaggingWeightCut = 4.20
 233 
 234 TaggedJetSelection_TaggedPtMin    = 30.0*GeV
 235 TaggedJetSelection_TaggedEtaMax   = 2.5
 236 TaggedJetSelection_UntaggedPtMin  = 15.0*GeV
 237 TaggedJetSelection_UntaggedEtaMax = 5.0
 238 
 239 TrackPtMin = 0.5*GeV
 240 
 241 
 242 # Raw Object selection parameters
 243 
 244 CellSelection_EtMin = 1.0*GeV
 245 TowerSelection_EMin = 1.0*GeV
 246 
 247 
 248 # Detector description version
 249 
 250 DetDescrVersion = "ATLAS-GEO-02-01-00"
 251 
 252 
 253 # New Jet generation and selection parameters
 254 
 255 NewJetNumber = 0
 256 NewJetName   = []
 257 NewJetPtMin  = []
 258 NewJetEtaMax = []
 259 
 260 NewJetNumber += 1
 261 NewJetName   += ["AntiKt4LCTopo"]
 262 NewJetPtMin  += [15.0*GeV]
 263 NewJetEtaMax += [5.0]
 264 
 265 #NewJetNumber += 1
 266 #NewJetName   += ["Kt4LCTopo"]
 267 #NewJetPtMin  += [15.0*GeV]
 268 #NewJetEtaMax += [5.0]
 269 
 270 #NewJetNumber += 1
 271 #NewJetName   += ["SISCone4LCTopo"]
 272 #NewJetPtMin  += [15.0*GeV]
 273 #NewJetEtaMax += [5.0]
 274 
 275 
 276 # Existing Jet tagging and selection parameters
 277 
 278 ReTaggedJetNumber = 0
 279 ReTaggedJetName   = []
 280 ReTaggedJetPtMin  = []
 281 ReTaggedJetEtaMax = []
 282 
 283 ReTaggedJetNumber += 1
 284 ReTaggedJetName   += ["Cone4H1Topo"]
 285 ReTaggedJetPtMin  += [15.0*GeV]
 286 ReTaggedJetEtaMax += [5.0]
 287 
 288 ReTaggedJetNumber += 1
 289 ReTaggedJetName   += ["Cone7H1Tower"]
 290 ReTaggedJetPtMin  += [15.0*GeV]
 291 ReTaggedJetEtaMax += [5.0]
 292 
 293 
 294 # Event selection parameters
 295 
 296 MissingEtMin       = 0.0*GeV
 297 
 298 LeptonNumberMin    = 0
 299 LeptonPtMin        = 30.0*GeV
 300 
 301 JetNumberMin       = 0
 302 JetPtMin           = 15.0*GeV
 303 
 304 TaggedJetNumberMin = 0
 305 TaggedJetPtMin     = 30.0*GeV
 306 
 307 
 308 # Sequence settings
 309 
 310 include("DESYsetSingleTopDPDMaker.py")  
 311 
 312   
 313 # Sequence execution
 314 			
 315 include("DESYexeSingleTopDPDMaker.py")