GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » General » Digitization and Reconstruction in Fairroot to new detectors
Digitization and Reconstruction in Fairroot to new detectors [message #13673] Wed, 27 June 2012 19:16 Go to next message
Raghav Kunnawalkam is currently offline  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 #13674 is a reply to message #13673] Wed, 27 June 2012 19:24 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 8.28.178*
You need only one data object and one task. I suggest to look at the emc hit and the emc hit producer. All the other classes are not needed.
Re: Digitization and Reconstruction in Fairroot to new detectors [message #13675 is a reply to message #13674] Wed, 27 June 2012 19:53 Go to previous messageGo to next message
Raghav Kunnawalkam is currently offline  Raghav Kunnawalkam
Messages: 63
Registered: February 2012
Location: Stony Brook, New York
continuous participant
From: *physics.sunysb.edu
Hi Stefano

Ok, so the hit producer calculates the hit for my detector and converts it into actual data based on the geometry defined in emcstructure.cxx. So if i just create these classes and also the ones they depend on, then i should be good to do some digi right? is this what Florian ment by doing simple smearing, or is it something else.

I also want to know what the nonuniformityfile is. In the beginning of emchitproducer, you create a constructor with that particular file in it /input/EmcDigiNoniformitypars.root

also is there a macro available that utilizes these classes to create the digi file?

Thanks a lot
Cheers
Raghav
Re: Digitization and Reconstruction in Fairroot to new detectors [message #13676 is a reply to message #13675] Wed, 27 June 2012 20:42 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 8.28.178*
You don't need structure or other params, you can delete all those things and edit the process hits as you wish. Maybe the mdt hit producer is easier to understand.
Re: Digitization and Reconstruction in Fairroot to new detectors [message #13681 is a reply to message #13676] Thu, 28 June 2012 16:57 Go to previous message
Raghav Kunnawalkam is currently offline  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
Previous Topic: switching off the information output
Next Topic: Cherenkov photons. Accuracy problem with Geant3
Goto Forum:
  


Current Time: Mon Oct 14 19:49:17 CEST 2024

Total time taken to generate the page: 0.00935 seconds