From 44f7ab75c822a75e1620f4aebfa9a293f8dcf462 Mon Sep 17 00:00:00 2001
From: adrien matta <matta@lpccaen.in2p3.fr>
Date: Mon, 20 Jul 2020 17:23:03 +0200
Subject: [PATCH] * Strasse outer barrel analysis
---
NPLib/Detectors/Strasse/TStrassePhysics.cxx | 188 ++++++++++++++++++--
NPLib/Detectors/Strasse/TStrassePhysics.h | 3 +-
Projects/Strasse/Control.cxx | 52 +++---
3 files changed, 203 insertions(+), 40 deletions(-)
diff --git a/NPLib/Detectors/Strasse/TStrassePhysics.cxx b/NPLib/Detectors/Strasse/TStrassePhysics.cxx
index 54fb6719c..a0274a76c 100644
--- a/NPLib/Detectors/Strasse/TStrassePhysics.cxx
+++ b/NPLib/Detectors/Strasse/TStrassePhysics.cxx
@@ -204,7 +204,100 @@ void TStrassePhysics::AddInnerDetector(double R, double Z, double Phi, double Sh
///////////////////////////////////////////////////////////////////////////
void TStrassePhysics::AddOuterDetector(double R, double Z, double Phi, double Shift){
-}
+ m_NumberOfOuterDetectors++;
+ double ActiveWidth = Outer_Wafer_Width-2.*Outer_Wafer_GuardRing;
+ double ActiveLength = Outer_Wafer_Length-Outer_Wafer_PADExternal-Outer_Wafer_PADInternal-2*Outer_Wafer_GuardRing;
+ double LongitudinalPitch = ActiveWidth/Outer_Wafer_LongitudinalStrips;
+ double TransversePitch = ActiveLength/Outer_Wafer_TransverseStrips;
+//cout << ActiveWidth << " " << ActiveLength << " " << LongitudinalPitch << " " << TransversePitch << endl;
+
+ // Vector C position of detector face center
+ TVector3 C(Shift,R,Z);// center of the whole detector, including PCB
+ C.RotateZ(-Phi);
+
+ // Vector W normal to detector face (pointing to the back)
+ TVector3 W(0,1,0);
+ W.RotateZ(-Phi);
+
+ // Vector U on detector face (parallel to Z axis/longitudinal strips)
+ TVector3 U = TVector3(0,0,1);
+ // Vector V on detector face (parallel to transverse strips)
+ TVector3 V = W.Cross(U);
+
+ // Adding position for downstream silicon:
+ // Moving to corner of the silicon
+ TVector3 P_1_1 = C
+ +U*0.5*(Outer_PCB_UpstreamWidth-Outer_PCB_DownstreamWidth) // In between wafer
+ -U*0.5*Outer_PCB_MidWidth // Internal wafer edge
+ -U*Outer_Wafer_Length // External wafer edge
+ +U*(Outer_Wafer_GuardRing+Outer_Wafer_PADExternal) // External active wafer edge
+ +U*0.5*TransversePitch // middle of strip
+ -V*0.5*(Outer_PCB_StarboardWidth-Outer_PCB_PortWidth)
+ -V*0.5*Outer_Wafer_Width
+ +V*Outer_Wafer_GuardRing
+ +V*0.5*LongitudinalPitch; // middle of strip
+
+ vector<double> lineX;
+ vector<double> lineY;
+ vector<double> lineZ;
+
+ vector<vector<double>> OneDetectorStripPositionX;
+ vector<vector<double>> OneDetectorStripPositionY;
+ vector<vector<double>> OneDetectorStripPositionZ;
+
+ TVector3 P;
+ for(int i=0; i<Outer_Wafer_TransverseStrips; i++){
+ lineX.clear();
+ lineY.clear();
+ lineZ.clear();
+ for(int j=0; j<Outer_Wafer_LongitudinalStrips; j++){
+ P = P_1_1 + i*U*TransversePitch + j*V*LongitudinalPitch;
+ lineX.push_back(P.X());
+ lineY.push_back(P.Y());
+ lineZ.push_back(P.Z());
+ }
+
+ OneDetectorStripPositionX.push_back(lineX);
+ OneDetectorStripPositionY.push_back(lineY);
+ OneDetectorStripPositionZ.push_back(lineZ);
+ }
+
+ m_OuterStripPositionX.push_back(OneDetectorStripPositionX);
+ m_OuterStripPositionY.push_back(OneDetectorStripPositionY);
+ m_OuterStripPositionZ.push_back(OneDetectorStripPositionZ);
+
+ // Adding position for upstream silicon:
+ // Moving to corner of the silicon
+ P_1_1 = C
+ +U*0.5*(Outer_PCB_UpstreamWidth-Outer_PCB_DownstreamWidth) // In between wafer
+ +U*0.5*Outer_PCB_MidWidth // Internal wafer edge
+ +U*(Outer_Wafer_GuardRing+Outer_Wafer_PADInternal) // Internal active wafer edge
+ +U*0.5*TransversePitch// middle of strip
+ -V*0.5*(Outer_PCB_StarboardWidth-Outer_PCB_PortWidth)
+ -V*0.5*Outer_Wafer_Width
+ +V*Outer_Wafer_GuardRing
+ +V*0.5*LongitudinalPitch; // middle of strip
+
+ for(int i=0; i<Outer_Wafer_TransverseStrips; i++){
+ lineX.clear();
+ lineY.clear();
+ lineZ.clear();
+
+ for(int j=0; j<Outer_Wafer_LongitudinalStrips; j++){
+ P = P_1_1 + i*U*TransversePitch + j*V*LongitudinalPitch;
+ lineX.push_back(P.X());
+ lineY.push_back(P.Y());
+ lineZ.push_back(P.Z());
+
+ }
+
+ m_OuterStripPositionX[m_NumberOfOuterDetectors-1].push_back(lineX);
+ m_OuterStripPositionY[m_NumberOfOuterDetectors-1].push_back(lineY);
+ m_OuterStripPositionZ[m_NumberOfOuterDetectors-1].push_back(lineZ);
+
+ }
+}
+
///////////////////////////////////////////////////////////////////////////
TVector3 TStrassePhysics::GetInnerPositionOfInteraction(const int i){
@@ -218,9 +311,9 @@ TVector3 TStrassePhysics::GetInnerPositionOfInteraction(const int i){
///////////////////////////////////////////////////////////////////////////
TVector3 TStrassePhysics::GetOuterPositionOfInteraction(const int i){
- TVector3 Position = TVector3(GetOuterStripPositionX(DetectorNumber[i], OuterStripL[i], OuterStripT[i]),
- GetOuterStripPositionY(DetectorNumber[i], OuterStripL[i], OuterStripT[i]),
- GetOuterStripPositionZ(DetectorNumber[i], OuterStripL[i], OuterStripT[i]));
+ TVector3 Position = TVector3(GetOuterStripPositionX(DetectorNumber[i], OuterStripT[i], OuterStripL[i]),
+ GetOuterStripPositionY(DetectorNumber[i], OuterStripT[i], OuterStripL[i]),
+ GetOuterStripPositionZ(DetectorNumber[i], OuterStripT[i], OuterStripL[i]));
return Position;
}
@@ -262,16 +355,17 @@ void TStrassePhysics::BuildPhysicalEvent() {
PreTreat();
if(1 /*CheckEvent() == 1*/){
- vector<TVector2> couple = Match_X_Y();
+ vector<TVector2> inner = MatchInner();
+ vector<TVector2> outer = MatchOuter();
- EventMultiplicity = couple.size();
- for(unsigned int i=0; i<couple.size(); i++){
- int N = m_PreTreatedData->GetInner_TE_DetectorNbr(couple[i].X());
- int T = m_PreTreatedData->GetInner_TE_StripNbr(couple[i].X());
- int L = m_PreTreatedData->GetInner_LE_StripNbr(couple[i].Y());
+ EventMultiplicity = inner.size();
+
+ for(unsigned int i=0; i<inner.size(); i++){
+ int N = m_PreTreatedData->GetInner_TE_DetectorNbr(inner[i].X());
+ int T = m_PreTreatedData->GetInner_TE_StripNbr(inner[i].X());
+ int L = m_PreTreatedData->GetInner_LE_StripNbr(inner[i].Y());
- double TE = m_PreTreatedData->GetInner_TE_Energy(couple[i].X());
- double LE = m_PreTreatedData->GetInner_LE_Energy(couple[i].Y());
+ double TE = m_PreTreatedData->GetInner_TE_Energy(inner[i].X());
DetectorNumber.push_back(N);
InnerStripT.push_back(T);
InnerStripL.push_back(L);
@@ -279,14 +373,41 @@ void TStrassePhysics::BuildPhysicalEvent() {
InnerPosX.push_back(GetInnerPositionOfInteraction(i).x());
InnerPosY.push_back(GetInnerPositionOfInteraction(i).y());
InnerPosZ.push_back(GetInnerPositionOfInteraction(i).z());
-
+ // look for outer
+ double outerE = 0;
+ int outerT=0;
+ int outerL=0;
+ for(unsigned int j=0; j<outer.size(); j++){
+ if(m_PreTreatedData->GetInner_TE_DetectorNbr(outer[j].X())==N){
+ outerE = m_PreTreatedData->GetOuter_TE_Energy(outer[j].X());
+ outerT = m_PreTreatedData->GetOuter_TE_StripNbr(outer[j].X());
+ outerL = m_PreTreatedData->GetOuter_LE_StripNbr(outer[j].Y());
+ }
+ }
+ if(outerE){
+ OuterStripT.push_back(outerT);
+ OuterStripL.push_back(outerL);
+ E.push_back(outerE);
+ OuterPosX.push_back(GetOuterPositionOfInteraction(i).x());
+ OuterPosY.push_back(GetOuterPositionOfInteraction(i).y());
+ OuterPosZ.push_back(GetOuterPositionOfInteraction(i).z());
+ }
+ else{
+ OuterStripT.push_back(-1000);
+ OuterStripL.push_back(-1000);
+ E.push_back(-1000);
+ OuterPosX.push_back(-1000);
+ OuterPosY.push_back(-1000);
+ OuterPosZ.push_back(-1000);
+ }
}
+
}
}
///////////////////////////////////////////////////////////////////////////
-vector<TVector2> TStrassePhysics::Match_X_Y(){
+vector<TVector2> TStrassePhysics::MatchInner(){
vector<TVector2> ArrayOfGoodCouple;
static unsigned int m_TEMult, m_LEMult;
@@ -308,7 +429,7 @@ vector<TVector2> TStrassePhysics::Match_X_Y(){
if(XDetNbr == YDetNbr){
// Declaration of variable for clarity
double TE = m_PreTreatedData->GetInner_TE_Energy(i);
- double LE = m_PreTreatedData->GetInner_LE_Energy(i);
+ double LE = m_PreTreatedData->GetInner_LE_Energy(j);
// look if energy matches
if(abs(TE-LE)/2.<m_StripEnergyMatching){
@@ -320,6 +441,42 @@ vector<TVector2> TStrassePhysics::Match_X_Y(){
return ArrayOfGoodCouple;
}
+///////////////////////////////////////////////////////////////////////////
+vector<TVector2> TStrassePhysics::MatchOuter(){
+ vector<TVector2> ArrayOfGoodCouple;
+
+ static unsigned int m_TEMult, m_LEMult;
+ m_TEMult = m_PreTreatedData->GetOuterMultTEnergy();
+ m_LEMult = m_PreTreatedData->GetOuterMultLEnergy();
+
+ if(m_TEMult>m_MaximumStripMultiplicityAllowed || m_LEMult>m_MaximumStripMultiplicityAllowed){
+ return ArrayOfGoodCouple;
+ }
+
+ for(unsigned int i=0; i<m_TEMult; i++){
+ for(unsigned int j=0; j<m_LEMult; j++){
+
+ // Declaration of variable for clarity
+ int XDetNbr = m_PreTreatedData->GetOuter_TE_DetectorNbr(i);
+ int YDetNbr = m_PreTreatedData->GetOuter_LE_DetectorNbr(j);
+
+ // if same detector check energy
+ if(XDetNbr == YDetNbr){
+ // Declaration of variable for clarity
+ double TE = m_PreTreatedData->GetOuter_TE_Energy(i);
+ double LE = m_PreTreatedData->GetOuter_LE_Energy(j);
+
+ // look if energy matches
+ if(abs(TE-LE)/2.<m_StripEnergyMatching){
+ ArrayOfGoodCouple.push_back(TVector2(i,j));
+ }
+ }
+ }
+ }
+
+ return ArrayOfGoodCouple;
+}
+
///////////////////////////////////////////////////////////////////////////
int TStrassePhysics::CheckEvent(){
@@ -460,7 +617,6 @@ void TStrassePhysics::ReadAnalysisConfig() {
///////////////////////////////////////////////////////////////////////////
void TStrassePhysics::Clear() {
EventMultiplicity = 0;
-
// DSSD
DetectorNumber.clear();
E.clear();
diff --git a/NPLib/Detectors/Strasse/TStrassePhysics.h b/NPLib/Detectors/Strasse/TStrassePhysics.h
index 901a92117..fb0f34902 100644
--- a/NPLib/Detectors/Strasse/TStrassePhysics.h
+++ b/NPLib/Detectors/Strasse/TStrassePhysics.h
@@ -58,7 +58,8 @@ class TStrassePhysics : public TObject, public NPL::VDetector {
void Clear(const Option_t*) {};
public:
- vector<TVector2> Match_X_Y();
+ vector<TVector2> MatchInner();
+ vector<TVector2> MatchOuter();
int CheckEvent();
//////////////////////////////////////////////////////////////
diff --git a/Projects/Strasse/Control.cxx b/Projects/Strasse/Control.cxx
index 2a5a4e1eb..0d0d9840f 100644
--- a/Projects/Strasse/Control.cxx
+++ b/Projects/Strasse/Control.cxx
@@ -1,33 +1,39 @@
void Control(){
TFile* file = TFile::Open("../../Outputs/Analysis/PhysicsTree.root");
TTree* PhysicsTree= (TTree*) file->FindObjectAny("PhysicsTree");
- TCanvas* c = new TCanvas("Control","Control",1000,1000);
- c->Divide(2,2);
- //string cond = "InnerPosX!=-60 && InnerStripL[0]==64 && InnerStripT[0]==152";
-
- string cond = "InnerPosX!=-60 && (DetectorNumber==3 || DetectorNumber==1)";
- c->cd(1);
-// PhysicsTree->Draw("InnerPosX:fDetected_Position_X",cond.c_str(),"col") ;
-//
- PhysicsTree->Draw("InnerPosX-fDetected_Position_X>>h(1000,-2,2)",cond.c_str(),"col") ;
- c->cd(2);
- //PhysicsTree->Draw("InnerPosY:fDetected_Position_Y",cond.c_str(),"col") ;
-
- PhysicsTree->Draw("InnerPosY-fDetected_Position_Y>>h2(1000,-2,2)",cond.c_str(),"col") ;
- c->cd(3);
- //PhysicsTree->Draw("InnerPosZ:fDetected_Position_Z",cond.c_str(),"col") ;
-
- PhysicsTree->Draw("InnerPosZ-fDetected_Position_Z>>h3(1000,-2,2)",cond.c_str(),"col") ;
- c->cd(3);
- c->cd(4);
+ string cond = "InnerPosX!=-1000";
+/*
+ TCanvas* cInner = new TCanvas("ControlInner","ControlInner",1000,1000);
+ cInner->Divide(2,2);
+ cInner->cd(1);
+ PhysicsTree->Draw("InnerPosX:fDetected_Position_X",cond.c_str(),"col") ;
+ cInner->cd(2);
+ PhysicsTree->Draw("InnerPosY:fDetected_Position_Y",cond.c_str(),"col") ;
+ cInner->cd(3);
+ PhysicsTree->Draw("InnerPosZ:fDetected_Position_Z",cond.c_str(),"col") ;
+ cInner->cd(4);
PhysicsTree->Draw("InnerPosY:InnerPosX:InnerPosZ", cond.c_str(),"");
+*/
+ TCanvas* cOuter = new TCanvas("ControlOuter","ControlOuter",1000,1000);
+ cOuter->Divide(2,2);
+ cond = "OuterPosX!=-1000";
+ cOuter->cd(1);
+ PhysicsTree->Draw("OuterPosX:fDetected_Position_X[3]",cond.c_str(),"col") ;
+ cOuter->cd(2);
+ PhysicsTree->Draw("OuterPosY:fDetected_Position_Y[3]",cond.c_str(),"col") ;
+ cOuter->cd(3);
+ PhysicsTree->Draw("OuterPosZ:fDetected_Position_Z[3]",cond.c_str(),"col") ;
+ cOuter->cd(4);
+ PhysicsTree->Draw("OuterPosY:OuterPosX:OuterPosZ", cond.c_str(),"");
+
- /* TCanvas* c2 = new TCanvas("Control 2", "Control2",500,500,2000,1000);
+ TCanvas* c2 = new TCanvas("Control 2", "Control2",500,500,2000,1000);
c2->Divide(2,1);
c2->cd(1);
- PhysicsTree->Draw("InnerPosY:InnerPosX","InnerPosX!=-60");
+ PhysicsTree->Draw("OuterPosY:OuterPosX",cond.c_str());
+ PhysicsTree->Draw("InnerPosY:InnerPosX",cond.c_str(),"same");
c2->cd(2);
- PhysicsTree->Draw("fDetected_Position_Y:fDetected_Position_X","InnerPosX!=-60","");
-*/
+ PhysicsTree->Draw("OuterPosY:OuterPosZ",cond.c_str());
+ PhysicsTree->Draw("InnerPosY:InnerPosZ",cond.c_str(),"same");
}
--
GitLab