| Digitization and Reconstruction in Fairroot to new detectors [message #13673] |
Wed, 27 June 2012 19:16  |
Raghav Kunnawalkam
Messages: 63
Registered: February 2012
Location: Stony Brook, New York
|
continuous participant |
From: *physics.sunysb.edu
|
|
Hi All
I just understood today that if you have to do digitization to look at the detector response then you have to create a bunch of classes for its digitization and reconstruction. (like how it is given in pandaroot)
If we look at EMC in panda for example, it has a bunch of classes under EmcDigi, EmcReco, etc.. which i think are necessary for doing reconstruction.
I just followed the procedure given in the fairroot.gsi.de website for creating my own detector and i only have few classes like fairemca.cxx, fairemcacontfact.cxx, fairemcageo.cxx, fairemcageopar.cxx, fairemcapoint.cxx.
(if we say that i named my detector as Emca)
So if i need to do any sort of reconstruction on this particular detector which is actually my ElectroMagnetic CAlorimeter (EMCA), then i need all the classes present in /pandaroot/emc/ ?
I will work on getting the classes and renaming them and all that but i want to know if it is standard procedure.
Has anyone done reconstruction/ digitization from the very beginning so that they can tell me the steps to follow. I also know that i need to create a parameter file which has information for digi and reco. I went through some of the parameter files in /pandaroot/macro/params and i see that there are certain parameters for certain detectors and some others for other detectors. I want to know if there is like a manual of some sort telling me all the possible parameters and if there is like a default configuration that i can use for the time being.
Also i learnt from Florian today that i can do a "simple smearing" of the MCPoint, as a first step. But i dont know how to go about doing that in my digi macro. Is there an example or something like that available online in the repository?
Thanks a lot guys.
Cheers
Raghav
|
|
|
|
|
|
|
|
|
|
| Re: Digitization and Reconstruction in Fairroot to new detectors [message #13681 is a reply to message #13676] |
Thu, 28 June 2012 16:57  |
Raghav Kunnawalkam
Messages: 63
Registered: February 2012
Location: Stony Brook, New York
|
continuous participant |
From: *dyn.optonline.net
|
|
Hi Stefano
I am kinda confused as to where and how to define my read outs. For example say that i have a barrel tracker that is a cylinder from z = -10 to z = +10, with radius say 5 to 10 with all having the same units.
Now say that in the real world, my barrel tracker (Batr) has read outs every 1 unit in radius and 5 degree in theta and phi. So i have like a little section with
r = 5 to 6, and theta = 0 to 5 and phi = 0 to 5 inside which i cannot differentiate between hits.
I am very confused as to which class actually does my digitization. from my reading since yesterday and what you told i see that it is the hit along with hit producer that creates the data that i need.
I think that these are the steps that i need to follow, Please correct me if i am wrong.
1. First i create the following classes for my detector
FairBatr.cxx, FairBatrContFact.cxx, FairBatrGeo.cxx, FairBatrGeoPar.cxx, FairBatrPoint.cxx, FairBatrHit.cxx, FairBatrHitProducer.cxx, FairBatrDigi.cxx, FairBatrDigiProducer.cxx, ( if i want to do reco after this then i also add some more classes)
2. create the readouts in FairBatrHit, just like there are mdtdigibox and mdtdigistrip
3. Inside FairBatrHitProducer, in the Exec function do the actual digitization.
this is exactly what i am confused about.
4. pretty much after part 2 i am confused.
Thanks a lot for all your help
Cheers
Raghav
|
|
|
|
GSI Forum
Search
Help
Members
Login
Home
