Hello. The code looks just fine. If you want to load other calibrations from ascii to do BrTofTimeOffset calibration, just load those calibration module before this module (as SetLoadAscii(kTRUE) or something like that). Of course, don't forget to set SetSaveAscii(kTRUE) for this calibration module. Hiro Bekele, Selemon wrote: >Hi Peter, > > thanks for your useful comment. I have been revisiting the >delta delay and effective speed of light calibrations. I would like >to continue with time offset calibratios. I want to do it reading from >ascii files instead of reading from the database. But I see the following >piece of code in "BrTofTimeOffsetCalModule.cxx" > >//____________________________________________________________________ >void BrTofTimeOffsetCalModule::Event(BrEventNode* inNode, BrEventNode* outNode) >{ > // Per event method > // Note the algortimh actually used is different if the t0 comes from either > // BB vertex or from the pp trigger start counters. > // > SetState(kEvent); > > > // ---- get BB stuff > > if (fCommitAscii || fLoadAscii) > return; >.... > >} > >//____________________________________________________________________ >void BrTofTimeOffsetCalModule::Finish() >{ > //----------------------------------------------------------- > // Fit histos with gaus > //----------------------------------------------------------- > > // Job-level finalisation > SetState(kFinish); > > // if load ascii mode > if (fLoadAscii) > return; > > >.... > > >} > >This seems to me one cannot do time offset calibration by reading >other calibration parameters from ascii files. Is that true? Do I >neccessarily have to commit values to the database in order to do >calibrations? I am just being careful. > >any comments are welcome, > >thanks, > >Selemon, > >-----Original Message----- >From: Peter H.L. Christiansen [mailto:pchristi@nbi.dk] >Sent: Fri 9/16/2005 11:11 AM >To: Bekele, Selemon >Cc: Hans Hjersing Dalsgaard; brahms-dev-l >Subject: RE: [Brahms-dev-l] fstree2datamap > >Hi, > >Let me give a comment. sigmaang, sigmamult, sigmat are used to >parameterise the overall sigmam^2, see e.g. my Ph.d. appendix D (and >probably Djamels thesis as well). This also describes how to extract the >values by fitting the data.... > >Cheers, > Peter > >On Fri, 16 Sep 2005, Bekele, Selemon wrote: > > > >>Hi Hans, >> >> I found the following piece of code in >> >>fstree2datamap.cxx,v >> >> which I assume is the official version of the code in CVS from which the excutable >> >>"fstree2datamap" is made. >> >>//______________________________________________________________________ >>void SetupTof2Pid(TString sett) >>{ >> Float_t pi = mass[0]; >> Float_t p = mass[2]; >> >> Float_t sigmaAng = 0; >> Float_t sigmaMult = 0; >> Float_t sigmaTof = 0; >> >> if (sett == "4B608") { >> pi = 0.0073251; p = 0.854308; >> sigmaAng = 0.0001; sigmaMult = 0.0223507; sigmaTof = 0.0015291; >> } >> else if (sett == "4A608") { >> pi = 0.00868369; p = 0.858988; >> sigmaAng = 0.0001; sigmaMult = 0.0319085; sigmaTof = 0.00151653; >> } >> else if (sett == "6B861") { >> pi = 0.00694888; p = 0.849384; >> sigmaAng = 0.0001; sigmaMult = 0.0191514; sigmaTof = 0.00156334; >> } >> else if (sett == "6A861") { >> pi = 0.00694888; p = 0.849384; >> sigmaAng = 0.0001; sigmaMult = 0.0191514; sigmaTof = 0.00156334; >> } >> else if (sett == "8B861") { >> pi = 0.00467955; p = 0.849995; >> sigmaAng = 0.00818456; sigmaMult = 0.000100182; sigmaTof = 0.00153155; >> } >> else if (sett == "8A861") { >> pi = 0.00467955; p = 0.849995; >> sigmaAng = 0.00818456; sigmaMult = 0.000100182; sigmaTof = 0.00153155; >> } >> >> >> fPidPi2 ->SetParameters(sigmaAng, sigmaMult, sigmaTof, pi); >> fPidP2 ->SetParameters(sigmaAng, sigmaMult, sigmaTof, p); >> >> // diagnostic >> fPiLow2 ->SetParameters(sigmaAng, sigmaMult, sigmaTof, pi, -kMaxSigm2); >> fPiHigh2->SetParameters(sigmaAng, sigmaMult, sigmaTof, pi, kMaxSigm2); >> >> fPLow2 ->SetParameters(sigmaAng, sigmaMult, sigmaTof, p, -kMaxSigm2); >> fPHigh2->SetParameters(sigmaAng, sigmaMult, sigmaTof, p, kMaxSigm2); >> >>} >> >> >>################################################ >> >> >>My question is, how are the hardwired numbers for >> >> pi >> p >> sigmaAng >> sigmaMult >> sigmaTof >> >>are calculated? How does one get these values for other field settings? >> >>Thanks for the help, >> >>Selemon >> >> >>-----Original Message----- >>From: Hans Hjersing Dalsgaard [mailto:canute@fys.ku.dk] >>Sent: Thu 9/15/2005 3:41 PM >>To: Bekele, Selemon >>Subject: RE: [Brahms-dev-l] fstree2datamap >> >>Hi Selemon, >> >>The problem is that when Truls revised the code, he focused on making it >>compile. The TOF PID in FS is handled by three functions (and a fourth >>that sets up hardcoded values for H2) , 'TofPid', 'TofPid2' and 'RichPid' >>(I think). The RICH pid is not implemented and the H1 pid is not set up so >>the only thing 'working' at the moment is the H2 pid, handled by >>'Tof2Pid'. There is, however, a problem with this function, namely that it >>contains a cut in the delta(dy) which is quite reasonable but not working, >>so by removing this cut you will get not-empty histograms. At the moment I >>am trying to write some TOF pid in the FS that makes sense but progress is >>slow. If you can help me I would be grateful. I have no idea whatsoever >>about the RICH pid, as it is apparently not generated by 'dst2tree'... >>Regards, >> >>Hans >> >> >> >>>Hi Casper, Hans >>> >>> thanks for your prompt reply. Your suggestion worked. Now I am able to >>> make >>>a microDst for the fs. Please have a look at the plot >>> >>>http://www4.rcf.bnl.gov/~tigist/BbVsZdcVertex.gif >>> >>>The above plot tells me the microDst is OK. Now when I try to make an >>>"fstree2datamap", >>>all the histograms are empty. I tried all I can to understand what is >>>going on but have >>>not succeded yet. This is what I do: >>> >>>~/>fstree2datamap -a FS_accMap -c 0 -C 90 -b fstree -o fsdataH1_ -p H1 >>>-S 4A2442 >>> >>>which is all one has to do to get a data map. Do you guys have any idea >>>why >>>the histograms turn out to be empty? >>> >>>Thanks, >>> >>>-----Original Message----- >>>From: Casper Nygaard [mailto:caz@fys.ku.dk] >>>Sent: Thu 9/15/2005 2:11 AM >>>To: Bekele, Selemon >>>Subject: Re: [Brahms-dev-l] fstree2datamap >>> >>>The fstree2datamap doesnt search for the .root part of the name. You >>>can easily circumvent the problem by making a symbolic link to the >>>file; e.g. ln -s fstree4B442.root fstree4B442 >>> >>>This should help. >>> >>>caz >>> >>>"Bekele, Selemon" <bekeleku@ku.edu> said: >>> >>> >>> >>>>Hi, >>>> >>>> I was trying to make data maps using the excutable >>>> "fstree2datamap". >>>>I have files the files >>>> >>>> accMap4B2442.root >>>> fstree4B2442.root >>>> >>>>in the same directory from which I run "fstree2datamap" as shown >>>>below: >>>> >>>>#################################3 >>>> >>>>tigist@rcas0032:~/MyBrahms/brahms_app/tml_app/analyze_app/anatools> >>>> >>>> >>>fstree2datamap -b fstree -S 4B2442 -p H1 -c 0 -C 20 -a FS_accMap -o >>>fsdataH1_ >>> >>> >>>> AcceptanceMaps::LoadMaps : Loading maps from FS_accMap4B2442.root >>>> Defining 11 vertex bin limits: -50 -40 -30 -20 -10 0 10 20 30 40 50 >>>> ...loading maps for h, pi, k, p, d, . >>>> file : FS_accMap4B2442.root contains geometry for run 13147. >>>>----------------------------------------------------- >>>> Selecting pi,K,p, setting 4 B 2442 >>>>----------------------------------------------------- >>>> --> found [1] Tree files >>>>Error in <TFile::TFile>: file fstree4B2442 does not exist >>>>Error in <TFile::TFile>: file fstree4B2442 does not exist >>>>Error in <TFile::TFile>: file fstree4B2442 does not exist >>>>Error in <TFile::TFile>: file fstree4B2442 does not exist >>>> --- Vertex binning from map: >>>> Nbins : 10 >>>> Range : -50 | 50 >>>>Looping over 0 events... >>>> pion efficiency map: bin [140] >>>> kaon efficiency map: bin [140] >>>> proton efficiency map: bin [140] >>>> >>>>################################################## >>>> >>>> I expect the file "fstree4B2442.root" to be picked up by the >>>> program. >>>> >>>> >>>I wanted to look at >>> >>> >>>>the source code to see what I might need to modify in my input but >>>>could >>>> >>>> >>>not find it easily. >>> >>> >>>>Any ideas? >>>> >>>>Thanks, >>>> >>>>Selemon >>>> >>>> >>>>_______________________________________________ >>>>Brahms-dev-l mailing list >>>>Brahms-dev-l@lists.bnl.gov >>>>http://lists.bnl.gov/mailman/listinfo/brahms-dev-l >>>> >>>> >>>> >>> >>>-- >>>*************************************************************************** >>>--------------------------------- >>>/ Casper Nygaard, stud.scient / "In the end it is not the amount >>>of / Frederiksborgvej 54, 3.tv / years in your life, but the >>>amount / 2400 Copenhagen NV / of life in your years >>>that counts" / tlf#: (+45)27125518) / - >>>Abraham Lincoln --------------------------------- >>>*************************************************************************** >>> >>> >>> >>> >>>_______________________________________________ >>>Brahms-dev-l mailing list >>>Brahms-dev-l@lists.bnl.gov >>>http://lists.bnl.gov/mailman/listinfo/brahms-dev-l >>> >>> >> >> >> >>_______________________________________________ >>Brahms-dev-l mailing list >>Brahms-dev-l@lists.bnl.gov >>http://lists.bnl.gov/mailman/listinfo/brahms-dev-l >> >> >> > > > _______________________________________________ Brahms-dev-l mailing list Brahms-dev-l@lists.bnl.gov http://lists.bnl.gov/mailman/listinfo/brahms-dev-lReceived on Tue Sep 20 14:28:56 2005
This archive was generated by hypermail 2.1.8 : Tue Sep 20 2005 - 14:29:01 EDT
