GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » Reproduce the result of Pos (Bormio2012) 018 (Analysis and simulation )
Reproduce the result of Pos (Bormio2012) 018 [message #15937] Mon, 10 March 2014 18:52 Go to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 117.198.193*
Hello all

I am new to PANDAroot. I have installed the external packages of apr13 and the pandaroot from the website https://subversion.gsi.de/fairroot/pandaroot/trunk pandaroot. Everything is working very well. I have also tried to understand tutorial of apr13 given in pandaroot/tutorial/apr13. which is also working well.

Now I want to do the simulation for x(3872) and using the master thesis Simulation of X(3872) Decays Using the PandaRoot Framework by Martin Johannes Galuska and I simply want to reproduce the data given in thesis and in paper Pos (Bormio2012) 018. Decay file is also ready. simulation file is also ready and to prepare the this file I have taken the help of tut_sim.C file in apr13 tutorial. Now what are the modification needed for the invariant mass spectroscopy and missing mass and lastly the resonance scan. I know some the Physics of it but little bit weak in simulation. Please help me out.

Thank You


Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #15938 is a reply to message #15937] Mon, 10 March 2014 19:09 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
From: *pools.arcor-ip.net
Hello,

I can send you all the files, macros and bash scripts which I used to run the simulation and to analyse the simulated data. However, the simulation was performed with an older version of PandaRoot and since that time a lot has changed in the software framework. As I currently have a lot of urgent deadlines to meet, I cannot make the changes necessary to the aforementioned files for you, so you might need the help of someone who is already experienced in PandaRoot simulations to help you.

For the resonance scan I ran a simulation for the background from the direct process and for the signal for each of the 20 points using a bash script, but still some manual work is needed to get to the final plot.

Kind regards,
Martin
Re: Reproduce the result of Pos (Bormio2012) 018 [message #15939 is a reply to message #15938] Tue, 11 March 2014 10:14 Go to previous messageGo to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
Hello

Thank you very much for your response.

I request you to send me the files I will try to modified it.


Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #15941 is a reply to message #15939] Wed, 12 March 2014 01:23 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
From: *pools.arcor-ip.net
Hello,

please send me your email address either as a reply here or as a private message. I will email you the files and my Master thesis which contains some documentation.

Kind regards,
Martin
Re: Reproduce the result of Pos (Bormio2012) 018 [message #15948 is a reply to message #15941] Thu, 13 March 2014 03:31 Go to previous messageGo to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
Hello

My email address is nakulphy@gmail.com

Please send the required files to this email id.

Thank you


Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #15957 is a reply to message #15948] Thu, 13 March 2014 12:10 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
From: 140.181.11*
Thank you. I will send you the files as soon as I have access to my backup harddisk which is in my office. As you know this week there is the PANDA CM taking place. Next week I will be at another conference, so it might take a couple of days until I can send you the files. Thank you for your patience.
Re: Reproduce the result of Pos (Bormio2012) 018 [message #15961 is a reply to message #15957] Fri, 14 March 2014 10:27 Go to previous messageGo to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
Thank you very much.

Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16007 is a reply to message #15961] Tue, 18 March 2014 20:30 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
From: *physik.uni-giessen.de
Hello,

I have just sent you the files you asked for.

Kind regards,
Martin
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16038 is a reply to message #16007] Tue, 25 March 2014 11:37 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
From: *pool.mediaWays.net
Quote:

Hello Dr. Martin

I was trying to run the run_sim_tpcorsttcombi.C but I am facing some problem it it.

in detector setup the file tpc.geo is not available in pandaroot/geometry in jan14 pandaroot with revision 24201.

therefore I tried to modified it but I may be wrong in it. The simulation is running but gives some error and that I am sending it to you with the file name "tpc". I am also sending you the another error in the file "root tpc".

I am also sending your simulation file where I have changed the detector and also the root file that is generated in the simulation.

Please help me out..

Thank you


Hello Nakul,

please use the forum for support requests as I am very busy currently and other people will be able to help you as well.

As I have written earlier, there have been many changes in the software since I did the X(3872) analysis. One of these changes is that the TPC option was decided against and therefore, this detector geometry is not included in PandaRoot anymore.

I suggest to use the macros from $VMCWORKDIR/macro/run as a basis for your files. You could do a diff of the macros I sent you with the macros in macro/run from the PandaRoot svn revision I used for my analysis (I quoted that number at some point in my Master thesis) to get an idea of the code I added to the standard macros and then you can apply similar changes to the current version of the macros in the PandaRoot revision you are using right now.

Kind regards,
Martin

PS: I do not have a Dr. title.
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16084 is a reply to message #15937] Fri, 28 March 2014 09:28 Go to previous messageGo to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
hello

I have tried to do as you suggested and I am sending you things which I got. the simulation and digi works but there is some problem with the reco. I am sending you the reco error.

nakul@linux-kqx4:~/pandaroot/x3872/x3872> root run_reco_tpcorsttcombi.C
*******************************************
* *
* W E L C O M E to R O O T *
* *
* Version 5.34/05 14 February 2013 *
* *
* You are welcome to visit our Web site *
* http://root.cern.ch *
* *
*******************************************

ROOT 5.34/05 (tags/v5-34-05@48582, Oct 03 2013, 16:35:25 on linuxx8664gcc)

CINT/ROOT C/C++ Interpreter version 5.18.00, July 2, 2010
Type ? for help. Commands must be C++ statements.
Enclose multiple statements between { }.
root [0]
Processing run_reco_tpcorsttcombi.C...
Error in <TClonesArray::SetClass>: called with a null pointer
FairRootManager::OpenOutFile("reco_sttcombi.root")
Use detector with STT setup.???
Info in (PndGeoHandling::Instance): Making a new instance using the framework.
[INFO ] The input consists out of the following trees and files:
[INFO ] - cbmsim
[INFO ] - digi_sttcombi.root
[INFO ] - FriendTree_1
[INFO ] - points_sttcombi.root
[INFO ] Parameter and input file are available, Assure that basic info is there for the run!
[INFO ] The number of entries in chain is 10

*************************************************************
initialisation for run id 1396013441
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet
Info in <TGeoManager::CloseGeometry>: Geometry loaded from file...
Info in <TGeoManager::SetTopVolume>: Top volume is cave. Master volume is cave
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
Info in <TGeoManager::Voxelize>: Voxelizing...
Info in <TGeoManager::CountLevels>: max level = 12, max placements = 4550
Info in <TGeoManager::CloseGeometry>: 2003695 nodes/ 5810 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
Container FairBaseParSet initialized from ROOT file.
-I container name PndGeoSttPar
-I container name PndGeoFtsPar

*************************************************************
initialisation for run id 1396013441
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet
-I- FairRunTimeDB::InitContainer() PndSensorNamePar
Container PndSensorNamePar initialized from ROOT file.
-I- FairRunTimeDB::InitContainer() PndGeoSttPar
Container PndGeoSttPar initialized from ROOT file.
-I- FairRunTimeDB::InitContainer() PndGemDetectors
Gem_Disk1_Gem1_Sensor_GEMmixture type 0 has 786 front and 4250 back channels
Gem_Disk1_Gem6_Sensor_GEMmixture type 2 has 4750 front and 9000 back channels
Gem_Disk2_Gem1_Sensor_GEMmixture type 0 has 786 front and 5350 back channels
Gem_Disk2_Gem6_Sensor_GEMmixture type 2 has 5850 front and 11200 back channels
Gem_Disk3_Gem1_Sensor_GEMmixture type 0 has 786 front and 7150 back channels
Gem_Disk3_Gem6_Sensor_GEMmixture type 2 has 7650 front and 14800 back channels
PndGemDetectors initialized from Ascii file
-I- FairRunTimeDB::InitContainer() PndGeoFtsPar
Container PndGeoFtsPar initialized from ROOT file.
PndFieldCreator::SetParm() 0x880db90

*************************************************************
initialisation for run id 1396013441
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet
-I- FairRunTimeDB::InitContainer() PndSensorNamePar
-I- FairRunTimeDB::InitContainer() PndGeoSttPar
-I- FairRunTimeDB::InitContainer() PndGemDetectors
PndGemDetectors initialized from Ascii file
-I- FairRunTimeDB::InitContainer() PndGeoFtsPar
-I- FairRunTimeDB::InitContainer() PndMultiFieldPar
Container PndMultiFieldPar initialized from ROOT file.
OBJ: PndTransPar PndTransPar Trans. Field parameter container
OBJ: PndDipole1Par PndDipole1Par Dipole Field parameter container
OBJ: PndDipole2Par PndDipole2Par Dipole Field parameter container
OBJ: PndSolenoid1Par PndSolenoid1Par Solenoid 1st region parameter container
OBJ: PndSolenoid2Par PndSolenoid2Par Solenoid 2nd region parameter container
OBJ: PndSolenoid3Par PndSolenoid3Par Solenoid 3rd region parameter container
OBJ: PndSolenoid4Par PndSolenoid4Par Solenoid 4th region parameter container
[INFO ] PndFieldMap: Reading field map from ROOT file /home/nakul/pandaroot/jan14/input/TransMap.0890.root
[INFO ] PndFieldMap: Reading field map from ROOT file /home/nakul/pandaroot/jan14/input/DipoleMap1.0890.root
[INFO ] PndFieldMap: Reading field map from ROOT file /home/nakul/pandaroot/jan14/input/DipoleMap2.0890.root
[INFO ] PndFieldMap: Reading field map from ROOT file /home/nakul/pandaroot/jan14/input/SolenoidMap1.root
[INFO ] PndFieldMap: Reading field map from ROOT file /home/nakul/pandaroot/jan14/input/SolenoidMap2.root
[INFO ] PndFieldMap: Reading field map from ROOT file /home/nakul/pandaroot/jan14/input/SolenoidMap3.root
[INFO ] PndFieldMap: Reading field map from ROOT file /home/nakul/pandaroot/jan14/input/SolenoidMap4.root
------------------------FairGeane::Init ()------------------------------
Loading Geant3 libraries ...
Loading Geant3 libraries ... finished

MZSTOR. ZEBRA table base TAB(0) in /MZCC/ at adr 255514687 F3AD83F HEX

MZSTOR. Initialize Store 0 in /GCBANK/
with Store/Table at absolute adrs 255623941 255514687
HEX F3C8305 F3AD83F
HEX 1AB2A 0
relative adrs 109354 0
with 1 Str. in 2 Links in 5300 Low words in 4999970 words.
This store has a fence of 16 words.

MZLOGL. Set Log Level 0 for store 0
1***** GEANT Version 3.21/11 Released on 100298
0***** Correction Cradle Version 0.1100

MZDIV. Initialize Division Constant in Store 0
NW/NWMAX= 20004000000, MODE/KIND= 1 2
Division 20 initialized.

MZLINK. Initialize Link Area /GCLINK/ for Store 0 NL/NS= 20 20

MZLINK. Initialize Link Area /GCSLNK/ for Store 0 NL/NS= 100 100
-I- G3Config: Geant3 with TGeo has been created for Geane.
-I- Geane.C: NOPRNT flag set to 1
-I- Geane.C: IERR flags are not printed. If you want to switch them on, please set fErtrio1->noprnt = 0 in Geane.C
Energy straggling area parameter from user set to: 0.999

Calculating cross section tables, see gphysi.dat for more information

Cross section calculation concluded successfully
I- FairGeane::FairGeane: Geane is Initialized
-W- PndMvdRiemannTrackFinderTask::Init: No Branch Names given with AddHitBranch(TString branchName)! Standard BranchNames taken!
-I- PndMvdRiemannTrackFinderTask: Initialisation successfull
-I- -------------------
-I- PndSttMvdGemTracking: using branches MVDHitsPixel MVDHitsStrip STTHit GEMHit
-I- to change one or more of these use PndSttMvdGemTracking:SetBranchName( TStrings ); the order of TStrings is mvd pixel name, mvd strip name, stt name, gem name
starting track for extrapolation SttMvdTrack SttMvdTrackCand
-I- -------------------
[INFO ] Branch: SttMvdTrackCand not found in Tree
[INFO ] Branch: SttMvdTrackCand not found in Tree
fStartTrackCandBranchName SttMvdTrackCand not found
Error in <PndSttMvdGemTracking::PndSttMvdGemTracking:Init>: stt + mvd trackcand - array not found!
-I- PndMCTrackAssociator::Init
[INFO ] Branch: SttMvdGemTrack not found in Tree
[INFO ] Branch: SttMvdGemTrack not found in Tree
-I- PndMCTrackAssociator::Init: No PndTrack array!
-I- PndRecoKalmanTask:Init :: Using GeaneTrackRep
#####Mapper:filltubearray######
fGeoType=1
*** PndRecoKalmanFit::Init MVDHitsStrip array found
*** PndRecoKalmanFit::Init MVDHitsPixel array found
[INFO ] Branch: SttHelixHit not found in Tree
[INFO ] Branch: SttHelixHit not found in Tree
*** PndRecoKalmanFit::Init SttHit array found
*** PndRecoKalmanFit::Init GEMHit array found
*** PndRecoKalmanFit::Init MdtHit array found
*** PndRecoKalmanFit::Init FtsHit array found
===PndRecoKalmanFit::Init() finished ===================================================
[INFO ] Branch: SttMvdGemTrack not found in Tree
[INFO ] Branch: SttMvdGemTrack not found in Tree
Error in <PndRecoKalmanTask::PndRecoKalmanTask::Init>: track-array not found!
-I- PndMCTrackAssociator::Init
[INFO ] Branch: SttMvdGemGenTrack not found in Tree
[INFO ] Branch: SttMvdGemGenTrack not found in Tree
-I- PndMCTrackAssociator::Init: No PndTrack array!
[INFO ] The number of entries in chain is 10
terminate called after throwing an instance of 'std::out_of_range'
what(): vector::_M_range_check
nakul@linux-kqx4:~/pandaroot/x3872/x3872>


thank you


Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16086 is a reply to message #16084] Fri, 28 March 2014 09:36 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Again,
please use the macros in macro/run, or the ones in the tutorial.
You need only to change the simulation macro sim_complete, setting your .dec file and the beam momentun. Apart from this, the rest should stay untouched, digi_complete reco_complete and pid_complete.
After pid, you can run your analysis macro as Martin did, even if I suppose you need to modify something. ana_complete is a good example, but you need to modify it maybe for your purposes.

The macros you are using are somehow obsolete, you use tasks which are not used anymore and they could have bugs making everything crash,
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16088 is a reply to message #16086] Fri, 28 March 2014 09:55 Go to previous messageGo to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
Thank you very much.

yes its work. actually that old files were also working with some small modifications. that is why I was using that files.

macro/run is now working properly
In sim, didi, rec, and pid_complete.C is now working properly and now I am focusing on the analysis file. but how to check that our simulations going on way that we are expecting ? in all the simulations at the end it gives the message "Test passed All ok". but how to verify this ?

during simulation the message comes
**** GTRIGI: IEVENT= 4 IDEVT= 4 Random Seeds = 123456 0
[INFO ] FairPrimaryGenerator: (Event 4) 6 primary tracks from vertex (-0.071961, 0.202136, 0.100250 ) with beam gradiant (0.000000, 0.000000) Event Time = 0.000000 (ns)
POINT EXECUTION *********************
in which what is the meaning of beam gradiant (,) ?

Thank you



Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16089 is a reply to message #16086] Fri, 28 March 2014 09:56 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
From: *wlan.uni-giessen.de
Yes, Stefano is right. I have not changed the digitization, reconstruction or pid macros. I only added some custom code in the simulation macro which you need to more or less copy over to the current version of the simulation macro from macro/run. The analysis macro will probably have to be rewritten due to many changes. The shell scripts should work fine with no or only minor adjustments.
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16090 is a reply to message #16089] Fri, 28 March 2014 10:12 Go to previous messageGo to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
yes sir I was trying to do the both older and the macro/run.
and there is no any problem in macro/run it is working properly. But how to check that our simulation is going as our expectations. and how to check the status of the ongoing simulations ???


Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16091 is a reply to message #16090] Fri, 28 March 2014 10:40 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
From: *wlan.uni-giessen.de
Sorry, I don't understand your question.
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16093 is a reply to message #16091] Fri, 28 March 2014 10:49 Go to previous messageGo to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
the root files that are created during the simulations such as sim_complete.root, simparams.root, etc...
how to study these files ? and what is the meaning of beam gradiant which is created during the sim_complete.C ? which is given below.

**** GTRIGI: IEVENT= 9 IDEVT= 9 Random Seeds = 123456 0
[INFO ] FairPrimaryGenerator: (Event 9) 8 primary tracks from vertex (-0.189898, 0.091527, 0.224391 ) with beam gradiant (0.000000, 0.000000) Event Time = 0.000000 (ns)
POINT EXECUTION *********************


Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16094 is a reply to message #16093] Fri, 28 March 2014 10:59 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *gsi.de
Hello nakulphy,

nakulphy wrote on Fri, 28 March 2014 10:49
the root files that are created during the simulations such as sim_complete.root, simparams.root, etc...
how to study these files ?

What do you want to study? Is is detector performance of the individual subsystems? Then you will have to dig into the sim, reco and pid files.
Is it particle physics, i.e. studying the detected results of a specific event generator input (e.g. "your" channel), then you don't need to worry about those files and go on with analysis.

nakulphy wrote on Fri, 28 March 2014 10:49
and what is the meaning of beam gradiant which is created during the sim_complete.C ? which is given below.

**** GTRIGI: IEVENT= 9 IDEVT= 9 Random Seeds = 123456 0
[INFO ] FairPrimaryGenerator: (Event 9) 8 primary tracks from vertex (-0.189898, 0.091527, 0.224391 ) with beam gradiant (0.000000, 0.000000) Event Time = 0.000000 (ns)
POINT EXECUTION *********************

This is output of the framework (event generator interface) during simulation. We don't use sophisticated beam properties, yet, so it's nothing to worry about.

Cheers
Ralf Kliemt
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16097 is a reply to message #16094] Fri, 28 March 2014 11:20 Go to previous messageGo to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
Thank you very much....

Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16381 is a reply to message #16097] Mon, 21 April 2014 11:13 Go to previous messageGo to next message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
Hello

simulation, digitization, reconstruction, particle identification (from macro/run) is completed and the calculation of no of event is also completed and for this calculations we have used Mathematica 9.
now the analysis part is left. what is the method that we have to follow ?

Thank you


Nakul Soni
Re: Reproduce the result of Pos (Bormio2012) 018 [message #16385 is a reply to message #16381] Mon, 21 April 2014 22:48 Go to previous messageGo to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
From: *pool.mediaWays.net
Hello,

I have not written many analysis macros recently, but I assume that the hint Stefano gave you in his post from March, 28th is still valid.

Quote:

ana_complete is a good example, but you need to modify it maybe for your purposes.


You can find the file in macro/run of your PandaRoot source folder or here. You can use it as an example of how to write an analysis macro. For the analysis strategy you can follow the analysis macro which I sent you in the email some time ago. The file is called run_ana_tpcorsttcombi.C, but don't expect the code to work without (numerous) modifications. There were lots of changes to the analysis framework after I wrote the code. I suggest to take advantage of the new PID functionality just to mention one thing that comes to mind.

Kind regards,
Martin

[Updated on: Mon, 21 April 2014 22:49]

Report message to a moderator

Re: Reproduce the result of Pos (Bormio2012) 018 [message #16389 is a reply to message #16385] Tue, 22 April 2014 04:46 Go to previous message
nakulphy is currently offline  nakulphy
Messages: 13
Registered: March 2014
Location: India
occasional visitor
From: 14.139.121*
Hello

Thank you. I am using the same macro.


Nakul Soni
Previous Topic: pi0 in fast simulation
Next Topic: [FIXED] Error in <PndAnaPidCombiner::Apply>
Goto Forum:
  


Current Time: Thu Dec 05 23:57:11 CET 2024

Total time taken to generate the page: 0.00742 seconds