<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
<div><br>
</div>
<div>Hi Frank, </div>
<div><br>
</div>
<div>
<blockquote type="cite">
<div>Maybe I can ask a follow-up question... if we<br>
extend the cfs acceptance, this would also affect the nstable>2 cut<br>
(which would now also count particles in the extended acceptance).<br>
Would that be correct?</div>
</blockquote>
<br>
</div>
<div>Right. "nstable" part should be in |\eta| < 2.5 . So, if we extend the rapidity acceptance in the CFS projection, then we should apply a rapidity cut in the "nstable" part.. you know, something like below.</div>
<div><br>
</div>
<div>
<div> int n_stable = 0;</div>
<div> foreach (const Particle& p, cfs.particles()) {</div>
<div> const PdgId pid = abs(p.pdgId());</div>
<div><b><i> const double eta = fabs(p.momentum().eta());</i></b></div>
<div>
<div><b><i><span class="Apple-tab-span" style="white-space:pre"></span>if (inRange(eta, -2.5, 2.5)) continue;</i></b></div>
</div>
<div> if (pid == 11 || pid == 13 || pid == 211 || pid == 321 || pid == 2212) n_stable++;</div>
<div> }</div>
</div>
<div><br>
</div>
<div>or we can keep the CFS projection as it is and another projection can be applied for the trigger requirement: <span class="Apple-style-span" style="font-family: Menlo; font-size: 13px; "><b> "ChargedFinalState trigfs(-2.09, 3.86, 100*MeV);"</b></span></div>
<div><span class="Apple-style-span" style="font-family: Menlo; font-size: 13px; "><b><br>
</b></span></div>
<div>By the way, I don't know if we need 100*MeV for ATLAS MBTS requirement (this is applied in the current code) ? Probably, we don't need it but this should be checked by someone from ATLAS..</div>
<div><br>
</div>
<div>Cheers, </div>
<div>Sercan</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div></div>
<br>
<div>
<div>On Oct 23, 2012, at 11:15 AM, Frank Siegert wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<div>Hi Sercan,<br>
<br>
Thanks for noticing this inconsistency. I take it, you're familiar<br>
with the analysis? Maybe I can ask a follow-up question... if we<br>
extend the cfs acceptance, this would also affect the nstable>2 cut<br>
(which would now also count particles in the extended acceptance).<br>
Would that be correct? If so, I'll make the change in the cfs<br>
constructor directly.<br>
<br>
Cheers,<br>
Frank<br>
<br>
On 23 October 2012 11:08, Sercan Sen <<a href="mailto:Sercan.Sen@cern.ch">Sercan.Sen@cern.ch</a>> wrote:<br>
<blockquote type="cite"><br>
</blockquote>
<blockquote type="cite">Hello,<br>
</blockquote>
<blockquote type="cite"><br>
</blockquote>
<blockquote type="cite">I think there is a small correction needed on ATLAS_2011_I944826 routine.<br>
</blockquote>
<blockquote type="cite"><br>
</blockquote>
<blockquote type="cite">It applies charged final states projection for the particles in |\eta| <<br>
</blockquote>
<blockquote type="cite">2.5, however, search for the hits in MBTS rapidity region 2.09 < |\eta| <<br>
</blockquote>
<blockquote type="cite">3.84. see below.. In this case, trigger cut is limited to the<br>
</blockquote>
<blockquote type="cite">2.09<|\eta|<2.5.<br>
</blockquote>
<blockquote type="cite"><br>
</blockquote>
<blockquote type="cite"> ChargedFinalState cfs(-2.5, 2.5, 100*MeV);<br>
</blockquote>
<blockquote type="cite"> addProjection(cfs, "CFS");<br>
</blockquote>
<blockquote type="cite"><br>
</blockquote>
<blockquote type="cite"> const ChargedFinalState& cfs =<br>
</blockquote>
<blockquote type="cite">applyProjection<ChargedFinalState>(event, "CFS");<br>
</blockquote>
<blockquote type="cite"><br>
</blockquote>
<blockquote type="cite"> // ATLAS MBTS trigger requirement of at least one hit in either<br>
</blockquote>
<blockquote type="cite">hemisphere<br>
</blockquote>
<blockquote type="cite"> const ChargedFinalState& cfs =<br>
</blockquote>
<blockquote type="cite">applyProjection<ChargedFinalState>(event, "CFS");<br>
</blockquote>
<blockquote type="cite"> int n_mbts = 0;<br>
</blockquote>
<blockquote type="cite"> foreach (const Particle& p, cfs.particles()) {<br>
</blockquote>
<blockquote type="cite"> const double eta = fabs(p.momentum().eta());<br>
</blockquote>
<blockquote type="cite"> if (inRange(eta, 2.09, 3.84)) n_mbts++;<br>
</blockquote>
<blockquote type="cite"> }<br>
</blockquote>
<blockquote type="cite"><br>
</blockquote>
<blockquote type="cite"> if (n_mbts < 1) {<br>
</blockquote>
<blockquote type="cite"> MSG_DEBUG("Failed trigger cut");<br>
</blockquote>
<blockquote type="cite"> vetoEvent;<br>
</blockquote>
<blockquote type="cite"> }<br>
</blockquote>
<blockquote type="cite"><br>
</blockquote>
<blockquote type="cite">Cheers,<br>
</blockquote>
<blockquote type="cite">Sercan<br>
</blockquote>
<blockquote type="cite"><br>
</blockquote>
<blockquote type="cite">_______________________________________________<br>
</blockquote>
<blockquote type="cite">Rivet mailing list<br>
</blockquote>
<blockquote type="cite"><a href="mailto:Rivet@projects.hepforge.org">Rivet@projects.hepforge.org</a><br>
</blockquote>
<blockquote type="cite"><a href="http://www.hepforge.org/lists/listinfo/rivet">http://www.hepforge.org/lists/listinfo/rivet</a><br>
</blockquote>
</div>
</blockquote>
</div>
<br>
</body>
</html>