Skip to content
Snippets Groups Projects
Commit 86c29954 authored by alberica.toia@cern.ch's avatar alberica.toia@cern.ch Committed by Pierre-Alain Loizeau
Browse files

Apply clang-format

parent daf4fc16
No related branches found
No related tags found
1 merge request!295Cosy2019 geometry par fixed
Expand all Hide whitespace changes
Inline Side-by-side
fles/cosy2019/unpacker/CbmStsCoincHodo.cxx
+ 486
− 541
View file @ 86c29954
This diff is collapsed.
fles/cosy2019/unpacker/CbmStsCoincHodo.h
+ 141
− 144
View file @ 86c29954
......@@ -10,157 +10,154 @@
#include "FairTask.h"
#include "TString.h"
#include "TH3.h"
#include "TClonesArray.h"
#include "TH3.h"
#include "TString.h"
//class TClonesArray;
class TH1;
class TH2;
class TH3;
class CbmStsCoincHodo : public FairTask
{
public:
/** Default constructor **/
CbmStsCoincHodo();
CbmStsCoincHodo(const CbmStsCoincHodo&) = delete;
CbmStsCoincHodo operator=(const CbmStsCoincHodo&) = delete;
/** Constructor with parameters (Optional) **/
// CbmStsCoincHodo(Int_t verbose);
/** Destructor **/
~CbmStsCoincHodo();
/** Initiliazation of task at the beginning of a run **/
virtual InitStatus Init();
/** ReInitiliazation of task when the runID changes **/
virtual InitStatus ReInit();
/** Executed for each event. **/
virtual void Exec(Option_t*);
/** Finish task called at the end of the run **/
virtual void Finish();
inline void SetOutFilename( TString sNameIn ) { fOutFileName = sNameIn; }
inline void SetTsStart( Int_t i ) { TsStart = i; }
inline void SetTsStop( Int_t i ) { TsStop = i; }
private:
void CreateHistos();
void WriteHistos();
/** Input array from previous already existing data level **/
TClonesArray* arrayClusters = nullptr;
TClonesArray* arrayHits = nullptr;
/// Constants
Double_t dPosZHodoA = 200.0;
Double_t dPosZSts = 229.0;
Double_t dPosZHodoB = 264.0;
Double_t dHodoDistZ = dPosZHodoB - dPosZHodoA;
Double_t dStsDistZ = dPosZSts - dPosZHodoA;
Double_t dMidStsHodoA = ( dPosZHodoA + dPosZSts ) / 2.0;
Double_t dMidStsHodoB = ( dPosZHodoB + dPosZSts ) / 2.0;
Int_t iCoincLimitClk = 100;
Double_t dClockCycle = 3.125; // ns
Double_t dCoincLimit = iCoincLimitClk * dClockCycle;
//
Int_t fNbTs = 0;
Int_t fNrOfStsDigis = 0;
//
Int_t TsStart;
Int_t TsStop;
// Histos
TH1 * phHitsStsTime;
TH1 * phHitsHodoATime;
TH1 * phHitsHodoBTime;
TH2 * phHitsPositionHodoA;
TH2 * phHitsPositionSts;
TH2 * phHitsPositionHodoB;
TH2 * phNbHitsCompHodo;
TH2 * phNbHitsCompStsHodoA;
TH2 * phNbHitsCompStsHodoB;
TH2 * phHitsCoincCorrXX;
TH2 * phHitsCoincCorrYY;
TH2 * phHitsCoincCorrXY;
TH2 * phHitsCoincCorrYX;
TH2 * phHitsPositionCoincA;
TH2 * phHitsPositionCoincB;
TH2 * phHitsPositionDiff;
TH1 * phHitsTimeDiff;
TH1 * phHitsCoincDist;
TH1 * phHitsCoincAngle;
TH2 * phHitsSingleCoincCorrXX;
TH2 * phHitsSingleCoincCorrYY;
TH2 * phHitsSingleCoincCorrXY;
TH2 * phHitsSingleCoincCorrYX;
TH2 * phHitsSinglePositionCoincA;
TH2 * phHitsSinglePositionCoincB;
TH2 * phHitsSinglePositionDiff;
TH1 * phHitsSingleTimeDiff;
TH1 * phHitsSingleCoincDist;
TH1 * phHitsSingleCoincAngle;
TH2 * phHitsBestCoincCorrXX;
TH2 * phHitsBestCoincCorrYY;
TH2 * phHitsBestCoincCorrXY;
TH2 * phHitsBestCoincCorrYX;
TH2 * phHitsBestPositionCoincA;
TH2 * phHitsBestPositionCoincB;
TH2 * phHitsBestPositionDiff;
TH1 * phHitsBestTimeDiff;
TH1 * phHitsBestCoincDist;
TH1 * phHitsBestCoincAngle;
TH2 * phHitsPositionCoincExtr;
TH2 * phHitsStsCoincCorrXX;
TH2 * phHitsStsCoincCorrYY;
TH2 * phHitsStsCoincCorrXY;
TH2 * phHitsStsCoincCorrYX;
TH2 * phHitsStsPositionCoincExtr;
TH2 * phHitsStsPositionCoinc;
TH1 * phHitsStsTimeDiff;
TH2 * phHitsStsPositionDiff;
TH2 * phHitsStsPositionDiffInv;
TH2 * phHitsStsBestCoincCorrXX;
TH2 * phHitsStsBestCoincCorrYY;
TH2 * phHitsStsBestCoincCorrXY;
TH2 * phHitsStsBestCoincCorrYX;
TH2 * phHitsStsBestPositionCoincExtr;
TH2 * phHitsStsBestPositionCoinc;
TH2 * phHitsStsBestPositionShiftCoinc;
TH1 * phHitsStsBestTimeDiff;
TH2 * phHitsStsBestPositionDiff;
TH2 * phHitsStsBestPositionDiffInv;
TH3 * phHitsStsBestDiff;
TH2 * phHitsStsEff;
TString fOutFileName{"testCoincHodo.root"};
ClassDef(CbmStsCoincHodo,1);
class CbmStsCoincHodo : public FairTask {
public:
/** Default constructor **/
CbmStsCoincHodo();
CbmStsCoincHodo(const CbmStsCoincHodo&) = delete;
CbmStsCoincHodo operator=(const CbmStsCoincHodo&) = delete;
/** Constructor with parameters (Optional) **/
// CbmStsCoincHodo(Int_t verbose);
/** Destructor **/
~CbmStsCoincHodo();
/** Initiliazation of task at the beginning of a run **/
virtual InitStatus Init();
/** ReInitiliazation of task when the runID changes **/
virtual InitStatus ReInit();
/** Executed for each event. **/
virtual void Exec(Option_t*);
/** Finish task called at the end of the run **/
virtual void Finish();
inline void SetOutFilename(TString sNameIn) { fOutFileName = sNameIn; }
inline void SetTsStart(Int_t i) { TsStart = i; }
inline void SetTsStop(Int_t i) { TsStop = i; }
private:
void CreateHistos();
void WriteHistos();
/** Input array from previous already existing data level **/
TClonesArray* arrayClusters = nullptr;
TClonesArray* arrayHits = nullptr;
/// Constants
Double_t dPosZHodoA = 200.0;
Double_t dPosZSts = 229.0;
Double_t dPosZHodoB = 264.0;
Double_t dHodoDistZ = dPosZHodoB - dPosZHodoA;
Double_t dStsDistZ = dPosZSts - dPosZHodoA;
Double_t dMidStsHodoA = (dPosZHodoA + dPosZSts) / 2.0;
Double_t dMidStsHodoB = (dPosZHodoB + dPosZSts) / 2.0;
Int_t iCoincLimitClk = 100;
Double_t dClockCycle = 3.125; // ns
Double_t dCoincLimit = iCoincLimitClk * dClockCycle;
//
Int_t fNbTs = 0;
Int_t fNrOfStsDigis = 0;
//
Int_t TsStart;
Int_t TsStop;
// Histos
TH1* phHitsStsTime;
TH1* phHitsHodoATime;
TH1* phHitsHodoBTime;
TH2* phHitsPositionHodoA;
TH2* phHitsPositionSts;
TH2* phHitsPositionHodoB;
TH2* phNbHitsCompHodo;
TH2* phNbHitsCompStsHodoA;
TH2* phNbHitsCompStsHodoB;
TH2* phHitsCoincCorrXX;
TH2* phHitsCoincCorrYY;
TH2* phHitsCoincCorrXY;
TH2* phHitsCoincCorrYX;
TH2* phHitsPositionCoincA;
TH2* phHitsPositionCoincB;
TH2* phHitsPositionDiff;
TH1* phHitsTimeDiff;
TH1* phHitsCoincDist;
TH1* phHitsCoincAngle;
TH2* phHitsSingleCoincCorrXX;
TH2* phHitsSingleCoincCorrYY;
TH2* phHitsSingleCoincCorrXY;
TH2* phHitsSingleCoincCorrYX;
TH2* phHitsSinglePositionCoincA;
TH2* phHitsSinglePositionCoincB;
TH2* phHitsSinglePositionDiff;
TH1* phHitsSingleTimeDiff;
TH1* phHitsSingleCoincDist;
TH1* phHitsSingleCoincAngle;
TH2* phHitsBestCoincCorrXX;
TH2* phHitsBestCoincCorrYY;
TH2* phHitsBestCoincCorrXY;
TH2* phHitsBestCoincCorrYX;
TH2* phHitsBestPositionCoincA;
TH2* phHitsBestPositionCoincB;
TH2* phHitsBestPositionDiff;
TH1* phHitsBestTimeDiff;
TH1* phHitsBestCoincDist;
TH1* phHitsBestCoincAngle;
TH2* phHitsPositionCoincExtr;
TH2* phHitsStsCoincCorrXX;
TH2* phHitsStsCoincCorrYY;
TH2* phHitsStsCoincCorrXY;
TH2* phHitsStsCoincCorrYX;
TH2* phHitsStsPositionCoincExtr;
TH2* phHitsStsPositionCoinc;
TH1* phHitsStsTimeDiff;
TH2* phHitsStsPositionDiff;
TH2* phHitsStsPositionDiffInv;
TH2* phHitsStsBestCoincCorrXX;
TH2* phHitsStsBestCoincCorrYY;
TH2* phHitsStsBestCoincCorrXY;
TH2* phHitsStsBestCoincCorrYX;
TH2* phHitsStsBestPositionCoincExtr;
TH2* phHitsStsBestPositionCoinc;
TH2* phHitsStsBestPositionShiftCoinc;
TH1* phHitsStsBestTimeDiff;
TH2* phHitsStsBestPositionDiff;
TH2* phHitsStsBestPositionDiffInv;
TH3* phHitsStsBestDiff;
TH2* phHitsStsEff;
TString fOutFileName {"testCoincHodo.root"};
ClassDef(CbmStsCoincHodo, 1);
};
#endif
macro/beamtime/cosy2019/check_coinc_Hodo.C
+ 10
− 14
View file @ 86c29954
void check_coinc_Hodo(
Int_t nEvents = 10,
UInt_t run=25,
UInt_t prefix=0)
void check_coinc_Hodo(Int_t nEvents = 10, UInt_t run = 25, UInt_t prefix = 0)
{
TString runId = TString::Format("%04u", run);
TString runId = TString::Format("%04u", run);
TString prefixId = TString::Format("%04u", prefix);
// TString fileName="/lustre/cbm/users/alberica/cosy2019/rec/TESTrec_cosy_";
// fileName = fileName + runId + "_" + prefixId + ".root";
TString fileName="0025.rec.root";
TString fileName = "0025.rec.root";
// ========================================================================
// Adjust this part according to your requirements
......@@ -22,8 +19,8 @@ void check_coinc_Hodo(
// ------------------------------------------------------------------------
// ----- Analysis run --------------------------------------------------
FairRunOnline *fRun= new FairRunOnline();
fRun->ActivateHttpServer( 100, 8080 ); // refresh each 100 events
FairRunOnline* fRun = new FairRunOnline();
fRun->ActivateHttpServer(100, 8080); // refresh each 100 events
FairFileSource* inputSource = new FairFileSource(fileName);
fRun->SetSource(inputSource);
......@@ -35,18 +32,17 @@ void check_coinc_Hodo(
CbmStsCoincHodo* timeChecker = new CbmStsCoincHodo();
timeChecker->SetTsStart(0);
timeChecker->SetTsStop(5000);
timeChecker->SetOutFilename(Form("testCoincHodo_%i_%i_2_8.root",run,prefix));
timeChecker->SetOutFilename(Form("testCoincHodo_%i_%i_2_8.root", run, prefix));
fRun->AddTask(timeChecker);
// ----- Intialise and run --------------------------------------------
fRun->Init();
cout << "Starting run" << endl;
fRun->Run(0,nEvents);
fRun->Run(0, nEvents);
// ------------------------------------------------------------------------
// ----- Finish -------------------------------------------------------
timer.Stop();
Double_t rtime = timer.RealTime();
......@@ -60,12 +56,12 @@ void check_coinc_Hodo(
// Extract the maximal used memory an add is as Dart measurement
// This line is filtered by CTest and the value send to CDash
FairSystemInfo sysInfo;
Float_t maxMemory=sysInfo.GetMaxMemory();
Float_t maxMemory = sysInfo.GetMaxMemory();
cout << "<DartMeasurement name=\"MaxMemory\" type=\"numeric/double\">";
cout << maxMemory;
cout << "</DartMeasurement>" << endl;
Float_t cpuUsage=ctime/rtime;
Float_t cpuUsage = ctime / rtime;
cout << "<DartMeasurement name=\"CpuLoad\" type=\"numeric/double\">";
cout << cpuUsage;
cout << "</DartMeasurement>" << endl;
......@@ -73,7 +69,7 @@ void check_coinc_Hodo(
FairMonitor* tempMon = FairMonitor::GetMonitor();
tempMon->Print();
// RemoveGeoManager();
// RemoveGeoManager();
cout << " Test passed" << endl;
cout << " All ok " << endl;
}
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment