B-tagging efficiency

[DIBYASHREE SENGUPTA]

Hi,

I am trying to modify the b-tagging efficiency as follows :

if the jet is a c-jet or a light quark and gluon jet then it should be mistagged as a bjet with the following probability :

if pt <= 15.0 then probability = 0

if eta <=3.0 and pt > 15 and pt <=100 then probability = 1/150

if eta <=3.0 and pt > 100 and pt <=250 then probability = 1/(-0.667*pt+216.67)

if eta <=3.0 and pt > 250 then probability = 1/50

if eta > 3.0 then probability = 0

Following is the code I have used in the delphes_card.dat for implementing the above conditions :
###########
# b-tagging
###########

module BTagging BTagging {

set JetInputArray JetEnergyScale/jets

set BitNumber 0

# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
# PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
# gluon's PDG code has the lowest priority

# default efficiency formula (misidentification rate)
add EfficiencyFormula {0}{ (pt <= 15.0) * (0.00) +

(abs(eta) <= 3.0) * (pt > 15.0 && pt <= 100.0) * (1.0/150.0) +

(abs(eta) <= 3.0) * (pt > 100.0 && pt <= 250.0) *(1.0/(-0.667*pt+216.67)) +
(abs(eta) <= 3.0) * (pt > 250.0) * (1.0/50.0)+
(abs(eta) > 3.0) * (0.00)}

# efficiency formula for c-jets (misidentification rate)

add EfficiencyFormula {4}{ (pt <= 15.0) * (0.00)+

(abs(eta) <= 3.0) * (pt > 15.0 && pt <= 100.0) * (1.0/150.0) +

(abs(eta) <= 3.0) * (pt > 100.0 && pt <= 250.0) *(1.0/(-0.667*pt+216.67)) +
(abs(eta) <= 3.0) * (pt > 250.0) * (1.0/50.0)+
(abs(eta) > 3.0) * (0.00)}

# efficiency formula for b-jets
add EfficiencyFormula {5} { (pt <= 15.0) * (0.00)+

(abs(eta) <= 3.0) * (pt > 15.0) * (0.60) +

(abs(eta) > 3.0) * (0.00)}

}

However, when I am running delphes using this card it is exiting without running.

What should I do?

[Michele Selvaggi]

Can you please attach your card?
Michele

[DIBYASHREE SENGUPTA]

Replying to mselvaggi:

Can you please attach your card?
Michele

Hi,

I have attached the delphes_card.dat.

Looking forward for your reply.

Thanks,
Dibyashree

[Michele Selvaggi]

You were missing empty spaces. Find the corrected card attached.

[DIBYASHREE SENGUPTA]

Replying to mselvaggi:

You were missing empty spaces. Find the corrected card attached.

Thanks Michele . The code now works perfectly.

Regards,
Dibyashree