Dear Torben,
Thanks for the answer.
Regards,
Xinchun
On 5/4/13 8:45 AM, "Torben Ferber (DESY)" <[log in to unmask]> wrote:
>Dear Xinchun,
>
>formation zone reweighting is only applied to DIS events in GENIE. If
>your events are exclusive QEL, the weight=1 assignment is correct since
>the formation zone model is not used for QEL events. You may test this
>by generating DIS events and reweight them using your analysis code and
>you will see weights!=1.
>
>best regards,
>Torben
>
>--
>Torben Ferber
>
>Deutsches Elektronen Synchrotron (DESY)
>Notkestraße 85, Bld. 1d, O2/412
>22607 Hamburg
>Germany
>
>On 05/04/2013 02:24 PM, Xinchun Tian wrote:
>> HI Costas,
>>
>> I tried grwght1scan, with the following command, but when I check with
>>the
>> output root file, all the weights array are 1.
>> shell>grwght1scan -f ../ccqe/gntp.21211200.ghep.root -n 5000 -s FormZone
>> -t 3
>> While
>> shell>grwght1scan -f ../ccqe/gntp.21211200.ghep.root -n 5000 -s MaCCQE
>>-t 3
>> Do give me none 1 weights.
>> The above behavior is consistent with my observation using the method
>> introduced by your 2009 tutorial.
>>
>> Is this expected or a bug?
>>
>> I am using GENIE v3665c, associated with NOvA software.
>>
>> Thanks,
>> Xinchun
>>
>> On 5/3/13 8:54 AM, "[log in to unmask]"
>> <[log in to unmask]> wrote:
>>
>>> Dear Xinchun.
>>>
>>> it seems to me that you don't use the reweighting package correctly.
>>> Please have a look at the example application shown here:
>>>
>>>https://genie.hepforge.org/trac/browser/trunk/src/support/rwght/gRwght1S
>>>ca
>>> n.cxx
>>> cheers
>>> Costas
>>>
>>> On 3 May 2013, at 03:30, Xinchun Tian wrote:
>>>
>>> Dear experts,
>>>
>>> I have couple of questions on GENIE Reweight.
>>>
>>> 1) From UserPhysicsOptions.xml, I can find the default values of
>>> different parameters, but where can I find the sigmas?
>>>
>>> 2) Formation zone parameter reweight does not work. I tried different
>>> values of QEL-Ma to make sure my code works, but with different
>>> FZONE-ct0, the reweight factor is always 1, am I missing something? I
>>> attach piece of relevant code for your reference.
>>> <!--
>>>
>>>
>>>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>>>~~
>>> ~~~~~~~~~~~~~~~~~~~~~~~
>>> Parameters relevant to formation zone simulation
>>> - ct0 is the formation time times the speed of light (given in fm)
>>> - KPt2 is the parameter multiplying pT2 in formation zone calc.
>>> -->
>>> <param type="double" name="FZONE-ct0"> 0.342 </param>
>>> <param type="double" name="FZONE-KPt2"> 0.0 </param>
>>>
>>> 3) What is the right way to apply this reweight factor to the event to
>>> get the reweight event?
>>>
>>> Thanks,
>>> Xinchun
>>> University of South carolina
>>>
>>> // access the algorithm factory and the configuration pool
>>> AlgFactory * algf = AlgFactory::Instance();
>>> AlgConfigPool * algc = AlgConfigPool::Instance();
>>>
>>> // get the physics parameters that GENIE developers find meaningfull
>>> for a user to tweak
>>> Registry * user_conf = algc->GlobalParameterList();
>>> user_conf->UnLock();
>>> user_conf->Set("FZONE-ct0", 2.0);
>>> // user_conf->Set("QEL-Ma", 0.99);
>>> algf->ForceReconfiguration();
>>> LOG("gRewght", pINFO) << "User options / current " << *user_conf;
>>>
>>> // create a weight calculator
>>> GReWeightNuXSecHelper wcalc;
>>>
>>> // Event loop
>>> for(Long64_t iev = 0; iev < nmax; iev++) {
>>> er_tree->GetEntry(iev);
>>>
>>> NtpMCRecHeader rec_header = mcrec->hdr;
>>> EventRecord & event = *(mcrec->event);
>>>
>>> double wght = wcalc.NewWeight(event);
>>> std::cout<<wght<<std::endl;
>>> }
>>>
>>> --
>>> Scanned by iCritical.
>
|