From 17dfddc8dbe69d20a1cf68d7cb03604aa8e3bd10 Mon Sep 17 00:00:00 2001
From: Unknown <unknown>
Date: Wed, 2 Oct 2013 10:00:25 +0000
Subject: [PATCH]
---
NPAnalysis/Template/src/Analysis.cc | 25 ++++---------------------
1 file changed, 4 insertions(+), 21 deletions(-)
diff --git a/NPAnalysis/Template/src/Analysis.cc b/NPAnalysis/Template/src/Analysis.cc
index 8916b88b8..a22fa769c 100644
--- a/NPAnalysis/Template/src/Analysis.cc
+++ b/NPAnalysis/Template/src/Analysis.cc
@@ -30,16 +30,6 @@ int main(int argc,char** argv)
DetectorManager* myDetector = new DetectorManager ;
myDetector -> ReadConfigurationFile(detectorfileName) ;
- // Attach more branch to the output
- double Ex = 0 ; double EE = 0 ; double TT = 0 ; double X = 0 ; double Y = 0 ; int det ;
- RootOutput::getInstance()->GetTree()->Branch("ExcitationEnergy",&Ex,"Ex/D") ;
- RootOutput::getInstance()->GetTree()->Branch("E",&EE,"EE/D") ;
- RootOutput::getInstance()->GetTree()->Branch("A",&TT,"TT/D") ;
- RootOutput::getInstance()->GetTree()->Branch("X",&X,"X/D") ;
- RootOutput::getInstance()->GetTree()->Branch("Y",&Y,"Y/D") ;
-
- // Get Must2 Pointer:
- MUST2Array* M2 = (MUST2Array*) myDetector -> m_Detector["MUST2"] ;
// Get the formed Chained Tree and Treat it
TChain* Chain = RootInput:: getInstance() -> GetChain() ;
int i;
@@ -51,18 +41,11 @@ int main(int argc,char** argv)
myDetector -> ClearEventPhysics() ;
myDetector -> BuildPhysicalEvent() ;
+ /************************************************
+
+ Put your code here
- double E = M2 -> GetEnergyDeposit();
- TVector3 A = M2 -> GetPositionOfInteraction();
- double Theta = ThetaCalculation ( A , TVector3(0,0,1) ) ;
- if(E>-1000) Ex = myReaction -> ReconstructRelativistic( E , Theta ) ;
- else Ex = -100 ;
- EE = E ; TT = Theta/deg ;
- if(E>-1000){
- X = A . X();
- Y = A . Y();}
- else{X = -1000 ; Y = -1000;}
-
+ ************************************************/
RootOutput::getInstance()->GetTree()->Fill() ;
}
cout << "A total of " << i << " event has been annalysed " << endl ;
--
GitLab