PamVMC_update/include/PamVMCSDMgr.h

#ifndef PAMVMCSDMGR_H
#define PAMVMCSDMGR_H
#include <iostream>

#include "TString.h"
#include "TObjString.h"
#include "TMap.h"
#include "PamVMCDetectorSD.h"

#define NTRPL 6 //number of TOF planes

class PamVMCSDMgr: public TObject {

 private:

  static PamVMCSDMgr * fsd;
  TMap fsdmap;

  //some trigger sd's for fast use
  PamVMCDetectorSD* ftrig_sd [NTRPL];
  TString ftrig_colls[NTRPL];
  Bool_t fisgoodfortrig; // boolean veriable to indicate that 
                         //all triggered sd's are present

 protected:
  PamVMCSDMgr() {
    TString tmp[NTRPL]={"S11Y","S12X","S21X","S22Y","S31Y","S32X"};
    for(Int_t i = 0; i<NTRPL; i++) ftrig_colls[i]=tmp[i];
  };
  
 public:
  
  static PamVMCSDMgr * Instance();

  ~PamVMCSDMgr() {fsdmap.DeleteAll();};

  PamVMCDetectorSD * GetSD(const char *name){
    return (PamVMCDetectorSD*)fsdmap(name);
  }

  void SetSD(const char *name, PamVMCDetectorSD *detSD){
    fsdmap.Add(new TObjString(name),detSD);
  }

  void Register(){

      TMapIter *n= (TMapIter *)fsdmap.MakeIterator(); 
      TObject *o; while(( o=(TObject *) n->Next())) { 
        ((PamVMCDetectorSD *)fsdmap.GetValue(o))->Register();
      }
      delete n;
      //check if all trigger SD's are registered. If some TOF layers
      //are missed, return false
      fisgoodfortrig = kTRUE;
      for(Int_t i = 0; i<NTRPL; i++){
        ftrig_sd[i]=(PamVMCDetectorSD *)fsdmap.GetValue(ftrig_colls[i]);
        if(!ftrig_sd[i]) { 
          cout<<"Warning, trigger SD: "<<ftrig_colls[i]<<" is missed, there are no triggers"<<endl;
          fisgoodfortrig = kFALSE;
        }
      }
  }
  void Clear(const Option_t* =""){
    TMapIter *n= (TMapIter *)fsdmap.MakeIterator(); 
    TObject *o; while( (o=(TObject *) n->Next())) { 
    ((PamVMCDetectorSD *)fsdmap.GetValue(o))->ClearHitColl();
    }
    delete n;
  }

  void Compress(const Option_t* =""){
    TMapIter *n= (TMapIter *)fsdmap.MakeIterator(); 
    TObject *o; while( (o=(TObject *) n->Next())) { 
    ((PamVMCDetectorSD *)fsdmap.GetValue(o))->Compress();
    }
    delete n;
  }


  void PreEvent(const Option_t* =""){
    TMapIter *n= (TMapIter *)fsdmap.MakeIterator(); 
    TObject *o; while( (o=(TObject *) n->Next())) { 
    ((PamVMCDetectorSD *)fsdmap.GetValue(o))->PreEvent();
    }
    delete n;
  }

  void Digitize(const Option_t* =""){
    TMapIter *n= (TMapIter *)fsdmap.MakeIterator(); 
    TObject *o; while( (o=(TObject *) n->Next())) { 
    ((PamVMCDetectorSD *)fsdmap.GetValue(o))->Digitize();
    }
    delete n;
  }


  void Print( const Option_t* ="") const { 
    cout << "PamVMCSDMgr, list of SD stored " << endl;
    cout<<"pointers: "<<endl;
    TMapIter *n= (TMapIter *)fsdmap.MakeIterator(); 
    TObject *o; while( (o=(TObject *) n->Next())) { 
      cout<<"->DETECTOR:"<<fsdmap.GetValue(o)<<endl;
      ((PamVMCDetectorSD *)fsdmap.GetValue(o))->Print();
    }
    delete n;
    fsdmap.Print();
  }

  Int_t GetTrigger(){
    if(!fisgoodfortrig) return 0;
    

    Int_t s1or = ftrig_sd[0]->GetNhits() + ftrig_sd[1]->GetNhits();
    Int_t s1an = ftrig_sd[0]->GetNhits() * ftrig_sd[1]->GetNhits();

    Int_t s2or = ftrig_sd[2]->GetNhits() + ftrig_sd[3]->GetNhits();
    Int_t s2an = ftrig_sd[2]->GetNhits() * ftrig_sd[3]->GetNhits();

    Int_t s3or = ftrig_sd[4]->GetNhits() + ftrig_sd[4]->GetNhits();
    Int_t s3an = ftrig_sd[4]->GetNhits() * ftrig_sd[4]->GetNhits();

    //TOF2
    if (s1an*s2an*s3an) return 2;
    //TOF1
    if (s1or*s2or*s3or) return 1;
    //TOF4
    if (s2an*s3an) return 4;

    return 0;
    
  }

};

#endif
1	#ifndef PAMVMCSDMGR_H
2	#define PAMVMCSDMGR_H
3	#include <iostream>
4
5	#include "TString.h"
6	#include "TObjString.h"
7	#include "TMap.h"
8	#include "PamVMCDetectorSD.h"
9
10	#define NTRPL 6 //number of TOF planes
11
12	class PamVMCSDMgr: public TObject {
13
14	private:
15
16	static PamVMCSDMgr * fsd;
17	TMap fsdmap;
18
19	//some trigger sd's for fast use
20	PamVMCDetectorSD* ftrig_sd [NTRPL];
21	TString ftrig_colls[NTRPL];
22	Bool_t fisgoodfortrig; // boolean veriable to indicate that
23	//all triggered sd's are present
24
25	protected:
26	PamVMCSDMgr() {
27	TString tmp[NTRPL]={"S11Y","S12X","S21X","S22Y","S31Y","S32X"};
28	for(Int_t i = 0; i<NTRPL; i++) ftrig_colls[i]=tmp[i];
29	};
30
31	public:
32
33	static PamVMCSDMgr * Instance();
34
35	~PamVMCSDMgr() {fsdmap.DeleteAll();};
36
37	PamVMCDetectorSD * GetSD(const char *name){
38	return (PamVMCDetectorSD*)fsdmap(name);
39	}
40
41	void SetSD(const char name, PamVMCDetectorSD detSD){
42	fsdmap.Add(new TObjString(name),detSD);
43	}
44
45	void Register(){
46
47	TMapIter n= (TMapIter )fsdmap.MakeIterator();
48	TObject o; while(( o=(TObject ) n->Next())) {
49	((PamVMCDetectorSD *)fsdmap.GetValue(o))->Register();
50	}
51	delete n;
52	//check if all trigger SD's are registered. If some TOF layers
53	//are missed, return false
54	fisgoodfortrig = kTRUE;
55	for(Int_t i = 0; i<NTRPL; i++){
56	ftrig_sd[i]=(PamVMCDetectorSD *)fsdmap.GetValue(ftrig_colls[i]);
57	if(!ftrig_sd[i]) {
58	cout<<"Warning, trigger SD: "<<ftrig_colls[i]<<" is missed, there are no triggers"<<endl;
59	fisgoodfortrig = kFALSE;
60	}
61	}
62	}
63	void Clear(const Option_t* =""){
64	TMapIter n= (TMapIter )fsdmap.MakeIterator();
65	TObject o; while( (o=(TObject ) n->Next())) {
66	((PamVMCDetectorSD *)fsdmap.GetValue(o))->ClearHitColl();
67	}
68	delete n;
69	}
70
71	void Compress(const Option_t* =""){
72	TMapIter n= (TMapIter )fsdmap.MakeIterator();
73	TObject o; while( (o=(TObject ) n->Next())) {
74	((PamVMCDetectorSD *)fsdmap.GetValue(o))->Compress();
75	}
76	delete n;
77	}
78
79
80	void PreEvent(const Option_t* =""){
81	TMapIter n= (TMapIter )fsdmap.MakeIterator();
82	TObject o; while( (o=(TObject ) n->Next())) {
83	((PamVMCDetectorSD *)fsdmap.GetValue(o))->PreEvent();
84	}
85	delete n;
86	}
87
88	void Digitize(const Option_t* =""){
89	TMapIter n= (TMapIter )fsdmap.MakeIterator();
90	TObject o; while( (o=(TObject ) n->Next())) {
91	((PamVMCDetectorSD *)fsdmap.GetValue(o))->Digitize();
92	}
93	delete n;
94	}
95
96
97	void Print( const Option_t* ="") const {
98	cout << "PamVMCSDMgr, list of SD stored " << endl;
99	cout<<"pointers: "<<endl;
100	TMapIter n= (TMapIter )fsdmap.MakeIterator();
101	TObject o; while( (o=(TObject ) n->Next())) {
102	cout<<"->DETECTOR:"<<fsdmap.GetValue(o)<<endl;
103	((PamVMCDetectorSD *)fsdmap.GetValue(o))->Print();
104	}
105	delete n;
106	fsdmap.Print();
107	}
108
109	Int_t GetTrigger(){
110	if(!fisgoodfortrig) return 0;
111
112
113	Int_t s1or = ftrig_sd[0]->GetNhits() + ftrig_sd[1]->GetNhits();
114	Int_t s1an = ftrig_sd[0]->GetNhits() * ftrig_sd[1]->GetNhits();
115
116	Int_t s2or = ftrig_sd[2]->GetNhits() + ftrig_sd[3]->GetNhits();
117	Int_t s2an = ftrig_sd[2]->GetNhits() * ftrig_sd[3]->GetNhits();
118
119	Int_t s3or = ftrig_sd[4]->GetNhits() + ftrig_sd[4]->GetNhits();
120	Int_t s3an = ftrig_sd[4]->GetNhits() * ftrig_sd[4]->GetNhits();
121
122	//TOF2
123	if (s1ans2ans3an) return 2;
124	//TOF1
125	if (s1ors2ors3or) return 1;
126	//TOF4
127	if (s2an*s3an) return 4;
128
129	return 0;
130
131	}
132
133	};
134
135	#endif