GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » General » Fieldmaps for the half current solenoid
Re: Fieldmaps for the half current solenoid [message #15973 is a reply to message #14938] Sat, 15 March 2014 00:19 Go to previous messageGo to previous message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *dip0.t-ipconnect.de
Hi Prome and Stefano,

The final issue for 1/2 half field map still remain due to combining two stuff between field map and accessor in my point of view.


I assume that total size of cubic(or vector) and structure must be same for half field map.
Then additional field map will be placed into the pandaroot/input/ accodring to transition region with 2 different beam momentum and one common solenoid
TransMap_Low.0150.root
TransMap_Low.0300.root
SolenoidMap_Low1.root
SolenoidMap_Low2.root
SolenoidMap_Low3.root
SolenoidMap_Low4.root


One can access reduced map via /field/PndMultiField by "Half" in user side.
I expect simply additional option "Half" with above naming scheme
Quote:

else if (Map=="Half") {
PndTransMap *map_t= new PndTransMap("TransMap_Low", "R", fBeamMom);
PndDipoleMap *map_d1= new PndDipoleMap("DipoleMap1", "R", fBeamMom);
PndDipoleMap *map_d2= new PndDipoleMap("DipoleMap2", "R", fBeamMom);
PndSolenoidMap *map_s1= new PndSolenoidMap("SolenoidMap_Low1", "R");
PndSolenoidMap *map_s2= new PndSolenoidMap("SolenoidMap_Low2", "R");
PndSolenoidMap *map_s3= new PndSolenoidMap("SolenoidMap_Low3", "R");
PndSolenoidMap *map_s4= new PndSolenoidMap("SolenoidMap_Low4", "R");
AddField(map_t);
AddField(map_d1);
AddField(map_d2);
AddField(map_s1);
AddField(map_s2);
AddField(map_s3);
AddField(map_s4);
}
....



Basically between 1.5 GeV and 3.0 GeV only one transition map should be enough.
However we have already two field map for transition region by producing from Prome.
Then field/PndTransMap.cxx have to be replaced by a piece of hanling for this two transition map.

Quote:

fType = 4;
TString Suffix="";
FairRunSim *fRun= FairRunSim::Instance();
if(fRun) fBeamMom= fRun->GetBeamMom();

//if(fBeamMom< 3.0)Suffix=".0150" ;//inactivate
if(fBeamMom>= 1.5 && fBeamMom <= 2.25)Suffix=".0150" ; //baem momentum close to the 1.5 GeV
if(fBeamMom>= 2.25 && fBeamMom <= 3.0)Suffix=".0300" ; //beam momentum clsoe to the 3.0 GeV
//actually between 1.5 and 3.0 GeV only one transition map.

else if (fBeamMom< 6.0 && fBeamMom >= 3.0)Suffix=".0406";
else if (fBeamMom< 10.0 && fBeamMom >= 6.0 )Suffix=".0890" ;
else if (fBeamMom< 13.0 && fBeamMom >= 10.0)Suffix=".1191";
else if (fBeamMom> 13.0) Suffix=".1500";


TString NewName=mapName;
NewName=mapName+Suffix;
SetName(NewName.Data());
TString dir = getenv("VMCWORKDIR");
fFileName = dir + "/input/" + NewName;
if ( fileType[0] == 'R' ) fFileName += ".root";
else fFileName += ".dat";




I think then we have no problem at all for this issue.
Do I miss somewhere else?

Best wishes,
Donghee







 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: [NOTSUPPORTED]Problem in Installation of external package dec 2013
Next Topic: Target Specifications
Goto Forum:
  


Current Time: Thu Apr 18 06:21:55 CEST 2024

Total time taken to generate the page: 0.01206 seconds