ILD card not working

[Zhao-Huan Yu]

Hi,

I was using the recently added card for the ILD detector (Delphes 3.3.0, delphes_card_ILD.tcl) in MadGraph5 2.3.0, but the resulting lhco file contained no event. It seemed that there were errors when the output root file was converted to the lhco file. I got the following error message in the delphes.log file:

Reading delphes.root
Input file contains 1000 events
WARNING: cannot access branch 'Track', return NULL pointer
WARNING: cannot access branch 'Tower', return NULL pointer
ERROR: ROOT file doesn't contain all required branches

Actually, I could not find any Track and Tower branches in the root file.

[Michele Selvaggi]

Hi,

the fact that you have a root file means that the card is working properly. The error you have is because the root2lhco converter is not adapted to convert the delphes output from any card.

You can either:

• analyse directly the root file
• add in the ILD card, in the TreeWriter section, the following lines:

add Branch AngularSmearing/tracks Track Track
add Branch TowerMerger/towers Tower Tower

Michele

[Zhao-Huan Yu]

Hi Michele,

Thanks for your reply! I added the lines in the TreeWriter section as you said:

module TreeWriter TreeWriter {
#  add Branch InputArray BranchName BranchClass
  add Branch Delphes/allParticles Particle GenParticle
  add Branch GenJetFinder/jets GenJet Jet

  add Branch AngularSmearing/tracks Track Track
  add Branch TowerMerger/towers Tower Tower

  add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
  add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
  add Branch ECal/eflowPhotons Photon Photon

  add Branch ElectronFilter/electrons Electron Electron
  add Branch MuonMomentumSmearing/muons Muon Muon
  add Branch JetEnergyScale/jets Jet Jet
  add Branch MissingET/momentum MissingET MissingET
  add Branch ScalarHT/energy ScalarHT ScalarHT
}

Then for the e+ e- -> e+ e- process, the resulting lhco file contained events like this:

   #  typ      eta      phi      pt    jmas   ntrk   btag  had/em   dum1   dum2
   0             1        0
   1    1   -2.199    0.313   54.73    0.00   -1.0    0.0    0.00    0.0    0.0
   2    1    2.199   -2.829   54.12    0.00    1.0    0.0    0.00    0.0    0.0
   3    4   -2.199    0.313   54.73    0.00    1.0    0.0  999.90    0.0    0.0
   4    4    2.199   -2.829   54.12    0.00    1.0    0.0  999.90    0.0    0.0
   5    6    0.000   -2.829    0.61    0.00    0.0    0.0    0.00    0.0    0.0

It seemed that both the e- and the e+ appeared twice, and falsely recognized as jets in the second time.

[Michele Selvaggi]

Hi,

I have attached a new card. can you please try with that?
Be aware that in order to treat properly the potential overlap between leptons and jets, one has to define Isolation criteria. I have taken them to be the same as in CMS, which is completely arbitrary. If you want different criteria you'll have to change accordingly the Isolation part of the card.

Michele

[Zhao-Huan Yu]

Hi Michele,

Now it seems that the new card has a problem with the muon finder. I generate e+ e- > mu+ mu- and the results look like this:

   #  typ      eta      phi      pt    jmas   ntrk   btag  had/em   dum1   dum2
   0             1        0
   1    4   -0.161    0.114  246.60    0.00    1.0    0.0  999.90    0.0    0.0
   2    4    0.161   -3.028  246.33    0.00    1.0    0.0  999.90    0.0    0.0
   3    6    0.000   -3.028    0.26    0.00    0.0    0.0    0.00    0.0    0.0
   0             2        0
   1    4    0.738   -0.117  194.44   -0.00    1.0    0.0  999.90    0.0    0.0
   2    4   -0.738    3.025  193.62   -0.00    1.0    0.0  999.90    0.0    0.0
   3    6    0.000    3.025    0.82    0.00    0.0    0.0    0.00    0.0    0.0
   0             3        0
   1    4   -0.998   -2.985  161.91    0.00    1.0    0.0  999.90    0.0    0.0
   2    4    0.998    0.156  161.72   -0.00    1.0    0.0  999.90    0.0    0.0
   3    6    0.000    0.156    0.19    0.00    0.0    0.0    0.00    0.0    0.0

[Michele Selvaggi]

Hi,

it should be fixed now,

Michele

[Zhao-Huan Yu]

Hi Michele,

The v3 card works well after I add the following line in the TreeWriter section:

add Branch TrackMerger/tracks Track Track

Thank you very much!