Skip to content
Snippets Groups Projects
Commit b6ffc01f authored by adrien-matta's avatar adrien-matta
Browse files

* Fixing small issue in Example1

parent 211954bb
No related branches found
No related tags found
No related merge requests found
......@@ -58,7 +58,7 @@ void Analysis::Init(){
Si_X_M2 = 0;
Si_Y_M2 = 0;
ZTarget = 0;
TargetThickness = m_DetectorManager->GetTargetThickness();
TargetThickness = m_DetectorManager->GetTargetThickness()*micrometer;
// Energy loss table: the G4Table are generated by the simulation
He3CD2 = EnergyLoss("He3_CD2.G4table","G4Table",100 );
He3Al = EnergyLoss("He3_Al.G4table","G4Table",10);
......@@ -80,6 +80,7 @@ void Analysis::TreatEvent(){
// Beam energy is measured using F3 and F2 plastic TOF
double BeamEnergy = Rand.Gaus(Initial->GetIncidentInitialKineticEnergy(),4.5);
BeamEnergy = Li11CD2.Slow(BeamEnergy,TargetThickness/2.,0);
He10Reaction->SetBeamEnergy(BeamEnergy);
//////////////////////////// LOOP on MUST2 + SSSD Hit //////////////////
for(unsigned int countSSSD = 0 ; countSSSD < SSSD->Energy.size() ; countSSSD++){
......@@ -158,13 +159,15 @@ void Analysis::TreatEvent(){
/************************************************/
// Part 3 : Excitation Energy Calculation
Ex = He10Reaction -> ReconstructRelativistic( ELab , ThetaLab );
ThetaLab=ThetaLab/deg;
/************************************************/
/************************************************/
// Part 4 : Theta CM Calculation
ThetaCM = He10Reaction -> EnergyLabToThetaCM( ELab , 0)/deg;
ThetaCM = He10Reaction -> EnergyLabToThetaCM( ELab , ThetaLab)/deg;
ThetaLab=ThetaLab/deg;
/************************************************/
}
} //end loop SSSD
}//end loop MUST2
......
......@@ -16,7 +16,7 @@ EDE= (TCutG*) File_EDE->FindObjectAny("EDE");
////////////////////////////////////////////////////////////////////////////////
void LoadChain(){
chain = new TChain("ResultTree");
chain = new TChain("PhysicsTree");
chain->Add("../../Outputs/Analysis/Example1.root");
}
......@@ -41,7 +41,8 @@ ETOF->Draw("same");
// Kinematical Line //
c1->cd(3);
chain->Draw("ELab:ThetaLab>>hKine(500,0,45,400,0,40)","MUST2.CsI_E<0 && MUST2.TelescopeNumber<5 && EDE && ETOF","colz");
//chain->Draw("ELab:ThetaLab>>hKine(500,0,45,400,0,40)","MUST2.CsI_E<0 && MUST2.TelescopeNumber<5 && EDE && ETOF","colz");
chain->Draw("ELab:ThetaLab","MUST2.CsI_E<0 && MUST2.TelescopeNumber<5 && EDE && ETOF","colz");
NPL::Reaction r("11Li(d,3He)10He@553");
r.SetExcitationHeavy(1.4);
......
ConfigMust2
MAX_STRIP_MULTIPLICITY 1
STRIP_ENERGY_MATCHING_NUMBER_OF_SIGMA 5
STRIP_ENERGY_MATCHING_SIGMA 0.02
STRIP_ENERGY_MATCHING_SIGMA 0.02
DISABLE_CHANNEL MM1STRY12
DISABLE_CHANNEL MM2STRY12
DISABLE_CHANNEL MM3STRY12
DISABLE_CHANNEL MM4STRY12
DISABLE_CHANNEL MM1STRX12
DISABLE_CHANNEL MM2STRX12
DISABLE_CHANNEL MM3STRX12
DISABLE_CHANNEL MM4STRX12
DISABLE_CHANNEL MM1STRY124
DISABLE_CHANNEL MM4STRX3
DISABLE_CHANNEL MM4STRX4
DISABLE_CHANNEL MM4STRX5
DISABLE_CHANNEL MM4STRX6
DISABLE_CHANNEL MM4STRX8
DISABLE_CHANNEL MM4STRX31
DISABLE_CHANNEL MM2STRX12
DISABLE_CHANNEL MM3STRX12
DISABLE_CHANNEL MM4STRX12
DISABLE_ALL MM5
DISABLE_ALL MM6
DISABLE_ALL MM7
DISABLE_ALL MM8
SI_X_E_RAW_THRESHOLD 0
CSI_E_RAW_THRESHOLD 0
CSI_SIZE 256
CSI_E_RAW_THRESHOLD 0
CSI_SIZE 256
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment